DOI: 10.25452/figshare.plus.26517577
title: Single cell ATACseq of mouse developing cranial motor neurons, spinal motor neurons, and surrounding neuronal tissue at e10.5 and e11.5.

description: Although Mendelian disorders are overwhelmingly attributed to protein-coding pathogenic variants, a majority of unsolved cases do not harbor obvious causal pathogenic variants in the coding sequence, suggesting a potential non-coding etiology. However, classification of pathogenicity in non-coding sequence remains prohibitive due to a vastly increased search space and the lack of a standardized rubric for interpretation. Here, we present an integrated single cell multiomic framework to nominate pathogenic non-coding variants for the congenital cranial dysinnervation disorders (CCDDs). The CCDDs are Mendelian neurodevelopmental disorders that result from aberrant development of cranial motor neurons in the embryonic brainstem. We created a non-coding reference atlas of single cell chromatin accessibility profiles for 86,089 embryonic mouse cranial motor neurons (cMNs). We found that high-quality single cell ATAC-seq (scATAC) profiles alone were a strong predictor of enhancement (64% in vivo validation rate). To further aid in interpretation, we integrated single cell histone modification and gene expression information to distinguish individual enhancers and their cognate genes. Relatively subtle differences in cellular composition of input data often led to substantial differences in predicted enhancer strength, cognate gene, and tissue of activity. Next, we mapped candidate non-coding variants from 899 whole genome sequences from 270 CCDD pedigrees to the murine cMN-specific regulatory elements and trained a machine learning classifier to accurately predict the functional effects of patient variants within these elements. We then performed high coverage scATACseq and site-specific footprinting analysis on an allelic series of CRISPR-humanised mice to validate our machine learning predictions and render important clues to the mode of pathogenicity. Finally, we performed peak- and gene-centric allelic aggregation to nominate non-coding variants, including those regulating MN1 and EBF3, respectively. Altogether this work extends non-coding variant analysis to Mendelian disease and presents a generalizable framework for nominating novel non-coding variants in other rare disorders.Fluorescence-assisted microdissection was performed to collect samples cMN3/4, cMN7, and sMN from <i>Isl1</i><sup><em>MN</em></sup>:GFP mice and likewise to collect samples of cMN6, cMN12, and sMN from <i>Hb9</i>:GFP mice, each at both e10.5 and e11.5. Nuclei were isolated in accordance with Low Cell Input Nuclei Isolation guidelines provided by ‘Demonstrated Protocol – Nuclei Isolation for Single Cell ATAC Sequencing Rev A (Protocol #CG000169) from 10x Genomics. performed scATAC transposition, droplet formation, and library construction as described in protocol CG000168 using v1 reagents (10x Genomics). scATAC libraries were sequenced on the Illumina NextSeq 500 system using standard Illumina chemistry. Paired inserts were minimum 2 x 34 bp in length excluding indices, and libraries were distributed to achieve an estimated coverage of ≥ 25,000 read pairs per cell in accordance with 10x Genomics guidelines. bcl outputs from NextSeq 500 were processed using 10x Cell Ranger ATAC mkfastq pipeline (https://www.10xgenomics.com/support/software/cell-ranger/) for demultiplexing and fastq file generation. ATAC fragment counts were tabulated using the Cell Ranger count function using default parameters. .bam and/or fragment files were processed for peak calling and differential accessibility analysis using R ArchR package (https://github.com/GreenleafLab/ArchR).