Ongoing evolution of KRAB zinc finger protein-coding genes in modern humans

Background Krüppel-associated box (KRAB) zinc finger proteins (KZFPs) constitute the largest and fastest evolving family of gene regulators encoded by the human genome. Recent data indicate that many KZFPs serve as repressors of transposable element-embedded regulatory sequences (TEeRS) and that the evolutionary turnover of KZFP genes is mainly attributable to the changing transposable element (TE) load of their hosts. However, how natural selection and genetic variation are shaping this process is still poorly defined. Methods Genetic information was collected from nine primate species and 138,500 human genomes. Gene-wide as well as functional amino acid position specific constraint was calculated across all human KZFPs. Results We found that the most conserved KZFPs, some of which go back close to 400 million years, have been subjected to marked negative selection in the evolutionarily recent past and are very homogeneous within the human population. In contrast, younger, largely primate-restricted family members present evidence of less negative selection than the rest of genome and lower levels of coding constraint, particularly within the sequences encoding the functional sites of their zinc finger (ZF) arrays. We defined 33 sets of KZFP paralogs, which pairwise displayed a broad range of coding constraints differentials, with more recently emerged paralogs usually displaying a higher frequency of putatively deleterious mutations and missense variants within the functional sites of their ZF arrays than their source gene. Finally, we identified three KZFP genes more constrained in the genomes of individuals of African ancestry than in Europeans, with their modes of expression or DNA targets pointing to possible links between these inter-populational genetic differences and regional differences in the prevalence of some diseases. Conclusions This work shows how the ongoing selection of KZFPs contributes to modern human genetic variation, in particular through the constraint of putatively deleterious- and missense variants in functional protein sites, and how ongoing interplays between environment and KZFP genes might be impacting the biology of modern humans.

[1]  R. Aebersold,et al.  The interactome of KRAB zinc finger proteins reveals the evolutionary history of their functional diversification , 2019, The EMBO journal.

[2]  R. Jaenisch,et al.  Hominoid-Specific Transposable Elements and KZFPs Facilitate Human Embryonic Genome Activation and Control Transcription in Naive Human ESCs , 2019, Cell stem cell.

[3]  Petra C. Schwalie,et al.  ZFP30 promotes adipogenesis through the KAP1-mediated activation of a retrotransposon-derived Pparg2 enhancer , 2019, Nature Communications.

[4]  Ryan L. Collins,et al.  The mutational constraint spectrum quantified from variation in 141,456 humans , 2020, Nature.

[5]  A. Ferguson-Smith,et al.  ZNF445 is a primary regulator of genomic imprinting , 2019, Genes & development.

[6]  David Haussler,et al.  The UCSC Genome Browser database: 2019 update , 2018, Nucleic Acids Res..

[7]  Silvio C. E. Tosatto,et al.  The Pfam protein families database in 2019 , 2018, Nucleic Acids Res..

[8]  Brent S. Pedersen,et al.  A map of constrained coding regions in the human genome , 2017, bioRxiv.

[9]  Yixuan Wang,et al.  The Role of KRAB-ZFPs in Transposable Element Repression and Mammalian Evolution. , 2017, Trends in genetics : TIG.

[10]  Nicola J. Rinaldi,et al.  Genetic effects on gene expression across human tissues , 2017, Nature.

[11]  Philip M. Kim,et al.  Non-base-contacting residues enable kaleidoscopic evolution of metazoan C2H2 zinc finger DNA binding , 2017, Genome Biology.

[12]  D. Trono,et al.  KRAB zinc finger proteins , 2017, Development.

[13]  M. Daly,et al.  Regional missense constraint improves variant deleteriousness prediction , 2017, bioRxiv.

[14]  Xiaodong Cheng,et al.  A placental growth factor is silenced in mouse embryos by the zinc finger protein ZFP568 , 2017, Science.

[15]  D. Trono,et al.  KRAB zinc-finger proteins contribute to the evolution of gene regulatory networks , 2017, Nature.

[16]  M. Esteller,et al.  Genome-wide DNA methylation pattern in visceral adipose tissue differentiates insulin-resistant from insulin-sensitive obese subjects. , 2016, Translational research : the journal of laboratory and clinical medicine.

[17]  F. Cunningham,et al.  The Ensembl Variant Effect Predictor , 2016, bioRxiv.

[18]  Ayal B. Gussow,et al.  The intolerance to functional genetic variation of protein domains predicts the localization of pathogenic mutations within genes , 2016, Genome Biology.

[19]  A. Wong,et al.  The evolution of gene expression and binding specificity of the largest transcription factor family in primates , 2016, Evolution; international journal of organic evolution.

[20]  Paul J. McLaren,et al.  The Characteristics of Heterozygous Protein Truncating Variants in the Human Genome , 2015, PLoS Comput. Biol..

[21]  James Y. Zou Analysis of protein-coding genetic variation in 60,706 humans , 2015, Nature.

[22]  A. Telenti,et al.  Identification of potential HIV restriction factors by combining evolutionary genomic signatures with functional analyses , 2015, Retrovirology.

[23]  C. Wood,et al.  Epidemiology and Transmission of Kaposi’s Sarcoma-Associated Herpesvirus , 2014, Viruses.

[24]  L. Stubbs,et al.  Deep Vertebrate Roots for Mammalian Zinc Finger Transcription Factor Subfamilies , 2014, Genome biology and evolution.

[25]  A. Levine,et al.  The diversity of zinc-finger genes on human chromosome 19 provides an evolutionary mechanism for defense against inherited endogenous retroviruses , 2013, Cell Death and Differentiation.

[26]  A. Lupo,et al.  KRAB-Zinc Finger Proteins: A Repressor Family Displaying Multiple Biological Functions , 2013, Current genomics.

[27]  Atul J. Butte,et al.  Ethnic Differences in the Relationship Between Insulin Sensitivity and Insulin Response , 2013, Diabetes Care.

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

[29]  Yan Zeng,et al.  Knockdown of ZNF268, which Is Transcriptionally Downregulated by GATA-1, Promotes Proliferation of K562 Cells , 2012, PloS one.

[30]  J. van Helden,et al.  RSAT peak-motifs: motif analysis in full-size ChIP-seq datasets , 2011, Nucleic acids research.

[31]  D. Trono,et al.  In Embryonic Stem Cells, ZFP57/KAP1 Recognize a Methylated Hexanucleotide to Affect Chromatin and DNA Methylation of Imprinting Control Regions , 2011, Molecular cell.

[32]  Cory Y. McLean,et al.  GREAT improves functional interpretation of cis-regulatory regions , 2010, Nature Biotechnology.

[33]  Luke J. Matthews,et al.  The 10kTrees website: A new online resource for primate phylogeny , 2010 .

[34]  H. Wen,et al.  The zinc finger DNA-binding domain of K-RBP plays an important role in regulating Kaposi's sarcoma-associated herpesvirus RTA-mediated gene expression. , 2009, Virology.

[35]  Bartek Wilczynski,et al.  Biopython: freely available Python tools for computational molecular biology and bioinformatics , 2009, Bioinform..

[36]  R. Emerson,et al.  Adaptive Evolution in Zinc Finger Transcription Factors , 2009, PLoS genetics.

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

[38]  Y. Matsui,et al.  Meisetz, A Novel Histone Tri-Methyltransferase, Regulates Meiosis-Specific Epigenesis , 2006, Cell cycle.

[39]  R. Urrutia KRAB-containing zinc-finger repressor proteins , 2003, Genome Biology.

[40]  D. Ballinger,et al.  A Broad Role for the Zinc Finger Protein ZNF202 in Human Lipid Metabolism* , 2000, The Journal of Biological Chemistry.