OSCR

MeCP2 gene dosage-dependent neurodevelopmentally restricted defects arise by aberrant activation of cell fate-determining bivalent genes.

Code ↔ Paper

18 matches between paragraphs of the paper and lines of its authors' code, computed by the harvester (lexical-v1). Click a colored paragraph or line to see its counterpart.

The 18 matches · 2 of them tie a paragraph to a whole file, not to given lines: weak matches, whose lines are not tinted
  1. [1] § Methods › Computational analysis › RNA-seq data analysis and integration ↔ 03_RNAseq_analysis/02_star_alignment.sh, lines 46–99 · score 0.98 · GeneCounts, SortedByCoordinate, alignIntronMax, alignSJoverhangMin, outFilterMismatchNoverReadLmax, outSAMtype
  2. [2] § Methods › Computational analysis › SAMMY-seq bioinformatics analysis ↔ 06_SAMMY_seq/01_process_fastq.sh, lines 193–275 · score 0.97 · MarkDuplicates, PCR duplicates, extendReads, normalizeUsing, mapping quality, binSize
  3. [3] § Methods › Computational analysis › SAMMY-seq bioinformatics analysis ↔ 06_SAMMY_seq/02b_sammy_analysis_blacklisted.py, lines 150–202 · score 0.89 · extendReads, binSize, bamCoverage, RPKM normalization, BigWig, SAMMY seq
  4. [4] § Methods › Computational analysis › CUT&Tag data processing and analysis ↔ 06_SAMMY_seq/01_process_fastq.sh, lines 193–275 · score 0.85 · PCR duplicates, mapping quality, bamCoverage, bigWig, Alignments, Trim
  5. [5] § Methods › Computational analysis › CUT&Tag data processing and analysis ↔ 02_CpG_enrichment/integrate_rnaseq.R, lines 41–97 · score 0.83 · transcript database, TxDb.Mmusculus.UCSC.mm10.knownGene, nearest genes, TSS regions, mapping, Peak
  6. [6] § Methods › Computational analysis › RNA-seq data analysis and integration ↔ 03_RNAseq_analysis/04_deseq2_analysis.R, lines 149–215 · score 0.82 · variance stabilizing transformation, DESeq2, log2FC, log2 fold change, VST, padj
  7. [7] § Methods › Computational analysis › CUT&Tag data processing and analysis ↔ 04_binding_expression_integration/integrate_genebody_rnaseq.py, lines 177–237 · score 0.80 · binding expression correlation, exogenous signal, endogenous signal, log2 fold changes, gene expression, MeCP2
  8. [8] § Methods › Computational analysis › CUT&Tag data processing and analysis ↔ 04_binding_expression_integration/peak_expression_correlation.py, lines 225–289 · score 0.79 · peak enrichment ratios, binding expression correlation, log2 fold changes, MeCP2 enrichment, aggregated, DEA
  9. [9] § Methods › Computational analysis › Smarcb1 binding analysis and integration ↔ 07_visualization/01_metaprofiles.sh, lines 47–124 · score 0.79 · plotHeatmap, computeMatrix, plotProfile, Metaprofiles, width, bp
  10. [10] § Methods › Computational analysis › CUT&Tag data processing and analysis ↔ 02_CpG_enrichment/cpg_enrichment_analysis.py, lines 1–24 · score 0.73 · pyBigWig, signal intensities, RPKM normalized, quantification, island, enrichment
  11. [11] § Methods › Computational analysis › CUT&Tag data processing and analysis ↔ 04_binding_expression_integration/integrate_rnaseq_base.py, lines 61–113 · score 0.70 · nearest genes, log2 fold change, endo signal, closest, RNA seq, padj
  12. [12] § Results › Mecp2 localization into CpG islands in both NPCs and Ns ↔ 06_SAMMY_seq/07_target_comparison_heatmap.py, lines 210–339 · score 0.57 · transcriptional end site, transcriptional start sites, located, heatmaps, upstream, genomic
  13. [13] § Results › Differential gene deregulation and phenotypic changes induced by Mecp2 overexpression in primary neuronal progenitor cells and neurons ↔ 03_RNAseq_analysis/05_functional_enrichment.R, lines 43–92 · score 0.57 · Gene Ontology, upregulated genes, downregulated genes, fold change, GO, enriched
  14. [14] § Results › Differential occupancy of CpG islands by Mecp2-endo and -exo in NPCs and Ns ↔ 03_RNAseq_analysis/05_functional_enrichment.R, lines 43–92 · score 0.55 · Gene Ontology, upregulated genes, downregulated genes, GO, enriched, enrichment
  15. [15] § Results › Differential gene deregulation and phenotypic changes induced by Mecp2 overexpression in primary neuronal progenitor cells and neurons ↔ 03_RNAseq_analysis/generate_de_table.R, lines 150–188 · score 0.54 · Benjamini Hochberg, detected genes, log2FC, log2 fold changes, protein, downregulated
  16. [16] § Results › Exogenous Mecp2 upregulated bivalent genes in NPCs by interacting with the SWI/SNF complex ↔ 04_binding_expression_integration/integrate_rnaseq_base.py, lines 268–322 · score 0.54 · Mecp2 enriched genes, Mecp2 enrichment, RNA seq, fold change, endo, CpG
  17. [17] § Results › Exogenous Mecp2 exerted a more pronounced effect on the chromatin state in NPCs than in Ns ↔ 06_SAMMY_seq/02_run_sammy_analysis.sh, the whole file · a weak match · score 0.53 · chromatin fraction, SAMMY seq, chromatin state, accessibility, S2S, S3
  18. [18] § Methods › Computational analysis › SAMMY-seq bioinformatics analysis ↔ 06_SAMMY_seq/02_run_sammy_analysis.sh, the whole file · a weak match · score 0.51 · chromatin fractions, SAMMY seq, S2S, S3

Paper

Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC

The paper is loaded when this pane is shown.

The authors' code

Shell · 275 lines · 9.2 KB · no license · 2 matches

  1. #!/bin/bash
  2. #SBATCH --job-name=process_samples
  3. #SBATCH --output=logs/process_samples_%a.out
  4. #SBATCH --error=logs/process_samples_%a.err
  5. #SBATCH --nodes=1
  6. #SBATCH --ntasks=1
  7. #SBATCH --cpus-per-task=16
  8. #SBATCH --mem=64G
  9. #SBATCH --time=24:00:00
  10. #SBATCH --partition=workq
  11. #SBATCH --array=1-34%10 # Update based on the number of files
  12. # Error handling
  13. set -e # Exit on error
  14. set -u # Exit on undefined variable
  15. # Activate conda environment with required tools (BWA, samtools, picard, bamCoverage)
  16. conda activate snakemake # Adjust to your environment name
  17. # Set WORKDIR to your working directory
  18. WORKDIR="" # e.g., "/path/to/sammy_seq/analysis"
  19. cd $WORKDIR
  20. # Directory setup
  21. FASTQ_DIR="${WORKDIR}/../DATA/Sammy_Seq_fastq"
  22. OUTPUT_DIR="${WORKDIR}/results"
  23. # Need to set up BWA index for mm10 genome
  24. BWA_GENOME_DIR="${WORKDIR}/mm10"
  25. BLACKLIST_FILE="${WORKDIR}/mm10-blacklist.v2.bed"
  26. THREADS=16 # Using all CPUs allocated to the job
  27. # Create a directory for trimmed reads
  28. mkdir -p $OUTPUT_DIR/trimmed
  29. mkdir -p $OUTPUT_DIR/bam/temp
  30. mkdir -p logs
  31. # Download blacklist file if it doesn't exist
  32. if [ ! -f "$BLACKLIST_FILE" ]; then
  33. echo "Downloading mm10 blacklist file..."
  34. mkdir -p $(dirname $BLACKLIST_FILE)
  35. wget -O $BLACKLIST_FILE https://www.encodeproject.org/files/ENCFF547MET/@@download/ENCFF547MET.bed.gz
  36. gunzip -f $BLACKLIST_FILE.gz
  37. fi
  38. # Check if the FASTQ directory exists
  39. if [ ! -d "$FASTQ_DIR" ]; then
  40. echo "Error: FASTQ directory $FASTQ_DIR does not exist"
  41. exit 1
  42. fi
  43. echo "Using FASTQ directory: $FASTQ_DIR"
  44. # Create output directories
  45. mkdir -p $OUTPUT_DIR/bam
  46. mkdir -p logs
  47. # Check if the filtered_sample_list.txt exists
  48. if [ ! -f "filtered_sample_list.txt" ]; then
  49. echo "Error: filtered_sample_list.txt not found"
  50. echo "Please run check_number_of_files.sh first to generate the file list"
  51. exit 1
  52. fi
  53. # Count how many files we have
  54. FILE_COUNT=$(wc -l < filtered_sample_list.txt)
  55. echo "Found $FILE_COUNT files in the filtered sample list"
  56. # Display the first few files to verify
  57. echo -e "\nFirst 10 files to process:"
  58. head -n 10 filtered_sample_list.txt
  59. # Check if the array size matches the number of files
  60. if [[ "$SLURM_ARRAY_TASK_MAX" ]]; then
  61. ARRAY_SIZE=$((SLURM_ARRAY_TASK_MAX))
  62. if [ "$ARRAY_SIZE" != "$FILE_COUNT" ]; then
  63. echo "Warning: Array size ($ARRAY_SIZE) does not match the number of files ($FILE_COUNT)"
  64. echo "For optimal processing, you should cancel this job and resubmit with:"
  65. echo "sbatch --array=1-${FILE_COUNT}%10 $0"
  66. fi
  67. else
  68. echo "Note: Running in non-array mode or array information not available"
  69. fi
  70. # Get the FASTQ file for this array job
  71. FASTQ=$(sed -n "${SLURM_ARRAY_TASK_ID}p" filtered_sample_list.txt)
  72. if [ -z "$FASTQ" ]; then
  73. echo "Error: No file found for task ID ${SLURM_ARRAY_TASK_ID}"
  74. exit 1
  75. fi
  76. FASTQ="$FASTQ_DIR/$FASTQ"
  77. if [ ! -f "$FASTQ" ]; then
  78. echo "Error: FASTQ file $FASTQ not found"
  79. exit 1
  80. fi
  81. # Get base name for output files
  82. BASE=$(basename $FASTQ _R1_001.fastq.gz)
  83. echo "Processing $BASE..."
  84. # Create a directory to store trimmed adapter files if it doesn't exist
  85. ADAPTER_DIR="$WORKDIR/adapters"
  86. mkdir -p $ADAPTER_DIR
  87. ADAPTER_FILE="$ADAPTER_DIR/TruSeq3-SE-2.fa"
  88. # Create adapter file if it doesn't exist
  89. if [ ! -f "$ADAPTER_FILE" ]; then
  90. echo "Creating TruSeq3-SE-2.fa adapter file..."
  91. cat > "$ADAPTER_FILE" << 'EOL'
  92. >TruSeq3_IndexedAdapter
  93. AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
  94. >TruSeq3_UniversalAdapter
  95. AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGTA
  96. EOL
  97. fi
  98. # Step 1: Trim reads using Trimmomatic as specified in the paper
  99. TRIMMED_FASTQ="$OUTPUT_DIR/trimmed/${BASE}.trimmed.fastq.gz"
  100. if [ -f "$TRIMMED_FASTQ" ]; then
  101. echo "Trimmed file already exists, skipping trimming step..."
  102. else
  103. echo "Trimming reads with Trimmomatic..."
  104. trimmomatic SE $FASTQ $TRIMMED_FASTQ \
  105. ILLUMINACLIP:$ADAPTER_FILE:2:30:10 \
  106. SLIDINGWINDOW:4:15 \
  107. MINLEN:36
  108. fi
  109. # Step 2: Align with BWA as specified in the paper
  110. echo "Aligning with BWA..."
  111. # Check if BWA index exists and is valid, if not create it
  112. BWA_INDEX_FILES=("$BWA_GENOME_DIR/mm10.amb" "$BWA_GENOME_DIR/mm10.ann" "$BWA_GENOME_DIR/mm10.bwt" "$BWA_GENOME_DIR/mm10.pac" "$BWA_GENOME_DIR/mm10.sa")
  113. INDEX_VALID=true
  114. # Check if all index files exist and have non-zero size
  115. for index_file in "${BWA_INDEX_FILES[@]}"; do
  116. if [ ! -s "$index_file" ]; then
  117. INDEX_VALID=false
  118. echo "BWA index file $index_file is missing or empty."
  119. break
  120. fi
  121. done
  122. # Even if all files exist, perform a more thorough validation by attempting to use the index
  123. if [ "$INDEX_VALID" = true ]; then
  124. echo "Testing BWA index validity..."
  125. # Create a small test FASTQ with a few reads
  126. TEST_FASTQ="${BWA_GENOME_DIR}/test_reads.fastq"
  127. echo -e "@test_read\nACGTACGTACGTACGTACGT\n+\nIIIIIIIIIIIIIIIIIIII" > "$TEST_FASTQ"
  128. # Try to align the test read - redirect stderr to capture any errors
  129. TEST_OUTPUT=$(bwa mem "$BWA_GENOME_DIR/mm10" "$TEST_FASTQ" 2>&1 || true)
  130. # Check if the output contains error messages related to index files
  131. if echo "$TEST_OUTPUT" | grep -q -E "(Can't|Cannot|failed|error|Unexpected end of file)"; then
  132. INDEX_VALID=false
  133. echo "BWA index validation failed with error:"
  134. echo "$TEST_OUTPUT" | grep -E "(Can't|Cannot|failed|error|Unexpected end of file)"
  135. else
  136. echo "BWA index validation successful."
  137. fi
  138. # Clean up test file
  139. rm -f "$TEST_FASTQ"
  140. fi
  141. if [ "$INDEX_VALID" = false ]; then
  142. echo "BWA index is invalid or incomplete. Removing old index files and recreating..."
  143. # Remove any existing index files that might be corrupted
  144. rm -f ${BWA_INDEX_FILES[@]}
  145. # Also remove any potential temporary files that might have been created during a failed indexing
  146. rm -f "$BWA_GENOME_DIR/mm10."*
  147. mkdir -p $BWA_GENOME_DIR
  148. # Check if reference genome exists
  149. if [ ! -f "$WORKDIR/references/mm10.fa" ]; then
  150. echo "Error: Reference genome file $WORKDIR/references/mm10.fa not found"
  151. exit 1
  152. fi
  153. # Check available disk space before indexing
  154. GENOME_SIZE=$(du -b "$WORKDIR/references/mm10.fa" | cut -f1)
  155. REQUIRED_SPACE=$((GENOME_SIZE * 5)) # BWA index typically requires ~5x the genome size
  156. AVAILABLE_SPACE=$(df -B1 "$BWA_GENOME_DIR" | awk 'NR==2 {print $4}')
  157. if [ "$AVAILABLE_SPACE" -lt "$REQUIRED_SPACE" ]; then
  158. echo "Error: Not enough disk space for BWA indexing. Need at least $((REQUIRED_SPACE/1024/1024/1024)) GB, but only $((AVAILABLE_SPACE/1024/1024/1024)) GB available."
  159. exit 1
  160. fi
  161. echo "Creating BWA index. This may take a while..."
  162. bwa index -p $BWA_GENOME_DIR/mm10 "$WORKDIR/references/mm10.fa"
  163. # Verify index was created successfully
  164. for index_file in "${BWA_INDEX_FILES[@]}"; do
  165. if [ ! -s "$index_file" ]; then
  166. echo "Error: Failed to create BWA index file $index_file"
  167. exit 1
  168. fi
  169. done
  170. else
  171. echo "BWA index is valid and complete."
  172. fi
  173. # Check if BAM file already exists
  174. if [ -f "$OUTPUT_DIR/bam/${BASE}.bam" ]; then
  175. echo "BAM file already exists, skipping alignment step..."
  176. else
  177. # Align with BWA using parameters from the paper
  178. echo "Performing BWA alignment..."
  179. bwa mem -t $THREADS $BWA_GENOME_DIR/mm10 $OUTPUT_DIR/trimmed/${BASE}.trimmed.fastq.gz > $OUTPUT_DIR/bam/${BASE}.sam
  180. # Convert SAM to BAM
  181. echo "Converting SAM to BAM..."
  182. samtools view -bS $OUTPUT_DIR/bam/${BASE}.sam > $OUTPUT_DIR/bam/${BASE}.bam
  183. rm $OUTPUT_DIR/bam/${BASE}.sam
  184. fi
  185. # Sort BAM if it exists but hasn't been sorted yet
  186. if [ -f "$OUTPUT_DIR/bam/${BASE}.bam" ]; then
  187. # Check if BAM is already sorted
  188. if ! samtools view -H "$OUTPUT_DIR/bam/${BASE}.bam" | grep -q 'SO:coordinate'; then
  189. echo "Sorting BAM..."
  190. samtools sort -@ $THREADS $OUTPUT_DIR/bam/${BASE}.bam -o $OUTPUT_DIR/bam/${BASE}.sorted.bam
  191. mv $OUTPUT_DIR/bam/${BASE}.sorted.bam $OUTPUT_DIR/bam/${BASE}.bam
  192. else
  193. echo "BAM file already sorted, skipping sorting step..."
  194. fi
  195. fi
  196. # Step 3: Mark and remove PCR duplicates using Picard
  197. if [ -f "$OUTPUT_DIR/bam/${BASE}.dedup.bam" ]; then
  198. echo "Deduplicated BAM already exists, skipping deduplication step..."
  199. else
  200. echo "Marking and removing PCR duplicates with Picard..."
  201. picard MarkDuplicates \
  202. I=$OUTPUT_DIR/bam/${BASE}.bam \
  203. O=$OUTPUT_DIR/bam/${BASE}.dedup.bam \
  204. M=$OUTPUT_DIR/bam/${BASE}.metrics.txt \
  205. REMOVE_DUPLICATES=true \
  206. TMP_DIR=$OUTPUT_DIR/bam/temp
  207. fi
  208. # Step 4: Filter by mapping quality (q>=1) as in the paper
  209. if [ -f "$OUTPUT_DIR/bam/${BASE}.filtered.bam" ]; then
  210. echo "Filtered BAM already exists, skipping filtering step..."
  211. else
  212. echo "Filtering by mapping quality..."
  213. samtools view -q 1 -b $OUTPUT_DIR/bam/${BASE}.dedup.bam > $OUTPUT_DIR/bam/${BASE}.filtered.bam
  214. fi
  215. # Index BAM if needed
  216. if [ ! -f "$OUTPUT_DIR/bam/${BASE}.filtered.bam.bai" ]; then
  217. echo "Indexing BAM..."
  218. samtools index $OUTPUT_DIR/bam/${BASE}.filtered.bam
  219. else
  220. echo "BAM index already exists, skipping indexing step..."
  221. fi
  222. # Step 5: Create bigWig file
  223. if [ -f "$OUTPUT_DIR/bigwig/${BASE}.bw" ]; then
  224. echo "BigWig file already exists, skipping bigWig creation step..."
  225. else
  226. echo "Creating bigWig for $BASE..."
  227. bamCoverage \
  228. --bam $OUTPUT_DIR/bam/${BASE}.filtered.bam \
  229. --outFileName $OUTPUT_DIR/bigwig/${BASE}.bw \
  230. --binSize 50 \
  231. --normalizeUsing RPKM \
  232. --effectiveGenomeSize 2652783500 \
  233. --blackListFileName $BLACKLIST_FILE \
  234. --extendReads 250 \
  235. --numberOfProcessors $THREADS
  236. fi
  237. echo "Processing complete for $BASE"

01_process_fastq.sh at commit a378313, no license · at the source

Overview

Authors: Mirko Luoni1,2, Michal Kubacki2, Serena Gea Giannelli2, Francesca Morosi1,2, Claudia Di Berardino2, Angelo Iannielli1,2, Alessandro Sessa3, Gaia Colasante2, Emanuele Di Patrizio Soldateschi4, Chiara Lanzuolo4,5, Vania Broccoli1,2
  1. Institute of Neuroscience, National Research Council (CNR),Milan, Italy
  2. Stem Cell and Neurogenesis Unit, Division of Neuroscience, IRCCS San Raffaele Scientific Institute,Milan, Italy
  3. Neuroepigenetics Unit, Division of Neuroscience, IRCCS San Raffaele Scientific Institute,Milan, Italy
  4. INGM Istituto Nazionale Genetica Molecolare Romeo ed Enrica Invernizzi,Milan, Italy
  5. Institute of Biomedical Technologies, National Research Council (CNR),Milan, Italy
Journal: Nature communications, volume 17, issue 1, article 3225
Dates: received 2 June 2025; accepted 4 March 2026; published online 23 April 2026
Type: Research article · Language: English
License: CC BY-NC-ND
Identifiers: DOI 10.1038/s41467-026-71432-w · PMID 42026056 · PMCID PMC13106735 · OpenAlex W7155387614
Open access: gold, a free copy (OpenAlex)
Status: code verified
Categories: human (organism), mouse (organism), other condition (population), cellular / molecular (subfield)
Methods: Spectral & time-frequency, Statistics, Smoothing, state filtering, decompositions, Preprocessing, Connectivity
Keywords: Epigenetics, Molecular neuroscience, Gene expression, Gene regulation
MeSH: Gene Dosage*, Methyl-CpG-Binding Protein 2*, Animals, CpG Islands, Humans, Mice, Neural Stem Cells, Neurodevelopment, Neurons, Rett Syndrome, Transcriptional Activation (* major topic)
Topic: Genetics and Neurodevelopmental Disorders (Genetics, Biochemistry, Genetics and Molecular Biology), according to OpenAlex
Citations: cited by 1 paper (Europe PMC); 89 references in the paper

Abstract

The abstract is not reproduced here: the paper's license (CC BY-NC-ND) does not allow it. Read it in the paper, at the publisher or on Europe PMC.

Repositories

Its files are read in the Code ↔ Paper reader above, with 18 matches between paragraphs and lines of code.

Zenodo 18257503

License: CC-BY-4.0
State: the link answers, verified on 29 September 2026
Evidence: files inventoried
Size: 1 file
Software Heritage: not checked
Found in: the text, “CUT&Tag data processing and analysis”
Not found: README, license file, CITATION.cff, environment file, tests, continuous integration, documentation
Tools: pandas (23 files), Matplotlib (21 files), NumPy (21 files), seaborn (18 files), SciPy (14 files), tidyverse (8 files), BEDTools (5 files), deepTools (5 files), SAMtools (5 files), clusterProfiler (3 files), DESeq2 (2 files), ggplot2 (2 files), pheatmap (2 files), FastQC (1 file), scikit-learn (1 file), STAR (1 file), statsmodels (1 file), Subread (featureCounts) (1 file)
Availability: 1 check, the latest on 29 September 2026: the link answers (HTTP 200)
  • 29 September 2026: the link answers (HTTP 200)
51 files
At the source:

u9856112956-prog/mecp2_analysis

License: none: the authors keep all their rights
State: the link answers, verified on 29 September 2026
Evidence: files inventoried
Commit: a378313b060fefaa5e141b10bffedc0a6e2f36bd, 15 January 2026
Languages: Python (23), R (14), Shell (13)
Size: 51 files, 50 scripts
Software Heritage: not archived
Found in: the Zenodo archive record
Holds: README
Not found: license file, CITATION.cff, environment file, tests, continuous integration, documentation
Tools: pandas (23 files), Matplotlib (21 files), NumPy (21 files), seaborn (18 files), SciPy (14 files), tidyverse (8 files), BEDTools (5 files), deepTools (5 files), SAMtools (5 files), clusterProfiler (3 files), DESeq2 (2 files), ggplot2 (2 files), pheatmap (2 files), FastQC (1 file), scikit-learn (1 file), STAR (1 file), statsmodels (1 file), Subread (featureCounts) (1 file)
Availability: 1 check, the latest on 29 September 2026: the link answers
  • 29 September 2026: the link answers
51 files

The paper's code and data availability statement is in the Data section.

Tracing map

Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.

What the map holds:

  • 2 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
  • 100 scripts, each with its path and the digest of its content;
  • 18 matches between paragraphs of the paper and lines of the code (method lexical-v1);
  • neither the text of the paper nor the code itself.

Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.

Data

Datasets cited

Code and data availability statement

The paper has a code and data availability statement. Its license (CC BY-NC-ND) does not allow reproducing it here; in short, from what the harvester recognized in it:

Read it in the paper: doi.org/10.1038/s41467-026-71432-w.

Versions

The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.

Version 1, 29 September 2026: the first record

Recorded: type, language, journal, volume, issue, pages, dates, 11 authors, 4 keywords, 11 MeSH terms, 1 funder, 88 references.

Cite

This paper

Luoni, M., Kubacki, M., Giannelli, S. G., Morosi, F., Di Berardino, C., Iannielli, A., Sessa, A., Colasante, G., Di Patrizio Soldateschi, E., Lanzuolo, C., & Broccoli, V. (2026). MeCP2 gene dosage-dependent neurodevelopmentally restricted defects arise by aberrant activation of cell fate-determining bivalent genes. Nature communications, 17(1), 3225. https://doi.org/10.1038/s41467-026-71432-w

BibTeX

@article{luoni2026mecp2,
author = {Luoni, Mirko and Kubacki, Michal and Giannelli, Serena Gea and Morosi, Francesca and Di Berardino, Claudia and Iannielli, Angelo and Sessa, Alessandro and Colasante, Gaia and Di Patrizio Soldateschi, Emanuele and Lanzuolo, Chiara and Broccoli, Vania},
title = {{MeCP2 gene dosage-dependent neurodevelopmentally restricted defects arise by aberrant activation of cell fate-determining bivalent genes}},
journal = {Nature communications},
year = {2026},
month = apr,
volume = {17},
number = {1},
pages = {3225},
publisher = {Nature Publishing Group},
issn = {2041-1723},
doi = {10.1038/s41467-026-71432-w},
url = {https://doi.org/10.1038/s41467-026-71432-w},
pmid = {42026056},
pmcid = {PMC13106735}
}

RIS

TY - JOUR
AU - Luoni, Mirko
AU - Kubacki, Michal
AU - Giannelli, Serena Gea
AU - Morosi, Francesca
AU - Di Berardino, Claudia
AU - Iannielli, Angelo
AU - Sessa, Alessandro
AU - Colasante, Gaia
AU - Di Patrizio Soldateschi, Emanuele
AU - Lanzuolo, Chiara
AU - Broccoli, Vania
TI - MeCP2 gene dosage-dependent neurodevelopmentally restricted defects arise by aberrant activation of cell fate-determining bivalent genes
T2 - Nature communications
J2 - Nat Commun
PY - 2026
DA - 2026/04/23
VL - 17
IS - 1
SP - 3225
SN - 2041-1723
PB - Nature Publishing Group
DO - 10.1038/s41467-026-71432-w
UR - https://doi.org/10.1038/s41467-026-71432-w
LA - en
ER -

CSL-JSON

{
"id": "10.1038/s41467-026-71432-w",
"type": "article-journal",
"title": "MeCP2 gene dosage-dependent neurodevelopmentally restricted defects arise by aberrant activation of cell fate-determining bivalent genes",
"container-title": "Nature communications",
"author": [
{
"family": "Luoni",
"given": "Mirko"
},
{
"family": "Kubacki",
"given": "Michal"
},
{
"family": "Giannelli",
"given": "Serena Gea"
},
{
"family": "Morosi",
"given": "Francesca"
},
{
"family": "Di Berardino",
"given": "Claudia"
},
{
"family": "Iannielli",
"given": "Angelo"
},
{
"family": "Sessa",
"given": "Alessandro"
},
{
"family": "Colasante",
"given": "Gaia"
},
{
"family": "Di Patrizio Soldateschi",
"given": "Emanuele"
},
{
"family": "Lanzuolo",
"given": "Chiara"
},
{
"family": "Broccoli",
"given": "Vania"
}
],
"container-title-short": "Nat Commun",
"volume": "17",
"issue": "1",
"page": "3225",
"DOI": "10.1038/s41467-026-71432-w",
"PMID": "42026056",
"PMCID": "PMC13106735",
"ISSN": "2041-1723",
"publisher": "Nature Publishing Group",
"URL": "https://doi.org/10.1038/s41467-026-71432-w",
"language": "en",
"issued": {
"date-parts": [
[
2026,
4,
23
]
]
}
}

The tracing map gets a citation of its own once an author has validated it and it has a DOI.

Similar papers

The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.

[1] doi:10.1038/s41586-026-10832-w [code]
Subnuclear genome compartmentalization controls bivalent chromatin activity.
Journal: Nature
In common: deepTools, FastQC, Subread (featureCounts), 9 other tools, 4 references
[2] doi:10.21203/rs.3.rs-9927928/v1 [code]
Genome-wide and allele-resolved maps of the radial architecture of the mouse genome
Journal: Research Square (preprint)
In common: FastQC, Subread (featureCounts), STAR, 11 other tools, mouse, 1 reference
[3] doi:10.1126/sciadv.aed2952 [code]
Activation of transposable elements is linked to a region- and cell type-specific interferon response in Parkinson's disease.
Journal: Science advances
In common: deepTools, Subread (featureCounts), STAR, 12 other tools, cellular / molecular
[4] doi:10.1093/bioinformatics/btag592 [code]
Network-based stratification of allele-specific expression reveals patient subgroups in Huntington's disease.
Journal: Bioinformatics (Oxford, England)
In common: FastQC, Subread (featureCounts), STAR, 12 other tools, other condition
[5] doi:10.1038/s41586-026-10512-9 [code]
Astrocyte glucocorticoid receptor signalling restricts neuronal plasticity.
Journal: Nature
In common: Subread (featureCounts), STAR, BEDTools, 10 other tools, mouse, cellular / molecular
[6] doi:10.1038/s41467-026-72598-z [code]
Functional impact of genetic background on variable expressivity in neurodevelopmental disorders.
Journal: Nature communications
In common: deepTools, BEDTools, SAMtools, 11 other tools, other condition
[7] doi:10.7554/elife.107393 [code]
Chromosome-scale genome assembly of the European common cuttlefish &lt;i&gt;Sepia officinalis&lt;/i&gt;.
Journal: eLife
In common: Subread (featureCounts), STAR, BEDTools, 10 other tools, cellular / molecular
[8] doi:10.1038/s41531-026-01287-x [code]
Faecalibacterium prausnitzii, depleted in the Parkinson's disease microbiome, improves motor deficits in α-synuclein overexpressing mice.
Journal: NPJ Parkinson's disease
In common: FastQC, Subread (featureCounts), STAR, 8 other tools, mouse, cellular / molecular
[9] doi:10.1016/j.xhgg.2026.100652 [code]
CRISPR-engineered deletion of POGZ alters transcription factor binding at promoters of genes involved in synaptic signaling.
Journal: HGG advances
In common: STAR, BEDTools, SAMtools, 8 other tools, cellular / molecular, 2 references
[10] doi:10.1038/s42003-026-10957-8 [code]
Brain defence by the extracellular matrix protein Cochlin.
Journal: Communications biology
In common: Subread (featureCounts), SAMtools, DESeq2, 11 other tools, mouse, cellular / molecular

Contribute

The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.

Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.

Request its removal

To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).

Discussion, reproductions, activity

Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.

Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.

Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.