Analysis of missense variants in the human genome reveals widespread gene-specific clustering and improves prediction of pathogenicity.

We used a machine learning approach to analyze the within-gene distribution of missense variants observed in hereditary conditions and cancer. When applied to 840 genes from the ClinVar database, this approach detected a significant non-random distribution of pathogenic and benign variants in 387 (46%) and 172 (20%) genes, respectively, revealing that variant clustering is widespread across the human exome. This clustering likely occurs as a consequence of mechanisms shaping pathogenicity at the protein level, as illustrated by the overlap of some clusters with known functional domains. We then took advantage of these findings to develop a pathogenicity predictor, MutScore, that integrates qualitative features of DNA substitutions with the new additional information derived from this positional clustering. Using a random forest approach, MutScore was able to identify pathogenic missense mutations with very high accuracy, outperforming existing predictive tools, especially for variants associated with autosomal-dominant disease and cancer. Thus, the within-gene clustering of pathogenic and benign DNA changes is an important and previously underappreciated feature of the human exome, which can be harnessed to improve the prediction of pathogenicity and disambiguation of DNA variants of uncertain significance.

Genomic and transcriptomic landscape of conjunctival melanoma.

Conjunctival melanoma (CJM) is a rare but potentially lethal and highly-recurrent cancer of the eye. Similar to cutaneous melanoma (CM), it originates from melanocytes. Unlike CM, however, CJM is relatively poorly characterized from a genomic point of view. To fill this knowledge gap and gain insight into the genomic nature of CJM, we performed whole-exome (WES) or whole-genome sequencing (WGS) of tumor-normal tissue pairs in 14 affected individuals, as well as RNA sequencing in a subset of 11 tumor tissues. Our results show that, similarly to CM, CJM is also characterized by a very high mutation load, composed of approximately 500 somatic mutations in exonic regions. This, as well as the presence of a UV light-induced mutational signature, are clear signs of the role of sunlight in CJM tumorigenesis. In addition, the genomic classification of CM proposed by TCGA seems to be well-applicable to CJM, with the presence of four typical subclasses defined on the basis of the most frequently mutated genes: BRAF, NF1, RAS, and triple wild-type. In line with these results, transcriptomic analyses revealed similarities with CM as well, namely the presence of a transcriptomic subtype enriched for immune genes and a subtype enriched for genes associated with keratins and epithelial functions. Finally, in seven tumors we detected somatic mutations in ACSS3, a possible new candidate oncogene. Transfected conjunctival melanoma cells overexpressing mutant ACSS3 showed higher proliferative activity, supporting the direct involvement of this gene in the tumorigenesis of CJM. Altogether, our results provide the first unbiased and complete genomic and transcriptomic classification of CJM.

Peripheral neuropathy and cognitive impairment associated with a novel monoallelic HARS variant.

BACKGROUND: A 49-year-old male presented with late-onset demyelinating peripheral neuropathy, cerebellar atrophy, and cognitive deficit. Nerve biopsy revealed intra-axonal inclusions suggestive of polyglucosan bodies, raising the suspicion of adult polyglucosan bodies disease (OMIM 263570). METHODS AND RESULTS: While known genes associated with polyglucosan bodies storage were negative, whole-exome sequencing identified an unreported monoallelic variant, c.397G>T (p.Val133Phe), in the histidyl-tRNA synthetase (HARS) gene. While we did not identify mutations in genes known to be associated with polygucosan body disease, whole-exome sequencing revealed an unreported monoallelic variant, c.397G>T in the histidyl-tRNA synthetase (HARS) gene, encoding a substitution (Val133Phe) in the catalytic domain. Expression of this variant in patient cells resulted in reduced aminoacylation activity in extracts obtained from dermal fibroblasts, without compromising overall protein synthesis. INTERPRETATION: Genetic variants in the genes coding for the different aminoacyl-tRNA synthases are associated with various clinical conditions. To date, a number of HARS variant have been associated with peripheral neuropathy, but not cognitive deficits. Further studies are needed to explore why HARS mutations confer a neuronal-specific phenotype.

Lamin B receptor-related disorder is associated with a spectrum of skeletal dysplasia phenotypes.

LBR (Lamin B Receptor) encodes a bifunctional protein important for cholesterol biosynthesis and heterochromatin organization on the inner nuclear membrane. Pathogenic variants in LBR are associated with marked phenotypic variability, ranging from the benign Pelger-Huët anomaly to lethal Greenberg Dysplasia. We performed trio exome sequencing (ES) on two patients with atypical variants of skeletal dysplasia and their unaffected parents. Patient 1 exhibited frontal bossing, mid-face hypoplasia, short stature with rhizomelic limb shortening, and relative macrocephaly at birth. Although remained short, Patient 1 later showed spontaneous improvement in her skeletal findings. Exome sequencing revealed two novel variants in LBR, c.1504C > G (p.Arg502Gly) in exon 12 and c.1748G > T (p.Arg583Leu) in exon 14, which were inherited from her unaffected father and mother, respectively. Sterol analysis revealed an increased level of cholesta­8,14­dien­3ß­ol to 2.9% of total sterols, consistent with a functional deficiency of 3ß­hydroxysterol Δ14­reductase. Patient 2 presented at birth with short stature and marked rhizomelic limb shortening but later exhibited decreasing severity of shortening of the long bones and improvement in the radiographic skeletal abnormalities although he continued to be significantly short at age 10 years. Exome sequencing revealed that Patient 2 is homozygous for a pathogenic variant c.1534C > T (p.Arg512Trp) in exon 12 of LBR, which was inherited from his unaffected consanguineous parents. This report provides further evidence for a phenotypic spectrum of LBR-associated disorders and expands the genotypic spectrum by describing 3 novel disease-causing variants that have not been previously associated with a disease. Moreover, our data on Patient 1 demonstrate that variants throughout the gene appear to influence both the sterol reductase and nuclear functions of LBR.

EXTL3 mutations cause skeletal dysplasia, immune deficiency, and developmental delay.

We studied three patients with severe skeletal dysplasia, T cell immunodeficiency, and developmental delay. Whole-exome sequencing revealed homozygous missense mutations affecting exostosin-like 3 (EXTL3), a glycosyltransferase involved in heparan sulfate (HS) biosynthesis. Patient-derived fibroblasts showed abnormal HS composition and altered fibroblast growth factor 2 signaling, which was rescued by overexpression of wild-type EXTL3 cDNA. Interleukin-2-mediated STAT5 phosphorylation in patients' lymphocytes was markedly reduced. Interbreeding of the extl3-mutant zebrafish (box) with Tg(rag2:green fluorescent protein) transgenic zebrafish revealed defective thymopoiesis, which was rescued by injection of wild-type human EXTL3 RNA. Targeted differentiation of patient-derived induced pluripotent stem cells showed a reduced expansion of lymphohematopoietic progenitor cells and defects of thymic epithelial progenitor cell differentiation. These data identify EXTL3 mutations as a novel cause of severe immune deficiency with skeletal dysplasia and developmental delay and underline a crucial role of HS in thymopoiesis and skeletal and brain development.

Comprehensive Genetic Landscape of Uveal Melanoma by Whole-Genome Sequencing.

Uveal melanoma (UM) is a rare intraocular tumor that, similar to cutaneous melanoma, originates from melanocytes. To gain insights into its genetics, we performed whole-genome sequencing at very deep coverage of tumor-control pairs in 33 samples (24 primary and 9 metastases). Genome-wide, the number of coding mutations was rather low (only 17 variants per tumor on average; range 7-28), thus radically different from cutaneous melanoma, where hundreds of exonic DNA insults are usually detected. Furthermore, no UV light-induced mutational signature was identified. Recurrent coding mutations were found in the known UM drivers GNAQ, GNA11, BAP1, EIF1AX, and SF3B1. Other genes, i.e., TP53BP1, CSMD1, TTC28, DLK2, and KTN1, were also found to harbor somatic mutations in more than one individual, possibly indicating a previously undescribed association with UM pathogenesis. De novo assembly of unmatched reads from non-coding DNA revealed peculiar copy-number variations defining specific UM subtypes, which in turn could be associated with metastatic transformation. Mutational-driven comparison with other tumor types showed that UM is very similar to pediatric tumors, characterized by very few somatic insults and, possibly, important epigenetic changes. Through the analysis of whole-genome sequencing data, our findings shed new light on the molecular genetics of uveal melanoma, delineating it as an atypical tumor of the adult for which somatic events other than mutations in exonic DNA shape its genetic landscape and define its metastatic potential.

NANS-mediated synthesis of sialic acid is required for brain and skeletal development.

We identified biallelic mutations in NANS, the gene encoding the synthase for N-acetylneuraminic acid (NeuNAc; sialic acid), in nine individuals with infantile-onset severe developmental delay and skeletal dysplasia. Patient body fluids showed an elevation in N-acetyl-D-mannosamine levels, and patient-derived fibroblasts had reduced NANS activity and were unable to incorporate sialic acid precursors into sialylated glycoproteins. Knockdown of nansa in zebrafish embryos resulted in abnormal skeletal development, and exogenously added sialic acid partially rescued the skeletal phenotype. Thus, NANS-mediated synthesis of sialic acid is required for early brain development and skeletal growth. Normal sialylation of plasma proteins was observed in spite of NANS deficiency. Exploration of endogenous synthesis, nutritional absorption, and rescue pathways for sialic acid in different tissues and developmental phases is warranted to design therapeutic strategies to counteract NANS deficiency and to shed light on sialic acid metabolism and its implications for human nutrition.

Mutations in the heat-shock protein A9 (HSPA9) gene cause the EVEN-PLUS syndrome of congenital malformations and skeletal dysplasia.

We and others have reported mutations in LONP1, a gene coding for a mitochondrial chaperone and protease, as the cause of the human CODAS (cerebral, ocular, dental, auricular and skeletal) syndrome (MIM 600373). Here, we delineate a similar but distinct condition that shares the epiphyseal, vertebral and ocular changes of CODAS but also included severe microtia, nasal hypoplasia, and other malformations, and for which we propose the name of EVEN-PLUS syndrome for epiphyseal, vertebral, ear, nose, plus associated findings. In three individuals from two families, no mutation in LONP1 was found; instead, we found biallelic mutations in HSPA9, the gene that codes for mHSP70/mortalin, another highly conserved mitochondrial chaperone protein essential in mitochondrial protein import, folding, and degradation. The functional relationship between LONP1 and HSPA9 in mitochondrial protein chaperoning and the overlapping phenotypes of CODAS and EVEN-PLUS delineate a family of "mitochondrial chaperonopathies" and point to an unexplored role of mitochondrial chaperones in human embryonic morphogenesis.

NBAS mutations cause a multisystem disorder involving bone, connective tissue, liver, immune system, and retina.

We report two unrelated patients with a multisystem disease involving liver, eye, immune system, connective tissue, and bone, caused by biallelic mutations in the neuroblastoma amplified sequence (NBAS) gene. Both presented as infants with recurrent episodes triggered by fever with vomiting, dehydration, and elevated transaminases. They had frequent infections, hypogammaglobulinemia, reduced natural killer cells, and the Pelger-Huët anomaly of their granulocytes. Their facial features were similar with a pointed chin and proptosis; loose skin and reduced subcutaneous fat gave them a progeroid appearance. Skeletal features included short stature, slender bones, epiphyseal dysplasia with multiple phalangeal pseudo-epiphyses, and small C1-C2 vertebrae causing cervical instability and myelopathy. Retinal dystrophy and optic atrophy were present in one patient. NBAS is a component of the synthaxin-18 complex and is involved in nonsense-mediated mRNA decay control. Putative loss-of-function mutations in NBAS are already known to cause disease in humans. A specific founder mutation has been associated with short stature, optic nerve atrophy and Pelger-Huët anomaly of granulocytes (SOPH) in the Siberian Yakut population. A more recent report associates NBAS mutations with recurrent acute liver failure in infancy in a group of patients of European descent. Our observations indicate that the phenotypic spectrum of NBAS deficiency is wider than previously known and includes skeletal, hepatic, metabolic, and immunologic aspects. Early recognition of the skeletal phenotype is important for preventive management of cervical instability.

Whole genome sequencing as a means to assess pathogenic mutations in medical genetics and cancer.

The past decade has seen the emergence of next-generation sequencing (NGS) technologies, which have revolutionized the field of human molecular genetics. With NGS, significant portions of the human genome can now be assessed by direct sequence analysis, highlighting normal and pathological variants of our DNA. Recent advances have also allowed the sequencing of complete genomes, by a method referred to as whole genome sequencing (WGS). In this work, we review the use of WGS in medical genetics, with specific emphasis on the benefits and the disadvantages of this technique for detecting genomic alterations leading to Mendelian human diseases and to cancer.

Functional chromatin features are associated with structural mutations in cancer.

BACKGROUND: Structural mutations (SMs) play a major role in cancer development. In some cancers, such as breast and ovarian, DNA double-strand breaks (DSBs) occur more frequently in transcribed regions, while in other cancer types such as prostate, there is a consistent depletion of breakpoints in transcribed regions. Despite such regularity, little is understood about the mechanisms driving these effects. A few works have suggested that protein binding may be relevant, e.g. in studies of androgen receptor binding and active chromatin in specific cell types. We hypothesized that this behavior might be general, i.e. that correlation between protein-DNA binding (and open chromatin) and breakpoint locations is common across divergent cancers. RESULTS: We investigated this hypothesis by comprehensively analyzing the relationship among 457 ENCODE protein binding ChIP-seq experiments, 125 DnaseI and 24 FAIRE experiments, and 14,600 SMs from 8 diverse cancer datasets covering 147 samples. In most cancers, including breast and ovarian, we found enrichment of protein binding and open chromatin in the vicinity of SM breakpoints at distances up to 200 kb. Furthermore, for all cancer types we observed an enhanced enrichment in regions distant from genes when compared to regions proximal to genes, suggesting that the SM-induction mechanism is independent from the bias of DSBs to occur near transcribed regions. We also observed a stronger effect for sites with more than one protein bound. CONCLUSIONS: Protein binding and open chromatin state are associated with nearby SM breakpoints in many cancer datasets. These observations suggest a consistent mechanism underlying SM locations across different cancers.