Lineage-specific expansion of KRAB zinc-finger transcription factor genes: implications for the evolution of vertebrate regulatory networks

On the dependent recognition of some long zinc finger proteins

The human genome contains about 800 C2H2 zinc finger proteins (ZFPs), and most of them are composed of long arrays of zinc fingers. Standard ZFP recognition model asserts longer finger arrays should recognize longer DNA-binding sites. However, recent experimental efforts to identify in vivo ZFP binding sites contradict this assumption, with many exhibiting short motifs. Here we use ZFY, CTCF, ZIM3, and ZNF343 as examples to address three closely related questions: What are the reasons that impede current motif discovery methods? What are the functions of those seemingly unused fingers and how can we improve the motif discovery algorithms based on long ZFPs' biophysical properties? Using ZFY, we employed a variety of methods and find evidence for 'dependent recognition' where downstream fingers can recognize some previously undiscovered motifs only in the presence of an intact core site. For CTCF, high-throughput measurements revealed its upstream specificity profile depends on the strength of its core. Moreover, the binding strength of the upstream site modulates CTCF's sensitivity to different epigenetic modifications within the core, providing new insight into how the previously identified intellectual disability-causing and cancer-related mutant R567W disrupts upstream recognition and deregulates the epigenetic control by CTCF. Our results establish that, because of irregular motif structures, variable spacing and dependent recognition between sub-motifs, the specificities of long ZFPs are significantly underestimated, so we developed an algorithm, ModeMap, to infer the motifs and recognition models of ZIM3 and ZNF343, which facilitates high-confidence identification of specific binding sites, including repeats-derived elements. With revised concept, technique, and algorithm, we can discover the overlooked specificities and functions of those 'extra' fingers, and therefore decipher their broader roles in human biology and diseases.

Structural Evaluation and Conformational Dynamics of ZNF141T474I Mutation Provoking Postaxial Polydactyly Type A

Postaxial Polydactyly (PAP) is a congenital disorder of limb abnormalities characterized by posterior extra digits. Mutations in the N-terminal region of the Zinc finger protein 141 (ZNF141) gene were recently linked with PAP type A. Zinc finger proteins exhibit similarity at their N-terminal regions due to C2-H2 type Zinc finger domains, but their functional preferences vary significantly by the binding patterns of DNA. Methods: This study delineates the pathogenic association, miss-fold aggregation, and conformational paradigm of a missense variant (c.1420C > T; p.T474I) in ZNF141 gene segregating PAP through a molecular dynamics simulations approach. Results: In ZNF141 protein, helices play a crucial role by attaching three specific target DNA base pairs. In ZNF141T474I protein, H1, H3, and H6 helices attain more flexibility by acquiring loop conformation. The outward disposition of the proximal portion of H9-helix in mutant protein occurs due to the loss of prior beta-hairpins at the C terminal region of the C2-H2 domain. The loss of hydrogen bonds and exposure of hydrophobic residues to solvent and helices turning to loops cause dysfunction of ZNF141 protein. These significant changes in the stability and conformation of the mutant protein were validated using essential dynamics and cross-correlation maps, which revealed that upon point mutation, the overall motion of the proteins and the correlation between them were completely different, resulting in Postaxial polydactyly type A. Conclusions: This study provides molecular insights into the structural association of ZNF141 protein with PAP type A. Identification of active site residues and legends offers new therapeutic targets for ZNF141 protein. Further, it reiterates the functional importance of the last residue of a protein.

KRAB family is involved in network shifts in response to osmotic stress in camels

A feature of the camel is its tolerance to osmotic stress. However, few studies of osmotic stress in vivo or comparative analyses between different tissues of the camel have been performed. Here, we report the roles of Krüppel-associated box domain containing zinc-finger repressor proteins (KRAB-ZFPs) in transcriptional networks under osmotic stress in camels by analyzing transcriptomes of four different tissues under various osmotic conditions. We found that 273 of 278 KRAB-ZFPs were expressed in our data set, being involved in all of the 65 identified networks and exhibiting their extensive functional diversity. We also found that 110 KRAB-ZFPs were hub genes involved in more than half of the networks. We demonstrated that the osmotic stress response is involved in network shifts and that KRAB-ZFPs mediate this process. Finally, we presented the diverse mechanisms of osmotic stress responses in different tissues. These results revealed the genetic architecture of systematic physiological response in vivo to osmotic stress in camels. Our work will lead to new directions for studying the mechanism of osmotic stress response in anti-arid mammals.

Genetic Variation in Human Gene Regulatory Factors Uncovers Regulatory Roles in Local Adaptation and Disease

Differences in gene regulation have been suggested to play essential roles in the evolution of phenotypic changes. Although DNA changes in cis-regulatory elements affect only the regulation of its corresponding gene, variations in gene regulatory factors (trans) can have a broader effect, because the expression of many target genes might be affected. Aiming to better understand how natural selection may have shaped the diversity of gene regulatory factors in human, we assembled a catalog of all proteins involved in controlling gene expression. We found that at least five DNA-binding transcription factor classes are enriched among genes located in candidate regions for selection, suggesting that they might be relevant for understanding regulatory mechanisms involved in human local adaptation. The class of KRAB-ZNFs, zinc-finger (ZNF) genes with a Krüppel-associated box, stands out by first, having the most genes located on candidate regions for positive selection. Second, displaying most nonsynonymous single nucleotide polymorphisms (SNPs) with high genetic differentiation between populations within these regions. Third, having 27 KRAB-ZNF gene clusters with high extended haplotype homozygosity. Our further characterization of nonsynonymous SNPs in ZNF genes located within candidate regions for selection, suggests regulatory modifications that might influence the expression of target genes at population level. Our detailed investigation of three candidate regions revealed possible explanations for how SNPs may influence the prevalence of schizophrenia, eye development, and fertility in humans, among other phenotypes. The genetic variation we characterized here may be responsible for subtle to rough regulatory changes that could be important for understanding human adaptation.

The mouse genome displays highly dynamic populations of KRAB-zinc finger protein genes and related genetic units

KRAB-containing poly-zinc finger proteins (KZFPs) constitute the largest family of transcription factors encoded by mammalian genomes, and growing evidence indicates that they fulfill functions critical to both embryonic development and maintenance of adult homeostasis. KZFP genes underwent broad and independent waves of expansion in many higher vertebrates lineages, yet comprehensive studies of members harbored by a given species are scarce. Here we present a thorough analysis of KZFP genes and related units in the murine genome. We first identified about twice as many elements than previously annotated as either KZFP genes or pseudogenes, notably by assigning to this family an entity formerly considered as a large group of Satellite repeats. We then could delineate an organization in clusters distributed throughout the genome, with signs of recombination, translocation, duplication and seeding of new sites by retrotransposition of KZFP genes and related genetic units (KZFP/rGUs). Moreover, we harvested evidence indicating that closely related paralogs had evolved through both drifting and shifting of sequences encoding for zinc finger arrays. Finally, we could demonstrate that the KAP1-SETDB1 repressor complex tames the expression of KZFP/rGUs within clusters, yet that the primary targets of this regulation are not the KZFP/rGUs themselves but enhancers contained in neighboring endogenous retroelements and that, underneath, KZFPs conserve highly individualized patterns of expression.

Genome-wide Regulatory Roles of the C2H2-type Zinc Finger Protein ZNF764 on the Glucocorticoid Receptor

The C2H2-type zinc finger protein ZNF764 acts as an enhancer for several steroid hormone receptors, and haploinsufficiency of this gene may be responsible for tissue resistance to multiple steroid hormones including glucocorticoids observed in a patient with 16p11.2 microdeletion. We examined genome-wide regulatory actions of ZNF764 on the glucocorticoid receptor (GR) in HeLa cells as a model system. ZNF764- and GR-binding sites demonstrated similar distribution in various genomic features. They positioned predominantly around 50-500 kbs from the transcription start sites of their nearby genes, and were closely localized with each other, overlapping in ~37% of them. ZNF764 demonstrated differential on/off effects on GR-binding and subsequent mRNA expression: some genes were highly dependent on the presence/absence of ZNF764, but others were not. Pathway analysis revealed that these 3 gene groups were involved in distinct cellular activities. ZNF764 physically interacted with GR at ligand-binding domain through its KRAB domain, and both its physical interaction to GR and zinc finger domain appear to be required for ZNF764 to regulate GR transcriptional activity. Thus, ZNF764 is a cofactor directing GR transcriptional activity toward specific biologic pathways by changing GR binding and transcriptional activity on the glucocorticoid-responsive genes.

Transcription factor ZNF25 is associated with osteoblast differentiation of human skeletal stem cells

Background

The differentiation of human bone marrow derived skeletal stem cells (known as human bone marrow stromal or mesenchymal stem cells, hMSCs) into osteoblasts involves the activation of a small number of well-described transcription factors. To identify additional osteoblastic transcription factors, we studied gene expression of hMSCs during ex vivo osteoblast differentiation.

Results

Clustering of gene expression, and literature investigation, revealed three transcription factors of interest – ZNF25, ZNF608 and ZBTB38. siRNA knockdown of ZNF25 resulted in significant suppression of alkaline phosphatase (ALP) activity. This effect was not present for ZNF608 and ZBTB38. To identify possible target genes of ZNF25, we analyzed gene expression following ZNF25 siRNA knockdown. This revealed a 23-fold upregulation of matrix metallopeptidase 1 and an 18-fold upregulation of leucine-rich repeat containing G protein-coupled receptor 5 and RAN-binding protein 3-like. We also observed enrichment in extracellular matrix organization, skeletal system development and regulation of ossification in the entire upregulated set of genes. Consistent with its function as a transcription factor during osteoblast differentiation of hMSC, we showed that the ZNF25 protein exhibits nuclear localization and is expressed in osteoblastic and osteocytic cells in vivo. ZNF25 is conserved in tetrapod vertebrates and contains a KRAB (Krueppel-associated box) transcriptional repressor domain.

Conclusions

This study shows that the uncharacterized transcription factor, ZNF25, is associated with differentiation of hMSC to osteoblasts.

The stem cell zinc finger 1 (SZF1)/ZNF589 protein has a human-specific evolutionary nucleotide DNA change and acts as a regulator of cell viability in the hematopoietic system

The stem cell zinc finger 1 (SZF1)/ZNF589 protein belongs to the large family of Krüppel-associated box domain-zinc finger (KRAB-ZNF) transcription factors, which are present only in higher vertebrates and epigenetically repress transcription by recruiting chromatin-modifying complexes to the promoter regions of their respective target genes. Although the distinct biological functions of most KRAB-ZNF proteins remain unknown, recent publications indicate their implication in fundamental processes, such as cell proliferation, apoptosis, differentiation, development, and tumorigenesis. SZF1/ZNF589 was first identified as a gene with SZF1-1 isoform specifically expressed in CD34(+) hematopoietic cells, strongly suggesting a role in epigenetic control of gene expression in hematopoietic stem/progenitor cells (HSPCs). However, the function of SZF1/ZNF589 in hematopoiesis has not yet been elucidated. Our study reveals SZF1/ZNF589 as a gene with a human-specific nucleotide DNA-change, conferring potential species-specific functional properties. Through shRNA-mediated loss-of-function experiments, we found that changes in expression of fundamental apoptosis-controlling genes are induced on SZF1/ZNF589 knockdown, resulting in inhibited growth of hematopoietic cell lines and decreased progenitor potential of primary human bone marrow CD34(+) cells. Moreover, we found that the SZF1/ZNF589 gene is differentially regulated during hypoxia in CD34(+) HSPCs in a cytokine-dependent manner, implicating its possible involvement in the maintenance of the hypoxic physiologic status of hematopoietic stem cells. Our results establish the role of SZF1/ZNF589 as a new functional regulator of the hematopoietic system.

Epigenetic regulation of puberty via Zinc finger protein-mediated transcriptional repression

In primates, puberty is unleashed by increased GnRH release from the hypothalamus following an interval of juvenile quiescence. GWAS implicates Zinc finger (ZNF) genes in timing human puberty. Here we show that hypothalamic expression of several ZNFs decreased in agonadal male monkeys in association with the pubertal reactivation of gonadotropin secretion. Expression of two of these ZNFs, GATAD1 and ZNF573, also decreases in peripubertal female monkeys. However, only GATAD1 abundance increases when gonadotropin secretion is suppressed during late infancy. Targeted delivery of GATAD1 or ZNF573 to the rat hypothalamus delays puberty by impairing the transition of a transcriptional network from an immature repressive epigenetic configuration to one of activation. GATAD1 represses transcription of two key puberty-related genes, KISS1 and TAC3, directly, and reduces the activating histone mark H3K4me2 at each promoter via recruitment of histone demethylase KDM1A. We conclude that GATAD1 epitomizes a subset of ZNFs involved in epigenetic repression of primate puberty.

Pervasive variation of transcription factor orthologs contributes to regulatory network evolution

Differences in transcriptional regulatory networks underlie much of the phenotypic variation observed across organisms. Changes to cis-regulatory elements are widely believed to be the predominant means by which regulatory networks evolve, yet examples of regulatory network divergence due to transcription factor (TF) variation have also been observed. To systematically ascertain the extent to which TFs contribute to regulatory divergence, we analyzed the evolution of the largest class of metazoan TFs, Cys2-His2 zinc finger (C2H2-ZF) TFs, across 12 Drosophila species spanning ~45 million years of evolution. Remarkably, we uncovered that a significant fraction of all C2H2-ZF 1-to-1 orthologs in flies exhibit variations that can affect their DNA-binding specificities. In addition to loss and recruitment of C2H2-ZF domains, we found diverging DNA-contacting residues in ~44% of domains shared between D. melanogaster and the other fly species. These diverging DNA-contacting residues, found in ~70% of the D. melanogaster C2H2-ZF genes in our analysis and corresponding to ~26% of all annotated D. melanogaster TFs, show evidence of functional constraint: they tend to be conserved across phylogenetic clades and evolve slower than other diverging residues. These same variations were rarely found as polymorphisms within a population of D. melanogaster flies, indicating their rapid fixation. The predicted specificities of these dynamic domains gradually change across phylogenetic distances, suggesting stepwise evolutionary trajectories for TF divergence. Further, whereas proteins with conserved C2H2-ZF domains are enriched in developmental functions, those with varying domains exhibit no functional enrichments. Our work suggests that a subset of highly dynamic and largely unstudied TFs are a likely source of regulatory variation in Drosophila and other metazoans.

Deep vertebrate roots for mammalian zinc finger transcription factor subfamilies

While many vertebrate transcription factor (TF) families are conserved, the C2H2 zinc finger (ZNF) family stands out as a notable exception. In particular, novel ZNF gene types have arisen, duplicated, and diverged independently throughout evolution to yield many lineage-specific TF genes. This evolutionary dynamic not only raises many intriguing questions but also severely complicates identification of those ZNF genes that remain functionally conserved. To address this problem, we searched for vertebrate “DNA binding orthologs” by mining ZNF loci from eight sequenced genomes and then aligning the patterns of DNA-binding amino acids, or “fingerprints,” extracted from the encoded ZNF motifs. Using this approach, we found hundreds of lineage-specific genes in each species and also hundreds of orthologous groups. Most groups of orthologs displayed some degree of fingerprint divergence between species, but 174 groups showed fingerprint patterns that have been very rigidly conserved. Focusing on the dynamic KRAB-ZNF subfamily—including nearly 400 human genes thought to possess potent KRAB-mediated epigenetic silencing activities—we found only three genes conserved between mammals and nonmammalian groups. These three genes, members of an ancient familial cluster, encode an unusual KRAB domain that functions as a transcriptional activator. Evolutionary analysis confirms the ancient provenance of this activating KRAB and reveals the independent expansion of KRAB-ZNFs in every vertebrate lineage. Most human ZNF genes, from the most deeply conserved to the primate-specific genes, are highly expressed in immune and reproductive tissues, indicating that they have been enlisted to regulate evolutionarily divergent biological traits.

Recruitment and biological consequences of histone modification of H3K27me3 and H3K9me3

Two histone marks, H3K27me3 and H3K9me3, are well known for their repressive roles in the genic and nongenic regions of metazoan genomes. Several protein complexes are known to be responsible for generating these marks, including polycomb repression complex 2 and several H3K9 methylases. Recent studies have shown that the targeting of these histone-modifying complexes within mammalian genomes may be mediated through several DNA-binding proteins, including AEBP2, JARID2, and YY1. In this review, we discuss the potential targeting mechanisms in light of the recent results that have been derived from genome-wide chromatin immunoprecipitation sequencing data and the in vivo functions of these two histone marks in light of the results derived from mouse and human genetic studies.

The KRAB zinc finger protein RSL1 modulates sex-biased gene expression in liver and adipose tissue to maintain metabolic homeostasis

Krüppel-associated box zinc finger proteins (KRAB-ZFPs) are a huge family of vertebrate-specific repressors that modify gene expression in an epigenetic manner. Despite a well-defined repression mechanism, few biological roles or gene targets of KRAB-ZFP are known. Regulator of sex-limitation 1 (RSL1) is a mouse KRAB-ZFP that enforces male-predominant expression in the liver, affecting body mass and pubertal timing. Here we show that female but not male Rsl1(-/-) mice gain more weight than wild-type mice on a high-fat diet (HFD) and that key liver and white adipose tissue (WAT) metabolic genes are altered in both Rsl1(-/-) sexes in response to dietary stress. Expression profiling of Rsl1-sensitive genes in liver and WAT indicates that RSL1 accentuates sex-biased gene expression in liver but greatly diminishes it in WAT. RSL1 expression solely in liver is sufficient to limit diet-induced weight gain and suppress lipogenic genes in WAT, indicating that RSL1 balances metabolism via liver-to-adipose-tissue communication. RSL1's effects on adult physiology exemplify a significant modulatory capacity of KRAB-ZFPs, in the absence of which there is widespread metabolic dysregulation. This ability to buffer against gene expression noise, coupled with extensive individual genetic variation, highlights the enormous potential of KRAB-Zfp genes as candidate risk factors for complex diseases.

A study on the distribution of 37 well conserved families of C2H2 zinc finger genes in eukaryotes

Background

The C2H2 zinc-finger (ZNF) containing gene family is one of the largest and most complex gene families in metazoan genomes. These genes are known to exist in almost all eukaryotes, and they constitute a major subset of eukaryotic transcription factors. The genes of this family usually occur as clusters in genomes and are thought to have undergone a massive expansion in vertebrates by multiple tandem duplication events (BMC Evol Biol 8:176, 2008).

Results

In this study, we combined two popular approaches for homolog detection, Reciprocal Best Hit (RBH) (Proc Natl Acad Sci USA 95:6239–6244, 1998) and Hidden–Markov model (HMM) profiles search (Bioinformatics 14:755-763, 1998), on a diverse set of complete genomes of 124 eukaryotic species ranging from excavates to humans to identify all detectable members of 37 C2H2 ZNF gene families. We succeeded in identifying 3,890 genes as distinct members of 37 C2H2 gene families. These 37 families are distributed among the eukaryotes as progressive additions of gene blocks with increasing complexity of the organisms. The first block featuring the protists had 7 families, the second block featuring plants had 2 families, the third block featuring the fungi had 2 families (one of which was also present in plants) and the final block consisted of metazoans with 25 families. Among the metazoans, the simpler unicellular metazoans had just 15 of the 25 families while most of the bilaterians had all 25 families making up a total of 37 families. Multiple potential examples of lineage-specific gene duplications and gene losses were also observed.

Conclusions

Our hybrid approach combines features of the both RBH and HMM methods for homolog detection. This largely automated technique is much faster than manual methods and is able to detect homologs accurately and efficiently among a diverse set of organisms. Our analysis of the 37 evolutionarily conserved C2H2 ZNF gene families revealed a stepwise appearance of ZNF families, agreeing well with the phylogenetic relationship of the organisms compared and their presumed stepwise increase in complexity (Science 300:1694, 2003).

Global analysis of target genes of 21 members of the ZAD transcription factor family in Drosophila melanogaster

The zinc-finger associated domain (ZAD) family is the largest transcription factor family in dipteran insects. Still, their functional significance is barely recognized in the literature due in part to their resistance to mutagenesis screens in genetic studies. Therefore, we employed in vitro techniques to identify the DNA-binding characteristics of several members of the Drosophila melanogaster ZAD family in an effort to study their target genes. In this comprehensive investigation, we constructed a panel of GST-Zinc finger (ZnF) array chimera from 21 selected ZAD proteins and used them to select binding sites from an oligonucleotide library by employing electrophoretic mobility shift assays (EMSA). Samples of the binding population were sequenced and used to derive DNA-binding consensus sequence for each member. These consensus sequences were tested for complex formation with their respective protein chimera and the specificity of binding ascertained by competition EMSA. Bioinformatics tools were used to identify potential genetic targets. The identified consensus sequences were distinct for each member and the putative genomic targets were clustered in the regulatory regions of specific genes. This appears to be consistent with a conservation of function between members and also suggests that the overlapping functions of ZAD proteins are the result of positive selection to maintain redundancy and not simply artifacts of recent expansion. Putative target genes suggest a major role of the ZAD family members in the regulation of several early developmental genes including homeobox transcription factors.

Coevolution of retroelements and tandem zinc finger genes

Vertebrate genomes encode large and highly variable numbers of tandem C2H2 zinc finger (tandem ZF) transcription factor proteins. In mammals, most tandem ZF genes also encode a KRAB domain (KZNF proteins). Very little is known about what forces have driven the number and diversity of tandem ZF genes. Recent studies suggest that one role of KZNF proteins is to bind and repress transcription of exogenous retroviruses and their endogenous counterpart LTR retroelements. We report a striking correlation across vertebrate genomes between the number of LTR retroelements and the number of host tandem ZF genes. This correlation is specific to LTR retroelements and ZF genes and was not explained by covariation in other genomic features. We further show that recently active LTR retroelements are correlated with recent tandem ZF gene duplicates across vertebrates. On branches of the primate phylogeny, we find that the appearance of new families of endogenous retroviruses is strongly predictive of the appearance of new duplicate KZNF genes. We hypothesize that retroviral and LTR retroelement burden drives evolution of host tandem ZF genes. This hypothesis is consistent with previously described molecular evolutionary patterns in duplicate ZF genes throughout vertebrates. To further explore these patterns, we investigated 34 duplicate human KZNF gene pairs, all of which underwent an early burst of divergence in the major nucleotide contact residues of their ZF domains, followed by purifying selection in both duplicates. Our results support a host-pathogen model for tandem ZF gene evolution, in which new LTR retroelement challenges drive duplication and divergence of host tandem ZF genes.

Gain, loss and divergence in primate zinc-finger genes: a rich resource for evolution of gene regulatory differences between species

The molecular changes underlying major phenotypic differences between humans and other primates are not well understood, but alterations in gene regulation are likely to play a major role. Here we performed a thorough evolutionary analysis of the largest family of primate transcription factors, the Krüppel-type zinc finger (KZNF) gene family. We identified and curated gene and pseudogene models for KZNFs in three primate species, chimpanzee, orangutan and rhesus macaque, to allow for a comparison with the curated set of human KZNFs. We show that the recent evolutionary history of primate KZNFs has been complex, including many lineage-specific duplications and deletions. We found 213 species-specific KZNFs, among them 7 human-specific and 23 chimpanzee-specific genes. Two human-specific genes were validated experimentally. Ten genes have been lost in humans and 13 in chimpanzees, either through deletion or pseudogenization. We also identified 30 KZNF orthologs with human-specific and 42 with chimpanzee-specific sequence changes that are predicted to affect DNA binding properties of the proteins. Eleven of these genes show signatures of accelerated evolution, suggesting positive selection between humans and chimpanzees. During primate evolution the most extensive re-shaping of the KZNF repertoire, including most gene additions, pseudogenizations, and structural changes occurred within the subfamily homininae. Using zinc finger (ZNF) binding predictions, we suggest potential impact these changes have had on human gene regulatory networks. The large species differences in this family of TFs stands in stark contrast to the overall high conservation of primate genomes and potentially represents a potent driver of primate evolution.

Genome structure analysis of molluscs revealed whole genome duplication and lineage specific repeat variation

Comparative genome structure analysis allows us to identify novel genes, repetitive sequences and gene duplications. To explore lineage-specific genomic changes of the molluscs that is good model for development of nervous system in invertebrate, we conducted comparative genome structure analyses of three molluscs, pygmy squid, nautilus and scallops using partial genome shotgun sequencing. Most effective elements on the genome structural changes are repetitive elements (REs) causing expansion of genome size and whole genome duplication producing large amount of novel functional genes. Therefore, we investigated variation and proportion of REs and whole genome duplication. We, first, identified variations of REs in the three molluscan genomes by homology-based and de novo RE detection. Proportion of REs were 9.2%, 4.0%, and 3.8% in the pygmy squid, nautilus and scallop, respectively. We, then, estimated genome size of the species as 2.1, 4.2 and 1.8 Gb, respectively, with 2× coverage frequency and DNA sequencing theory. We also performed a gene duplication assay based on coding genes, and found that large-scale duplication events occurred after divergence from the limpet Lottia, an out-group of the three molluscan species. Comparison of all the results suggested that RE expansion did not relate to the increase in genome size of nautilus. Despite close relationships to nautilus, the squid has the largest portion of REs and smaller genome size than nautilus. We also identified lineage-specific RE and gene-family expansions, possibly relate to acquisition of the most complicated eye and brain systems in the three species.

Function and Evolution of C2H2 Zinc Finger Arrays

Krüppel-type or C2H2 zinc fingers represent a dominant DNA-binding motif in eukaryotic transcription factor (TF) proteins. In Krüppel-type (KZNF) TFs, KZNF motifs are arranged in arrays of three to as many as 40 tandem units, which cooperate to define the unique DNA recognition properties of the protein. Each finger contains four amino acids located at specific positions, which are brought into direct contact with adjacent nucleotides in the DNA sequence as the KZNF array winds around the major groove of the alpha helix. This arrangement creates an intimate and potentially predictable relationship between the amino acid sequence of KZNF arrays and the nucleotide sequence of target binding sites. The large number of possible combinations and arrangements of modular KZNF motifs, and the increasing lengths of KZNF arrays in vertebrate species, has created huge repertoires of functionally unique TF proteins. The properties of this versatile DNA-binding motif have been exploited independently many times over the course of evolution, through attachment to effector motifs that confer activating, repressing or other activities to the proteins. Once created, some of these novel inventions have expanded in specific evolutionary clades, creating large families of TFs that are lineage- or species-unique. This chapter reviews the properties and their remarkable evolutionary history of eukaryotic KZNF TF proteins, with special focus on large families that dominate the TF landscapes in different metazoan species.

The ancient mammalian KRAB zinc finger gene cluster on human chromosome 8q24.3 illustrates principles of C2H2 zinc finger evolution associated with unique expression profiles in human tissues

Background

Expansion of multi-C2H2 domain zinc finger (ZNF) genes, including the Krüppel-associated box (KRAB) subfamily, paralleled the evolution of tetrapodes, particularly in mammalian lineages. Advances in their cataloging and characterization suggest that the functions of the KRAB-ZNF gene family contributed to mammalian speciation.

Results

Here, we characterized the human 8q24.3 ZNF cluster on the genomic, the phylogenetic, the structural and the transcriptome level. Six (ZNF7, ZNF34, ZNF250, ZNF251, ZNF252, ZNF517) of the seven locus members contain exons encoding KRAB domains, one (ZNF16) does not. They form a paralog group in which the encoded KRAB and ZNF protein domains generally share more similarities with each other than with other members of the human ZNF superfamily. The closest relatives with respect to their DNA-binding domain were ZNF7 and ZNF251. The analysis of orthologs in therian mammalian species revealed strong conservation and purifying selection of the KRAB-A and zinc finger domains. These findings underscore structural/functional constraints during evolution. Gene losses in the murine lineage (ZNF16, ZNF34, ZNF252, ZNF517) and potential protein truncations in primates (ZNF252) illustrate ongoing speciation processes. Tissue expression profiling by quantitative real-time PCR showed similar but distinct patterns for all tested ZNF genes with the most prominent expression in fetal brain. Based on accompanying expression signatures in twenty-six other human tissues ZNF34 and ZNF250 revealed the closest expression profiles. Together, the 8q24.3 ZNF genes can be assigned to a cerebellum, a testis or a prostate/thyroid subgroup. These results are consistent with potential functions of the ZNF genes in morphogenesis and differentiation. Promoter regions of the seven 8q24.3 ZNF genes display common characteristics like missing TATA-box, CpG island-association and transcription factor binding site (TFBS) modules. Common TFBS modules partly explain the observed expression pattern similarities.

Conclusions

The ZNF genes at human 8q24.3 form a relatively old mammalian paralog group conserved in eutherian mammals for at least 130 million years. The members persisted after initial duplications by undergoing subfunctionalizations in their expression patterns and target site recognition. KRAB-ZNF mediated repression of transcription might have shaped organogenesis in mammalian ontogeny.

Evolution of C2H2-zinc finger genes and subfamilies in mammals: species-specific duplication and loss of clusters, genes and effector domains

Background

C2H2 zinc finger genes (C2H2-ZNF) constitute the largest class of transcription factors in humans and one of the largest gene families in mammals. Often arranged in clusters in the genome, these genes are thought to have undergone a massive expansion in vertebrates, primarily by tandem duplication. However, this view is based on limited datasets restricted to a single chromosome or a specific subset of genes belonging to the large KRAB domain-containing C2H2-ZNF subfamily.

Results

Here, we present the first comprehensive study of the evolution of the C2H2-ZNF family in mammals. We assembled the complete repertoire of human C2H2-ZNF genes (718 in total), about 70% of which are organized into 81 clusters across all chromosomes. Based on an analysis of their N-terminal effector domains, we identified two new C2H2-ZNF subfamilies encoding genes with a SET or a HOMEO domain. We searched for the syntenic counterparts of the human clusters in other mammals for which complete gene data are available: chimpanzee, mouse, rat and dog. Cross-species comparisons show a large variation in the numbers of C2H2-ZNF genes within homologous mammalian clusters, suggesting differential patterns of evolution. Phylogenetic analysis of selected clusters reveals that the disparity in C2H2-ZNF gene repertoires across mammals not only originates from differential gene duplication but also from gene loss. Further, we discovered variations among orthologs in the number of zinc finger motifs and association of the effector domains, the latter often undergoing sequence degeneration. Combined with phylogenetic studies, physical maps and an analysis of the exon-intron organization of genes from the SCAN and KRAB domains-containing subfamilies, this result suggests that the SCAN subfamily emerged first, followed by the SCAN-KRAB and finally by the KRAB subfamily.

Conclusion

Our results are in agreement with the "birth and death hypothesis" for the evolution of C2H2-ZNF genes, but also show that this hypothesis alone cannot explain the considerable evolutionary variation within the subfamilies of these genes in mammals. We, therefore, propose a new model involving the interdependent evolution of C2H2-ZNF gene subfamilies.

The Ajuba LIM domain protein is a corepressor for SNAG domain mediated repression and participates in nucleocytoplasmic Shuttling

The SNAG repression domain is comprised of a highly conserved 21-amino acid sequence, is named for its presence in the Snail/growth factor independence-1 class of zinc finger transcription factors, and is present in a variety of proto-oncogenic transcription factors and developmental regulators. The prototype SNAG domain containing oncogene, growth factor independence-1, is responsible for the development of T cell thymomas. The SNAIL proteins also encode the SNAG domain and play key roles in epithelial mesenchymal differentiation events during development and metastasis. Significantly, these oncogenic functions require a functional SNAG domain. The molecular mechanisms of SNAG domain-mediated transcriptional repression are largely unknown. Using a yeast two-hybrid strategy, we identified Ajuba, a multiple LIM domain protein that can function as a corepressor for the SNAG domain. Ajuba interacts with the SNAG domain in vitro and in vivo, colocalizes with it, and enhances SNAG-mediated transcriptional repression. Ajuba shuttles between the cytoplasm and the nucleus and may form a novel intracellular signaling system. Using an integrated reporter gene combined with chromatin immunoprecipitation, we observed rapid, SNAG-dependent assembly of a multiprotein complex that included Ajuba, SNAG, and histone modifications consistent with the repressed state. Thus, SNAG domain proteins may bind Ajuba, trapping it in the nucleus where it functions as an adapter or molecular scaffold for the assembly of macromolecular repression complexes at target promoters.

Genomic organization and imprinting of the Peg3 domain in bovine

Using multiple mammalian genomic sequences, we have analyzed the evolution and imprinting of several genes located in the Peg3 domain, including Mim1 (approved name, Mimt1), Usp29, Zim3, and Zfp264. A series of comparative analyses shows that the overall genomic structure of this 500-kb imprinted domain has been well maintained throughout mammalian evolution but that several lineage-specific changes have also occurred in each species. In the bovine domain, Usp29 has lost its protein-coding capability, Zim3 has been duplicated, and the expression of Zfp264 has become biallelic in brain and testis, which differs from paternal expression of mouse Zfp264 in brain. In contrast, the two transcript genes of cow, Mim1 and Usp29, both lacking protein-coding capability, are still expressed mainly from the paternal allele, indicating the imprinting of these two genes in cow. The imprinting of Mim1 and Usp29 along with Peg3 is the most evolutionarily selected feature in this imprinted domain, suggesting significant function of these three genes, either as protein-coding or as untranslated transcript genes.

A novel KRAB-Zinc finger protein interacts with latency-associated nuclear antigen of Kaposi's sarcoma-associated herpesvirus and activates transcription via terminal repeat sequences

Kaposi's sarcoma-associated herpesvirus (KSHV) establishes latent infection in various cells in vitro as well as KSHV-associated tumor cells in vivo. The latency-associated nuclear antigen (LANA) of KSHV is one of a small number of genes expressed in the latent phase of KSHV infection. This antigen is crucial for establishment of the latent infection, such as replication of KSHV genomic DNA and maintenance of infection via direct interaction with terminal repeats (TRs) in the viral genome. Using a yeast two-hybrid screening method, we isolated a novel LANA-interacting protein (designated as KZLP; KRAB Zinc finger LANA interacting Protein) from a human peripheral leukocyte cDNA library. KZLP encodes a KRAB domain and 12 Kruppel-type zinc fingers. Reverse transcription polymerase chain reaction showed that KZLP was expressed ubiquitously in various cell lines including those infected with KSHV. A luciferase assay showed that KZLP could activate the KSHV open reading frame K1 promoter containing TRs in 293T cells, and that such activation required multiple TR sequences. In contrast, LANA repressed the activity of the K1 promoter through TRs, and again this repression required multiple TR units. Moreover, LANA almost completely abrogated the KZLP-mediated transcriptional activation. Our results suggest that KZLP and LANA regulate gene expression through TRs in the KSHV viral genome, including the K1 gene in latent KSHV-infected cells.

Expansion of symmetric exon-bordering domains does not explain evolution of lineage specific genes in mammals

In order to examine the evolution of lineage specific genes, we analyzed intron phase distributions and exon-bordering domains in primate and rodent specific genes. We found that the expansion of symmetric exon-bordering domains could not explain the evolution of lineage specific genes. Rather internal intron loss of a domain can partially explain the excess of class 1-1 intron phases in the lineage specific genes. We suggest the event that led to excess of symmetric exons in lineage specific genes had little bearing on shaping the phenotypes specific to the individual lineage. Instead, Kruppel-associated box (KRAB) proteins associated with zinc finger C2H2 (zf-C2H2) type are likely to be responsible for the lineage specific function.

Survey of microarray technologies suitable to elucidate transcriptional networks as exemplified by studying KRAB zinc finger gene families

Current microarray systems are suitable to monitor genome-wide expression patterns, to detect single-nucleotide polymorphisms (SNP), to identify target genes of transcription factors and DNA-protein interaction sites thereof as well as to determine genomic sites that are modified by methylation of CpG islands. In this review, advantages and limitations of individual microarray technologies are presented as well as experiences from ongoing studies on KRAB zinc finger gene families are taken to exemplify how different microarray approaches are applicable to elucidate complex transcriptional networks of gene regulation. However, bioinformaticians should be aware that each microarray technology has limitations in its sensitivity and selectivity that has to be taken into account once data mining on comprehensive genome-wide microarray data is conducted. In many cases, microarray results are the initial step to identify target genes of interest and to study the molecular regulation of biological processes thereof followed and validated by complementary proteome, metabolome or toponome analysis. Thus, microarray technologies can be considered a reliable approach for determining gene functions that might be modulated by electromagnetic fields.

Comparative genomics modeling of the NRSF/REST repressor network: from single conserved sites to genome-wide repertoire

We constructed and applied an open source informatic framework called Cistematic in an effort to predict the target gene repertoire for transcription factors with large binding sites. Cistematic uses two different evolutionary conservation-filtering algorithms in conjunction with several analysis modules. Beginning with a single conserved and biologically tested site for the neuronal repressor NRSF/REST, Cistematic generated a refined PSFM (position specific frequency matrix) based on conserved site occurrences in mouse, human, and dog genomes. Predictions from this model were validated by chromatin immunoprecipitation (ChIP) followed by quantitative PCR. The combination of transfection assays and ChIP enrichment data provided an objective basis for setting a threshold for membership and rank-ordering a final gene cohort model consisting of 842 high-confidence sites in the human genome associated with 733 genes. Statistically significant enrichment of NRSE-associated genes was found for neuron-specific Gene Ontology (GO) terms and neuronal mRNA expression profiles. A more extensive evolutionary survey showed that NRSE sites matching the PSFM model exist in roughly similar numbers in all fully sequenced vertebrate genomes but are notably absent from invertebrate and protochordate genomes, as is NRSF itself. Some NRSF/REST sites reside in repeats, which suggests a mechanism for both ancient and modern dispersal of NRSEs through vertebrate genomes. Multiple predicted sites are located near neuronal microRNA and splicing-factor genes, and these tested positive for NRSF/REST occupancy in vivo. The resulting network model integrates post-transcriptional and translational controllers, including candidate feedback loops on NRSF and its corepressor, CoREST.

Evolutionary expansion and divergence in the ZNF91 subfamily of primate-specific zinc finger genes

Most genes are conserved in mammals, but certain gene families have acquired large numbers of lineage-specific loci through repeated rounds of gene duplication, divergence, and loss that have continued in each mammalian group. One such family encodes KRAB-zinc finger (KRAB-ZNF) proteins, which function as transcriptional repressors. One particular subfamily of KRAB-ZNF genes, including ZNF91, has expanded specifically in primates to comprise more than 110 loci in the human genome. Genes of the ZNF91 subfamily reside in large gene clusters near centromeric regions of human chromosomes 19 and 7 with smaller clusters or isolated copies in other locations. Phylogenetic analysis indicates that many of these genes arose before the split between the New and Old World monkeys, but the ZNF91 subfamily has continued to expand and diversify throughout the evolution of apes and humans. Paralogous loci are distinguished by divergence within their zinc finger arrays, indicating selection for proteins with different regulatory targets. In addition, many loci produce multiple alternatively spliced transcripts encoding proteins that may serve separate and perhaps even opposing regulatory roles because of the modular motif structure of KRAB-ZNF genes. The tissue-specific expression patterns and rapid structural divergence of ZNF91 subfamily genes suggest a role in determining gene expression differences between species and the evolution of novel primate traits.

A comprehensive catalog of human KRAB-associated zinc finger genes: insights into the evolutionary history of a large family of transcriptional repressors

Krüppel-type zinc finger (ZNF) motifs are prevalent components of transcription factor proteins in all eukaryotes. KRAB-ZNF proteins, in which a potent repressor domain is attached to a tandem array of DNA-binding zinc-finger motifs, are specific to tetrapod vertebrates and represent the largest class of ZNF proteins in mammals. To define the full repertoire of human KRAB-ZNF proteins, we searched the genome sequence for key motifs and then constructed and manually curated gene models incorporating those sequences. The resulting gene catalog contains 423 KRAB-ZNF protein-coding loci, yielding alternative transcripts that altogether predict at least 742 structurally distinct proteins. Active rounds of segmental duplication, involving single genes or larger regions and including both tandem and distributed duplication events, have driven the expansion of this mammalian gene family. Comparisons between the human genes and ZNF loci mined from the draft mouse, dog, and chimpanzee genomes not only identified 103 KRAB-ZNF genes that are conserved in mammals but also highlighted a substantial level of lineage-specific change; at least 136 KRAB-ZNF coding genes are primate specific, including many recent duplicates. KRAB-ZNF genes are widely expressed and clustered genes are typically not coregulated, indicating that paralogs have evolved to fill roles in many different biological processes. To facilitate further study, we have developed a Web-based public resource with access to gene models, sequences, and other data, including visualization tools to provide genomic context and interaction with other public data sets.