MeCP2 gene dosage-dependent neurodevelopmentally restricted defects arise by aberrant activation of cell fate-determining bivalent genes.
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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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] § 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
- #!/bin/bash
- #SBATCH --job-name=process_samples
- #SBATCH --output=logs/process_samples_%a.out
- #SBATCH --error=logs/process_samples_%a.err
- #SBATCH --nodes=1
- #SBATCH --ntasks=1
- #SBATCH --cpus-per-task=16
- #SBATCH --mem=64G
- #SBATCH --time=24:00:00
- #SBATCH --partition=workq
- #SBATCH --array=1-34%10 # Update based on the number of files
- # Error handling
- set -e # Exit on error
- set -u # Exit on undefined variable
- # Activate conda environment with required tools (BWA, samtools, picard, bamCoverage)
- conda activate snakemake # Adjust to your environment name
- # Set WORKDIR to your working directory
- WORKDIR="" # e.g., "/path/to/sammy_seq/analysis"
- cd $WORKDIR
- # Directory setup
- FASTQ_DIR="${WORKDIR}/../DATA/Sammy_Seq_fastq"
- OUTPUT_DIR="${WORKDIR}/results"
- # Need to set up BWA index for mm10 genome
- BWA_GENOME_DIR="${WORKDIR}/mm10"
- BLACKLIST_FILE="${WORKDIR}/mm10-blacklist.v2.bed"
- THREADS=16 # Using all CPUs allocated to the job
- # Create a directory for trimmed reads
- mkdir -p $OUTPUT_DIR/trimmed
- mkdir -p $OUTPUT_DIR/bam/temp
- mkdir -p logs
- # Download blacklist file if it doesn't exist
- if [ ! -f "$BLACKLIST_FILE" ]; then
- echo "Downloading mm10 blacklist file..."
- mkdir -p $(dirname $BLACKLIST_FILE)
- wget -O $BLACKLIST_FILE https://www.encodeproject.org/files/ENCFF547MET/@@download/ENCFF547MET.bed.gz
- gunzip -f $BLACKLIST_FILE.gz
- fi
- # Check if the FASTQ directory exists
- if [ ! -d "$FASTQ_DIR" ]; then
- echo "Error: FASTQ directory $FASTQ_DIR does not exist"
- exit 1
- fi
- echo "Using FASTQ directory: $FASTQ_DIR"
- # Create output directories
- mkdir -p $OUTPUT_DIR/bam
- mkdir -p logs
- # Check if the filtered_sample_list.txt exists
- if [ ! -f "filtered_sample_list.txt" ]; then
- echo "Error: filtered_sample_list.txt not found"
- echo "Please run check_number_of_files.sh first to generate the file list"
- exit 1
- fi
- # Count how many files we have
- FILE_COUNT=$(wc -l < filtered_sample_list.txt)
- echo "Found $FILE_COUNT files in the filtered sample list"
- # Display the first few files to verify
- echo -e "\nFirst 10 files to process:"
- head -n 10 filtered_sample_list.txt
- # Check if the array size matches the number of files
- if [[ "$SLURM_ARRAY_TASK_MAX" ]]; then
- ARRAY_SIZE=$((SLURM_ARRAY_TASK_MAX))
- if [ "$ARRAY_SIZE" != "$FILE_COUNT" ]; then
- echo "Warning: Array size ($ARRAY_SIZE) does not match the number of files ($FILE_COUNT)"
- echo "For optimal processing, you should cancel this job and resubmit with:"
- echo "sbatch --array=1-${FILE_COUNT}%10 $0"
- fi
- else
- echo "Note: Running in non-array mode or array information not available"
- fi
- # Get the FASTQ file for this array job
- FASTQ=$(sed -n "${SLURM_ARRAY_TASK_ID}p" filtered_sample_list.txt)
- if [ -z "$FASTQ" ]; then
- echo "Error: No file found for task ID ${SLURM_ARRAY_TASK_ID}"
- exit 1
- fi
- FASTQ="$FASTQ_DIR/$FASTQ"
- if [ ! -f "$FASTQ" ]; then
- echo "Error: FASTQ file $FASTQ not found"
- exit 1
- fi
- # Get base name for output files
- BASE=$(basename $FASTQ _R1_001.fastq.gz)
- echo "Processing $BASE..."
- # Create a directory to store trimmed adapter files if it doesn't exist
- ADAPTER_DIR="$WORKDIR/adapters"
- mkdir -p $ADAPTER_DIR
- ADAPTER_FILE="$ADAPTER_DIR/TruSeq3-SE-2.fa"
- # Create adapter file if it doesn't exist
- if [ ! -f "$ADAPTER_FILE" ]; then
- echo "Creating TruSeq3-SE-2.fa adapter file..."
- cat > "$ADAPTER_FILE" << 'EOL'
- >TruSeq3_IndexedAdapter
- AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
- >TruSeq3_UniversalAdapter
- AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGTA
- EOL
- fi
- # Step 1: Trim reads using Trimmomatic as specified in the paper
- TRIMMED_FASTQ="$OUTPUT_DIR/trimmed/${BASE}.trimmed.fastq.gz"
- if [ -f "$TRIMMED_FASTQ" ]; then
- echo "Trimmed file already exists, skipping trimming step..."
- else
- echo "Trimming reads with Trimmomatic..."
- trimmomatic SE $FASTQ $TRIMMED_FASTQ \
- ILLUMINACLIP:$ADAPTER_FILE:2:30:10 \
- SLIDINGWINDOW:4:15 \
- MINLEN:36
- fi
- # Step 2: Align with BWA as specified in the paper
- echo "Aligning with BWA..."
- # Check if BWA index exists and is valid, if not create it
- 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")
- INDEX_VALID=true
- # Check if all index files exist and have non-zero size
- for index_file in "${BWA_INDEX_FILES[@]}"; do
- if [ ! -s "$index_file" ]; then
- INDEX_VALID=false
- echo "BWA index file $index_file is missing or empty."
- break
- fi
- done
- # Even if all files exist, perform a more thorough validation by attempting to use the index
- if [ "$INDEX_VALID" = true ]; then
- echo "Testing BWA index validity..."
- # Create a small test FASTQ with a few reads
- TEST_FASTQ="${BWA_GENOME_DIR}/test_reads.fastq"
- echo -e "@test_read\nACGTACGTACGTACGTACGT\n+\nIIIIIIIIIIIIIIIIIIII" > "$TEST_FASTQ"
- # Try to align the test read - redirect stderr to capture any errors
- TEST_OUTPUT=$(bwa mem "$BWA_GENOME_DIR/mm10" "$TEST_FASTQ" 2>&1 || true)
- # Check if the output contains error messages related to index files
- if echo "$TEST_OUTPUT" | grep -q -E "(Can't|Cannot|failed|error|Unexpected end of file)"; then
- INDEX_VALID=false
- echo "BWA index validation failed with error:"
- echo "$TEST_OUTPUT" | grep -E "(Can't|Cannot|failed|error|Unexpected end of file)"
- else
- echo "BWA index validation successful."
- fi
- # Clean up test file
- rm -f "$TEST_FASTQ"
- fi
- if [ "$INDEX_VALID" = false ]; then
- echo "BWA index is invalid or incomplete. Removing old index files and recreating..."
- # Remove any existing index files that might be corrupted
- rm -f ${BWA_INDEX_FILES[@]}
- # Also remove any potential temporary files that might have been created during a failed indexing
- rm -f "$BWA_GENOME_DIR/mm10."*
- mkdir -p $BWA_GENOME_DIR
- # Check if reference genome exists
- if [ ! -f "$WORKDIR/references/mm10.fa" ]; then
- echo "Error: Reference genome file $WORKDIR/references/mm10.fa not found"
- exit 1
- fi
- # Check available disk space before indexing
- GENOME_SIZE=$(du -b "$WORKDIR/references/mm10.fa" | cut -f1)
- REQUIRED_SPACE=$((GENOME_SIZE * 5)) # BWA index typically requires ~5x the genome size
- AVAILABLE_SPACE=$(df -B1 "$BWA_GENOME_DIR" | awk 'NR==2 {print $4}')
- if [ "$AVAILABLE_SPACE" -lt "$REQUIRED_SPACE" ]; then
- 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."
- exit 1
- fi
- echo "Creating BWA index. This may take a while..."
- bwa index -p $BWA_GENOME_DIR/mm10 "$WORKDIR/references/mm10.fa"
- # Verify index was created successfully
- for index_file in "${BWA_INDEX_FILES[@]}"; do
- if [ ! -s "$index_file" ]; then
- echo "Error: Failed to create BWA index file $index_file"
- exit 1
- fi
- done
- else
- echo "BWA index is valid and complete."
- fi
- # Check if BAM file already exists
- if [ -f "$OUTPUT_DIR/bam/${BASE}.bam" ]; then
- echo "BAM file already exists, skipping alignment step..."
- else
- # Align with BWA using parameters from the paper
- echo "Performing BWA alignment..."
- bwa mem -t $THREADS $BWA_GENOME_DIR/mm10 $OUTPUT_DIR/trimmed/${BASE}.trimmed.fastq.gz > $OUTPUT_DIR/bam/${BASE}.sam
- # Convert SAM to BAM
- echo "Converting SAM to BAM..."
- samtools view -bS $OUTPUT_DIR/bam/${BASE}.sam > $OUTPUT_DIR/bam/${BASE}.bam
- rm $OUTPUT_DIR/bam/${BASE}.sam
- fi
- # Sort BAM if it exists but hasn't been sorted yet
- if [ -f "$OUTPUT_DIR/bam/${BASE}.bam" ]; then
- # Check if BAM is already sorted
- if ! samtools view -H "$OUTPUT_DIR/bam/${BASE}.bam" | grep -q 'SO:coordinate'; then
- echo "Sorting BAM..."
- samtools sort -@ $THREADS $OUTPUT_DIR/bam/${BASE}.bam -o $OUTPUT_DIR/bam/${BASE}.sorted.bam
- mv $OUTPUT_DIR/bam/${BASE}.sorted.bam $OUTPUT_DIR/bam/${BASE}.bam
- else
- echo "BAM file already sorted, skipping sorting step..."
- fi
- fi
- # Step 3: Mark and remove PCR duplicates using Picard
- if [ -f "$OUTPUT_DIR/bam/${BASE}.dedup.bam" ]; then
- echo "Deduplicated BAM already exists, skipping deduplication step..."
- else
- echo "Marking and removing PCR duplicates with Picard..."
- picard MarkDuplicates \
- I=$OUTPUT_DIR/bam/${BASE}.bam \
- O=$OUTPUT_DIR/bam/${BASE}.dedup.bam \
- M=$OUTPUT_DIR/bam/${BASE}.metrics.txt \
- REMOVE_DUPLICATES=true \
- TMP_DIR=$OUTPUT_DIR/bam/temp
- fi
- # Step 4: Filter by mapping quality (q>=1) as in the paper
- if [ -f "$OUTPUT_DIR/bam/${BASE}.filtered.bam" ]; then
- echo "Filtered BAM already exists, skipping filtering step..."
- else
- echo "Filtering by mapping quality..."
- samtools view -q 1 -b $OUTPUT_DIR/bam/${BASE}.dedup.bam > $OUTPUT_DIR/bam/${BASE}.filtered.bam
- fi
- # Index BAM if needed
- if [ ! -f "$OUTPUT_DIR/bam/${BASE}.filtered.bam.bai" ]; then
- echo "Indexing BAM..."
- samtools index $OUTPUT_DIR/bam/${BASE}.filtered.bam
- else
- echo "BAM index already exists, skipping indexing step..."
- fi
- # Step 5: Create bigWig file
- if [ -f "$OUTPUT_DIR/bigwig/${BASE}.bw" ]; then
- echo "BigWig file already exists, skipping bigWig creation step..."
- else
- echo "Creating bigWig for $BASE..."
- bamCoverage \
- --bam $OUTPUT_DIR/bam/${BASE}.filtered.bam \
- --outFileName $OUTPUT_DIR/bigwig/${BASE}.bw \
- --binSize 50 \
- --normalizeUsing RPKM \
- --effectiveGenomeSize 2652783500 \
- --blackListFileName $BLACKLIST_FILE \
- --extendReads 250 \
- --numberOfProcessors $THREADS
- fi
- echo "Processing complete for $BASE"
01_process_fastq.sh at commit a378313, no license · at the source
Overview
- Institute of Neuroscience, National Research Council (CNR),Milan, Italy
- Stem Cell and Neurogenesis Unit, Division of Neuroscience, IRCCS San Raffaele Scientific Institute,Milan, Italy
- Neuroepigenetics Unit, Division of Neuroscience, IRCCS San Raffaele Scientific Institute,Milan, Italy
- INGM Istituto Nazionale Genetica Molecolare Romeo ed Enrica Invernizzi,Milan, Italy
- Institute of Biomedical Technologies, National Research Council (CNR),Milan, Italy
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
Availability: 1 check, the latest on 29 September 2026: the link answers (HTTP 200)
- 29 September 2026: the link answers (HTTP 200)
51 files
- 01_CUTandTAG_processing/
SEACR_1.3.R , R, 153 lines - 01_CUTandTAG_processing/
SEACR_1.3.sh , Shell, 199 lines - 01_CUTandTAG_processing/
annotate_peaks.R , R, 40 lines - 01_CUTandTAG_processing/
calculate_frip.R , R, 21 lines - 01_CUTandTAG_processing/
calculate_tss_enrichment , R, 37 lines.R - 01_CUTandTAG_processing/
combine_peaks_broad.py , Python, 158 lines - 01_CUTandTAG_processing/
combine_peaks_narrow.py , Python, 119 lines - 01_CUTandTAG_processing/
diffbind_analysis.R , R, 47 lines - 01_CUTandTAG_processing/
filter_peaks.R , R, 22 lines - 02_CpG_enrichment/
cpg_enrichment_analysis. , Python, 364 linespy - 02_CpG_enrichment/
cpg_enrichment_in_island , Python, 363 liness.py - 02_CpG_enrichment/
cpg_gene_analysis.R , R, 246 lines - 02_CpG_enrichment/
integrate_rnaseq.R , R, 113 lines - 02_CpG_enrichment/
peak_localization.R , R, 308 lines - 03_RNAseq_analysis/
01_fastqc.sh , Shell, 57 lines - 03_RNAseq_analysis/
02_star_alignment.sh , Shell, 99 lines - 03_RNAseq_analysis/
03_featurecounts.sh , Shell, 77 lines - 03_RNAseq_analysis/
04_deseq2_analysis.R , R, 215 lines - 03_RNAseq_analysis/
05_functional_enrichment , R, 256 lines.R - 03_RNAseq_analysis/
06_visualizations.R , R, 632 lines - 03_RNAseq_analysis/
generate_de_table.R , R, 249 lines - 03_RNAseq_analysis/
generate_volcano_plots.R , R, 310 lines - 04_binding_expression_in
tegration/ , Python, 644 linesbinding_heatmap.py - 04_binding_expression_in
tegration/ , Python, 798 linesgene_body_analysis.py - 04_binding_expression_in
tegration/ , Python, 585 linesintegrate_genebody_rnase q.py - 04_binding_expression_in
tegration/ , Python, 325 linesintegrate_rnaseq_base.py - 04_binding_expression_in
tegration/ , Python, 765 linespeak_expression_correlat ion.py - 04_binding_expression_in
tegration/ , Python, 641 linespromoter_expression_corr elation.py - 05_Smarcb1_analysis/
01_cpg_metaprofile.py , Python, 511 lines - 05_Smarcb1_analysis/
02_cpg_targets.py , Python, 132 lines - 05_Smarcb1_analysis/
03_smarcb1_comparison.py , Python, 151 lines - 05_Smarcb1_analysis/
04_smarcb1_comparison_ex , Python, 176 linestended.py - 05_Smarcb1_analysis/
05_smarcb1_dea_analysis. , Python, 166 linespy - 06_SAMMY_seq/
00_build_index.sh , Shell, 39 lines - 06_SAMMY_seq/
01_process_fastq.sh , Shell, 275 lines - 06_SAMMY_seq/
02_run_sammy_analysis.sh , Shell, 33 lines - 06_SAMMY_seq/
02_sammy_analysis.py , Python, 674 lines - 06_SAMMY_seq/
02b_run_sammy_blackliste , Shell, 35 linesd.sh - 06_SAMMY_seq/
02b_sammy_analysis_black , Python, 712 lineslisted.py - 06_SAMMY_seq/
03_chromatin_state_analy , Python, 761 linessis.py - 06_SAMMY_seq/
03_run_chromatin_analysi , Shell, 39 liness.sh - 06_SAMMY_seq/
04_chromatin_comparison. , Python, 713 linespy - 06_SAMMY_seq/
04_run_comparison_neu.sh , Shell, 39 lines - 06_SAMMY_seq/
04_run_comparison_nsc.sh , Shell, 39 lines - 06_SAMMY_seq/
05_metaprofile_heatmap.p , Python, 391 linesy - 06_SAMMY_seq/
05_run_metaprofile_heatm , Shell, 49 linesap.sh - 06_SAMMY_seq/
06_associate_regions_gen , Python, 280 lineses.py - 06_SAMMY_seq/
07_target_comparison_hea , Python, 735 linestmap.py - 06_SAMMY_seq/
sammy_analysis_extended. , Python, 880 linespy - 07_visualization/
01_metaprofiles.sh , Shell, 124 lines - README.md, Text, 1 line
u9856112956-prog/mecp2_analysis
a378313b060fefaa5e141b10bffedc0a6e2f36bd, 15 January 2026Availability: 1 check, the latest on 29 September 2026: the link answers
- 29 September 2026: the link answers
51 files
- 01_CUTandTAG_processing/
SEACR_1.3.R , R, 153 lines - 01_CUTandTAG_processing/
SEACR_1.3.sh , Shell, 199 lines - 01_CUTandTAG_processing/
annotate_peaks.R , R, 40 lines - 01_CUTandTAG_processing/
calculate_frip.R , R, 21 lines - 01_CUTandTAG_processing/
calculate_tss_enrichment , R, 37 lines.R - 01_CUTandTAG_processing/
combine_peaks_broad.py , Python, 158 lines - 01_CUTandTAG_processing/
combine_peaks_narrow.py , Python, 119 lines - 01_CUTandTAG_processing/
diffbind_analysis.R , R, 47 lines - 01_CUTandTAG_processing/
filter_peaks.R , R, 22 lines - 02_CpG_enrichment/
cpg_enrichment_analysis. , Python, 364 lines, 1 matchpy - 02_CpG_enrichment/
cpg_enrichment_in_island , Python, 363 liness.py - 02_CpG_enrichment/
cpg_gene_analysis.R , R, 246 lines - 02_CpG_enrichment/
integrate_rnaseq.R , R, 113 lines, 1 match - 02_CpG_enrichment/
peak_localization.R , R, 308 lines - 03_RNAseq_analysis/
01_fastqc.sh , Shell, 57 lines - 03_RNAseq_analysis/
02_star_alignment.sh , Shell, 99 lines, 1 match - 03_RNAseq_analysis/
03_featurecounts.sh , Shell, 77 lines - 03_RNAseq_analysis/
04_deseq2_analysis.R , R, 215 lines, 1 match - 03_RNAseq_analysis/
05_functional_enrichment , R, 256 lines, 2 matches.R - 03_RNAseq_analysis/
06_visualizations.R , R, 632 lines - 03_RNAseq_analysis/
generate_de_table.R , R, 249 lines, 1 match - 03_RNAseq_analysis/
generate_volcano_plots.R , R, 310 lines - 04_binding_expression_in
tegration/ , Python, 644 linesbinding_heatmap.py - 04_binding_expression_in
tegration/ , Python, 798 linesgene_body_analysis.py - 04_binding_expression_in
tegration/ , Python, 585 lines, 1 matchintegrate_genebody_rnase q.py - 04_binding_expression_in
tegration/ , Python, 325 lines, 2 matchesintegrate_rnaseq_base.py - 04_binding_expression_in
tegration/ , Python, 765 lines, 1 matchpeak_expression_correlat ion.py - 04_binding_expression_in
tegration/ , Python, 641 linespromoter_expression_corr elation.py - 05_Smarcb1_analysis/
01_cpg_metaprofile.py , Python, 511 lines - 05_Smarcb1_analysis/
02_cpg_targets.py , Python, 132 lines - 05_Smarcb1_analysis/
03_smarcb1_comparison.py , Python, 151 lines - 05_Smarcb1_analysis/
04_smarcb1_comparison_ex , Python, 176 linestended.py - 05_Smarcb1_analysis/
05_smarcb1_dea_analysis. , Python, 166 linespy - 06_SAMMY_seq/
00_build_index.sh , Shell, 39 lines - 06_SAMMY_seq/
01_process_fastq.sh , Shell, 275 lines, 2 matches - 06_SAMMY_seq/
02_run_sammy_analysis.sh , Shell, 33 lines, 2 matches - 06_SAMMY_seq/
02_sammy_analysis.py , Python, 674 lines - 06_SAMMY_seq/
02b_run_sammy_blackliste , Shell, 35 linesd.sh - 06_SAMMY_seq/
02b_sammy_analysis_black , Python, 712 lines, 1 matchlisted.py - 06_SAMMY_seq/
03_chromatin_state_analy , Python, 761 linessis.py - 06_SAMMY_seq/
03_run_chromatin_analysi , Shell, 39 liness.sh - 06_SAMMY_seq/
04_chromatin_comparison. , Python, 713 linespy - 06_SAMMY_seq/
04_run_comparison_neu.sh , Shell, 39 lines - 06_SAMMY_seq/
04_run_comparison_nsc.sh , Shell, 39 lines - 06_SAMMY_seq/
05_metaprofile_heatmap.p , Python, 391 linesy - 06_SAMMY_seq/
05_run_metaprofile_heatm , Shell, 49 linesap.sh - 06_SAMMY_seq/
06_associate_regions_gen , Python, 280 lineses.py - 06_SAMMY_seq/
07_target_comparison_hea , Python, 735 lines, 1 matchtmap.py - 06_SAMMY_seq/
sammy_analysis_extended. , Python, 880 linespy - 07_visualization/
01_metaprofiles.sh , Shell, 124 lines, 1 match - README.md, Text, 1 line
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
- geo:GSE299957, at NCBI GEO; found in “Data availability”
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:
- it points to a dataset: NCBI GEO GSE299957
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://
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/
url = {https://
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/
VL - 17
IS - 1
SP - 3225
SN - 2041-1723
PB - Nature Publishing Group
DO - 10.1038/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1038/
"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":
"volume": "17",
"issue": "1",
"page": "3225",
"DOI": "10.1038/
"PMID": "42026056",
"PMCID": "PMC13106735",
"ISSN": "2041-1723",
"publisher": "Nature Publishing Group",
"URL": "https://
"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: NatureIn 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 genomeJournal: 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 advancesIn 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: NatureIn 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 communicationsIn 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: eLifeIn 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 diseaseIn 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 advancesIn 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 biologyIn 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.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 2 repositories of the authors' code, each at its verified commit and with its license, 100 scripts, and 18 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:e603aa1df2c836ca…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[, paste the snippet at the top, then “Commit changes…” and, to review it first, “Create a new branch and start a pull request”. You open the pull request; OSCR asks for no permission.
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.
