Whole-Genome Sequencing of Individuals from a Founder Population Identifies Candidate Genes for Asthma

Asthma is a complex genetic disease caused by a combination of genetic and environmental risk factors. We sought to test classes of genetic variants largely missed by genome-wide association studies (GWAS), including copy number variants (CNVs) and low-frequency variants, by performing whole-genome sequencing (WGS) on 16 individuals from asthma-enriched and asthma-depleted families. The samples were obtained from an extended 13-generation Hutterite pedigree with reduced genetic heterogeneity due to a small founding gene pool and reduced environmental heterogeneity as a result of a communal lifestyle. We sequenced each individual to an average depth of 13-fold, generated a comprehensive catalog of genetic variants, and tested the most severe mutations for association with asthma. We identified and validated 1960 CNVs, 19 nonsense or splice-site single nucleotide variants (SNVs), and 18 insertions or deletions that were out of frame. As follow-up, we performed targeted sequencing of 16 genes in 837 cases and 540 controls of Puerto Rican ancestry and found that controls carry a significantly higher burden of mutations in IL27RA (2.0% of controls; 0.23% of cases; nominal p = 0.004; Bonferroni p = 0.21). We also genotyped 593 CNVs in 1199 Hutterite individuals. We identified a nominally significant association (p = 0.03; Odds ratio (OR) = 3.13) between a 6 kbp deletion in an intron of NEDD4L and increased risk of asthma. We genotyped this deletion in an additional 4787 non-Hutterite individuals (nominal p = 0.056; OR = 1.69). NEDD4L is expressed in bronchial epithelial cells, and conditional knockout of this gene in the lung in mice leads to severe inflammation and mucus accumulation. Our study represents one of the early instances of applying WGS to complex disease with a large environmental component and demonstrates how WGS can identify risk variants, including CNVs and low-frequency variants, largely untested in GWAS

Mechanistic dissection of chromatin topology disruption in the 5q14.3 MEF2C locus as an indirect driver of neurodevelopmental disorders

Structural variants have the potential to create long-range positional effects, uncouple genes from regulatory elements, and facilitate aberrant 3D chromatin folding. In an independent study, we revealed a significant enrichment of intergenic balanced chromosomal abnormality (BCA) breakpoints from congenital anomaly cases at chromosome 5q14.3. All 11 5q14.3 BCA carriers had breakpoints disrupting the topologically associated domain (TAD) housing MEF2C, a known driver of neurodevelopmental disorders. In a second study, we showed, using 4C-seq, that MEF2C forms proximal and distal interaction loops within this TAD, contacting multiple neuronal enhancers. To understand the functional impact of these BCAs, we performed a mechanistic dissection of the 3D regulatory network at the 5q14.3 locus and its constituent functional elements using Cas9-based genome editing. We generated >200 cell lines, representing deletions of MEF2C alongside five MEF2C TAD boundary and interaction loop targets in iPS-derived neural stem cells (NSCs) and cortical induced neurons (iNs). We profiled changes in expression and 3D chromatin interactions within these models using RNA-seq and 4C-seq. MEF2C was variably differentially expressed upon deletion of proximal and distal MEF2C loop regions, with effects depending on cell type and variant class. Underlying chromatin interaction patterns revealed evidence of loop maintenance in these models, possibly via CTCF buffering, highlighting compensatory mechanisms against 3D chromatin disruption. In contrast, deletion of the proximal TAD boundary facilitated increased contacts with predicted enhancers in the adjacent TAD. Our results suggest novel regulatory mechanisms driving phenotypic outcomes for the 5q14.3 region, with significant implications for interpretation of pathogenic structural variation

Biallelic mutation of FBXL7 suggests a novel form of Hennekam syndrome

Hennekam lymphangiectasia-lymphedema syndrome is an autosomal recessive disorder characterized by congenital lymphedema, intestinal lymphangiectasia, facial dysmorphism, and variable intellectual disability. Known disease genes include CCBE1, FAT4, and ADAMTS3. In a patient with clinically diagnosed Hennekam syndrome but without mutations or copy-number changes in the three known disease genes, we identified a homozygous single-exon deletion affecting FBXL7. Specifically, exon 3, which encodes the F-box domain and several leucine-rich repeats of FBXL7, is eliminated. Our analyses of databases representing >100,000 control individuals failed to identify biallelic loss-of-function variants in FBXL7. Published studies in Drosophila indicate Fbxl7 interacts with Fat, of which human FAT4 is an ortholog, and mutation of either gene yields similar morphological consequences. These data suggest that FBXL7 may be the fourth gene for Hennekam syndrome, acting via a shared pathway with FAT4

Interchromosomal core duplicons drive both evolutionary instability and disease susceptibility of the Chromosome 8p23.1 region

Recurrent rearrangements of Chromosome 8p23.1 are associated with congenital heart defects and developmental delay. The complexity of this region has led to inconsistencies in the current reference assembly, confounding studies of genetic variation. Using comparative sequence-based approaches, we generated a high-quality 6.3-Mbp alternate reference assembly of an inverted Chromosome 8p23.1 haplotype. Comparison with nonhuman primates reveals a 746-kbp duplicative transposition and two separate inversion events that arose in the last million years of human evolution. The breakpoints associated with these rearrangements map to an ape-specific interchromosomal core duplicon that clusters at sites of evolutionary inversion (P = 7.8 × 10(-5)). Refinement of microdeletion breakpoints identifies a subgroup of patients that map to the same interchromosomal core involved in the evolutionary formation of the duplication blocks. Our results define a higher-order genomic instability element that has shaped the structure of specific chromosomes during primate evolution contributing to rearrangements associated with inversion and disease

Evolution and diversity of copy number variation in the great ape lineage

Copy number variation (CNV) contributes to disease and has restructured the genomes of great apes. The diversity and rate of this process, however, have not been extensively explored among great ape lineages. We analyzed 97 deeply sequenced great ape and human genomes and estimate 16% (469 Mb) of the hominid genome has been affected by recent CNV. We identify a comprehensive set of fixed gene deletions (n = 340) and duplications (n = 405) as well as >13.5 Mb of sequence that has been specifically lost on the human lineage. We compared the diversity and rates of copy number and single nucleotide variation across the hominid phylogeny. We find that CNV diversity partially correlates with single nucleotide diversity (r2 = 0.5) and recapitulates the phylogeny of apes with few exceptions. Duplications significantly outpace deletions (2.8-fold). The load of segregating duplications remains significantly higher in bonobos, Western chimpanzees, and Sumatran orangutans—populations that have experienced recent genetic bottlenecks (P = 0.0014, 0.02, and 0.0088, respectively). The rate of fixed deletion has been more clocklike with the exception of the chimpanzee lineage, where we observe a twofold increase in the chimpanzee–bonobo ancestor (P = 4.79 × 10−9) and increased deletion load among Western chimpanzees (P = 0.002). The latter includes the first genomic disorder in a chimpanzee with features resembling Smith-Magenis syndrome mediated by a chimpanzee-specific increase in segmental duplication complexity. We hypothesize that demographic effects, such as bottlenecks, have contributed to larger and more gene-rich segments being deleted in the chimpanzee lineage and that this effect, more generally, may account for episodic bursts in CNV during hominid evolution.P.H.S. is supported by a Howard Hughes International Student Fellowship.T.M.B. is supported by an ERC Starting Grant (260372). T.M.B. is an ICREA Research Investigator (Institut Catala d’Estudis i Recerca Avancats de la Generalitat de Catalunya). This work was supported, in part, by U.S. National Institutes of Health (NIH) grant HG002385 to E.E.E., BFU2009-13409-C02-02 to J.P.M., and MICINN (Spain) BFU2011-28549 to T.M.B. E.E.E. is an investigator of the Howard Hughes Medical Institut

Evolution and diversity of copy number variation in the great ape lineage

Copy number variation (CNV) contributes to disease and has restructured the genomes of great apes. The diversity and rate of this process, however, have not been extensively explored among great ape lineages. We analyzed 97 deeply sequenced great ape and human genomes and estimate 16% (469 Mb) of the hominid genome has been affected by recent CNV. We identify a comprehensive set of fixed gene deletions (n = 340) and duplications (n = 405) as well as >13.5 Mb of sequence that has been specifically lost on the human lineage. We compared the diversity and rates of copy number and single nucleotide variation across the hominid phylogeny. We find that CNV diversity partially correlates with single nucleotide diversity (r2 = 0.5) and recapitulates the phylogeny of apes with few exceptions. Duplications significantly outpace deletions (2.8-fold). The load of segregating duplications remains significantly higher in bonobos, Western chimpanzees, and Sumatran orangutans—populations that have experienced recent genetic bottlenecks (P = 0.0014, 0.02, and 0.0088, respectively). The rate of fixed deletion has been more clocklike with the exception of the chimpanzee lineage, where we observe a twofold increase in the chimpanzee–bonobo ancestor (P = 4.79 × 10−9) and increased deletion load among Western chimpanzees (P = 0.002). The latter includes the first genomic disorder in a chimpanzee with features resembling Smith-Magenis syndrome mediated by a chimpanzee-specific increase in segmental duplication complexity. We hypothesize that demographic effects, such as bottlenecks, have contributed to larger and more gene-rich segments being deleted in the chimpanzee lineage and that this effect, more generally, may account for episodic bursts in CNV during hominid evolution.P.H.S. is supported by a Howard Hughes International Student Fellowship.T.M.B. is supported by an ERC Starting Grant (260372). T.M.B. is an ICREA Research Investigator (Institut Catala d’Estudis i Recerca Avancats de la Generalitat de Catalunya). This work was supported, in part, by U.S. National Institutes of Health (NIH) grant HG002385 to E.E.E., BFU2009-13409-C02-02 to J.P.M., and MICINN (Spain) BFU2011-28549 to T.M.B. E.E.E. is an investigator of the Howard Hughes Medical Institut

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization

Abstract Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF’s RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted

A cross-disorder dosage sensitivity map of the human genome

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the prop-erties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 trip-losensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will pro-vide broad utility for human disease research and clinical genetics