Non-coding deletions identify Maenli lncRNA as a limb-specific En1 regulator.

Long non-coding RNAs (lncRNAs) can be important components in gene-regulatory networks1, but the exact nature and extent of their involvement in human Mendelian disease is largely unknown. Here we show that genetic ablation of a lncRNA locus on human chromosome 2 causes a severe congenital limb malformation. We identified homozygous 27-63-kilobase deletions located 300 kilobases upstream of the engrailed-1 gene (EN1) in patients with a complex limb malformation featuring mesomelic shortening, syndactyly and ventral nails (dorsal dimelia). Re-engineering of the human deletions in mice resulted in a complete loss of En1 expression in the limb and a double dorsal-limb phenotype that recapitulates the human disease phenotype. Genome-wide transcriptome analysis in the developing mouse limb revealed a four-exon-long non-coding transcript within the deleted region, which we named Maenli. Functional dissection of the Maenli locus showed that its transcriptional activity is required for limb-specific En1 activation in cis, thereby fine-tuning the gene-regulatory networks controlling dorso-ventral polarity in the developing limb bud. Its loss results in the En1-related dorsal ventral limb phenotype, a subset of the full En1-associated phenotype. Our findings demonstrate that mutations involving lncRNA loci can result in human Mendelian disease.

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.

DOMINO: Using Machine Learning to Predict Genes Associated with Dominant Disorders.

In contrast to recessive conditions with biallelic inheritance, identification of dominant (monoallelic) mutations for Mendelian disorders is more difficult, because of the abundance of benign heterozygous variants that act as massive background noise (typically, in a 400:1 excess ratio). To reduce this overflow of false positives in next-generation sequencing (NGS) screens, we developed DOMINO, a tool assessing the likelihood for a gene to harbor dominant changes. Unlike commonly-used predictors of pathogenicity, DOMINO takes into consideration features that are the properties of genes, rather than of variants. It uses a machine-learning approach to extract discriminant information from a broad array of features (N = 432), including: genomic data, intra-, and interspecies conservation, gene expression, protein-protein interactions, protein structure, etc. DOMINO's iterative architecture includes a training process on 985 genes with well-established inheritance patterns for Mendelian conditions, and repeated cross-validation that optimizes its discriminant power. When validated on 99 newly-discovered genes with pathogenic mutations, the algorithm displays an excellent final performance, with an area under the curve (AUC) of 0.92. Furthermore, unsupervised analysis by DOMINO of real sets of NGS data from individuals with intellectual disability or epilepsy correctly recognizes known genes and predicts 9 new candidates, with very high confidence. In summary, DOMINO is a robust and reliable tool that can infer dominance of candidate genes with high sensitivity and specificity, making it a useful complement to any NGS pipeline dealing with the analysis of the morbid human genome.

Childhood neurodegeneration associated with a specific UBTF variant: a new case report and review of the literature.

BACKGROUND: A new monogenic neurodegenerative disease affecting ribosomal metabolism has recently been identified in association with a monoallelic UBTF putative gain of function variant (NM_001076683.1:c.628G>A, hg19). Phenotype is consistent among these probands with progressive motor, cognitive, and behavioural regression in early to middle childhood. CASE PRESENTATION: We report on a child with this monoallelic UBTF variant who presented with progressive disease including regression, episodes of subacute deterioration during febrile illnesses and a remarkable EEG pattern with a transient pattern of semi-periodic slow waves. CONCLUSIONS: This case further supports the phenotype-genotype correlation of neurodegeneration associated with UBTF c.628G>A. Moreover, it brings new insights into the clinical features and EEG that could possibly serve as diagnostic markers of this otherwise nonspecific phenotype.

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.

Mutations in CEP78 Cause Cone-Rod Dystrophy and Hearing Loss Associated with Primary-Cilia Defects.

Cone-rod degeneration (CRD) belongs to the disease spectrum of retinal degenerations, a group of hereditary disorders characterized by an extreme clinical and genetic heterogeneity. It mainly differentiates from other retinal dystrophies, and in particular from the more frequent disease retinitis pigmentosa, because cone photoreceptors degenerate at a higher rate than rod photoreceptors, causing severe deficiency of central vision. After exome analysis of a cohort of individuals with CRD, we identified biallelic mutations in the orphan gene CEP78 in three subjects from two families: one from Greece and another from Sweden. The Greek subject, from the island of Crete, was homozygous for the c.499+1G>T (IVS3+1G>T) mutation in intron 3. The Swedish subjects, two siblings, were compound heterozygotes for the nearby mutation c.499+5G>A (IVS3+5G>A) and for the frameshift-causing variant c.633delC (p.Trp212Glyfs(∗)18). In addition to CRD, these three individuals had hearing loss or hearing deficit. Immunostaining highlighted the presence of CEP78 in the inner segments of retinal photoreceptors, predominantly of cones, and at the base of the primary cilium of fibroblasts. Interaction studies also showed that CEP78 binds to FAM161A, another ciliary protein associated with retinal degeneration. Finally, analysis of skin fibroblasts derived from affected individuals revealed abnormal ciliary morphology, as compared to that of control cells. Altogether, our data strongly suggest that mutations in CEP78 cause a previously undescribed clinical entity of a ciliary nature characterized by blindness and deafness but clearly distinct from Usher syndrome, a condition for which visual impairment is due to retinitis pigmentosa.

Mutations in the polyglutamylase gene TTLL5, expressed in photoreceptor cells and spermatozoa, are associated with cone-rod degeneration and reduced male fertility.

Hereditary retinal degenerations encompass a group of genetic diseases characterized by extreme clinical variability. Following next-generation sequencing and autozygome-based screening of patients presenting with a peculiar, recessive form of cone-dominated retinopathy, we identified five homozygous variants [p.(Asp594fs), p.(Gln117*), p.(Met712fs), p.(Ile756Phe), and p.(Glu543Lys)] in the polyglutamylase-encoding gene TTLL5, in eight patients from six families. The two male patients carrying truncating TTLL5 variants also displayed a substantial reduction in sperm motility and infertility, whereas those carrying missense changes were fertile. Defects in this polyglutamylase in humans have recently been associated with cone photoreceptor dystrophy, while mouse models carrying truncating mutations in the same gene also display reduced fertility in male animals. We examined the expression levels of TTLL5 in various human tissues and determined that this gene has multiple viable isoforms, being highly expressed in testis and retina. In addition, antibodies against TTLL5 stained the basal body of photoreceptor cells in rat and the centrosome of the spermatozoon flagellum in humans, suggesting a common mechanism of action in these two cell types. Taken together, our data indicate that mutations in TTLL5 delineate a novel, allele-specific syndrome causing defects in two as yet pathogenically unrelated functions, reproduction and vision.

[Next generation sequencing : a diagnostic tool for inherited immune defects]. / Séquençage à haut débit : outil de diagnostic des déficits immunitaires héréditaires.

Establishing the definitive diagnosis in the case of inherited immune defects (IID) is often challenging because the clinical features can be heterogeneous, atypical and overlapping different disease entities. The next generation sequencing technology (NGS) allows identifying genetic variants that are responsible for the observed clinical presentations. The use of NGS applied to the genes mutated in IIDs or known to be involved in the development, differentiation and regulation of the immune system allows to target hundreds of relevant genes in well characterized patients suspected of carrying inherited immune defects. This approach answers both diagnostic and research needs, facilitates the understanding of the mechanisms that underlie IIDs, and ultimately leads to the discovery of new therapeutic targets.

Etablir un diagnostic précis est souvent difficile dans le cas des déficits immunitaires héréditaires (DIH) car les manifestations cliniques sont hétérogènes et parfois atypiques ou communes à des maladies distinctes. Le séquençage à haut débit, aussi appelé next generation sequencing (NGS), permet d'identifier des variants génétiques responsables de caractéristiques cliniques. Le recours au NGS pour découvrir les gènes mutés dans les DIH ou connus pour être impliqués dans le développement, la différenciation et la régulation du système immunitaire, permet de cibler des centaines de gènes d'intérêt sur des patients suspectés de DIH pour lesquels le tableau clinique est documenté. Cette approche répond à la fois à des besoins de diagnostic et de recherche pour comprendre les mécanismes qui régissent les DIH et découvrir de nouveaux outils ou cibles thérapeutiques.

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.

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.

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.

An atypical form of 60S ribosomal subunit in Diamond-Blackfan anemia linked to RPL17 variants.

Diamond-Blackfan anemia syndrome (DBA) is a ribosomopathy associated with loss-of-function variants in more than 20 ribosomal protein (RP) genes. Here, we report the genetic, functional, and biochemical dissection of 2 multigenerational pedigrees with variants in RPL17, a large ribosomal subunit protein-encoding gene. Affected individuals had clinical features and erythroid proliferation defects consistent with DBA. Further, RPL17/uL22 depletion resulted in anemia and micrognathia in zebrafish larvae, and in vivo complementation studies indicated that RPL17 variants were pathogenic. Lymphoblastoid cell lines (LCLs) derived from patients displayed a ribosomal RNA maturation defect reflecting haploinsufficiency of RPL17. The proteins encoded by RPL17 variants were not incorporated into ribosomes, but 10%-20% of 60S ribosomal subunits contained a short form of 5.8S rRNA (5.8SC), a species that is marginal in normal cells. These atypical 60S subunits were actively engaged in translation. Ribosome profiling showed changes of the translational profile, but those are similar to LCLs bearing RPS19 variants. These results link an additional RP gene to DBA. They show that ribosomes can be modified substantially by RPL17 haploinsufficiency but support the paradigm that translation alterations in DBA are primarily related to insufficient ribosome production rather than to changes in ribosome structure or composition.

CNV Detection from Exome Sequencing Data in Routine Diagnostics of Rare Genetic Disorders: Opportunities and Limitations.

To assess the potential of detecting copy number variations (CNVs) directly from exome sequencing (ES) data in diagnostic settings, we developed a CNV-detection pipeline based on ExomeDepth software and applied it to ES data of 450 individuals. Initially, only CNVs affecting genes in the requested diagnostic gene panels were scored and tested against arrayCGH results. Pathogenic CNVs were detected in 18 individuals. Most detected CNVs were larger than 400 kb (11/18), but three individuals had small CNVs impacting one or a few exons only and were thus not detectable by arrayCGH. Conversely, two pathogenic CNVs were initially missed, as they impacted genes not included in the original gene panel analysed, and a third one was missed as it was in a poorly covered region. The overall combined diagnostic rate (SNVs + CNVs) in our cohort was 36%, with wide differences between clinical domains. We conclude that (1) the ES-based CNV pipeline detects efficiently large and small pathogenic CNVs, (2) the detection of CNV relies on uniformity of sequencing and good coverage, and (3) in patients who remain unsolved by the gene panel analysis, CNV analysis should be extended to all captured genes, as diagnostically relevant CNVs may occur everywhere in the genome.

De novo variants in CACNA1E found in patients with intellectual disability, developmental regression and social cognition deficit but no seizures.

BACKGROUND: De novo variants in the voltage-gated calcium channel subunit α1 E gene (CACNA1E) have been described as causative of epileptic encephalopathy with contractures, macrocephaly and dyskinesias. METHODS: Following the observation of an index patient with developmental delay and autism spectrum disorder (ASD) without seizures who had a de novo deleterious CACNA1E variant, we screened GeneMatcher for other individuals with CACNA1E variants and neurodevelopmental phenotypes without epilepsy. The spectrum of pathogenic CACNA1E variants was compared to the mutational landscape of variants in the gnomAD control population database. RESULTS: We identified seven unrelated individuals with intellectual disability, developmental regression and ASD-like behavioral profile, and notably without epilepsy, who had de novo heterozygous putatively pathogenic variants in CACNA1E. Age of onset of clinical manifestation, presence or absence of regression and degree of severity were variable, and no clear-cut genotype-phenotype association could be recognized. The analysis of disease-associated variants and their comparison to benign variants from the control population allowed for the identification of regions in the CACNA1E protein that seem to be intolerant to substitutions and thus more likely to harbor pathogenic variants. As in a few reported cases with CACNA1E variants and epilepsy, one patient showed a positive clinical behavioral response to topiramate, a specific calcium channel modulator. LIMITATIONS: The significance of our study is limited by the absence of functional experiments of the effect of identified variants, the small sample size and the lack of systematic ASD assessment in all participants. Moreover, topiramate was given to one patient only and for a short period of time. CONCLUSIONS: Our results indicate that CACNA1E variants may result in neurodevelopmental disorders without epilepsy and expand the mutational and phenotypic spectrum of this gene. CACNA1E deserves to be included in gene panels for non-specific developmental disorders, including ASD, and not limited to patients with seizures, to improve diagnostic recognition and explore the possible efficacy of topiramate.

AutoMap is a high performance homozygosity mapping tool using next-generation sequencing data.

Homozygosity mapping is a powerful method for identifying mutations in patients with recessive conditions, especially in consanguineous families or isolated populations. Historically, it has been used in conjunction with genotypes from highly polymorphic markers, such as DNA microsatellites or common SNPs. Traditional software performs rather poorly with data from Whole Exome Sequencing (WES) and Whole Genome Sequencing (WGS), which are now extensively used in medical genetics. We develop AutoMap, a tool that is both web-based or downloadable, to allow performing homozygosity mapping directly on VCF (Variant Call Format) calls from WES or WGS projects. Following a training step on WES data from 26 consanguineous families and a validation procedure on a matched cohort, our method shows higher overall performances when compared with eight existing tools. Most importantly, when tested on real cases with negative molecular diagnosis from an internal set, AutoMap detects three gene-disease and multiple variant-disease associations that were previously unrecognized, projecting clear benefits for both molecular diagnosis and research activities in medical genetics.

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.