Adaptive evolution in zinc finger transcription factors

Suboptimal Intermediates Underlie Evolution of the Bicoid Homeodomain

Changes in regulatory networks generate materials for evolution to create phenotypic diversity. For transcription networks, multiple studies have shown that alterations in binding sites of cis-regulatory elements correlate well with the gain or loss of specific features of the body plan. Less is known about alterations in the amino acid sequences of the transcription factors (TFs) that bind these elements. Here we study the evolution of Bicoid (Bcd), a homeodomain (HD) protein that is critical for anterior embryo patterning in Drosophila. The ancestor of Bcd (AncBcd) emerged after a duplication of a Zerknullt (Zen)-like ancestral protein (AncZB) in a suborder of flies. AncBcd diverged from AncZB, gaining novel transcriptional and translational activities. We focus on the evolution of the HD of AncBcd, which binds to DNA and RNA, and is comprised of four subdomains: an N-terminal arm (NT) and three helices; H1, H2, and Recognition Helix (RH). Using chimeras of subdomains and gene rescue assays in Drosophila, we show that robust patterning activity of the Bcd HD (high frequency rescue to adulthood) is achieved only when amino acid substitutions in three separate subdomains (NT, H1, and RH) are combined. Other combinations of subdomains also yield full rescue, but with lower penetrance, suggesting alternative suboptimal activities. Our results suggest a multistep pathway for the evolution of the Bcd HD that involved intermediate HD sequences with suboptimal activities, which constrained and enabled further evolutionary changes. They also demonstrate critical epistatic forces that contribute to the robust function of a DNA-binding domain.

TFs for TEs: the transcription factor repertoire of mammalian transposable elements

In this review, Hermant and Torres-Padilla summarize and discuss the transcription factors known to be involved in the sequence-specific recognition and transcriptional activation of specific transposable element families or subfamilies.

Transposable elements (TEs) are genetic elements capable of changing position within the genome. Although their mobilization can constitute a threat to genome integrity, nearly half of modern mammalian genomes are composed of remnants of TE insertions. The first critical step for a successful transposition cycle is the generation of a full-length transcript. TEs have evolved cis-regulatory elements enabling them to recruit host-encoded factors driving their own, selfish transcription. TEs are generally transcriptionally silenced in somatic cells, and the mechanisms underlying their repression have been extensively studied. However, during germline formation, preimplantation development, and tumorigenesis, specific TE families are highly expressed. Understanding the molecular players at stake in these contexts is of utmost importance to establish the mechanisms regulating TEs, as well as the importance of their transcription to the biology of the host. Here, we review the transcription factors known to be involved in the sequence-specific recognition and transcriptional activation of specific TE families or subfamilies. We discuss the diversity of TE regulatory elements within mammalian genomes and highlight the importance of TE mobilization in the dispersal of transcription factor-binding sites over the course of evolution.

Mice Deficient in the RNA-Binding Protein Zfp871 Are Prone to Early Death and Steatohepatitis in Part through the p53-Mdm2 Axis

p53 transcription factor is activated upon exposure to various cellular stresses, leading to growth suppression. However, aberrant activation of p53 can lead to defects in embryonic development and other abnormalities. Here, we identified zinc finger protein Zfp871 as a p53 target gene. We showed that as an RNA-binding protein, Zfp871 binds to Mdm2 3'UTR and stabilizes Mdm2 mRNA, which in turn suppresses p53 expression through increased expression of Mdm2 E3 ubiquitin ligase. Consistently, Zfp871 deficiency increases p53 expression, leading to growth suppression in a p53-dependent manner. To determine the role of Zfp871 in the p53 pathway, we used Zfp871-deficient mouse model and found that Zfp871-null mice were prone to embryonic/pre-weaning lethality, which can be partially rescued by simultaneous deletion of Trp53. We also found that mice heterozygous for Zfp871 had a short lifespan and were susceptible to steatohepatitis but not to spontaneous tumors. To determine the underlying mechanism, RNA-seq analysis was performed and showed that an array of genes involved in development, lipid metabolism, and inflammation is regulated by Zfp871 in conjunction with p53. Taken together, we conclude that the Zfp871-Mdm2-p53 pathway plays a critical role in tumor-free survival and development. IMPLICATIONS: A fine equilibrium of p53 is required for preventing damaging effects of aberrant p53 expression. We identify the Zfp871-Mdm2-p53 pathway that plays a critical role in development of mice and steatohepatitis.

Inflammasome, the Constitutive Heterochromatin Machinery, and Replication of an Oncogenic Herpesvirus

The success of long-term host–virus partnerships is predicated on the ability of the host to limit the destructive potential of the virus and the virus’s skill in manipulating its host to persist undetected yet replicate efficiently when needed. By mastering such skills, herpesviruses persist silently in their hosts, though perturbations in this host–virus equilibrium can result in disease. The heterochromatin machinery that tightly regulates endogenous retroviral elements and pericentromeric repeats also silences invading genomes of alpha-, beta-, and gammaherpesviruses. That said, how these viruses disrupt this constitutive heterochromatin machinery to replicate and spread, particularly in response to disparate lytic triggers, is unclear. Here, we review how the cancer-causing gammaherpesvirus Epstein–Barr virus (EBV) uses the inflammasome as a security system to alert itself of threats to its cellular home as well as to flip the virus-encoded lytic switch, allowing it to replicate and escape in response to a variety of lytic triggers. EBV provides the first example of an infectious agent able to actively exploit the inflammasome to spark its replication. Revealing an unexpected link between the inflammasome and the epigenome, this further brings insights into how the heterochromatin machinery uses differential strategies to maintain the integrity of the cellular genome whilst guarding against invading pathogens. These recent insights into EBV biology and host–viral epigenetic regulation ultimately point to the NLRP3 inflammasome as an attractive target to thwart herpesvirus reactivation.

Conservation of Aging and Cancer Epigenetic Signatures across Human and Mouse

Aging and cancer are two interrelated processes, with aging being a major risk factor for the development of cancer. Parallel epigenetic alterations have been described for both, although differences, especially within the DNA hypomethylation scenario, have also been recently reported. Although many of these observations arise from the use of mouse models, there is a lack of systematic comparisons of human and mouse epigenetic patterns in the context of disease. However, such comparisons are significant as they allow to establish the extent to which some of the observed similarities or differences arise from pre-existing species-specific epigenetic traits. Here, we have used reduced representation bisulfite sequencing to profile the brain methylomes of young and old, tumoral and nontumoral brain samples from human and mouse. We first characterized the baseline epigenomic patterns of the species and subsequently focused on the DNA methylation alterations associated with cancer and aging. Next, we described the functional genomic and epigenomic context associated with the alterations, and finally, we integrated our data to study interspecies DNA methylation levels at orthologous CpG sites. Globally, we found considerable differences between the characteristics of DNA methylation alterations in cancer and aging in both species. Moreover, we describe robust evidence for the conservation of the specific cancer and aging epigenomic signatures in human and mouse. Our observations point toward the preservation of the functional consequences of these alterations at multiple levels of genomic regulation. Finally, our analyses reveal a role for the genomic context in explaining disease- and species-specific epigenetic traits.

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.

Innovation of heterochromatin functions drives rapid evolution of essential ZAD-ZNF genes in Drosophila

Contrary to dogma, evolutionarily young and dynamic genes can encode essential functions. We find that evolutionarily dynamic ZAD-ZNF genes, which encode the most abundant class of insect transcription factors, are more likely to encode essential functions in Drosophila melanogaster than ancient, conserved ZAD-ZNF genes. We focus on the Nicknack ZAD-ZNF gene, which is evolutionarily young, poorly retained in Drosophila species, and evolves under strong positive selection. Yet we find that it is necessary for larval development in D. melanogaster. We show that Nicknack encodes a heterochromatin-localizing protein like its paralog Oddjob, also an evolutionarily dynamic yet essential ZAD-ZNF gene. We find that the divergent D. simulans Nicknack protein can still localize to D. melanogaster heterochromatin and rescue viability of female but not male Nicknack-null D. melanogaster. Our findings suggest that innovation for rapidly changing heterochromatin functions might generally explain the essentiality of many evolutionarily dynamic ZAD-ZNF genes in insects.

Effects of Naturally Occurring Mutations in Bovine Leukemia Virus 5'-LTR and Tax Gene on Viral Transcriptional Activity

Bovine leukemia virus (BLV) is a deltaretrovirus infecting bovine B cells and causing enzootic bovine leucosis (EBL). The long terminal repeat (LTR) plays an indispensable role in viral gene expression. The BLV Tax protein acts as the main transactivator of LTR-driven transcription of BLV viral genes. The aim of this study was to analyze mutations in the BLV LTR region and tax gene to determine their association with transcriptional activity. LTRs were obtained from one hundred and six BLV isolates and analyzed for their genetic variability. Fifteen variants were selected and characterized based on mutations in LTR regulatory elements, and further used for in vitro transcription assays. Reporter vectors containing the luciferase gene under the control of each variant BLV promoter sequence, in addition to variant Tax expression vectors, were constructed. Both types of plasmids were used for cotransfection of HeLa cells and the level of luciferase activity was measured as a proxy of transcriptional activity. Marked differences in LTR promoter activity and Tax transactivation activity were observed amongst BLV variants. These results demonstrate that mutations in both the BLV LTR and tax gene can affect the promoter activity, which may have important consequences on proviral load, viral fitness, and transmissibility in BLV-infected cattle.

Characterization and evolutionary origin of novel C2H2 zinc finger protein (ZNF648) required for both erythroid and megakaryocyte differentiation in humans

Human ZNF648 is a novel poly C-terminal C2H2 zinc finger protein identified amongst the most dysregulated proteins in erythroid cells differentiated from iPSC. Its nuclear localisation and structure indicate it is likely a DNA-binding protein. Using a combination of ZNF648 overexpression in an iPSC line and primary adult erythroid cells, ZNF648 knockdown in primary adult erythroid cells and megakaryocytes, comparative proteomics and transcriptomics we show that ZNF648 is required for both erythroid and megakaryocyte differentiation. Orthologues of ZNF648 were detected across Mammals, Reptilia, Actinopterygii, in some Aves, Amphibia and Coelacanthiformes suggesting the gene originated in the common ancestor of Osteichthyes (Euteleostomi or bony fish). Conservation of the C-terminal zinc finger domain is higher, with some variation in zinc finger number but a core of at least six zinc fingers conserved across all groups, with the N-terminus recognisably similar within but not between major lineages. This suggests the N-terminus of ZNF648 evolves faster than the C-terminus, however this is not due to exon-shuffling as the entire coding region of ZNF648 is within a single exon. As for other such transcription factors, the N-terminus likely carries out regulatory functions, but showed no sequence similarity to any known domains. The greater functional constraint on the zinc finger domain suggests ZNF648 binds at least some similar regions of DNA in the different organisms. However, divergence of the N-terminal region may enable differential expression, allowing adaptation of function in the different organisms.

Regulation of Expression and Latency in BLV and HTLV

Human T-lymphotrophic virus type 1 (HTLV-1) and Bovine leukemia virus (BLV) belong to the Deltaretrovirus genus. HTLV-1 is the etiologic agent of the highly aggressive and currently incurable cancer adult T-cell leukemia (ATL) and a neurological disease HTLV-1-associated myelopathy (HAM)/tropical spastic paraparesis (TSP). BLV causes neoplastic proliferation of B cells in cattle: enzootic bovine leucosis (EBL). Despite the severity of these conditions, infection by HTLV-1 and BLV appear in most cases clinically asymptomatic. These viruses can undergo latency in their hosts. The silencing of proviral gene expression and maintenance of latency are central for the establishment of persistent infection, as well as for pathogenesis in vivo. In this review, we will present the mechanisms that control proviral activation and retroviral latency in deltaretroviruses, in comparison with other exogenous retroviruses. The 5′ long terminal repeats (5′-LTRs) play a main role in controlling viral gene expression. While the regulation of transcription initiation is a major mechanism of silencing, we discuss topics that include (i) the epigenetic control of the provirus, (ii) the cis-elements present in the LTR, (iii) enhancers with cell-type specific regulatory functions, (iv) the role of virally-encoded transactivator proteins, (v) the role of repressors in transcription and silencing, (vi) the effect of hormonal signaling, (vii) implications of LTR variability on transcription and latency, and (viii) the regulatory role of non-coding RNAs. Finally, we discuss how a better understanding of these mechanisms may allow for the development of more effective treatments against Deltaretroviruses.

Strain-Specific Epigenetic Regulation of Endogenous Retroviruses: The Role of Trans-Acting Modifiers

Approximately 10 percent of the mouse genome consists of endogenous retroviruses (ERVs), relics of ancient retroviral infections that are classified based on their relatedness to exogenous retroviral genera. Because of the ability of ERVs to retrotranspose, as well as their cis-acting regulatory potential due to functional elements located within the elements, mammalian ERVs are generally subject to epigenetic silencing by DNA methylation and repressive histone modifications. The mobilisation and expansion of ERV elements is strain-specific, leading to ERVs being highly polymorphic between inbred mouse strains, hinting at the possibility of the strain-specific regulation of ERVs. In this review, we describe the existing evidence of mouse strain-specific epigenetic control of ERVs and discuss the implications of differential ERV regulation on epigenetic inheritance models. We consider Krüppel-associated box domain (KRAB) zinc finger proteins as likely candidates for strain-specific ERV modifiers, drawing on insights gained from the study of the strain-specific behaviour of transgenes. We conclude by considering the coevolution of KRAB zinc finger proteins and actively transposing ERV elements, and highlight the importance of cross-strain studies in elucidating the mechanisms and consequences of strain-specific ERV regulation.

The genetic factors of bilaterian evolution

The Cambrian explosion was a unique animal radiation ~540 million years ago that produced the full range of body plans across bilaterians. The genetic mechanisms underlying these events are unknown, leaving a fundamental question in evolutionary biology unanswered. Using large-scale comparative genomics and advanced orthology evaluation techniques, we identified 157 bilaterian-specific genes. They include the entire Nodal pathway, a key regulator of mesoderm development and left-right axis specification; components for nervous system development, including a suite of G-protein-coupled receptors that control physiology and behaviour, the Robo-Slit midline repulsion system, and the neurotrophin signalling system; a high number of zinc finger transcription factors; and novel factors that previously escaped attention. Contradicting the current view, our study reveals that genes with bilaterian origin are robustly associated with key features in extant bilaterians, suggesting a causal relationship.

Independent Whole-Genome Duplications Define the Architecture of the Genomes of the Devastating West African Cacao Black Pod Pathogen Phytophthora megakarya and Its Close Relative Phytophthora palmivora

Phytophthora megakarya and P. palmivora are oomycete pathogens that cause black pod rot of cacao (Theobroma cacao), the most economically important disease on cacao globally. While P. palmivora is a cosmopolitan pathogen, P. megakarya, which is more aggressive on cacao than P. palmivora, has been reported only in West and Central Africa where it has been spreading and devastating cacao farms since the 1950s. In this study, we reconstructed the complete diploid genomes of multiple isolates of both species using single-molecule real-time sequencing. Thirty-one additional genotypes were sequenced to analyze inter- and intra-species genomic diversity. The P. megakarya genome is exceptionally large (222 Mbp) and nearly twice the size of P. palmivora (135 Mbp) and most known Phytophthora species (∼100 Mbp on average). Previous reports pointed toward a whole-genome duplication (WGD) in P. palmivora In this study, we demonstrate that both species underwent independent and relatively recent WGD events. In P. megakarya we identified a unique combination of WGD and large-scale transposable element driven genome expansion, which places this genome in the upper range of Phytophthora genome sizes, as well as effector pools with 1,382 predicted RxLR effectors. Finally, this study provides evidence of adaptive evolution of effectors like RxLRs and Crinklers, and discusses the implications of effector expansion and diversification.

The geometric influence on the Cys2His2 zinc finger domain and functional plasticity

The Cys₂His₂ zinc finger is the most common DNA-binding domain expanding in metazoans since the fungi human split. A proposed catalyst for this expansion is an arms race to silence transposable elements yet it remains poorly understood how this domain is able to evolve the required specificities. Likewise, models of its DNA binding specificity remain error prone due to a lack of understanding of how adjacent fingers influence each other's binding specificity. Here, we use a synthetic approach to exhaustively investigate binding geometry, one of the dominant influences on adjacent finger function. By screening over 28 billion protein–DNA interactions in various geometric contexts we find the plasticity of the most common natural geometry enables more functional amino acid combinations across all targets. Further, residues that define this geometry are enriched in genomes where zinc fingers are prevalent and specificity transitions would be limited in alternative geometries. Finally, these results demonstrate an exhaustive synthetic screen can produce an accurate model of domain function while providing mechanistic insight that may have assisted in the domains expansion.

Comparative Transcriptomics Reveals Gene Families Associated with Predatory Behavior in Photuris femme fatale Fireflies

Identifying the basis of phenotypic variation is a key objective of genetics. This work has been mostly limited to model systems with a plethora of genetic manipulation and functional characterization tools. With the development of high-throughput sequencing and new computational tools, it is possible to identify candidate genes related to phenotypic variation in non-model organisms. Fireflies are excellent for studying phenotypic variation because of their diverse and well-characterized behaviors. Most adult fireflies emit a single mating flash pattern and do not eat. In contrast, adult females of many species in the genus Photuris employ multiple flash patterns and prey upon mate-seeking males of other firefly species. To investigate the genetic basis for this variation, we used comparative transcriptomics to identify positively selected genes between a predatory firefly, Photuris sp., and a non-predatory relative, Photuris frontalis, controlling for genes generally under selection in fireflies by comparing to a Photinus firefly. Nine gene families were identified under positive selection in the predatory versus non-predatory Photuris comparison, including genes involved in digestion, detoxification, vision, reproduction, and neural processes. These results generate intriguing hypotheses about the genetic basis for insect behavior and highlight the utility of comparative transcriptomic tools to investigate complex behaviors in non-model systems.

Structural genome analysis in cultivated potato taxa

KEY MESSAGE

Twelve potato accessions were selected to represent two principal views on potato taxonomy. The genomes were sequenced and analyzed for structural variation (copy number variation) against three published potato genomes. The common potato (Solanum tuberosum L.) is an important staple crop with a highly heterozygous and complex tetraploid genome. The other taxa of cultivated potato contain varying ploidy levels (2X-5X), and structural variations are common in the genomes of these species, likely contributing to the diversification or agronomic traits during domestication. Increased understanding of the genomes and genomic variation will aid in the exploration of novel agronomic traits. Thus, sequencing data from twelve potato landraces, representing the four ploidy levels, were used to identify structural genomic variation compared to the two currently available reference genomes, a double monoploid potato genome and a diploid inbred clone of S. chacoense. The results of a copy number variation analysis showed that in the majority of the genomes, while the number of deletions is greater than the number of duplications, the number of duplicated genes is greater than the number of deleted ones. Specific regions in the twelve potato genomes have a high density of CNV events. Further, the auxin-induced SAUR genes (involved in abiotic stress), disease resistance genes and the 2-oxoglutarate/Fe(II)-dependent oxygenase superfamily proteins, among others, had increased copy numbers in these sequenced genomes relative to the references.

C2H2-Type Zinc Finger Proteins in Brain Development, Neurodevelopmental, and Other Neuropsychiatric Disorders: Systematic Literature-Based Analysis

Neurodevelopmental disorders (NDDs) are multifaceted pathologic conditions manifested with intellectual disability, autistic features, psychiatric problems, motor dysfunction, and/or genetic/chromosomal abnormalities. They are associated with skewed neurogenesis and brain development, in part through dysfunction of the neural stem cells (NSCs) where abnormal transcriptional regulation on key genes play significant roles. Recent accumulated evidence highlights C₂H₂-type zinc finger proteins (C₂H₂-ZNFs), the largest transcription factor family in humans, as important targets for the pathologic processes associated with NDDs. In this review, we identified their significant accumulation (74 C₂H₂-ZNFs: ~10% of all human member proteins) in brain physiology and pathology. Specifically, we discuss their physiologic contribution to brain development, particularly focusing on their actions in NSCs. We then explain their pathologic implications in various forms of NDDs, such as morphological brain abnormalities, intellectual disabilities, and psychiatric disorders. We found an important tendency that poly-ZNFs and KRAB-ZNFs tend to be involved in the diseases that compromise gross brain structure and human-specific higher-order functions, respectively. This may be consistent with their characteristic appearance in the course of species evolution and corresponding contribution to these brain activities.

Ovarian Transcriptomic Analyses in the Urban Human Health Pest, the Western Black Widow Spider

Due to their abundance and ability to invade diverse environments, many arthropods have become pests of economic and health concern, especially in urban areas. Transcriptomic analyses of arthropod ovaries have provided insight into life history variation and fecundity, yet there are few studies in spiders despite their diversity within arthropods. Here, we generated a de novo ovarian transcriptome from 10 individuals of the western black widow spider (Latrodectus hesperus), a human health pest of high abundance in urban areas, to conduct comparative ovarian transcriptomic analyses. Biological processes enriched for metabolism—specifically purine, and thiamine metabolic pathways linked to oocyte development—were significantly abundant in L. hesperus. Functional and pathway annotations revealed overlap among diverse arachnid ovarian transcriptomes for highly-conserved genes and those linked to fecundity, such as oocyte maturation in vitellogenin and vitelline membrane outer layer proteins, hormones, and hormone receptors required for ovary development, and regulation of fertility-related genes. Comparative studies across arachnids are greatly needed to understand the evolutionary similarities of the spider ovary, and here, the identification of ovarian proteins in L. hesperus provides potential for understanding how increased fecundity is linked to the success of this urban pest.

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 Arms Race Between KRAB–Zinc Finger Proteins and Endogenous Retroelements and Its Impact on Mammals

Nearly half of the human genome consists of endogenous retroelements (EREs) and their genetic remnants, a small fraction of which carry the potential to propagate in the host genome, posing a threat to genome integrity and cell/organismal survival. The largest family of transcription factors in tetrapods, the Krüppel-associated box domain zinc finger proteins (KRAB-ZFPs), binds to specific EREs and represses their transcription. Since their first appearance over 400 million years ago, KRAB-ZFPs have undergone dramatic expansion and diversification in mammals, correlating with the invasions of new EREs. In this article we review our current understanding of the structure, function, and evolution of KRAB-ZFPs and discuss growing evidence that the arms race between KRAB-ZFPs and the EREs they target is a major driving force for the evolution of new traits in mammals, often accompanied by domestication of EREs themselves.

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.

Mechanisms of noncanonical binding dynamics in multivalent protein-protein interactions

Protein multivalency can provide increased affinity and specificity relative to monovalent counterparts, but these emergent biochemical properties and their mechanistic underpinnings are difficult to predict as a function of the biophysical properties of the multivalent binding partners. Here, we present a mathematical model that accurately simulates binding kinetics and equilibria of multivalent protein-protein interactions as a function of the kinetics of monomer-monomer binding, the structure and topology of the multidomain interacting partners, and the valency of each partner. These properties are all experimentally or computationally estimated a priori, including approximating topology with a worm-like chain model applicable to a variety of structurally disparate systems, thus making the model predictive without parameter fitting. We conceptualize multivalent binding as a protein-protein interaction network: ligand and receptor valencies determine the number of interacting species in the network, with monomer kinetics and structural properties dictating the dynamics of each species. As predicted by the model and validated by surface plasmon resonance experiments, multivalent interactions can generate several noncanonical macroscopic binding dynamics, including a transient burst of high-energy configurations during association, biphasic equilibria resulting from interligand competition at high concentrations, and multiexponential dissociation arising from differential lifetimes of distinct network species. The transient burst was only uncovered when extending our analysis to trivalent interactions due to the significantly larger network, and we were able to predictably tune burst magnitude by altering linker rigidity. This study elucidates mechanisms of multivalent binding and establishes a framework for model-guided analysis and engineering of such interactions.

Host-transposon interactions: conflict, cooperation, and cooption

Transposable elements (TEs) are mobile DNA sequences that colonize genomes and threaten genome integrity. As a result, several mechanisms appear to have emerged during eukaryotic evolution to suppress TE activity. However, TEs are ubiquitous and account for a prominent fraction of most eukaryotic genomes. We argue that the evolutionary success of TEs cannot be explained solely by evasion from host control mechanisms. Rather, some TEs have evolved commensal and even mutualistic strategies that mitigate the cost of their propagation. These coevolutionary processes promote the emergence of complex cellular activities, which in turn pave the way for cooption of TE sequences for organismal function.

A conserved regulatory program initiates lateral plate mesoderm emergence across chordates

Cardiovascular lineages develop together with kidney, smooth muscle, and limb connective tissue progenitors from the lateral plate mesoderm (LPM). How the LPM initially emerges and how its downstream fates are molecularly interconnected remain unknown. Here, we isolate a pan-LPM enhancer in the zebrafish-specific draculin (drl) gene that provides specific LPM reporter activity from early gastrulation. In toto live imaging and lineage tracing of drl-based reporters captures the dynamic LPM emergence as lineage-restricted mesendoderm field. The drl pan-LPM enhancer responds to the transcription factors EomesoderminA, FoxH1, and MixL1 that combined with Smad activity drive LPM emergence. We uncover specific activity of zebrafish-derived drl reporters in LPM-corresponding territories of several chordates including chicken, axolotl, lamprey, Ciona, and amphioxus, revealing a universal upstream LPM program. Altogether, our work provides a mechanistic framework for LPM emergence as defined progenitor field, possibly representing an ancient mesodermal cell state that predates the primordial vertebrate embryo.

Numerous tissues are derived from the lateral plate mesoderm (LPM) but how this is specified is unclear. Here, the authors identify a pan-LPM reporter activity found in the zebrafish draculin (drl) gene that also shows transgenic activity in LPM-corresponding territories of several chordates, including chicken, axolotl, lamprey, Ciona, and amphioxus.

ZFP57 regulation of transposable elements and gene expression within and beyond imprinted domains

BACKGROUND

KRAB zinc finger proteins (KZFPs) represent one of the largest families of DNA-binding proteins in vertebrate genomes and appear to have evolved to silence transposable elements (TEs) including endogenous retroviruses through sequence-specific targeting of repressive chromatin states. ZFP57 is required to maintain the post-fertilization DNA methylation memory of parental origin at genomic imprints. Here we conduct RNA-seq and ChIP-seq analyses in normal and ZFP57 mutant mouse ES cells to understand the relative importance of ZFP57 at imprints, unique and repetitive regions of the genome.

RESULTS

Over 80% of ZFP57 targets are TEs, however, ZFP57 is not essential for their repression. The remaining targets lie within unique imprinted and non-imprinted sequences. Though the loss of ZFP57 influences imprinted genes as expected, the majority of unique gene targets lose H3K9me3 with little effect on DNA methylation and very few exhibit alterations in expression. Comparison of ZFP57 mutants with DNA methyltransferase-deleted ES cells (TKO) identifies a remarkably similar pattern of H3K9me3 loss across the genome. These data define regions where H3K9me3 is secondary to DNA methylation and we propose that ZFP57 is the principal if not sole methylation-sensitive KZFP in mouse ES cells. Finally, we examine dynamics of DNA and H3K9 methylation during pre-implantation development and show that sites bound by ZFP57 in ES cells maintain DNA methylation and H3K9me3 at imprints and at non-imprinted regions on the maternally inherited chromosome throughout preimplantation development.

CONCLUSION

Our analyses suggest the evolution of a rare DNA methylation-sensitive KZFP that is not essential for repeat silencing, but whose primary function is to maintain DNA methylation and repressive histone marks at germline-derived imprinting control regions.

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.

Plasticity of Promoter-Core Sequences Allows Bacteria to Compensate for the Loss of a Key Global Regulatory Gene

Transcription regulatory networks (TRNs) are of central importance for both short-term phenotypic adaptation in response to environmental fluctuations and long-term evolutionary adaptation, with global regulatory genes often being targets of natural selection in laboratory experiments. Here, we combined evolution experiments, whole-genome resequencing, and molecular genetics to investigate the driving forces, genetic constraints, and molecular mechanisms that dictate how bacteria can cope with a drastic perturbation of their TRNs. The crp gene, encoding a major global regulator in Escherichia coli, was deleted in four different genetic backgrounds, all derived from the Long-Term Evolution Experiment (LTEE) but with different TRN architectures. We confirmed that crp deletion had a more deleterious effect on growth rate in the LTEE-adapted genotypes; and we showed that the ptsG gene, which encodes the major glucose-PTS transporter, gained CRP (cyclic AMP receptor protein) dependence over time in the LTEE. We then further evolved the four crp-deleted genotypes in glucose minimal medium, and we found that they all quickly recovered from their growth defects by increasing glucose uptake. We showed that this recovery was specific to the selective environment and consistently relied on mutations in the cis-regulatory region of ptsG, regardless of the initial genotype. These mutations affected the interplay of transcription factors acting at the promoters, changed the intrinsic properties of the existing promoters, or produced new transcription initiation sites. Therefore, the plasticity of even a single promoter region can compensate by three different mechanisms for the loss of a key regulatory hub in the E. coli TRN.

The KRAB-zinc-finger protein ZFP708 mediates epigenetic repression at RMER19B retrotransposons

Global epigenetic reprogramming is vital to purge germ cell-specific epigenetic features to establish the totipotent state of the embryo. This process transpires to be carefully regulated and is not an undirected, radical erasure of parental epigenomes. The TRIM28 complex has been shown to be crucial in embryonic epigenetic reprogramming by regionally opposing DNA demethylation to preserve vital parental information to be inherited from germline to soma. Yet the DNA-binding factors guiding this complex to specific targets are largely unknown. Here, we uncover and characterize a novel, maternally expressed, TRIM28-interacting KRAB zinc-finger protein: ZFP708. It recruits the repressive TRIM28 complex to RMER19B retrotransposons to evoke regional heterochromatin formation. ZFP708 binding to these hitherto unknown TRIM28 targets is DNA methylation and H3K9me3 independent. ZFP708 mutant mice are viable and fertile, yet embryos fail to inherit and maintain DNA methylation at ZFP708 target sites. This can result in activation of RMER19B-adjacent genes, while ectopic expression of ZFP708 results in transcriptional repression. Finally, we describe the evolutionary conservation of ZFP708 in mice and rats, which is linked to the conserved presence of the targeted RMER19B retrotransposons in these species.

Summary: Analysis of the function and targets of a maternal KRAB-zinc-finger protein, ZFP708, found to specifically mediate maintenance of DNA methylation at a subset of LTR retrotransposons during embryonic epigenetic reprogramming.

CRISPR/Cas: An intriguing genomic editing tool with prospects in treating neurodegenerative diseases

The CRISPR/Cas genome editing tool has led to a revolution in biological research. Its ability to target multiple genomic loci simultaneously allows its application in gene function and genomic manipulation studies. Its involvement in the sequence specific gene editing in different backgrounds has changed the scenario of treating genetic diseases. By unravelling the mysteries behind complex neuronal circuits, it not only paved way in understanding the pathogenesis of the disease but helped in the development of large animal models of different neuronal diseases; thereby opened the gateways of successfully treating different neuronal diseases. This review explored the possibility of using of CRISPR/Cas in engineering DNA at the embryonic stage, as well as during the functioning of different cell types in the brain, to delineate implications related to the use of this super-specialized genome editing tool to overcome various neurodegenerative diseases that arise as a result of genetic mutations.

SETDB1-dependent heterochromatin stimulates alternative lengthening of telomeres

Alternative lengthening of telomeres, or ALT, is a recombination-based process that maintains telomeres to render some cancer cells immortal. The prevailing view is that ALT is inhibited by heterochromatin because heterochromatin prevents recombination. To test this model, we used telomere-specific quantitative proteomics on cells with heterochromatin deficiencies. In contrast to expectations, we found that ALT does not result from a lack of heterochromatin; rather, ALT is a consequence of heterochromatin formation at telomeres, which is seeded by the histone methyltransferase SETDB1. Heterochromatin stimulates transcriptional elongation at telomeres together with the recruitment of recombination factors, while disrupting heterochromatin had the opposite effect. Consistently, loss of SETDB1, disrupts telomeric heterochromatin and abrogates ALT. Thus, inhibiting telomeric heterochromatin formation in ALT cells might offer a new therapeutic approach to cancer treatment.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

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.

Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation

BACKGROUND

Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes.

FINDINGS

Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated.

CONCLUSION

The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

The Human Transcription Factors

Transcription factors (TFs) recognize specific DNA sequences to control chromatin and transcription, forming a complex system that guides expression of the genome. Despite keen interest in understanding how TFs control gene expression, it remains challenging to determine how the precise genomic binding sites of TFs are specified and how TF binding ultimately relates to regulation of transcription. This review considers how TFs are identified and functionally characterized, principally through the lens of a catalog of over 1,600 likely human TFs and binding motifs for two-thirds of them. Major classes of human TFs differ markedly in their evolutionary trajectories and expression patterns, underscoring distinct functions. TFs likewise underlie many different aspects of human physiology, disease, and variation, highlighting the importance of continued effort to understand TF-mediated gene regulation.

Sequences encoding C2H2 zinc fingers inhibit polyadenylation and mRNA export in human cells

The large C2H2-Zinc Finger (C2H2-ZNF) gene family has rapidly expanded in primates through gene duplication. There is consequently considerable sequence homology between family members at both the nucleotide and amino acid level, allowing for coordinated regulation and shared functions. Here we show that multiple C2H2-ZNF mRNAs experience differential polyadenylation resulting in populations with short and long poly(A) tails. Furthermore, a significant proportion of C2H2-ZNF mRNAs are retained in the nucleus. Intriguingly, both short poly(A) tails and nuclear retention can be specified by the repeated elements that encode zinc finger motifs. These Zinc finger Coding Regions (ZCRs) appear to restrict polyadenylation of nascent RNAs and at the same time impede their export. However, the polyadenylation process is not necessary for nuclear retention of ZNF mRNAs. We propose that inefficient polyadenylation and export may allow C2H2-ZNF mRNAs to moonlight as non-coding RNAs or to be stored for later use.

C2H2-Type Zinc Finger Proteins: Evolutionarily Old and New Partners of the Nuclear Hormone Receptors

Nuclear hormone receptors (NRs) are evolutionarily conserved ligand-dependent transcription factors. They are essential for human life, mediating the actions of lipophilic molecules, such as steroid hormones and metabolites of fatty acid, cholesterol, and external toxic compounds. The C2H2-type zinc finger proteins (ZNFs) form the largest family of the transcription factors in humans and are characterized by multiple, tandemly arranged zinc fingers. Many of the C2H2-type ZNFs are conserved throughout evolution, suggesting their involvement in preserved biological activities, such as general transcriptional regulation and development/differentiation of organs/tissues observed in the early embryonic phase. However, some C2H2-type ZNFs, such as those with the Krüppel-associated box (KRAB) domain, appeared relatively late in evolution and have significantly increased family members in mammals including humans, possibly modulating their complicated transcriptional network and/or supporting the morphological development/functions specific to them. Such evolutional characteristics of the C2H2-type ZNFs indicate that these molecules influence the NR functions conserved through evolution, whereas some also adjust them to meet with specific needs of higher organisms. We review the interaction between NRs and C2H2-type ZNFs by focusing on some of the latter molecules.

Sequence Determinants for Nuclear Retention and Cytoplasmic Export of mRNAs and lncRNAs

Eukaryotes are divided into two major compartments: the nucleus where RNA is synthesized and processed, and the cytoplasm, where mRNA is translated into proteins. Although many different RNAs are made, only a subset is allowed access to the cytoplasm, primarily RNAs involved in protein synthesis (mRNA, tRNA, and rRNA). In contrast, nuclear retained transcripts are mostly long non-coding RNAs (lncRNAs) whose role in cell physiology has been a source of much investigation in the past few years. In addition, it is likely that many non-functional RNAs, which arise by spurious transcription and misprocessing of functional RNAs, are also retained in the nucleus and degraded. In this review, the main sequence features that dictate whether any particular mRNA or lncRNA is a substrate for retention in the nucleus, or export to the cytoplasm, are discussed. Although nuclear export is promoted by RNA-splicing due to the fact that the spliceosome can help recruit export factors to the mature RNA, nuclear export does not require splicing. Indeed, most stable unspliced transcripts are well exported and associate with these same export factors in a splicing-independent manner. In contrast, nuclear retention is promoted by specialized cis-elements found in certain RNAs. This new understanding of the determinants of nuclear retention and cytoplasmic export provides a deeper understanding of how information flow is regulated in eukaryotic cells. Ultimately these processes promote the evolution of complexity in eukaryotes by shaping the genomic content through constructive neutral evolution.

Genome-wide survey, characterization, and expression analysis of RING finger protein genes in Brassica oleracea and their syntenic comparison to Brassica rapa and Arabidopsis thaliana

The ubiquitin-mediated post-translational regulatory pathway regulates a broad range of cell functions in all eukaryotes. It requires the involvement of a large number of E3 ligases, of which more than one third belongs to the RING protein family as in Arabidopsis thaliana. In this study, a total of 756 RING domains in 734 predicted proteins were identified in Brassica oleracea. Their encoding genes were characterized by RING domain type, additional domain, and expression pattern, and mapped on the nine chromosomes of B. oleracea. Comparison of these results with B. rapa and A. thaliana revealed some common as well as species-specific features. Our results showed that the differential gene loss following the whole genome triplication has largely contributed to the RING protein gene number variation among these species, although other factors such as tandem duplication, RING domain loss, or modification had also contributed to this variation. Analysis of RNA-seq data showed that these RING protein genes were functionally diversified and involved in all the stages of plant growth and development, and that the triplicated members were also diverged in expression with one member often more dominantly expressed over the two others in the majority of cases. Our study lays the foundation for further functional determination of each RING protein gene among species of the genus Brassica.

PRDM9, a driver of the genetic map

During meiosis, maternal and paternal chromosomes undergo exchanges by homologous recombination. This is essential for fertility and contributes to genome evolution. In many eukaryotes, sites of meiotic recombination, also called hotspots, are regions of accessible chromatin, but in many vertebrates, their location follows a distinct pattern and is specified by PR domain-containing protein 9 (PRDM9). The specification of meiotic recombination hotspots is achieved by the different activities of PRDM9: DNA binding, histone methyltransferase, and interaction with other proteins. Remarkably, PRDM9 activity leads to the erosion of its own binding sites and the rapid evolution of its DNA-binding domain. PRDM9 may also contribute to reproductive isolation, as it is involved in hybrid sterility potentially due to a reduction of its activity in specific heterozygous contexts.

Expression and functional evaluation of CaZNF830 during pepper response to Ralstonia solanacearum or high temperature and humidity

Extensive transcriptional reprogramming after pathogen attack determines immunity to these invaders and plant development. Zinc finger (ZNF) transcription factors regulate important processes in plants such as development, vegetative activities and plant immunity. Despite their immense significance, majority of ZNF transcription factors (TF) involved in pepper immunity and resistance to heat stress have not been focused much. Herein, we identified and functionally characterized CaZNF830 in pepper defense against Ralstonia solanacearum inoculation (RSI) and tolerance to high temperature and high humidity (HTHH). Transient expression analysis of CaZNF830-GFP fusion protein in tobacco leaves revealed its localization to the nucleus. Transcription of CaZNF830 is induced in pepper plants upon RSI or HTHH. Consistent with this, fluorometric GUS enzymatic assay driven by pCaZNF830 presented significantly enhanced activity under RSI and HTHH in comparison with the control plants. The silencing of CaZNF830 by virus induced gene silencing (VIGS) significantly compromised pepper immunity against RSI with enhanced growth of Ralstonia solanacearum in pepper plants. Silencing of CaZNF830 also impaired tolerance to HTHH coupled with decreased expression levels of immunity and thermo-tolerance associated marker genes including CaHIR1, CaNPR1, CaPR1, CaABR1 and CaHSP24. By contrast, the transient over-expression of CaZNF830 in pepper leaves by infiltration of GV3101 cells containing 35S::CaZNF830-HA induced HR mimic cell death, H2O2 accumulation and activated the transcriptions of the tested defense-relative or thermo-tolerance associated marker genes. RT-PCR and immune-blotting assay confirmed the stable expression of HA-tagged CaZNF830 mRNA and protein in pepper. All these results suggest that CaZNF830 acts as a positive regulator of plant immunity against RSI or tolerance to HTHH, it is induced by RSI or HTHH and consequently activate pepper immunity against RSI or tolerance to HTHH by directly or indirectly transcriptional modulation of many defense-linked genes.

Mammalian transposable elements and their impacts on genome evolution

Transposable elements (TEs) are genetic elements with the ability to mobilize and replicate themselves in a genome. Mammalian genomes are dominated by TEs, which can reach copy numbers in the hundreds of thousands. As a result, TEs have had significant impacts on mammalian evolution. Here we summarize the current understanding of TE content in mammal genomes and find that, with a few exceptions, most fall within a predictable range of observations. First, one third to one half of the genome is derived from TEs. Second, most mammalian genomes are dominated by LINE and SINE retrotransposons, more limited LTR retrotransposons, and minimal DNA transposon accumulation. Third, most mammal genome contains at least one family of actively accumulating retrotransposon. Finally, horizontal transfer of TEs among lineages is rare. TE exaptation events are being recognized with increasing frequency. Despite these beneficial aspects of TE content and activity, the majority of TE insertions are neutral or deleterious. To limit the deleterious effects of TE proliferation, the genome has evolved several defense mechanisms that act at the epigenetic, transcriptional, and post-transcriptional levels. The interaction between TEs and these defense mechanisms has led to an evolutionary arms race where TEs are suppressed, evolve to escape suppression, then are suppressed again as the defense mechanisms undergo compensatory change. The result is complex and constantly evolving interactions between TEs and host genomes.

Integrated analysis of motif activity and gene expression changes of transcription factors

The ability to predict transcription factors based on sequence information in regulatory elements is a key step in systems-level investigation of transcriptional regulation. Here, we have developed a novel tool, IMAGE, for precise prediction of causal transcription factors based on transcriptome profiling and genome-wide maps of enhancer activity. High precision is obtained by combining a near-complete database of position weight matrices (PWMs), generated by compiling public databases and systematic prediction of PWMs for uncharacterized transcription factors, with a state-of-the-art method for PWM scoring and a novel machine learning strategy, based on both enhancers and promoters, to predict the contribution of motifs to transcriptional activity. We applied IMAGE to published data obtained during 3T3-L1 adipocyte differentiation and showed that IMAGE predicts causal transcriptional regulators of this process with higher confidence than existing methods. Furthermore, we generated genome-wide maps of enhancer activity and transcripts during human mesenchymal stem cell commitment and adipocyte differentiation and used IMAGE to identify positive and negative transcriptional regulators of this process. Collectively, our results demonstrate that IMAGE is a powerful and precise method for prediction of regulators of gene expression.

Waking up dormant tumor suppressor genes with zinc fingers, TALEs and the CRISPR/dCas9 system

The aberrant epigenetic silencing of tumor suppressor genes (TSGs) plays a major role during carcinogenesis and regaining these dormant functions by engineering of sequence-specific epigenome editing tools offers a unique opportunity for targeted therapies. However, effectively normalizing the expression and regaining tumor suppressive functions of silenced TSGs by artificial transcription factors (ATFs) still remains a major challenge. Herein we describe novel combinatorial strategies for the potent reactivation of two class II TSGs, MASPIN and REPRIMO, in cell lines with varying epigenetic states, using the CRISPR/dCas9 associated system linked to a panel of effector domains (VP64, p300, VPR and SAM complex), as well as with protein-based ATFs, Zinc Fingers and TALEs. We found that co-delivery of multiple effector domains using a combination of CRISPR/dCas9 and TALEs or SAM complex maximized activation in highly methylated promoters. In particular, CRISPR/dCas9 VPR with SAM upregulated MASPIN mRNA (22,145-fold change) in H157 lung cancer cells, with accompanying re-expression of MASPIN protein, which led to a concomitant inhibition of cell proliferation and induction of apoptotic cell death. Consistently, CRISPR/dCas9 VP64 with SAM upregulated REPRIMO (680-fold change), which led to phenotypic reprogramming in AGS gastric cancer cells. Altogether, our results outlined novel sequence-specific, combinatorial epigenome editing approaches to reactivate highly methylated TSGs as a promising therapy for cancer and other diseases.

ZFP57 and the Targeted Maintenance of Postfertilization Genomic Imprints

Epigenetic modifications play an important role in modulating genome function. In mammals, inappropriate epigenetic states can cause embryonic lethality and various acquired and inherited diseases; hence, it is important to understand how such states are formed and maintained in particular genomic contexts. Genomic imprinting is a process in which epigenetic states provide a sustained memory of parental origin and cause gene expression/repression from only one of the two parental chromosomes. Genomic imprinting is therefore a valuable model to decipher the principles and processes associated with the targeting and maintenance of epigenetic states in general. Krüppel-associated box zinc finger proteins (KRAB-ZFPs) are proteins that have the potential to mediate this. ZFP57, one of the best characterized proteins in this family, has been shown to target and maintain epigenetic states at imprinting control regions after fertilization. Its role in imprinting through the use of ZFP57 mutants in mouse and the wider implications of KRAB-ZFPs for the targeted maintenance of epigenetic states are discussed here.

Comparison of ChIP-Seq Data and a Reference Motif Set for Human KRAB C2H2 Zinc Finger Proteins

KRAB C2H2 zinc finger proteins (KZNFs) are the largest and most diverse family of human transcription factors, likely due to diversifying selection driven by novel endogenous retroelements (EREs), but the vast majority lack binding motifs or functional data. Two recent studies analyzed a majority of the human KZNFs using either ChIP-seq (60 proteins) or ChIP-exo (221 proteins) in the same cell type (HEK293). The ChIP-exo paper did not describe binding motifs, however. Thirty-nine proteins are represented in both studies, enabling the systematic comparison of the data sets presented here. Typically, only a minority of peaks overlap, but the two studies nonetheless display significant similarity in ERE binding for 32/39, and yield highly similar DNA binding motifs for 23 and related motifs for 34 (MoSBAT similarity score >0.5 and >0.2, respectively). Thus, there is overall (albeit imperfect) agreement between the two studies. For the 242 proteins represented in at least one study, we selected a highest-confidence motif for each protein, utilizing several motif-derivation approaches, and evaluating motifs within and across data sets. Peaks for the majority (158) are enriched (96% with AUC >0.6 predicting peak vs. nonpeak) for a motif that is supported by the C2H2 “recognition code,” consistent with intrinsic sequence specificity driving DNA binding in cells. An additional 63 yield motifs enriched in peaks, but not supported by the recognition code, which could reflect indirect binding. Altogether, these analyses validate both data sets, and provide a reference motif set with associated quality metrics.

Zinc Fingers, TALEs, and CRISPR Systems: A Comparison of Tools for Epigenome Editing

The completion of genome, epigenome, and transcriptome mapping in multiple cell types has created a demand for precision biomolecular tools that allow researchers to functionally manipulate DNA, reconfigure chromatin structure, and ultimately reshape gene expression patterns. Epigenetic editing tools provide the ability to interrogate the relationship between epigenetic modifications and gene expression. Importantly, this information can be exploited to reprogram cell fate for both basic research and therapeutic applications. Three different molecular platforms for epigenetic editing have been developed: zinc finger proteins (ZFs), transcription activator-like effectors (TALEs), and the system of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR-associated (Cas) proteins. These platforms serve as custom DNA-binding domains (DBDs), which are fused to epigenetic modifying domains to manipulate epigenetic marks at specific sites in the genome. The addition and/or removal of epigenetic modifications reconfigures local chromatin structure, with the potential to provoke long-lasting changes in gene transcription. Here we summarize the molecular structure and mechanism of action of ZF, TALE, and CRISPR platforms and describe their applications for the locus-specific manipulation of the epigenome. The advantages and disadvantages of each platform will be discussed with regard to genomic specificity, potency in regulating gene expression, and reprogramming cell phenotypes, as well as ease of design, construction, and delivery. Finally, we outline potential applications for these tools in molecular biology and biomedicine and identify possible barriers to their future clinical implementation.

TRIM28 and Interacting KRAB-ZNFs Control Self-Renewal of Human Pluripotent Stem Cells through Epigenetic Repression of Pro-differentiation Genes

Reprogramming to induced pluripotent stem cells (iPSCs) and differentiation of pluripotent stem cells (PSCs) are regulated by epigenetic machinery. Tripartite motif protein 28 (TRIM28), a universal mediator of Krüppel-associated box domain zinc fingers (KRAB-ZNFs), is known to regulate both processes; however, the exact mechanism and identity of participating KRAB-ZNF genes remain unknown. Here, using a reporter system, we show that TRIM28/KRAB-ZNFs alter DNA methylation patterns in addition to H3K9me3 to cause stable gene repression during reprogramming. Using several expression datasets, we identified KRAB-ZNFs (ZNF114, ZNF483, ZNF589) in the human genome that maintain pluripotency. Moreover, we identified target genes repressed by these KRAB-ZNFs. Mechanistically, we demonstrated that these KRAB-ZNFs directly alter gene expression of important developmental genes by modulating H3K9me3 and DNA methylation of their promoters. In summary, TRIM28 employs KRAB-ZNFs to evoke epigenetic silencing of its target differentiation genes via H3K9me3 and DNA methylation.

•

Upon reprogramming KRAB-repressor evokes stable silencing of its target genes

•

KRAB-ZNFs repress target genes required for differentiation of pluripotent cells

•

KRAB-ZNFs are crucial for the maintenance of pluripotency of human stem cells

•

TRIM28/KRAB-ZNFs repress developmental genes through H3K9 and de novo DNA methylation

Comparative analysis of weighted gene co-expression networks in human and mouse

The application of complex network modeling to analyze large co-expression data sets has gained traction during the last decade. In particular, the use of the weighted gene co-expression network analysis framework has allowed an unbiased and systems-level investigation of genotype-phenotype relationships in a wide range of systems. Since mouse is an important model organism for biomedical research on human disease, it is of great interest to identify similarities and differences in the functional roles of human and mouse orthologous genes. Here, we develop a novel network comparison approach which we demonstrate by comparing two gene-expression data sets from a large number of human and mouse tissues. The method uses weighted topological overlap alongside the recently developed network-decomposition method of s-core analysis, which is suitable for making gene-centrality rankings for weighted networks. The aim is to identify globally central genes separately in the human and mouse networks. By comparing the ranked gene lists, we identify genes that display conserved or diverged centrality-characteristics across the networks. This framework only assumes a single threshold value that is chosen from a statistical analysis, and it may be applied to arbitrary network structures and edge-weight distributions, also outside the context of biology. When conducting the comparative network analysis, both within and across the two species, we find a clear pattern of enrichment of transcription factors, for the homeobox domain in particular, among the globally central genes. We also perform gene-ontology term enrichment analysis and look at disease-related genes for the separate networks as well as the network comparisons. We find that gene ontology terms related to regulation and development are generally enriched across the networks. In particular, the genes FOXE3, RHO, RUNX2, ALX3 and RARA, which are disease genes in either human or mouse, are on the top-10 list of globally central genes in the human and mouse networks.

The emergence of piRNAs against transposon invasion to preserve mammalian genome integrity

Transposable elements (TEs) contribute to the large amount of repetitive sequences in mammalian genomes and have been linked to species-specific genome innovations by rewiring regulatory circuitries. However, organisms need to restrict TE activity to ensure genome integrity, especially in germline cells to protect the transmission of genetic information to the next generation. This review features our current understandings of mammalian PIWI-interacting RNAs (piRNAs) and their role in TE regulation in spermatogenesis. Here we discuss functional implication and explore additional molecular mechanisms that inhibit transposon activity and altogether illustrate the paradoxical arms race between genome evolution and stability.