Where are the missing gene defects in inherited retinal disorders? Intronic and synonymous variants contribute at least to 4% of CACNA1F-mediated inherited retinal disorders

Inherited retinal disorders (IRD) represent clinically and genetically heterogeneous diseases. To date, pathogenic variants have been identified in similar to 260 genes. Albeit that many genes are implicated in IRD, for 30-50% of the cases, the gene defect is unknown. These cases may be explained by novel gene defects, by overlooked structural variants, by variants in intronic, promoter or more distant regulatory regions, and represent synonymous variants of known genes contributing to the dysfunction of the respective proteins. Patients with one subgroup of IRD, namely incomplete congenital stationary night blindness (icCSNB), show a very specific phenotype. The major cause of this condition is the presence of a hemizygous pathogenic variant in CACNA1F. A comprehensive study applying direct Sanger sequencing of the gene-coding regions, exome and genome sequencing applied to a large cohort of patients with a clinical diagnosis of icCSNB revealed indeed that seven of the 189 CACNA1F-related cases have intronic and synonymous disease-causing variants leading to missplicing as validated by minigene approaches. These findings highlight that gene-locus sequencing may be a very efficient method in detecting disease-causing variants in clinically well-characterized patients with a diagnosis of IRD, like icCSNB.

Computational analysis of ultraconserved non-coding DNA elements in vertebrate genomes

Slavica Dimitrieva Janeva

Recent advancements in DNA sequencing technologies, accompanied by parallel development of sophisticated computational methods, offer an unprecedented opportunity to explore and study the genetic material of many organisms. However, despite the tremendous progress in the field of genomics, our understanding of the genetic information encoded by the DNA is far from complete. Whole-genome comparisons of vertebrate species have uncovered the existence of non-coding DNA sequences that exhibit exceptionally high levels of similarity over hundreds of base pairs across distant species, referred to as ultraconserved non-coding elements (UCNEs). The existence of such sequences represents one of the biggest mysteries in current biology. Many UCNEs exhibit even stronger conservation than protein coding regions, but to date no molecular mechanism has been described that would require such a high degree of conservation over such long sequences. This thesis focuses on computational analysis of UCNEs across vertebrate genomes, with an ultimate goal to provide insights into the functioning of these elements and to unravel the reasons for their extreme conservation. Using computational approaches we studied the salient characteristics of UCNEs, and explored their genomic environment and their organization across genomes. Like previous studies, we observed that UCNEs are organized as large clusters around essential developmental genes, forming genomic regulatory blocks. Then we sought to understand the reasons behind this strong clustering of UCNEs. The clustering could reflect functional cooperativity between UCNEs. Alternatively, if each UCNE acts independently, their high concentration near developmental genes could merely reflect the extreme regulatory complexity of these genes. In a special setting of genomic context analysis, we analyzed the fate of UCNEs in teleost fish genomes that were subjected to an additional round of whole-genome duplication. We found that in most cases all UCNEs of a block were retained in one copy only, but together on the same chromosome. Conversely, the corresponding target genes were often retained in two copies, one completely devoid of UCNEs. Our results suggested that UCNEs of a cluster function in a highly cooperative manner. We propose that a multitude of cooperative cis-interactions between UCNEs is the reason for their extreme sequence conservation. Our work presents a novel hypothesis about the reasons for UCNE conservation and their mode of action, which is yet to be confirmed experimentally. The results from our study are made accessible through a new web resource called UCNEbase (http://ccg.vital-it.ch/UCNEbase). UCNEbase provides informa- tion about the genomic organization of UCNEs and their association with developmental regulatory genes, and enables all data to be explored in a genome browser environment. As part of a side project, we developed a new sparsified version of algorithms for computational prediction of RNA secondary structures and introduced a new program sibRNAfold. We performed a thorough analysis of the time complexity of sparsified RNA folding algorithms covering a large parameter space and empirically demonstrated that sparse RNA folding has cubic time complexity rather than quadratic, as claimed previously.

EPFL2013

Analysis of Stop-Gain and Frameshift Variants in Human Innate Immunity Genes

Istvan Bartha, Jacques Fellay, Amalio Telenti, Ioannis Xenarios

Loss-of-function variants in innate immunity genes are associated with Mendelian disorders in the form of primary immunodeficiencies. Recent resequencing projects report that stop-gains and frameshifts are collectively prevalent in humans and could be responsible for some of the inter-individual variability in innate immune response. Current computational approaches evaluating loss-of-function in genes carrying these variants rely on gene-level characteristics such as evolutionary conservation and functional redundancy across the genome. However, innate immunity genes represent a particular case because they are more likely to be under positive selection and duplicated. To create a ranking of severity that would be applicable to innate immunity genes we evaluated 17,764 stop-gain and 13,915 frameshift variants from the NHLBI Exome Sequencing Project and 1,000 Genomes Project. Sequence-based features such as loss of functional domains, isoform-specific truncation and nonsense-mediated decay were found to correlate with variant allele frequency and validated with gene expression data. We integrated these features in a Bayesian classification scheme and benchmarked its use in predicting pathogenic variants against Online Mendelian Inheritance in Man (OMIM) disease stop-gains and frameshifts. The classification scheme was applied in the assessment of 335 stop-gains and 236 frameshifts affecting 227 interferon-stimulated genes. The sequence-based score ranks variants in innate immunity genes according to their potential to cause disease, and complements existing gene-based pathogenicity scores. Specifically, the sequence-based score improves measurement of functional gene impairment, discriminates across different variants in a given gene and appears particularly useful for analysis of less conserved genes.

Public Library Science2014

Computational analysis of ultraconserved non-coding DNA elements in vertebrate genomes

Slavica Dimitrieva Janeva

EPFL2013