Person:

Collins, Ryan

Background: Structural 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. Results: We 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. Conclusions: These 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 material The online version of this article (doi:10.1186/s13059-017-1158-6) contains supplementary material, which is available to authorized users.

Recurrent, reciprocal genomic disorders resulting from non-allelic homologous recombination (NAHR) between near-identical segmental duplications (SDs) are a major cause of human disease, often producing phenotypically distinct syndromes. The genomic architecture of flanking SDs presents a significant challenge for modeling these syndromes; however, the capability to efficiently generate reciprocal copy number variants (CNVs) that mimic NAHR would represent an invaluable modeling tool. We describe here a CRISPR/Cas9 genome engineering method, Single-guide-CRISPR/Cas-targeting-Of-Repetitive-Elements (SCORE), to model reciprocal genomic disorders and demonstrate its capabilities by generating reciprocal CNVs of 16p11.2 and 15q13.3, including alteration of one copy-equivalent of the SDs that mediate NAHR in vivo. The method is reproducible and RNAseq reliably clusters transcriptional signatures from human subjects with in vivo CNV and their corresponding in vitro models. This new approach will provide broad applicability for the study of genomic disorders and, with further development, may also permit efficient correction of these defects.

Abstract The BAM and CRAM formats provide a supplementary linear index that facilitates rapid access to sequence alignments in arbitrary genomic regions. Comparing consecutive entries in a BAM or CRAM index allows one to infer the number of alignment records per genomic region for use as an effective proxy of sequence depth in each genomic region. Based on these properties, we have developed indexcov, an efficient estimator of whole-genome sequencing coverage to rapidly identify samples with aberrant coverage profiles, reveal large-scale chromosomal anomalies, recognize potential batch effects, and infer the sex of a sample. Indexcov is available at https://github.com/brentp/goleft under the MIT license.