Rare variants in CAPN2 increase risk for isolated hypoplastic left heart syndrome.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart defect (CHD) characterized by hypoplasia of the left ventricle and aorta along with stenosis or atresia of the aortic and mitral valves. HLHS represents only ∼4%-8% of all CHDs but accounts for ∼25% of deaths. HLHS is an isolated defect (i.e., iHLHS) in 70% of families, the vast majority of which are simplex. Despite intense investigation, the genetic basis of iHLHS remains largely unknown. We performed exome sequencing on 331 families with iHLHS aggregated from four independent cohorts. A Mendelian-model-based analysis demonstrated that iHLHS was not due to single, large-effect alleles in genes previously reported to underlie iHLHS or CHD in >90% of families in this cohort. Gene-based association testing identified increased risk for iHLHS associated with variation in CAPN2 (p = 1.8 × 10-5), encoding a protein involved in functional adhesion. Functional validation studies in a vertebrate animal model (Xenopus laevis) confirmed CAPN2 is essential for cardiac ventricle morphogenesis and that in vivo loss of calpain function causes hypoplastic ventricle phenotypes and suggest that human CAPN2707C>T and CAPN21112C>T variants, each found in multiple individuals with iHLHS, are hypomorphic alleles. Collectively, our findings show that iHLHS is typically not a Mendelian condition, demonstrate that CAPN2 variants increase risk of iHLHS, and identify a novel pathway involved in HLHS pathogenesis.

Phenotypic and mutational spectrum of ROR2-related Robinow syndrome.

Robinow syndrome is characterized by a triad of craniofacial dysmorphisms, disproportionate-limb short stature, and genital hypoplasia. A significant degree of phenotypic variability seems to correlate with different genes/loci. Disturbances of the noncanonical WNT-pathway have been identified as the main cause of the syndrome. Biallelic variants in ROR2 cause an autosomal recessive form of the syndrome with distinctive skeletal findings. Twenty-two patients with a clinical diagnosis of autosomal recessive Robinow syndrome were screened for variants in ROR2 using multiple molecular approaches. We identified 25 putatively pathogenic ROR2 variants, 16 novel, including single nucleotide variants and exonic deletions. Detailed phenotypic analyses revealed that all subjects presented with a prominent forehead, hypertelorism, short nose, abnormality of the nasal tip, brachydactyly, mesomelic limb shortening, short stature, and genital hypoplasia in male patients. A total of 19 clinical features were present in more than 75% of the subjects, thus pointing to an overall uniformity of the phenotype. Disease-causing variants in ROR2, contribute to a clinically recognizable autosomal recessive trait phenotype with multiple skeletal defects. A comprehensive quantitative clinical evaluation of this cohort delineated the phenotypic spectrum of ROR2-related Robinow syndrome. The identification of exonic deletion variant alleles further supports the contention of a loss-of-function mechanism in the etiology of the syndrome.

Further delineation of van den Ende-Gupta syndrome: Genetic heterogeneity and overlap with congenital heart defects and skeletal malformations syndrome.

Van den Ende-Gupta syndrome (VDEGS) is a rare autosomal recessive condition characterized by distinctive facial and skeletal features, and in most affected persons, by biallelic pathogenic variants in SCARF2. We review the type and frequency of the clinical features in 36 reported individuals with features of VDEGS, 15 (42%) of whom had known pathogenic variants in SCARF2, 6 (16%) with negative SCARF2 testing, and 15 (42%) not tested. We also report three new individuals with pathogenic variants in SCARF2 and clinical features of VDEGS. Of the six persons without known pathogenic variants in SCARF2, three remain unsolved despite extensive genetic testing. Three were found to have pathogenic ABL1 variants using whole exome sequencing (WES) or whole genome sequencing (WGS). Their phenotype was consistent with the congenital heart disease and skeletal malformations syndrome (CHDSKM), which has been associated with ABL1 variants. Of the three unsolved cases, two were brothers who underwent WGS and targeted long-range sequencing of both SCARF2 and ABL1, and the third person who underwent WES and RNA sequencing for SCARF2. Because these affected individuals with classical features of VDEGS lacked a detectable pathogenic SCARF2 variant, genetic heterogeneity is likely. Our study shows the importance of performing genetic testing on individuals with the VDEGS "phenotype," either as a targeted gene analysis (SCARF2, ABL1) or WES/WGS. Additionally, individuals with the combination of arachnodactyly and blepharophimosis should undergo echocardiography while awaiting results of molecular testing due to the overlapping physical features of VDEGS and CHDSKM.

Paralog Studies Augment Gene Discovery: DDX and DHX Genes.

Members of a paralogous gene family in which variation in one gene is known to cause disease are eight times more likely to also be associated with human disease. Recent studies have elucidated DHX30 and DDX3X as genes for which pathogenic variant alleles are involved in neurodevelopmental disorders. We hypothesized that variants in paralogous genes encoding members of the DExD/H-box RNA helicase superfamily might also underlie developmental delay and/or intellectual disability (DD and/or ID) disease phenotypes. Here we describe 15 unrelated individuals who have DD and/or ID, central nervous system (CNS) dysfunction, vertebral anomalies, and dysmorphic features and were found to have probably damaging variants in DExD/H-box RNA helicase genes. In addition, these individuals exhibit a variety of other tissue and organ system involvement including ocular, outer ear, hearing, cardiac, and kidney tissues. Five individuals with homozygous (one), compound-heterozygous (two), or de novo (two) missense variants in DHX37 were identified by exome sequencing. We identified ten total individuals with missense variants in three other DDX/DHX paralogs: DHX16 (four individuals), DDX54 (three individuals), and DHX34 (three individuals). Most identified variants are rare, predicted to be damaging, and occur at conserved amino acid residues. Taken together, these 15 individuals implicate the DExD/H-box helicases in both dominantly and recessively inherited neurodevelopmental phenotypes and highlight the potential for more than one disease mechanism underlying these disorders.

Exome Sequencing of a Primary Ovarian Insufficiency Cohort Reveals Common Molecular Etiologies for a Spectrum of Disease.

CONTEXT: Primary ovarian insufficiency (POI) encompasses a spectrum of premature menopause, including both primary and secondary amenorrhea. For 75% to 90% of individuals with hypergonadotropic hypogonadism presenting as POI, the molecular etiology is unknown. Common etiologies include chromosomal abnormalities, environmental factors, and congenital disorders affecting ovarian development and function, as well as syndromic and nonsyndromic single gene disorders suggesting POI represents a complex trait. OBJECTIVE: To characterize the contribution of known disease genes to POI and identify molecular etiologies and biological underpinnings of POI. DESIGN, SETTING, AND PARTICIPANTS: We applied exome sequencing (ES) and family-based genomics to 42 affected female individuals from 36 unrelated Turkish families, including 31 with reported parental consanguinity. RESULTS: This analysis identified likely damaging, potentially contributing variants and molecular diagnoses in 16 families (44%), including 11 families with likely damaging variants in known genes and five families with predicted deleterious variants in disease genes (IGSF10, MND1, MRPS22, and SOHLH1) not previously associated with POI. Of the 16 families, 2 (13%) had evidence for potentially pathogenic variants at more than one locus. Absence of heterozygosity consistent with identity-by-descent mediated recessive disease burden contributes to molecular diagnosis in 15 of 16 (94%) families. GeneMatcher allowed identification of additional families from diverse genetic backgrounds. CONCLUSIONS: ES analysis of a POI cohort further characterized locus heterogeneity, reaffirmed the association of genes integral to meiotic recombination, demonstrated the likely contribution of genes involved in hypothalamic development, and documented multilocus pathogenic variation suggesting the potential for oligogenic inheritance contributing to the development of POI.

Insights into genetics, human biology and disease gleaned from family based genomic studies.

Identifying genes and variants contributing to rare disease phenotypes and Mendelian conditions informs biology and medicine, yet potential phenotypic consequences for variation of >75% of the ~20,000 annotated genes in the human genome are lacking. Technical advances to assess rare variation genome-wide, particularly exome sequencing (ES), enabled establishment in the United States of the National Institutes of Health (NIH)-supported Centers for Mendelian Genomics (CMGs) and have facilitated collaborative studies resulting in novel "disease gene" discoveries. Pedigree-based genomic studies and rare variant analyses in families with suspected Mendelian conditions have led to the elucidation of hundreds of novel disease genes and highlighted the impact of de novo mutational events, somatic variation underlying nononcologic traits, incompletely penetrant alleles, phenotypes with high locus heterogeneity, and multilocus pathogenic variation. Herein, we highlight CMG collaborative discoveries that have contributed to understanding both rare and common diseases and discuss opportunities for future discovery in single-locus Mendelian disorder genomics. Phenotypic annotation of all human genes; development of bioinformatic tools and analytic methods; exploration of non-Mendelian modes of inheritance including reduced penetrance, multilocus variation, and oligogenic inheritance; construction of allelic series at a locus; enhanced data sharing worldwide; and integration with clinical genomics are explored. Realizing the full contribution of rare disease research to functional annotation of the human genome, and further illuminating human biology and health, will lay the foundation for the Precision Medicine Initiative.

Identifying Genes Whose Mutant Transcripts Cause Dominant Disease Traits by Potential Gain-of-Function Alleles.

Premature termination codon (PTC)-bearing transcripts are often degraded by nonsense-mediated decay (NMD) resulting in loss-of-function (LoF) alleles. However, not all PTCs result in LoF mutations, i.e., some such transcripts escape NMD and are translated to truncated peptide products that result in disease due to gain-of-function (GoF) effects. Since the location of the PTC is a major factor determining transcript fate, we hypothesized that depletion of protein-truncating variants (PTVs) within the gene region predicted to escape NMD in control databases could provide a rank for genic susceptibility for disease through GoF versus LoF. We developed an NMD escape intolerance score to rank genes based on the depletion of PTVs that would render them able to escape NMD using the Atherosclerosis Risk in Communities Study (ARIC) and the Exome Aggregation Consortium (ExAC) control databases, which was further used to screen the Baylor-Center for Mendelian Genomics disease database. This analysis revealed 1,996 genes significantly depleted for PTVs that are predicted to escape from NMD, i.e., PTVesc; further studies provided evidence that revealed a subset as candidate genes underlying Mendelian phenotypes. Importantly, these genes have characteristically low pLI scores, which can cause them to be overlooked as candidates for dominant diseases. Collectively, we demonstrate that this NMD escape intolerance score is an effective and efficient tool for gene discovery in Mendelian diseases due to production of truncated or altered proteins. More importantly, we provide a complementary analytical tool to aid identification of genes associated with dominant traits through a mechanism distinct from LoF.

CHRNA7 copy number gains are enriched in adolescents with major depressive and anxiety disorders.

OBJECTIVE: Neuronal nicotinic acetylcholine receptors (nAChRs), specifically the α7 nAChR encoded by the gene CHRNA7, have been implicated in behavior regulation in animal models. In humans, copy number variants (CNVs) of CHRNA7 are found in a range of neuropsychiatric disorders, including mood and anxiety disorders. Here, we aimed to determine the prevalence of CHRNA7 CNVs among adolescents and young adults with major depressive disorder (MDD) and anxiety disorders. METHODS: Twelve to 21 year-old participants with MDD and/or anxiety disorders (34% males, mean ±â€¯std age: 18.9 ±â€¯1.8 years) were assessed for CHRNA7 copy number state using droplet digital PCR (ddPCR) and genomic quantitative PCR (qPCR). Demographic, anthropometric, and clinical data, including the Beck Anxiety Index (BAI), Beck Depression Inventory (BDI), and the Inventory of Depressive Symptoms (IDS) were collected and compared across individuals with and without a CHRNA7 CNV. RESULTS: Of 205 individuals, five (2.4%) were found to carry a CHRNA7 gain, significantly higher than the general population. No CHRNA7 deletions were identified. Clinically, the individuals carrying CHRNA7 duplications did not differ significantly from copy neutral individuals with MDD and/or anxiety disorders. CONCLUSIONS: CHRNA7 gains are relatively prevalent among young individuals with MDD and anxiety disorders (odds ratio = 4.032) without apparent distinguishing clinical features. Future studies should examine the therapeutic potential of α7 nAChR targeting drugs to ameliorate depressive and anxiety disorders.

WNT Signaling Perturbations Underlie the Genetic Heterogeneity of Robinow Syndrome.

Locus heterogeneity characterizes a variety of skeletal dysplasias often due to interacting or overlapping signaling pathways. Robinow syndrome is a skeletal disorder historically refractory to molecular diagnosis, potentially stemming from substantial genetic heterogeneity. All current known pathogenic variants reside in genes within the noncanonical Wnt signaling pathway including ROR2, WNT5A, and more recently, DVL1 and DVL3. However, ∼70% of autosomal-dominant Robinow syndrome cases remain molecularly unsolved. To investigate this missing heritability, we recruited 21 families with at least one family member clinically diagnosed with Robinow or Robinow-like phenotypes and performed genetic and genomic studies. In total, four families with variants in FZD2 were identified as well as three individuals from two families with biallelic variants in NXN that co-segregate with the phenotype. Importantly, both FZD2 and NXN are relevant protein partners in the WNT5A interactome, supporting their role in skeletal development. In addition to confirming that clustered -1 frameshifting variants in DVL1 and DVL3 are the main contributors to dominant Robinow syndrome, we also found likely pathogenic variants in candidate genes GPC4 and RAC3, both linked to the Wnt signaling pathway. These data support an initial hypothesis that Robinow syndrome results from perturbation of the Wnt/PCP pathway, suggest specific relevant domains of the proteins involved, and reveal key contributors in this signaling cascade during human embryonic development. Contrary to the view that non-allelic genetic heterogeneity hampers gene discovery, this study demonstrates the utility of rare disease genomic studies to parse gene function in human developmental pathways.

CHRNA7 Deletions are Enriched in Risperidone-Treated Children and Adolescents.

OBJECTIVE: Aggression is among the most common indications for referral to child and adolescent mental health services and is often challenging to treat. Understanding the biological underpinnings of aggression could help optimize treatment efficacy. Neuronal nicotinic acetylcholine receptors (nAChRs), specifically the α7 nAChR, encoded by the gene CHRNA7, have been implicated in aggressive behaviors in animal models as well as humans. Copy number variants (CNVs) of CHRNA7 are found in individuals with neuropsychiatric disorders, often with comorbid aggression. In this study, we aimed to determine the prevalence of CHRNA7 CNVs among individuals treated with risperidone, predominantly for irritability and aggression. METHODS: Risperidone-treated children and adolescents were assessed for CHRNA7 copy number state using droplet digital PCR and genomic quantitative PCR. Demographic, anthropometric, and clinical data, including the Child Behavior Checklist (CBCL), were collected and compared across individuals with and without the CHRNA7 deletion. RESULTS: Of 218 individuals (90% males, mean age: 12.3 ± 2.3 years), 7 (3.2%) were found to carry a CHRNA7 deletion and one proband carried a CHRNA7 duplication (0.46%). T-scores for rule breaking, aggression, and externalizing behavior factors of the CBCL were higher in the deletion group, despite taking 58% higher dose of risperidone. CONCLUSIONS: CHRNA7 loss may contribute to a phenotype of severe aggression. Given the high prevalence of the deletion among risperidone-treated youth, future studies should examine the therapeutic potential of α7 nAChR-targeting drugs to target aggression associated with CHRNA7 deletions.

REST Final-Exon-Truncating Mutations Cause Hereditary Gingival Fibromatosis.

Hereditary gingival fibromatosis (HGF) is the most common genetic form of gingival fibromatosis that develops as a slowly progressive, benign, localized or generalized enlargement of keratinized gingiva. HGF is a genetically heterogeneous disorder and can be transmitted either as an autosomal-dominant or autosomal-recessive trait or appear sporadically. To date, four loci (2p22.1, 2p23.3-p22.3, 5q13-q22, and 11p15) have been mapped to autosomes and one gene (SOS1) has been associated with the HGF trait observed to segregate in a dominant inheritance pattern. Here we report 11 individuals with HGF from three unrelated families. Whole-exome sequencing (WES) revealed three different truncating mutations including two frameshifts and one nonsense variant in RE1-silencing transcription factor (REST) in the probands from all families and further genetic and genomic analyses confirmed the WES-identified findings. REST is a transcriptional repressor that is expressed throughout the body; it has different roles in different cellular contexts, such as oncogenic and tumor-suppressor functions and hematopoietic and cardiac differentiation. Here we show the consequences of germline final-exon-truncating mutations in REST for organismal development and the association with the HGF phenotype.

Dual molecular diagnosis contributes to atypical Prader-Willi phenotype in monozygotic twins.

We describe monozygotic twin girls with genetic variation at two separate loci resulting in a blended phenotype of Prader-Willi syndrome and Pitt-Hopkins syndrome. These girls were diagnosed in early infancy with Prader-Willi syndrome, but developed an atypical phenotype, with apparent intellectual deficiency and lack of obesity. Array-comparative genomic hybridization confirmed a de novo paternal deletion of the 15q11.2q13 region and exome sequencing identified a second mutational event in both girls, which was a novel variant c.145+1G>A affecting a TCF4 canonical splicing site inherited from the mosaic mother. RNA studies showed that the variant abolished the donor splicing site, which was accompanied by activation of an alternative non-canonical splicing-site which then predicts a premature stop codon in the following exon. Clinical re-evaluation of the twins indicated that both variants are likely contributing to the more severe phenotypic presentation. Our data show that atypical clinical presentations may actually be the expression of blended clinical phenotypes arising from independent pathogenic events at two loci.

Identification of a RAI1-associated disease network through integration of exome sequencing, transcriptomics, and 3D genomics.

BACKGROUND: Smith-Magenis syndrome (SMS) is a developmental disability/multiple congenital anomaly disorder resulting from haploinsufficiency of RAI1. It is characterized by distinctive facial features, brachydactyly, sleep disturbances, and stereotypic behaviors. METHODS: We investigated a cohort of 15 individuals with a clinical suspicion of SMS who showed neither deletion in the SMS critical region nor damaging variants in RAI1 using whole exome sequencing. A combination of network analysis (co-expression and biomedical text mining), transcriptomics, and circularized chromatin conformation capture (4C-seq) was applied to verify whether modified genes are part of the same disease network as known SMS-causing genes. RESULTS: Potentially deleterious variants were identified in nine of these individuals using whole-exome sequencing. Eight of these changes affect KMT2D, ZEB2, MAP2K2, GLDC, CASK, MECP2, KDM5C, and POGZ, known to be associated with Kabuki syndrome 1, Mowat-Wilson syndrome, cardiofaciocutaneous syndrome, glycine encephalopathy, mental retardation and microcephaly with pontine and cerebellar hypoplasia, X-linked mental retardation 13, X-linked mental retardation Claes-Jensen type, and White-Sutton syndrome, respectively. The ninth individual carries a de novo variant in JAKMIP1, a regulator of neuronal translation that was recently found deleted in a patient with autism spectrum disorder. Analyses of co-expression and biomedical text mining suggest that these pathologies and SMS are part of the same disease network. Further support for this hypothesis was obtained from transcriptome profiling that showed that the expression levels of both Zeb2 and Map2k2 are perturbed in Rai1 -/- mice. As an orthogonal approach to potentially contributory disease gene variants, we used chromatin conformation capture to reveal chromatin contacts between RAI1 and the loci flanking ZEB2 and GLDC, as well as between RAI1 and human orthologs of the genes that show perturbed expression in our Rai1 -/- mouse model. CONCLUSIONS: These holistic studies of RAI1 and its interactions allow insights into SMS and other disorders associated with intellectual disability and behavioral abnormalities. Our findings support a pan-genomic approach to the molecular diagnosis of a distinctive disorder.

De novo missense variants in PPP1CB are associated with intellectual disability and congenital heart disease.

Intellectual disabilities are genetically heterogeneous and can be associated with congenital anomalies. Using whole-exome sequencing (WES), we identified five different de novo missense variants in the protein phosphatase-1 catalytic subunit beta (PPP1CB) gene in eight unrelated individuals who share an overlapping phenotype of dysmorphic features, macrocephaly, developmental delay or intellectual disability (ID), congenital heart disease, short stature, and skeletal and connective tissue abnormalities. Protein phosphatase-1 (PP1) is a serine/threonine-specific protein phosphatase involved in the dephosphorylation of a variety of proteins. The PPP1CB gene encodes a PP1 subunit that regulates the level of protein phosphorylation. All five altered amino acids we observed are highly conserved among the PP1 subunit family, and all are predicted to disrupt PP1 subunit binding and impair dephosphorylation. Our data suggest that our heterozygous de novo PPP1CB pathogenic variants are associated with syndromic intellectual disability.

DVL3 Alleles Resulting in a -1 Frameshift of the Last Exon Mediate Autosomal-Dominant Robinow Syndrome.

Robinow syndrome is a rare congenital disorder characterized by mesomelic limb shortening, genital hypoplasia, and distinctive facial features. Recent reports have identified, in individuals with dominant Robinow syndrome, a specific type of variant characterized by being uniformly located in the penultimate exon of DVL1 and resulting in a -1 frameshift allele with a premature termination codon that escapes nonsense-mediated decay. Here, we studied a cohort of individuals who had been clinically diagnosed with Robinow syndrome but who had not received a molecular diagnosis from variant studies of DVL1, WNT5A, and ROR2. Because of the uniform location of frameshift variants in DVL1-mediated Robinow syndrome and the functional redundancy of DVL1, DVL2, and DVL3, we elected to pursue direct Sanger sequencing of the penultimate exon of DVL1 and its paralogs DVL2 and DVL3 to search for potential disease-associated variants. Remarkably, targeted sequencing identified five unrelated individuals harboring heterozygous, de novo frameshift variants in DVL3, including two splice acceptor mutations and three 1 bp deletions. Similar to the variants observed in DVL1-mediated Robinow syndrome, all variants in DVL3 result in a -1 frameshift, indicating that these highly specific alterations might be a common cause of dominant Robinow syndrome. Here, we review the current knowledge of these peculiar variant alleles in DVL1- and DVL3-mediated Robinow syndrome and further elucidate the phenotypic features present in subjects with DVL1 and DVL3 frameshift mutations.

Hutterite-type cataract maps to chromosome 6p21.32-p21.31, cosegregates with a homozygous mutation in LEMD2, and is associated with sudden cardiac death.

BACKGROUND: Juvenile-onset cataracts are known among the Hutterites of North America. Despite being identified over 30 years ago, this autosomal recessive condition has not been mapped, and the disease gene is unknown. METHODS: We performed whole exome sequencing of three Hutterite-type cataract trios and follow-up genotyping and mapping in four extended kindreds. RESULTS: Trio exomes enabled genome-wide autozygosity mapping, which localized the disease gene to a 9.5-Mb region on chromosome 6p. This region contained two candidate variants, LEMD2 c.T38G and MUC21 c.665delC. Extended pedigrees recruited for variant genotyping revealed multiple additional relatives with juvenile-onset cataract, as well as six deceased relatives with both cataracts and sudden cardiac death. The candidate variants were genotyped in 84 family members, including 17 with cataracts; only the variant in LEMD2 cosegregated with cataracts (LOD = 9.62). SNP-based fine mapping within the 9.5 Mb linked region supported this finding by refining the cataract locus to a 0.5- to 2.9-Mb subregion (6p21.32-p21.31) containing LEMD2 but not MUC21. LEMD2 is expressed in mouse and human lenses and encodes a LEM domain-containing protein; the c.T38G missense mutation is predicted to mutate a highly conserved residue within this domain (p.Leu13Arg). CONCLUSION: We performed a genetic and genomic study of Hutterite-type cataract and found evidence for an association of this phenotype with sudden cardiac death. Using combined genetic and genomic approaches, we mapped cataracts to a small portion of chromosome 6 and propose that they result from a homozygous missense mutation in LEMD2.