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.

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.

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.

GNB1-Related Rod-Cone Dystrophy: A Case Report

Introduction

The GNB1 (guanine nucleotide-binding protein, β1) gene encodes for the ubiquitous β1 subunit of heterotrimeric G proteins, which are associated with G-protein-coupled receptors (GPCRs). GNB1 mutations cause a neurodevelopmental disorder characterized by a broad clinical spectrum. A novel variant has recently been confirmed in a case of rod-cone dystrophy.

Case Presentation

We describe the second confirmed case of a classical rod-cone dystrophy associated with a mutation located in exon 6 of GNB1 [NM_002074.5:c.217G>C, p.(Ala73Pro)] in a 56-year-old patient also presenting mild intellectual disability, attention deficit/hyperactivity disorder, and truncal obesity.

Conclusion

This paper confirms the role of GNB1 in the pathogenesis of a classic rod-cone dystrophy and highlights the importance of including this gene in the genetic analysis panel for inherited retinal diseases.

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.

TBC1D32 variants disrupt retinal ciliogenesis and cause retinitis pigmentosa

Retinitis pigmentosa (RP) is the most common inherited retinal disease (IRD) and is characterized by photoreceptor degeneration and progressive vision loss. We report 4 patients presenting with RP from 3 unrelated families with variants in TBC1D32, which to date has never been associated with an IRD. To validate TBC1D32 as a putative RP causative gene, we combined Xenopus in vivo approaches and human induced pluripotent stem cell-derived (iPSC-derived) retinal models. Our data showed that TBC1D32 was expressed during retinal development and that it played an important role in retinal pigment epithelium (RPE) differentiation. Furthermore, we identified a role for TBC1D32 in ciliogenesis of the RPE. We demonstrated elongated ciliary defects that resulted in disrupted apical tight junctions, loss of functionality (delayed retinoid cycling and altered secretion balance), and the onset of an epithelial-mesenchymal transition-like phenotype. Last, our results suggested photoreceptor differentiation defects, including connecting cilium anomalies, that resulted in impaired trafficking to the outer segment in cones and rods in TBC1D32 iPSC-derived retinal organoids. Overall, our data highlight a critical role for TBC1D32 in the retina and demonstrate that TBC1D32 mutations lead to RP. We thus identify TBC1D32 as an IRD-causative gene.

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.

In-depth molecular profiling of an intronic GNAO1 mutant as the basis for personalized high-throughput drug screening

BACKGROUND

The GNAO1 gene, encoding the major neuronal G protein Gαo, is mutated in a subset of pediatric encephalopathies. Most such mutations consist of missense variants.

METHODS

In this study, we present a precision medicine workflow combining next-generation sequencing (NGS) diagnostics, molecular etiology analysis, and personalized drug discovery.

FINDINGS

We describe a patient carrying a de novo intronic mutation (NM_020988.3:c.724-8G>A), leading to epilepsy-negative encephalopathy with motor dysfunction from the second decade. Our data show that this mutation creates a novel splice acceptor site that in turn causes an in-frame insertion of two amino acid residues, Pro-Gln, within the regulatory switch III region of Gαo. This insertion misconfigures the switch III loop and creates novel interactions with the catalytic switch II region, resulting in increased GTP uptake, defective GTP hydrolysis, and aberrant interactions with effector proteins. In contrast, intracellular localization, Gβγ interactions, and G protein-coupled receptor (GPCR) coupling of the Gαo[insPQ] mutant protein remain unchanged.

CONCLUSIONS

This in-depth analysis characterizes the heterozygous c.724-8G>A mutation as partially dominant negative, providing clues to the molecular etiology of this specific pathology. Further, this analysis allows us to establish and validate a high-throughput screening platform aiming at identifying molecules that could correct the aberrant biochemical functions of the mutant Gαo.

FUNDING

This work was supported by the Joint Seed Money Funding scheme between the University of Geneva and the Hebrew University of Jerusalem.

The first genetic landscape of inherited retinal dystrophies in Portuguese patients identifies recurrent homozygous mutations as a frequent cause of pathogenesis

Inherited retinal diseases (IRDs) are a group of ocular conditions characterized by an elevated genetic and clinical heterogeneity. They are transmitted almost invariantly as monogenic traits. However, with more than 280 disease genes identified so far, association of clinical phenotypes with genotypes can be very challenging, and molecular diagnosis is essential for genetic counseling and correct management of the disease. In addition, the prevalence and the assortment of IRD mutations are often population-specific. In this work, we examined 230 families from Portugal, with individuals suffering from a variety of IRD diagnostic classes (270 subjects in total). Overall, we identified 157 unique mutations (34 previously unreported) in 57 distinct genes, with a diagnostic rate of 76%. The IRD mutational landscape was, to some extent, different from those reported in other European populations, including Spanish cohorts. For instance, the EYS gene appeared to be the most frequently mutated, with a prevalence of 10% among all IRD cases. This was, in part, due to the presence of a recurrent and seemingly founder mutation involving the deletion of exons 13 and 14 of this gene. Moreover, our analysis highlighted that as many as 51% of our cases had mutations in a homozygous state. To our knowledge, this is the first study assessing a cross-sectional genotype-phenotype landscape of IRDs in Portugal. Our data reveal a rather unique distribution of mutations, possibly shaped by a small number of rare ancestral events that have now become prevalent alleles in patients.

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.

Biallelic variants in coenzyme Q10 biosynthesis pathway genes cause a retinitis pigmentosa phenotype

The aim of this study was to investigate coenzyme Q10 (CoQ₁₀) biosynthesis pathway defects in inherited retinal dystrophy. Individuals affected by inherited retinal dystrophy (IRD) underwent exome or genome sequencing for molecular diagnosis of their condition. Following negative IRD gene panel analysis, patients carrying biallelic variants in CoQ₁₀ biosynthesis pathway genes were identified. Clinical data were collected from the medical records. Haplotypes harbouring the same missense variant were characterised from family genome sequencing (GS) data and direct Sanger sequencing. Candidate splice variants were characterised using Oxford Nanopore Technologies single molecule sequencing. The CoQ₁₀ status of the human plasma was determined in some of the study patients. 13 individuals from 12 unrelated families harboured candidate pathogenic genotypes in the genes: PDSS1, COQ2, COQ4 and COQ5. The PDSS1 variant c.589 A > G was identified in three affected individuals from three unrelated families on a possible ancestral haplotype. Three variants (PDSS1 c.468-25 A > G, PDSS1 c.722-2 A > G, COQ5 c.682-7 T > G) were shown to lead to cryptic splicing. 6 affected individuals were diagnosed with non-syndromic retinitis pigmentosa and 7 had additional clinical findings. This study provides evidence of CoQ₁₀ biosynthesis pathway gene defects leading to non-syndromic retinitis pigmentosa in some cases. Intronic variants outside of the canonical splice-sites represent an important cause of disease. RT-PCR nanopore sequencing is effective in characterising these splice defects.

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.

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.

Fast and highly sensitive full-length single-cell RNA sequencing using FLASH-seq

We present FLASH-seq (FS), a full-length single-cell RNA sequencing (scRNA-seq) method with increased sensitivity and reduced hands-on time compared to Smart-seq3. The entire FS protocol can be performed in ~4.5 hours, is simple to automate and can be easily miniaturized to decrease resource consumption. The FS protocol can also use unique molecular identifiers (UMIs) for molecule counting while displaying reduced strand-invasion artifacts. FS will be especially useful for characterizing gene expression at high resolution across multiple samples.

FLASH-seq speeds up high-sensitivity scRNA-seq while sparing resources.

Identification of New Vulnerabilities in Conjunctival Melanoma Using Image-Based High Content Drug Screening

Recent evidence suggests that numerous similarities exist between the genomic landscapes of both conjunctival and cutaneous melanoma. Since alterations of several components of the MAP kinases, PI3K/mTOR, and cell cycle pathways have been reported in conjunctival melanoma, we decided to assess the sensitivity of conjunctival melanoma to targeted inhibition mostly of kinase inhibitors. A high content drug screening assay based on automated fluorescence microscopy was performed in three conjunctival melanoma cell lines with different genomic backgrounds with 489 kinase inhibitors and 53 other inhibitors. IC50 and apoptosis induction were respectively assessed for 53 and 48 compounds. The genomic background influenced the response to MAK and PI3K/mTOR inhibition, more specifically cell lines with BRAF V600E mutations were more sensitive to BRAF/MEK inhibition, while CRMM2 bearing the NRASQ61L mutation was more sensitive to PI3k/mTOR inhibition. All cell lines demonstrated sensitivity to cell cycle inhibition, being more pronounced in CRMM2, especially with polo-like inhibitors. Our data also revealed new vulnerabilities to Hsp90 and Src inhibition. This study demonstrates that the genomic background partially influences the response to targeted therapy and uncovers a large panel of potential vulnerabilities in conjunctival melanoma that may expand available options for the management of this tumor.

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.

Genetic analysis of consanguineous Pakistani families with congenital stationary night blindness

AIM

Congenital stationary night blindness (CSNB) is a rare, largely non progressive, inherited retinal disorder that can be clinically classified on the basis of fundus and electroretinogram (ERG) abnormalities.

METHODS

We analyzed four large consanguineous families from the Southern Punjab region of Pakistan including multiple individuals affected with CSNB. Exome sequencing (ES) 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 variants 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.

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.

Whole exome sequencing in 17 consanguineous Iranian pedigrees expands the mutational spectrum of inherited retinal dystrophies

Inherited retinal dystrophies (IRDs) constitute one of the most heterogeneous groups of Mendelian human disorders. Using autozygome-guided next-generation sequencing methods in 17 consanguineous pedigrees of Iranian descent with isolated or syndromic IRD, we identified 17 distinct genomic variants in 11 previously-reported disease genes. Consistent with a recessive inheritance pattern, as suggested by pedigrees, variants discovered in our study were exclusively bi-allelic and mostly in a homozygous state (in 15 families out of 17, or 88%). Out of the 17 variants identified, 5 (29%) were never reported before. Interestingly, two mutations (GUCY2D:c.564dup, p.Ala189ArgfsTer130 and TULP1:c.1199G > A, p.Arg400Gln) were also identified in four separate pedigrees (two pedigrees each). In addition to expanding the mutational spectrum of IRDs, our findings confirm that the traditional practice of endogamy in the Iranian population is a prime cause for the appearance of IRDs.

NGS and phenotypic ontology-based approaches increase the diagnostic yield in syndromic retinal diseases

Syndromic retinal diseases (SRDs) are a group of complex inherited systemic disorders, with challenging molecular underpinnings and clinical management. Our main goal is to improve clinical and molecular SRDs diagnosis, by applying a structured phenotypic ontology and next-generation sequencing (NGS)-based pipelines. A prospective and retrospective cohort study was performed on 100 probands with an a priori diagnosis of non-Usher SRDs, using available clinical data, including Human Phenotype Ontology annotation, and further classification into seven clinical categories (ciliopathies, specific syndromes and five others). Retrospective molecular diagnosis was assessed using different molecular and bioinformatic methods depending on availability. Subsequently, uncharacterized probands were prospectively screened using other NGS approaches to extend the number of analyzed genes. After phenotypic classification, ciliopathies were the most common SRD (35%). A global characterization rate of 52% was obtained, with six cases incompletely characterized for a gene that partially explained the phenotype. An improved characterization rate was achieved addressing prospective cases (83%) and well-recognizable syndrome (62%) subgroups. The 27% of the fully characterized cases were reclassified into a different clinical category after identification of the disease-causing gene. Clinical-exome sequencing is the most appropriate first-tier approach for prospective cases, whereas whole-exome sequencing and bioinformatic reanalysis increases the diagnosis of uncharacterized retrospective cases to 45%, mostly those with unspecific symptoms. Our study describes a comprehensive approach to SRDs in daily clinical practice and the importance of thorough clinical assessment and selection of the most appropriate molecular test to be used to solve these complex cases and elucidate novel associations.

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.

CEP78 functions downstream of CEP350 to control biogenesis of primary cilia by negatively regulating CP110 levels

CEP78 is a centrosomal protein implicated in ciliogenesis and ciliary length control, and mutations in the CEP78 gene cause retinal cone-rod dystrophy associated with hearing loss. However, the mechanism by which CEP78 affects cilia formation is unknown. Based on a recently discovered disease-causing CEP78 p.L150S mutation, we identified the disease-relevant interactome of CEP78. We confirmed that CEP78 interacts with the EDD1-DYRK2-DDB1^VPRBP E3 ubiquitin ligase complex, which is involved in CP110 ubiquitination and degradation, and identified a novel interaction between CEP78 and CEP350 that is weakened by the CEP78^L150S mutation. We show that CEP350 promotes centrosomal recruitment and stability of CEP78, which in turn leads to centrosomal recruitment of EDD1. Consistently, cells lacking CEP78 display significantly increased cellular and centrosomal levels of CP110, and depletion of CP110 in CEP78-deficient cells restored ciliation frequency to normal. We propose that CEP78 functions downstream of CEP350 to promote ciliogenesis by negatively regulating CP110 levels via an EDD1-dependent mechanism.

Transcriptomic Signature Differences Between SARS-CoV-2 and Influenza Virus Infected Patients

The reason why most individuals with COVID-19 have relatively limited symptoms while other develop respiratory distress with life-threatening complications remains unknown. Increasing evidence suggests that COVID-19 associated adverse outcomes mainly rely on dysregulated immunity. Here, we compared transcriptomic profiles of blood cells from 103 patients with different severity levels of COVID-19 with that of 27 healthy and 22 influenza-infected individuals. Data provided a complete overview of SARS-CoV-2-induced immune signature, including a dramatic defect in IFN responses, a reduction of toxicity-related molecules in NK cells, an increased degranulation of neutrophils, a dysregulation of T cells, a dramatic increase in B cell function and immunoglobulin production, as well as an important over-expression of genes involved in metabolism and cell cycle in patients infected with SARS-CoV-2 compared to those infected with influenza viruses. These features also differed according to COVID-19 severity. Overall and specific gene expression patterns across groups can be visualized on an interactive website (https://bix.unil.ch/covid/). Collectively, these transcriptomic host responses to SARS-CoV-2 infection are discussed in the context of current studies, thereby improving our understanding of COVID-19 pathogenesis and shaping the severity level of COVID-19.

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.

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.

Homozygosity mapping is a useful tool for identifying candidate mutations in recessive conditions, however application to next generation sequencing data has been sub-optimal. Here, the authors present AutoMap, which efficiently identifies runs of homozygosity in whole exome/genome sequencing data.

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.

Management of Full-Thickness Macular Hole in A Genetically Confirmed Case with Usher Syndrome

INTRODUCTION

Full-thickness macular hole (FTMH) formation is rarely seen in patients with retinitis pigmentosa (RP) and can have an adverse impact on their residual visual function. The underlying mechanisms are unknown, and clinical experience is limited regarding surgical outcomes. Here, we describe the surgical management of FTMH in a young patient with genetically confirmed Usher syndrome, the most common form of syndromic RP.

CASE REPORT

A 28-year-old woman presented with blurred vision in her right eye (RE). She had a history of RP and bilateral hearing impairment since childhood. Fundoscopy and spectral-domain optical coherence tomography revealed a FTMH in the RE along with typical RP features bilaterally. After pars plana vitrectomy (PPV) with internal limiting membrane peel and gas tamponade, the FTMH closed. Six months after PPV the patient underwent cataract surgery in the affected eye, and the visual acuity remained stable compared to baseline. The clinical diagnosis of Usher syndrome was genetically confirmed by whole exome sequencing (WES), which revealed the presence of two pathogenic nucleotide variants in trans (compound heterozygosity) in the gene USH2A.

CONCLUSION

We report a rare case of successful closure of a FTMH in a patient with Usher syndrome. Surgical treatment of FTMH can help preserve the central vision in RP patients, whose peripheral vision is severely affected.

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.

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.

Human herpes virus type 1 (HSV-1) infects a large number of individuals during their life, with manifestations usually limited to mild and self-limiting inflammation of the oral mucosa (cold sore). However, HSV-1 can cause a very severe disease of the brain called Herpes simplex encephalitis (HSE) in 1 out of 250’000–500’000 individuals per year. The reasons why HSV-1 can cause such a devastating disease in a very limited number of individuals are unknown. Increasing evidence suggests that susceptibility to HSE in children can results from genetic variations in the immune system, in particular in a viral detection pathway called the Toll-like receptor 3 (TLR3)–interferon (IFN) axis. Fewer data are available to explain HSE in adult patients. Here, we describe two adult patients with HSE who carry mutations in a gene called mannan-binding lectin serine protease 2 (MASP2), which is part of an immune pathway different from the TLR3-IFN axis, called the lectin pathway of the complement system. We demonstrate that MASP2 mutations induce functional defects in immune defense against HSV-1 that prevent viral replication. Mice deficient in the lectin pathway have higher mortality compared to wild-type mice after HSV-1 infection. Altogether, our study suggests that susceptibility to HSE in adults relies of immune deficiencies that are different from those causing HSE in children.

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.

The flagellated tails of sperm cells require a stringent developmental process that is essential for motility and fertility. The components that comprise the sperm tail assemble in regulated steps with protein processing, transport, and structural assembly dependent on each other for sperm tail maturity. In this work, we have identified ARL2BP, a previously retinal-associated protein, to be essential for sperm tail development and assembly. We show that without functional ARL2BP in humans or mice, sperm tails fail to develop, starting with the assembly of the core microtubular structure within the tail. Loss of ARL2BP also effects other ciliated cells, indicating a unique role for ARL2BP in ciliary microtubule formation. This research on ARL2BP provides further understanding on the links between vision and fertility. This work also demonstrates how genomic studies for human patients and murine models can coincide to provide greater insight into disease.

A frequent variant in the Japanese population determines quasi-Mendelian inheritance of rare retinal ciliopathy

Hereditary retinal degenerations (HRDs) are Mendelian diseases characterized by progressive blindness and caused by ultra-rare mutations. In a genomic screen of 331 unrelated Japanese patients, we identify a disruptive Alu insertion and a nonsense variant (p.Arg1933*) in the ciliary gene RP1, neither of which are rare alleles in Japan. p.Arg1933* is almost polymorphic (frequency = 0.6%, amongst 12,000 individuals), does not cause disease in homozygosis or heterozygosis, and yet is significantly enriched in HRD patients (frequency = 2.1%, i.e., a 3.5-fold enrichment; p-value = 9.2 × 10-5). Familial co-segregation and association analyses show that p.Arg1933* can act as a Mendelian mutation in trans with the Alu insertion, but might also associate with disease in combination with two alleles in the EYS gene in a non-Mendelian pattern of heredity. Our results suggest that rare conditions such as HRDs can be paradoxically determined by relatively common variants, following a quasi-Mendelian model linking monogenic and complex inheritance.

Reciprocal modulation of mesenchymal stem cells and tumor cells promotes lung cancer metastasis

Metastasis is a multi-step process in which direct crosstalk between cancer cells and their microenvironment plays a key role. Here, we assessed the effect of paired tumor-associated and normal lung tissue mesenchymal stem cells (MSCs) on the growth and dissemination of primary human lung carcinoma cells isolated from the same patients. We show that the tumor microenvironment modulates MSC gene expression and identify a four-gene MSC signature that is functionally implicated in promoting metastasis. We also demonstrate that tumor-associated MSCs induce the expression of genes associated with an aggressive phenotype in primary lung cancer cells and selectively promote their dissemination rather than local growth. Our observations provide insight into mechanisms by which the stroma promotes lung cancer metastasis.

•

Distinct gene expression profiles distinguish normal lung and tumor-associated MSCs.

•

MSCs induce EMT- and hypoxia-related genes in primary tumor cells and promote their metastatic potential.

•

A 4-gene T-MSC signature is involved in MSC-induced metastasis promotion.

The tumor microenvironment, which includes mesenchymal stem cells (MSCs) among many other stromal cell types, plays a fundamental role in cancer metastasis. Although MSCs are suggested to participate in tumor progression, most studies thus far have been performed on bone marrow-derived MSCs and cancer cell lines. Using primary human pulmonary MSCs and paired lung cancer cells, we show that tumor cells modulate MSCs to acquire properties, including a four-gene signature, which allow them to promote tumor dissemination. Our results provide insight into the mutual cancer cell-stromal cell modulation that drives tumor dissemination.

Identification of Inherited Retinal Disease-Associated Genetic Variants in 11 Candidate Genes

Inherited retinal diseases (IRDs) display an enormous genetic heterogeneity. Whole exome sequencing (WES) recently identified genes that were mutated in a small proportion of IRD cases. Consequently, finding a second case or family carrying pathogenic variants in the same candidate gene often is challenging. In this study, we searched for novel candidate IRD gene-associated variants in isolated IRD families, assessed their causality, and searched for novel genotype-phenotype correlations. Whole exome sequencing was performed in 11 probands affected with IRDs. Homozygosity mapping data was available for five cases. Variants with minor allele frequencies ≤ 0.5% in public databases were selected as candidate disease-causing variants. These variants were ranked based on their: (a) presence in a gene that was previously implicated in IRD; (b) minor allele frequency in the Exome Aggregation Consortium database (ExAC); (c) in silico pathogenicity assessment using the combined annotation dependent depletion (CADD) score; and (d) interaction of the corresponding protein with known IRD-associated proteins. Twelve unique variants were found in 11 different genes in 11 IRD probands. Novel autosomal recessive and dominant inheritance patterns were found for variants in Small Nuclear Ribonucleoprotein U5 Subunit 200 (SNRNP200) and Zinc Finger Protein 513 (ZNF513), respectively. Using our pathogenicity assessment, a variant in DEAH-Box Helicase 32 (DHX32) was the top ranked novel candidate gene to be associated with IRDs, followed by eight medium and lower ranked candidate genes. The identification of candidate disease-associated sequence variants in 11 single families underscores the notion that the previously identified IRD-associated genes collectively carry > 90% of the defects implicated in IRDs. To identify multiple patients or families with variants in the same gene and thereby provide extra proof for pathogenicity, worldwide data sharing is needed.

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.

Likelihood-Free Inference in High-Dimensional Models

Methods that bypass analytical evaluations of the likelihood function have become an indispensable tool for statistical inference in many fields of science. These so-called likelihood-free methods rely on accepting and rejecting simulations based on summary statistics, which limits them to low-dimensional models for which the value of the likelihood is large enough to result in manageable acceptance rates. To get around these issues, we introduce a novel, likelihood-free Markov chain Monte Carlo (MCMC) method combining two key innovations: updating only one parameter per iteration and accepting or rejecting this update based on subsets of statistics approximately sufficient for this parameter. This increases acceptance rates dramatically, rendering this approach suitable even for models of very high dimensionality. We further derive that for linear models, a one-dimensional combination of statistics per parameter is sufficient and can be found empirically with simulations. Finally, we demonstrate that our method readily scales to models of very high dimensionality, using toy models as well as by jointly inferring the effective population size, the distribution of fitness effects (DFE) of segregating mutations, and selection coefficients for each locus from data of a recent experiment on the evolution of drug resistance in influenza.

Characteristic bimodal profiles of RNA polymerase II at thousands of active mammalian promoters

Background

In mammals, ChIP-seq studies of RNA polymerase II (PolII) occupancy have been performed to reveal how recruitment, initiation and pausing of PolII may control transcription rates, but the focus is rarely on obtaining finely resolved profiles that can portray the progression of PolII through sequential promoter states.

Results

Here, we analyze PolII binding profiles from high-coverage ChIP-seq on promoters of actively transcribed genes in mouse and humans. We show that the enrichment of PolII near transcription start sites exhibits a stereotypical bimodal structure, with one peak near active transcription start sites and a second peak 110 base pairs downstream from the first. Using an empirical model that reliably quantifies the spatial PolII signal, gene by gene, we show that the first PolII peak allows for refined positioning of transcription start sites, which is corroborated by mRNA sequencing. This bimodal signature is found both in mouse and humans. Analysis of the pausing-related factors NELF and DSIF suggests that the downstream peak reflects widespread pausing at the +1 nucleosome barrier. Several features of the bimodal pattern are correlated with sequence features such as CpG content and TATA boxes, as well as the histone mark H3K4me3.

Conclusions

We thus show how high coverage DNA sequencing experiments can reveal as-yet unnoticed bimodal spatial features of PolII accumulation that are frequent at individual mammalian genes and reminiscent of transcription initiation and pausing. The initiation-pausing hypothesis is corroborated by evidence from run-on sequencing and immunoprecipitation in other cell types and species.

Barium sulphate and soft-tissue radiology: allying the old and the new for the investigation of animal cutaneous microcirculation

The study of microcirculation using angiography is essential to the advancement of flap and angiogenesis research in plastic surgery. Until the mid-1980s, barium sulphate was the most commonly used contrast material, although it did not provide optimal visualisation of the vascular tree. In 1986, a new technique using lead oxide was proposed, which permitted very high-quality imaging and rapidly became the technique of choice, despite its high toxicity. We reconsider the former technique of barium-sulphate injection and combine it with soft-tissue radiology using mammographic film to achieve a radiological definition similar to that obtained with lead oxide, and discuss the advantages and disadvantages of the two methods. We conclude that barium sulphate and the use of mammographic film is an accurate, simple and non-toxic method of analysing the cutaneous circulation in small animals.