An evolutionary arms race between KRAB zinc finger genes 91/93 and SVA/L1 retrotransposons

Throughout evolution primate genomes have been modified by waves of retrotransposon insertions. For each wave, the host eventually finds a way to repress retrotransposon transcription and prevent further insertions. In mouse embryonic stem cells, transcriptional silencing of retrotransposons requires KAP1 (also known as TRIM28) and its repressive complex, which can be recruited to target sites by KRAB zinc-finger (KZNF) proteins such as murine-specific ZFP809 which binds to integrated murine leukaemia virus DNA elements and recruits KAP1 to repress them. KZNF genes are one of the fastest growing gene families in primates and this expansion is hypothesized to enable primates to respond to newly emerged retrotransposons. However, the identity of KZNF genes battling retrotransposons currently active in the human genome, such as SINE-VNTR-Alu (SVA) and long interspersed nuclear element 1 (L1), is unknown. Here we show that two primate-specific KZNF genes rapidly evolved to repress these two distinct retrotransposon families shortly after they began to spread in our ancestral genome. ZNF91 underwent a series of structural changes 8–12 million years ago that enabled it to repress SVA elements. ZNF93 evolved earlier to repress the primate L1 lineage until ∼12.5 million years ago when the L1PA3-subfamily of retrotransposons escaped ZNF93’s restriction through the removal of the ZNF93-binding site. Our data support a model where KZNF gene expansion limits the activity of newly emerged retrotransposon classes, and this is followed by mutations in these retrotransposons to evade repression, a cycle of events that could explain the rapid expansion of lineage-specific KZNF genes.

[1]  J. V. Moran,et al.  Initial sequencing and analysis of the human genome. , 2001, Nature.

[2]  J. V. Moran,et al.  Dynamic interactions between transposable elements and their hosts , 2011, Nature Reviews Genetics.

[3]  M. Oshimura,et al.  Enhanced apoptosis during early neuronal differentiation in mouse ES cells with autosomal imbalance , 2009, Cell Research.

[4]  G. Swergold Identification, characterization, and cell specificity of a human LINE-1 promoter , 1990, Molecular and cellular biology.

[5]  Cole Trapnell,et al.  TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions , 2013, Genome Biology.

[6]  W. Huber,et al.  which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. MAnorm: a robust model for quantitative comparison of ChIP-Seq data sets , 2011 .

[7]  Aaron R. Quinlan,et al.  Bioinformatics Applications Note Genome Analysis Bedtools: a Flexible Suite of Utilities for Comparing Genomic Features , 2022 .

[8]  D. C. Hancks,et al.  Active human retrotransposons: variation and disease. , 2012, Current opinion in genetics & development.

[9]  A. Löytynoja,et al.  Phylogeny-Aware Gap Placement Prevents Errors in Sequence Alignment and Evolutionary Analysis , 2008, Science.

[10]  H. Kazazian,et al.  The Minimal Active Human SVA Retrotransposon Requires Only the 5′-Hexamer and Alu-Like Domains , 2012, Molecular and Cellular Biology.

[11]  Chris P. Ponting,et al.  Bioinformatics Discovery Note Meisetz and the Birth of the Krab Motif , 2022 .

[12]  Stéphane Boissinot,et al.  Molecular evolution and tempo of amplification of human LINE-1 retrotransposons since the origin of primates. , 2005, Genome research.

[13]  J. V. Moran,et al.  An actively retrotransposing, novel subfamily of mouse L1 elements , 1998, The EMBO journal.

[14]  Austin G Smith,et al.  Conversion of embryonic stem cells into neuroectodermal precursors in adherent monoculture , 2003, Nature Biotechnology.

[15]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[16]  W. Berger,et al.  Full-length human L1 insertions retain the capacity for high frequency retrotransposition in cultured cells. , 1999, Human molecular genetics.

[17]  Koichiro Tamura,et al.  MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. , 2013, Molecular biology and evolution.

[18]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[19]  D. Trono,et al.  Interplay of TRIM28 and DNA methylation in controlling human endogenous retroelements , 2014, Genome research.

[20]  A Klug,et al.  Design of polyzinc finger peptides with structured linkers. , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[21]  T. Borodina,et al.  Transcriptome analysis by strand-specific sequencing of complementary DNA , 2009, Nucleic acids research.

[22]  David Haussler,et al.  Thousands of human mobile element fragments undergo strong purifying selection near developmental genes , 2007, Proceedings of the National Academy of Sciences.

[23]  Paul Flicek,et al.  Latent Regulatory Potential of Human-Specific Repetitive Elements , 2013, Molecular cell.

[24]  J. V. Moran,et al.  Determination of L1 retrotransposition kinetics in cultured cells. , 2000, Nucleic acids research.

[25]  P. Bray-Ward,et al.  Emergence of the ZNF91 Krüppel-associated box-containing zinc finger gene family in the last common ancestor of anthropoidea. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[26]  Laurie Gordon,et al.  A comprehensive catalog of human KRAB-associated zinc finger genes: insights into the evolutionary history of a large family of transcriptional repressors. , 2006, Genome research.

[27]  G. Benson,et al.  Tandem repeats finder: a program to analyze DNA sequences. , 1999, Nucleic acids research.

[28]  Helen M. Rowe,et al.  KAP1 controls endogenous retroviruses in embryonic stem cells , 2010, Nature.

[29]  K. Katoh,et al.  MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability , 2013, Molecular biology and evolution.

[30]  Mona Singh,et al.  Predicting DNA recognition by Cys2His2 zinc finger proteins , 2009, Bioinform..

[31]  D. Trono,et al.  Evolutionally dynamic L1 regulation in embryonic stem cells , 2014, Genes & development.

[32]  H. Kazazian Mobile Elements: Drivers of Genome Evolution , 2004, Science.

[33]  Tom H. Pringle,et al.  The human genome browser at UCSC. , 2002, Genome research.

[34]  S. Goff,et al.  Embryonic stem cells use ZFP809 to silence retroviral DNAs , 2009, Nature.

[35]  Koichiro Tamura,et al.  Estimating divergence times in large molecular phylogenies , 2012, Proceedings of the National Academy of Sciences.

[36]  M. Batzer,et al.  The impact of retrotransposons on human genome evolution , 2009, Nature Reviews Genetics.

[37]  Jerilyn A. Walker,et al.  SVA elements: a hominid-specific retroposon family. , 2005, Journal of molecular biology.

[38]  Clifford A. Meyer,et al.  Model-based Analysis of ChIP-Seq (MACS) , 2008, Genome Biology.

[39]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[40]  David Haussler,et al.  The UCSC Known Genes , 2006, Bioinform..

[41]  M. Nei,et al.  A Simple Method for Estimating and Testing Minimum-Evolution Trees , 1992 .

[42]  S. Goff,et al.  TRIM28 Mediates Primer Binding Site-Targeted Silencing of Murine Leukemia Virus in Embryonic Cells , 2007, Cell.

[43]  James H. Thomas,et al.  Coevolution of retroelements and tandem zinc finger genes. , 2011, Genome research.

[44]  D. Haussler,et al.  Gene Isoform Specificity through Enhancer-Associated Antisense Transcription , 2012, PloS one.

[45]  S. Pfaff,et al.  Transposable elements as genetic regulatory substrates in early development. , 2013, Trends in cell biology.