Patterns of genic intolerance of rare copy number variation in 59,898 human exomes

Increased burden of rare protein-truncating variants in constrained, brain-specific and synaptic genes in extremely impulsively violent males with antisocial personality disorder

The genetic correlates of extreme impulsive violence are poorly understood, and there have been few studies that have characterized a large group of affected individuals both clinically and genetically. We performed whole exome sequencing (WES) in 290 males with the life-course-persistent, extremely impulsively violent form of antisocial personality disorder (APD) and analyzed the spectrum of rare protein-truncating variants (rPTVs). Comparisons were made with 314 male controls and publicly available genotype data. Functional annotation tools were used for biological interpretation. Participants were significantly more likely to harbor rPTVs in genes that are intolerant to loss-of-function variants (odds ratio [OR] 2.06; p < 0.001), specifically expressed in brain (OR 2.80; p = 0.036) and enriched for those involved in neurotransmitter transport and synaptic processes. In 60 individuals (20%), we identified rPTVs that we classified as clinically relevant based on their clinical associations, biological function and gene expression patterns. Of these, 37 individuals harbored rPTVs in 23 genes that are associated with a monogenic neurological disorder, and 23 individuals harbored rPTVs in 20 genes reportedly intolerant to loss-of-function variants. The analysis presents evidence in support of a model where presence of either one or several private, functionally relevant mutations contribute significantly to individual risk of life-course-persistent APD and reveals multiple individuals who could be affected by clinically unrecognized neuropsychiatric Mendelian disease. Thus, Mendelian diseases and increased rPTV burden may represent important factors for the development of extremely impulsive violent life-course-persistent forms of APD irrespective of their clinical presentation.

Deep multiple-instance learning accurately predicts gene haploinsufficiency and deletion pathogenicity

Copy number losses (deletions) are a major contributor to the etiology of severe genetic disorders. Although haploinsufficient genes play a critical role in deletion pathogenicity, current methods for deletion pathogenicity prediction fail to integrate multiple lines of evidence for haploinsufficiency at the gene level, limiting their power to pinpoint deleterious deletions associated with genetic disorders. Here we introduce DosaCNV, a deep multiple-instance learning framework that, for the first time, models deletion pathogenicity jointly with gene haploinsufficiency. By integrating over 30 gene-level features potentially predictive of haploinsufficiency, DosaCNV shows unmatched performance in prioritizing pathogenic deletions associated with a broad spectrum of genetic disorders. Furthermore, DosaCNV outperforms existing methods in predicting gene haploinsufficiency even though it is not trained on known haploinsufficient genes. Finally, DosaCNV leverages a state-of-the-art technique to quantify the contributions of individual gene-level features to haploinsufficiency, allowing for human-understandable explanations of model predictions. Altogether, DosaCNV is a powerful computational tool for both fundamental and translational research.

Chromosomal Aberrations in Pediatric Patients With Moderate/Severe Developmental Delay/Intellectual Disability With Abundant Phenotypic Heterogeneities: A Single-Center Study

BACKGROUND

This study aimed to examine the clinical usefulness of chromosome microarray (CMA) for selective implementation in patients with unexplained moderate or severe developmental delay/intellectual disability (DD/ID) and/or combined with different dysphonic features in the Han Chinese population.

METHODS

We retrospectively analyzed data on 122 pediatric patients with unexplained isolated moderate/severe DD/ID with or without autism spectrum disorders, epilepsy, dystonia, and congenital abnormalities from a single-center neurorehabilitation clinic in southern China.

RESULTS

A total of 46 probands (37.7%) had abnormal CMA results among the 122 study patients. With the exclusion of aneuploidies, uniparental disomies, and multiple homozygotes, 37 patients harbored 39 pathogenic copy number variations (pCNVs) (median [interquartile range] size: 3.57 [1.6 to 7.1] Mb; 33 deletions and 6 duplications), enriched in chromosomes 5, 7, 15, 17, and 22, with a markedly high prevalence of Angelman/Prader-Willi syndrome (24.3% [nine of 37]). Three rare deletions in the regions 5q33.2q34, 17p13.2, and 13q33.2 were reported, with specific delineation of clinical phenotypes. The frequencies of pCNVs were 18%, 33.3%, 38.89%, 41.67%, and 100% for patients with 1, 2, 3, 4, and 5 study phenotypes, respectively; patients with more concomitant abnormalities in the heart, brain, craniofacial region, and/or other organs had a higher CMA diagnostic yield and pCNV prevalence (P < 0.05).

CONCLUSIONS

Clinical application of CMA as a first-tier test among patients with moderate/severe DD/ID combined with congenital structural anomalies improved diagnostic yields and the quality of clinical management in this series of patients.

GATK-gCNV enables the discovery of rare copy number variants from exome sequencing data

Copy number variants (CNVs) are major contributors to genetic diversity and disease. While standardized methods, such as the genome analysis toolkit (GATK), exist for detecting short variants, technical challenges have confounded uniform large-scale CNV analyses from whole-exome sequencing (WES) data. Given the profound impact of rare and de novo coding CNVs on genome organization and human disease, we developed GATK-gCNV, a flexible algorithm to discover rare CNVs from sequencing read-depth information, complete with open-source distribution via GATK. We benchmarked GATK-gCNV in 7,962 exomes from individuals in quartet families with matched genome sequencing and microarray data, finding up to 95% recall of rare coding CNVs at a resolution of more than two exons. We used GATK-gCNV to generate a reference catalog of rare coding CNVs in WES data from 197,306 individuals in the UK Biobank, and observed strong correlations between per-gene CNV rates and measures of mutational constraint, as well as rare CNV associations with multiple traits. In summary, GATK-gCNV is a tunable approach for sensitive and specific CNV discovery in WES data, with broad applications.

Copy Number Variations and Schizophrenia

Schizophrenia is a neurodevelopmental disorder with genetic and environmental factors involved in its aetiology. Genetic liability contributing to the development of schizophrenia is a subject of extensive research activity, as reliable data regarding its aetiology would enable the improvement of its therapy and the development of new methods of treatment. A multitude of studies in this field focus on genetic variants, such as copy number variations (CNVs) or single-nucleotide variants (SNVs). Certain genetic disorders caused by CNVs including 22q11.2 microdeletion syndrome, Burnside-Butler syndrome (15q11.2 BP1-BP2 microdeletion) or 1q21.1 microduplication/microdeletion syndrome are associated with a higher risk of developing schizophrenia. In this article, we provide a unifying framework linking these CNVs and their associated genetic disorders with schizophrenia and its various neural and behavioural abnormalities.

Diverse monogenic subforms of human spermatogenic failure

Non-obstructive azoospermia (NOA) is the most severe form of male infertility and typically incurable. Defining the genetic basis of NOA has proven challenging, and the most advanced classification of NOA subforms is not based on genetics, but simple description of testis histology. In this study, we exome-sequenced over 1000 clinically diagnosed NOA cases and identified a plausible recessive Mendelian cause in 20%. We find further support for 21 genes in a 2-stage burden test with 2072 cases and 11,587 fertile controls. The disrupted genes are primarily on the autosomes, enriched for undescribed human "knockouts", and, for the most part, have yet to be linked to a Mendelian trait. Integration with single-cell RNA sequencing data shows that azoospermia genes can be grouped into molecular subforms with synchronized expression patterns, and analogs of these subforms exist in mice. This analysis framework identifies groups of genes with known roles in spermatogenesis but also reveals unrecognized subforms, such as a set of genes expressed across mitotic divisions of differentiating spermatogonia. Our findings highlight NOA as an understudied Mendelian disorder and provide a conceptual structure for organizing the complex genetics of male infertility, which may provide a rational basis for disease classification.

Complex genomic rearrangements: an underestimated cause of rare diseases

Complex genomic rearrangements (CGRs) are known contributors to disease but are often missed during routine genetic screening. Identifying CGRs requires (i) identifying copy number variants (CNVs) concurrently with inversions, (ii) phasing multiple breakpoint junctions incis, as well as (iii) detecting and resolving structural variants (SVs) within repeats. We demonstrate how combining cytogenetics and new sequencing methodologies is being successfully applied to gain insights into the genomic architecture of CGRs. In addition, we review CGR patterns and molecular features revealed by studying constitutional genomic disorders. These data offer invaluable lessons to individuals interested in investigating CGRs, evaluating their clinical relevance and frequency, as well as assessing their impact(s) on rare genetic diseases.

Neurodevelopmental copy-number variants: A roadmap to improving outcomes by uniting patient advocates, researchers, and clinicians for collective impact

Copy-number variants and structural variants (CNVs/SVs) drive many neurodevelopmental-related disorders. While many neurodevelopmental-related CNVs/SVs give rise to complex phenotypes, the overlap in phenotypic presentation between independent CNVs can be extensive and provides a motivation for shared approaches. This confluence at the level of clinical phenotype implies convergence in at least some aspects of the underlying genomic mechanisms. With this perspective, our Commission on Novel Technologies for Neurodevelopmental CNVs asserts that the time has arrived to approach neurodevelopmental-related CNVs/SVs as a class of disorders that can be identified, investigated, and treated on the basis of shared mechanisms and/or pathways (e.g., molecular, neurological, or developmental). To identify common etiologic mechanisms among uncommon neurodevelopmental-related disorders and to potentially identify common therapies, it is paramount for teams of scientists, clinicians, and patients to unite their efforts. We bring forward novel, collaborative, and integrative strategies to translational CNV/SV research that engages diverse stakeholders to help expedite therapeutic outcomes. We articulate a clear vision for piloted roadmap strategies to reduce patient/caregiver burden and redundancies, increase efficiency, avoid siloed data, and accelerate translational discovery across CNV/SV-based syndromes.

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 properties 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 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Large-scale discovery of novel neurodevelopmental disorder-related genes through a unified analysis of single-nucleotide and copy number variants

Background

Previous large-scale studies of de novo variants identified a number of genes associated with neurodevelopmental disorders (NDDs); however, it was also predicted that many NDD-associated genes await discovery. Such genes can be discovered by integrating copy number variants (CNVs), which have not been fully considered in previous studies, and increasing the sample size.

Methods

We first constructed a model estimating the rates of de novo CNVs per gene from several factors such as gene length and number of exons. Second, we compiled a comprehensive list of de novo single-nucleotide variants (SNVs) in 41,165 individuals and de novo CNVs in 3675 individuals with NDDs by aggregating our own and publicly available datasets, including denovo-db and the Deciphering Developmental Disorders study data. Third, summing up the de novo CNV rates that we estimated and SNV rates previously established, gene-based enrichment of de novo deleterious SNVs and CNVs were assessed in the 41,165 cases. Significantly enriched genes were further prioritized according to their similarity to known NDD genes using a deep learning model that considers functional characteristics (e.g., gene ontology and expression patterns).

Results

We identified a total of 380 genes achieving statistical significance (5% false discovery rate), including 31 genes affected by de novo CNVs. Of the 380 genes, 52 have not previously been reported as NDD genes, and the data of de novo CNVs contributed to the significance of three genes (GLTSCR1, MARK2, and UBR3). Among the 52 genes, we reasonably excluded 18 genes [a number almost identical to the theoretically expected false positives (i.e., 380 × 0.05 = 19)] given their constraints against deleterious variants and extracted 34 “plausible” candidate genes. Their validity as NDD genes was consistently supported by their similarity in function and gene expression patterns to known NDD genes. Quantifying the overall similarity using deep learning, we identified 11 high-confidence (> 90% true-positive probabilities) candidate genes: HDAC2, SUPT16H, HECTD4, CHD5, XPO1, GSK3B, NLGN2, ADGRB1, CTR9, BRD3, and MARK2.

Conclusions

We identified dozens of new candidates for NDD genes. Both the methods and the resources developed here will contribute to the further identification of novel NDD-associated genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01042-w.

Copy number variants calling from WES data through eXome hidden Markov model (XHMM) identifies additional 2.5% pathogenic genomic imbalances smaller than 30 kb undetected by array-CGH

It has been estimated that Copy Number Variants (CNVs) account for 10%-20% of patients affected by Developmental Disorder (DD)/Intellectual Disability (ID). Although array comparative genomic hybridization (array-CGH) represents the gold-standard for the detection of genomic imbalances, common Agilent array-CGH 4 × 180 kb arrays fail to detect CNVs smaller than 30 kb. Whole Exome sequencing (WES) is becoming the reference application for the detection of gene variants and makes it possible also to infer genomic imbalances at single exon resolution. However, the contribution of small CNVs in DD/ID is still underinvestigated. We made use of the eXome Hidden Markov Model (XHMM) software, a tool utilized by the ExAC consortium, to detect CNVs from whole exome sequencing data, in a cohort of 200 unsolved DD/DI patients after array-CGH and WES-based single nucleotide/indel variant analyses. In five out of 200 patients (2.5%), we identified pathogenic CNV(s) smaller than 30 kb, ranging from one to six exons. They included two heterozygous deletions in TCF4 and STXBP1 and three homozygous deletions in PPT1, CLCN2, and PIGN. After reverse phenotyping, all variants were reported as causative. This study shows the interest in applying sequencing-based CNV detection, from available WES data, to reduce the diagnostic odyssey of additional patients unsolved DD/DI patients and compare the CNV-detection yield of Agilent array-CGH 4 × 180kb versus whole exome sequencing.

Chromosomal microarray analysis of 410 Han Chinese patients with autism spectrum disorder or unexplained intellectual disability and developmental delay

Copy number variants (CNVs) are recognized as a crucial genetic cause of neurodevelopmental disorders (NDDs). Chromosomal microarray analysis (CMA), the first-tier diagnostic test for individuals with NDDs, has been utilized to detect CNVs in clinical practice, but most reports are still from populations of European ancestry. To contribute more worldwide clinical genomics data, we investigated the genetic etiology of 410 Han Chinese patients with NDDs (151 with autism and 259 with unexplained intellectual disability (ID) and developmental delay (DD)) using CMA (Affymetrix) after G-banding karyotyping. Among all the NDD patients, 109 (26.6%) carried clinically relevant CNVs or uniparental disomies (UPDs), and 8 (2.0%) had aneuploidies (6 with trisomy 21 syndrome, 1 with 47,XXY, 1 with 47,XYY). In total, we found 129 clinically relevant CNVs and UPDs, including 32 CNVs in 30 ASD patients, and 92 CNVs and 5 UPDs in 79 ID/DD cases. When excluding the eight patients with aneuploidies, the diagnostic yield of pathogenic and likely pathogenic CNVs and UPDs was 20.9% for all NDDs (84/402), 3.3% in ASD (5/151), and 31.5% in ID/DD (79/251). When aneuploidies were included, the diagnostic yield increased to 22.4% for all NDDs (92/410), and 33.6% for ID/DD (87/259). We identified a de novo CNV in 14.9% (60/402) of subjects with NDDs. Interestingly, a higher diagnostic yield was observed in females (31.3%, 40/128) compared to males (16.1%, 44/274) for all NDDs (P = 4.8 × 10^-4), suggesting that a female protective mechanism exists for deleterious CNVs and UPDs.

Overcoming constraints on the detection of recessive selection in human genes from population frequency data

The identification of genes that evolve under recessive natural selection is a long-standing goal of population genetics research that has important applications to the discovery of genes associated with disease. We found that commonly used methods to evaluate selective constraint at the gene level are highly sensitive to genes under heterozygous selection but ubiquitously fail to detect recessively evolving genes. Additionally, more sophisticated likelihood-based methods designed to detect recessivity similarly lack power for a human gene of realistic length from current population sample sizes. However, extensive simulations suggested that recessive genes may be detectable in aggregate. Here, we offer a method informed by population genetics simulations designed to detect recessive purifying selection in gene sets. Applying this to empirical gene sets produced significant enrichments for strong recessive selection in genes previously inferred to be under recessive selection in a consanguineous cohort and in genes involved in autosomal recessive monogenic disorders.

Using induced pluripotent stem cells to investigate human neuronal phenotypes in 1q21.1 deletion and duplication syndrome

Copy Number Variation (CNV) at the 1q21.1 locus is associated with a range of neurodevelopmental and psychiatric disorders in humans, including abnormalities in head size and motor deficits. Yet, the functional consequences of these CNVs (both deletion and duplication) on neuronal development remain unknown. To determine the impact of CNV at the 1q21.1 locus on neuronal development, we generated induced pluripotent stem cells from individuals harbouring 1q21.1 deletion or duplication and differentiated them into functional cortical neurons. We show that neurons with 1q21.1 deletion or duplication display reciprocal phenotype with respect to proliferation, differentiation potential, neuronal maturation, synaptic density and functional activity. Deletion of the 1q21.1 locus was also associated with an increased expression of lower cortical layer markers. This difference was conserved in the mouse model of 1q21.1 deletion, which displayed altered corticogenesis. Importantly, we show that neurons with 1q21.1 deletion and duplication are associated with differential expression of calcium channels and demonstrate that physiological deficits in neurons with 1q21.1 deletion or duplication can be pharmacologically modulated by targeting Ca²⁺ channel activity. These findings provide biological insight into the neuropathological mechanism underlying 1q21.1 associated brain disorder and indicate a potential target for therapeutic interventions.

Trio-Based Low-Pass Genome Sequencing Reveals Characteristics and Significance of Rare Copy Number Variants in Prenatal Diagnosis

Background: Low-pass genome sequencing (GS) detects clinically significant copy number variants (CNVs) in prenatal diagnosis. However, detection at improved resolutions leads to an increase in the number of CNVs identified, increasing the difficulty of clinical interpretation and management.

Methods: Trio-based low-pass GS was performed in 315 pregnancies undergoing invasive testing. Rare CNVs detected in the fetuses were investigated. The characteristics of rare CNVs were described and compared to curated CNVs in other studies.

Results: A total of 603 rare CNVs, namely, 597 constitutional and 6 mosaic CNVs, were detected in 272 fetuses (272/315, 86.3%), providing 1.9 rare CNVs per fetus (603/315). Most CNVs were smaller than 1 Mb (562/603, 93.2%), while 1% (6/603) were mosaic. Forty-six de novo (7.6%, 46/603) CNVs were detected in 11.4% (36/315) of the cases. Eighty-four CNVs (74 fetuses, 23.5%) involved disease-causing genes of which the mode of inheritance was crucial for interpretation and assessment of recurrence risk. Overall, 31 pathogenic/likely pathogenic CNVs were detected, among which 25.8% (8/31) were small (<100 kb; n = 3) or mosaic CNVs (n = 5).

Conclusion: We examined the landscape of rare CNVs with parental inheritance assignment and demonstrated that they occur frequently in prenatal diagnosis. This information has clinical implications regarding genetic counseling and consideration for trio-based CNV analysis.

Prevalence and phenotypic impact of rare potentially damaging variants in autism spectrum disorder

BACKGROUND

The Autism Sequencing Consortium identified 102 high-confidence autism spectrum disorder (ASD) genes, showing that individuals with ASD and with potentially damaging single nucleotide variation (pdSNV) in these genes had lower cognitive levels and delayed age at walking, when compared to ASD participants without pdSNV. Here, we made use of a Swedish sample of individuals with ASD (called PAGES, for Population-Based Autism Genetics & Environment Study) to evaluate the frequency of pdSNV and their impact on medical and psychiatric phenotypes, using an epidemiological frame and universal health reporting. We then combine findings with those for potentially damaging copy number variation (pdCNV).

METHODS

SNV and CNV calls were generated from whole-exome sequencing and chromosome microarray data, respectively. Birth and medical register data were used to collect phenotypes.

RESULTS

Of 808 individuals assessed by sequencing, 69 (9%) had pdSNV in the 102 ASC genes, and 144 (18%) had pdSNV in the 102 ASC genes or in a larger set of curated neurodevelopmental genes (from the Deciphering Developmental Disorders study, the gene2phenotype database, and the Radboud University gene lists). Three or more individuals had pdSNV in GRIN2B, POGZ, SATB1, DYNC1H1, SCN8A, or CREBBP. In comparison, out of the 996 individuals from whom CNV were called, 105 (11%) carried one or more pdCNV, including four or more individuals with CNV in the recurrent 15q11q13, 22q11.2, and 16p11.2 loci. Carriers of pdSNV were more likely to have intellectual disability (ID) and epilepsy, while carriers of pdCNV showed increased rates of congenital anomalies and scholastic skill disorders. Carriers of either pdSNV or pdCNV were more likely to have ID, scholastic skill disorders, and epilepsy.

LIMITATIONS

The cohort only included individuals with autistic disorder, the more severe form of ASD, and phenotypes are defined from medical registers. Not all genes studied are definitively ASD genes, and we did not have de novo information to aid in classification.

CONCLUSIONS

In this epidemiological sample, rare pdSNV were more common than pdCNV and the combined yield of potentially damaging variation was substantial at 27%. The results provide compelling rationale for the use of high-throughout sequencing as part of routine clinical workup for ASD and support the development of precision medicine in ASD.

Size matters: Large copy number losses in Hirschsprung disease patients reveal genes involved in enteric nervous system development

Hirschsprung disease (HSCR) is a complex genetic disease characterized by absence of ganglia in the intestine. HSCR etiology can be explained by a unique combination of genetic alterations: rare coding variants, predisposing haplotypes and Copy Number Variation (CNV). Approximately 18% of patients have additional anatomical malformations or neurological symptoms (HSCR-AAM). Pinpointing the responsible culprits within a CNV is challenging as often many genes are affected. Therefore, we selected candidate genes based on gene enrichment strategies using mouse enteric nervous system transcriptomes and constraint metrics. Next, we used a zebrafish model to investigate whether loss of these genes affects enteric neuron development in vivo.

This study included three groups of patients, two groups without coding variants in disease associated genes: HSCR-AAM and HSCR patients without associated anomalies (HSCR-isolated). The third group consisted of all HSCR patients in which a confirmed pathogenic rare coding variant was identified. We compared these patient groups to unaffected controls. Predisposing haplotypes were determined, confirming that every HSCR subgroup had increased contributions of predisposing haplotypes, but their contribution was highest in isolated HSCR patients without RET coding variants. CNV profiling proved that specifically HSCR-AAM patients had larger Copy Number (CN) losses. Gene enrichment strategies using mouse enteric nervous system transcriptomes and constraint metrics were used to determine plausible candidate genes located within CN losses. Validation in zebrafish using CRISPR/Cas9 targeting confirmed the contribution of UFD1L, TBX2, SLC8A1, and MAPK8 to ENS development. In addition, we revealed epistasis between reduced Ret and Gnl1 expression and between reduced Ret and Tubb5 expression in vivo. Rare large CN losses—often de novo—contribute to HSCR in HSCR-AAM patients. We proved the involvement of six genes in enteric nervous system development and Hirschsprung disease.

Hirschsprung disease is a congenital disorder characterized by the absence of intestinal neurons in the distal part of the intestine. It is a complex genetic disorder in which multiple variations in our genome combined, result in disease. One of these variations are Copy Number Variations (CNVs): large segments of our genome that are duplicated or deleted. Patients often have Hirschsprung disease without other symptoms. However, a proportion of patients has additional associated anatomical malformations and neurological symptoms. We found that CNVs, present in patients with associated anomalies, are more often larger compared to unaffected controls or Hirschsprung patients without other symptoms. Furthermore, Copy Number (CN) losses are enriched for constrained coding regions (CCR; genes usually not impacted by genomic alterations in unaffected controls) of which the expression is higher in the developing intestinal neurons compared to the intestine. We modelled loss of these candidate genes in zebrafish by disrupting the zebrafish orthologues by genome editing. For several genes this resulted in changes in intestinal neuron development, reminiscent of HSCR observed in patients. The results presented here highlight the importance of Copy Number profiling, zebrafish validation and evaluating all CCR expressed in developing intestinal neurons during diagnostic evaluation.

CNVIntegrate: the first multi-ethnic database for identifying copy number variations associated with cancer

Human copy number variations (CNVs) and copy number alterations (CNAs) are DNA segments (>1000 base pairs) of duplications or deletions with respect to the reference genome, potentially causing genomic imbalance leading to diseases such as cancer. CNVs further cause genetic diversity in healthy populations and are predominant drivers of gene/genome evolution. Initiatives have been taken by the research community to establish large-scale databases to comprehensively characterize CNVs in humans. Exome Aggregation Consortium (ExAC) is one such endeavor that catalogs CNVs, of nearly 60 000 healthy individuals across five demographic clusters. Furthermore, large projects such as the Catalogue of Somatic Mutations in Cancer (COSMIC) and the Cancer Cell Line Encyclopedia (CCLE) combine CNA data from cancer-affected individuals and large panels of human cancer cell lines, respectively. However, we lack a structured and comprehensive CNV/CNA resource including both healthy individuals and cancer patients across large populations. CNVIntegrate is the first web-based system that hosts CNV and CNA data from both healthy populations and cancer patients, respectively, and concomitantly provides statistical comparisons between copy number frequencies of multiple ethnic populations. It further includes, for the first time, well-cataloged CNV and CNA data from Taiwanese healthy individuals and Taiwan Breast Cancer data, respectively, along with imported resources from ExAC, COSMIC and CCLE. CNVIntegrate offers a CNV/CNA-data hub for structured information retrieval for clinicians and scientists towards important drug discoveries and precision treatments. Database URL: http://cnvintegrate.cgm.ntu.edu.tw/.

The physiological and pathological functions of VEGFR3 in cardiac and lymphatic development and related diseases

Vascular endothelial growth factor receptors (VEGFRs) are part of the evolutionarily conserved VEGF signalling pathways that regulate the development and maintenance of the body's cardiovascular and lymphovascular systems. VEGFR3, encoded by the FLT4 gene, has an indispensable and well-characterized function in development and establishment of the lymphatic system. Autosomal dominant VEGFR3 mutations, that prevent the receptor functioning as a homodimer, cause one of the major forms of hereditary primary lymphoedema; Milroy disease. Recently, we and others have shown that FLT4 variants, distinct to those observed in Milroy disease cases, predispose individuals to Tetralogy of Fallot, the most common cyanotic congenital heart disease, demonstrating a novel function for VEGFR3 in early cardiac development. Here, we examine the familiar and emerging roles of VEGFR3 in the development of both lymphovascular and cardiovascular systems, respectively, compare how distinct genetic variants in FLT4 lead to two disparate human conditions, and highlight the research still required to fully understand this multifaceted receptor.

Origins of human genetics. A personal perspective

Genetics evolved as a field of science after 1900 with new theories being derived from experiments obtained in fruit flies, bacteria, and viruses. This personal account suggests that the origins of human genetics can best be traced to the years 1949 to 1959. Several genetic scientific advances in genetics in 1949 yielded results directly relating to humans for the first time, except for a few earlier observations. In 1949 the first textbook of human genetics was published, the American Journal of Human Genetics was founded, and in the previous year the American Society of Human Genetics. In 1940 in Britain a textbook entitled Introduction to Medical Genetics served as a foundation for introducing genetic aspects into medicine. The introduction of new methods for analyzing chromosomes and new biochemical assays using cultured cells in 1959 and subsequent years revealed that many human diseases, including cancer, have genetic causes. It became possible to arrive at a precise cause-related genetic diagnosis. As a result the risk of occurrence or re-occurrence of a disease within a family could be assessed correctly. Genetic counseling as a new concept became a basis for improved patient care. Taken together the advances in medically orientated genetic research and patient care since 1949 have resulted in human genetics being both, a basic medical and a basic biological science. Prior to 1949 genetics was not generally viewed in a medical context. Although monogenic human diseases were recognized in 1902, their occurrence and distribution were considered mainly at the population level.

Voltage-Gated Sodium Channel β1 Gene: An Overview

BACKGROUND

Voltage-gated sodium channels are protein complexes composed of 2 subunits, namely, pore-forming α- and regulatory β-subunits. A β-subunit consists of 5 proteins encoded by 4 genes (i.e., SCN1B-SCN4B).

SUMMARY

β1-Subunits regulate sodium ion channel functions, including gating properties, subcellular localization, and kinetics. Key Message: Sodium channel β1- and its variant β1B-subunits are encoded by SCN1B. These variants are associated with many human diseases, such as epilepsy, Brugada syndrome, Dravet syndrome, and cancers. On the basis of previous research, we aimed to provide an overview of the structure, expression, and involvement of SCN1B in physiological processes and focused on its role in diseases.

Novel MYT1 variants expose the complexity of oculo-auriculo-vertebral spectrum genetic mechanisms

Oculo-auriculo-vertebral spectrum (OAVS) is a developmental disorder characterized by anomalies mainly involving the structures derived from the first and second pharyngeal arches. The spectrum presents with heterogeneous clinical features and complex etiology with genetic factors not yet completely understood. To date, MYT1 is the most important gene unambiguously associated with the spectrum and with functional data confirmation. In this work, we aimed to identify new single nucleotide variants (SNVs) affecting MYT1 in a cohort of 73 Brazilian patients diagnosed with OAVS. In addition, we investigated copy number variations (CNVs) encompassing this gene or its cis-regulatory elements and compared the frequency of these events in patients versus a cohort of 455 Brazilian control individuals. A new SNV, predicted as likely deleterious, was identified in five unrelated patients with OAVS. All five patients presented hearing impairment and orbital asymmetry suggesting an association with the variant. CNVs near MYT1, located in its neighboring topologically associating domain (TAD), were found to be enriched in patients when compared to controls, indicating a possible involvement of this region with OAVS pathogenicity. Our findings highlight the genetic complexity of the spectrum that seems to involve more than one variant type and inheritance patterns.

A pediatric case of primary ciliary dyskinesia caused by novel copy number variation in PIH1D3

An 11-month-old boy with productive cough was referred to our hospital. He had nasal obstruction immediately after birth, and wheezing, wet cough, and rhinorrhea were observed daily after the neonatal period. Clinical and imaging findings revealed secretory otitis media, chronic sinusitis, and bronchiectasis. Primary ciliary dyskinesia was suspected. Transmission electron microscopy of nasal cilia showed defects of the outer and inner dynein arms. Genetic examinations of the family revealed copy number variation in PIH1 domain-containing 3 (PIH1D3) in the proband and mother. This is the first report of a Japanese patient with primary ciliary dyskinesia caused by copy number variation in PIH1D3.

Dosage-sensitive functions in embryonic development drove the survival of genes on sex-specific chromosomes in snakes, birds, and mammals

Different ancestral autosomes independently evolved into sex chromosomes in snakes, birds, and mammals. In snakes and birds, females are ZW and males are ZZ; in mammals, females are XX and males are XY. Although X and Z Chromosomes retain nearly all ancestral genes, sex-specific W and Y Chromosomes suffered extensive genetic decay. In both birds and mammals, the genes that survived on sex-specific chromosomes are enriched for broadly expressed, dosage-sensitive regulators of gene expression, subject to strong purifying selection. To gain deeper insight into the processes that govern survival on sex-specific chromosomes, we carried out a meta-analysis of survival across 41 species-three snakes, 24 birds, and 14 mammals-doubling the number of ancestral genes under investigation and increasing our power to detect enrichments among survivors relative to nonsurvivors. Of 2564 ancestral genes, representing an eighth of the ancestral amniote genome, only 324 survive on present-day sex-specific chromosomes. Survivors are enriched for dosage-sensitive developmental processes, particularly development of neural crest-derived structures, such as the face. However, there was no enrichment for expression in sex-specific tissues, involvement in sex determination or gonadogenesis pathways, or conserved sex-biased expression. Broad expression and dosage sensitivity contributed independently to gene survival, suggesting that pleiotropy imposes additional constraints on the evolution of dosage compensation. We propose that maintaining the viability of the heterogametic sex drove gene survival on amniote sex-specific chromosomes, and that subtle modulation of the expression of survivor genes and their autosomal orthologs has disproportionately large effects on development and disease.

Effect Sizes of Deletions and Duplications on Autism Risk Across the Genome

OBJECTIVE

Deleterious copy number variants (CNVs) are identified in up to 20% of individuals with autism. However, levels of autism risk conferred by most rare CNVs remain unknown. The authors recently developed statistical models to estimate the effect size on IQ of all CNVs, including undocumented ones. In this study, the authors extended this model to autism susceptibility.

METHODS

The authors identified CNVs in two autism populations (Simons Simplex Collection and MSSNG) and two unselected populations (IMAGEN and Saguenay Youth Study). Statistical models were used to test nine quantitative variables associated with genes encompassed in CNVs to explain their effects on IQ, autism susceptibility, and behavioral domains.

RESULTS

The "probability of being loss-of-function intolerant" (pLI) best explains the effect of CNVs on IQ and autism risk. Deleting 1 point of pLI decreases IQ by 2.6 points in autism and unselected populations. The effect of duplications on IQ is threefold smaller. Autism susceptibility increases when deleting or duplicating any point of pLI. This is true for individuals with high or low IQ and after removing de novo and known recurrent neuropsychiatric CNVs. When CNV effects on IQ are accounted for, autism susceptibility remains mostly unchanged for duplications but decreases for deletions. Model estimates for autism risk overlap with previously published observations. Deletions and duplications differentially affect social communication, behavior, and phonological memory, whereas both equally affect motor skills.

CONCLUSIONS

Autism risk conferred by duplications is less influenced by IQ compared with deletions. The model applied in this study, trained on CNVs encompassing >4,500 genes, suggests highly polygenic properties of gene dosage with respect to autism risk and IQ loss. These models will help to interpret CNVs identified in the clinic.

Performance of In Silico Prediction Tools for the Detection of Germline Copy Number Variations in Cancer Predisposition Genes in 4208 Female Index Patients with Familial Breast and Ovarian Cancer

Simple Summary

The identification of germline copy number variants (CNVs) by targeted nextgeneration sequencing frequently relies on in silico prediction tools with unknown sensitivities. We investigated the performances of four in silico CNV prediction tools in 17 cancer predisposition genes in a large series of 4208 female index patients with familial breast and/or ovarian cancer. We identified 77 CNVs in 76 out of 4208 patients; six CNVs were missed by at least one of the prediction tools. Experimental verification of in silico predicted CNVs is required due to high frequencies of false positive predictions. For female index patients with familial breast and/or ovarian cancer, CNV detection should not be restricted to BRCA1/2 due to the relevant proportion of CNVs in further cancer predisposition genes.

Abstract

The identification of germline copy number variants (CNVs) by targeted next-generation sequencing (NGS) frequently relies on in silico CNV prediction tools with unknown sensitivities. We investigated the performances of four in silico CNV prediction tools, including one commercial (Sophia Genetics DDM) and three non-commercial tools (ExomeDepth, GATK gCNV, panelcn.MOPS) in 17 cancer predisposition genes in 4208 female index patients with familial breast and/or ovarian cancer (BC/OC). CNV predictions were verified via multiplex ligation-dependent probe amplification. We identified 77 CNVs in 76 out of 4208 patients (1.81%); 33 CNVs were identified in genes other than BRCA1/2, mostly in ATM, CHEK2, and RAD51C and less frequently in BARD1, MLH1, MSH2, PALB2, PMS2, RAD51D, and TP53. The Sophia Genetics DDM software showed the highest sensitivity; six CNVs were missed by at least one of the non-commercial tools. The positive predictive values ranged from 5.9% (74/1249) for panelcn.MOPS to 79.1% (72/91) for ExomeDepth. Verification of in silico predicted CNVs is required due to high frequencies of false positive predictions, particularly affecting target regions at the extremes of the GC content or target length distributions. CNV detection should not be restricted to BRCA1/2 due to the relevant proportion of CNVs in further BC/OC predisposition genes.

NCKAP1 Disruptive Variants Lead to a Neurodevelopmental Disorder with Core Features of Autism

NCKAP1/NAP1 regulates neuronal cytoskeletal dynamics and is essential for neuronal differentiation in the developing brain. Deleterious variants in NCKAP1 have been identified in individuals with autism spectrum disorder (ASD) and intellectual disability; however, its clinical significance remains unclear. To determine its significance, we assemble genotype and phenotype data for 21 affected individuals from 20 unrelated families with predicted deleterious variants in NCKAP1. This includes 16 individuals with de novo (n = 8), transmitted (n = 6), or inheritance unknown (n = 2) truncating variants, two individuals with structural variants, and three with potentially disruptive de novo missense variants. We report a de novo and ultra-rare deleterious variant burden of NCKAP1 in individuals with neurodevelopmental disorders which needs further replication. ASD or autistic features, language and motor delay, and variable expression of intellectual or learning disability are common clinical features. Among inherited cases, there is evidence of deleterious variants segregating with neuropsychiatric disorders. Based on available human brain transcriptomic data, we show that NCKAP1 is broadly and highly expressed in both prenatal and postnatal periods and demostrate enriched expression in excitatory neurons and radial glias but depleted expression in inhibitory neurons. Mouse in utero electroporation experiments reveal that Nckap1 loss of function promotes neuronal migration during early cortical development. Combined, these data support a role for disruptive NCKAP1 variants in neurodevelopmental delay/autism, possibly by interfering with neuronal migration early in cortical development.

Developmental Gene Expression Differences between Humans and Mammalian Models

Identifying the molecular programs underlying human organ development and how they differ from model species is key for understanding human health and disease. Developmental gene expression profiles provide a window into the genes underlying organ development and a direct means to compare them across species. We use a transcriptomic resource covering the development of seven organs to characterize the temporal profiles of human genes associated with distinct disease classes and to determine, for each human gene, the similarity of its spatiotemporal expression with its orthologs in rhesus macaque, mouse, rat, and rabbit. We find clear associations between spatiotemporal profiles and the phenotypic manifestations of diseases. We also find that half of human genes differ from their mouse orthologs in their temporal trajectories in at least one of the organs. These include more than 200 genes associated with brain, heart, and liver disease for which mouse models should undergo extra scrutiny.

A Practical Guide for Structural Variation Detection in the Human Genome

Profiling genetic variants-including single nucleotide variants, small insertions and deletions, copy number variations, and structural variations (SVs)-from both healthy individuals and individuals with disease is a key component of genetic and biomedical research. SVs are large-scale changes in the genome and involve breakage and rejoining of DNA fragments. They may affect thousands to millions of nucleotides and can lead to loss, gain, and reshuffling of genes and regulatory elements. SVs are known to impact gene expression and potentially result in altered phenotypes and diseases. Therefore, identifying SVs from the human genomes is particularly important. In this review, I describe advantages and disadvantages of the available high-throughput assays for the discovery of SVs, which are the most challenging genetic alterations to detect. A practical guide is offered to suggest the most suitable strategies for discovering different types of SVs including common germline, rare, somatic, and complex variants. I also discuss factors to be considered, such as cost and performance, for different strategies when designing experiments. Last, I present several approaches to identify potential SV artifacts caused by samples, experimental procedures, and computational analysis. © 2020 Wiley Periodicals LLC.

Clinician-centric diagnosis of rare genetic diseases: performance of a gene pertinence metric in decision support for clinicians

Background

In diagnosis of rare genetic diseases we face a decision as to the degree to which the sequencing lab offers one or more diagnoses based on clinical input provided by the clinician, or the clinician reaches a diagnosis based on the complete set of variants provided by the lab. We tested a software approach to assist the clinician in making the diagnosis based on clinical findings and an annotated genomic variant table, using cases already solved using less automated processes.

Results

For the 81 cases studied (involving 216 individuals), 70 had genetic abnormalities with phenotypes previously described in the literature, and 11 were not described in the literature at the time of analysis (“discovery genes”). These included cases beyond a trio, including ones with different variants in the same gene. In 100% of cases the abnormality was recognized. Of the 70, the abnormality was ranked #1 in 94% of cases, with an average rank 1.1 for all cases. Large CNVs could be analyzed in an integrated analysis, performed in 24 of the cases. The process is rapid enough to allow for periodic reanalysis of unsolved cases.

Conclusions

A clinician-friendly environment for clinical correlation can be provided to clinicians who are best positioned to have the clinical information needed for this interpretation.

The Consequences of Abnormal Gene Dosage: Lessons from Chromosome 18

Accurate interpretation of genomic copy number variation (CNV) remains a challenge and has important consequences for both congenital and late-onset disease. Hemizygosity dosage characterization of the genes on chromosome 18 reveals a spectrum of outcomes ranging from no clinical effect, to risk factors for disease, to both low- and high-penetrance disease. These data are important for accurate and predictive clinical management. Additionally, the potential mechanisms of reduced penetrance due to dosage compensation are discussed as a key to understanding avenues for potential treatment. This review describes the chromosome 18 findings, and discusses the molecular mechanisms that allow haploinsufficiency, reduced penetrance, and dosage compensation.

Mapping and characterization of structural variation in 17,795 human genomes

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline1 to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

Detection of copy-number variations from NGS data using read depth information: a diagnostic performance evaluation

The detection of copy-number variations (CNVs) from NGS data is underexploited as chip-based or targeted techniques are still commonly used. We assessed the performances of a workflow centered on CANOES, a bioinformatics tool based on read depth information. We applied our workflow to gene panel (GP) and whole-exome sequencing (WES) data, and compared CNV calls to quantitative multiplex PCR of short fluorescent fragments (QMSPF) or array comparative genomic hybridization (aCGH) results. From GP data of 3776 samples, we reached an overall positive predictive value (PPV) of 87.8%. This dataset included a complete comprehensive QMPSF comparison of four genes (60 exons) on which we obtained 100% sensitivity and specificity. From WES data, we first compared 137 samples with aCGH and filtered comparable events (exonic CNVs encompassing enough aCGH probes) and obtained an 87.25% sensitivity. The overall PPV was 86.4% following the targeted confirmation of candidate CNVs from 1056 additional WES. In addition, our CANOES-centered workflow on WES data allowed the detection of CNVs with a resolution of single exons, allowing the detection of CNVs that were missed by aCGH. Overall, switching to an NGS-only approach should be cost-effective as it allows a reduction in overall costs together with likely stable diagnostic yields. Our bioinformatics pipeline is available at: https://gitlab.bioinfo-diag.fr/nc4gpm/canoes-centered-workflow .

Functional annotation of rare structural variation in the human brain

Structural variants (SVs) contribute to many disorders, yet, functionally annotating them remains a major challenge. Here, we integrate SVs with RNA-sequencing from human post-mortem brains to quantify their dosage and regulatory effects. We show that genic and regulatory SVs exist at significantly lower frequencies than intergenic SVs. Functional impact of copy number variants (CNVs) stems from both the proportion of genic and regulatory content altered and loss-of-function intolerance of the gene. We train a linear model to predict expression effects of rare CNVs and use it to annotate regulatory disruption of CNVs from 14,891 independent genome-sequenced individuals. Pathogenic deletions implicated in neurodevelopmental disorders show significantly more extreme regulatory disruption scores and if rank ordered would be prioritized higher than using frequency or length alone. This work shows the deleteriousness of regulatory SVs, particularly those altering CTCF sites and provides a simple approach for functionally annotating the regulatory consequences of CNVs.

Properties of structural variants and short tandem repeats associated with gene expression and complex traits

Structural variants (SVs) and short tandem repeats (STRs) comprise a broad group of diverse DNA variants which vastly differ in their sizes and distributions across the genome. Here, we identify genomic features of SV classes and STRs that are associated with gene expression and complex traits, including their locations relative to eGenes, likelihood of being associated with multiple eGenes, associated eGene types (e.g., coding, noncoding, level of evolutionary constraint), effect sizes, linkage disequilibrium with tagging single nucleotide variants used in GWAS, and likelihood of being associated with GWAS traits. We identify a set of high-impact SVs/STRs associated with the expression of three or more eGenes via chromatin loops and show that they are highly enriched for being associated with GWAS traits. Our study provides insights into the genomic properties of structural variant classes and short tandem repeats that are associated with gene expression and human traits.

Mapping and characterization of structural variation in 17,795 human genomes

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline¹ to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

A structural variation reference for medical and population genetics

Structural variants (SVs) rearrange large segments of DNA¹ and can have profound consequences in evolution and human disease^2,3. As national biobanks, disease-association studies, and clinical genetic testing have grown increasingly reliant on genome sequencing, population references such as the Genome Aggregation Database (gnomAD)⁴ have become integral in the interpretation of single-nucleotide variants (SNVs)⁵. However, there are no reference maps of SVs from high-coverage genome sequencing comparable to those for SNVs. Here we present a reference of sequence-resolved SVs constructed from 14,891 genomes across diverse global populations (54% non-European) in gnomAD. We discovered a rich and complex landscape of 433,371 SVs, from which we estimate that SVs are responsible for 25–29% of all rare protein-truncating events per genome. We found strong correlations between natural selection against damaging SNVs and rare SVs that disrupt or duplicate protein-coding sequence, which suggests that genes that are highly intolerant to loss-of-function are also sensitive to increased dosage⁶. We also uncovered modest selection against noncoding SVs in cis-regulatory elements, although selection against protein-truncating SVs was stronger than all noncoding effects. Finally, we identified very large (over one megabase), rare SVs in 3.9% of samples, and estimate that 0.13% of individuals may carry an SV that meets the existing criteria for clinically important incidental findings⁷. This SV resource is freely distributed via the gnomAD browser⁸ and will have broad utility in population genetics, disease-association studies, and diagnostic screening.

A large empirical assessment of sequence-resolved structural variants from 14,891 genomes across diverse global populations in the Genome Aggregation Database (gnomAD) provides a reference map for disease-association studies, population genetics, and diagnostic screening.

Genomically Aided Diagnosis of Severe Developmental Disorders

Our ability to make accurate and specific genetic diagnoses in individuals with severe developmental disorders has been transformed by data derived from genomic sequencing technologies. These data reveal both the patterns and rates of different mutational mechanisms and identify regions of the human genome with fewer mutations than would be expected. In outbred populations, the most common identifiable cause of severe developmental disorders is de novo mutation affecting the coding region in one of approximately 500 different genes, almost universally showing constraint. Simply combining the location of a de novo genomic event with its predicted consequence on the gene product gives significant diagnostic power. Our knowledge of the diversity of phenotypic consequences associated with comparable diagnostic genotypes at each locus is improving. Computationally useful phenotype data will improve diagnostic interpretation of ultrarare genetic variants and, in the long run, indicate which specific embryonic processes have been perturbed.

Infantile hypertrophic pyloric stenosis in patients with esophageal atresia

BACKGROUND

Patients born with esophageal atresia (EA) have a higher incidence of infantile hypertrophic pyloric stenosis (IHPS), suggestive of a relationship. A shared etiology makes sense from a developmental perspective as both affected structures are foregut derived. A genetic component has been described for both conditions as single entities and EA and IHPS are variable components in several monogenetic syndromes. We hypothesized that defects disturbing foregut morphogenesis are responsible for this combination of malformations.

METHODS

We investigated the genetic variation of 15 patients with both EA and IHPS with unaffected parents using exome sequencing and SNP array-based genotyping, and compared the results to mouse transcriptome data of the developing foregut.

RESULTS

We did not identify putatively deleterious de novo mutations or recessive variants. However, we detected rare inherited variants in EA or IHPS disease genes or in genes important in foregut morphogenesis, expressed at the proper developmental time-points. Two pathways were significantly enriched (p < 1 × 10^-5 ): proliferation and differentiation of smooth muscle cells and self-renewal of satellite cells.

CONCLUSIONS

None of our findings could fully explain the combination of abnormalities on its own, which makes complex inheritance the most plausible genetic explanation, most likely in combination with mechanical and/or environmental factors. As we did not find one defining monogenetic cause for the EA/IHPS phenotype, the impact of the corrective surgery could should be further investigated.

Characterization of Single Gene Copy Number Variants in Schizophrenia

BACKGROUND

Genetic studies of schizophrenia have implicated numerous risk loci including several copy number variants (CNVs) of large effect and hundreds of loci of small effect. In only a few cases has a specific gene been clearly identified. Rare CNVs affecting a single gene offer a potential avenue to discovering schizophrenia risk genes.

METHODS

CNVs were generated from exome sequencing of 4913 schizophrenia cases and 6188 control subjects from Sweden. We integrated two CNV calling methods (XHMM and ExomeDepth) to expand our set of single-gene CNVs and leveraged two different approaches for validating these variants (quantitative polymerase chain reaction and NanoString).

RESULTS

We found a significant excess of all rare CNVs (deletions: p = .0004, duplications: p = .0006) and single-gene CNVs (deletions: p = .04, duplications: p = .03) in schizophrenia cases compared with control subjects. An expanded set of CNVs generated from integrating multiple approaches showed a significant burden of deletions in 11 of 21 gene sets previously implicated in schizophrenia and across all genes in those sets (p = .008), although no tests survived correction. We performed an extensive validation of all deletions in the significant set of voltage-gated calcium channels among CNVs called from both exome sequencing and genotyping arrays. In total, 4 exonic, single-gene deletions were validated in schizophrenia cases and none in control subjects (p = .039), of which all were identified by exome sequencing.

CONCLUSIONS

These results point to the potential contribution of single-gene CNVs to schizophrenia, indicate that the utility of exome sequencing for CNV calling has yet to be maximized, and note that single-gene CNVs should be included in gene-focused studies using other classes of variation.

ATP7B variant penetrance explains differences between genetic and clinical prevalence estimates for Wilson disease

Wilson disease (WD) is a genetic disorder of copper metabolism caused by variants in the copper transporting P-type ATPase gene ATP7B. Estimates for WD population prevalence vary with 1 in 30,000 generally quoted. However, some genetic studies have reported much higher prevalence rates. The aim of this study was to estimate the population prevalence of WD and the pathogenicity/penetrance of WD variants by determining the frequency of ATP7B variants in a genomic sequence database. A catalogue of WD-associated ATP7B variants was constructed, and then, frequency information for these was extracted from the gnomAD data set. Pathogenicity of variants was assessed by (a) comparing gnomAD allele frequencies against the number of reports for variants in the WD literature and (b) using variant effect prediction algorithms. 231 WD-associated ATP7B variants were identified in the gnomAD data set, giving an initial estimated population prevalence of around 1 in 2400. After exclusion of WD-associated ATP7B variants with predicted low penetrance, the revised estimate showed a prevalence of around 1 in 20,000, with higher rates in the Asian and Ashkenazi Jewish populations. Reanalysis of other recent genetic studies using our penetrance criteria also predicted lower population prevalences for WD in the UK and France than had been reported. Our results suggest that differences in variant penetrance can explain the discrepancy between reported epidemiological and genetic prevalences of WD. They also highlight the challenge in defining penetrance when assigning causality to some ATP7B variants.

Revisiting the Population Genetics of Human Height

CONTEXT

Recent advances in genetics and genomics present unique opportunities for enhancing knowledge of human physiology and disease susceptibility. An outstanding example of these new insights may be seen in the study of human height, of which it has been estimated that approximately 80% is genetically determined. Over the past decade, large-scale population analyses have led to the identification of novel variation in genes and loci individually associated with changes in adult height of as much as 2 cm.

OBJECTIVE

To assess these same variants in the genomes of 213 158 individuals compiled by the Genome Aggregation Database (GnomAD) consortium, representing different population groups from around the world.

RESULTS

The majority of these height-changing alleles are substantially less prevalent in GnomAD than found previously in other cohorts, with 4 of 5 amino acid substitution variants with the largest impact on adult height being more frequent in the European population than in other groups.

CONCLUSIONS

A larger-scale analysis of individuals from diverse backgrounds will be necessary to ensure a full and accurate understanding of the genetic underpinnings of human height throughout the world, and additional studies will be needed to discern the biochemical and molecular mechanisms governing the physiological processes that explain how these variant proteins might selectively impact the biology of the growth plate. Broader understanding of the genetics of height also should set the stage for more comprehensive investigation into the causes of prevalent polygenic human diseases.

Sensitivity to gene dosage and gene expression affects genes with copy number variants observed among neuropsychiatric diseases

BACKGROUND

Copy number variants (CNVs) have been reported to be associated with diseases, traits, and evolution. However, it is hard to determine which gene should have priority as a target for further functional experiments if a CNV is rare or a singleton. In this study, we attempted to overcome this issue by using two approaches: by assessing the influences of gene dosage sensitivity and gene expression sensitivity. Dosage sensitive genes derived from two-round whole-genome duplication in previous studies. In addition, we proposed a cross-sectional omics approach that utilizes open data from GTEx to assess the effect of whole-genome CNVs on gene expression.

METHODS

Affymetrix Genome-Wide SNP Array 6.0 was used to detect CNVs by PennCNV and CNV Workshop. After quality controls for population stratification, family relationship and CNV detection, 287 patients with narcolepsy, 133 patients with essential hypersomnia, 380 patients with panic disorders, 164 patients with autism, 784 patients with Alzheimer disease and 1280 healthy individuals remained for the enrichment analysis.

RESULTS

Overall, significant enrichment of dosage sensitive genes was found across patients with narcolepsy, panic disorders and autism. Particularly, significant enrichment of dosage-sensitive genes in duplications was observed across all diseases except for Alzheimer disease. For deletions, less or no enrichment of dosage-sensitive genes with deletions was seen in the patients when compared to the healthy individuals. Interestingly, significant enrichments of genes with expression sensitivity in brain were observed in patients with panic disorder and autism. While duplications presented a higher burden, deletions did not cause significant differences when compared to the healthy individuals. When we assess the effect of sensitivity to genome dosage and gene expression at the same time, the highest ratio of enrichment was observed in the group including dosage-sensitive genes and genes with expression sensitivity only in brain. In addition, shared CNV regions among the five neuropsychiatric diseases were also investigated.

CONCLUSIONS

This study contributed the evidence that dosage-sensitive genes are associated with CNVs among neuropsychiatric diseases. In addition, we utilized open data from GTEx to assess the effect of whole-genome CNVs on gene expression. We also investigated shared CNV region among neuropsychiatric diseases.

Structural variation at the CYP2C locus: Characterization of deletion and duplication alleles

The human CYP2C locus harbors the polymorphic CYP2C18, CYP2C19, CYP2C9, and CYP2C8 genes, and of these, CYP2C19 and CYP2C9 are directly involved in the metabolism of ~15% of all medications. All variant CYP2C19 and CYP2C9 star (*) allele haplotypes currently cataloged by the Pharmacogene Variation (PharmVar) Consortium are defined by sequence variants. To determine if structural variation also occurs at the CYP2C locus, the 10q23.33 region was interrogated across deidentified clinical chromosomal microarray (CMA) data from 20,642 patients tested at two academic medical centers. Fourteen copy number variants that affected the coding region of CYP2C genes were detected in the clinical CMA cohorts, which ranged in size from 39.2 to 1,043.3 kb. Selected deletions and duplications were confirmed by MLPA or ddPCR. Analysis of the clinical CMA and an additional 78,839 cases from the Database of Genomic Variants (DGV) and ClinGen (total n = 99,481) indicated that the carrier frequency of a CYP2C structural variant is ~1 in 1,000, with ~1 in 2,000 being a CYP2C19 full gene or partial-gene deletion carrier, designated by PharmVar as CYP2C19*36 and *37, respectively. Although these structural variants are rare in the general population, their detection will likely improve metabolizer phenotype prediction when interrogated for research and/or clinical testing.

Copy number variation in DRC1 is the major cause of primary ciliary dyskinesia in the Japanese population

BACKGROUND

Primary ciliary dyskinesia (PCD) is a rare genetic disorder caused by functional impairment of cilia throughout the body. The involvement of copy number variation (CNV) in the development of PCD is largely unknown.

METHODS

We examined 93 Japanese patients with clinically suspected PCD from 84 unrelated families. CNV was examined either by exome sequencing of a PCD gene panel or by whole-exome sequencing (WES). The identified alterations were validated by PCR and Sanger sequencing. Nasal ciliary ultrastructure was examined by electron microscopy.

RESULTS

Analysis of CNV by the panel or WES revealed a biallelic deletion in the dynein regulatory complex subunit 1 (DRC1) gene in 21 patients, which accounted for 49% of the PCD patients in whom a disease-causing gene was found. Sanger sequencing of the PCR product revealed a 27,748-bp biallelic deletion including exons 1-4 of DRC1 with identical breakpoints in all 21 patients. The ciliary ultrastructure of the patients with this CNV showed axonemal disorganization and the loss or gain of central microtubules.

CONCLUSION

The deletion of DRC1 is the major cause of PCD in Japan and this alteration can cause various ciliary ultrastructural abnormalities.

Exome Sequencing in Individuals with Isolated Biliary Atresia

Biliary atresia (BA) is a severe pediatric liver disease resulting in necroinflammatory obliteration of the extrahepatic biliary tree. BA presents within the first few months of life as either an isolated finding or with additional syndromic features. The etiology of isolated BA is unknown, with evidence for infectious, environmental, and genetic risk factors described. However, to date, there are no definitive causal genes identified for isolated BA in humans, and the question of whether single gene defects play a major role remains open. We performed exome-sequencing in 101 North American patients of European descent with isolated BA (including 30 parent-child trios) and considered several experimental designs to identify potentially deleterious protein-altering variants that may be involved in the disease. In a case-only analysis, we did not identify genes with variants shared among more than two probands, and burden tests of rare variants using a case-case control design did not yield significant results. In the trio analysis of 30 simplex families (patient and parent trios), we identified 66 de novo variants in 66 genes including potentially deleterious variants in STIP1 and REV1. STIP1 is a co-chaperone for the heat-shock protein, HSP90, and has been shown to have diverse functions in yeast, flies and mammals, including stress-responses. REV1 is known to be a key player in DNA repair pathway and to interact with HSP90. In conclusion, our results do not support the hypothesis that a simple genetic model is responsible for the majority of cases of isolated BA. Our finding of de novo variants in genes linked to evolutionarily conserved stress responses (STIP1 and REV1) suggests that exploration of how genetic susceptibility and environmental exposure may interact to cause BA is warranted.