SUMMARY Balanced chromosomal abnormalities (BCAs) represent a reservoir of single gene disruptions in neurodevelopmental disorders (NDD). We sequenced BCAs in autism and related NDDs, revealing disruption of 33 loci in four general categories: 1) genes associated with abnormal neurodevelopment (e.g., AUTS2, FOXP1, CDKL5), 2) single gene contributors to microdeletion syndromes (MBD5, SATB2, EHMT1, SNURF-SNRPN), 3) novel risk loci (e.g., CHD8, KIRREL3, ZNF507), and 4) genes associated with later onset psychiatric disorders (e.g., TCF4, ZNF804A, PDE10A, GRIN2B, ANK3). We also discovered profoundly increased burden of copy number variants among 19,556 neurodevelopmental cases compared to 13,991 controls (p = 2.07×10−47) and enrichment of polygenic risk alleles from autism and schizophrenia genome-wide association studies (p = 0.0018 and 0.0009, respectively). Our findings suggest a polygenic risk model of autism incorporating loci of strong effect and indicate that some neurodevelopmental genes are sensitive to perturbation by multiple mutational mechanisms, leading to variable phenotypic outcomes that manifest at different life stages.

Despite their clinical significance, characterization of balanced chromosomal abnormalities (BCAs) has largely been restricted to cytogenetic resolution. We explored the landscape of BCAs at nucleotide resolution in 273 subjects with a spectrum of congenital anomalies. Whole-genome sequencing revised 93% of karyotypes and revealed complexity that was cryptic to karyotyping in 21% of BCAs, highlighting the limitations of conventional cytogenetic approaches. At least 33.9% of BCAs resulted in gene disruption that likely contributed to the developmental phenotype, 5.2% were associated with pathogenic genomic imbalances, and 7.3% disrupted topologically associated domains (TADs) encompassing known syndromic loci. Remarkably, BCA breakpoints in eight subjects altered a single TAD encompassing MEF2C, a known driver of 5q14.3 microdeletion syndrome, resulting in decreased MEF2C expression. This study proposes that sequence-level resolution dramatically improves prediction of clinical outcomes for balanced rearrangements, and provides insight into novel pathogenic mechanisms such as altered regulation due to changes in chromosome topology.

We defined the genetic landscape of balanced chromosomal rearrangements at nucleotide resolution by sequencing 141 breakpoints from cytogenetically-interpreted translocations and inversions. We confirm that the recently described phenomenon of “chromothripsis” (massive chromosomal shattering and reorganization) is not unique to cancer cells but also occurs in the germline where it can resolve to a karyotypically balanced state with frequent inversions. We detected a high incidence of complex rearrangements (19.2%) and substantially less reliance on microhomology (31%) than previously observed in benign CNVs. We compared these results to experimentally-generated DNA breakage-repair by sequencing seven transgenic animals, and revealed extensive rearrangement of the transgene and host genome with similar complexity to human germline alterations. Inversion is the most common rearrangement, suggesting that a combined mechanism involving template switching and non-homologous repair mediates the formation of balanced complex rearrangements that are viable, stably replicated and transmitted unaltered to subsequent generations.

Persons with neurodevelopmental disorders or autism spectrum disorder (ASD) often harbor chromosomal microdeletions, yet the individual genetic contributors within these regions have not been systematically evaluated. We established a consortium of clinical diagnostic and research laboratories to accumulate a large cohort with genetic alterations of chromosomal region 2q23.1 and acquired 65 subjects with microdeletion or translocation. We sequenced translocation breakpoints; aligned microdeletions to determine the critical region; assessed effects on mRNA expression; and examined medical records, photos, and clinical evaluations. We identified a single gene, methyl-CpG-binding domain 5 (MBD5), as the only locus that defined the critical region. Partial or complete deletion of MBD5 was associated with haploinsufficiency of mRNA expression, intellectual disability, epilepsy, and autistic features. Fourteen alterations, including partial deletions of noncoding regions not typically captured or considered pathogenic by current diagnostic screening, disrupted MBD5 alone. Expression profiles and clinical characteristics were largely indistinguishable between MBD5-specific alteration and deletion of the entire 2q23.1 interval. No copy-number alterations of MBD5 were observed in 7878 controls, suggesting MBD5 alterations are highly penetrant. We surveyed MBD5 coding variations among 747 ASD subjects compared to 2043 non-ASD subjects analyzed by whole-exome sequencing and detected an association with a highly conserved methyl-CpG-binding domain missense variant, p.79Gly>Glu (c.236G>A) (p = 0.012). These results suggest that genetic alterations of MBD5 cause features of 2q23.1 microdeletion syndrome and that this epigenetic regulator significantly contributes to ASD risk, warranting further consideration in research and clinical diagnostic screening and highlighting the importance of chromatin remodeling in the etiology of these complex disorders.

BackgroundStructural variation (SV) influences genome organization and contributes to human disease. However, the complete mutational spectrum of SV has not been routinely captured in disease association studies.ResultsWe sequenced 689 participants with autism spectrum disorder (ASD) and other developmental abnormalities to construct a genome-wide map of large SV. Using long-insert jumping libraries at 105X mean physical coverage and linked-read whole-genome sequencing from 10X Genomics, we document seven major SV classes at ~5 kb SV resolution. Our results encompass 11,735 distinct large SV sites, 38.1% of which are novel and 16.8% of which are balanced or complex. We characterize 16 recurrent subclasses of complex SV (cxSV), revealing that: (1) cxSV are larger and rarer than canonical SV; (2) each genome harbors 14 large cxSV on average; (3) 84.4% of large cxSVs involve inversion; and (4) most large cxSV (93.8%) have not been delineated in previous studies. Rare SVs are more likely to disrupt coding and regulatory non-coding loci, particularly when truncating constrained and disease-associated genes. We also identify multiple cases of catastrophic chromosomal rearrangements known as chromoanagenesis, including somatic chromoanasynthesis, and extreme balanced germline chromothripsis events involving up to 65 breakpoints and 60.6 Mb across four chromosomes, further defining rare categories of extreme cxSV.ConclusionsThese data provide a foundational map of large SV in the morbid human genome and demonstrate a previously underappreciated abundance and diversity of cxSV that should be considered in genomic studies of human disease.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-017-1158-6) contains supplementary material, which is available to authorized users.

The contribution of balanced chromosomal rearrangements to complex disorders remains unclear because they are not detected routinely by genome-wide microarrays and clinical localization is imprecise. Failure to consider these events bypasses a potentially powerful complement to single nucleotide polymorphism and copy-number association approaches to complex disorders, where much of the heritability remains unexplained. To capitalize on this genetic resource, we have applied optimized sequencing and analysis strategies to test whether these potentially high-impact variants can be mapped at reasonable cost and throughput. By using a whole-genome multiplexing strategy, rearrangement breakpoints could be delineated at a fraction of the cost of standard sequencing. For rearrangements already mapped regionally by karyotyping and fluorescence in situ hybridization, a targeted approach enabled capture and sequencing of multiple breakpoints simultaneously. Importantly, this strategy permitted capture and unique alignment of up to 97% of repeat-masked sequences in the targeted regions. Genome-wide analyses estimate that only 3.7% of bases should be routinely omitted from genomic DNA capture experiments. Illustrating the power of these approaches, the rearrangement breakpoints were rapidly defined to base pair resolution and revealed unexpected sequence complexity, such as co-occurrence of inversion and translocation as an underlying feature of karyotypically balanced alterations. These findings have implications ranging from genome annotation to de novo assemblies and could enable sequencing screens for structural variations at a cost comparable to that of microarrays in standard clinical practice.

X-linked Dystonia-Parkinsonism (XDP) is a Mendelian neurodegenerative disease that is endemic to the Philippines and is associated with a founder haplotype. We integrated multiple genome and transcriptome assembly technologies to narrow the causal mutation to the TAF1 locus, which included a SINE-VNTR-Alu (SVA) retrotransposition into intron 32 of the gene. Transcriptome analyses identified decreased expression of the canonical cTAF1 transcript among XDP probands, and de novo assembly across multiple pluripotent stem-cell-derived neuronal lineages discovered aberrant TAF1 transcription that involved alternative splicing and intron retention (IR) in proximity to the SVA that was anti-correlated with overall TAF1 expression. CRISPR/Cas9 excision of the SVA rescued this XDP-specific transcriptional signature and normalized TAF1 expression in probands. These data suggest an SVA-mediated aberrant transcriptional mechanism associated with XDP and may provide a roadmap for layered technologies and integrated assembly-based analyses for other unsolved Mendelian disorders.

SUMMARY Conventional cytogenetic testing offers low-resolution detection of balanced karyotypic abnormalities but cannot provide the precise, gene-level knowledge required to predict outcomes. The use of high-resolution whole-genome deep sequencing is currently impractical for the purpose of routine clinical care. We show here that whole-genome “jumping libraries” can offer an immediately applicable, nucleotide-level complement to conventional genetic diagnostics within a time frame that allows for clinical action. We performed large-insert sequencing of DNA extracted from amniotic-fluid cells with a balanced de novo translocation. The amniotic-fluid sample was from a patient in the third trimester of pregnancy who underwent amniocentesis because of severe polyhydramnios after multiple fetal anomalies had been detected on ultrasonography. Using a 13-day sequence and analysis pipeline, we discovered direct disruption of CHD7, a causal locus in the CHARGE syndrome (coloboma of the eye, heart anomaly, atresia of the choanae, retardation, and genital and ear anomalies). Clinical findings at birth were consistent with the CHARGE syndrome, a diagnosis that could not have been reliably inferred from the cytogenetic breakpoint. This case study illustrates the potential power of customized whole-genome jumping libraries when used to augment prenatal karyotyping.