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

Interplay of Transcriptomic Regulation, Microbiota, and Signaling Pathways in Lung and Gut Inflammation-Induced Tumorigenesis

Inflammation can positively and negatively affect tumorigenesis based on the duration, scope, and sequence of related events through the regulation of signaling pathways. A transcriptomic analysis of five pulmonary arterial hypertension, twelve Crohn's disease, and twelve ulcerative colitis high throughput sequencing datasets using R language specialized libraries and gene enrichment analyses identified a regulatory network in each inflammatory disease. IRF9 and LINC01089 in pulmonary arterial hypertension are related to the regulation of signaling pathways like MAPK, NOTCH, human papillomavirus, and hepatitis c infection. ZNF91 and TP53TG1 in Crohn's disease are related to the regulation of PPAR, MAPK, and metabolic signaling pathways. ZNF91, VDR, DLEU1, SATB2-AS1, and TP53TG1 in ulcerative colitis are related to the regulation of PPAR, AMPK, and metabolic signaling pathways. The activation of the transcriptomic network and signaling pathways might be related to the interaction of the characteristic microbiota of the inflammatory disease, with the lung and gut cell receptors present in membrane rafts and complexes. The transcriptomic analysis highlights the impact of several coding and non-coding RNAs, suggesting their relationship with the unlocking of cell phenotypic plasticity for the acquisition of the hallmarks of cancer during lung and gut cell adaptation to inflammatory phenotypes.

Gene expression networks regulated by human personality

Genome-wide association studies of human personality have been carried out, but transcription of the whole genome has not been studied in relation to personality in humans. We collected genome-wide expression profiles of adults to characterize the regulation of expression and function in genes related to human personality. We devised an innovative multi-omic approach to network analysis to identify the key control elements and interactions in multi-modular networks. We identified sets of transcribed genes that were co-expressed in specific brain regions with genes known to be associated with personality. Then we identified the minimum networks for the co-localized genes using bioinformatic resources. Subjects were 459 adults from the Young Finns Study who completed the Temperament and Character Inventory and provided peripheral blood for genomic and transcriptomic analysis. We identified an extrinsic network of 45 regulatory genes from seed genes in brain regions involved in self-regulation of emotional reactivity to extracellular stimuli (e.g., self-regulation of anxiety) and an intrinsic network of 43 regulatory genes from seed genes in brain regions involved in self-regulation of interpretations of meaning (e.g., production of concepts and language). We discovered that interactions between the two networks were coordinated by a control hub of 3 miRNAs and 3 protein-coding genes shared by both. Interactions of the control hub with proteins and ncRNAs identified more than 100 genes that overlap directly with known personality-related genes and more than another 4000 genes that interact indirectly. We conclude that the six-gene hub is the crux of an integrative network that orchestrates information-transfer throughout a multi-modular system of over 4000 genes enriched in liquid-liquid-phase-separation (LLPS)-related RNAs, diverse transcription factors, and hominid-specific miRNAs and lncRNAs. Gene expression networks associated with human personality regulate neuronal plasticity, epigenesis, and adaptive functioning by the interactions of salience and meaning in self-awareness.

Emergence of the ZNF675 Gene During Primate Evolution-Influenced Human Neurodevelopment Through Changing HES1 Autoregulation

In this study, we investigated recurrent copy number variations (CNVs) in the 19p12 locus, which are associated with neurodevelopmental disorders. The two genes in this locus, ZNF675 and ZNF681, arose via gene duplication in primates, and their presence in several pathological CNVs in the human population suggests that either or both of these genes are required for normal human brain development. ZNF675 and ZNF681 are members of the Krüppel-associated box zinc finger (KZNF) protein family, a class of transcriptional repressors important for epigenetic silencing of specific genomic regions. About 170 primate-specific KZNFs are present in the human genome. Although KZNFs are primarily associated with repressing retrotransposon-derived DNA, evidence is emerging that they can be co-opted for other gene regulatory processes. We show that genetic deletion of ZNF675 causes developmental defects in cortical organoids, and our data suggest that part of the observed neurodevelopmental phenotype is mediated by a gene regulatory role of ZNF675 on the promoter of the neurodevelopmental gene Hes family BHLH transcription factor 1 (HES1). We also find evidence for the recently evolved regulation of genes involved in neurological disorders, microcephalin 1 and sestrin 3. We show that ZNF675 interferes with HES1 auto-inhibition, a process essential for the maintenance of neural progenitors. As a striking example of how some KZNFs have integrated into preexisting gene expression networks, these findings suggest the emergence of ZNF675 has caused a change in the balance of HES1 autoregulation. The association of ZNF675 CNV with human developmental disorders and ZNF675-mediated regulation of neurodevelopmental genes suggests that it evolved into an important factor for human brain development.

Hominoid SVA-lncRNA AK057321 targets human-specific SVA retrotransposons in SCN8A and CDK5RAP2 to initiate neuronal maturation

SINE-VNTR-Alu (SVA) retrotransposons arose and expanded in the genome of hominoid primates concurrent with the slowing of brain maturation. We report genes with intronic SVA transposons are enriched for neurodevelopmental disease and transcribed into long non-coding SVA-lncRNAs. Human-specific SVAs in microcephaly CDK5RAP2 and epilepsy SCN8A gene introns repress their expression via transcription factor ZNF91 to delay neuronal maturation. Deleting the SVA in CDK5RAP2 initiates multi-dimensional and in SCN8A selective sodium current neuronal maturation by upregulating these genes. SVA-lncRNA AK057321 forms RNA:DNA heteroduplexes with the genomic SVAs and upregulates these genes to initiate neuronal maturation. SVA-lncRNA AK057321 also promotes species-specific cortex and cerebellum-enriched expression upregulating human genes with intronic SVAs (e.g., HTT, CHAF1B and KCNJ6) but not mouse orthologs. The diversity of neuronal genes with intronic SVAs suggest this hominoid-specific SVA transposon-based gene regulatory mechanism may act at multiple steps to specialize and achieve neoteny of the human brain.

Primate protein-ligand interfaces exhibit significant conservation and unveil human-specific evolutionary drivers

Despite the vast phenotypic differences observed across primates, their protein products are largely similar to each other at the sequence level. We hypothesized that, since proteins accomplish all their functions via interactions with other molecules, alterations in the sites that participate in these interactions may be of critical importance. To uncover the extent to which these sites evolve across primates, we built a structurally-derived dataset of ~4,200 one-to-one orthologous sequence groups across 18 primate species, consisting of ~68,000 ligand-binding sites that interact with DNA, RNA, small molecules, ions, or peptides. Using this dataset, we identify functionally important patterns of conservation and variation within the amino acid residues that facilitate protein-ligand interactions across the primate phylogeny. We uncover that interaction sites are significantly more conserved than other sites, and that sites binding DNA and RNA further exhibit the lowest levels of variation. We also show that the subset of ligand-binding sites that do vary are enriched in components of gene regulatory pathways and uncover several instances of human-specific ligand-binding site changes within transcription factors. Altogether, our results suggest that ligand-binding sites have experienced selective pressure in primates and propose that variation in these sites may have an outsized effect on phenotypic variation in primates through pleiotropic effects on gene regulation.

Identification of genes contributing to cisplatin resistance in osteosarcoma cells

Osteosarcomas are prevalent in children and young adults and have a high recurrence rate. Cisplatin, doxorubicin, and methotrexate are common adjuvant chemotherapy drugs for treatment of osteosarcoma, but multidrug resistance is a growing problem. Therefore, understanding the molecular mechanisms of chemotherapy resistance in osteosarcoma cells is crucial for developing new therapeutic approaches and ultimately improving the prognosis of osteosarcoma patients. To identify genes associated with cisplatin resistance in osteosarcoma, we screened a large-scale mutant library generated by transfecting human osteosarcoma cells with a piggyBac (PB) transposon-based gene activation vector. Several candidate genes were identified by using Splinkerette-PCR paired with Next Generation Sequencing. We created a disease-free survival predictor model, which includes ZNF720, REEP3, CNNM2, and CGREF1, using TARGET (Therapeutically Applicable Research to Generate Effective Treatments) datasets. Additionally, the results of our enrichment analysis between the Four_genes_high group and Low_group suggested that these four genes may participate in cisplatin resistance in osteosarcoma through cross talk between various signaling pathways, especially the signaling pathway related to bone formation. These data may help guide future studies into chemotherapy for osteosarcoma.

Prognostic Relevance of ZNF844 and Chr 19p13.2 KRAB-Zinc Finger Proteins in Clear Cell Renal Carcinoma

BACKGROUND/AIM

Clear-cell renal cell carcinoma (ccRCC) is the most common and aggressive form of all urological cancers, with poor prognosis and high mortality. Despite growing evidence of involvement in carcinogenesis, the role of KRAB-ZFP in ccRCC has not been fully explored. KRAB Zinc finger proteins (KRAB-ZFPs) are the largest family of mammalian transcription regulators. They are differentially expressed in various tissues during cellular development and phenotypic differentiation.

MATERIALS AND METHODS

In this study, the levels of transcripts of ccRCC from The Cancer Genome Atlas (TCGA) dataset were used to identify prognostic biomarkers in this disease.

RESULTS

Using bioinformatics techniques, we demonstrate that approximately 60% of KRAB zinc finger proteins located on chromosome 19p13.2 are differentially expressed, with all but two being down-regulated in ccRCC. Moreover, ZNF844, a paralog of ZNF433, was the most down-regulated across all histological grades and pathological stages (p<0.001). In addition, the decrease in ZNF844 expression was associated with poor patient survival (HR=0.41; 95% CI=0.3-0.56; p<0.0001). Gene Set Enrichment Analysis of genes inversely co-expressed with ZNF844 revealed that enriched pathways were consistently related to immune and translation processes (p<0.05, FDR <0.05). Lastly, ZNF844 expression showed moderate, inverse correlation to Helper T-cell (CD4 or Th1) subtype 1 (R=-0.558, p=5.15×10^-39) infiltration and with the exhausted T-cell phenotype (R=-0.37; p=4.1×10^-21).

CONCLUSION

Down-regulation of KRAB-ZFPs at 19p13.2 may represent a signature for ccRCC. Moreover, ZNF844 is a prognostic marker for ccRCC and may serve as a putative immune-related tumor suppressor gene.

Lack of Molecular Mimicry between Nonhuman Primates and Infectious Pathogens: The Possible Genetic Bases

Recently, it was found that proteomes from poliovirus, measles virus, dengue virus, and severe acute respiratory syndrome-related Coronavirus 2 (SARS-CoV-2) have high molecular mimicry at the heptapeptide level with the human proteome, while heptapeptide commonality is minimal or absent with proteomes from nonhuman primates, that is, gorilla, chimpanzee, and rhesus macaque. To acquire more data on the issue, analyses here have been expanded to Ebola virus, Francisella tularensis , human immunodeficiency virus-1 (HIV-1), Toxoplasma gondii , Variola virus, and Yersinia pestis . Results confirm that heptapeptide overlap is high between pathogens and Homo sapiens , but not between pathogens and primates. Data are discussed in light of the possible genetic bases that differently model primate phenomes, thus possibly underlying the zero/low level of molecular mimicry between infectious agents and primates. Notably, this study might help address preclinical vaccine tests that currently utilize primates as animal models, since autoimmune cross-reactions and the consequent adverse events cannot occur in absentia of shared sequences.

Conserved paradoxical relationships among the evolutionary, structural and expressional features of KRAB zinc-finger proteins reveal their special functional characteristics

Background

One striking feature of the large KRAB domain-containing zinc finger protein (KZFP) family is its rapid evolution, leading to hundreds of member genes with various origination time in a certain mammalian genome. However, a comprehensive genome-wide and across-taxa analysis of the structural and expressional features of KZFPs with different origination time is lacking. This type of analysis will provide valuable clues about the functional characteristics of this special family.

Results

In this study, we found several conserved paradoxical phenomena about this issue. 1) Ordinary young domains/proteins tend to be disordered, but most of KRAB domains are completely structured in 64 representative species across the superclass of Sarcopterygii and most of KZFPs are also highly structured, indicating their rigid and unique structural and functional characteristics; as exceptions, old-zinc-finger-containing KZFPs have relatively disordered KRAB domains and linker regions, contributing to diverse interacting partners and functions. 2) In general, young or highly structured proteins tend to be spatiotemporal specific and have low abundance. However, by integrated analysis of 29 RNA-seq datasets, including 725 samples across early embryonic development, embryonic stem cell differentiation, embryonic and adult organs, tissues in 7 mammals, we found that KZFPs tend to express ubiquitously with medium abundance regardless of evolutionary age and structural disorder degree, indicating the wide functional requirements of KZFPs in various states. 3) Clustering and correlation analysis reveal that there are differential expression patterns across different spatiotemporal states, suggesting the specific-high-expression KZFPs may play important roles in the corresponding states. In particular, part of young-zinc-finger-containing KZFPs are highly expressed in early embryonic development and ESCs differentiation into endoderm or mesoderm. Co-expression analysis revealed that young-zinc-finger-containing KZFPs are significantly enriched in five co-expression modules. Among them, one module, including 13 young-zinc-finger-containing KZFPs, showed an ‘early-high and late-low’ expression pattern. Further functional analysis revealed that they may function in early embryonic development and ESC differentiation via participating in cell cycle related processes.

Conclusions

This study shows the conserved and special structural, expressional features of KZFPs, providing new clues about their functional characteristics and potential causes of their rapid evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00346-w.

Few SINEs of life: Alu elements have little evidence for biological relevance despite elevated translation

Transposable elements colonize genomes and with time may end up being incorporated into functional regions. SINE Alu elements, which appeared in the primate lineage, are ubiquitous in the human genome and more than a thousand overlap annotated coding exons. Although almost all Alu-derived coding exons appear to be in alternative transcripts, they have been incorporated into the main coding transcript in at least 11 genes. The extent to which Alu regions are incorporated into functional proteins is unclear, but we detected reliable peptide evidence to support the translation to protein of 33 Alu-derived exons. All but one of the Alu elements for which we detected peptides were frame-preserving and there was proportionally seven times more peptide evidence for Alu elements as for other primate exons. Despite this strong evidence for translation to protein we found no evidence of selection, either from cross species alignments or human population variation data, among these Alu-derived exons. Overall, our results confirm that SINE Alu elements have contributed to the expansion of the human proteome, and this contribution appears to be stronger than might be expected over such a relatively short evolutionary timeframe. Despite this, the biological relevance of these modifications remains open to question.

The Role of SINE-VNTR-Alu (SVA) Retrotransposons in Shaping the Human Genome

Retrotransposons can alter the regulation of genes both transcriptionally and post-transcriptionally, through mechanisms such as binding transcription factors and alternative splicing of transcripts. SINE-VNTR-Alu (SVA) retrotransposons are the most recently evolved class of retrotransposable elements, found solely in primates, including humans. SVAs are preferentially found at genic, high GC loci, and have been termed "mobile CpG islands". We hypothesise that the ability of SVAs to mobilise, and their non-random distribution across the genome, may result in differential regulation of certain pathways. We analysed SVA distribution patterns across the human reference genome and identified over-representation of SVAs at zinc finger gene clusters. Zinc finger proteins are able to bind to and repress SVA function through transcriptional and epigenetic mechanisms, and the interplay between SVAs and zinc fingers has been proposed as a major feature of genome evolution. We describe observations relating to the clustering patterns of both reference SVAs and polymorphic SVA insertions at zinc finger gene loci, suggesting that the evolution of this network may be ongoing in humans. Further, we propose a mechanism to direct future research and validation efforts, in which the interplay between zinc fingers and their epigenetic modulation of SVAs may regulate a network of zinc finger genes, with the potential for wider transcriptional consequences.

Identification of a nine-gene panel as a prognostic indicator for recurrence with muscle-invasive bladder cancer

BACKGROUND AND OBJECTIVES

Bladder cancer is one of the most common and highly recurrent cancers worldwide. Recurrence-associated genes may potentially predict cancer recurrence. We aimed to construct a recurrence-associated gene panel to improve the prognostic prediction of bladder cancer.

METHODS

Based on DNA sequencing and clinical data from the TCGA-BLCA project, we identified 10 potential driver genes significantly associated with recurrence of bladder cancer. We performed multivariable logistic regression analysis to construct an optimized recurrence prediction model with nine recurrence-associated genes (EME1, AKAP9, ZNF91, PARD3, STAG2, ZFP36L2, METTL3, POLR3B, and MUC7) and clinical information as the independent variables.

RESULTS

The area under the receiver operating characteristic (ROC) curve was 0.80 in this model, much higher than that of the baseline model (AUC = 0.73) and the same trend was also validated in its subset. Decision curve analysis also revealed that there is a significant net benefit gained by adding nine genes mutation to the baseline model. Furthermore, Kaplan-Meier survival analysis showed that eight out of the nine genes (excluding MUC7) had good effects on the overall prognosis of patients.

CONCLUSIONS

This nine-gene panel will most likely be a useful tool for prognostic evaluation and will facilitate the personalized management of patients with bladder cancer.

Endogenous Retroviruses Function as Gene Expression Regulatory Elements During Mammalian Pre-implantation Embryo Development

Pre-implantation embryo development encompasses several key developmental events, especially the activation of zygotic genome activation (ZGA)-related genes. Endogenous retroviruses (ERVs), which are regarded as “deleterious genomic parasites”, were previously considered to be “junk DNA”. However, it is now known that ERVs, with limited conservatism across species, mediate conserved developmental processes (e.g., ZGA). Transcriptional activation of ERVs occurs during the transition from maternal control to zygotic genome control, signifying ZGA. ERVs are versatile participants in rewiring gene expression networks during epigenetic reprogramming. Particularly, a subtle balance exists between ERV activation and ERV repression in host–virus interplay, which leads to stage-specific ERV expression during pre-implantation embryo development. A large portion of somatic cell nuclear transfer (SCNT) embryos display developmental arrest and ZGA failure during pre-implantation embryo development. Furthermore, because of the close relationship between ERV activation and ZGA, exploring the regulatory mechanism underlying ERV activation may also shed more light on the enigma of SCNT embryo development in model animals.

Non-base-contacting residues enable kaleidoscopic evolution of metazoan C2H2 zinc finger DNA binding

Background

The C2H2 zinc finger (C2H2-ZF) is the most numerous protein domain in many metazoans, but is not as frequent or diverse in other eukaryotes. The biochemical and evolutionary mechanisms that underlie the diversity of this DNA-binding domain exclusively in metazoans are, however, mostly unknown.

Results

Here, we show that the C2H2-ZF expansion in metazoans is facilitated by contribution of non-base-contacting residues to DNA binding energy, allowing base-contacting specificity residues to mutate without catastrophic loss of DNA binding. In contrast, C2H2-ZF DNA binding in fungi, plants, and other lineages is constrained by reliance on base-contacting residues for DNA-binding functionality. Reconstructions indicate that virtually every DNA triplet was recognized by at least one C2H2-ZF domain in the common progenitor of placental mammals, but that extant C2H2-ZF domains typically bind different sequences from these ancestral domains, with changes facilitated by non-base-contacting residues.

Conclusions

Our results suggest that the evolution of C2H2-ZFs in metazoans was expedited by the interaction of non-base-contacting residues with the DNA backbone. We term this phenomenon “kaleidoscopic evolution,” to reflect the diversity of both binding motifs and binding motif transitions and the facilitation of their diversification.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1287-y) contains supplementary material, which is available to authorized users.

The pioneering role of PRDM9 indel mutations in tarsier evolution

PRDM9 is currently the sole speciation gene found in vertebrates causing hybrid sterility probably due to incompatible alleles. Its role in defining the double strand break loci during the meiotic prophase I is crucial for proper chromosome segregation. Therefore, the rapid turnover of the loci determining zinc finger array seems to be causative for incompatibilities. We here investigated the zinc finger domain-containing exon of PRDM9 in 23 tarsiers. Tarsiers, the most basal extant haplorhine primates, exhibit two frameshifting indels at the 5'-end of the array. The first mutation event interrupts the reading frame and function while the second compensates both. The fixation of this allele variant in tarsiers led to hypothesize that de- and reactivation of the zinc finger domain drove the speciation in early haplorhine or tarsiiform primates. Moreover, the high allelic diversity within Tarsius points to multiple effects of genetic drift reflecting their phylogeographic history since the Miocene.

Transcriptional repression of p27 is essential for murine embryonic development

The Nczf gene has been identified as one of Ncx target genes and encodes a novel KRAB zinc-finger protein, which functions as a sequence specific transcriptional repressor. In order to elucidate Nczf functions, we generated Nczf knockout (Nczf−/−) mice. Nczf−/− mice died around embryonic day 8.5 (E8.5) with small body size and impairment of axial rotation. Histopathological analysis revealed that the cell number decreased and pyknotic cells were occasionally observed. We examined the expression of cell cycle related genes in Nczf−/− mice. p27 expression was increased in E8.0 Nczf−/− mice compared to that of wild type mice. Nczf knockdown by siRNA resulted in increased expression of p27 in mouse embryonic fibroblasts (MEFs). Furthermore, p27 promoter luciferase reporter gene analysis confirmed the regulation of p27 mRNA expression by Nczf. Nczf−/−; p27−/− double knockout mice survived until E11.5 and the defect of axial rotation was restored. These data suggest that p27 repression by Nczf is essential in the developing embryo.

Human Lineage-Specific Transcriptional Regulation through GA-Binding Protein Transcription Factor Alpha (GABPa)

A substantial fraction of phenotypic differences between closely related species are likely caused by differences in gene regulation. While this has already been postulated over 30 years ago, only few examples of evolutionary changes in gene regulation have been verified. Here, we identified and investigated binding sites of the transcription factor GA-binding protein alpha (GABPa) aiming to discover cis-regulatory adaptations on the human lineage. By performing chromatin immunoprecipitation-sequencing experiments in a human cell line, we found 11,619 putative GABPa binding sites. Through sequence comparisons of the human GABPa binding regions with orthologous sequences from 34 mammals, we identified substitutions that have resulted in 224 putative human-specific GABPa binding sites. To experimentally assess the transcriptional impact of those substitutions, we selected four promoters for promoter-reporter gene assays using human and African green monkey cells. We compared the activities of wild-type promoters to mutated forms, where we have introduced one or more substitutions to mimic the ancestral state devoid of the GABPa consensus binding sequence. Similarly, we introduced the human-specific substitutions into chimpanzee and macaque promoter backgrounds. Our results demonstrate that the identified substitutions are functional, both in human and nonhuman promoters. In addition, we performed GABPa knock-down experiments and found 1,215 genes as strong candidates for primary targets. Further analyses of our data sets link GABPa to cognitive disorders, diabetes, KRAB zinc finger (KRAB-ZNF), and human-specific genes. Thus, we propose that differences in GABPa binding sites played important roles in the evolution of human-specific phenotypes.

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.

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.

Spotting the enemy within: Targeted silencing of foreign DNA in mammalian genomes by the Krüppel-associated box zinc finger protein family

Tandem C2H2-type zinc finger proteins (ZFPs) constitute the largest transcription factor family in animals. Tandem-ZFPs bind DNA in a sequence-specific manner through arrays of multiple zinc finger domains that allow high flexibility and specificity in target recognition. In tetrapods, a large proportion of tandem-ZFPs contain Krüppel-associated-box (KRAB) repression domains, which are able to induce epigenetic silencing through the KAP1 corepressor. The KRAB-ZFP family continuously amplified in tetrapods through segmental gene duplications, often accompanied by deletions, duplications, and mutations of the zinc finger domains. As a result, tetrapod genomes contain unique sets of KRAB-ZFP genes, consisting of ancient and recently evolved family members. Although several hundred human and mouse KRAB-ZFPs have been identified or predicted, the biological functions of most KRAB-ZFP family members have gone unexplored. Furthermore, the evolutionary forces driving the extraordinary KRAB-ZFP expansion and diversification have remained mysterious for decades. In this review, we highlight recent studies that associate KRAB-ZFPs with the repression of parasitic DNA elements in the mammalian germ line and discuss the hypothesis that the KRAB-ZFP family primarily evolved as an adaptive genomic surveillance system against foreign DNA. Finally, we comment on the computational, genetic, and biochemical challenges of studying KRAB-ZFPs and attempt to predict how these challenges may be soon overcome.

Retrotransposons shape species-specific embryonic stem cell gene expression

Over half of our genome is composed of retrotransposons, which are mobile elements that can readily amplify their copy number by replicating through an RNA intermediate. Most of these elements are no longer mobile but still contain regulatory sequences that can serve as promoters, enhancers or repressors for cellular genes. Despite dominating our genetic content, little is known about the precise functions of retrotransposons, which include both endogenous retroviruses (ERVs) and non-LTR elements like long interspersed nuclear element 1 (LINE-1). However, a few recent cutting-edge publications have illustrated how retrotransposons shape species-specific stem cell gene expression by two opposing mechanisms, involving their recruitment of stem cell-enriched transcription factors (TFs): firstly, they can activate expression of genes linked to naïve pluripotency, and secondly, they can induce repression of proximal genes. The paradox that different retrotransposons are active or silent in embryonic stem cells (ESCs) can be explained by differences between retrotransposon families, between individual copies within the same family, and between subpopulations of ESCs. Since they have coevolved with their host genomes, some of them have been co-opted to perform species-specific beneficial functions, while others have been implicated in genetic disease. In this review, we will discuss retrotransposon functions in ESCs, focusing on recent mechanistic advances of how HERV-H has been adopted to preserve human naïve pluripotency and how particular LINE-1, SVA and ERV family members recruit species-specific transcriptional repressors. This review highlights the fine balance between activation and repression of retrotransposons that exists to harness their ability to drive evolution, while minimizing the risk they pose to genome integrity.

C2H2 zinc finger proteins greatly expand the human regulatory lexicon

Cys2-His2 zinc finger (C2H2-ZF) proteins represent the largest class of putative human transcription factors. However, for most C2H2-ZF proteins it is unknown whether they even bind DNA or, if they do, to which sequences. Here, by combining data from a modified bacterial one-hybrid system with protein-binding microarray and chromatin immunoprecipitation analyses, we show that natural C2H2-ZFs encoded in the human genome bind DNA both in vitro and in vivo, and we infer the DNA recognition code using DNA-binding data for thousands of natural C2H2-ZF domains. In vivo binding data are generally consistent with our recognition code and indicate that C2H2-ZF proteins recognize more motifs than all other human transcription factors combined. We provide direct evidence that most KRAB-containing C2H2-ZF proteins bind specific endogenous retroelements (EREs), ranging from currently active to ancient families. The majority of C2H2-ZF proteins, including KRAB proteins, also show widespread binding to regulatory regions, indicating that the human genome contains an extensive and largely unstudied adaptive C2H2-ZF regulatory network that targets a diverse range of genes and pathways.

Phylogenetic analyses provide the first insights into the evolution of OVATE family proteins in land plants

BACKGROUND AND AIMS

The OVATE gene encodes a nuclear-localized regulatory protein belonging to a distinct family of plant-specific proteins known as the OVATE family proteins (OFPs). OVATE was first identified as a key regulator of fruit shape in tomato, with nonsense mutants displaying pear-shaped fruits. However, the role of OFPs in plant development has been poorly characterized.

METHODS

Public databases were searched and a total of 265 putative OVATE protein sequences were identified from 13 sequenced plant genomes that represent the major evolutionary lineages of land plants. A phylogenetic analysis was conducted based on the alignment of the conserved OVATE domain from these 13 selected plant genomes. The expression patterns of tomato SlOFP genes were analysed via quantitative real-time PCR. The pattern of OVATE gene duplication resulting in the expansion of the gene family was determined in arabidopsis, rice and tomato.

KEY RESULTS

Genes for OFPs were found to be present in all the sampled land plant genomes, including the early-diverged lineages, mosses and lycophytes. Phylogenetic analysis based on the amino acid sequences of the conserved OVATE domain defined 11 sub-groups of OFPs in angiosperms. Different evolutionary mechanisms are proposed for OVATE family evolution, namely conserved evolution and divergent expansion. Characterization of the AtOFP family in arabidopsis, the OsOFP family in rice and the SlOFP family in tomato provided further details regarding the evolutionary framework and revealed a major contribution of tandem and segmental duplications towards expansion of the OVATE gene family.

CONCLUSIONS

This first genome-wide survey on OFPs provides new insights into the evolution of the OVATE protein family and establishes a solid base for future functional genomics studies on this important but poorly characterized regulatory protein family in plants.

Accumulation of interspersed and sex-specific repeats in the non-recombining region of papaya sex chromosomes

Background

The papaya Y chromosome has undergone a degenerative expansion from its ancestral autosome, as a consequence of recombination suppression in the sex determining region of the sex chromosomes. The non-recombining feature led to the accumulation of repetitive sequences in the male- or hermaphrodite-specific regions of the Y or the Y^h chromosome (MSY or HSY). Therefore, repeat composition and distribution in the sex determining region of papaya sex chromosomes would be informative to understand how these repetitive sequences might be involved in the early stages of sex chromosome evolution.

Results

Detailed composition of interspersed, sex-specific, and tandem repeats was analyzed from 8.1 megabases (Mb) HSY and 5.3 Mb corresponding X chromosomal regions. Approximately 77% of the HSY and 64% of the corresponding X region were occupied by repetitive sequences. Ty3-gypsy retrotransposons were the most abundant interspersed repeats in both regions. Comparative analysis of repetitive sequences between the sex determining region of papaya X chromosome and orthologous autosomal sequences of Vasconcellea monoica, a close relative of papaya lacking sex chromosomes, revealed distinctive differences in the accumulation of Ty3-Gypsy, suggesting that the evolution of the papaya sex determining region may accompany Ty3-Gypsy element accumulation. In total, 21 sex-specific repeats were identified from the sex determining region; 20 from the HSY and one from the X. Interestingly, most HSY-specific repeats were detected in two regions where the HSY expansion occurred, suggesting that the HSY expansion may result in the accumulation of sex-specific repeats or that HSY-specific repeats might play an important role in the HSY expansion. The analysis of simple sequence repeats (SSRs) revealed that longer SSRs were less abundant in the papaya sex determining region than the other chromosomal regions.

Conclusion

Major repetitive elements were Ty3-gypsy retrotransposons in both the HSY and the corresponding X. Accumulation of Ty3-Gypsy retrotransposons in the sex determining region of papaya X chromosome was significantly higher than that in the corresponding region of V. monoica, suggesting that Ty3-Gypsy could be crucial for the expansion and evolution of the sex determining region in papaya. Most sex-specific repeats were located in the two HSY expansion regions.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-335) contains supplementary material, which is available to authorized users.

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.

ZNF-mediated resistance to imatinib mesylate in gastrointestinal stromal tumor

Although imatinib mesylate (IM) has transformed the treatment of gastrointestinal stromal tumors (GIST), many patients experience primary/secondary drug resistance. In a previous study, we identified a gene signature, consisting mainly of Kruppel-associated box (KRAB) domain containing zinc finger (ZNF) transcriptional repressors that predict short-term response to IM. To determine if these genes have functional significance, a siRNA library targeting these genes was constructed and applied to GIST cells in vitro. These screens identified seventeen “IM sensitizing genes” in GIST cells (sensitization index (SI) <0.85 ratio of drug/vehicle) with a false discovery rate (FDR) <15%, including twelve ZNF genes, the majority of which are located within the HSA19p12–13.1 locus. These genes were shown to be highly specific to IM and another tyrosine kinase inhibitor (TKI), sunitinib, in GIST cells. In order to determine mechanistically how these ZNFs might be modulating response to IM, RNAi approaches were used to individually silence genes within the predictive signature in GIST cells and expression profiling was performed. Knockdown of the 14 IM-sensitizing genes (10 ZNFs) universally led to downregulation of six genes, including TGFb3, periostin, and NEDD9. These studies implicate a role of KRAB-ZNFs in modulating response to TKIs in GIST.

A common mode of recognition for methylated CpG

Much is known about vertebrate DNA methylation, however it is not known how methylated CpG within particular sequences is recognized. Two recent structures of C2H2 zinc finger (ZnF) proteins in complex with methylated DNA reveal a common recognition mode for 5-methylcytosine (5mC) that involves a 5mC-Arg-G triad. In the two ZnF proteins, an arginine that precedes the first Zn-binding histidine (RH motif) can interact with a 5mCpG or TpG dinucleotide. Among a family of >300 human Krüppel-associated box (KRAB) domain containing ZnF proteins examined, two-thirds contained at least one ZnF that included an RH motif. We propose that the RH-ZnF motifs provide specificity for 5mCpG, whereas the neighboring Zn fingers recognize the surrounding DNA sequence context.

A prominent role of KRAB-ZNF transcription factors in mammalian speciation?

The mechanisms of speciation have been one of the most debated topics in evolutionary biology. Among all reproductive barriers, postzygotic reproductive isolation is perhaps the one that has attracted the most attention from geneticists. Despite remarkable advances in the identification of loci involved in Drosophila speciation, little is known about the genes, functions, and biochemical interactions of the molecules underlying hybrid sterility and inviability in mammals. Here, we discuss the main evolutionary and molecular features that make transcription factors (TFs), especially the family of zinc finger proteins with a Krüppel-associated box domain (KRAB-ZNF), strong candidates to play an important role in postzygotic reproductive isolation. Motivated by the recent identification of the gene encoding PR domain zinc finger protein 9 (Prdm9; a KRAB-ZNF gene) as the first hybrid sterility gene identified in mammals, we further propose integrative approaches to study KRAB-ZNF genes with the main goal of characterizing the molecular pathways and interactions involved in hybrid incompatibilities.

Gypsy and the birth of the SCAN domain

SCAN is a protein domain frequently found at the N termini of proteins encoded by mammalian tandem zinc finger (ZF) genes, whose structure is known to be similar to that of retroviral gag capsid domains and whose multimerization has been proposed as a model for retroviral assembly. We report that the SCAN domain is derived from the C-terminal portion of the gag capsid (CA) protein from the Gmr1-like family of Gypsy/Ty3-like retrotransposons. On the basis of sequence alignments and phylogenetic distributions, we show that the ancestral host SCAN domain (ESCAN for extended SCAN) was exapted from a full-length CA gene from a Gmr1-like retrotransposon at or near the root of the tetrapod animal branch. A truncated variant of ESCAN that corresponds to the annotated SCAN domain arose shortly thereafter and appears to be the only form extant in mammals. The Anolis lizard has a large number of tandem ZF genes with N-terminal ESCAN or SCAN domains. We predict DNA binding sites for all Anolis ESCAN-ZF and SCAN-ZF proteins and demonstrate several highly significant matches to Anolis Gmr1-like sequences, suggesting that at least some of these proteins target retroelements. SCAN is known to mediate protein dimerization, and the CA protein multimerizes to form the core retroviral and retrotransposon capsid structure. We speculate that the SCAN domain originally functioned to target host ZF proteins to retroelement capsids.

Reassessing domain architecture evolution of metazoan proteins: major impact of errors caused by confusing paralogs and epaktologs

In the accompanying paper (Nagy, Szláma, Szarka, Trexler, Bányai, Patthy, Reassessing Domain Architecture Evolution of Metazoan Proteins: Major Impact of Gene Prediction Errors) we showed that in the case of UniProtKB/TrEMBL, RefSeq, EnsEMBL and NCBI's GNOMON predicted protein sequences of Metazoan species the contribution of erroneous (incomplete, abnormal, mispredicted) sequences to domain architecture (DA) differences of orthologous proteins might be greater than those of true gene rearrangements. Based on these findings, we suggest that earlier genome-scale studies based on comparison of predicted (frequently mispredicted) protein sequences may have led to some erroneous conclusions about the evolution of novel domain architectures of multidomain proteins. In this manuscript we examine the impact of confusing paralogous and epaktologous multidomain proteins (i.e., those that are related only through the independent acquisition of the same domain types) on conclusions drawn about DA evolution of multidomain proteins in Metazoa. To estimate the contribution of this type of error we have used as reference UniProtKB/Swiss-Prot sequences from protein families with well-characterized evolutionary histories. We have used two types of paralogy-group construction procedures and monitored the impact of various parameters on the separation of true paralogs from epaktologs on correctly annotated Swiss-Prot entries of multidomain proteins. Our studies have shown that, although public protein family databases are contaminated with epaktologs, analysis of the structure of sequence similarity networks of multidomain proteins provides an efficient means for the separation of epaktologs and paralogs. We have also demonstrated that contamination of protein families with epaktologs increases the apparent rate of DA change and introduces a bias in DA differences in as much as it increases the proportion of terminal over internal DA differences. We have shown that confusing paralogous and epaktologous multidomain proteins significantly increases the apparent rate of DA change in Metazoa and introduces a positional bias in favor of terminal over internal DA changes. Our findings caution that earlier studies based on analysis of datasets of protein families that were contaminated with epaktologs may have led to some erroneous conclusions about the evolution of novel domain architectures of multidomain proteins. A reassessment of the DA evolution of multidomain proteins is presented in an accompanying paper [1].

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.

Unusually effective microRNA targeting within repeat-rich coding regions of mammalian mRNAs

MicroRNAs (miRNAs) regulate numerous biological processes by base-pairing with target messenger RNAs (mRNAs), primarily through sites in 3' untranslated regions (UTRs), to direct the repression of these targets. Although miRNAs have sometimes been observed to target genes through sites in open reading frames (ORFs), large-scale studies have shown such targeting to be generally less effective than 3' UTR targeting. Here, we show that several miRNAs each target significant groups of genes through multiple sites within their coding regions. This ORF targeting, which mediates both predictable and effective repression, arises from highly repeated sequences containing miRNA target sites. We show that such sequence repeats largely arise through evolutionary duplications and occur particularly frequently within families of paralogous C(2)H(2) zinc-finger genes, suggesting the potential for their coordinated regulation. Examples of ORFs targeted by miR-181 include both the well-known tumor suppressor RB1 and RBAK, encoding a C(2)H(2) zinc-finger protein and transcriptional binding partner of RB1. Our results indicate a function for repeat-rich coding sequences in mediating post-transcriptional regulation and reveal circumstances in which miRNA-mediated repression through ORF sites can be reliably predicted.

Sequence specificity is obtained from the majority of modular C2H2 zinc-finger arrays

C2H2 zinc fingers (C2H2-ZFs) are the most prevalent type of vertebrate DNA-binding domain, and typically appear in tandem arrays (ZFAs), with sequential C2H2-ZFs each contacting three (or more) sequential bases. C2H2-ZFs can be assembled in a modular fashion, providing one explanation for their remarkable evolutionary success. Given a set of modules with defined three-base specificities, modular assembly also presents a way to construct artificial proteins with specific DNA-binding preferences. However, a recent survey of a large number of three-finger ZFAs engineered by modular assembly reported high failure rates (∼70%), casting doubt on the generality of modular assembly. Here, we used protein-binding microarrays to analyze 28 ZFAs that failed in the aforementioned study. Most (17) preferred specific sequences, which in all but one case resembled the intended target sequence. Like natural ZFAs, the engineered ZFAs typically yielded degenerate motifs, binding dozens to hundreds of related individual sequences. Thus, the failure of these proteins in previous assays is not due to lack of sequence-specific DNA-binding activity. Our findings underscore the relevance of individual C2H2-ZF sequence specificities within tandem arrays, and support the general ability of modular assembly to produce ZFAs with sequence-specific DNA-binding activity.

Gains, losses and changes of function after gene duplication: study of the metallothionein family

Metallothioneins (MT) are small proteins involved in heavy metal detoxification and protection against oxidative stress and cancer. The mammalian MT family originated through a series of duplication events which generated four major genes (MT1 to MT4). MT1 and MT2 encode for ubiquitous proteins, while MT3 and MT4 evolved to accomplish specific roles in brain and epithelium, respectively. Herein, phylogenetic, transcriptional and polymorphic analyses are carried out to expose gains, losses and diversification of functions that characterize the evolutionary history of the MT family. The phylogenetic analyses show that all four major genes originated through a single duplication event prior to the radiation of mammals. Further expansion of the MT1 gene has occurred in the primate lineage reaching in humans a total of 13 paralogs, five of which are pseudogenes. In humans, the reading frame of all five MT1 pseudogenes is reconstructed by sequence homology with a functional duplicate revealing that loss of invariant cysteines is the most frequent event accounting for pseudogeneisation. Expression analyses based on EST counts and RT-PCR experiments show that, as for MT1 and MT2, human MT3 is also ubiquitously expressed while MT4 transcripts are present in brain, testes, esophagus and mainly in thymus. Polymorphic variation reveals two deleterious mutations (Cys30Tyr and Arg31Trp) in MT4 with frequencies reaching about 30% in African and Asian populations suggesting the gene is inactive in some individuals and physiological compensation for its loss must arise from a functional equivalent. Altogether our findings provide novel data on the evolution and diversification of MT gene duplicates, a valuable resource for understanding the vast set of biological processes in which these proteins are involved.

Evidence for somatic gene conversion and deletion in bipolar disorder, Crohn's disease, coronary artery disease, hypertension, rheumatoid arthritis, type-1 diabetes, and type-2 diabetes

BACKGROUND

During gene conversion, genetic information is transferred unidirectionally between highly homologous but non-allelic regions of DNA. While germ-line gene conversion has been implicated in the pathogenesis of some diseases, somatic gene conversion has remained technically difficult to investigate on a large scale.

METHODS

A novel analysis technique is proposed for detecting the signature of somatic gene conversion from SNP microarray data. The Wellcome Trust Case Control Consortium has gathered SNP microarray data for two control populations and cohorts for bipolar disorder (BD), cardiovascular disease (CAD), Crohn's disease (CD), hypertension (HT), rheumatoid arthritis (RA), type-1 diabetes (T1D) and type-2 diabetes (T2D). Using the new analysis technique, the seven disease cohorts are analyzed to identify cohort-specific SNPs at which conversion is predicted. The quality of the predictions is assessed by identifying known disease associations for genes in the homologous duplicons, and comparing the frequency of such associations with background rates.

RESULTS

Of 28 disease/locus pairs meeting stringent conditions, 22 show various degrees of disease association, compared with only 8 of 70 in a mock study designed to measure the background association rate (P < 10-9). Additional candidate genes are identified using less stringent filtering conditions. In some cases, somatic deletions appear likely. RA has a distinctive pattern of events relative to other diseases. Similarities in patterns are apparent between BD and HT.

CONCLUSIONS

The associations derived represent the first evidence that somatic gene conversion could be a significant causative factor in each of the seven diseases. The specific genes provide potential insights about disease mechanisms, and are strong candidates for further study.

Dynamic control of endogenous retroviruses during development

Close to half of the human genome encompasses mobile genetic elements, most of which are retrotransposons. These genetic invaders are formidable evolutionary forces that have shaped the architecture of the genomes of higher organisms, with some conserving the ability to induce new integrants within their hosts' genome. Expectedly, the control of endogenous retroviruses is tight and multi-pronged. It is most crucially established in the germ line and during the first steps of embryogenesis, primarily through transcriptional mechanisms that have likely evolved under their very pressure, but are now engaged in controlling gene expression at large, notably during early development.

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.

Gene inactivation and its implications for annotation in the era of personal genomics

The first wave of personal genomes documents how no single individual genome contains the full complement of functional genes. Here, we describe the extent of variation in gene and pseudogene numbers between individuals arising from inactivation events such as premature termination or aberrant splicing due to single-nucleotide polymorphisms. This highlights the inadequacy of the current reference sequence and gene set. We present a proposal to define a reference gene set that will remain stable as more individuals are sequenced. In particular, we recommend that the ancestral allele be used to define the reference sequence from which a core human reference gene annotation set can be derived. In addition, we call for the development of an expanded gene set to include human-specific genes that have arisen recently and are absent from the ancestral set.

ZNF274 recruits the histone methyltransferase SETDB1 to the 3' ends of ZNF genes

Only a small percentage of human transcription factors (e.g. those associated with a specific differentiation program) are expressed in a given cell type. Thus, cell fate is mainly determined by cell type-specific silencing of transcription factors that drive different cellular lineages. Several histone modifications have been associated with gene silencing, including H3K27me3 and H3K9me3. We have previously shown that genes for the two largest classes of mammalian transcription factors are marked by distinct histone modifications; homeobox genes are marked by H3K27me3 and zinc finger genes are marked by H3K9me3. Several histone methyltransferases (e.g. G9a and SETDB1) may be involved in mediating the H3K9me3 silencing mark. We have used ChIP-chip and ChIP-seq to demonstrate that SETDB1, but not G9a, is associated with regions of the genome enriched for H3K9me3. One current model is that SETDB1 is recruited to specific genomic locations via interaction with the corepressor TRIM28 (KAP1), which is in turn recruited to the genome via interaction with zinc finger transcription factors that contain a Kruppel-associated box (KRAB) domain. However, specific KRAB-ZNFs that recruit TRIM28 (KAP1) and SETDB1 to the genome have not been identified. We now show that ZNF274 (a KRAB-ZNF that contains 5 C2H2 zinc finger domains), can interact with KAP1 both in vivo and in vitro and, using ChIP-seq, we show that ZNF274 binding sites co-localize with SETDB1, KAP1, and H3K9me3 at the 3′ ends of zinc finger genes. Knockdown of ZNF274 with siRNAs reduced the levels of KAP1 and SETDB1 recruitment to the binding sites. These studies provide the first identification of a KRAB domain-containing ZNF that is involved in recruitment of the KAP1 and SETDB1 to specific regions of the human genome.

Rapid sequence and expression divergence suggest selection for novel function in primate-specific KRAB-ZNF genes

Recent segmental duplications (SDs), arising from duplication events that occurred within the past 35-40 My, have provided a major resource for the evolution of proteins with primate-specific functions. KRAB zinc finger (KRAB-ZNF) transcription factor genes are overrepresented among genes contained within these recent human SDs. Here, we examine the structural and functional diversity of the 70 human KRAB-ZNF genes involved in the most recent primate SD events including genes that arose in the hominid lineage. Despite their recent advent, many parent-daughter KRAB-ZNF gene pairs display significant differences in zinc finger structure and sequence, expression, and splicing patterns, each of which could significantly alter the regulatory functions of the paralogous genes. Paralogs that emerged on the lineage to humans and chimpanzees have undergone more evolutionary changes per unit of time than genes already present in the common ancestor of rhesus macaques and great apes. Taken together, these data indicate that a substantial fraction of the recently evolved primate-specific KRAB-ZNF gene duplicates have acquired novel functions that may possibly define novel regulatory pathways and suggest an active ongoing selection for regulatory diversity in primates.

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.

KRAB-zinc finger proteins and KAP1 can mediate long-range transcriptional repression through heterochromatin spreading

Krüppel-associated box domain-zinc finger proteins (KRAB–ZFPs) are tetrapod-specific transcriptional repressors encoded in the hundreds by the human genome. In order to explore their as yet ill-defined impact on gene expression, we developed an ectopic repressor assay, allowing the study of KRAB–mediated transcriptional regulation at hundreds of different transcriptional units. By targeting a drug-controllable KRAB–containing repressor to gene-trapping lentiviral vectors, we demonstrate that KRAB and its corepressor KAP1 can silence promoters located several tens of kilobases (kb) away from their DNA binding sites, with an efficiency which is generally higher for promoters located within 15 kb or less. Silenced promoters exhibit a loss of histone H3-acetylation, an increase in H3 lysine 9 trimethylation (H3K9me3), and a drop in RNA Pol II recruitment, consistent with a block of transcriptional initiation following the establishment of silencing marks. Furthermore, we reveal that KRAB–mediated repression is established by the long-range spreading of H3K9me3 and heterochromatin protein 1 β (HP1β) between the repressor binding site and the promoter. We confirm the biological relevance of this phenomenon by documenting KAP1–dependent transcriptional repression at an endogenous KRAB–ZFP gene cluster, where KAP1 binds to the 3′ end of genes and mediates propagation of H3K9me3 and HP1β towards their 5′ end. Together, our data support a model in which KRAB/KAP1 recruitment induces long-range repression through the spread of heterochromatin. This finding not only suggests auto-regulatory mechanisms in the control of KRAB–ZFP gene clusters, but also provides important cues for interpreting future genome-wide DNA binding data of KRAB–ZFPs and KAP1.

The regulation of gene activity by transcription factors is crucial to the function of all cells. Here, we studied the mechanisms of action of the largest family of gene regulators encoded by the human genome, the so-called KRAB–containing zinc finger proteins (KRAB–ZFPs), which in concert with their universal cofactor KAP1 act as transcriptional repressors. For this, we used two parallel approaches. First, by targeting an ectopic KRAB domain to hundreds of different genes, we found that KRAB/KAP1 can repress promoters located several tens of kilobases from the repressor DNA docking site. We further could show that KRAB induces such long-range effects by mediating the spread of repressive chromatin marks along the body of the gene, resulting in a block of transcriptional initiation at the promoter. In a second set of experiments, we analyzed an endogenous KRAB–ZFP gene cluster, where we could also document KAP1–dependent heterochromatin spreading and transcriptional repression. Together, these results support a model whereby KRAB–ZFPs and KAP1 can mediate long-range transcriptional repression through the spread of silencing chromatin marks. This study thus provides insight into KRAB/KAP1–induced gene regulation at KRAB–ZFP gene clusters, and will further help interpret genome-wide studies of KRAB–ZFPs and KAP1 DNA binding patterns.

Phylogenomics reveals subfamilies of fungal nonribosomal peptide synthetases and their evolutionary relationships

BACKGROUND

Nonribosomal peptide synthetases (NRPSs) are multimodular enzymes, found in fungi and bacteria, which biosynthesize peptides without the aid of ribosomes. Although their metabolite products have been the subject of intense investigation due to their life-saving roles as medicinals and injurious roles as mycotoxins and virulence factors, little is known of the phylogenetic relationships of the corresponding NRPSs or whether they can be ranked into subgroups of common function. We identified genes (NPS) encoding NRPS and NRPS-like proteins in 38 fungal genomes and undertook phylogenomic analyses in order to identify fungal NRPS subfamilies, assess taxonomic distribution, evaluate levels of conservation across subfamilies, and address mechanisms of evolution of multimodular NRPSs. We also characterized relationships of fungal NRPSs, a representative sampling of bacterial NRPSs, and related adenylating enzymes, including alpha-aminoadipate reductases (AARs) involved in lysine biosynthesis in fungi.

RESULTS

Phylogenomic analysis identified nine major subfamilies of fungal NRPSs which fell into two main groups: one corresponds to NPS genes encoding primarily mono/bi-modular enzymes which grouped with bacterial NRPSs and the other includes genes encoding primarily multimodular and exclusively fungal NRPSs. AARs shared a closer phylogenetic relationship to NRPSs than to other acyl-adenylating enzymes. Phylogenetic analyses and taxonomic distribution suggest that several mono/bi-modular subfamilies arose either prior to, or early in, the evolution of fungi, while two multimodular groups appear restricted to and expanded in fungi. The older mono/bi-modular subfamilies show conserved domain architectures suggestive of functional conservation, while multimodular NRPSs, particularly those unique to euascomycetes, show a diversity of architectures and of genetic mechanisms generating this diversity.

CONCLUSIONS

This work is the first to characterize subfamilies of fungal NRPSs. Our analyses suggest that mono/bi-modular NRPSs have more ancient origins and more conserved domain architectures than most multimodular NRPSs. It also demonstrates that the alpha-aminoadipate reductases involved in lysine biosynthesis in fungi are closely related to mono/bi-modular NRPSs. Several groups of mono/bi-modular NRPS metabolites are predicted to play more pivotal roles in cellular metabolism than products of multimodular NRPSs. In contrast, multimodular subfamilies of NRPSs are of more recent origin, are restricted to fungi, show less stable domain architectures, and biosynthesize metabolites which perform more niche-specific functions than mono/bi-modular NRPS products. The euascomycete-only NRPS subfamily, in particular, shows evidence for extensive gain and loss of domains suggestive of the contribution of domain duplication and loss in responding to niche-specific pressures.

Lineage-specific transcription factors and the evolution of gene regulatory networks

Nature is replete with examples of diverse cell types, tissues and body plans, forming very different creatures from genomes with similar gene complements. However, while the genes and the structures of proteins they encode can be highly conserved, the production of those proteins in specific cell types and at specific developmental time points might differ considerably between species. A full understanding of the factors that orchestrate gene expression will be essential to fully understand evolutionary variety. Transcription factor (TF) proteins, which form gene regulatory networks (GRNs) to act in cooperative or competitive partnerships to regulate gene expression, are key components of these unique regulatory programs. Although many TFs are conserved in structure and function, certain classes of TFs display extensive levels of species diversity. In this review, we highlight families of TFs that have expanded through gene duplication events to create species-unique repertoires in different evolutionary lineages. We discuss how the hierarchical structures of GRNs allow for flexible small to large-scale phenotypic changes. We survey evidence that explains how newly evolved TFs may be integrated into an existing GRN and how molecular changes in TFs might impact the GRNs. Finally, we review examples of traits that evolved due to lineage-specific TFs and species differences in GRNs.

Accelerated evolution of the Prdm9 speciation gene across diverse metazoan taxa

The onset of prezygotic and postzygotic barriers to gene flow between populations is a hallmark of speciation. One of the earliest postzygotic isolating barriers to arise between incipient species is the sterility of the heterogametic sex in interspecies' hybrids. Four genes that underlie hybrid sterility have been identified in animals: Odysseus, JYalpha, and Overdrive in Drosophila and Prdm9 (Meisetz) in mice. Mouse Prdm9 encodes a protein with a KRAB motif, a histone methyltransferase domain and several zinc fingers. The difference of a single zinc finger distinguishes Prdm9 alleles that cause hybrid sterility from those that do not. We find that concerted evolution and positive selection have rapidly altered the number and sequence of Prdm9 zinc fingers across 13 rodent genomes. The patterns of positive selection in Prdm9 zinc fingers imply that rapid evolution has acted on the interface between the Prdm9 protein and the DNA sequences to which it binds. Similar patterns are apparent for Prdm9 zinc fingers for diverse metazoans, including primates. Indeed, allelic variation at the DNA–binding positions of human PRDM9 zinc fingers show significant association with decreased risk of infertility. Prdm9 thus plays a role in determining male sterility both between species (mouse) and within species (human). The recurrent episodes of positive selection acting on Prdm9 suggest that the DNA sequences to which it binds must also be evolving rapidly. Our findings do not identify the nature of the underlying DNA sequences, but argue against the proposed role of Prdm9 as an essential transcription factor in mouse meiosis. We propose a hypothetical model in which incompatibilities between Prdm9-binding specificity and satellite DNAs provide the molecular basis for Prdm9-mediated hybrid sterility. We suggest that Prdm9 should be investigated as a candidate gene in other instances of hybrid sterility in metazoans.

Speciation, the process by which one species splits into two, involves reproductive barriers between previously interbreeding populations. The question of how speciation occurs has rightly occupied the attention of biologists since before Darwin's “On the Origin of Species.” Studies of recently diverged species have revealed the presence of hybrid sterility genes (colloquially referred to as “speciation genes”), alleles of which are associated with sterility of interspecies hybrids. Mouse Prdm9 is the only known such gene in vertebrate animals. Here we report that the Prdm9 protein has evolved extremely rapidly in its DNA-binding domain, comprising an array of “zinc fingers.” This suggests that hybrid sterility may arise from a mismatch between the DNA-binding specificity of Prdm9 and rapidly evolving DNA. We propose that Prdm9 binds to satellite-DNA repeats evolving rapidly within and between different species. Prdm9 evolution is unusual because other hybrid sterility genes appear only to evolve rapidly in isolated bursts, whereas Prdm9 has evolved rapidly over 700 million years, in many rodent species, diverse primates and other metazoans. This leads to the tantalizing possibility that Prdm9 may have served as a “speciation gene” on other occasions in metazoan evolution, a possibility that will now need to be investigated.

Gene expression signatures and response to imatinib mesylate in gastrointestinal stromal tumor

Despite initial efficacy of imatinib mesylate in most gastrointestinal stromal tumor (GIST) patients, many experience primary/secondary drug resistance. Therefore, clinical management of GIST may benefit from further molecular characterization of tumors before and after imatinib mesylate treatment. As part of a recent phase II trial of neoadjuvant/adjuvant imatinib mesylate treatment for advanced primary and recurrent operable GISTs (Radiation Therapy Oncology Group S0132), gene expression profiling using oligonucleotide microarrays was done on tumor samples obtained before and after imatinib mesylate therapy. Patients were classified according to changes in tumor size after treatment based on computed tomography scan measurements. Gene profiling data were evaluated with Statistical Analysis of Microarrays to identify differentially expressed genes (in pretreatment GIST samples). Based on Statistical Analysis of Microarrays [False Discovery Rate (FDR), 10%], 38 genes were expressed at significantly lower levels in the pretreatment biopsy samples from tumors that significantly responded to 8 to 12 weeks of imatinib mesylate, that is, >25% tumor reduction. Eighteen of these genes encoded Krüppel-associated box (KRAB) domain containing zinc finger (ZNF) transcriptional repressors. Importantly, 10 KRAB-ZNF genes mapped to a single locus on chromosome 19p, and a subset predicted likely response to imatinib mesylate-based therapy in a naïve panel of GIST. Furthermore, we found that modifying expression of genes within this predictive signature can enhance the sensitivity of GIST cells to imatinib mesylate. Using clinical pretreatment biopsy samples from a prospective neoadjuvant phase II trial, we have identified a gene signature that includes KRAB-ZNF 91 subfamily members that may be both predictive of and functionally associated with likely response to short-term imatinib mesylate treatment.

Genomic evidence for non-random endemic populations of decaying exons from mammalian genes

Background

Functional diversification of genes in mammalian genomes is engendered by a number of processes, e.g., gene duplication and alternative splicing. Gene duplication is classically discussed as leading to neofunctionalization (generation of new functions), subfunctionalization (generation of a varied function), or pseudogenization (loss of the gene and its function).

Results

Here, we focus on the process of pseudogenization, but specifically for individual exons from genes. It is at present unclear to what extent pseudogenization of individual exon duplications affects gene evolution, i.e., is it a random phenomenon, or is it associated with specific types of genes and encoded proteins, and positions in gene structures? We gathered genomic evidence for pseudogenic exons (ΨEs, i.e., exons disabled by frameshifts and premature stop codons), to examine for significant trends in their distribution across four mammalian genomes (specifically human, cow, mouse and rat). Across these four genomes, we observed a consistent population of ΨEs, associated with 0.4–1.0% of genes. These ΨE populations exhibit codon substitution patterns that are typical of an endemic population of decaying sequences. In human, ΨEs have significant over-representation for functional categories related to 'ion binding' and 'nucleic-acid binding', compared to duplicated exons in general. Also, ΨEs tend to be associated with some protein domains that are abundant generally, e.g., Zinc-finger and immunoglobulin protein domains, but not others, e.g., EGF-like domains. Positionally, ΨEs are also significantly associated with the 5' end of genes, but despite this, individual stop codons are positioned so that there is significant avoidance of potential targeting to nonsense-mediated decay. In human, ΨEs are often associated with alternative splicing (in 22 out of 284 genes with ΨEs in their milieu), and can have different parts of their sequence differentially spliced in alternative transcripts. Some unusual cases of ΨEs embedded within 5' and 3' non-coding exons are observed.

Conclusion

Our results indicate the types of genes that harbour ΨEs, and demonstrate that ΨEs have non-random distribution within gene structures. These ΨEs may function in gene regulation through generation of transcribed pseudogenes, or regulatory alternate transcripts.

TFCat: the curated catalog of mouse and human transcription factors

TFCat is a catalog of mouse and human transcription factors based on a reliable core collection of annotations obtained by expert review of the scientific literature

Unravelling regulatory programs governed by transcription factors (TFs) is fundamental to understanding biological systems. TFCat is a catalog of mouse and human TFs based on a reliable core collection of annotations obtained by expert review of the scientific literature. The collection, including proven and homology-based candidate TFs, is annotated within a function-based taxonomy and DNA-binding proteins are organized within a classification system. All data and user-feedback mechanisms are available at the TFCat portal .