The evolution of gene expression and binding specificity of the largest transcription factor family in primates

The KRAB-containing zinc finger (KRAB-ZF) proteins represent the largest family of transcription factors (TFs) in humans, yet for the great majority, their function and specific genomic target remain unknown. However, it has been shown that a large fraction of these genes arose from segmental duplications, and that they have expanded in gene and zinc finger number throughout vertebrate evolution. To determine whether this expansion is linked to selective pressures acting on different domains, we have manually curated all KRAB-ZF genes present in the human genome together with their orthologous genes in three closely related species and assessed the evolutionary forces acting at the sequence level as well as on their expression profiles. We provide evidence that KRAB-ZFs can be separated into two categories according to the polymorphism present in their DNA-contacting residues. Those carrying a nonsynonymous single nucleotide polymorphism (SNP) in their DNA-contacting amino acids exhibit significantly reduced expression in all tissues, have emerged in a recent lineage, and seem to be less strongly constrained evolutionarily than those without such a polymorphism. This work provides evidence for a link between age of the TF, as well as polymorphism in their DNA-contacting residues and expression levels-both of which may be jointly affected by selection.

Computational study of transcription factor binding sites

Romain Fernand Pietro Groux

Any living organism contains a whole set of instructions encoded as genes on the DNA. This set of instructions contains all the necessary information that the organism will ever need, from its development to a mature individual to environment specific responses. Since all these instructions are not needed at the same time, the gene expression needs to be regulated. Eukaryotic genomes are stored inside nuclei as chromatin. The chromatin is the association of DNA with dedicated storage proteins - the histones - and the necessary machinery to regulate and express genes (RNA polymerases or RNAPs). In the nuclei, histones are assembled into octamers around which are wrapped ~148bp of DNA. This structure is known as the nucleosome. The repetition of nucleosomes along the genome allows to drastically compact the genome, eventually allowing to fit it inside the nucleus. However, this comes at the cost of rendering the DNA sequence inaccessible to DNA readers, such as the RNAPs and transcription factors (TFs). TFs are a class of proteins that have the remarkable property of recognizing and binding specific DNA sequences. More striking, each TF can recognize a multitude of different - but similar - DNA sequences providing TFs with a wide sequence specificity range. Eventually, this allows the cell to recruit TFs at dedicated locations in the genome called regulatory elements (REs). The action of TFs at REs is crucial to gene expression. Indeed, TFs are involved in many processes such as the opening of the chromatin structure or the recruitment of RNAPs. However if TFs can influence the chromatin structure, the opposite is also true as histones impede TF binding on DNA. Thus the regulation of genes relies on a subtle and complex crosstalk between the chromatin and TFs. To better understand how TFs and chromatin interact together to regulate gene expression, I lead several projects prospecting TF binding specificity and the chromatin structure at REs in human. First, I used ENCODE next generation sequencing (NGS) data to explore how TF binding influences the nearby nucleosome organization and the propensity of some TFs to bind together. The results suggest that regular nucleosome arrays are found near all TFs. They also point out two special cases. When CTCF binds with the cohesin complex, it seems to drive the nucleosome organization, which is a unique feature among all TFs investigated. Additionally I present evidence supporting that EBF1 is a pioneer factor - a special class of TFs able to bind nucleosome. Secondly, I developed several clustering algorithms and software to partition genomic regions according to NGS data and/or on their DNA sequences. These methods allow to discover important trends, for instance different nucleosome architectures . I illustrated the usefulness of these methods for the study of chromatin accessibility data and the identification of REs. Thirdly, I participated to the assessment of SMiLE-seq, a new microfluidic device that generates TF specificity data. The creation of TF specificity models and their comparison with other publicly available models demonstrated the value of SMiLE-seq to study TF specificity. Finally, I participated in the development of a software that predicts TF binding sites. A careful benchmarking suggested that this software is - at the time of writing - the best available software in terms of speed while showing other performances similar to its competitors.

EPFL2020

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

The Interplay between Cytidine Deaminases, HIV and Domesticated Genetic Parasites

"There is more virus in us than us in us". John Coffin's famous sentence illustrates that particular nucleic acid sequences related to exogenous viruses, called retrotransposons, constitute almost half of the human genome and largely exceeds the amount of protein-coding DNA (section 1.4). At the beginning of the past century, biologists realized the size of the genomes of various eukaryotes did not correlate with the level of complexity of these species. In extreme cases, certain species such as the fish L.paradoxa had 30-times more genomic DNA than do primates. In 1980, Francis Crick and Carmen Sapienza independently postulated that part of this once-called "junk DNA" must stem from self-replicating units, categorized as "selfish DNA" to distinguish them from more "altruistic" cellular genes. Indeed, while genes amplify in a population because of the selective advantage they confer to their host, the selfish DNA amplifies without making any contribution to the host phenotype and survives rather as a not-too-deleterious parasite. In early 2000's, the accomplishment of the human genome sequencing projects gave formal credence to Crick's and Sapienza's initial hypothesis. Retrotransposons replicate in the genome of their host by a copy-and-paste mechanism. Accordingly, any new retrotransposition event causes a genetic alteration potentially detrimental to the host. On an evolutionary perspective however, accumulation of retrotransposon-derived sequences can be beneficial to the species after these once "selfish" elements become "domesticated" to deserve more physiological roles. In human or mouse for instance, retrotransposons are known to cover a large spectrum of cellular functions, from gene regulation to antiviral defense (section 1.6). In order to maintain under tight control the balance between the good and adverse effects of retroelements, eukaryotes have evolved cellular defenses aimed at inhibiting the replication of such elements (section 1.7). Gene silencing by DNA methylation is a mechanism widely used in plants, fungi or vertebrates to control the expression of self-replicating elements. For instance in mammals, specific KRAB-zinc finger proteins can target retroviral-derived sequences and mediate their transcriptional silencing (section 3). The first aim of this study was to characterize the sensitivity of retroviruses to KRAB-mediated epigenetic silencing in varied chromosomal contexts (section 5). Since some retroviruses like HIV favor integration within transcriptionally active regions, these viruses may escape blockade by integrating chromosomal area refractory to epigenetic silencing. We found that instead, KRAB-mediated repression mechanism was not dependent upon integration site of a retroviral vector, hence the emergence of viruses escaping KRAB-mediated repression would be unlikely. Mammalian cells are also able to control the replication of retrotransposons at a post-transcriptional level via the APOBEC3 family of cytidine deaminases (section 2). Since certain retrotransposons are related to exogenous retroviruses, it is perhaps not surprising that APOBEC3 proteins can also inhibit a range of such elements including HIV. While there is only one APOBEC3 gene in the mouse genome, the family considerably expanded in other species such as primates were it comprises seven members (section 2.1). Although it is tempting to postulate that the expansion of the family in primates was driven by the strong selective pressure imposed by various genetic threats growing in number and complexity, it is not understood how only a few homologous proteins can target such a disparate group of replicating elements. The second aim of this study was to understand the mechanism by which different APOBEC3 proteins can target and inactivate particular self-perpetuating elements (section 4). We launched an integrative study focused on two human APOBEC3 homologues with distinct antiviral functions. APOBEC3A, a single-domain cytidine deaminase, can block efficiently the parvovirus adeno-associated virus type 2 (AAV-2) or the retrotransposon LINE-1 but cannot inhibit HIV replication. APOBEC3G is formed of two independent cytidine deaminase domains, only the C-terminal of which is catalytically active. APOBEC3G is endowed with restriction activity against HIV, but is largely ineffective against AAV-2 or LINE-1. We combined functional studies with structural modeling and phylogenetic analyses and found that restriction in each case involved similar regions at the surface of the protein and nearby the catalytic center, yet with distinct shapes and nucleic-acid interacting properties between the homologous proteins. While on APOBEC3G, the amino-acids in this region mediate intersubunit contacts between monomers and the predicted dimer interface shapes a groove that may facilitate binding to single-stranded RNA, on APOBEC3A however, it is corresponding residues on the monomer that form a DNA-binding groove important for editing. In summary, this study describes the nucleic-acid-interacting domains of APOBEC3A and APOBEC3G essential for their intrinsic restriction activities, and further suggests that structural variations in this domain may contribute to define the sets of targets inhibited by each APOBEC3 protein. APOBEC3 proteins may be considered as modular units, whose restriction function can evolve with emerging selective pressures by specific variations in their nucleic-acid binding domain, alternatively by fusion of two APOBEC3 proteins endowed with different functionalities. APOBEC3 proteins would efficiently protect the genome integrity of their host against highly variable genetic threats or "selfish DNA", in combination with other restriction factors such as the KRAB-zinc finger proteins. Contrary to epigenetic silencing however, post-transcriptional restriction as provided by APOBEC3 proteins may prohibit the deleterious effect of retrotransposition, without altering other retroelement functions putatively important for the physiology of the cell.

EPFL2010