Bi-allelic variants in COQ8B, a gene involved in the biosynthesis of coenzyme Q10, lead to non-syndromic retinitis pigmentosa.

Retinitis pigmentosa (RP) is a Mendelian disease characterized by gradual loss of vision, due to the progressive degeneration of retinal cells. Genetically, it is highly heterogeneous, with pathogenic variants identified in more than 100 genes so far. Following a large-scale sequencing screening, we identified five individuals (four families) with recessive and non-syndromic RP, carrying as well bi-allelic DNA changes in COQ8B, a gene involved in the biosynthesis of coenzyme Q10. Specifically, we detected compound heterozygous assortments of five disease-causing variants (c.187C>T [p.Arg63Trp], c.566G>A [p.Trp189Ter], c.1156G>A [p.Asp386Asn], c.1324G>A [p.Val442Met], and c.1560G>A [p.Trp520Ter]), all segregating with disease according to a recessive pattern of inheritance. Cell-based analysis of recombinant proteins deriving from these genotypes, performed by target engagement assays, showed in all cases a significant decrease in ligand-protein interaction compared to the wild type. Our results indicate that variants in COQ8B lead to recessive non-syndromic RP, possibly by impairing the biosynthesis of coenzyme Q10, a key component of oxidative phosphorylation in the mitochondria.

Mutations in SAMD7 cause autosomal-recessive macular dystrophy with or without cone dysfunction.

Sterile alpha motif domain containing 7 (SAMD7) is a component of the Polycomb repressive complex 1, which inhibits transcription of many genes, including those activated by the transcription factor Cone-Rod Homeobox (CRX). Here we report bi-allelic mutations in SAMD7 as a cause of autosomal-recessive macular dystrophy with or without cone dysfunction. Four of these mutations affect splicing, while another mutation is a missense variant that alters the repressive effect of SAMD7 on CRX-dependent promoter activity, as shown by in vitro assays. Immunostaining of human retinal sections revealed that SAMD7 is localized in the nuclei of both rods and cones, as well as in those of cells belonging to the inner nuclear layer. These results place SAMD7 as a gene crucial for human retinal function and demonstrate a significant difference in the role of SAMD7 between the human and the mouse retina.

Detection of elusive DNA copy-number variations in hereditary disease and cancer through the use of noncoding and off-target sequencing reads.

Copy-number variants (CNVs) play a substantial role in the molecular pathogenesis of hereditary disease and cancer, as well as in normal human interindividual variation. However, they are still rather difficult to identify in mainstream sequencing projects, especially involving exome sequencing, because they often occur in DNA regions that are not targeted for analysis. To overcome this problem, we developed OFF-PEAK, a user-friendly CNV detection tool that builds on a denoising approach and the use of "off-target" DNA reads, which are usually discarded by sequencing pipelines. We benchmarked OFF-PEAK on data from targeted sequencing of 96 cancer samples, as well as 130 exomes of individuals with inherited retinal disease from three different populations. For both sets of data, OFF-PEAK demonstrated excellent performance (>95% sensitivity and >80% specificity vs. experimental validation) in detecting CNVs from in silico data alone, indicating its immediate applicability to molecular diagnosis and genetic research.

Substitution of a single non-coding nucleotide upstream of TMEM216 causes non-syndromic retinitis pigmentosa and is associated with reduced TMEM216 expression.

Genome analysis of individuals affected by retinitis pigmentosa (RP) identified two rare nucleotide substitutions at the same genomic location on chromosome 11 (g.61392563 [GRCh38]), 69 base pairs upstream of the start codon of the ciliopathy gene TMEM216 (c.-69G>A, c.-69G>T [GenBank: NM_001173991.3]), in individuals of South Asian and African ancestry, respectively. Genotypes included 71 homozygotes and 3 mixed heterozygotes in trans with a predicted loss-of-function allele. Haplotype analysis showed single-nucleotide variants (SNVs) common across families, suggesting ancestral alleles within the two distinct ethnic populations. Clinical phenotype analysis of 62 available individuals from 49 families indicated a similar clinical presentation with night blindness in the first decade and progressive peripheral field loss thereafter. No evident systemic ciliopathy features were noted. Functional characterization of these variants by luciferase reporter gene assay showed reduced promotor activity. Nanopore sequencing confirmed the lower transcription of the TMEM216 c.-69G>T allele in blood-derived RNA from a heterozygous carrier, and reduced expression was further recapitulated by qPCR, using both leukocytes-derived RNA of c.-69G>T homozygotes and total RNA from genome-edited hTERT-RPE1 cells carrying homozygous TMEM216 c.-69G>A. In conclusion, these variants explain a significant proportion of unsolved cases, specifically in individuals of African ancestry, suggesting that reduced TMEM216 expression might lead to abnormal ciliogenesis and photoreceptor degeneration.

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.

Multidisciplinary team directed analysis of whole genome sequencing reveals pathogenic non-coding variants in molecularly undiagnosed inherited retinal dystrophies.

The purpose of this paper is to identify likely pathogenic non-coding variants in inherited retinal dystrophy (IRD) genes, using genome sequencing (GS). Patients with IRD were recruited to the study and underwent comprehensive ophthalmological evaluation and GS. The results of GS were investigated through virtual gene panel analysis, and plausible pathogenic variants and clinical phenotype evaluated by the multidisciplinary team (MDT) discussion. For unsolved patients in whom a specific gene was suspected to harbor a missed pathogenic variant, targeted re-analysis of non-coding regions was performed on GS data. Candidate variants were functionally tested by messenger RNA analysis, minigene or luciferase reporter assays. Previously unreported, likely pathogenic, non-coding variants in 7 genes (PRPF31, NDP, IFT140, CRB1, USH2A, BBS10 and GUCY2D), were identified in 11 patients. These were shown to lead to mis-splicing (PRPF31, IFT140, CRB1 and USH2A) or altered transcription levels (BBS10 and GUCY2D). MDT-led, phenotype-driven, non-coding variant re-analysis of GS is effective in identifying the missing causative alleles.

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.

Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative.

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

Mutations in the ribosome biogenesis factor gene LTV1 are linked to LIPHAK syndrome, a novel poikiloderma-like disorder.

In the framework of the UK 100 000 Genomes Project, we investigated the genetic origin of a previously undescribed recessive dermatological condition, which we named LIPHAK (LTV1-associated Inflammatory Poikiloderma with Hair abnormalities and Acral Keratoses), in four affected individuals from two UK families of Pakistani and Indian origins, respectively. Our analysis showed that only one gene, LTV1, carried rare biallelic variants that were shared in all affected individuals, and specifically they bore the NM_032860.5:c.503A > G, p.(Asn168Ser) change, found homozygously in all of them. In addition, high-resolution homozygosity mapping revealed the presence of a small 652-kb stretch on chromosome 6, encompassing LTV1, that was haploidentical and common to all affected individuals. The c.503A > G variant was predicted by in silico tools to affect the correct splicing of LTV1's exon 5. Minigene-driven splicing assays in HEK293T cells and in a skin sample from one of the patients confirmed that this variant was indeed responsible for the creation of a new donor splice site, resulting in aberrant splicing and in a premature termination codon in exon 6 of this gene. LTV1 encodes one of the ribosome biogenesis factors that promote the assembly of the small (40S) ribosomal subunit. In yeast, defects in LTV1 alter the export of nascent ribosomal subunits to the cytoplasm; however, the role of this gene in human pathology is unknown to date. Our data suggest that LIPHAK could be a previously unrecognized ribosomopathy.

Loss-of-function variants in UBAP1L cause autosomal recessive retinal degeneration.

PURPOSE: Inherited retinal diseases (IRDs) are a group of monogenic conditions that can lead to progressive blindness. Their missing heritability is still considerable, due in part to the presence of disease genes that await molecular identification. The purpose of this work was to identify novel genetic associations with IRDs. METHODS: Patients underwent a comprehensive ophthalmological evaluation using standard-of-care tests, such as detailed retinal imaging (macular optical coherence tomography and short-wavelength fundus autofluorescence) and electrophysiological testing. Exome and genome sequencing, as well as computer-assisted data analysis were used for genotyping and detection of DNA variants. A minigene-driven splicing assay was performed to validate the deleterious effects of 1 of such variants. RESULTS: We identified 8 unrelated families from Hungary, the United States, Israel, and The Netherlands with members presenting with a form of autosomal recessive and nonsyndromic retinal degeneration, predominantly described as rod-cone dystrophy but also including cases of cone/cone-rod dystrophy. Age of disease onset was very variable, with some patients experiencing first symptoms during their fourth decade of life or later. Myopia greater than 5 diopters was present in 5 of 7 cases with available refractive data, and retinal detachment was reported in 2 cases. All ascertained patients carried biallelic loss-of-function variants in UBAP1L (HGNC: 40028), a gene with unknown function and with homologies to UBAP1, encoding a protein involved in ubiquitin metabolism. One of these pathogenic variants, the intronic NM_001163692.2:c.910-7G>A substitution, was identified in 5 unrelated families. Minigene-driven splicing assays in HEK293T cells confirmed that this DNA change is responsible for the creation of a new acceptor splice site, resulting in aberrant splicing. CONCLUSION: We identified UBAP1L as a novel IRD gene. Although its function is currently unknown, UBAP1L is almost exclusively expressed in photoreceptors and the retinal pigment epithelium, hence possibly explaining the link between pathogenic variants in this gene and an ocular phenotype.

The p.C759F Variant in USH2A Is a Pathogenic Mutation: Systematic Literature Review and Meta-Analysis of 667 Genotypes.

BACKGROUND: Although the p.C759F (c.2276G>T, p.Cys759Phe) variant in the USH2A gene has been identified in association with retinal degeneration by several authors, its pathogenicity has been questioned once by the publication of two unaffected homozygotes from a single family. OBJECTIVES: The objective of the study was to ascertain the role of p.C759F in hereditary retinal disease. METHODS: We examined 87 research articles reporting on patients carrying this variant and then used this information as primary data for a series of meta-analytical tests. RESULTS: Independent statistical analyses showed that p.C759F (i) is highly enriched in patients with respect to healthy individuals, (ii) represents a clear-cut recessive allele causing disease when it is in trans with other mutations, (iii) is pathogenic in homozygotes. CONCLUSIONS: Our results confirm that p.C759F is a bona fide mutation, leading to retinal blindness according to a recessive pattern of inheritance.

Biallelic variants in ZNF526 cause a severe neurodevelopmental disorder with microcephaly, bilateral cataract, epilepsy and simplified gyration.

BACKGROUND: Next-generation sequencing, combined with international pooling of cases, has impressively enhanced the discovery of genes responsible for Mendelian neurodevelopmental disorders, particularly in individuals affected by clinically undiagnosed diseases. To date, biallelic missense variants in ZNF526 gene, encoding a Krüppel-type zinc-finger protein, have been reported in three families with non-syndromic intellectual disability. METHODS: Here, we describe five individuals from four unrelated families with an undiagnosed neurodevelopmental disorder in which we performed exome sequencing, on a combination of trio-based (4 subjects) or single probands (1 subject). RESULTS: We identified five patients from four unrelated families with homozygous ZNF526 variants by whole exome sequencing. Four had variants resulting in truncation of ZNF526; they were affected by severe prenatal and postnatal microcephaly (ranging from -4 SD to -8 SD), profound psychomotor delay, hypertonic-dystonic movements, epilepsy and simplified gyral pattern on MRI. All of them also displayed bilateral progressive cataracts. A fifth patient had a homozygous missense variant and a slightly less severe disorder, with postnatal microcephaly (-2 SD), progressive bilateral cataracts, severe intellectual disability and unremarkable brain MRI.Mutant znf526 zebrafish larvae had notable malformations of the eye and central nervous system, resembling findings seen in the human holoprosencephaly spectrum. CONCLUSION: Our findings support the role of ZNF526 biallelic variants in a complex neurodevelopmental disorder, primarily affecting brain and eyes, resulting in severe microcephaly, simplified gyral pattern, epileptic encephalopathy and bilateral cataracts.

A Novel Intronic Deletion in PDE6B Causes Autosomal Recessive Retinitis Pigmentosa by Interfering with RNA Splicing.

INTRODUCTION: Retinitis pigmentosa (RP) is a rare degenerative retinal disease caused by mutations in approximately seventy genes. Currently, despite the availability of large-scale DNA sequencing technologies, ∼30-40% of patients still cannot be diagnosed at the molecular level. In this study, we investigated a novel intronic deletion of PDE6B, encoding the beta subunit of phosphodiesterase 6 in association with recessive RP. METHODS: Three unrelated consanguineous families were recruited from the northwestern part of Pakistan. Whole exome sequencing was performed for the proband of each family, and the data were analyzed according to an in-house computer pipeline. Relevant DNA variants in all available members of these families were assessed through Sanger sequencing. A minigene-based splicing assay was also performed. RESULTS: The clinical phenotype for all patients was compatible with rod cone degeneration, with the onset during childhood. Whole exome sequencing revealed a homozygous 18 bp intronic deletion (NM_000283.3:c.1921-20_1921-3del) in PDE6B, which co-segregated with disease in 10 affected individuals. In vitro splicing tests showed that this deletion causes aberrant RNA splicing of the gene, leading to the in-frame deletion of 6 codons and, likely, to disease. CONCLUSION: Our findings further expand the mutational spectrum of the PDE6B gene.

An Alu-mediated duplication in NMNAT1, involved in NAD biosynthesis, causes a novel syndrome, SHILCA, affecting multiple tissues and organs.

We investigated the genetic origin of the phenotype displayed by three children from two unrelated Italian families, presenting with a previously unrecognized autosomal recessive disorder that included a severe form of spondylo-epiphyseal dysplasia, sensorineural hearing loss, intellectual disability and Leber congenital amaurosis (SHILCA), as well as some brain anomalies that were visible at the MRI. Autozygome-based analysis showed that these children shared a 4.76 Mb region of homozygosity on chromosome 1, with an identical haplotype. Nonetheless, whole-exome sequencing failed to identify any shared rare coding variants, in this region or elsewhere. We then determined the transcriptome of patients' fibroblasts by RNA sequencing, followed by additional whole-genome sequencing experiments. Gene expression analysis revealed a 4-fold downregulation of the gene NMNAT1, residing indeed in the shared autozygous interval. Short- and long-read whole-genome sequencing highlighted a duplication involving 2 out of the 5 exons of NMNAT1 main isoform (NM_022787.3), leading to the production of aberrant mRNAs. Pathogenic variants in NMNAT1 have been previously shown to cause non-syndromic Leber congenital amaurosis (LCA). However, no patient with null biallelic mutations has ever been described, and murine Nmnat1 knockouts show embryonic lethality, indicating that complete absence of NMNAT1 activity is probably not compatible with life. The rearrangement found in our cases, presumably causing a strong but not complete reduction of enzymatic activity, may therefore result in an intermediate syndromic phenotype with respect to LCA and lethality.

Genetic Analysis of Consanguineous Pakistani Families with Congenital Stationary Night Blindness.

INTRODUCTION: Congenital stationary night blindness (CSNB) is a rare, largely nonprogressive, inherited retinal disorder that can be clinically classified on the basis of fundus and electroretinogram abnormalities. METHODS: We analyzed four large consanguineous families from the Southern Punjab region of Pakistan including multiple individuals affected with CSNB. Exome sequencing was performed in probands of all four families; Sanger sequencing was performed in additional members to test co-segregation of the variants identified. RESULTS: We identified two novel and likely pathogenic variants in two pedigrees, namely, NM_002905.4:c.668A>C (p.Gln223Pro) in RDH5 and NM_022567.2:c.908del (p.Gly303ValfsTer45) in NYX. In the two other families, the variants NM_002905.4:c.319G>C (p.Gly107Arg) in RDH5 and NM_000541.5:c.874C>T (p.Arg292Ter) in SAG were identified. These latter mutations have been reported previously, but not in the Pakistani population. CONCLUSIONS: Our findings expand the mutational spectrum of CSNB, in particular within the population of Southern Punjab.

Mutations in ARL2BP, a protein required for ciliary microtubule structure, cause syndromic male infertility in humans and mice.

Cilia are evolutionarily conserved hair-like structures with a wide spectrum of key biological roles, and their dysfunction has been linked to a growing class of genetic disorders, known collectively as ciliopathies. Many strides have been made towards deciphering the molecular causes for these diseases, which have in turn expanded the understanding of cilia and their functional roles. One recently-identified ciliary gene is ARL2BP, encoding the ADP-Ribosylation Factor Like 2 Binding Protein. In this study, we have identified multiple ciliopathy phenotypes associated with mutations in ARL2BP in human patients and in a mouse knockout model. Our research demonstrates that spermiogenesis is impaired, resulting in abnormally shaped heads, shortened and mis-assembled sperm tails, as well as in loss of axonemal doublets. Additional phenotypes in the mouse included enlarged ventricles of the brain and situs inversus. Mouse embryonic fibroblasts derived from knockout animals revealed delayed depolymerization of primary cilia. Our results suggest that ARL2BP is required for the structural maintenance of cilia as well as of the sperm flagellum, and that its deficiency leads to syndromic ciliopathy.

New clinical and molecular evidence linking mutations in ARSG to Usher syndrome type IV.

In murine and canine animal models, mutations in the Arylsulfatase G gene (ARSG) cause a particular lysosomal storage disorder characterized by neurological phenotypes. Recently, two variants in the same gene were found to be associated with an atypical form of Usher syndrome in humans, leading to visual and auditory impairment without the involvement of the central nervous system. In this study, we identified three novel pathogenic variants in ARSG, which segregated recessively with the disease in two families from Portugal. The probands were affected with retinitis pigmentosa and sensorineural hearing loss, generally with an onset of symptoms in their fourth decade of life. Functional experiments showed that these pathogenic variants abolish the sulfatase activity of the Arylsulfatase G enzyme and impede the appropriate lysosomal localization of the protein product, which appears to be retained in the endoplasmic reticulum. Our data enable to definitely confirm that different biallelic variants in ARSG cause a specific deaf-blindness syndrome, by abolishing the activity of the enzyme it encodes.

Herpes simplex encephalitis in adult patients with MASP-2 deficiency.

We report here two cases of Herpes simplex virus encephalitis (HSE) in adult patients with very rare, previously uncharacterized, non synonymous heterozygous G634R and R203W substitution in mannan-binding lectin serine protease 2 (MASP2), a gene encoding a key protease of the lectin pathway of the complement system. None of the 2 patients had variants in genes involved in the TLR3-interferon signaling pathway. Both MASP2 variants induced functional defects in vitro, including a reduced (R203W) or abolished (G634R) protein secretion, a lost capability to cleave MASP-2 precursor into its active form (G634R) and an in vivo reduced antiviral activity (G634R). In a murine model of HSE, animals deficient in mannose binding lectins (MBL, the main pattern recognition molecule associated with MASP-2) had a decreased survival rate and an increased brain burden of HSV-1 compared to WT C57BL/6J mice. Altogether, these data suggest that MASP-2 deficiency can increase susceptibility to adult HSE.

Phenotypic expansion of CACNA1C-associated disorders to include isolated neurological manifestations.

PURPOSE: CACNA1C encodes the alpha-1-subunit of a voltage-dependent L-type calcium channel expressed in human heart and brain. Heterozygous variants in CACNA1C have previously been reported in association with Timothy syndrome and long QT syndrome. Several case reports have suggested that CACNA1C variation may also be associated with a primarily neurological phenotype. METHODS: We describe 25 individuals from 22 families with heterozygous variants in CACNA1C, who present with predominantly neurological manifestations. RESULTS: Fourteen individuals have de novo, nontruncating variants and present variably with developmental delays, intellectual disability, autism, hypotonia, ataxia, and epilepsy. Functional studies of a subgroup of missense variants via patch clamp experiments demonstrated differential effects on channel function in vitro, including loss of function (p.Leu1408Val), neutral effect (p.Leu614Arg), and gain of function (p.Leu657Phe, p.Leu614Pro). The remaining 11 individuals from eight families have truncating variants in CACNA1C. The majority of these individuals have expressive language deficits, and half have autism. CONCLUSION: We expand the phenotype associated with CACNA1C variants to include neurodevelopmental abnormalities and epilepsy, in the absence of classic features of Timothy syndrome or long QT syndrome.

Immune deficiency, autoimmune disease and intellectual disability: A pleiotropic disorder caused by biallelic variants in the TPP2 gene.

Four individuals from two families presented with a multisystemic condition of suspected genetic origin that was diagnosed only after genome analysis. The main phenotypic features were immune system dysregulation (severe immunodeficiency with autoimmunity) and intellectual disability. The four individuals were found to be homozygous for a 4.4 Kb deletion removing exons 20-23 (NM_003291.4) of the TPP2 gene, predicting a frameshift with premature termination of the protein. The deletion was located on a shared chromosome 13 haplotype indicating a Swiss founder mutation. Tripeptidyl peptidase 2 (TPP2) is a protease involved in HLA/antigen complex processing and amino acid homeostasis. Biallelic variants in TPP2 have been described in 10 individuals with variable features including immune deficiency, autoimmune cytopenias, and intellectual disability or chronic sterile brain inflammation mimicking multiple sclerosis. Our observations further delineate this severe condition not yet included in the OMIM catalog. Timely recognition of TPP2 deficiency is crucial since (1) immune surveillance is needed and hematopoietic stem cell transplantation may be necessary, and (2) for provision of genetic counselling. Additionally, enzyme replacement therapy, as already established for TPP1 deficiency, might be an option in the future.