Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects.

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

Gene expression associations with body mass index in the Multi-Ethnic Study of Atherosclerosis.

BACKGROUND/OBJECTIVES: Obesity, defined as excessive fat accumulation that represents a health risk, is increasing in adults and children, reaching global epidemic proportions. Body mass index (BMI) correlates with body fat and future health risk, yet differs in prediction by fat distribution, across populations and by age. Nonetheless, few genetic studies of BMI have been conducted in ancestrally diverse populations. Gene expression association with BMI was assessed in the Multi-Ethnic Study of Atherosclerosis (MESA) in four self-identified race and ethnicity (SIRE) groups to identify genes associated with obesity. SUBJECTS/METHODS: RNA-sequencing was performed on 1096 MESA participants (37.8% white, 24.3% Hispanic, 28.4% African American, and 9.5% Chinese American) and linear models were used to assess the association of expression from each gene for its effect on BMI, adjusting for age, sex, sequencing center, study site, five expression and four genetic principal components in each self-identified race group. Sample-size-weighted meta-analysis was performed to identify genes with BMI-associated expression across ancestry groups. RESULTS: Within individual SIRE groups, there were zero to three genes whose expression is significantly (p < 1.97 × 10-6) associated with BMI. Across all groups, 45 genes were identified by meta-analysis whose expression was significantly associated with BMI, explaining 29.7% of BMI variation. The 45 genes are expressed in a variety of tissues and cell types and are enriched for obesity-related processes including erythrocyte function, oxygen binding and transport, and JAK-STAT signaling. CONCLUSIONS: We have identified genes whose expression is significantly associated with obesity in a multi-ethnic cohort. We have identified novel genes associated with BMI as well as confirmed previously identified genes from earlier genetic analyses. These novel genes and their biological pathways represent new targets for understanding the biology of obesity as well as new therapeutic intervention to reduce obesity and improve global public health.

Benchmarking association analyses of continuous exposures with RNA-seq in observational studies.

Large datasets of hundreds to thousands of individuals measuring RNA-seq in observational studies are becoming available. Many popular software packages for analysis of RNA-seq data were constructed to study differences in expression signatures in an experimental design with well-defined conditions (exposures). In contrast, observational studies may have varying levels of confounding transcript-exposure associations; further, exposure measures may vary from discrete (exposed, yes/no) to continuous (levels of exposure), with non-normal distributions of exposure. We compare popular software for gene expression-DESeq2, edgeR and limma-as well as linear regression-based analyses for studying the association of continuous exposures with RNA-seq. We developed a computation pipeline that includes transformation, filtering and generation of empirical null distribution of association P-values, and we apply the pipeline to compute empirical P-values with multiple testing correction. We employ a resampling approach that allows for assessment of false positive detection across methods, power comparison and the computation of quantile empirical P-values. The results suggest that linear regression methods are substantially faster with better control of false detections than other methods, even with the resampling method to compute empirical P-values. We provide the proposed pipeline with fast algorithms in an R package Olivia, and implemented it to study the associations of measures of sleep disordered breathing with RNA-seq in peripheral blood mononuclear cells in participants from the Multi-Ethnic Study of Atherosclerosis.

Use of neuromuscular blockade for neck dissection and association with iatrogenic nerve injury.

BACKGROUND: Cranial nerve injury is an uncommon but significant complication of neck dissection. We examined the association between the use of intraoperative neuromuscular blockade and iatrogenic cranial nerve injury during neck dissection. METHODS: This was a single-center, retrospective, electronic health record review. Study inclusion criteria stipulated patients > 18 years who had ≥ 2 neck lymphatic levels dissected for malignancy under general anesthesia with a surgery date between 2008 - 2018. Use of neuromuscular blockade during neck dissection was the primary independent variable. This was defined as any use of rocuronium, cisatracurium, or vecuronium upon anesthesia induction without reversal with sugammadex prior to surgical incision. Univariate tests were used to compare variables between those patients with, and those without, iatrogenic cranial nerve injury. Multivariable logistic regression determined predictors of cranial nerve injury and was performed incorporating Firth's estimation given low prevalence of the primary outcome. RESULTS: Our cohort consisted of 925 distinct neck dissections performed in 897 patients. Neuromuscular blockade was used during 285 (30.8%) neck dissections. Fourteen instances (1.5% of surgical cases) of nerve injury were identified. On univariate logistic regression, use of neuromuscular blockade was not associated with iatrogenic cranial nerve injury (OR: 1.73, 95% CI: 0.62 - 4.86, p = 0.30). There remained no significant association on multivariable logistic regression controlling for patient age, sex, weight, ASA class, paralytic dose, history of diabetes, stroke, coronary artery disease, carotid atherosclerosis, myocardial infarction, and cardiac arrythmia (OR: 1.87, 95% CI: 0.63 - 5.51, p = 0.26). CONCLUSIONS: In this study, use of neuromuscular blockade intraoperatively during neck dissection was not associated with increased rates of iatrogenic cranial nerve injury. While this investigation provides early support for safe use of neuromuscular blockade during neck dissection, future investigation with greater power remains necessary.

De Novo SOX4 Variants Cause a Neurodevelopmental Disease Associated with Mild Dysmorphism.

SOX4, together with SOX11 and SOX12, forms group C of SRY-related (SOX) transcription factors. They play key roles, often in redundancy, in multiple developmental pathways, including neurogenesis and skeletogenesis. De novo SOX11 heterozygous mutations have been shown to cause intellectual disability, growth deficiency, and dysmorphic features compatible with mild Coffin-Siris syndrome. Using trio-based exome sequencing, we here identify de novo SOX4 heterozygous missense variants in four children who share developmental delay, intellectual disability, and mild facial and digital morphological abnormalities. SOX4 is highly expressed in areas of active neurogenesis in human fetuses, and sox4 knockdown in Xenopus embryos diminishes brain and whole-body size. The SOX4 variants cluster in the highly conserved, SOX family-specific HMG domain, but each alters a different residue. In silico tools predict that each variant affects a distinct structural feature of this DNA-binding domain, and functional assays demonstrate that these SOX4 proteins carrying these variants are unable to bind DNA in vitro and transactivate SOX reporter genes in cultured cells. These variants are not found in the gnomAD database of individuals with presumably normal development, but 12 other SOX4 HMG-domain missense variants are recorded and all demonstrate partial to full activity in the reporter assay. Taken together, these findings point to specific SOX4 HMG-domain missense variants as the cause of a characteristic human neurodevelopmental disorder associated with mild facial and digital dysmorphism.

Mutations in the Epithelial Cadherin-p120-Catenin Complex Cause Mendelian Non-Syndromic Cleft Lip with or without Cleft Palate.

Non-syndromic cleft lip with or without cleft palate (NS-CL/P) is one of the most common human birth defects and is generally considered a complex trait. Despite numerous loci identified by genome-wide association studies, the effect sizes of common variants are relatively small, with much of the presumed genetic contribution remaining elusive. We report exome-sequencing results in 209 people from 72 multi-affected families with pedigree structures consistent with autosomal-dominant inheritance and variable penetrance. Herein, pathogenic variants are described in four genes encoding components of the p120-catenin complex (CTNND1, PLEKHA7, PLEKHA5) and an epithelial splicing regulator (ESRP2), in addition to the known CL/P-associated gene, CDH1, which encodes E-cadherin. The findings were also validated in a second cohort of 497 people with NS-CL/P, comprising small families and singletons with pathogenic variants in these genes identified in 14% of multi-affected families and 2% of the replication cohort of smaller families. Enriched expression of each gene/protein in human and mouse embryonic oro-palatal epithelia, demonstration of functional impact of CTNND1 and ESRP2 variants, and recapitulation of the CL/P spectrum in Ctnnd1 knockout mice support a causative role in CL/P pathogenesis. These data show that primary defects in regulators of epithelial cell adhesion are the most significant contributors to NS-CL/P identified to date and that inherited and de novo single gene variants explain a substantial proportion of NS-CL/P.

Common α-globin variants modify hematologic and other clinical phenotypes in sickle cell trait and disease.

Co-inheritance of α-thalassemia has a significant protective effect on the severity of complications of sickle cell disease (SCD), including stroke. However, little information exists on the association and interactions for the common African ancestral α-thalassemia mutation (-α3.7 deletion) and ß-globin traits (HbS trait [SCT] and HbC trait) on important clinical phenotypes such as red blood cell parameters, anemia, and chronic kidney disease (CKD). In a community-based cohort of 2,916 African Americans from the Jackson Heart Study, we confirmed the expected associations between SCT, HbC trait, and the -α3.7 deletion with lower mean corpuscular volume/mean corpuscular hemoglobin and higher red blood cell count and red cell distribution width. In addition to the recently recognized association of SCT with lower estimated glomerular filtration rate and glycated hemoglobin (HbA1c), we observed a novel association of the -α3.7 deletion with higher HbA1c levels. Co-inheritance of each additional copy of the -α3.7 deletion significantly lowered the risk of anemia and chronic kidney disease among individuals with SCT (P-interaction = 0.031 and 0.019, respectively). Furthermore, co-inheritance of a novel α-globin regulatory variant was associated with normalization of red cell parameters in individuals with the -α3.7 deletion and significantly negated the protective effect of α-thalassemia on stroke in 1,139 patients with sickle cell anemia from the Cooperative Study of Sickle Cell Disease (CSSCD) (P-interaction = 0.0049). Functional assays determined that rs11865131, located in the major alpha-globin enhancer MCS-R2, was the most likely causal variant. These findings suggest that common α- and ß-globin variants interact to influence hematologic and clinical phenotypes in African Americans, with potential implications for risk-stratification and counseling of individuals with SCD and SCT.

The contribution of de novo coding mutations to autism spectrum disorder.

Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

Defect-Directed Growth of Symmetrically Branched Metal Nanocrystals.

Branched plasmonic nanocrystals (NCs) have attracted much attention due to electric field enhancements at their tips. Seeded growth provides routes to NCs with defined branching patterns and, in turn, near-field distributions with defined symmetries. Here, a systematic analysis was undertaken in which seeds containing different distributions of planar defects were used to grow branched NCs in order to understand how their distributions direct the branching. Characterization of the products by multimode electron tomography and analysis of the NC morphologies at different overgrowth stages indicate that the branching patterns are directed by the seed defects, with the emergence of branches from the seed faces consistent with minimizing volumetric strain energy at the expense of surface energy. These results contrast with growth of branched NCs from single-crystalline seeds and provide a new platform for the synthesis of symmetrically branched plasmonic NCs.

Mutations in GDF11 and the extracellular antagonist, Follistatin, as a likely cause of Mendelian forms of orofacial clefting in humans.

Cleft lip with or without cleft palate (CL/P) is generally viewed as a complex trait with multiple genetic and environmental contributions. In 70% of cases, CL/P presents as an isolated feature and/or deemed nonsyndromic. In the remaining 30%, CL/P is associated with multisystem phenotypes or clinically recognizable syndromes, many with a monogenic basis. Here we report the identification, via exome sequencing, of likely pathogenic variants in two genes that encode interacting proteins previously only linked to orofacial clefting in mouse models. A variant in GDF11 (encoding growth differentiation factor 11), predicting a p.(Arg298Gln) substitution at the Furin protease cleavage site, was identified in one family that segregated with CL/P and both rib and vertebral hypersegmentation, mirroring that seen in Gdf11 knockout mice. In the second family in which CL/P was the only phenotype, a mutation in FST (encoding the GDF11 antagonist, Follistatin) was identified that is predicted to result in a p.(Cys56Tyr) substitution in the region that binds GDF11. Functional assays demonstrated a significant impact of the specific mutated amino acids on FST and GDF11 function and, together with embryonic expression data, provide strong evidence for the importance of GDF11 and Follistatin in the regulation of human orofacial development.

Exploring by whole exome sequencing patients with initial diagnosis of Rubinstein-Taybi syndrome: the interconnections of epigenetic machinery disorders.

Rubinstein-Taybi syndrome (RSTS) is an autosomal-dominant neurodevelopmental disease affecting 1:125,000 newborns characterized by intellectual disability, growth retardation, facial dysmorphisms and skeletal abnormalities. RSTS is caused by mutations in genes encoding for writers of the epigenetic machinery: CREBBP (~ 60%) or its homologous EP300 (~ 10%). No causative mutation is identified in up to 30% of patients. We performed whole-exome sequencing (WES) on eight RSTS-like individuals who had normal high-resolution array CGH testing and were CREBBP- and EP300-mutation -negative, to identify the molecular cause. In four cases, we identified putatively causal variants in three genes (ASXL1, KMT2D and KMT2A) encoding members of the epigenetic machinery known to be associated with the Bohring-Opitz, Kabuki and Wiedemann-Steiner syndromes. Each variant is novel, de novo, fulfills the ACMG criteria and is predicted to result in loss-of-function leading to haploinsufficiency of the epi-gene. In two of the remaining cases, homozygous/compound heterozygous variants in XYLT2 and PLCB4 genes, respectively, associated with spondyloocular and auriculocondylar 2 syndromes and in the latter an additional candidate variant in XRN2, a gene yet unrelated to any disease, were detected, but their pathogenicity remains uncertain. These results underscore the broad clinical spectrum of Mendelian disorders of the epigenetic apparatus and the high rate of WES disclosure of the genetic basis in cases which may pose a challenge for phenotype encompassing distinct syndromes. The overlapping features of distinct intellectual disability syndromes reflect common pathogenic molecular mechanisms affecting the complex regulation of balance between open and closed chromatin.

Genome Sequencing of Autism-Affected Families Reveals Disruption of Putative Noncoding Regulatory DNA.

We performed whole-genome sequencing (WGS) of 208 genomes from 53 families affected by simplex autism. For the majority of these families, no copy-number variant (CNV) or candidate de novo gene-disruptive single-nucleotide variant (SNV) had been detected by microarray or whole-exome sequencing (WES). We integrated multiple CNV and SNV analyses and extensive experimental validation to identify additional candidate mutations in eight families. We report that compared to control individuals, probands showed a significant (p = 0.03) enrichment of de novo and private disruptive mutations within fetal CNS DNase I hypersensitive sites (i.e., putative regulatory regions). This effect was only observed within 50 kb of genes that have been previously associated with autism risk, including genes where dosage sensitivity has already been established by recurrent disruptive de novo protein-coding mutations (ARID1B, SCN2A, NR3C2, PRKCA, and DSCAM). In addition, we provide evidence of gene-disruptive CNVs (in DISC1, WNT7A, RBFOX1, and MBD5), as well as smaller de novo CNVs and exon-specific SNVs missed by exome sequencing in neurodevelopmental genes (e.g., CANX, SAE1, and PIK3CA). Our results suggest that the detection of smaller, often multiple CNVs affecting putative regulatory elements might help explain additional risk of simplex autism.

Autosomal-Recessive Hearing Impairment Due to Rare Missense Variants within S1PR2.

The sphingosine-1-phosphate receptors (S1PRs) are a well-studied class of transmembrane G protein-coupled sphingolipid receptors that mediate multiple cellular processes. However, S1PRs have not been previously reported to be involved in the genetic etiology of human traits. S1PR2 lies within the autosomal-recessive nonsyndromic hearing impairment (ARNSHI) locus DFNB68 on 19p13.2. From exome sequence data we identified two pathogenic S1PR2 variants, c.323G>C (p.Arg108Pro) and c.419A>G (p.Tyr140Cys). Each of these variants co-segregates with congenital profound hearing impairment in consanguineous Pakistani families with maximum LOD scores of 6.4 for family DEM4154 and 3.3 for family PKDF1400. Neither S1PR2 missense variant was reported among ∼120,000 chromosomes in the Exome Aggregation Consortium database, in 76 unrelated Pakistani exomes, or in 720 Pakistani control chromosomes. Both DNA variants affect highly conserved residues of S1PR2 and are predicted to be damaging by multiple bioinformatics tools. Molecular modeling predicts that these variants affect binding of sphingosine-1-phosphate (p.Arg108Pro) and G protein docking (p.Tyr140Cys). In the previously reported S1pr2(-/-) mice, stria vascularis abnormalities, organ of Corti degeneration, and profound hearing loss were observed. Additionally, hair cell defects were seen in both knockout mice and morphant zebrafish. Family PKDF1400 presents with ARNSHI, which is consistent with the lack of gross malformations in S1pr2(-/-) mice, whereas family DEM4154 has lower limb malformations in addition to hearing loss. Our findings suggest the possibility of developing therapies against hair cell damage (e.g., from ototoxic drugs) through targeted stimulation of S1PR2.

Recurrent De Novo and Biallelic Variation of ATAD3A, Encoding a Mitochondrial Membrane Protein, Results in Distinct Neurological Syndromes.

ATPase family AAA-domain containing protein 3A (ATAD3A) is a nuclear-encoded mitochondrial membrane protein implicated in mitochondrial dynamics, nucleoid organization, protein translation, cell growth, and cholesterol metabolism. We identified a recurrent de novo ATAD3A c.1582C>T (p.Arg528Trp) variant by whole-exome sequencing (WES) in five unrelated individuals with a core phenotype of global developmental delay, hypotonia, optic atrophy, axonal neuropathy, and hypertrophic cardiomyopathy. We also describe two families with biallelic variants in ATAD3A, including a homozygous variant in two siblings, and biallelic ATAD3A deletions mediated by nonallelic homologous recombination (NAHR) between ATAD3A and gene family members ATAD3B and ATAD3C. Tissue-specific overexpression of borR534W, the Drosophila mutation homologous to the human c.1582C>T (p.Arg528Trp) variant, resulted in a dramatic decrease in mitochondrial content, aberrant mitochondrial morphology, and increased autophagy. Homozygous null bor larvae showed a significant decrease of mitochondria, while overexpression of borWT resulted in larger, elongated mitochondria. Finally, fibroblasts of an affected individual exhibited increased mitophagy. We conclude that the p.Arg528Trp variant functions through a dominant-negative mechanism that results in small mitochondria that trigger mitophagy, resulting in a reduction in mitochondrial content. ATAD3A variation represents an additional link between mitochondrial dynamics and recognizable neurological syndromes, as seen with MFN2, OPA1, DNM1L, and STAT2 mutations.

Analysis of tumor-infiltrating CD103 resident memory T-cell content in recurrent laryngeal squamous cell carcinoma.

BACKGROUND: Recurrent laryngeal squamous cell carcinomas (LSCCs) are associated with poor outcomes, without reliable biomarkers to identify patients who may benefit from adjuvant therapies. Given the emergence of tumor-infiltrating lymphocytes (TIL) as a biomarker in head and neck squamous cell carcinoma, we generated predictive models to understand the utility of CD4+, CD8+ and/or CD103+ TIL status in patients with advanced LSCC. METHODS: Tissue microarrays were constructed from salvage laryngectomy specimens of 183 patients with recurrent/persistent LSCC and independently stained for CD4+, CD8+, and CD103+ TIL content. Cox proportional hazards regression analysis was employed to assess combinations of CD4+, CD8+, and CD103+ TIL levels for prediction of overall survival (OS), disease-specific survival (DSS), and disease-free survival (DFS) in patients with recurrent/persistent LSCC. RESULTS: High tumor CD103+ TIL content was associated with significantly improved OS, DSS, and DFS and was a stronger predictor of survival in recurrent/persistent LSCC than either high CD8+ or CD4+ TIL content. On multivariate analysis, an "immune-rich" phenotype, in which tumors were enriched for both CD103+ and CD4+ TILs, conferred a survival benefit (OS hazard ratio: 0.28, p = 0.0014; DSS hazard ratio: 0.09, p = 0.0015; DFS hazard ratio: 0.18, p = 0.0018) in recurrent/persistent LSCC. CONCLUSIONS: An immune profile driven by CD103+ TIL content, alone and in combination with CD4+ TIL content, is a prognostic biomarker of survival in patients with recurrent/persistent LSCC. Predictive models described herein may thus prove valuable in prognostic stratification and lead to personalized treatment paradigms for this patient population.

Elective Paratracheal Lymph Node Dissection in Salvage Laryngectomy.

BACKGROUND: Indications for and efficacy of paratracheal nodal dissection (PTND) in patients undergoing laryngectomy (salvage) for persistent or recurrent laryngeal squamous cell carcinoma are not well-defined. METHODS: A retrospective cohort study was performed for patients undergoing salvage laryngectomy with clinically and radiographically negative neck disease between 1998 and 2015 (n = 210). Univariate and multivariate Cox regression analyses were performed. RESULTS: PTND was performed on 77/210 patients (36%). The PTND cohort had a greater proportion of advanced T classification (rT3/rT4) tumors (78%) than subjects without PTND (55%; p = 0.001). There was a 14% rate of occult nodal metastases in the paratracheal basin; of these, 55% did not have pathologic lateral neck disease. Multivariate analysis controlling for tumor site, tumor stage, and pathologic lateral neck disease demonstrated that PTND was associated with improved overall survival [OS] (p = 0.03; hazard ratio [HR] 0.60, 95% confidence interval [CI] 0.38-0.96), disease-free survival [DFS] (p = 0.03; HR 0.55, 95% CI 0.31-0.96), and distant DFS survival (p = 0.01; HR 0.29, 95% CI 0.11-0.77). The rate of hypocalcemia did not differ between subjects who underwent bilateral PTND, unilateral PTND, or no PTND (p = 0.19 at discharge, p = 0.17 at last follow-up). CONCLUSIONS: PTND at the time of salvage laryngectomy was more common in patients with rT3/rT4 tumors and was associated with improved OS and DFS, with no effect on hypocalcemia. In patients undergoing PTND, the finding of occult paratracheal metastases was often independent of lateral neck metastases.

A statin-dependent QTL for GATM expression is associated with statin-induced myopathy.

Statins are prescribed widely to lower plasma low-density lipoprotein (LDL) concentrations and cardiovascular disease risk and have been shown to have beneficial effects in a broad range of patients. However, statins are associated with an increased risk, albeit small, of clinical myopathy and type 2 diabetes. Despite evidence for substantial genetic influence on LDL concentrations, pharmacogenomic trials have failed to identify genetic variations with large effects on either statin efficacy or toxicity, and have produced little information regarding mechanisms that modulate statin response. Here we identify a downstream target of statin treatment by screening for the effects of in vitro statin exposure on genetic associations with gene expression levels in lymphoblastoid cell lines derived from 480 participants of a clinical trial of simvastatin treatment. This analysis identified six expression quantitative trait loci (eQTLs) that interacted with simvastatin exposure, including rs9806699, a cis-eQTL for the gene glycine amidinotransferase (GATM) that encodes the rate-limiting enzyme in creatine synthesis. We found this locus to be associated with incidence of statin-induced myotoxicity in two separate populations (meta-analysis odds ratio = 0.60). Furthermore, we found that GATM knockdown in hepatocyte-derived cell lines attenuated transcriptional response to sterol depletion, demonstrating that GATM may act as a functional link between statin-mediated lowering of cholesterol and susceptibility to statin-induced myopathy.

The Genetic Basis of Mendelian Phenotypes: Discoveries, Challenges, and Opportunities.

Discovering the genetic basis of a Mendelian phenotype establishes a causal link between genotype and phenotype, making possible carrier and population screening and direct diagnosis. Such discoveries also contribute to our knowledge of gene function, gene regulation, development, and biological mechanisms that can be used for developing new therapeutics. As of February 2015, 2,937 genes underlying 4,163 Mendelian phenotypes have been discovered, but the genes underlying ∼50% (i.e., 3,152) of all known Mendelian phenotypes are still unknown, and many more Mendelian conditions have yet to be recognized. This is a formidable gap in biomedical knowledge. Accordingly, in December 2011, the NIH established the Centers for Mendelian Genomics (CMGs) to provide the collaborative framework and infrastructure necessary for undertaking large-scale whole-exome sequencing and discovery of the genetic variants responsible for Mendelian phenotypes. In partnership with 529 investigators from 261 institutions in 36 countries, the CMGs assessed 18,863 samples from 8,838 families representing 579 known and 470 novel Mendelian phenotypes as of January 2015. This collaborative effort has identified 956 genes, including 375 not previously associated with human health, that underlie a Mendelian phenotype. These results provide insight into study design and analytical strategies, identify novel mechanisms of disease, and reveal the extensive clinical variability of Mendelian phenotypes. Discovering the gene underlying every Mendelian phenotype will require tackling challenges such as worldwide ascertainment and phenotypic characterization of families affected by Mendelian conditions, improvement in sequencing and analytical techniques, and pervasive sharing of phenotypic and genomic data among researchers, clinicians, and families.

De novo mutations in NALCN cause a syndrome characterized by congenital contractures of the limbs and face, hypotonia, and developmental delay.

Freeman-Sheldon syndrome, or distal arthrogryposis type 2A (DA2A), is an autosomal-dominant condition caused by mutations in MYH3 and characterized by multiple congenital contractures of the face and limbs and normal cognitive development. We identified a subset of five individuals who had been putatively diagnosed with "DA2A with severe neurological abnormalities" and for whom congenital contractures of the limbs and face, hypotonia, and global developmental delay had resulted in early death in three cases; this is a unique condition that we now refer to as CLIFAHDD syndrome. Exome sequencing identified missense mutations in the sodium leak channel, non-selective (NALCN) in four families affected by CLIFAHDD syndrome. We used molecular-inversion probes to screen for NALCN in a cohort of 202 distal arthrogryposis (DA)-affected individuals as well as concurrent exome sequencing of six other DA-affected individuals, thus revealing NALCN mutations in ten additional families with "atypical" forms of DA. All 14 mutations were missense variants predicted to alter amino acid residues in or near the S5 and S6 pore-forming segments of NALCN, highlighting the functional importance of these segments. In vitro functional studies demonstrated that NALCN alterations nearly abolished the expression of wild-type NALCN, suggesting that alterations that cause CLIFAHDD syndrome have a dominant-negative effect. In contrast, homozygosity for mutations in other regions of NALCN has been reported in three families affected by an autosomal-recessive condition characterized mainly by hypotonia and severe intellectual disability. Accordingly, mutations in NALCN can cause either a recessive or dominant condition characterized by varied though overlapping phenotypic features, perhaps based on the type of mutation and affected protein domain(s).