The Characteristics of Heterozygous Protein Truncating Variants in the Human Genome

Sequencing projects have identified large numbers of rare stop-gain and frameshift variants in the human genome. As most of these are observed in the heterozygous state, they test a gene's tolerance to haploinsufficiency and dominant loss of function. We analyzed the distribution of truncating variants across 16,260 autosomal protein coding genes in 11,546 individuals. We observed 39,893 truncating variants affecting 12,062 genes, which significantly differed from an expectation of 12,916 genes under a model of neutral de novo mutation (p

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

KRAB zinc-finger proteins and their transposable element targets: between antagonism and cooperation

Jonas Caspar De Tribolet-Hardy

There are 377 Krüppel-associated box (KRAB) domain-containing zinc finger proteins (KZFPs) in the human genome, making them the largest family of transcription factors. KZFPs are defined by a N-terminal KRAB domain and several zinc-finger domains arranged in an array at the C-terminus of the proteins. The zinc-finger domains each form sequence specific interactions with double stranded DNA, allowing the zinc-finger array to target specific genomic sequences. The KRAB domain, through its interaction with the protein TRIM28 (also known as KAP1) allows for the stable silencing of transcription in a genomic region. Together these two domains allow KZFPs to bind specific regions of the genome and generally lead to the formation of heterochromatin and silencing of transcription allthough other functions for some KZFPs have been recently reported. KZFPs usually target transposable elements (TEs), mobile genomic elements that can move through the genome either by cut-and-paste or copy-paste mechanisms and make up almost half of the human genome. Different KZFPs bind to specific regions of TEs, limiting their expression and thus mitigating some of the threat they pose to human development, allowing both the TE and its host to survive. In recent years it became apparent that both certain TEs and KZFPs have roles beyond the previously described threat and mitigation of that threat. TEs harbor gene regulatory regions allowing, through their spread, for a dissemination of those regions through the genome and for new gene regulatory networks to form. KZFPs can affect the transcription of genes located in proximity to their binding sites leading to changes in gene regulation as well. A multitude of regulatory roles involving KZFPs and TEs have been described in recent years making them an exciting field of study. A major hurdle in understanding both KZFPs and TEs is to identify the binding sites of KZFPs, as they reveal both the targets of the KZFP and how, if at all, a TE is regulated by KZFPs. To do so, we expanded on previous efforts and aimed to perform chromatin immunoprecipitation followed by massively parallel sequencing (ChIP-seq) on every member of the human KZFP family. Here we present ChIP-seq experiments for 110 new KZFPs, which together with previously published data, allow us to study the binding of almost all human KZFPs (95%). The entirety of the data was analysed together to generate a coherent dataset and results for each KZFP are made available to the public on our web portal (https://tronoapps.epfl.ch/web/krabopedia/). The identified binding sites allowed us to witness the adaptation of the KZFP family to new TEs and showed how targeted sequences shift after segmental gene duplication events. Furthermore, we could corroborate that several KZFPs target the same TE subfamilies in a seemingly redundant fashion and show that unexpectedly these KZFPs arose independently in different genomic locations. These results represent a valuable tool for anyone studying KZFPs and should aid in the pursuit of both understanding KZFPs and TEs.

EPFL2022

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.