Young KRAB-zinc finger gene clusters are highly dynamic incubators of ERV-driven genetic heterogeneity in mice

KRAB-zinc finger proteins (KZFPs) comprise the largest family of mammalian transcription factors, rapidly evolving within and between species. Most KZFPs repress endogenous retroviruses (ERVs) and other retrotransposons, with KZFP gene numbers correlating with the ERV load across species, suggesting coevolution. How new KZFPs emerge in response to ERV invasions is currently unknown. Using a combination of long-read sequencing technologies and genome assembly, we present a first detailed comparative analysis of young KZFP gene clusters in the mouse lineage, which has undergone recent KZFP gene expansion and ERV infiltration. Detailed annotation of KZFP genes in a cluster on Mus musculus Chromosome 4 revealed parallel expansion and diversification of this locus in different mouse strains (C57BL/6J, 129S1/SvImJ and CAST/EiJ) and species (Mus spretus and Mus pahari). Our data supports a model by which new ERV integrations within young KZFP gene clusters likely promoted recombination events leading to the emergence of new KZFPs that repress them. At the same time, ERVs also increased their numbers by duplication instead of retrotransposition alone, unraveling a new mechanism for ERV enrichment at these loci. ChIP-seq experiments for H3K3me3 and H3K9me3 in pure C57BL/6J mESCs and ChIP-seq experiments for individual mouse KRAB-zinc finger proteins overexpressed in F9 embryo carcinoma cells.