Transcription factors: specific DNA binding and specific gene regulation

Receptor-activated transcription factors and beyond: multiple modes of Smad2/3-dependent transmission of TGF-β signaling

Transforming growth factor β (TGF-β) is a pleiotropic cytokine that is widely distributed throughout the body. Its receptor proteins, TGF-β type I and type II receptors, are also ubiquitously expressed. Therefore, the regulation of various signaling outputs in a context-dependent manner is a critical issue in this field. Smad proteins were originally identified as signal-activated transcription factors similar to signal transducer and activator of transcription proteins. Smads are activated by serine phosphorylation mediated by intrinsic receptor dual specificity kinases of the TGF-β family, indicating that Smads are receptor-restricted effector molecules downstream of ligands of the TGF-β family. Smad proteins have other functions in addition to transcriptional regulation, including post-transcriptional regulation of micro-RNA processing, pre-mRNA splicing, and m6A methylation. Recent technical advances have identified a novel landscape of Smad-dependent signal transduction, including regulation of mitochondrial function without involving regulation of gene expression. Therefore, Smad proteins are receptor-activated transcription factors and also act as intracellular signaling modulators with multiple modes of function. In this review, we discuss the role of Smad proteins as receptor-activated transcription factors and beyond. We also describe the functional differences between Smad2 and Smad3, two receptor-activated Smad proteins downstream of TGF-β, activin, myostatin, growth and differentiation factor (GDF) 11, and Nodal.

Transcription factor FOXM1 specifies chromatin DNA to extracellular vesicles

Extracellular vesicle DNAs (evDNAs) hold significant diagnostic value for various diseases and facilitate transcellular transfer of genetic material. Our study identifies transcription factor FOXM1 as a mediator for directing chromatin genes or DNA fragments (termed FOXM1-chDNAs) to extracellular vesicles (EVs). FOXM1 binds to MAP1LC3/LC3 in the nucleus, and FOXM1-chDNAs, such as the DUX4 gene and telomere DNA, are designated by FOXM1 binding and translocated to the cytoplasm before being released to EVs through the secretory autophagy during lysosome inhibition (SALI) process involving LC3. Disrupting FOXM1 expression or the SALI process impairs FOXM1-chDNAs incorporation into EVs. FOXM1-chDNAs can be transmitted to recipient cells via EVs and expressed in recipient cells when they carry functional genes. This finding provides an example of how chromatin DNA fragments are specified to EVs by transcription factor FOXM1, revealing its contribution to the formation of evDNAs from nuclear chromatin. It provides a basis for further exploration of the roles of evDNAs in biological processes, such as horizontal gene transfer.Abbreviation: ATG5: autophagy related 5; CCFs: cytoplasmic chromatin fragments; ChIP: chromatin immunoprecipitation; cytoDNA: cytoplasmic DNA; CQ: chloroquine; FOXM1-DBD: FOXM1 DNA binding domain; DUX4:double homeobox 4; EVs: extracellular vesicles; evDNAs: extracellular vesicle DNAs; FOXM1: forkhead box M1; FOXM1-chDNAs: chromatin DNA fragments directed by FOXM1 to EVs; HGT: horizontal gene transfer; LC3-II: lipid modified LC3; LMNB1: lamin B1; LIR: LC3-interacting region; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; MVBs: multivesicular bodies; M1-binding DNA: a linear DNA containing 72× FOXM1 binding sites; SALI: secretory autophagy during lysosome inhibition; siRNA: small interfering RNA; TetO-DUX4: TetO array-containing DUX4 DNA; TetO: tet operator; TetR: tet repressor.

Functional characterisation of Dof gene family and expression analysis under abiotic stresses and melatonin-mediated tolerance in pitaya (Selenicereus undatus)

DNA binding proteins with one finger (Dof ) transcription factors are essential for seed development and defence against various biotic and abiotic stresses in plants. Genomic analysis of Dof has not been determined yet in pitaya (Selenicereus undatus ). In this study, we have identified 26 Dof gene family members, renamed as HuDof-1 to HuDof-26 , and clustered them into seven subfamilies based on conserved motifs, domains, and phylogenetic analysis. The gene pairs of Dof family members were duplicated by segmental duplications that faced purifying selection, as indicated by the K a /K s ratio values. Promoter regions of HuDof genes contain many cis -acting elements related to phytohormones including abscisic acid, jasmonic acid, gibberellin, temperature, and light. We exposed pitaya plants to different environmental stresses and examined melatonin's influence on Dof gene expression levels. Signifcant expression of HuDof -2 and HuDof -6 were observed in different developmental stages of flower buds, flowers, pericarp, and pulp. Pitaya plants were subjected to abiotic stresses, and transcriptome analysis was carried out to identify the role of Dof gene family members. RNA-sequencing data and reverse transcription quantitative PCR-based expression analysis revealed three putative candidate genes (HuDof -1, HuDof -2, and HuDof -8), which might have diverse roles against the abiotic stresses. Our study provides a theoretical foundation for functional analysis through traditional and modern biotechnological tools for pitaya trait improvement.

Utilizing peptide-anchored DNA templates for novel programmable nanoparticle assemblies in biological macromolecules: A review

Biological macromolecules such as proteins and DNA are known to self-assemble into various structural moieties with distinct functions. While nucleic acids are the structural building blocks, peptides exemplify diversity as tailorable biochemical units. Thus, combining the scaffold properties of the biomacromolecule DNA and the functionality of peptides could evolve into a powerful method to obtain tailorable nano assemblies. In this review, we discuss the assembly of non-DNA-coated colloidal NPs on DNA/peptide templates using functional anchors. We begin with strategies for directly attaching metallic NPs to DNA templates to ascertain the functional role of DNA as a scaffold. Followed by methods to assemble peptides onto DNA templates to emphasize the functional versatility of biologically abundant DNA-binding peptides. Next, we focus on studies corroborating peptide self-assembling into macromolecular templates onto which NPs can attach to emphasize the properties of NP-binding peptides. Finally, we discuss the assembly of NPs on a DNA template with a focus on the bifunctional DNA-binding peptides with NP-binding affinity (peptide anchors). This review aims to highlight the immense potential of combining the functional power of DNA scaffolds and tailorable functionalities of peptides for NP assembly and the need to utilize them effectively to obtain tailorable hierarchical NP assemblies.

The forkhead DNA-binding domain binds specific G2-rich RNA sequences

Transcription factors contain a DNA-binding domain ensuring specific recognition of DNA target sequences. The family of forkhead (FOX) transcription factors is composed of dozens of paralogs in mammals. The forkhead domain (FHD) is a segment of about 100 amino acids that binds an A-rich DNA sequence. Using DNA and RNA PCR-SELEX, we show that recombinant FOXL2 proteins, either wild-type or carrying the oncogenic variant C134W, recognize similar DNA-binding sites. This suggests that the oncogenic variant does not alter the intrinsic sequence-specificity of FOXL2. Most importantly, we show that FOXL2 binds G2-rich RNA sequences whereas it virtually fails to bind similar sequences in DNA chemistry. Interestingly, a statistically significant subset of genes responding to the knock-down of FOXL2/Foxl2 harbor such G2-rich sequences and are involved in crucial signaling pathways and cellular processes. In addition, we show that FOXA1, FOXO3a and chimeric FOXL2 proteins containing the FHD of the former are also able to interact with some of the preferred FOXL2-binding sequences. Our results point to an unexpected and novel characteristic of the forkhead domain, the biological relevance of which remains to be explored.

International Union of Basic and Clinical Pharmacology CXIII: Nuclear Receptor Superfamily-Update 2023

The NR superfamily comprises 48 transcription factors in humans that control a plethora of gene network programs involved in a wide range of physiologic processes. This review will summarize and discuss recent progress in NR biology and drug development derived from integrating various approaches, including biophysical techniques, structural studies, and translational investigation. We also highlight how defective NR signaling results in various diseases and disorders and how NRs can be targeted for therapeutic intervention via modulation via binding to synthetic lipophilic ligands. Furthermore, we also review recent studies that improved our understanding of NR structure and signaling. SIGNIFICANCE STATEMENT: Nuclear receptors (NRs) are ligand-regulated transcription factors that are critical regulators of myriad physiological processes. NRs serve as receptors for an array of drugs, and in this review, we provide an update on recent research into the roles of these drug targets.

KaScape: a sequencing-based method for global characterization of protein‒DNA binding affinity

It is difficult to exhaustively screen all possible DNA binding sequences for a given transcription factor (TF). Here, we developed the KaScape method, in which TFs bind to all possible DNA sequences in the same DNA pool where DNA sequences are prepared by randomized oligo synthesis and the random length can be adjusted to a length such as 4, 5, 6, or 7. After separating bound from unbound double-stranded DNAs (dsDNAs), their sequences are determined by next-generation sequencing. To demonstrate the relative binding affinities of all possible DNA sequences determined by KaScape, we developed three-dimensional KaScape viewing software based on a K-mer graph. We applied KaScape to 12 plant TF family AtWRKY proteins and found that all AtWRKY proteins bound to the core sequence GAC with similar profiles. KaScape can detect not only binding sequences consistent with the consensus W-box "TTGAC(C/T)" but also other sequences with weak affinity. KaScape provides a high-throughput, easy-to-operate, sensitive, and exhaustive method for quantitatively characterizing the relative binding strength of a TF with all possible binding sequences, allowing us to comprehensively characterize the specificity and affinity landscape of transcription factors, particularly for moderate- and low-affinity binding sites.

Cooperative assembly confers regulatory specificity and long-term genetic circuit stability

A ubiquitous feature of eukaryotic transcriptional regulation is cooperative self-assembly between transcription factors (TFs) and DNA cis-regulatory motifs. It is thought that this strategy enables specific regulatory connections to be formed in gene networks between otherwise weakly interacting, low-specificity molecular components. Here, using synthetic gene circuits constructed in yeast, we find that high regulatory specificity can emerge from cooperative, multivalent interactions among artificial zinc-finger-based TFs. We show that circuits "wired" using the strategy of cooperative TF assembly are effectively insulated from aberrant misregulation of the host cell genome. As we demonstrate in experiments and mathematical models, this mechanism is sufficient to rescue circuit-driven fitness defects, resulting in genetic and functional stability of circuits in long-term continuous culture. Our naturally inspired approach offers a simple, generalizable means for building high-fidelity, evolutionarily robust gene circuits that can be scaled to a wide range of host organisms and applications.

Acute impact of salinity and C/N ratio on the formation and properties of soluble microbial products from activated sludge

The present study investigated the shock of NaCl and C/N ratio on properties of soluble microbial products (SMPs), focusing on their sized fractions. The results indicated that the NaCl stress increased the content of biopolymers, humic substances, building blocks, and LMW substances in SMPs, while the addition of 40 g NaCl L^-1 significantly changed their relative abundance in SMPs. The acute impact of both N-rich and N-deficient conditions accelerated the secretion of SMPs, but the characteristics of LMW substances differed. Meanwhile, the bio-utilization of SMPs has been enhanced with the increase of NaCl dosage but decreased with the increase of the C/N ratio. The mass balance of sized fractions in SMPs + EPS could be set up when NaCl dosage <10 g/L and C/N ratio >5, which indicates the hydrolysis of sized fractions in EPS mainly compensated for their increase/reduction in SMPs. Besides, the results of the toxic assessment indicated that the oxidative damage caused by the NaCl shock was an important factor affecting the property of SMPs, and the abnormal expression of DNA transcription cannot be neglected for bacteria metabolisms with the change of C/N ratio.

The molecular grammar of protein disorder guiding genome-binding locations

Intrinsically disordered regions (IDRs) direct transcription factors (TFs) towards selected genomic occurrences of their binding motif, as exemplified by budding yeast's Msn2. However, the sequence basis of IDR-directed TF binding selectivity remains unknown. To reveal this sequence grammar, we analyze the genomic localizations of >100 designed IDR mutants, each carrying up to 122 mutations within this 567-AA region. Our data points at multivalent interactions, carried by hydrophobic-mostly aliphatic-residues dispersed within a disordered environment and independent of linear sequence motifs, as the key determinants of Msn2 genomic localization. The implications of our results for the mechanistic basis of IDR-based TF binding preferences are discussed.

Complementary strategies for directing in vivo transcription factor binding through DNA binding domains and intrinsically disordered regions

DNA binding domains (DBDs) of transcription factors (TFs) recognize DNA sequence motifs that are highly abundant in genomes. Within cells, TFs bind a subset of motif-containing sites as directed by either their DBDs or DBD-external (nonDBD) sequences. To define the relative roles of DBDs and nonDBDs in directing binding preferences, we compared the genome-wide binding of 48 (∼30%) budding yeast TFs with their DBD-only, nonDBD-truncated, and nonDBD-only mutants. With a few exceptions, binding locations differed between DBDs and TFs, resulting from the cumulative action of multiple determinants mapped mostly to disordered nonDBD regions. Furthermore, TFs' preferences for promoters of the fuzzy nucleosome architecture were lost in DBD-only mutants, whose binding spread across promoters, implicating nonDBDs' preferences in this hallmark of budding yeast regulatory design. We conclude that DBDs and nonDBDs employ complementary DNA-targeting strategies, whose balance defines TF binding specificity along genomes.

Intra-Domain Residue Coevolution in Transcription Factors Contributes to DNA Binding Specificity

Understanding the basis of the DNA-binding specificity of transcription factors (TFs) has been of long-standing interest. Despite extensive efforts to map millions of putative TF binding sequences, identifying the critical determinants for DNA binding specificity remains a major challenge. The coevolution of residues in proteins occurs due to a shared evolutionary history. However, it is unclear how coevolving residues in TFs contribute to DNA binding specificity. Here, we systematically collected publicly available data sets from multiple large-scale high-throughput TF-DNA interaction screening experiments for the major TF families with large numbers of TF members. These families included the Homeobox, HLH, bZIP_1, Ets, HMG_box, ZF-C4, and Zn_clus TFs. We detected TF subclass-determining sites (TSDSs) and showed that the TSDSs were more likely to coevolve with other TSDSs than with non-TSDSs, particularly for the Homeobox, HLH, Ets, bZIP_1, and HMG_box TF families. By in silico modeling, we showed that mutation of the highly coevolving residues could significantly reduce the stability of the TF-DNA complex. The distant residues from the DNA interface also contributed to TF-DNA binding activity. Overall, our study gave evidence that coevolved residues relate to transcriptional regulation and provided insights into the potential application of engineered DNA-binding domains and proteins. IMPORTANCE While unraveling DNA-binding specificity of TFs is the key to understanding the basis and molecular mechanism of gene expression regulation, identifying the critical determinants that contribute to DNA binding specificity remains a major challenge. In this study, we provided evidence showing that coevolving residues in TF domains contributed to DNA binding specificity. We demonstrated that the TSDSs were more likely to coevolve with other TSDSs than with non-TSDSs. Mutation of the coevolving residue pairs (CRPs) could significantly reduce the stability of THE TF-DNA complex, and even the distant residues from the DNA interface contribute to TF-DNA binding activity. Collectively, our study expands our knowledge of the interactions among coevolved residues in TFs, tertiary contacting, and functional importance in refined transcriptional regulation. Understanding the impact of coevolving residues in TFs will help understand the details of transcription of gene regulation and advance the application of engineered DNA-binding domains and protein.

The Oncogenic FOXL2 C134W Mutation Is a Key Driver of Granulosa Cell Tumors

Adult-type granulosa cell tumors (AGCT) are the most common type of malignant ovarian sex cord-stromal tumors. Most AGCTs carry the somatic variant c.402C>G (p.C134W) affecting the transcription factor FOXL2. Germline dominant variants in FOXL2 are responsible for blepharophimosis syndrome, which is characterized by underdevelopment of the eyelid. In this work, we generated a mouse model harboring the C134W variant of FOXL2 to evaluate in vivo the poorly understood oncogenic role of FOXL2. The mutation was dominant regarding eyelid hypoplasia, reminiscent of blepharophimosis syndrome. Interestingly, Foxl2+/C134W female mice had reduced fertility and developed AGCTs through a progression from abnormal ovaries with aberrant granulosa cells to ovaries with stromal hyperplasia and atypia and on to tumors in adut mice. The genes dysregulated in mouse AGCTs exhibited the hallmarks of cancer and were consistent with a gain-of-function of the mutated allele affecting TGFβ signaling. A comparison of these data with previous results on human AGCTs indicated similar deregulated pathways. Finally, a mutational analysis of mouse AGCT transcriptomic data suggested the absence of additional driver mutations apart from FOXL2-C134W. These results provide a clear in vivo example in which a single mutational hit triggers tumor development associated with profound transcriptomic alterations.

SIGNIFICANCE

A newly generated mouse model carrying a FOXL2 mutation characteristic of adult-type granulosa cell tumors shows that FOXL2 C134W shifts the transcriptome towards a signature of granulosa cell cancer and drives tumorigenesis.

Sequence-Specific Structural Features and Solvation Properties of Transcription Factor Binding DNA Motifs: Insights from Molecular Dynamics Simulation

Sequence-specific recognition of transcription factor (TF) binding motifs in the target site of DNA over the vast amount of non-target DNA is of primary importance for the transcriptional regulation of gene expression by the TFs. Binding of TFs to the target site of DNA relies not only on the direct contact formation but also on the structural and conformational features of DNA. Recognition of DNA structural features or shape readout by proteins is an important factor in the context of TF-DNA interaction. Based on the atomistic molecular simulation, here we report the sequence-dependent unique structural features, solvation, and ion-binding properties of biologically relevant AT- and GC-rich human TF binding motifs in DNA. Counterion and water distribution around the motif is found to be sensitive to the motif sequence, which is accompanied with the DNA shape features. The motif sequence affects the electrostatic potential along the grooves, and cytosine methylation alters the DNA shape features. Characteristic solvation properties of TF binding motif DNA fragments infer that an ionic environment and hydration influences are essential to describe TF-DNA interactions.

M. tuberculosis Transcription Machinery: A Review on the Mycobacterial RNA Polymerase and Drug Discovery Efforts

Mycobacterium tuberculosis (MTB) is the main source of tuberculosis (TB), one of the oldest known diseases in the human population. Despite the drug discovery efforts of past decades, TB is still one of the leading causes of mortality and claimed more than 1.5 million lives worldwide in 2020. Due to the emergence of drug-resistant strains and patient non-compliance during treatments, there is a pressing need to find alternative therapeutic agents for TB. One of the important areas for developing new treatments is in the inhibition of the transcription step of gene expression; it is the first step to synthesize a copy of the genetic material in the form of mRNA. This further translates to functional protein synthesis, which is crucial for the bacteria living processes. MTB contains a bacterial DNA-dependent RNA polymerase (RNAP), which is the key enzyme for the transcription process. MTB RNAP has been targeted for designing and developing antitubercular agents because gene transcription is essential for the mycobacteria survival. Initiation, elongation, and termination are the three important sequential steps in the transcription process. Each step is complex and highly regulated, involving multiple transcription factors. This review is focused on the MTB transcription machinery, especially in the nature of MTB RNAP as the main enzyme that is regulated by transcription factors. The mechanism and conformational dynamics that occur during transcription are discussed and summarized. Finally, the current progress on MTB transcription inhibition and possible drug target in mycobacterial RNAP are also described to provide insight for future antitubercular drug design and development.

Molecular basis for lethal cross-talk between two unrelated bacterial transcription factors - the regulatory protein of a restriction-modification system and the repressor of a defective prophage

Bacterial gene expression depends on the efficient functioning of global transcriptional networks, however their interconnectivity and orchestration rely mainly on the action of individual DNA binding proteins called transcription factors (TFs). TFs interact not only with their specific target sites, but also with secondary (off-target) sites, and vary in their promiscuity. It is not clear yet what mechanisms govern the interactions with secondary sites, and how such rewiring affects the overall regulatory network, but this could clearly constrain horizontal gene transfer. Here, we show the molecular mechanism of one such off-target interaction between two unrelated TFs in Escherichia coli: the C regulatory protein of a Type II restriction-modification system, and the RacR repressor of a defective prophage. We reveal that the C protein interferes with RacR repressor expression, resulting in derepression of the toxic YdaT protein. These results also provide novel insights into regulation of the racR-ydaST operon. We mapped the C regulator interaction to a specific off-target site, and also visualized C protein dynamics, revealing intriguing differences in single molecule dynamics in different genetic contexts. Our results demonstrate an apparent example of horizontal gene transfer leading to adventitious TF cross-talk with negative effects on the recipient's viability. More broadly, this study represents an experimentally-accessible model of a regulatory constraint on horizontal gene transfer.

Transcription factor-nucleosome dynamics from plasma cfDNA identifies ER-driven states in breast cancer

Genome-wide binding profiles of estrogen receptor (ER) and FOXA1 reflect cancer state in ER+ breast cancer. However, routine profiling of tumor transcription factor (TF) binding is impractical in the clinic. Here, we show that plasma cell-free DNA (cfDNA) contains high-resolution ER and FOXA1 tumor binding profiles for breast cancer. Enrichment of TF footprints in plasma reflects the binding strength of the TF in originating tissue. We defined pure in vivo tumor TF signatures in plasma using ER+ breast cancer xenografts, which can distinguish xenografts with distinct ER states. Furthermore, state-specific ER-binding signatures can partition human breast tumors into groups with significantly different ER expression and mortality. Last, TF footprints in human plasma samples can identify the presence of ER+ breast cancer. Thus, plasma TF footprints enable minimally invasive mapping of the regulatory landscape of breast cancer in humans and open vast possibilities for clinical applications across multiple tumor types.

Dimerization processes for light-regulated transcription factor Photozipper visualized by high-speed atomic force microscopy

Dimerization is critical for transcription factors (TFs) to bind DNA and regulate a wide variety of cellular functions; however, the molecular mechanisms remain to be completely elucidated. Here, we used high-speed atomic force microscopy (HS-AFM) to observe the dimerization process for a photoresponsive TF Photozipper (PZ), which consists of light–oxygen–voltage-sensing (LOV) and basic-region-leucine-zipper (bZIP) domains. HS-AFM visualized not only the oligomeric states of PZ molecules forming monomers and dimers under controlled dark–light conditions but also the domain structures within each molecule. Successive AFM movies captured the dimerization process for an individual PZ molecule and the monomer–dimer reversible transition during dark–light cycling. Detailed AFM images of domain structures in PZ molecules demonstrated that the bZIP domain entangled under dark conditions was loosened owing to light illumination and fluctuated around the LOV domain. These observations revealed the role of the bZIP domain in the dimerization processes of a TF.

Chromatin dynamics controls epigenetic domain formation

In multicellular organisms, nucleosomes carry epigenetic information that defines distinct patterns of gene expression, which are inherited over multiple generations. The enhanced capacity for information storage arises by nucleosome modifications, which are triggered by enzymes. Modified nucleosomes can transfer the mark to others that are in proximity by a positive-feedback (modification begets modification) mechanism. We created a generic polymer model, referred to as 3DSpreader, in which each bead, representing a nucleosome, stochastically switches between unmodified (U) and modified (M) states depending on the states of the neighbors. Modification begins at a specific nucleation site (NS) that is permanently in the M state, and could spread to other loci that is dictated by chromatin dynamics. Transfer of marks among the non-nucleation loci occurs stochastically as chromatin evolves in time. If the spreading rate is slower than the chromatin relaxation rate, which is biologically pertinent, then finite-sized domains form, driven by contacts between nucleosomes through a three-dimensional looping mechanism. Surprisingly, simulations based on the 3DSpreader model result in finite bounded domains that arise without the need for any boundary elements. Maintenance of spatially and temporally stable domains requires the presence of the NS, whose removal eliminates finite-sized modified domains. The theoretical predictions are in excellent agreement with experimental data for H3K9me3 spreading in mouse embryonic stem cells.

Genome-Wide Identification of MDH Family Genes and Their Association with Salt Tolerance in Rice

Malate dehydrogenase (MDH) is widely present in nature and regulates plant growth and development, as well as playing essential roles, especially in abiotic stress responses. Nevertheless, there is no comprehensive knowledge to date on MDH family members in rice. In this study, a total of 12 MDH members in rice were identified through genome-wide analysis and divided into three groups on the basis of their phylogenetic relationship and protein-conserved motifs. Evolutionary analysis showed that MDH proteins from rice, maize and wheat shared a close phylogenetic relationship, and the MDH family was conserved in the long-term process of domestication. We identified two segmental duplication events involving four genes, which could be the major force driving the expansion of the OsMDH family. The expression profile, cis-regulatory elements and qRT-PCR results of these genes revealed that a few OsMDH showed high tissue specificity, almost all of which had stress response elements in the promoter region, and ten MDH members were significantly induced by salt stress. Through gene-based association analysis, we found a significant correlation between salt tolerance at the seedling stage and the genetic variation of OsMDH8.1 and OsMDH12.1. Additionally, we found that the polymorphism in the promoter region of OsMDH8.1 might be related to the salt tolerance of rice. This study aimed to provide valuable information on the functional study of the rice MDH gene family related to salt stress response and revealed that OsMDH8.1 might be an important gene for the cultivar improvement of salt tolerance in rice.

LcNAC13 Is Involved in the Reactive Oxygen Species-Dependent Senescence of the Rudimentary Leaves in Litchi chinensis

Litchi is an important evergreen fruit tree. Floral formation in litchi is induced by low temperatures (LTs). However, unstable flowering is a challenge for litchi production in times of global warming and climate change. Previous studies have shown that the methyl viologen dichloride hydrate-generated reactive oxygen species (ROS) could promote flowering. Leaves in the panicles may affect the development of the inflorescence in litchi under high-temperature condition. In this study, potted litchi trees were transferred to growth chambers at LT and high temperature (HT). From a previous dataset of the RNA sequencing of the ROS-treated rudimentary leaves, a NAC transcription factor-encoding gene LcNAC13 was identified. By genetic transformation of LcNAC13 to Arabidopsis thaliana and tobacco, it was found that the ROS-induced senescence of the leaves was accelerated. Silencing LcNAC13 by virus-induced gene silencing (VIGS) delayed ROS-dependent senescence. Our results suggested that LcNAC13 regulates rudimentary leaf senescence. Our study provided a new target gene for the future molecular breeding of new cultivars that could flower under global warming conditions.

Genome-Wide Identification of DNA Binding with One Finger (Dof) Gene Family in Tartary Buckwheat (Fagopyrum tataricum) and Analysis of Its Expression Pattern after Exogenous Hormone Stimulation

Simple Summary

A number of studies have demonstrated that DNA binding with one finger (Dof) proteins are involved in multiple biological processes. In the present study, Dof genes or proteins in Tartary buckwheat (FtDofs) were systematically analysed, including their physical properties, phylogenetic relationships, structure, motif composition, cis-acting elements present in promoter regions, chromosomal distribution, gene duplication events, syntenic relationships, expression patterns in different tissues and different fruit developmental stages and responses to exogenous hormone stimulation. The results indicated that the expansion of FtDofs was mainly due to segmental duplication. The tissue-specific expression patterns of FtDofs and their positive responses to exogenous hormone stimulation suggest that they play important roles in the growth and development of Tartary buckwheat as well as in the adaptation to environmental changes. Collectively, this study lays a foundation for further exploration of the function of FtDof genes in Tartary buckwheat.

Abstract

DNA binding with one finger (Dof) proteins have been proven to be involved in multiple biological processes. However, genome-wide identification of the Dof gene family has not been reported for Tartary buckwheat (Fagopyrum tataricum). In this study, 35 FtDof proteins were identified, and they could be divided into nine phylogenetic subgroups. Proteins within the same subgroup had similar gene structure and motif composition. Moreover, abundant cis-acting elements were present in the promoter regions of FtDof genes. Segmental duplication was the primary driving force for the evolution of the FtDof gene family. Synteny analysis indicated that Tartary buckwheat was closer to dicotyledons, and more orthologous Dof genes existed among them. The expression pattern of FtDofs in different tissues and at different fruit developmental stages varied. Different tissues contained several genes that were specifically expressed. FtDof expression was mainly upregulated under methyl jasmonate treatment and downregulated under other hormone treatments. Taken together, FtDofs may play important roles in the growth and development of Tartary buckwheat and in response to abiotic and biotic stresses. Therefore, the genome-wide identification and expression pattern analysis of the Tartary buckwheat Dof gene family lays a foundation for further exploration of the functional characteristics of FtDofs in the future.

Trends in biological data integration for the selection of enzymes and transcription factors related to cellulose and hemicellulose degradation in fungi

Fungi are key players in biotechnological applications. Although several studies focusing on fungal diversity and genetics have been performed, many details of fungal biology remain unknown, including how cellulolytic enzymes are modulated within these organisms to allow changes in main plant cell wall compounds, cellulose and hemicellulose, and subsequent biomass conversion. With the advent and consolidation of DNA/RNA sequencing technology, different types of information can be generated at the genomic, structural and functional levels, including the gene expression profiles and regulatory mechanisms of these organisms, during degradation-induced conditions. This increase in data generation made rapid computational development necessary to deal with the large amounts of data generated. In this context, the origination of bioinformatics, a hybrid science integrating biological data with various techniques for information storage, distribution and analysis, was a fundamental step toward the current state-of-the-art in the postgenomic era. The possibility of integrating biological big data has facilitated exciting discoveries, including identifying novel mechanisms and more efficient enzymes, increasing yields, reducing costs and expanding opportunities in the bioprocess field. In this review, we summarize the current status and trends of the integration of different types of biological data through bioinformatics approaches for biological data analysis and enzyme selection.

Proteome-Wide Profiling of Readers for DNA Modification

DNA modifications, represented by 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC), play important roles in epigenetic regulation of biological processes. The specific recognition of DNA modifications by the transcriptional protein machinery is thought to be a potential mechanism for epigenetic-driven gene regulation, and many modified DNA-specific binding proteins have been uncovered. However, the panoramic view of the roles of DNA modification readers at the proteome level remains largely unclear. Here, a recently developed concatenated tandem array of consensus transcription factor (TF) response elements (catTFREs) approach is employed to profile the binding activity of TFs at DNA modifications. Modified DNA-binding activity is quantified for 1039 TFs, representing 70% of the TFs in the human genome. Additionally, the modified DNA-binding activity of 600 TFs is monitored during the mouse brain development from the embryo to the adult stages. Readers of these DNA modifications are predicted, and the hierarchical networks between the transcriptional protein machinery and modified DNA are described. It is further demonstrated that ZNF24 and ZSCAN21 are potential readers of 5fC-modified DNA. This study provides a landscape of TF-DNA modification interactions that can be used to elucidate the epigenetic-related transcriptional regulation mechanisms under physiological conditions.

Genome-Wide Analysis and Expression Profiling of the Phospholipase D Gene Family in Solanum tuberosum

Phospholipase D (PLD) is the most important phospholipid hydrolase in plants, which can hydrolyze phospholipids into phosphatidic acid (PA) and choline. When plants encounter low temperature, drought and high salt stress, phospholipase D and its products play an important role in regulating plant growth and development and coping with stress. In this study, 16 members of StPLD gene family were identified in potato genome, which were distributed in α, β, δ, and ζ subfamilies, and their expression patterns under salt, high temperature, drought, and ABA stress were detected by qRT-PCR method. Gene expression analysis showed that the expression of StPLD genes in potato was upregulated and downregulated to varying degrees under the four stresses, indicating that the PLD gene family is involved in the interaction of potato plant hormones and abiotic stress signals. Chromosome distribution showed that StPLD gene was unevenly distributed on 8 chromosomes, and only one pair of tandem repeat genes was found. All StPLD promoters contain hormone and stress-related cis-regulatory elements to respond to different stresses. Structural analysis showed that StPLD genes in the same subgroup had a similar exon-intron structure. Our study provides a valuable reference for further research of the function and structure of PLD gene.

Identification of consistent post-translational regulatory triplets related to oncogenic and tumour suppressive modulators in childhood acute lymphoblastic leukemia

Background

Acute lymphoblastic leukemia (ALL) is the most common type of childhood cancer. It can be caused by mutations that turn on oncogenes or turn off tumour suppressor genes. For instance, changes in certain genes including Rb and p53 are common in ALL cells. Oncogenes and TSGs may serve as a modulator gene to regulate the gene expression level via their respective target genes. To investigate the regulatory relationship between oncogenes, tumour suppressor genes and transcription factors at the post translational level in childhood ALL, we performed an integrative network analysis on the gene regulation in the post-translational level for childhood ALL based on many publicly available cancer gene expression data including TARGET and GEO database.

Methods

We collected 259 childhood ALL-related genes from the latest online leukemia database, Leukemia Gene Literature Database. These 259 genes were selected from a comprehensive systematic literature with experimental evidences. The identified and curated genes were also associated with patient survival cases and we incorporated this pediatric ALL-related gene list into our analysis. We extracted the known human TFs from the TRRUST database. Among 259 childhood ALL-related genes, 101 unique regulators were mapped to the list of oncogene and tumour suppressor genes (TSGs) from the ONGene and the TSGene databases, and these included 74 TSGs, 62 oncogenes and 46 TF genes.

Results

The resulted regulation was presented as a hierarchical regulatory network with transcription factors (TFs) as intermediate regulators connecting the top modulators (oncogene and TSGs) to the common target genes. Cross-validation was applied to the results from the TARGET dataset by identifying the consistent regulatory motifs based on three independent ALL expression datasets. A three-layer regulatory network of consistent positive modulators in childhood ALL was constructed in which 74 modulators (40 oncogenes, 34 TSGs) are considered as the most important regulators. The middle layer and the bottom layer contain 34 TFs and 176 target genes, respectively. Oncogenes mostly participated in positive regulation of gene expression and the transcription process of RNA II polymerase, while TSGs were mainly involved in the negative regulation of gene expression. In addition, the oncogene-specific targets were enriched with regulators of the MAPK cascade while tumour suppressor-specific targets were associated with cell death.

Conclusion

The results revealed that oncogenes and TSGs possess a different functional regulatory pattern with regard to not only their biological functions but also their specific target genes in childhood ALL cancer progression. Taken together, our findings could contribute to a better understanding of the important regulatory mechanisms and this method could be used to analyse the targeted genes at the post-translational level in childhood ALL through integrative network analysis.

Time-Based Systems Biology Approaches to Capture and Model Dynamic Gene Regulatory Networks

All aspects of transcription and its regulation involve dynamic events. However, capturing these dynamic events in gene regulatory networks (GRNs) offers both a promise and a challenge. The promise is that capturing and modeling the dynamic changes in GRNs will allow us to understand how organisms adapt to a changing environment. The ability to mount a rapid transcriptional response to environmental changes is especially important in nonmotile organisms such as plants. The challenge is to capture these dynamic, genome-wide events and model them in GRNs. In this review, we cover recent progress in capturing dynamic interactions of transcription factors with their targets-at both the local and genome-wide levels-and how they are used to learn how GRNs operate as a function of time. We also discuss recent advances that employ time-based machine learning approaches to forecast gene expression at future time points, a key goal of systems biology.

NF-YB-Mediated Active Responses of Plant Growth under Salt and Temperature Stress in Eucalyptus grandis

The transcription factor NF-YB (nuclear factor-YB) family is a subfamily of the nuclear factor Y (NF-Y), which plays an important role in regulating plant growth, development and participates in various stress responses. Although the NF-Y family has been studied in many species, it is still obscure in Eucalyptus grandis. In this study, 23 EgNF-YB genes in eucalyptus were identified and unevenly distributed on 11 chromosomes. Phylogenetic analysis showed the EgNF-YB genes were divided into two clades, LEC-1 type and non-LEC1 type. The evolution of distinct clades was relatively conservative, the gene structures were analogous, and the differences of genetic structures among clades were small. The expression profiles showed that the distinct EgNF-YB genes were highly expressed in diverse tissues, and EgNF-YB4/6/13/19/23 functioned in response to salinity, heat and cold stresses. Our study characterized the phylogenetic relationship, gene structures and expression patterns of EgNF-YB gene family and investigated their potential roles in abiotic stress responses, which provides solid foundations for further functional analysis of NF-YB genes in eucalyptus.

Role of a local transcription factor in governing cellular carbon/nitrogen homeostasis in Pseudomonas fluorescens

Autoactivation of two-component systems (TCSs) can increase the sensitivity to signals but inherently cause a delayed response. Here, we describe a unique negative feedback mechanism enabling the global NtrB/NtrC regulator to rapidly respond to nitrogen starvation over the course of histidine utilization (hut) in Pseudomonas fluorescens. NtrBC directly activates transcription of hut genes, but overexpression will produce excess ammonium leading to NtrBC inactivation. To prevent this from occurring, the histidine-responsive repressor HutC fine-tunes ntrBC autoactivation: HutC and NtrC bind to the same operator site in the ntrBC promoter. This newly discovered low-affinity binding site shows little sequence similarity with the consensus sequence that HutC recognizes for substrate-specific induction of hut operons. A combination of genetic and transcriptomic analysis indicated that both ntrBC and hut promoter activities cannot be stably maintained in the ΔhutC background when histidine fluctuates at high concentrations. Moreover, the global carbon regulator CbrA/CbrB is involved in directly activating hut transcription while de-repressing hut translation via the CbrAB-CrcYZ-Crc/Hfq regulatory cascade. Together, our data reveal that the local transcription factor HutC plays a crucial role in governing NtrBC to maintain carbon/nitrogen homeostasis through the complex interactions between two TCSs (NtrBC and CbrAB) at the hut promoter.

Genome-wide identification and expression analysis of ADP-ribosylation factors associated with biotic and abiotic stress in wheat (Triticum aestivum L.)

The ARF gene family plays important roles in intracellular transport in eukaryotes and is involved in conferring tolerance to biotic and abiotic stresses in plants. To explore the role of these genes in the development of wheat (Triticum aestivum L.), 74 wheat ARF genes (TaARFs; including 18 alternate transcripts) were identified and clustered into seven sub-groups. Phylogenetic analysis revealed that TaARFA1 sub-group genes were strongly conserved. Numerous cis-elements functionally associated with the stress response and hormones were identified in the TaARFA1 sub-group, implying that these TaARFs are induced in response to abiotic and biotic stresses in wheat. According to available transcriptome data and qRT-PCR analysis, the TaARFA1 genes displayed tissue-specific expression patterns and were regulated by biotic stress (powdery mildew and stripe rust) and abiotic stress (cold, heat, ABA, drought and NaCl). Protein interaction network analysis further indicated that TaARFA1 proteins may interact with protein phosphatase 2C (PP2C), which is a key protein in the ABA signaling pathway. This comprehensive analysis will be useful for further functional characterization of TaARF genes and the development of high-quality wheat varieties.

Optimized combination methods for exploring and verifying disease-resistant transcription factors in melon

A large amount of omics data and number of bioinformatics tools has been produced. However, the methods for further exploring omics data are simple, in particular, to mine key regulatory genes, which are a priority concern in biological systems, and most of the specific functions are still unknown. First, raw data of two genotypes of melon (susceptible and resistant) were obtained by transcriptome analysis. Second, 391 transcription factors (TFs) were identified from the plant transcription factor database and cucurbit genomics database. Then, functional enrichment analysis indicated that these genes were mainly annotated in the process of transcription regulation. Third, 243 and 230 module-specific TFs were screened by weighted gene coexpression network analysis and short time series expression miner, respectively. Several TF genes, such as WRKYs and bHLHs, were regarded as key regulatory genes according to the values of significantly different modules. The coexpression network showed that these TF genes were significant correlated with resistance (R) genes, such as DRP2, RGA3, DRP1 and NB-ARC. Fourth, cis-acting element analysis illustrated that these R genes may bind to WRKY and bHLH. Finally, the expression of WRKY genes was verified by quantitative reverse transcription PCR (RT-qPCR). Phylogenetic analysis was carried out to further confirm that these TFs may play a critical role in Curcurbitaceae disease resistance. This study provides a new optimized combination strategy to explore the functions of TFs in a wide spectrum of biological processes. This strategy may also effectively predict potential relationships in the interactions of essential genes.

Design of the Crosslinking Reactions for Nucleic Acids-Binding Protein and Evaluation of the Reactivity

Selective chemical reactions of biomolecules are some of the important tools for investigations by biological studies. We have developed the selective crosslinking reactions to form covalent bonds to DNA or RNA using crosslinking oligonucleotides (CFO) bearing reactive bases. In this study, we designed the cross-linkable 4-amino-6-oxo-2-vinyltriazine derivative with an acyclic linker (acyAOVT) to react with the nucleic acids-binding protein based on our previous results. We hypothesized that the acyAOVT base would form a stable base pair with guanine by three hydrogen bonds at the positions of the vinyl group in the duplex DNA major groove, and the vinyl group can react with the nucleophilic species in the proximity, for example, the cysteine or lysine residue in the nucleic acids-binding protein. The synthesized oligonucleotides bearing the acyAOVT derivative showed a higher reactivity than that of the corresponding pyrimidine derivative without one nitrogen. The duplex containing acyAOVT-guanine (G) formed complexes with Hha1 DNMT even in the presence of 2-mercaptoethanol. We expect that our system will provide a useful tool for the molecular study of nucleic acids-binding proteins.

Mapping Regulatory Determinants in Plants

The domestication and improvement of many plant species have frequently involved modulation of transcriptional outputs and continue to offer much promise for targeted trait engineering. The cis-regulatory elements (CREs) controlling these trait-associated transcriptional variants however reside within non-coding regions that are currently poorly annotated in most plant species. This is particularly true in large crop genomes where regulatory regions constitute only a small fraction of the total genomic space. Furthermore, relatively little is known about how CREs function to modulate transcription in plants. Therefore understanding where regulatory regions are located within a genome, what genes they control, and how they are structured are important factors that could be used to guide both traditional and synthetic plant breeding efforts. Here, we describe classic examples of regulatory instances as well as recent advances in plant regulatory genomics. We highlight valuable molecular tools that are enabling large-scale identification of CREs and offering unprecedented insight into how genes are regulated in diverse plant species. We focus on chromatin environment, transcription factor (TF) binding, the role of transposable elements, and the association between regulatory regions and target genes.

Mir214-3p and Hnf4a/Hnf4α reciprocally regulate Ulk1 expression and autophagy in nonalcoholic hepatic steatosis

Macroautophagy/autophagy, a self-degradative process, regulates metabolic homeostasis in response to various stress conditions and is a therapeutic target for nonalcoholic fatty liver disease. We found that autophagic activity was inhibited as a result of a significant reduction in the expression of autophagy-related genes such as Ulk1 in a mouse model and patients with fatty liver. This downregulation was caused by increased Mir214-3p levels and decreased Hnf4a/Hnf4α mRNA levels in hepatocytes. Mir214-3p suppressed Ulk1 expression through direct binding at a 3' untranslated region sequence. Hnf4a directly activated transcription of Ulk1. We investigated lipid accumulation and the expression of autophagy-related genes in the livers of mice treated with anti-Mir214-3p. Hepatic steatosis was alleviated, and Ulk1 mRNA levels were significantly increased by locked nucleic acid-mediated Mir214-3p silencing. Additionally, autophagosome formation and MAP1LC3/LC3-II protein levels were increased, indicating an increase in autophagic activity. Interestingly, suppression of Mir214-3p did not ameliorate fatty liver under Ulk1 suppression, suggesting that reduced Mir214-3p levels mitigate hepatic steatosis through upregulation of Ulk1. These results demonstrate that inhibition of Mir214-3p expression ameliorated fatty liver disease through increased autophagic activity by increasing the expression of Ulk1. Thus, Mir214-3p is a potential therapeutic target for nonalcoholic fatty disease.Abbreviations: AMPK: adenosine monophosphate-activated protein kinase; ATG: autophagy-related; ChIP: chromatin immunoprecipitation; CTSB: cathepsin B; CTSL: cathepsin L; CQ: chloroquine; HFD: high-fat diet; HNF4A: hepatocyte nuclear factor 4, alpha; IF: immunofluorescence; IHC: immunohistochemistry; LDs: lipid droplets; Leup: leupeptin; LFD: low-fat diet; LNA: locked nucleic acid; MAP1LC3B/LC3B: microtubule-associated protein 1 light chain 3 beta; miRNA: microRNA; MTOR: mechanistic target of rapamycin kinase; NAFLD: non-alcoholic fatty liver disease; NASH: non-alcoholic steatohepatitis; PCR: polymerase chain reaction; TEM: transmission electron microscopy; TF: transcription factor; TLDA: TaqMan low-density array; ULK1: unc-51 like kinase 1; UTR: untranslated region.

Transcriptomic Analysis and Specific Expression of Transcription Factor Genes in the Root and Sporophyll of Dryopteris fragrans (L.) Schott

Background: Dryopteris fragrans, which is densely covered with glandular trichomes, is considered to be one of the ferns with the most medicinal potential. The transcriptomes from selected tissues of D. fragrans were collected and analyzed for functional and comparative genomic studies. The aim of this study was to determine the transcriptomic characteristics of wild D. fragrans sporangium in tissues from the SR (root), SL (sporophyll), and TRL (sporophyll with glandular trichomes removed). Results: Cluster analysis identified genes that were highly expressed in an organ-specific manner according to read mapping, feature counting, and normalization. The functional map identified gene clusters that can uniquely describe the function of each tissue. We identified a group of three tissue-specific transcription factors targeting the SL, SR, and TRL. In addition, highly expressed transcription factors (TFs) were found in each tissue-specific gene cluster, where ERF and bHLH transcription factors were the two types showing the most distinct expression patterns between the three different tissues. The specific expression of transcription factor genes varied between the different types of tissues. The numbers of transcription factors specifically expressed in the roots and sporophylls were 60 and 30, respectively, while only seven were found for the sporophylls with glandular trichomes removed. The expression of genes known to be associated with the development of glandular trichomes in flowering plants, including MIXTA, ATML1, and MYB106, were also validated and are discussed. In particular, a unigene encoding MIXTA was identified and exhibited the highest expression level in SL in D. fragrans. Conclusions: This study is the first report of global transcriptomic analysis in different tissues of D. fragrans, and the first to discuss these findings in the context of the development of homologous glandular trichomes. These results set the stage for further research on the development, stress resistance, and secondary metabolism of D. fragrans glandular trichomes.

Genome-wide identification and expression profiles of ERF subfamily transcription factors in Zea mays

The Ethylene-Response Factor (ERF) subfamily transcription factors (TFs) belong to the APETALA2/Ethylene-Responsive Factor (AP2/ERF) superfamily and play a vital role in plant growth and development. However, identification and analysis of the ERF subfamily genes in maize have not yet been performed at genome-wide level. In this study, a total of 76 ERF subfamily TFs were identified and were found to be unevenly distributed on the maize chromosomes. These maize ERF (ZmERF) TFs were classified into six groups, namely groups B1 to B6, based on phylogenetic analysis. Synteny analysis showed that 50, 54, and 58 of the ZmERF genes were orthologous to those in rice, Brachypodium, and Sorghum, respectively. Cis-element analysis showed that elements related to plant growth and development, hormones, and abiotic stress were identified in the promoter region of ZmERF genes. Expression profiles suggested that ZmERF genes might participate in plant development and in response to salinity and drought stresses. Our findings lay a foundation and provide clues for understanding the biological functions of ERF TFs in maize.

IGAP-integrative genome analysis pipeline reveals new gene regulatory model associated with nonspecific TF-DNA binding affinity

The human genome is regulated in a multi-dimensional way. While biophysical factors like Non-specific Transcription factor Binding Affinity (nTBA) act at DNA sequence level, other factors act above sequence levels such as histone modifications and 3-D chromosomal interactions. This multidimensionality of regulation requires many of these factors for a proper understanding of the regulatory landscape of the human genome. Here, we propose a new biophysical model for estimating nTBA. Integration of nTBA with chromatin modifications and chromosomal interactions, using a new Integrative Genome Analysis Pipeline (IGAP), reveals additive effects of nTBA to regulatory DNA sequences and identifies three types of genomic zones in the human genome (Inactive Genomic Zones, Poised Genomic Zones, and Active Genomic Zones). It also unveils a novel long distance gene regulatory model: chromosomal interactions reduce the physical distance between the high occupancy target (HOT) regions that results in high nTBA to DNA in the area, which in turn attract TFs to such regions with higher binding potential. These findings will help to elucidate the three-dimensional diffusion process that TFs use during their search for the right targets.

A report on DNA sequence determinants in gene expression

The biased usage of nucleotides in coding sequence and its correlation with gene expression has been observed in several studies. A complex set of interactions between genes and other components of the expression system determine the amount of proteins produced from coding sequences. It is known that the elongation rate of polypeptide chain is affected by both codon usage bias and specific amino acid compositional constraints. Therefore, it is of interest to review local DNA-sequence elements and other positional as well as combinatorial constraints that play significant role in gene expression.

A Conformational Switch in the Zinc Finger Protein Kaiso Mediates Differential Readout of Specific and Methylated DNA Sequences

Recognition of the epigenetic mark 5-methylcytosine (mC) at CpG sites in DNA has emerged as a novel function of many eukaryotic transcription factors (TFs). It remains unclear why the sequence specificity of these TFs differs for CpG-methylated motifs and consensus motifs. Here, we dissect the structural and dynamic basis for this differential DNA binding specificity in the human zinc finger TF Kaiso, which exhibits high affinity for two consecutive mCpG sites in variable contexts and also for a longer, sequence-specific Kaiso binding site (KBS). By integrating structural analysis and DNA binding studies with targeted protein mutagenesis and nucleotide substitutions, we identify distinct mechanisms for readout of methylated and KBS motifs by Kaiso. We show that a key glutamate residue (E535), critical for mCpG site recognition, adopts different conformations in complexes with specific and methylated DNA. These conformational differences, together with intrinsic variations in DNA flexibility and/or solvation at TpG versus mCpG sites, contribute to the different DNA affinity and sequence specificity. With methylated DNA, multiple direct contacts between E535 and the 5' mCpG site dominate the binding affinity, allowing for tolerance of different flanking DNA sequences. With KBS, Kaiso employs E535 as part of an indirect screen of the 5' flanking sequence, relying on key tyrosine-DNA interactions to stabilize an optimal DNA conformation and select against noncognate sites. These findings demonstrate how TFs use conformational adaptation and exploit variations in DNA flexibility to achieve distinct DNA readout outcomes and target a greater variety of regulatory and epigenetic sites than previously appreciated.

Retrons and their applications in genome engineering

Precision genome editing technologies have transformed modern biology. These technologies have arisen from the redirection of natural biological machinery, such as bacteriophage lambda proteins for recombineering and CRISPR nucleases for eliciting site-specific double-strand breaks. Less well-known is a widely distributed class of bacterial retroelements, retrons, that employ specialized reverse transcriptases to produce noncoding intracellular DNAs. Retrons' natural function and mechanism of genetic transmission have remained enigmatic. However, recent studies have harnessed their ability to produce DNA in situ for genome editing and evolution. This review describes retron biology and function in both natural and synthetic contexts. We also highlight areas that require further study to advance retron-based precision genome editing platforms.

Genomics, molecular and evolutionary perspective of NAC transcription factors

NAC (NAM, ATAF1,2, and CUC2) transcription factors are one of the largest transcription factor families found in the plants and are involved in diverse developmental and signalling events. Despite the availability of comprehensive genomic information from diverse plant species, the basic genomic, biochemical, and evolutionary details of NAC TFs have not been established. Therefore, NAC TFs family proteins from 160 plant species were analyzed in the current study. Study revealed, Brassica napus (410) encodes highest number and Klebsormidium flaccidum (3) encodes the lowest number of TFs. The study further revealed the presence of NAC TF in the Charophyte algae K. flaccidum. On average, the monocot plants encode higher number (141.20) of NAC TFs compared to the eudicots (125.04), gymnosperm (75), and bryophytes (22.66). Furthermore, our analysis revealed that several NAC TFs are membrane bound and contain monopartite, bipartite, and multipartite nuclear localization signals. NAC TFs were also found to encode several novel chimeric proteins and regulate a complex interactome network. In addition to the presence of NAC domain, several NAC proteins were found to encode other functional signature motifs as well. Relative expression analysis of NAC TFs in A. thaliana revealed root tissue treated with urea and ammonia showed higher level of expression and leaf tissues treated with urea showed lower level of expression. The synonymous codon usage is absent in the NAC TFs and it appears that they have evolved from orthologous ancestors and undergone vivid duplications to give rise to paralogous NAC TFs. The presence of novel chimeric NAC TFs are of particular interest and the presence of chimeric NAC domain with other functional signature motifs in the NAC TF might encode novel functional properties in the plants.

Identification of pancreatic cancer type related factors by Weighted Gene Co-Expression Network Analysis

This study aims to identify the core modules associated with pancreatic cancer (PC) types and the ncRNAs and transcription factors (TFs) that regulate core module genes by weighted gene co-expression network analysis (WGCNA). WGCNA was used to analyze the union of genes related to PC in NCBI and OMIM databases and the differentially expressed genes screened by TCGA-PAAD database. Samples were clustered according to gene expression in gene modules and Fisher exact method was performed. GO and KEGG were used for enrichment analysis to visually display module genes and screen driver genes. Hypergeometric test method was used to calculate pivot nodes among ncRNAs, TFs and mRNA based on RAID 2.0 and TRRUST v2 databases. The blue and yellow modules were identified as the core modules associated with PC types. MST1R, TMPRSS, MIR198, SULF1, COL1A1 and FAP were the core genes in the modules. Hypergeometric test results showed that ANCR, miR-3134, MT1DP, LOC154449, LOC28329 and other ncRNAs were key factors driving blue module genes, while LINC-ROR, UCA1, SNORD114-4, HEIH, SNORD114-6 and other ncRNAs were key factors driving yellow module genes. TFs with significant regulatory effect on blue module included LCOR, PIAS4, ZEB1, SNAI2, SMARCA4, etc. and on yellow module included HOXC6, PER2, HOXD3, TWIST2, VHL, etc. The core modules associated with PC types were proved as yellow and blue modules, and important ncRNAs and TFs regulating yellow and blue modules were found. This study provides relevant evidence for further identification of PC types.

Parsimonious evolutionary scenario for the origin of allostery and coevolution patterns in proteins

Proteins display generic properties that are challenging to explain by direct selection, notably allostery, the capacity to be regulated through long-range effects, and evolvability, the capacity to adapt to new selective pressures. An evolutionary scenario is proposed where proteins acquire these two features indirectly as a by-product of their selection for a more fundamental property, exquisite discrimination, the capacity to bind discriminatively very similar ligands. Achieving this task is shown to typically require proteins to undergo a conformational change. We argue that physical and evolutionary constraints impel this change to be controlled by a group of sites extending from the binding site. Proteins can thus acquire a latent potential for allosteric regulation and evolutionary adaptation because of long-range effects that initially arise as evolutionary spandrels. This scenario accounts for the groups of conserved and coevolving residues observed in multiple sequence alignments. However, we propose that most pairs of coevolving and contacting residues inferred from such alignments have a different origin, related to thermal stability. A physical model is presented that illustrates this evolutionary scenario and its implications. The scenario can be implemented in experiments of protein evolution to directly test its predictions.

An improved bind-n-seq strategy to determine protein-DNA interactions validated using the bacterial transcriptional regulator YipR

Background

Interactions between transcription factors and DNA lie at the centre of many biological processes including DNA recombination, replication, repair and transcription. Most bacteria encode diverse proteins that act as transcription factors to regulate various traits. Several technologies for identifying protein–DNA interactions at the genomic level have been developed. Bind-n-seq is a high-throughput in vitro method first deployed to analyse DNA interactions associated with eukaryotic zinc-finger proteins. The method has three steps (i) binding protein to a randomised oligonucleotide DNA target library, (ii) deep sequencing of bound oligonucleotides, and (iii) a computational algorithm to define motifs among the sequences. The classical Bind-n-seq strategy suffers from several limitations including a lengthy wet laboratory protocol and a computational algorithm that is difficult to use. We introduce here an improved, rapid, and simplified Bind-n-seq protocol coupled with a user-friendly downstream data analysis and handling algorithm, which has been optimized for bacterial target proteins. We validate this new protocol by showing the successful characterisation of the DNA-binding specificities of YipR (YajQ interacting protein regulator), a well-known transcriptional regulator of virulence genes in the bacterial phytopathogen Xanthomonas campestris pv. campestris (Xcc).

Results

The improved Bind-n-seq approach identified several DNA binding motif sequences for YipR, in particular the CCCTCTC motif, which were located in the promoter regions of 1320 Xcc genes. Informatics analysis revealed that many of these genes regulate functions associated with virulence, motility, and biofilm formation and included genes previously found involved in virulence. Additionally, electromobility shift assays show that YipR binds to the promoter region of XC_2633 in a CCCTCTC motif-dependent manner.

Conclusion

We present a new and rapid Bind-n-seq protocol that should be useful to investigate DNA-binding proteins in bacteria. The analysis of YipR DNA binding using this protocol identifies a novel DNA sequence motif in the promoter regions of target genes that define the YipR regulon.

Characterization of Human Dosage-Sensitive Transcription Factor Genes

Copy number changes in protein-coding genes are detrimental if the consequent changes in protein concentrations disrupt essential cellular functions. The dosage sensitivity of transcription factor (TF) genes is particularly interesting because their products are essential in regulating the expression of genetic information. From four recently curated data sets of dosage-sensitive genes (genes with conserved copy numbers across mammals, ohnologs, and two data sets of haploinsufficient genes), we compiled a data set of the most reliable dosage-sensitive (MRDS) genes and a data set of the most reliable dosage-insensitive (MRDIS) genes. The MRDS genes were those present in all four data sets, while the MRDIS genes were those absent from any one of the four data sets and with the probability of being loss of function-intolerant (pLI) values < 0.5 in both of the haploinsufficient gene data sets. Enrichment analysis of TF genes among the MRDS and MRDIS gene data sets showed that TF genes are more likely to be dosage-sensitive than other genes in the human genome. The nuclear receptor family was the most enriched TF family among the dosage-sensitive genes. TF families with very few members were also deemed more likely to be dosage-sensitive than TF families with more members. In addition, we found a certain number of dosage-insensitive TFs. The most typical were the Krüppel-associated box domain-containing zinc-finger proteins (KZFPs). Gene ontology (GO) enrichment analysis showed that the MRDS TFs were enriched for many more terms than the MRDIS TFs; however, the proteins interacting with these two groups of TFs did not show such sharp differences. Furthermore, we found that the MRDIS KZFPs were not significantly enriched for any GO terms, whereas their interacting proteins were significantly enriched for thousands of GO terms. Further characterizations revealed significant differences between MRDS TFs and MRDIS TFs in the lengths and nucleotide compositions of DNA-binding sites as well as in expression level, protein size, and selective force.

Transcriptome analyses of cells carrying the Type II Csp231I restriction-modification system reveal cross-talk between two unrelated transcription factors: C protein and the Rac prophage repressor

Restriction-modification (R–M) systems represent an effective mechanism of defence against invading bacteriophages, and are widely spread among bacteria and archaea. In acquiring a Type II R–M system via horizontal gene transfer, the new hosts become more resistant to phage infection, through the action of a restriction endonuclease (REase), which recognizes and cleaves specific target DNAs. To protect the host cell's DNA, there is also a methyltransferase (MTase), which prevents DNA cleavage by the cognate REase. In some R–M systems, the host also accepts a cis-acting transcription factor (C protein), which regulates the counteracting activities of REase and MTase to avoid host self-restriction. Our study characterized the unexpected phenotype of Escherichia coli cells, which manifested as extensive cell filamentation triggered by acquiring the Csp231I R–M system from Citrobacter sp. Surprisingly, we found that the cell morphology defect was solely dependent on the C regulator. Our transcriptome analysis supported by in vivo and in vitro assays showed that C protein directly silenced the expression of the RacR repressor to affect the Rac prophage-related genes. The rac locus ydaST genes, when derepressed, exerted a toxicity indicated by cell filamentation through an unknown mechanism. These results provide an apparent example of transcription factor cross-talk, which can have significant consequences for the host, and may represent a constraint on lateral gene transfer.

Short bZIP homologue of sulfur regulator Met4 from Ogataea parapolymorpha does not depend on DNA-binding cofactors for activating genes in sulfur starvation

The acquisition of sulfur from environment and its assimilation is essential for fungal growth and activities. Here, we describe novel features of the regulatory network of sulfur metabolism in Ogataea parapolymorpha, a thermotolerant methylotrophic yeast with high resistance to harsh environmental conditions. A short bZIP protein (OpMet4p) of O. parapolymorpha, displaying the combined structural characteristics of yeast and filamentous fungal Met4 homologues, plays a key role as a master regulator of cell homeostasis during sulfur limitation, but also its function is required for the tolerance of various stresses. Domain swapping analysis, combined with deletion analysis of the regulatory domains and genes encoding OpCbf1p, OpMet28p, and OpMet32p, indicated that OpMet4p does not require the interaction with these DNA-binding cofactors to induce the expression of sulfur genes, unlike the Saccharomyces cerevisiae Met4p. ChIP analysis confirmed the notion that OpMet4p, which contains a canonical bZIP domain, can bind the target DNA in the absence of cofactors, similar to homologues in other filamentous fungi. Collectively, the identified unique features of the O. parapolymorpha regulatory network, as the first report on the sulfur regulation by a short yeast Met4 homologue, provide insights into conservation and divergence of the sulfur regulatory networks among diverse ascomycetous fungi.

Recent insights into the genotype-phenotype relationship from massively parallel genetic assays

With the molecular revolution in Biology, a mechanistic understanding of the genotype-phenotype relationship became possible. Recently, advances in DNA synthesis and sequencing have enabled the development of deep mutational scanning assays, capable of scoring comprehensive libraries of genotypes for fitness and a variety of phenotypes in massively parallel fashion. The resulting empirical genotype-fitness maps pave the way to predictive models, potentially accelerating our ability to anticipate the behaviour of pathogen and cancerous cell populations from sequencing data. Besides from cellular fitness, phenotypes of direct application in industry (e.g. enzyme activity) and medicine (e.g. antibody binding) can be quantified and even selected directly by these assays. This review discusses the technological basis of and recent developments in massively parallel genetics, along with the trends it is uncovering in the genotype-phenotype relationship (distribution of mutation effects, epistasis), their possible mechanistic bases and future directions for advancing towards the goal of predictive genetics.