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

Transcription factor binding specificities of the oomycete Phytophthora infestans reflect conserved and divergent evolutionary patterns and predict function

BACKGROUND

Identifying the DNA-binding specificities of transcription factors (TF) is central to understanding gene networks that regulate growth and development. Such knowledge is lacking in oomycetes, a microbial eukaryotic lineage within the stramenopile group. Oomycetes include many important plant and animal pathogens such as the potato and tomato blight agent Phytophthora infestans, which is a tractable model for studying life-stage differentiation within the group.

RESULTS

Mining of the P. infestans genome identified 197 genes encoding proteins belonging to 22 TF families. Their chromosomal distribution was consistent with family expansions through unequal crossing-over, which were likely ancient since each family had similar sizes in most oomycetes. Most TFs exhibited dynamic changes in RNA levels through the P. infestans life cycle. The DNA-binding preferences of 123 proteins were assayed using protein-binding oligonucleotide microarrays, which succeeded with 73 proteins from 14 families. Binding sites predicted for representatives of the families were validated by electrophoretic mobility shift or chromatin immunoprecipitation assays. Consistent with the substantial evolutionary distance of oomycetes from traditional model organisms, only a subset of the DNA-binding preferences resembled those of human or plant orthologs. Phylogenetic analyses of the TF families within P. infestans often discriminated clades with canonical and novel DNA targets. Paralogs with similar binding preferences frequently had distinct patterns of expression suggestive of functional divergence. TFs were predicted to either drive life stage-specific expression or serve as general activators based on the representation of their binding sites within total or developmentally-regulated promoters. This projection was confirmed for one TF using synthetic and mutated promoters fused to reporter genes in vivo.

CONCLUSIONS

We established a large dataset of binding specificities for P. infestans TFs, representing the first in the stramenopile group. This resource provides a basis for understanding transcriptional regulation by linking TFs with their targets, which should help delineate the molecular components of processes such as sporulation and host infection. Our work also yielded insight into TF evolution during the eukaryotic radiation, revealing both functional conservation as well as diversification across kingdoms.

Chromatin accessibility in the Drosophila embryo is determined by transcription factor pioneering and enhancer activation

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation.

An RFX transcription factor regulates ciliogenesis in the closest living relatives of animals

Cilia allowed our protistan ancestors to sense and explore their environment, avoid predation, and capture bacterial prey.1,2,3 Regulated ciliogenesis was likely critical for early animal evolution,2,4,5,6 and in modern animals, deploying cilia in the right cells at the right time is crucial for development and physiology. Two transcription factors, RFX and FoxJ1, coordinate ciliogenesis in animals7,8,9 but are absent from the genomes of many other ciliated eukaryotes, raising the question of how the regulation of ciliogenesis in animals evolved.10,11 By comparing the genomes of animals with those of their closest living relatives, the choanoflagellates, we found that the genome of their last common ancestor encoded at least three RFX paralogs and a FoxJ1 homolog. Disruption of the RFX homolog cRFXa in the model choanoflagellate Salpingoeca rosetta resulted in delayed cell proliferation and aberrant ciliogenesis, marked by the collapse and resorption of nascent cilia. In cRFXa mutants, ciliogenesis genes and foxJ1 were significantly downregulated. Moreover, the promoters of S. rosetta ciliary genes are enriched for DNA motifs matching those bound by the cRFXa protein in vitro. These findings suggest that an ancestral cRFXa homolog coordinated ciliogenesis in the progenitors of animals and choanoflagellates and that the selective deployment of the RFX regulatory module may have been necessary to differentiate ciliated from non-ciliated cell types during early animal evolution.

Use of gene therapy for optic nerve protection: Current concepts

Gene therapy has become an essential treatment for optic nerve injury (ONI) in recent years, and great strides have been made using animal models. ONI, which is characterized by the loss of retinal ganglion cells (RGCs) and axons, can induce abnormalities in the pupil light reflex, visual field defects, and even vision loss. The eye is a natural organ to target with gene therapy because of its high accessibility and certain immune privilege. As such, numerous gene therapy trials are underway for treating eye diseases such as glaucoma. The aim of this review was to cover research progress made in gene therapy for ONI. Specifically, we focus on the potential of gene therapy to prevent the progression of neurodegenerative diseases and protect both RGCs and axons. We cover the basic information of gene therapy, including the classification of gene therapy, especially focusing on genome editing therapy, and then we introduce common editing tools and vector tools such as Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) -Cas9 and adeno-associated virus (AAV). We also summarize the progress made on understanding the roles of brain derived neurotrophic factor (BDNF), ciliary neurotrophic factor (CNTF), phosphatase-tensin homolog (PTEN), suppressor of cytokine signal transduction 3 (SOCS3), histone acetyltransferases (HATs), and other important molecules in optic nerve protection. However, gene therapy still has many challenges, such as misalignment and mutations, immunogenicity of AAV, time it takes and economic cost involved, which means that these issues need to be addressed before clinical trials can be considered.

MtTRC-1, a Novel Transcription Factor, Regulates Cellulase Production via Directly Modulating the Genes Expression of the Mthac-1 and Mtcbh-1 in Myceliophthora thermophila

The thermophilic fungus Myceliophthora thermophila has been used to produce industrial enzymes and biobased chemicals. In saprotrophic fungi, the mechanisms regulating cellulase production have been studied, which revealed the involvement of multiple transcription factors. However, in M. thermophila, the transcription factors influencing cellulase gene expression and secretion remain largely unknown. In this study, we identified and characterized a novel cellulase regulator (MtTRC-1) in M. thermophila through a combination of functional genomics and genetic analyses. Deletion of Mttrc-1 resulted in significantly decreased cellulase production and activities. Transcriptome analysis revealed downregulation of not only the encoding genes of main cellulases but also the transcriptional regulator MtHAC-1 of UPR pathway after disruption of MtTRC-1 under cellulolytic induction conditions. Herein, we also characterized the ortholog of the yeast HAC1p in M. thermophila. We show that Mthac-1 mRNA undergoes an endoplasmic reticulum (ER) stress-induced splicing by removing a 23-nucleotide (nt) intron. Notably, the protein secretion on cellulose was dramatically impaired by the deletion of MtHAC-1. Moreover, the colonial growth on various carbon sources was defective in the absence of MtHAC-1. Electrophoretic mobility shift assays and chromatin immunoprecipitation assays verified MtTRC-1 regulates the transcription of Mthac-1 and the major cellulase gene Mtcbh-1 by binding directly to the promoters in vitro and in vivo. Furthermore, DNase I footprinting assays identified the putative consensus binding site (5′-GNG/C-3′). These results revealed the importance of MtTRC-1 for positively regulating cellulase production. This finding has clarified the complex regulatory pathways involved in cellulolytic enzyme production.

IMPORTANCE In the present study, we characterized a novel regulator MtTRC-1 in M. thermophila, which regulated cellulase production through direct transcriptional regulation of the Mthac-1 and Mtcbh-1 genes. Our data demonstrated that MtHAC-1 is a key factor for the cellulase secretion capacity of M. thermophila. Our data indicate that this thermophilic fungus regulates cellulase production through a multilevels network, in which the protein secretory pathway is modulated by MtHAC-1-dependent UPR pathway and the cellulase gene expression is directly regulated in parallel by transcription factors. The conservation of Mttrc1 in filamentous fungi suggests this mechanism may be exploited to engineer filamentous fungal cell factories capable of producing proteins on an industrial scale.

Application of Gene Editing Technology in Resistance Breeding of Livestock

As a new genetic engineering technology, gene editing can precisely modify the specific gene sequence of the organism’s genome. In the last 10 years, with the rapid development of gene editing technology, zinc-finger nucleases (ZFNs), transcription activator-like endonucleases (TALENs), and CRISPR/Cas9 systems have been applied to modify endogenous genes in organisms accurately. Now, gene editing technology has been used in mice, zebrafish, pigs, cattle, goats, sheep, rabbits, monkeys, and other species. Breeding for disease-resistance in agricultural animals tends to be a difficult task for traditional breeding, but gene editing technology has made this easier. In this work, we overview the development and application of gene editing technology in the resistance breeding of livestock. Also, we further discuss the prospects and outlooks of gene editing technology in disease-resistance breeding.

Distinct structural bases for sequence-specific DNA binding by mammalian BEN domain proteins

The BEN domain is a recently recognized DNA binding module that is present in diverse metazoans and certain viruses. Several BEN domain factors are known as transcriptional repressors, but, overall, relatively little is known of how BEN factors identify their targets in humans. In particular, X-ray structures of BEN domain:DNA complexes are only known for Drosophila factors bearing a single BEN domain, which lack direct vertebrate orthologs. Here, we characterize several mammalian BEN domain (BD) factors, including from two NACC family BTB-BEN proteins and from BEND3, which has four BDs. In vitro selection data revealed sequence-specific binding activities of isolated BEN domains from all of these factors. We conducted detailed functional, genomic, and structural studies of BEND3. We show that BD4 is a major determinant for in vivo association and repression of endogenous BEND3 targets. We obtained a high-resolution structure of BEND3-BD4 bound to its preferred binding site, which reveals how BEND3 identifies cognate DNA targets and shows differences with one of its non-DNA-binding BEN domains (BD1). Finally, comparison with our previous invertebrate BEN structures, along with additional structural predictions using AlphaFold2 and RoseTTAFold, reveal distinct strategies for target DNA recognition by different types of BEN domain proteins. Together, these studies expand the DNA recognition activities of BEN factors and provide structural insights into sequence-specific DNA binding by mammalian BEN proteins.

Altering the Double-Stranded DNA Specificity of the bZIP Domain of Zta with Site-Directed Mutagenesis at N182

Zta, the Epstein-Barr virus bZIP transcription factor (TF), binds both unmethylated and methylated double-stranded DNA (dsDNA) in a sequence-specific manner. We studied the contribution of a conserved asparagine (N182) to sequence-specific dsDNA binding to four types of dsDNA: (i) dsDNA with cytosine in both strands ((DNA(C|C)), (ii, iii) dsDNA with 5-methylcytosine (5mC, M) or 5-hydroxymethylcytosine (5hmC, H) in one strand and cytosine in the second strand ((DNA(5mC|C) and DNA(5hmC|C)), and (iv) dsDNA with methylated cytosine in both strands in all CG dinucleotides ((DNA(5mCG)). We replaced asparagine with five similarly sized amino acids (glutamine (Q), serine (S), threonine (T), isoleucine (I), or valine (V)) and used protein binding microarrays to evaluate sequence-specific dsDNA binding. Zta preferentially binds the pseudo-palindrome TRE (AP1) motif (T^-4G^-3A^-2G/_C ⁰T²C³A⁴ ). Zta (N182Q) changes binding to A³ in only one half-site. Zta(N182S) changes binding to G³ in one or both halves of the motif. Zta(N182S) and Zta(N182Q) have 34- and 17-fold weaker median dsDNA binding, respectively. Zta(N182V) and Zta(N182I) have increased binding to dsDNA(5mC|C). Molecular dynamics simulations rationalize some of these results, identifying hydrogen bonds between glutamine and A³ , but do not reveal why serine preferentially binds G³ , suggesting that entropic interactions may mediate this new binding specificity.

REL Domain of NFATc2 Binding to Five Types of DNA Using Protein Binding Microarrays

NFATc2 is a DNA binding protein in the Rel family transcription factors, which binds a CGGAA motif better when both cytosines in the CG dinucleotide are methylated. Using protein binding microarrays (PBMs), we examined the DNA binding of NFATc2 to three additional types of DNA: single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) with either 5-methylcytosine (5mC, M) or 5-hydroxymethylcytosine (5hmC, H) in one strand and a cytosine in the second strand. ATTTCCAC, the complement of the core GGAA motif, is better bound as ssDNA compared to dsDNA. dsDNA containing the 5-mer CGGAA with either 5mC or 5hmC in one DNA strand is bound stronger than CGGAA. In contrast, the reverse complement TTCCG is bound weaker when it contains 5mC. Analysis of the available NFATc2:dsDNA complexes rationalizes these PBM data.

Crystal structure of the BRPF2 PWWP domain in complex with DNA reveals a different binding mode than the HDGF family of PWWP domains

The PWWP domain was first identified in the HDGF protein family and named after the conserved Proline-Tryptophan-Tryptophan-Proline motif in WHSC1. The PWWP domain-containing proteins play important roles in different biological processes, such as DNA replication, transcription, DNA repair, pre-mRNA processing by recognizing methylated histone and dsDNA simultaneously. Recently, how the HDGF family of PWWP domains recognize histone H3K36me3-modified nucleosome has been reported. In order to better understand the interactions between the PWWP domain and dsDNA, we carried out family-wide characterization of dsDNA binding abilities of human PWWP domains. Our binding assays confirmed that PWWP domains bind to dsDNA without sequence selectivity. Our crystal structure of the BRPF2 PWWP domain in complex with a 12-mer dsDNA reveals that the PWWP domain interacts with dsDNA by binding to its major groove, instead of the minor groove observed in the HDGF family of PWWP domains. Our study indicates that PWWP domains could bind to dsDNA in different modes.

bZIP Dimers CREB1, ATF2, Zta, ATF3|cJun, and cFos|cJun Prefer to Bind to Some Double-Stranded DNA Sequences Containing 5-Formylcytosine and 5-Carboxylcytosine

In mammalian cells, 5-methylcytosine (5mC) occurs in genomic double-stranded DNA (dsDNA) and is enzymatically oxidized to 5-hydroxymethylcytosine (5hmC), then to 5-formylcytosine (5fC), and finally to 5-carboxylcytosine (5caC). These cytosine modifications are enriched in regulatory regions of the genome. The effect of these oxidative products on five bZIP dimers (CREB1, ATF2, Zta, ATF3|cJun, and cFos|cJun) binding to five types of dsDNA was measured using protein binding microarrays. The five dsDNAs contain either cytosine in both DNA strands or cytosine in one strand and either 5mC, 5hmC, 5fC, or 5caC in the second strand. Some sequences containing the CEBP half-site GCAA are bound more strongly by all five bZIP domains when dsDNA contains 5mC, 5hmC, or 5fC. dsDNA containing 5caC in some TRE (AP-1)-like sequences, e.g., TGACTAA, is better bound by Zta, ATF3|cJun, and cFos|cJun.

Balance between asymmetry and abundance in multi-domain DNA-binding proteins may regulate the kinetics of their binding to DNA

DNA sequences are often recognized by multi-domain proteins that may have higher affinity and specificity than single-domain proteins. However, the higher affinity to DNA might be coupled with slower recognition kinetics. In this study, we address this balance between stability and kinetics for multi-domain Cys₂His₂- (C₂H₂-) type zinc-finger (ZF) proteins. These proteins are the most prevalent DNA-binding domain in eukaryotes and C₂H₂ type zinc-finger proteins (C₂H₂-ZFPs) constitute nearly one-half of all known and predicted transcription factors in human. Extensive contact with DNA via tandem ZF domains confers high stability on the sequence-specific complexes. However, this can limit target search efficiency, especially for low abundance ZFPs. Earlier, we found that asymmetrical distribution of electrostatic charge among the three ZF domains of the low abundance transcription factor Egr-1 facilitates its DNA search process. Here, on a diverse set of 273 human C₂H₂-ZFP comprised of 3–15 tandem ZF domains, we find that, in many cases, electrostatic charge and binding specificity are asymmetrically distributed among the ZF domains so that neighbouring domains have different DNA-binding properties. For proteins containing 3–6 ZF domains, we show that the low abundance proteins possess a higher degree of non-specific asymmetry and vice versa. Our findings suggest that where the electrostatics of tandem ZF domains are similar (i.e., symmetrical), the ZFPs are more abundant to optimize their DNA search efficiency. This study reveals new insights into the fundamental determinants of recognition by C₂H₂-ZFPs of their DNA binding sites in the cellular landscape. The importance of electrostatic asymmetry with respect to binding site recognition by C₂H₂-ZFPs suggests the possibility that it may also be important in other ZFP systems and reveals a new design feature for zinc finger engineering.

Optimal recognition of proteins to DNA is governed by various factors among them the thermodynamics, kinetics and specificity of the protein-DNA complex. Multi-domain DNA-binding proteins are expected to have higher affinity and specificity due to the extensive interface they form with DNA. However, larger interface may result with higher friction when these proteins scan the DNA for the target site via the sliding mechanism. A way to overcome this drawback is to have asymmetry in the protein so that the interface with DNA is smaller. Alternatively, higher abundance can also increase the search speed. Here, using computational analysis of large data set of multi-domain zinc finger DNA-binding proteins, we report a trade-off between asymmetry and abundance.

Applications of multiplex genome editing in higher plants

Multiplex genome editing involves the simultaneous targeting of multiple related or unrelated targets. The latter is most straightforward using the CRISPR/Cas9 system because multiple gRNAs can be delivered either as independent expression cassettes with their own promoters or as polycistronic transcripts processed into mature gRNAs by endogenous or introduced nucleases. Multiplex genome editing in plants initially focused on input traits such as herbicide resistance, but has recently expanded to include hormone biosynthesis and perception, metabolic engineering, plant development and molecular farming, with more than 100 simultaneous targeting events reported. Usually the coding region is targeted but recent examples also include promoter modifications to generate mutants with varying levels of gene expression.

Genome-wide identification and expression profiling of the C2H2-type zinc finger protein genes in the silkworm Bombyx mori

Cys2-His2 zinc finger (C2H2-ZF) proteins comprise the largest class of putative eukaryotic transcription factors. The zinc finger motif array is highly divergent, indicating that most proteins will have distinctive binding sites and perform different functions. However, the binding sites and functions of the majority of C2H2-ZF proteins remain unknown. In this study, we identified 327 C2H2-ZF protein genes in the silkworm, 290 in the monarch butterfly, 243 in the fruit fly, 107 in elegans, 673 in mouse, and 1,082 in human. The C2H2-ZF protein genes of the silkworm were classified into three main grouping clades according to a phylogenetic classification, and 312 of these genes could be mapped onto 27 chromosomes. Most silkworm C2H2-ZF protein genes exhibited specific expression in larval tissues. Furthermore, several C2H2-ZF protein genes had sex-specific expression during metamorphosis. In addition, we found that some C2H2-ZF protein genes are involved in metamorphosis and female reproduction by using expression clustering and gene annotation analysis. Among them, five genes were selected, BGIBMGA002091 (CTCF), BGIBMGA006492 (fru), BGIBMGA006230 (wor), BGIBMGA004640 (lola), and BIGBMGA004569, for quantitative real-time PCR analysis from larvae to adult ovaries. The results showed that the five genes had different expression patterns in ovaries, among which BGIBMGA002091 (CTCF) gene expression level was the highest, and its expression level increased rapidly in late pupae and adult stages. These findings provide a basis for further investigation of the functions of C2H2-ZF protein genes in the silkworm, and the results offer clues for further research into the development of metamorphosis and female reproduction in the silkworm.

Similarity regression predicts evolution of transcription factor sequence specificity

Transcription factor (TF) binding specificities (motifs) are essential for the analysis of gene regulation. Accurate prediction of TF motifs is critical, because it is infeasible to assay all TFs in all sequenced eukaryotic genomes. There is ongoing controversy regarding the degree of motif diversification among related species that is, in part, because of uncertainty in motif prediction methods. Here we describe similarity regression, a significantly improved method for predicting motifs, which we use to update and expand the Cis-BP database. Similarity regression inherently quantifies TF motif evolution, and shows that previous claims of near-complete conservation of motifs between human and Drosophila are inflated, with nearly half of the motifs in each species absent from the other, largely due to extensive divergence in C2H2 zinc finger proteins. We conclude that diversification in DNA-binding motifs is pervasive, and present a new tool and updated resource to study TF diversity and gene regulation across eukaryotes.

The bZIP mutant CEBPB (V285A) has sequence specific DNA binding propensities similar to CREB1

The bZIP homodimers CEBPB and CREB1 bind DNA containing methylated cytosines differently. CREB1 binds stronger to the C/EBP half-site GCAA when the cytosine is methylated. For CEBPB, methylation of the same cytosine does not affect DNA binding. The X-ray structure of CREB1 binding the half site GTCA identifies an alanine in the DNA binding region interacting with the methyl group of T, structurally analogous to the methyl group of methylated C. This alanine is replaced with a valine in CEBPB. To explore the contribution of this amino acid to binding with methylated cytosine of the GCAA half-site, we made the reciprocal mutants CEBPB(V285A) and CREB1(A297V) and used protein binding microarrays (PBM) to examine binding to four types of double-stranded DNA (dsDNA): 1) DNA with cytosine in both strands (DNA(C|C)), 2) DNA with 5-methylcytosine (M) in one strand and cytosine in the second strand (DNA(M|C)), 3) DNA with 5-hydroxymethylcytosine (H) in one strand and cytosine in the second strand (DNA(H|C)), and 4) DNA with both cytosines in all CG dinucleotides containing 5-methylcytosine (DNA(5mCG)). When binding to DNA(C|C), CEBPB (V285A) preferentially binds the CRE consensus motif (TGACGTCA), similar to CREB1. The reciprocal mutant, CREB1(A297V) binds DNA with some similarity to CEBPB, with strongest binding to the methylated PAR site 8-mer TTACGTAA. These data demonstrate that V285 residue inhibits CEBPB binding to methylated cytosine of the GCAA half-site.

Functional Analysis and Fine Mapping of the 9p22.2 Ovarian Cancer Susceptibility Locus

Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. SIGNIFICANCE: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene.See related commentary by Choi and Brown, p. 439.

Genome Editing of Pigs for Agriculture and Biomedicine

Pigs serve as an important agricultural resource and animal model in biomedical studies. Efficient and precise modification of pig genome by using recently developed gene editing tools has significantly broadened the application of pig models in various research areas. The three types of site-specific nucleases, namely, zinc-finger nucleases, transcription activator-like effector nucleases, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein, are the main gene editing tools that can efficiently introduce predetermined modifications, including knockouts and knockins, into the pig genome. These modifications can confer desired phenotypes to pigs to improve production traits, such as optimal meat production, enhanced feed digestibility, and disease resistance. Besides, given their genetic, anatomic, and physiologic similarities to humans, pigs can also be modified to model human diseases or to serve as an organ source for xenotransplantation to save human lives. To date, many genetically modified pig models with agricultural or biomedical values have been established by using gene editing tools. These pig models are expected to accelerate research progress in related fields and benefit humans.

Cys₂/His₂ Zinc-Finger Proteins in Transcriptional Regulation of Flower Development

Flower development is the core of higher-plant ontogenesis and is controlled by complex gene regulatory networks. Cys₂/His₂ zinc-finger proteins (C2H2-ZFPs) constitute one of the largest transcription factor families and are highly involved in transcriptional regulation of flowering induction, floral organ morphogenesis, and pollen and pistil maturation. Nevertheless, the molecular mechanism of C2H2-ZFPs has been gradually revealed only in recent years. During flowering induction, C2H2-ZFPs can modify the chromatin of FLOWERING LOCUS C, thereby providing additional insights into the quantification of transcriptional regulation caused by chromatin regulation. C2H2-ZFPs are involved in cell division and proliferation in floral organ development and are associated with hormonal regulation, thereby revealing how a flower is partitioned into four developmentally distinct whorls. The studies reviewed in this work integrate the information from the endogenous, hormonal, and environmental regulation of flower development. The structure of C2H2-ZFPs determines their function as transcriptional regulators. The findings indicate that C2H2-ZFPs play a crucial role in flower development. In this review, we summarize the current understanding of the structure, expression, and function of C2H2-ZFPs and discuss their molecular mechanism in flower development.

Replacing C189 in the bZIP domain of Zta with S, T, V, or A changes DNA binding specificity to four types of double-stranded DNA

Zta is a bZIP transcription factor (TF) in the Epstein-Barr virus that binds unmethylated and methylated DNA sequences. Substitution of cysteine 189 of Zta to serine (Zta(C189S)) results in a virus that is unable to execute the lytic cycle, which was attributed to a change in binding to methylated DNA sequences. To learn more about the role of this position in defining sequence-specific DNA binding, we mutated cysteine 189 to four other amino acids, producing Zta(C189S), Zta(C189T), Zta(C189A), and Zta(C189V) mutants. Zta and mutants were used in protein binding microarray (PBM) experiments to evaluate sequence-specific DNA binding to four types of double-stranded DNA (dsDNA): 1) with cytosine in both strands (DNA(C|C)), 2) with 5-methylcytosine (5mC) in one strand and cytosine in the second strand (DNA(5mC|C)), 3) with 5-hydroxymethylcytosine (5hmC) in one strand and cytosine in the second strand (DNA(5hmC|C)), and 4) with both cytosines in all CG dinucleotides containing 5mC (DNA(5mCG)). Zta(C189S) and Zta(C189T) bound the TRE (AP-1) motif (TGA^G/_CTCA) more strongly than wild-type Zta, while binding to other sequences, including the C/EBP half site GCAA was reduced. Binding of Zta(C189S) and Zta(C189T) to DNA containing modified cytosines (DNA(5mC|C), DNA(5hmC|C), and DNA(5mCG)) was reduced compared to Zta. Zta(C189A) and Zta(C189V) had higher non-specific binding to all four types of DNA. Our data suggests that position C189 in Zta impacts sequence-specific binding to DNA containing modified and unmodified cytosine.

An incoherent feedforward loop facilitates adaptive tuning of gene expression

We studied adaptive evolution of gene expression using long-term experimental evolution of Saccharomyces cerevisiae in ammonium-limited chemostats. We found repeated selection for non-synonymous variation in the DNA binding domain of the transcriptional activator, GAT1, which functions with the repressor, DAL80 in an incoherent type-1 feedforward loop (I1-FFL) to control expression of the high affinity ammonium transporter gene, MEP2. Missense mutations in the DNA binding domain of GAT1 reduce its binding to the GATAA consensus sequence. However, we show experimentally, and using mathematical modeling, that decreases in GAT1 binding result in increased expression of MEP2 as a consequence of properties of I1-FFLs. Our results show that I1-FFLs, one of the most commonly occurring network motifs in transcriptional networks, can facilitate adaptive tuning of gene expression through modulation of transcription factor binding affinities. Our findings highlight the importance of gene regulatory architectures in the evolution of gene expression.

Structural basis for the ability of MBD domains to bind methyl-CG and TG sites in DNA

Cytosine methylation is a well-characterized epigenetic mark and occurs at both CG and non-CG sites in DNA. Both methylated CG (mCG)- and mCH (H = A, C, or T)-containing DNAs, especially mCAC-containing DNAs, are recognized by methyl-CpG-binding protein 2 (MeCP2) to regulate gene expression in neuron development. However, the molecular mechanism involved in the binding of methyl-CpG-binding domain (MBD) of MeCP2 to these different DNA motifs is unclear. Here, we systematically characterized the DNA-binding selectivities of the MBD domains in MeCP2 and MBD1-4 with isothermal titration calorimetry-based binding assays, mutagenesis studies, and X-ray crystallography. We found that the MBD domains of MeCP2 and MBD1-4 bind mCG-containing DNAs independently of the sequence identity outside the mCG dinucleotide. Moreover, some MBD domains bound to both methylated and unmethylated CA dinucleotide-containing DNAs, with a preference for the CAC sequence motif. We also found that the MBD domains bind to mCA or nonmethylated CA DNA by recognizing the complementary TG dinucleotide, which is consistent with an overlooked ligand of MeCP2, i.e. the matrix/scaffold attachment regions (MARs/SARs) with a consensus sequence of 5'-GGTGT-3' that was identified in early 1990s. Our results also explain why MeCP2 exhibits similar binding affinity to both mCA- and hmCA-containing dsDNAs. In summary, our results suggest that in addition to mCG sites, unmethylated CA or TG sites also serve as DNA-binding sites for MeCP2 and other MBD-containing proteins. This discovery expands the genome-wide activity of MBD-containing proteins in gene regulation.

The intracellular immune receptor Rx1 regulates the DNA-binding activity of a Golden2-like transcription factor

Plant nucleotide-binding leucine-rich repeat (NLR) proteins enable the immune system to recognize and respond to pathogen attack. An early consequence of immune activation is transcriptional reprogramming, and some NLRs have been shown to act in the nucleus and interact with transcription factors. The Rx1 NLR protein of potato is further able to bind and distort double-stranded DNA. However, Rx1 host targets that support a role for Rx1 in transcriptional reprogramming at DNA are unknown. Here, we report a functional interaction between Rx1 and NbGlk1, a Golden2-like transcription factor. Rx1 binds to NbGlk1 in vitro and in planta. NbGlk1 binds to known Golden2-like consensus DNA sequences. Rx1 reduces the binding affinity of NbGlk1 for DNA in vitro. NbGlk1 activates cellular responses to potato virus X, whereas Rx1 associates with NbGlk1 and prevents its assembly on DNA in planta unless activated by PVX. This study provides new mechanistic insight into how an NLR can coordinate an immune signaling response at DNA following pathogen perceptions.

The Epstein-Barr Virus B-ZIP Protein Zta Recognizes Specific DNA Sequences Containing 5-Methylcytosine and 5-Hydroxymethylcytosine

The Epstein-Barr virus (EBV) B-ZIP transcription factor Zta binds to many DNA sequences containing methylated CG dinucleotides. Using protein binding microarrays (PBMs), we analyzed the sequence specific DNA binding of Zta to four kinds of double-stranded DNA (dsDNA): (1) DNA containing cytosine in both strands, (2) DNA with 5-methylcytosine (5mC) in one strand and cytosine in the second strand, (3) DNA with 5-hydroxymethylcytosine (5hmC) in one strand and cytosine in the second strand, and (4) DNA in which both cytosines in all CG dinucleotides contain 5mC. We compared these data to PBM data for three additional B-ZIP proteins (CREB1 and CEBPB homodimers and cJun|cFos heterodimers). With cytosine, Zta binds the TRE motif TGA^C/_GTCA as previously reported. With CG dinucleotides containing 5mC on both strands, many TRE motif variants containing a methylated CG dinucleotide at two positions in the motif, such as MGAGTCA and TGAGMGA (where M = 5mC), were preferentially bound. 5mC inhibits binding of Zta to both TRE motif half-sites GTCA and CTCA. Like the CREB1 homodimer, the Zta homodimer and the cJun|cFos heterodimer more strongly bind the C/EBP half-site tetranucleotide GCAA when it contains 5mC. Zta also binds dsDNA sequences containing 5hmC in one strand, although the effect is less dramatic than that observed for 5mC. Our results identify new DNA sequences that are well-bound by the viral B-ZIP protein Zta only when they contain 5mC or 5hmC, uncovering the potential for discovery of new viral and host regulatory programs controlled by EBV.

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.

Characteristics of Genome Editing Mutations in Cereal Crops

Designer nucleases allow the creation of new plant genotypes by introducing precisely-targeted double-strand breaks that are resolved by endogenous repair pathways. The major nuclease technologies are meganucleases, zinc-finger nucleases, transcription activator-like effector nucleases, and the CRISPR/Cas9 system. Each comprises a promiscuous endonuclease guided by protein-DNA or RNA-DNA interactions. A great deal is known about the principles of designer nucleases but much remains to be learned about their detailed behavioral characteristics in different plant species. The outcome of genome engineering reflects the intrinsic properties of each nuclease and target genome, causing variations in efficiency, accuracy, and mutation structure. In this article, we critically discuss the activities of designer nucleases in different cereals representing a broad range of genome characteristics.

5-Methylcytosine (5mC) and 5-Hydroxymethylcytosine (5hmC) Enhance the DNA Binding of CREB1 to the C/EBP Half-Site Tetranucleotide GCAA

In human and mouse stem cells and brain, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) can occur outside of CG dinucleotides. Using protein binding microarrays (PBMs) containing 60-mer DNA probes, we evaluated the effect of 5mC and 5hmC on one DNA strand on the double-stranded DNA binding of the mouse B-ZIP transcription factors (TFs) CREB1, ATF1, and JUND. 5mC inhibited binding of CREB1 to the canonical CRE half-site |GTCA but enhanced binding to the C/EBP half-site |GCAA. 5hmC inhibited binding of CREB1 to all 8-mers except TGAT|GCAA, where binding is enhanced. We observed similar DNA binding patterns with ATF1, a closely related B-ZIP domain. In contrast, both 5mC and 5hmC inhibited binding of JUND. These results identify new DNA sequences that are well-bound by CREB1 and ATF1 only when they contain 5mC or 5hmC. Analysis of two X-ray structures examines the consequences of 5mC and 5hmC on DNA binding by CREB and FOS|JUN.

Transcription Factors Encoded on Core and Accessory Chromosomes of Fusarium oxysporum Induce Expression of Effector Genes

Proteins secreted by pathogens during host colonization largely determine the outcome of pathogen-host interactions and are commonly called 'effectors'. In fungal plant pathogens, coordinated transcriptional up-regulation of effector genes is a key feature of pathogenesis and effectors are often encoded in genomic regions with distinct repeat content, histone code and rate of evolution. In the tomato pathogen Fusarium oxysporum f. sp. lycopersici (Fol), effector genes reside on one of four accessory chromosomes, known as the 'pathogenicity' chromosome, which can be exchanged between strains through horizontal transfer. The three other accessory chromosomes in the Fol reference strain may also be important for virulence towards tomato. Expression of effector genes in Fol is highly up-regulated upon infection and requires Sge1, a transcription factor encoded on the core genome. Interestingly, the pathogenicity chromosome itself contains 13 predicted transcription factor genes and for all except one, there is a homolog on the core genome. We determined DNA binding specificity for nine transcription factors using oligonucleotide arrays. The binding sites for homologous transcription factors were highly similar, suggesting that extensive neofunctionalization of DNA binding specificity has not occurred. Several DNA binding sites are enriched on accessory chromosomes, and expression of FTF1, its core homolog FTF2 and SGE1 from a constitutive promoter can induce expression of effector genes. The DNA binding sites of only these three transcription factors are enriched among genes up-regulated during infection. We further show that Ftf1, Ftf2 and Sge1 can activate transcription from their binding sites in yeast. RNAseq analysis revealed that in strains with constitutive expression of FTF1, FTF2 or SGE1, expression of a similar set of plant-responsive genes on the pathogenicity chromosome is induced, including most effector genes. We conclude that the Fol pathogenicity chromosome may be partially transcriptionally autonomous, but there are also extensive transcriptional connections between core and accessory chromosomes.

G9a and ZNF644 Physically Associate to Suppress Progenitor Gene Expression during Neurogenesis

Proliferating progenitor cells undergo changes in competence to give rise to post-mitotic progeny of specialized function. These cell-fate transitions typically involve dynamic regulation of gene expression by histone methyltransferase (HMT) complexes. However, the composition, roles, and regulation of these assemblies in regulating cell-fate decisions in vivo are poorly understood. Using unbiased affinity purification and mass spectrometry, we identified the uncharacterized C2H2-like zinc finger protein ZNF644 as a G9a/GLP-interacting protein and co-regulator of histone methylation. In zebrafish, functional characterization of ZNF644 orthologs, znf644a and znf644b, revealed complementary roles in regulating G9a/H3K9me2-mediated gene silencing during neurogenesis. The non-overlapping requirements for znf644a and znf644b during retinal differentiation demarcate critical aspects of retinal differentiation programs regulated by differential G9a-ZNF644 associations, such as transitioning proliferating progenitor cells toward differentiation. Collectively, our data point to ZNF644 as a critical co-regulator of G9a/H3K9me2-mediated gene silencing during neuronal differentiation.

Gene editing and its application for hematological diseases

The use of precise, rationally designed gene-editing nucleases allows for targeted genome and transcriptome modification, and at present, four major classes of nucleases are being employed: zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), meganucleases (MNs), and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9. Each reagent shares the ability to recognize and bind a target sequence of DNA. Depending on the properties of the reagent, the DNA can be cleaved on one or both strands, or epigenetic changes can be mediated. These novel properties can impact hematological disease by allowing for: (1) direct modification of hematopoietic stem/progenitor cells (HSPCs), (2) gene alteration of hematopoietic lineage committed terminal effectors, (3) genome engineering in non-hematopoietic cells with reprogramming to a hematopoietic phenotype, and (4) transcriptome modulation for gene regulation, modeling, and discovery.

Artificial zinc finger DNA binding domains: versatile tools for genome engineering and modulation of gene expression

Genome editing and alteration of gene expression by synthetic DNA binding activities gained a lot of momentum over the last decade. This is due to the development of new DNA binding molecules with enhanced binding specificity. The most commonly used DNA binding modules are zinc fingers (ZFs), TALE-domains, and the RNA component of the CRISPR/Cas9 system. These binding modules are fused or linked to either nucleases that cut the DNA and induce DNA repair processes, or to protein domains that activate or repress transcription of genes close to the targeted site in the genome. This review focuses on the structure, design, and applications of ZF DNA binding domains (ZFDBDs). ZFDBDs are relatively small and have been shown to penetrate the cell membrane without additional tags suggesting that they could be delivered to cells without a DNA or RNA intermediate. Advanced algorithms that are based on extensive knowledge of the mode of ZF/DNA interactions are used to design the amino acid composition of ZFDBDs so that they bind to unique sites in the genome. Off-target binding has been a concern for all synthetic DNA binding molecules. Thus, increasing the specificity and affinity of ZFDBDs will have a significant impact on their use in analytical or therapeutic settings.

Moonlighting transcriptional activation function of a fungal sulfur metabolism enzyme

Moonlighting proteins, including metabolic enzymes acting as transcription factors (TF), are present in a variety of organisms but have not been described in higher fungi so far. In a previous genome-wide analysis of the TF repertoire of the plant-symbiotic fungus Tuber melanosporum, we identified various enzymes, including the sulfur-assimilation enzyme phosphoadenosine-phosphosulfate reductase (PAPS-red), as potential transcriptional activators. A functional analysis performed in the yeast Saccharomyces cerevisiae, now demonstrates that a specific variant of this enzyme, PAPS-red A, localizes to the nucleus and is capable of transcriptional activation. TF moonlighting, which is not present in the other enzyme variant (PAPS-red B) encoded by the T. melanosporum genome, relies on a transplantable C-terminal polypeptide containing an alternating hydrophobic/hydrophilic amino acid motif. A similar moonlighting activity was demonstrated for six additional proteins, suggesting that multitasking is a relatively frequent event. PAPS-red A is sulfur-state-responsive and highly expressed, especially in fruitbodies, and likely acts as a recruiter of transcription components involved in S-metabolism gene network activation. PAPS-red B, instead, is expressed at low levels and localizes to a highly methylated and silenced region of the genome, hinting at an evolutionary mechanism based on gene duplication, followed by epigenetic silencing of this non-moonlighting gene variant.

The CRISPR/Cas9 system for gene editing and its potential application in pain research

The CRISPR/Cas9 system is a research hotspot in genome editing and regulation. Currently, it is used in genomic silencing and knock-in experiments as well as transcriptional activation and repression. This versatile system consists of two components: a guide RNA (gRNA) and a Cas9 nuclease. Recognition of a genomic DNA target is mediated through base pairing with a 20-base gRNA. The latter further recruits the Cas9 endonuclease protein to the target site and creates double-stranded breaks in the target DNA. Compared with traditional genome editing directed by DNA-binding protein domains, this short RNA-directed Cas9 endonuclease system is simple and easily programmable. Although this system may have off-target effects and in vivo delivery and immune challenges, researchers have employed this system in vivo to establish disease models, study specific gene functions under certain disease conditions, and correct genomic information for disease treatment. In regards to pain research, the CRISPR/Cas9 system may act as a novel tool in gene correction therapy for pain-associated hereditary diseases and may be a new approach for RNA-guided transcriptional activation or repression of pain-related genes. In addition, this system is also applied to loss-of-function mutations in pain-related genes and knockin of reporter genes or loxP tags at pain-related genomic loci. The CRISPR/Cas9 system will likely be carried out widely in both bench work and clinical settings in the pain field.

Integrating gene synthesis and microfluidic protein analysis for rapid protein engineering

The capability to rapidly design proteins with novel functions will have a significant impact on medicine, biotechnology and synthetic biology. Synthetic genes are becoming a commodity, but integrated approaches have yet to be developed that take full advantage of gene synthesis. We developed a solid-phase gene synthesis method based on asymmetric primer extension (APE) and coupled this process directly to high-throughput, on-chip protein expression, purification and characterization (via mechanically induced trapping of molecular interactions, MITOMI). By completely circumventing molecular cloning and cell-based steps, APE-MITOMI reduces the time between protein design and quantitative characterization to 3–4 days. With APE-MITOMI we synthesized and characterized over 400 zinc-finger (ZF) transcription factors (TF), showing that although ZF TFs can be readily engineered to recognize a particular DNA sequence, engineering the precise binding energy landscape remains challenging. We also found that it is possible to engineer ZF–DNA affinity precisely and independently of sequence specificity and that in silico modeling can explain some of the observed affinity differences. APE-MITOMI is a generic approach that should facilitate fundamental studies in protein biophysics, and protein design/engineering.

Motif co-regulation and co-operativity are common mechanisms in transcriptional, post-transcriptional and post-translational regulation

A substantial portion of the regulatory interactions in the higher eukaryotic cell are mediated by simple sequence motifs in the regulatory segments of genes and (pre-)mRNAs, and in the intrinsically disordered regions of proteins. Although these regulatory modules are physicochemically distinct, they share an evolutionary plasticity that has facilitated a rapid growth of their use and resulted in their ubiquity in complex organisms. The ease of motif acquisition simplifies access to basal housekeeping functions, facilitates the co-regulation of multiple biomolecules allowing them to respond in a coordinated manner to changes in the cell state, and supports the integration of multiple signals for combinatorial decision-making. Consequently, motifs are indispensable for temporal, spatial, conditional and basal regulation at the transcriptional, post-transcriptional and post-translational level. In this review, we highlight that many of the key regulatory pathways of the cell are recruited by motifs and that the ease of motif acquisition has resulted in large networks of co-regulated biomolecules. We discuss how co-operativity allows simple static motifs to perform the conditional regulation that underlies decision-making in higher eukaryotic biological systems. We observe that each gene and its products have a unique set of DNA, RNA or protein motifs that encode a regulatory program to define the logical circuitry that guides the life cycle of these biomolecules, from transcription to degradation. Finally, we contrast the regulatory properties of protein motifs and the regulatory elements of DNA and (pre-)mRNAs, advocating that co-regulation, co-operativity, and motif-driven regulatory programs are common mechanisms that emerge from the use of simple, evolutionarily plastic regulatory modules.

Engineering large animal models of human disease

The recent development of gene editing tools and methodology for use in livestock enables the production of new animal disease models. These tools facilitate site‐specific mutation of the genome, allowing animals carrying known human disease mutations to be produced. In this review, we describe the various gene editing tools and how they can be used for a range of large animal models of diseases. This genomic technology is in its infancy but the expectation is that through the use of gene editing tools we will see a dramatic increase in animal model resources available for both the study of human disease and the translation of this knowledge into the clinic. Comparative pathology will be central to the productive use of these animal models and the successful translation of new therapeutic strategies. © 2015 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Can genome engineering be used to target cancer-associated enhancers?

Transcriptional misregulation is involved in the development of many diseases, especially neoplastic transformation. Distal regulatory elements, such as enhancers, play a major role in specifying cell-specific transcription patterns in both normal and diseased tissues, suggesting that enhancers may be prime targets for therapeutic intervention. By focusing on modulating gene regulation mediated by cell type-specific enhancers, there is hope that normal epigenetic patterning in an affected tissue could be restored with fewer side effects than observed with treatments employing relatively nonspecific inhibitors such as epigenetic drugs. New methods employing genomic nucleases and site-specific epigenetic regulators targeted to specific genomic regions, using either artificial DNA-binding proteins or RNA-DNA interactions, may allow precise genome engineering at enhancers. However, this field is still in its infancy and further refinements that increase specificity and efficiency are clearly required.

A Novel AT-Rich DNA Recognition Mechanism for Bacterial Xenogeneic Silencer MvaT

Bacterial xenogeneic silencing proteins selectively bind to and silence expression from many AT rich regions of the chromosome. They serve as master regulators of horizontally acquired DNA, including a large number of virulence genes. To date, three distinct families of xenogeneic silencers have been identified: H-NS of Proteobacteria, Lsr2 of the Actinomycetes, and MvaT of Pseudomonas sp. Although H-NS and Lsr2 family proteins are structurally different, they all recognize the AT-rich DNA minor groove through a common AT-hook-like motif, which is absent in the MvaT family. Thus, the DNA binding mechanism of MvaT has not been determined. Here, we report the characteristics of DNA sequences targeted by MvaT with protein binding microarrays, which indicates that MvaT prefers binding flexible DNA sequences with multiple TpA steps. We demonstrate that there are clear differences in sequence preferences between MvaT and the other two xenogeneic silencer families. We also determined the structure of the DNA-binding domain of MvaT in complex with a high affinity DNA dodecamer using solution NMR. This is the first experimental structure of a xenogeneic silencer in complex with DNA, which reveals that MvaT recognizes the AT-rich DNA both through base readout by an "AT-pincer" motif inserted into the minor groove and through shape readout by multiple lysine side chains interacting with the DNA sugar-phosphate backbone. Mutations of key MvaT residues for DNA binding confirm their importance with both in vitro and in vivo assays. This novel DNA binding mode enables MvaT to better tolerate GC-base pair interruptions in the binding site and less prefer A tract DNA when compared to H-NS and Lsr2. Comparison of MvaT with other bacterial xenogeneic silencers provides a clear picture that nature has evolved unique solutions for different bacterial genera to distinguish foreign from self DNA.

Mapping and analysis of Caenorhabditis elegans transcription factor sequence specificities

Caenorhabditis elegans is a powerful model for studying gene regulation, as it has a compact genome and a wealth of genomic tools. However, identification of regulatory elements has been limited, as DNA-binding motifs are known for only 71 of the estimated 763 sequence-specific transcription factors (TFs). To address this problem, we performed protein binding microarray experiments on representatives of canonical TF families in C. elegans, obtaining motifs for 129 TFs. Additionally, we predict motifs for many TFs that have DNA-binding domains similar to those already characterized, increasing coverage of binding specificities to 292 C. elegans TFs (∼40%). These data highlight the diversification of binding motifs for the nuclear hormone receptor and C2H2 zinc finger families and reveal unexpected diversity of motifs for T-box and DM families. Motif enrichment in promoters of functionally related genes is consistent with known biology and also identifies putative regulatory roles for unstudied TFs.

DOI:http://dx.doi.org/10.7554/eLife.06967.001

Many scientists use ‘model’ species—such as the fruit fly or a nematode worm called Caenorhabditis elegans—in their research because these organisms have useful features that make it easier to carry out many experiments. For example, C. elegans has a smaller genome compared to many other animals, which is useful for studying the roles of individual genes or stretches of DNA.

Transcription factors are a type of protein that can bind to specific stretches of DNA and help to switch certain genes on or off. These ‘motifs’ may be close to the gene or further away in the genome, and therefore, must stand out amongst the rest of the DNA, like lights on a landing strip. However, the motifs for only 10% of the estimated 763 transcription factors in C. elegans have been identified so far.

In this study, Narasimhan, Lambert, Yang et al. used a technique called a ‘protein binding microarray’ to identify the motifs for many more of the C. elegans transcription factors. These findings were then used to predict motifs for other transcription factors. Together, these methods increased the proportion of C. elegans transcription factors with known DNA-binding motifs from 10% to around 40%.

Now that more DNA motifs have been identified, it is possible to look for similarities and differences between them. For example, Narasimhan, Lambert, Yang et al. found that transcription factors with similar sequences can bind to very varied motifs. On the other hand, some transcription factors that are very different are able to recognize very similar motifs. The experiments also indicate that motifs found very close to genes—in sequences known as ‘promoters’—may be able to interact with many proteins to influence the activity of genes.

Narasimhan, Lambert, Yang et al.'s findings increase the number of C. elegans transcription factors with a motif, bringing the knowledge of these proteins more in line with the better-studied transcription factors of humans and fruit flies. The next challenge is to identify DNA motifs for the remaining 60% of transcription factors.

DOI:http://dx.doi.org/10.7554/eLife.06967.002

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.

A systematic survey of the Cys2His2 zinc finger DNA-binding landscape

Cys2His2 zinc fingers (C2H2-ZFs) comprise the largest class of metazoan DNA-binding domains. Despite this domain's well-defined DNA-recognition interface, and its successful use in the design of chimeric proteins capable of targeting genomic regions of interest, much remains unknown about its DNA-binding landscape. To help bridge this gap in fundamental knowledge and to provide a resource for design-oriented applications, we screened large synthetic protein libraries to select binding C2H2-ZF domains for each possible three base pair target. The resulting data consist of >160 000 unique domain-DNA interactions and comprise the most comprehensive investigation of C2H2-ZF DNA-binding interactions to date. An integrated analysis of these independent screens yielded DNA-binding profiles for tens of thousands of domains and led to the successful design and prediction of C2H2-ZF DNA-binding specificities. Computational analyses uncovered important aspects of C2H2-ZF domain-DNA interactions, including the roles of within-finger context and domain position on base recognition. We observed the existence of numerous distinct binding strategies for each possible three base pair target and an apparent balance between affinity and specificity of binding. In sum, our comprehensive data help elucidate the complex binding landscape of C2H2-ZF domains and provide a foundation for efforts to determine, predict and engineer their DNA-binding specificities.

Recent advances in developing molecular tools for targeted genome engineering of mammalian cells

Various biological molecules naturally existing in diversified species including fungi, bacteria, and bacteriophage have functionalities for DNA binding and processing. The biological molecules have been recently actively engineered for use in customized genome editing of mammalian cells as the molecule-encoding DNA sequence information and the underlying mechanisms how the molecules work are unveiled. Excitingly, multiple novel methods based on the newly constructed artificial molecular tools have enabled modifications of specific endogenous genetic elements in the genome context at efficiencies that are much higher than that of the conventional homologous recombination based methods. This minireview introduces the most recently spotlighted molecular genome engineering tools with their key features and ongoing modifications for better performance. Such ongoing efforts have mainly focused on the removal of the inherent DNA sequence recognition rigidity from the original molecular platforms, the addition of newly tailored targeting functions into the engineered molecules, and the enhancement of their targeting specificity. Effective targeted genome engineering of mammalian cells will enable not only sophisticated genetic studies in the context of the genome, but also widely-applicable universal therapeutics based on the pinpointing and correction of the disease-causing genetic elements within the genome in the near future.

Engineering the control of mosquito-borne infectious diseases

Recent advances in genetic engineering are bringing new promise for controlling mosquito populations that transmit deadly pathogens. Here we discuss past and current efforts to engineer mosquito strains that are refractory to disease transmission or are suitable for suppressing wild disease-transmitting populations.

Determination and inference of eukaryotic transcription factor sequence specificity

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

Rule-based design of synthetic transcription factors in eukaryotes

To design and build living systems, synthetic biologists have at their disposal an increasingly large library of naturally derived and synthetic parts. These parts must be combined together in particular orders, orientations, and spacings to achieve desired functionalities. These structural constraints can be viewed as grammatical rules describing how to assemble parts together into larger functional units. Here, we develop a grammar for the design of synthetic transcription factors (sTFs) in eukaryotic cells and implement it within GenoCAD, a Computer-Aided Design (CAD) software for synthetic biology. Knowledge derived from experimental evidence was captured in this grammar to guide the user to create designer transcription factors that should operate as intended. The grammar can be easily updated and refined as our experience with using sTFs in different contexts increases. In combination with grammars that define other synthetic systems, we anticipate that this work will enable the more reliable, efficient, and automated design of synthetic cells with rich functionalities.

Role of stem cells in large animal genetic engineering in the TALENs-CRISPR era

The establishment of embryonic stem cells (ESCs) and gene targeting technologies in mice has revolutionised the field of genetics. The relative ease with which genes can be knocked out, and exogenous sequences introduced, has allowed the mouse to become the prime model for deciphering the genetic code. Not surprisingly, the lack of authentic ESCs has hampered the livestock genetics field and has forced animal scientists into adapting alternative technologies for genetic engineering. The recent discovery of the creation of induced pluripotent stem cells (iPSCs) by upregulation of a handful of reprogramming genes has offered renewed enthusiasm to animal geneticists. However, much like ESCs, establishing authentic iPSCs from the domestic animals is still beset with problems, including (but not limited to) the persistent expression of reprogramming genes and the lack of proven potential for differentiation into target cell types both in vitro and in vivo. Site-specific nucleases comprised of zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs) and clustered regulated interspaced short palindromic repeats (CRISPRs) emerged as powerful genetic tools for precisely editing the genome, usurping the need for ESC-based genetic modifications even in the mouse. In this article, in the aftermath of these powerful genome editing technologies, the role of pluripotent stem cells in livestock genetics is discussed.