Cross-Family Transcription Factor Interactions: An Additional Layer of Gene Regulation

Comprehensive integration of single-cell transcriptomic data illuminates the regulatory network architecture of plant cell fate specification

Plant morphogenesis relies on precise gene expression programs at the proper time and position which is orchestrated by transcription factors (TFs) in intricate regulatory networks in a cell-type specific manner. Here we introduced a comprehensive single-cell transcriptomic atlas of Arabidopsis seedlings. This atlas is the result of meticulous integration of 63 previously published scRNA-seq datasets, addressing batch effects and conserving biological variance. This integration spans a broad spectrum of tissues, including both below- and above-ground parts. Utilizing a rigorous approach for cell type annotation, we identified 47 distinct cell types or states, largely expanding our current view of plant cell compositions. We systematically constructed cell-type specific gene regulatory networks and uncovered key regulators that act in a coordinated manner to control cell-type specific gene expression. Taken together, our study not only offers extensive plant cell atlas exploration that serves as a valuable resource, but also provides molecular insights into gene-regulatory programs that varies from different cell types.

Decoding plant specialized metabolism: new mechanistic insights

Secondary metabolite (SM) production provides biotic and abiotic stress resistance and enables plants to adapt to the environment. Biosynthesis of these metabolites involves a complex interplay between transcription factors (TFs) and regulatory elements, with emerging evidence suggesting an integral role for chromatin dynamics. Here we review key TFs and epigenetic regulators that govern SM biosynthesis in different contexts. We summarize relevant emerging technologies and results from the model species arabidopsis (Arabidopsis thaliana) and outline aspects of regulation that may also function in food, feed, fiber, oil, or industrial crop plants. Finally, we highlight how effective translation of fundamental knowledge from model to non-model species can benefit understanding of SM production in a variety of ecological, agricultural, and pharmaceutical contexts.

Integrative systems biology of wheat susceptibility to Fusarium graminearum uncovers a conserved gene regulatory network and identifies master regulators targeted by fungal core effectors

BACKGROUND

Plant diseases are driven by an intricate set of defense mechanisms counterbalanced by the expression of host susceptibility factors promoted through the action of pathogen effectors. In spite of their central role in the establishment of the pathology, the primary components of plant susceptibility are still poorly understood and challenging to trace especially in plant-fungal interactions such as in Fusarium head blight (FHB) of bread wheat. Designing a system-level transcriptomics approach, we leveraged the analysis of wheat responses from a susceptible cultivar facing Fusarium graminearum strains of different aggressiveness and examined their constancy in four other wheat cultivars also developing FHB.

RESULTS

In this study, we describe unexpected differential expression of a conserved set of transcription factors and an original subset of master regulators were evidenced using a regulation network approach. The dual-integration with the expression data of pathogen effector genes combined with database mining, demonstrated robust connections with the plant molecular regulators and identified relevant candidate genes involved in plant susceptibility, mostly able to suppress plant defense mechanisms. Furthermore, taking advantage of wheat cultivars of contrasting susceptibility levels, a refined list of 142 conserved susceptibility gene candidates was proposed to be necessary host's determinants for the establishment of a compatible interaction.

CONCLUSIONS

Our findings emphasized major FHB determinants potentially controlling a set of conserved responses associated with susceptibility in bread wheat. They provide new clues for improving FHB control in wheat and also could conceivably leverage further original researches dealing with a broader spectrum of plant pathogens.

A meta-analysis of Boolean network models reveals design principles of gene regulatory networks

Gene regulatory networks (GRNs) play a central role in cellular decision-making. Understanding their structure and how it impacts their dynamics constitutes thus a fundamental biological question. GRNs are frequently modeled as Boolean networks, which are intuitive, simple to describe, and can yield qualitative results even when data are sparse. We assembled the largest repository of expert-curated Boolean GRN models. A meta-analysis of this diverse set of models reveals several design principles. GRNs exhibit more canalization, redundancy, and stable dynamics than expected. Moreover, they are enriched for certain recurring network motifs. This raises the important question why evolution favors these design mechanisms.

Identification and characterization of CYC2-like genes related to floral symmetric development in Tagetes erecta (Asteraceae)

Marigold (Tagetes erecta) is an annual herbaceous flower belonging to Asteraceae, whose capitulum is composed of bilateral symmetry ray florets on the outer periphery and radial symmetry disk florets on the inside. The flower symmetry evolution from radial symmetry to bilateral symmetry has changed the morphology, inflorescence architecture and function of florets among several lineages in Asteraceae. Several studies have identified that CYC2 genes in TCP transcription factor family are the key genes regulating the flower morphogenesis, such as corolla symmetry and stamen development. Here, seven TeCYC2 genes were cloned and phylogenetically grouped into the CYC2 branch of TCP transcription family. TeCYC2c and TeCYC2d were found to be expressed specifically in ray florets, TeCYC2b was strongly expressed in both ray and disk florets, TeCYC2g was significantly higher expressed in ray florets than in disk florets, while TeCYC2a, TeCYC2e1 and TeCYC2e2 were significantly expressed in disk florets, according to an examination of the expression profile. Among the ectopic expression lines of seven TeCYC2 genes in Arabidopsis thaliana, the flower symmetry of all transgenic lines was changed from radial symmetry to bilateral symmetry, and only the reproductive growth of TeCYC2c lines was affected. In TeCYC2c transgenic Arabidopsis, the pollen sac was difficult to crack, and the filaments were shorter than the pistils, resulting in a significant decrease in the seed setting rate. All TeCYC2 proteins were localized in the nucleus. Eight pairs of interactions between TeCYC2 proteins were validated by Y2H and BiFC assays, indicating the possibility of TeCYC2 proteins forming homodimers or heterodimers to improve functional specificity. Our findings verified the main regulatory role of TeCYC2c on the development of corollas and stamen in marigold, and analyzed the interaction network of the formation mechanism of floral symmetry in two florets, which provided more insights into the expansion of CYC2 genes in the evolution of Asteraceae inflorescence and contributed to elucidate the complex regulatory network, as well as the molecular breeding concerning flower form diversity in marigold.

Phylogenomics of transcriptionally active AP2/ERF and bHLH transcription factors and study of their promoter regions in Nothapodytes nimmoniana (J.Graham) Mabb

Nothapodytes nimmoniana is a medicinally important plant producing anticancer monoterpene indole alkaloid (MIA), camptothecin (CPT). The CPT is synthesised through the strictosidine intermediate following the MIA pathway; however, transcriptional regulation of CPT pathway is still elusive in N. nimmoniana. Biosynthesis of MIA is regulated by various transcription factors (TFs) belonging to AP2/ERF, bHLH, MYB, and WRKY families. The present study identified transcriptionally active full-length 105 AP2/ERF and 68 bHLH family TFs from the N. nimmoniana. AP2/ERF TFs were divided into three subfamilies along with a soloist, while bHLH TFs were divided into 10 subfamilies according to their phylogenetic similarities. Three group IXa ERFs, Nn-ERF22, Nn-ERF29, and Nn-ERF41, one subfamily IVa TF Nn-bHLH7, and three subfamilies IIIe Nn-bHLH33, Nn-bHLH51, and Nn-bHLH52 clustered with the TFs regulating alkaloid biosynthesis in Catharanthus roseus, tomato, tobacco, and Artemisia annua. Expression of these TFs in N. nimmoniana was higher in roots, which is a primary CPT accumulating tissue. Moreover, genome skimming approach was used to reconstruct the promoter regions of candidate ERF genes to identify the cis-regulatory elements. The presence of G-boxes and other jasmonic acid-responsive elements in the promoter suggests the regulation of ERFs by bHLHs. The present study effectively generated and used genomics resource for characterisation of regulatory TFs from non-model medicinal plant.

Transcription factor bZIP52 modulates Arabidopsis seed oil biosynthesis through interaction with WRINKLED1

Transcriptional regulation mediated by combinatorial interaction of transcription factors (TFs) is a key molecular mechanism modulating plant development and metabolism. Basic leucine zipper (bZIP) TFs play important roles in various plant developmental and physiological processes. However, their involvement in fatty acid biosynthesis is largely unknown. Arabidopsis (Arabidopsis thaliana) WRINKLED1 (WRI1) is a pivotal TF in regulation of plant oil biosynthesis and interacts with other positive and negative regulators. In this study, we identified two bZIP TFs, bZIP21 and bZIP52, as interacting partners of AtWRI1 by yeast-two-hybrid (Y2H)-based screening of an Arabidopsis TF library. We found that coexpression of bZIP52, but not bZIP21, with AtWRI1 reduced AtWRI1-mediated oil biosynthesis in Nicotiana benthamiana leaves. The AtWRI1-bZIP52 interaction was further verified by Y2H, in vitro pull-down, and bimolecular fluorescence complementation assays. Transgenic Arabidopsis plants overexpressing bZIP52 showed reduced seed oil accumulation, while the CRISPR/Cas9-edited bzip52 knockout mutant exhibited increased seed oil accumulation. Further analysis revealed that bZIP52 represses the transcriptional activity of AtWRI1 on the fatty acid biosynthetic gene promoters. Together, our findings suggest that bZIP52 represses fatty acid biosynthesis genes through interaction with AtWRI1, resulting in a reduction of oil production. Our work reports a previously uncharacterized regulatory mechanism that enables fine-tuning of seed oil biosynthesis.

Transcriptional repressor AGL79 positively regulates flowering time in Arabidopsis

The MADS-box gene family is widely distributed in higher plants and the members of the angiosperm-specific APETALA1/FRUITFULL (AP1/FUL) subfamily plays important roles in the regulation of plant reproductive development. Recent studies revealed that the AP1/FUL subfamily member Dt2, VEGETATIVE1/PsFRUITFULc (VEG1/PsFULc) and MtFRUITFULc (MtFULc) are essential for the stem growth, branching and inflorescence development in legume species soybean (Glycine max), pea (Pisum sativum) and Medicago truncatula. However, the biological function of their homologue in Arabidopsis thaliana, AGAMOUS-LIKE 79 (AGL79), has not been well elucidated. In this study, we investigated the developmental roles of Arabidopsis AGL79 by CRISPR/Cas9-mutagenesis and molecular and physiological analyses. We found that AGL79 mainly acts as a transcriptional repressor and positively regulates Arabidopsis flowering time. We further revealed that AGL79 interacts with SUPPRESSOR OF OVEREXPRESSION OF CO1 (SOC1) and represses the expression of TERMINAL FLOWER 1 (TFL1). Our results demonstrated the AGL79-mediated flowering regulation in Arabidopsis and added an additional layer of complexity to the understanding of flowering time regulation in dicot plants.

Kiwifruit bZIP transcription factor AcePosF21 elicits ascorbic acid biosynthesis during cold stress

Cold stress seriously affects plant development, resulting in heavy agricultural losses. L-ascorbic acid (AsA, vitamin C) is an antioxidant implicated in abiotic stress tolerance and metabolism of reactive oxygen species (ROS). Understanding whether and how cold stress elicits AsA biosynthesis to reduce oxidative damage is important for developing cold-resistant plants. Here, we show that the accumulation of AsA in response to cold stress is a common mechanism conserved across the plant kingdom, from single-cell algae to angiosperms. We identified a basic leucine zipper domain (bZIP) transcription factor (TF) of kiwifruit (Actinidia eriantha Benth.), AcePosF21, which was triggered by cold and is involved in the regulation of kiwifruit AsA biosynthesis and defense responses against cold stress. AcePosF21 interacted with the R2R3-MYB TF AceMYB102 and directly bound to the promoter of the gene encoding GDP-L-galactose phosphorylase 3 (AceGGP3), a key conduit for regulating AsA biosynthesis, to up-regulate AceGGP3 expression and produce more AsA, which neutralized the excess ROS induced by cold stress. On the contrary, VIGS or CRISPR-Cas9-mediated editing of AcePosF21 decreased AsA content and increased the generation of ROS in kiwifruit under cold stress. Taken together, we illustrated a model for the regulatory mechanism of AcePosF21-mediated regulation of AceGGP3 expression and AsA biosynthesis to reduce oxidative damage by cold stress, which provides valuable clues for manipulating the cold resistance of kiwifruit.

FaMYB123 interacts with FabHLH3 to regulate the late steps of anthocyanin and flavonol biosynthesis during ripening

In this work, we identified and functionally characterized the strawberry (Fragaria × ananassa) R2R3 MYB transcription factor FaMYB123. As in most genes associated with organoleptic properties of ripe fruit, FaMYB123 expression is ripening-related, receptacle-specific, and antagonistically regulated by ABA and auxin. Knockdown of FaMYB123 expression by RNAi in ripe strawberry fruit receptacles downregulated the expression of enzymes involved in the late steps of anthocyanin/flavonoid biosynthesis. Transgenic fruits showed a parallel decrease in the contents of total anthocyanin and flavonoid, especially malonyl derivatives of pelargonidin and cyanidins. The decrease was concomitant with accumulation of proanthocyanin, propelargonidins, and other condensed tannins associated mainly with green receptacles. Potential coregulation between FaMYB123 and FaMYB10, which may act on different sets of genes for the enzymes involved in anthocyanin production, was explored. FaMYB123 and FabHLH3 were found to interact and to be involved in the transcriptional activation of FaMT1, a gene responsible for the malonylation of anthocyanin components during ripening. Taken together, these results demonstrate that FaMYB123 regulates the late steps of the flavonoid pathway in a specific manner. In this study, a new function for an R2R3 MYB transcription factor, regulating the expression of a gene that encodes a malonyltransferase, has been elucidated.

The protein-protein interaction landscape of transcription factors during gynoecium development in Arabidopsis

Flowers are composed of organs whose identity is defined by the combinatorial activity of transcription factors (TFs). The interactions between MADS-box TFs and protein complex formation have been schematized in the floral quartet model of flower development. The gynoecium is the flower's female reproductive part, crucial for fruit and seed production and, hence, for reproductive success. After the establishment of carpel identity, many tissues arise to form a mature gynoecium. TFs have been described as regulators of gynoecium development, and some interactions and complexes have been identified. However, broad knowledge about the interactions among these TFs and their participation during development remains scarce. In this study, we used a systems biology approach to understand the formation of a complex reproductive unit-as the gynoecium-by mapping binary interactions between well-characterized TFs. We analyzed almost 4500 combinations and detected more than 250 protein-protein interactions (PPIs), resulting in a process-specific interaction map. Topological analyses suggest hidden functions and novel roles for many TFs. In addition, we observed a close relationship between TFs involved in auxin and cytokinin-signaling pathways and other TFs. Furthermore, we analyzed the network by combining PPI data, expression, and genetic data, which helped us to dissect it into several dynamic spatio-temporal subnetworks related to gynoecium development processes. Finally, we generated an extended PPI network that predicts new players in gynoecium development. Taken together, all these results serve as a valuable resource for the plant community.

Expression patterns and correlation analyses of muscle-specific genes in the process of sheep myoblast differentiation

The purpose of this study was to establish a system for the isolation, culture, and differentiation of sheep myoblasts, and to explore the expression patterns as well as mutual relationships of muscle-specific genes. Sheep fetal myoblasts (SFMs) were isolated by two-step enzymatic digestion, purified by differential adhesion and identified using immunofluorescence techniques. Two percent horse serum was used to induce differentiation in SFMs. Real-time quantitative and Western blot analyses were respectively used to detect the mRNA and protein expressions of muscle-specific genes including MyoD, MyoG, Myf5, Myf6, PAX3, PAX7, myomaker, desmin, MYH1, MYH2, MYH4, MYH7, and MSTN during the differentiation of SFMs. Finally, the correlation between muscle-specific genes was analyzed by the Pearson correlation coefficient method. The results showed that the isolated and purified SFMs could form myotubes after the induction for differentiation. The marker factors including MyoD, MyoG, myomaker, desmin, and MyHC were positively stained in SFMs. The mRNA expressions of MyoD, MyoG, and myomaker increased and then decreased, while Myf5, PAX3, and PAX7 decreased; Myf6, desmin, MYH1, MYH2, MYH4, and MYH7 increased; and MSTN fluctuated up and down during the differentiation of SFMs. The expression patterns of protein were basically consistent with those of mRNA except MSTN. There existed significant or highly significant correlations at mRNA or protein level among some genes. Some transcription factor proteins (MyoD, Myf5, Myf6, PAX3, PAX7) showed significant or highly significant correlations with the mRNA level of some other genes and/or themselves. In conclusion, SFMs with good myogenic differentiation ability were successfully isolated, and the expression patterns and correlations of muscle-specific genes during SFM differentiation were revealed, which laid an important foundation for elucidating the gene regulation mechanism of sheep myogenesis.

Transcriptional regulation of oil biosynthesis in seed plants: Current understanding, applications, and perspectives

Plants produce and accumulate triacylglycerol (TAG) in their seeds as an energy reservoir to support the processes of seed germination and seedling development. Plant seed oils are vital not only for the human diet but also as renewable feedstocks for industrial use. TAG biosynthesis consists of two major steps: de novo fatty acid biosynthesis in the plastids and TAG assembly in the endoplasmic reticulum. The latest advances in unraveling transcriptional regulation have shed light on the molecular mechanisms of plant oil biosynthesis. We summarize recent progress in understanding the regulatory mechanisms of well-characterized and newly discovered transcription factors and other types of regulators that control plant fatty acid biosynthesis. The emerging picture shows that plant oil biosynthesis responds to developmental and environmental cues that stimulate a network of interacting transcriptional activators and repressors, which in turn fine-tune the spatiotemporal regulation of the pathway genes.

The TCP transcription factor HvTB2 heterodimerizes with VRS5 and controls spike architecture in barley

Understanding the molecular network, including protein-protein interactions, of VRS5 provide new routes towards the identification of other key regulators of plant architecture in barley. The TCP transcriptional regulator TEOSINTE BRANCHED 1 (TB1) is a key regulator of plant architecture. In barley, an important cereal crop, HvTB1 (also referred to as VULGARE SIX-ROWED spike (VRS) 5), inhibits the outgrowth of side shoots, or tillers, and grains. Despite its key role in barley development, there is limited knowledge on the molecular network that is utilized by VRS5. In this work, we performed protein-protein interaction studies of VRS5. Our analysis shows that VRS5 potentially interacts with a diverse set of proteins, including other class II TCP's, NF-Y TF, but also chromatin remodelers. Zooming in on the interaction capacity of VRS5 with other TCP TFs shows that VRS5 preferably interacts with other class II TCP TFs in the TB1 clade. Induced mutagenesis through CRISPR-Cas of one of the putative VRS5 interactors, HvTB2 (also referred to as COMPOSITUM 1 and BRANCHED AND INDETERMINATE SPIKELET 1), resulted in plants that have lost their characteristic unbranched spike architecture. More specifically, hvtb2 mutants exhibited branches arising at the main spike, suggesting that HvTB2 acts as inhibitor of branching. Our protein-protein interaction studies of VRS5 resulted in the identification of HvTB2 as putative interactor of VRS5, another key regulator of spike architecture in barley. The study presented here provides a first step to underpin the protein-protein interactome of VRS5 and to identify other, yet unknown, key regulators of barley plant architecture.

Transcriptional Regulation of the Chicken ASMT Gene- A Preliminary Analysis

1. Melatonin is an indole hormone that, among its myriad biological functions, regulates circadian and seasonal rhythms in animals. The ASMT gene plays an essential role in melatonin synthesis. However, in chickens, little is known about the regulatory elements governing its transcription.2. The following study identified the transcription start site of the chicken ASMT gene by 5'-RACE. Then, the proximal minimal promoter was identified using a series of 5' truncations of the ASMT promoter (e.g., -3502/+17, -2698/+17, -2003/+17, -1378/+17, and -254/+17). Site-directed mutagenesis, overexpression, and electrophoretic mobility shift assay (EMSA) were applied to show that the transcription factor Oct-1 binds to the promoter region of ASMT.3. The translation start site was located 19 bp upstream from the translational start site. The luciferase reporter assay confirmed that the core promoter of chicken ASMT gene was in the -254/+17 region. Using site-directed mutagenesis, overexpression, and EMSA, Oct-1 bound the promoter of ASMT.4. Overall, Oct1 plays an important role in the transcriptional regulation of chicken ASMT gene.

TCP15 interacts with GOLDEN2-LIKE 1 to control cotyledon opening in Arabidopsis

After germination, exposure to light promotes the opening and expansion of the cotyledons and the development of the photosynthetic apparatus in a process called de-etiolation. This process is crucial for seedling establishment and photoautotrophic growth. TEOSINTE BRANCHED 1, CYCLOIDEA, and PROLIFERATING CELL FACTORS (TCP) transcription factors are important developmental regulators of plant responses to internal and external signals that are grouped into two main classes. In this study, we identified GOLDEN2-LIKE 1 (GLK1), a key transcriptional regulator of photomorphogenesis, as a protein partner of class I TCPs during light-induced cotyledon opening and expansion in Arabidopsis. The class I TCP TCP15 and GLK1 are mutually required for cotyledon opening and the induction of SAUR and EXPANSIN genes, involved in cell expansion. TCP15 also participates in the expression of photosynthesis-associated genes regulated by GLK1, like LHCB1.4 and LHCB2.2. Furthermore, GLK1 and TCP15 bind to the same promoter regions of different target genes containing either GLK or TCP binding motifs and binding of TCP15 is affected in a GLK1-deficient background, suggesting that a complex between TCP15 and GLK1 participates in the induction of these genes. We postulate that GLK1 helps to recruit TCP15 for the modulation of cell expansion genes in cotyledons and that the functional interaction between these transcription factors may serve to coordinate the expression of cell expansion genes with that of genes involved in the development of the photosynthetic apparatus.

The Arabidopsis transcription factors AtPHL1 and AtHB23 act together promoting carbohydrate transport from pedicel-silique nodes to seeds

Carbohydrates are produced in green tissues through photosynthesis and then transported to sink tissues. Carbon partitioning is a strategic process, fine regulated, involving specific sucrose transporters in each connecting tissue. Here we report that a screening of an Arabidopsis transcription factor (TF) library using the homeodomain-leucine zipper I member AtHB23 as bait, allowed identifying the TF AtPHL1 interacting with the former. An independent Y2H assay, and in planta by BiFC, confirmed such interaction. AtHB23 and AtPHL1 coexpressed in the pedicel-silique nodes and the funiculus. Mutant plants (phl1, and amiR23) showed a marked reduction of lipid content in seeds, although lipid composition did not change compared to the wild type. While protein and carbohydrate contents were not significantly different between mutants and control mature seeds, we observed a reduced carbohydrate content in mutant plants young siliques (7 days after pollination). Moreover, using a CFDA probe, we revealed an impaired transport to the seeds, and the gene encoding the carbohydrate transporters SWEET10 and SWEET11, usually expressed in connecting tissues, was repressed in the amiR23 and phl1 mutant plants. Altogether, the results indicated that AtHB23 and AtPHL1 act together, promoting sucrose transport, and the lack of any of them provoked a reduction in seeds lipid content.

Coordinated resource allocation to plant growth-defense tradeoffs

Plant resource allocation patterns often reveal tradeoffs that favor growth (G) over defense (D), or vice versa. Ecologists most often explain G-D tradeoffs through principles of economic optimality, in which negative trait correlations are attributed to the reconciliation of fitness costs. Recently, researchers in molecular biology have developed 'big data' resources including multi-omic (e.g. transcriptomic, proteomic and metabolomic) studies that describe the cellular processes controlling gene expression in model species. In this synthesis, we bridge ecological theory with discoveries in multi-omics biology to better understand how selection has shaped the mechanisms of G-D tradeoffs. Multi-omic studies reveal strategically coordinated patterns in resource allocation that are enabled by phytohormone crosstalk and transcriptional signal cascades. Coordinated resource allocation justifies the framework of optimality theory, while providing mechanistic insight into the feedbacks and control hubs that calibrate G-D tradeoff commitments. We use the existing literature to describe the coordinated resource allocation hypothesis (CoRAH) that accounts for balanced cellular controls during the expression of G-D tradeoffs, while sustaining stored resource pools to buffer the impacts of future stresses. The integrative mechanisms of the CoRAH unify the supply- and demand-side perspectives of previous G-D tradeoff theories.

Both Binding Strength and Evolutionary Accessibility Affect the Population Frequency of Transcription Factor Binding Sequences in Arabidopsis thaliana

Mutations in DNA sequences that bind transcription factors and thus modulate gene expression are a source of adaptive variation in gene expression. To understand how transcription factor binding sequences evolve in natural populations of the thale cress Arabidopsis thaliana, we integrated genomic polymorphism data for loci bound by transcription factors with in vitro data on binding affinity for these transcription factors. Specifically, we studied 19 different transcription factors, and the allele frequencies of 8,333 genomic loci bound in vivo by these transcription factors in 1,135 A. thaliana accessions. We find that transcription factor binding sequences show very low genetic diversity, suggesting that they are subject to purifying selection. High frequency alleles of such binding sequences tend to bind transcription factors strongly. Conversely, alleles that are absent from the population tend to bind them weakly. In addition, alleles with high frequencies also tend to be the endpoints of many accessible evolutionary paths leading to these alleles. We show that both high affinity and high evolutionary accessibility contribute to high allele frequency for at least some transcription factors. Although binding sequences with stronger affinity are more frequent, we did not find them to be associated with higher gene expression levels. Epistatic interactions among individual mutations that alter binding affinity are pervasive and can help explain variation in accessibility among binding sequences. In summary, combining in vitro binding affinity data with in vivo binding sequence data can help understand the forces that affect the evolution of transcription factor binding sequences in natural populations.

A rapid, efficient, and low-cost BiFC protocol and its application in studying in vivo interaction of seed-specific transcription factors, RISBZ and RPBF

Proteins regulate cellular and biological processes in all living organisms. More than 80% of the proteins interact with one another to perform their respective functions; therefore, studying the protein-protein-interaction has gained attention in functional characterization studies. Bimolecular fluorescence complement (BiFC) assay is widely adopted to determine the physical interaction of two proteins in vivo. Here, we developed a simple, yet effective BiFC assay for protein-protein-interaction using transient Agrobacterium-mediated-transformation of onion epidermal cells by taking case study of Rice-P-box-Binding-Factor (RPBF) and rice-seed-specific-bZIP (RISBZ) in vivo interaction. Our result revealed that both the proteins, i.e., RISBZ and RPBF, interacted in the nucleus and cytosol. These two transcription factors are known for their coordinate/synergistic regulation of seed-protein content via concurrent binding to the promoter region of the seed storage protein (SSP) encoding genes. We further validated our results with BiFC assay in Nicotiana by agroinfiltration method, which exhibited similar results as Agrobacterium-mediated-transformation of onion epidermal cells. We also examined the subcellular localization of RISBZ and RPBF to assess the efficacy of the protocol. The subcellular localization and BiFC assay presented here is quite easy-to-follow, reliable, and reproducible, which can be completed within 2-3 days without using costly instruments and technologies that demand a high skill set.

Overexpression of HOXD8 inhibits the proliferation, migration and invasion of breast cancer cells by downregulating ILP2 expression

Breast cancer is one of the most common malignant tumors in women. Although a number of homeobox (HOX) genes are known to serve an important role in breast cancer, the role of HOXD8 in breast cancer remains unclear. The aim of the present study was to investigate the role of HOXD8 in the physiological behaviors of breast cancer cells. The Gene Expression Profiling Interactive Analysis database was used to analyze the expression of HOXD8 in patients with breast cancer and in healthy subjects. Western blotting was performed to determine the expression levels of HOXD8 in several breast cancer cell lines; subsequently, HOXD8 expression was knocked down and overexpressed in MCF-7 cells. Cell Counting Kit-8, colony formation, wound healing and Transwell assays were used to evaluate the effects of HOXD8 on breast cancer cell viability, proliferation, migration and invasion, respectively. Chromatin immunoprecipitation and dual-luciferase reporter assays were conducted to identify the binding sites between HOXD8 and inhibitor of apoptosis-like protein-2 (ILP2). In addition, ILP2 expression levels were knocked down in MCF-7 cells. The results demonstrated that the expression levels of HOXD8 were significantly downregulated in breast cancer tissues and cell lines, and that the overexpression of HOXD8 inhibited the proliferation, invasion and migration of cancer cells. HOXD8 was shown to bind to the ILP2 promoter to regulate the expression of ILP2. Furthermore, ILP2 knockdown reversed the effects of HOXD8 knockdown on breast cancer cell proliferation, invasion and migration. In conclusion, the findings of the present study suggested that HOXD8 may inhibit the proliferation, migration and invasion of breast cancer cells by downregulating ILP2 expression.

The interaction of MYB, bHLH and WD40 transcription factors in red pear (Pyrus pyrifolia) peel

Sunlight enhanced peel color and significantly up-regulated the expression of PyMYB10 and PybHLH genes. MYB-bHLH-WD40 transcriptional complex forms in the light and is involved in regulating anthocyanin accumulation in the peel. Anthocyanin is the major pigment in the peel of Yunnan red pear (Pyrus pyrifolia (Burm.) Nak.). A transcriptional activation protein complex, involving members of the transcription factor classes of MYB, bHLH and WD40, regulates anthocyanin biosynthesis. This complex was examined in the peel of red pear. In order to clarify the interaction of PyMYB10, PybHLH and PyWD40, fruit were bagged then peel samples collected 0, 3, 5, and 7 days after bag removal. Samples were used for Western blotting and protein interaction analysis. The results showed that sunlight enhanced peel color and significantly up-regulated the expression of both PyMYB10 and PybHLH genes. Co-immunoprecipitation (Co-IP) analysis showed that PybHLH interacted with PyMYB10 or PyWD40, and PyMYB10 interacted with PyWD40. Using onion cells as a model system, bimolecular fluorescence complementation (BiFC) confirmed these interactions and showed that the interaction localized to the nuclei. GST Pull down and Far-Western blotting assays demonstrated that PybHLH interacted with PyMYB10 or PyWD40, respectively, and PyMYB10 interacted with PyWD40 in vitro. In addition, EMSA assay showed that PyMYB10 can directly bind to the promoter of the gene encoding the anthocyanin biosynthesis enzyme anthocyanidin synthase (PyANS). Taken together, these results showed that the ternary complex of PyMYB10, PybHLH and PyWD40 transcription factors forms to regulate anthocyanin biosynthesis and accumulation in Yunnan red pear.

Comparative analyses of saprotrophy in Salisapilia sapeloensis and diverse plant pathogenic oomycetes reveal lifestyle-specific gene expression

Oomycetes include many devastating plant pathogens. Across oomycete diversity, plant-infecting lineages are interspersed by non-pathogenic ones. Unfortunately, our understanding of the evolution of lifestyle switches is hampered by a scarcity of data on the molecular biology of saprotrophic oomycetes, ecologically important primary colonizers of dead tissue that can serve as informative reference points for understanding the evolution of pathogens. Here, we established Salisapilia sapeloensis as a tractable system for the study of saprotrophic oomycetes. We generated multiple transcriptomes from S. sapeloensis and compared them with (i) 22 oomycete genomes and (ii) the transcriptomes of eight pathogenic oomycetes grown under 13 conditions. We obtained a global perspective on gene expression signatures of oomycete lifestyles. Our data reveal that oomycete saprotrophs and pathogens use similar molecular mechanisms for colonization but exhibit distinct expression patterns. We identify a S. sapeloensis-specific array and expression of carbohydrate-active enzymes and putative regulatory differences, highlighted by distinct expression levels of transcription factors. Salisapilia sapeloensis expresses only a small repertoire of candidates for virulence-associated genes. Our analyses suggest lifestyle-specific gene regulatory signatures and that, in addition to variation in gene content, shifts in gene regulatory networks underpin the evolution of oomycete lifestyles.

ZmOrphan94 Transcription Factor Downregulates ZmPEPC1 Gene Expression in Maize Bundle Sheath Cells

Spatial separation of the photosynthetic reactions is a key feature of C₄ metabolism. In most C₄ plants, this separation requires compartmentation of photosynthetic enzymes between mesophyll (M) and bundle sheath (BS) cells. The upstream region of the gene encoding the maize PHOSPHOENOLPYRUVATE CARBOXYLASE 1 (ZmPEPC1) has been shown sufficient to drive M-specific ZmPEPC1 gene expression. Although this region has been well characterized, to date, only few trans-factors involved in the ZmPEPC1 gene regulation were identified. Here, using a yeast one-hybrid approach, we have identified three novel maize transcription factors ZmHB87, ZmCPP8, and ZmOrphan94 as binding to the ZmPEPC1 upstream region. Bimolecular fluorescence complementation assays in maize M protoplasts unveiled that ZmOrphan94 forms homodimers and interacts with ZmCPP8 and with two other ZmPEPC1 regulators previously reported, ZmbHLH80 and ZmbHLH90. Trans-activation assays in maize M protoplasts unveiled that ZmHB87 does not have a clear transcriptional activity, whereas ZmCPP8 and ZmOrphan94 act as activator and repressor, respectively. Moreover, we observed that ZmOrphan94 reduces the trans-activation activity of both activators ZmCPP8 and ZmbHLH90. Using the electromobility shift assay, we showed that ZmOrphan94 binds to several cis-elements present in the ZmPEPC1 upstream region and one of these cis-elements overlaps with the ZmbHLH90 binding site. Gene expression analysis revealed that ZmOrphan94 is preferentially expressed in the BS cells, suggesting that ZmOrphan94 is part of a transcriptional regulatory network downregulating ZmPEPC1 transcript level in the BS cells. Based on both this and our previous work, we propose a model underpinning the importance of a regulatory mechanism within BS cells that contributes to the M-specific ZmPEPC1 gene expression.

BMP2-dependent gene regulatory network analysis reveals Klf4 as a novel transcription factor of osteoblast differentiation

Transcription factors (TFs) regulate the expression of target genes, inducing changes in cell morphology or activities needed for cell fate determination and differentiation. The BMP signaling pathway is widely regarded as one of the most important pathways in vertebrate skeletal biology, of which BMP2 is a potent inducer, governing the osteoblast differentiation of bone marrow stromal cells (BMSCs). However, the mechanism by which BMP2 initiates its downstream transcription factor cascade and determines the direction of differentiation remains largely unknown. In this study, we used RNA-seq, ATAC-seq, and animal models to characterize the BMP2-dependent gene regulatory network governing osteoblast lineage commitment. Sp7-Cre; Bmp2fx/fx mice (BMP2-cKO) were generated and exhibited decreased bone density and lower osteoblast number (n > 6). In vitro experiments showed that BMP2-cKO mouse bone marrow stromal cells (mBMSCs) had an impact on osteoblast differentiation and deficient cell proliferation. Osteogenic medium induced mBMSCs from BMP2-cKO mice and control were subjected to RNA-seq and ATAC-seq analysis to reveal differentially expressed TFs, along with their target open chromatin regions. Combined with H3K27Ac CUT&Tag during osteoblast differentiation, we identified 2338 BMP2-dependent osteoblast-specific active enhancers. Motif enrichment assay revealed that over 80% of these elements were directly targeted by RUNX2, DLX5, MEF2C, OASIS, and KLF4. We deactivated Klf4 in the Sp7 + lineage to validate the role of KLF4 in osteoblast differentiation of mBMSCs. Compared to the wild-type, Sp7-Cre; Klf4fx/+ mice (KLF4-Het) were smaller in size and had abnormal incisors resembling BMP2-cKO mice. Additionally, KLF4-Het mice had fewer osteoblasts and decreased osteogenic ability. RNA-seq and ATAC-seq revealed that KLF4 mainly "co-bound" with RUNX2 to regulate downstream genes. Given the significant overlap between KLF4- and BMP2-dependent NFRs and enriched motifs, our findings outline a comprehensive BMP2-dependent gene regulatory network specifically governing osteoblast differentiation of the Sp7 + lineage, in which Klf4 is a novel transcription factor.

Evolution of A bHLH Interaction Motif

Intrinsically disordered proteins and regions with their associated short linear motifs play key roles in transcriptional regulation. The disordered MYC-interaction motif (MIM) mediates interactions between MYC and MYB transcription factors in Arabidopsis thaliana that are critical for constitutive and induced glucosinolate (GLS) biosynthesis. GLSs comprise a class of plant defense compounds that evolved in the ancestor of the Brassicales order. We used a diverse set of search strategies to discover additional occurrences of the MIM in other proteins and in other organisms and evaluate the findings by means of structural predictions, interaction assays, and biophysical experiments. Our search revealed numerous MIM instances spread throughout the angiosperm lineage. Experiments verify that several of the newly discovered MIM-containing proteins interact with MYC TFs. Only hits found within the same transcription factor family and having similar characteristics could be validated, indicating that structural predictions and sequence similarity are good indicators of whether the presence of a MIM mediates interaction. The experimentally validated MIMs are found in organisms outside the Brassicales order, showing that MIM function is broader than regulating GLS biosynthesis.

Determinants of correlated expression of transcription factors and their target genes

While transcription factors (TFs) are known to regulate the expression of their target genes (TGs), only a weak correlation of expression between TFs and their TGs has generally been observed. As lack of correlation could be caused by additional layers of regulation, the overall correlation distribution may hide the presence of a subset of regulatory TF-TG pairs with tight expression coupling. Using reported regulatory pairs in the plant Arabidopsis thaliana along with comprehensive gene expression information and testing a wide array of molecular features, we aimed to discern the molecular determinants of high expression correlation of TFs and their TGs. TF-family assignment, stress-response process involvement, short genomic distances of the TF-binding sites to the transcription start site of their TGs, few required protein-protein-interaction connections to establish physical interactions between the TF and polymerase-II, unambiguous TF-binding motifs, increased numbers of miRNA target-sites in TF-mRNAs, and a young evolutionary age of TGs were found particularly indicative of high TF-TG correlation. The modulating roles of post-transcriptional, post-translational processes, and epigenetic factors have been characterized as well. Our study reveals that regulatory pairs with high expression coupling are associated with specific molecular determinants.

Processing of coding and non-coding RNAs in plant development and environmental responses

Precursor RNAs undergo extensive processing to become mature RNAs. RNA transcripts are subjected to 5′ capping, 3′-end processing, splicing, and modification; they also form dynamic secondary structures during co-transcriptional and post-transcriptional processing. Like coding RNAs, non-coding RNAs (ncRNAs) undergo extensive processing. For example, secondary small interfering RNA (siRNA) transcripts undergo RNA processing, followed by further cleavage to become mature siRNAs. Transcriptome studies have revealed roles for co-transcriptional and post-transcriptional RNA processing in the regulation of gene expression and the coordination of plant development and plant–environment interactions. In this review, we present the latest progress on RNA processing in gene expression and discuss phased siRNAs (phasiRNAs), a kind of germ cell-specific secondary small RNA (sRNA), focusing on their functions in plant development and environmental responses.

The function of the WRI1-TCP4 regulatory module in lipid biosynthesis

The plant-specific TCP transcription factors play pivotal roles in various processes of plant growth and development. However, little is known regarding the functions of TCPs in plant oil biosynthesis. Our recent work showed that TCP4 mediates oil production via interaction with WRINKLED1 (WRI1), an essential transcription factor governing plant fatty acid biosynthesis. Arabidopsis WRI1 (AtWRI1) physically interacts with multiple TCPs, including TCP4, TCP10, and TCP24. Transient co-expression of AtWRI1 with TCP4, but not TCP10 or TCP24, represses oil accumulation in Nicotiana benthamiana leaves. Increased TCP4 in transgenic plants overexpressing a miR319-resistant TCP4 (rTCP4) decreased the expression of AtWRI1 target genes. The tcp4 knockout mutant, the jaw-D mutant with significant reduction of TCP4 expression, and a tcp2 tcp4 tcp10 triple mutant, display increased seed oil contents compared to the wild-type Arabidopsis. The APETALA2 (AP2) transcription factor WRI1 is characterized by regulating fatty acid biosynthesis through cross-family interactions with multiple transcriptional, post-transcriptional, and post-translational regulators. The interacting regulator modules control the range of AtWRI1 transcriptional activity, allowing spatiotemporal modulation of lipid production. Interaction of TCP4 with AtWRI1, which results in a reduction of AtWRI1 activity, represents a newly discovered mechanism that enables the fine-tuning of plant oil biosynthesis.

Multi-Omics Strategies for Decoding Smoke-Assisted Germination Pathways and Seed Vigour

The success of seed germination and the successful establishment of seedlings across diverse environmental conditions depends on seed vigour, which is of both economic and ecologic importance. The smoke-derived exogenous compound karrikins (KARs) and the endogenous plant hormone strigolactone (SL) are two classes of butanolide-containing molecules that follow highly similar signalling pathways to control diverse biological activities in plants. Unravelling the precise mode-of-action of these two classes of molecules in model species has been a key research objective. However, the specific and dynamic expression of biomolecules upon stimulation by these signalling molecules remains largely unknown. Genomic and post-genomic profiling approaches have enabled mining and association studies across the vast genetic diversity and phenotypic plasticity. Here, we review the background of smoke-assisted germination and vigour and the current knowledge of how plants perceive KAR and SL signalling and initiate the crosstalk with the germination-associated hormone pathways. The recent advancement of 'multi-omics' applications are discussed in the context of KAR signalling and with relevance to their adoption for superior agronomic trait development. The remaining challenges and future opportunities for integrating multi-omics datasets associated with their application in KAR-dependent seed germination and abiotic stress tolerance are also discussed.

CmCYC2-like transcription factors may interact with each other or bind to the promoter to regulate floral symmetry development in Chrysanthemum morifolium

The complex capitulum of Chrysanthemum morifolium is often comprised of bilaterally symmetrical ray florets and radially symmetrical disc florets. The TCP transcription factor clade CYCLOIDEA2 (CYC2) appears to play a vital role in determining floral symmetry and in regulating floral organ development in Asteraceae. Our previous study identified six CmCYC2 genes from chrysanthemum and showed that CmCYC2c participated in the regulation of ray floret identity. However, the functions of other CmCYC2 genes and the underlying molecular mechanism of CmCYC2-mediated floral development regulation in chrysanthemums have not been elucidated. In this study, we analysed the function of CmCYC2 genes by ectopic expression of CmCYC2 in Arabidopsis. Then, we examined the protein-protein interaction using yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays. Finally, we analysed the protein-DNA interaction using yeast one-hybrid (Y1H) and dual-luciferase reporter assays. We found that ectopic expression of CmCYC2 genes in the Arabidopsis tcp1 mutant changed its floral symmetry and flowering time. Y2H and BiFC assays confirmed three pairs of interactions between CmCYC2 proteins, that is, CmCYC2b-CmCYC2d, CmCYC2b-CmCYC2e and CmCYC2c-CmCYC2d, suggesting that heterodimeric complexes may form between CmCYC2 proteins to increase their functional specificity. The results of Y1H and dual-luciferase reporter assays indicate that CmCYC2c can bind to the promoter of ClCYC2f. Our findings provided clues that CmCYC2-like transcription factors may interact with each other or bind to the promoter to regulate floral symmetry development in C. morifolium. KEY MESSAGE: CmCYC2-like transcription factors may interact with each other or bind to the promoter to regulate floral symmetry development in Chrysanthemum morifolium.

The barley stripe mosaic virus expression system reveals the wheat C2H2 zinc finger protein TaZFP1B as a key regulator of drought tolerance

BACKGROUND

Drought stress is one of the major factors limiting wheat production globally. Improving drought tolerance is important for agriculture sustainability. Although various morphological, physiological and biochemical responses associated with drought tolerance have been documented, the molecular mechanisms and regulatory genes that are needed to improve drought tolerance in crops require further investigation. We have used a novel 4-component version (for overexpression) and a 3-component version (for underexpression) of a barley stripe mosaic virus-based (BSMV) system for functional characterization of the C2H2-type zinc finger protein TaZFP1B in wheat. These expression systems avoid the need to produce transgenic plant lines and greatly speed up functional gene characterization.

RESULTS

We show that overexpression of TaZFP1B stimulates plant growth and up-regulates different oxidative stress-responsive genes under well-watered conditions. Plants that overexpress TaZFP1B are more drought tolerant at critical periods of the plant's life cycle. Furthermore, RNA-Seq analysis revealed that plants overexpressing TaZFP1B reprogram their transcriptome, resulting in physiological and physical modifications that help wheat to grow and survive under drought stress. In contrast, plants transformed to underexpress TaZFP1B are significantly less tolerant to drought and growth is negatively affected.

CONCLUSIONS

This study clearly shows that the two versions of the BSMV system can be used for fast and efficient functional characterization of genes in crops. The extent of transcriptome reprogramming in plants that overexpress TaZFP1B indicates that the encoded transcription factor is a key regulator of drought tolerance in wheat.

Genome-wide MNase hypersensitivity assay unveils distinct classes of open chromatin associated with H3K27me3 and DNA methylation in Arabidopsis thaliana

BACKGROUND

Regulation of transcription depends on interactions between cis-regulatory elements (CREs) and regulatory proteins. Active CREs are imbedded in open chromatin that are accessible to nucleases. Several techniques, including DNase-seq, which is based on nuclease DNase I, and ATAC-seq, which is based on transposase Tn5, have been widely used to identify genomic regions associated with open chromatin. These techniques have played a key role in dissecting the regulatory networks in gene expression in both animal and plant species.

RESULTS

We develop a technique, named MNase hypersensitivity sequencing (MH-seq), to identify genomic regions associated with open chromatin in Arabidopsis thaliana. Genomic regions enriched with MH-seq reads are referred as MNase hypersensitive sites (MHSs). MHSs overlap with the majority (~ 90%) of the open chromatin identified previously by DNase-seq and ATAC-seq. Surprisingly, 22% MHSs are not covered by DNase-seq or ATAC-seq reads, which are referred to "specific MHSs" (sMHSs). sMHSs tend to be located away from promoters, and a substantial portion of sMHSs are derived from transposable elements. Most interestingly, genomic regions containing sMHSs are enriched with epigenetic marks, including H3K27me3 and DNA methylation. In addition, sMHSs show a number of distinct characteristics including association with transcriptional repressors. Thus, sMHSs span distinct classes of open chromatin that may not be accessible to DNase I or Tn5. We hypothesize that the small size of the MNase enzyme relative to DNase I or Tn5 allows its access to relatively more condensed chromatin domains.

CONCLUSION

MNase can be used to identify open chromatin regions that are not accessible to DNase I or Tn5. Thus, MH-seq provides an important tool to identify and catalog all classes of open chromatin in plants.

Combinatorial interactions of the LEC1 transcription factor specify diverse developmental programs during soybean seed development

The LEAFY COTYLEDON1 (LEC1) transcription factor is a central regulator of seed development, because it controls diverse biological programs during seed development, such as embryo morphogenesis, photosynthesis, and seed maturation. To understand how LEC1 regulates different gene sets during development, we explored the possibility that LEC1 acts in combination with other transcription factors. We identified and compared genes that are directly transcriptionally regulated by ABA-RESPONSIVE ELEMENT BINDING PROTEIN3 (AREB3), BASIC LEUCINE ZIPPER67 (bZIP67), and ABA INSENSITIVE3 (ABI3) with those regulated by LEC1. We showed that LEC1 operates with specific sets of transcription factors to regulate different gene sets and, therefore, distinct developmental processes. Thus, LEC1 controls diverse processes through its combinatorial interactions with other transcription factors. DNA binding sites for the transcription factors are closely clustered in genomic regions upstream of target genes, defining cis-regulatory modules that are enriched for DNA sequence motifs that resemble sequences known to be bound by these transcription factors. Moreover, cis-regulatory modules for genes regulated by distinct transcription factor combinations are enriched for different sets of DNA motifs. Expression assays with embryo cells indicate that the enriched DNA motifs are functional cis elements that regulate transcription. Together, the results suggest that combinatorial interactions between LEC1 and other transcription factors are mediated by cis-regulatory modules containing clustered cis elements and by physical interactions that are documented to occur between the transcription factors.

IDDomainSpotter: Compositional bias reveals domains in long disordered protein regions-Insights from transcription factors

Protein domains constitute regions of distinct structural properties and molecular functions that are retained when removed from the rest of the protein. However, due to the lack of tertiary structure, the identification of domains has been largely neglected for long (>50 residues) intrinsically disordered regions. Here we present a sequence-based approach to assess and visualize domain organization in long intrinsically disordered regions based on compositional sequence biases. An online tool to find putative intrinsically disordered domains (IDDomainSpotter) in any protein sequence or sequence alignment using any particular sequence trait is available at http://www.bio.ku.dk/sbinlab/IDDomainSpotter. Using this tool, we have identified a putative domain enriched in hydrophilic and disorder-promoting residues (Pro, Ser, and Thr) and depleted in positive charges (Arg and Lys) bordering the folded DNA-binding domains of several transcription factors (p53, GCR, NAC46, MYB28, and MYB29). This domain, from two different MYB transcription factors, was characterized biophysically to determine its properties. Our analyses show the domain to be extended, dynamic and highly disordered. It connects the DNA-binding domain to other disordered domains and is present and conserved in several transcription factors from different families and domains of life. This example illustrates the potential of IDDomainSpotter to predict, from sequence alone, putative domains of functional interest in otherwise uncharacterized disordered proteins.

Specificity of MYB interactions relies on motifs in ordered and disordered contexts

Physical interactions between members of the MYB and bHLH transcription factor (TF) families regulate many important biological processes in plants. Not all reported MYB-bHLH interactions can be explained by the known binding sites in the R3 repeat of the MYB DNA-binding domain. Noteworthy, most of the sequence diversity of MYB TFs lies in their non-MYB regions, which contain orphan small subgroup-defining motifs not yet linked to molecular functions. Here, we identified the motif mediating interaction between MYB TFs from subgroup 12 and their bHLH partners. Unlike other known MYB-bHLH interactions, the motif locates to the centre of the predicted disordered non-MYB region. We characterised the core motif, which enabled accurate prediction of previously unknown bHLH-interacting MYB TFs in Arabidopsis thaliana, and we confirmed its functional importance in planta. Our results indicate a correlation between the MYB-bHLH interaction affinity and the phenotypic output controlled by the TF complex. The identification of an interaction motif outside R3 indicates that MYB-bHLH interactions must have arisen multiple times, independently and suggests many more motifs of functional relevance to be harvested from subgroup-specific studies.

R2R3 MYB Transcription Factors - Functions outside the DNA-Binding Domain

Several transcription factor (TF) families, including the MYB family, regulate a wide array of biological processes. TFs contain DNA-binding domains (DBDs) and regulatory regions; although information on protein structure is scarce for plant MYB TFs, various in silico methods suggest that the non-MYB regions contain extensive intrinsically disordered regions (IDRs). Although IDRs do not fold into stable globular structures, they comprise functional regions including interaction motifs, and recent research has shown that IDRs perform crucial biological roles. We map here domain organization, disorder predictions, and functional regions across the entire Arabidopsis thaliana R2R3 MYB TF family, and highlight where an increased research focus will be necessary to shape a new understanding of structure-function relationships in plant TFs.

Specific TCP transcription factors interact with and stabilize PRR2 within different nuclear sub-domains

Plants possess a large set of transcription factors both involved in the control of plant development or in plant stress responses coordination. We previously identified PRR2, a Pseudo-Response Regulator, as a plant-specific CML-interacting partner. We reported that PRR2 acts as a positive actor of plant defense by regulating the production of antimicrobial compounds. Here, we report new data on the interaction between PRR2 and transcription factors belonging to the Teosinte branched Cycloidea and PCF (TCP) family. TCPs have been described to be involved in plant development and immunity. We evaluated the ability of PRR2 to interact with seven TCPs representative of the different subclades of the family. PRR2 is able to interact with TCP13, TCP15, TCP19 and TCP20 in yeast two-hybrid system and in planta interactions were validated for TCP19 and TCP20. Transient expression in tobacco highlighted that PRR2 protein is more easily detected when co-expressed with TCP19 or TC20. This stabilization is associated with a specific sub-nuclear localization of the complex in Cajal bodies or in nuclear speckles according to the interaction of PRR2 with TCP19 or TCP20 respectively. The interaction between PRR2 and TCP19 or TCP20 would contribute to the biological function in specific nuclear compartments.

Transcription factor VqERF114 regulates stilbene synthesis in Chinese wild Vitis quinquangularis by interacting with VqMYB35

VqERF114 regulates stilbene synthesis by interacting with VqMYB35. Resveratrol is a stilbene, an important class of secondary metabolites that accumulates in some plant species, including grapevine. In the plant, these are involved in the response to attack by plant pathogens and, as a component of the human diet, they offer a range of significant health benefits. Stilbene synthase (STS), the key enzyme responsible for resveratrol synthesis, has been characterised in a small number of plant species. However, the regulatory mechanisms for stilbene synthesis are uncertain. Here, an ERF family transcription factor from Chinese wild Vitis quinquangularis, VqERF114, was characterised as an indirect regulator of stilbene synthesis. A transient overexpression assay of VqERF114 in grapevine leaves led to increased STS expression and stilbene accumulation. However, VqERF114 did not bind to the promoters of VqSTSs but the MYB transcription factor, VqMYB35, did interact with VqERF114. This interaction was confirmed by a yeast two-hybrid assay and bimolecular fluorescence complementation. Furthermore, VqMYB35 showed activation effects on the expressions of VqSTS15, VqSTS28, VqSTS42 and VqSTS46 by binding directly to the MBS elements in their promoters. Co-overexpression of VqERF114 and VqMYB35 resulted in higher VqSTSs expression and more stilbene synthesis. These results demonstrate that VqERF114 regulates stilbene synthesis by interacting with VqMYB35.

TF2Network: predicting transcription factor regulators and gene regulatory networks in Arabidopsis using publicly available binding site information

A gene regulatory network (GRN) is a collection of regulatory interactions between transcription factors (TFs) and their target genes. GRNs control different biological processes and have been instrumental to understand the organization and complexity of gene regulation. Although various experimental methods have been used to map GRNs in Arabidopsis thaliana, their limited throughput combined with the large number of TFs makes that for many genes our knowledge about regulating TFs is incomplete. We introduce TF2Network, a tool that exploits the vast amount of TF binding site information and enables the delineation of GRNs by detecting potential regulators for a set of co-expressed or functionally related genes. Validation using two experimental benchmarks reveals that TF2Network predicts the correct regulator in 75–92% of the test sets. Furthermore, our tool is robust to noise in the input gene sets, has a low false discovery rate, and shows a better performance to recover correct regulators compared to other plant tools. TF2Network is accessible through a web interface where GRNs are interactively visualized and annotated with various types of experimental functional information. TF2Network was used to perform systematic functional and regulatory gene annotations, identifying new TFs involved in circadian rhythm and stress response.

The transcriptional activity of a temperature-sensitive transcription factor module is associated with pollen shedding time in pine

It has long been known that the pollen shedding time in pine trees is correlated with temperature, but the molecular basis for this has remained largely unknown. To better understand the mechanisms driving temperature response and to identify the hub regulators of pollen shedding time regulation in Pinus tabuliformis Carr., we identified a set of temperature-sensitive genes by carrying out a comparative transcriptome analysis using six early pollen shedding trees (EPs) and six late pollen shedding trees (LPs) during mid-winter and at three consecutive time points in early spring. We carried out a weighted gene co-expression network analysis and constructed a transcription factor (TF) collaborative network, merging the common but differentially expressed TFs of the EPs and LPs into a joint network. We found five hub genes in the core TF module whose expression was rapidly induced by low temperatures. The transcriptional activity of this TF module was strongly associated with pollen shedding time, and likely to produce the fine balance between cold hardiness and growth activity in early spring. We confirmed the key role of temperature in regulating flowering time and identified a transcription factor module associated with pollen shedding time in P. tabuliformis. This suggests that repression of growth activity by repressors is the main mechanism balancing growth and cold hardiness in pine trees in early spring. Our results provide new insights into the molecular mechanisms regulating seasonal flowering time in pines.

A pepper RING-type E3 ligase, CaASRF1, plays a positive role in drought tolerance via modulation of CaAIBZ1 stability

Plants have evolved complex defense mechanisms to adapt and survive under adverse growth conditions. Abscisic acid (ABA) is a phytohormone that plays a pivotal role in the stress response, especially regulation of the stomatal aperture in response to drought. Here, we identified the pepper CaASRF1 (Capsicum annuum ABA Sensitive RING Finger E3 ligase 1) gene, which modulates drought stress tolerance via ABA-mediated signaling. CaASRF1 contains a C3H2C3-type RING finger domain, which functions as an E3 ligase by attaching ubiquitins to the target proteins. CaASRF1 expression was enhanced after exposure to ABA, drought and NaCl. Loss-of-function in pepper plants and gain-of-function in Arabidopsis plants revealed that CaASRF1 positively modulates ABA signaling and the drought stress response. Moreover, CaASRF1 interacted with and was associated with degradation of the bZIP transcription factor CaAIBZ1 (Capsicum annuum ASRF1-Interacting bZIP transcription factor 1). Contrary to CaASRF1 phenotypes, CaAIBZ1-silenced pepper and CaAIBZ1-overexpressing Arabidopsis exhibited drought-tolerant and drought-sensitive phenotypes, respectively. Taken together, our data indicate that CaASRF1 positively modulates ABA signaling and the drought stress response via modulation of CaAIBZ1 stability.

Gynoecium development: networks in Arabidopsis and beyond

Life has always found a way to preserve itself. One strategy that has been developed for this purpose is sexual reproduction. In land plants, the gynoecium is considered to be at the top of evolutionary innovation, since it has been a key factor in the success of the angiosperms. The gynoecium is composed of carpels with different tissues that need to develop and differentiate in the correct way. In order to control and guide gynoecium development, plants have adapted elements of pre-existing gene regulatory networks (GRNs) but new ones have also evolved. The GRNs can interact with internal factors (e.g. hormones and other metabolites) and external factors (e.g. mechanical signals and temperature) at different levels, giving robustness and flexibility to gynoecium development. Here, we review recent findings regarding the role of cytokinin-auxin crosstalk and the genes that connect these hormonal pathways during early gynoecium development. We also discuss some examples of internal and external factors that can modify GRNs. Finally, we make a journey through the flowering plant lineage to determine how conserved are these GRNs that regulate gynoecium and fruit development.

New roles of NO TRANSMITTING TRACT and SEEDSTICK during medial domain development in Arabidopsis fruits

The gynoecium, the female reproductive part of the flower, is key for plant sexual reproduction. During its development, inner tissues such as the septum and the transmitting tract tissue, important for pollen germination and guidance, are formed. In Arabidopsis, several transcription factors are known to be involved in the development of these tissues. One of them is NO TRANSMITTING TRACT (NTT), essential for transmitting tract formation. We found that the NTT protein can interact with several gynoecium-related transcription factors, including several MADS-box proteins, such as SEEDSTICK (STK), known to specify ovule identity. Evidence suggests that NTT and STK control enzyme and transporter-encoding genes involved in cell wall polysaccharide and lipid distribution in gynoecial medial domain cells. The results indicate that the simultaneous loss of NTT and STK activity affects polysaccharide and lipid deposition and septum fusion, and delays entry of septum cells to their normal degradation program. Furthermore, we identified KAWAK, a direct target of NTT and STK, which is required for the correct formation of fruits in Arabidopsis These findings position NTT and STK as important factors in determining reproductive competence.

Identification of Putative UL54 (ICP27) Transcription Regulatory Sequences Binding to Oct-1, v-Myb, Pax-6 and Hairy in Herpes Simplex Viruses

An oncolytic herpes simplex virus (oHSV) has proven amenable in oncolytic virotherapy and was approved to treat melanoma. The immediate-early (IE) protein ICP27 encoded by gene UL54 is essential for HSV infection. Post-transcriptional modification of UL54 would increase tumor targeting of oHSVs. However, UL54 gene transcription regulatory sequences and factors were not reported yet. Here we isolated a new strain LXMW of type 1 HSV (HSV-1-LXMW) in China and found it's closely related to HSV-1 strains Patton and H129 in the US by the first and next generation DNA sequencing viral DNA phylogenetic analysis. Using a weight matrix-based program Match, we found the UL54 transcription regulatory sequences binding to the transcription factors Oct-1, v-Myb and Pax-6 in HSV-1-LXMW, while the sequences binding to Oct-1 and Hairy in a HSV-2 strain. Further validation showed that HSV-1 and HSV-2 shared the common sequence binding to Oct-1, but had unique sequences to bind v-Myb and Pax-6, or Hairy, respectively, by DNA sequence alignment of total 11 HSV strains. The published results howed that the expression of transcription factors is consistent with the tissue tropism of HSV-1 and HSV-2. In the current article a new HSV-1 strain LXMW was isolated and its putative HSV UL54 transcription regulatory sequences and factors were identified for the first time. Our findings highlight the new understanding of the principles of transcriptional regulation in HSV biology and oncolytic virotherapy.

The SAUR gene family: the plant's toolbox for adaptation of growth and development

The family of small auxin up-regulated RNA (SAUR) genes is a family of auxin-responsive genes with ~60-140 members in most higher plant species. Despite the early discovery of their auxin responsiveness, their function and mode of action remained unknown for a long time. In recent years, the importance of SAUR genes in the regulation of dynamic and adaptive growth, and the molecular mechanisms by which SAUR proteins act are increasingly well understood. SAURs play a central role in auxin-induced acid growth, but can also act independently of auxin, tissue specifically regulated by various other hormone pathways and transcription factors. In this review, we summarize recent advances in the characterization of the SAUR genes in Arabidopsis and other plant species. We particularly elaborate on their capacity to fine-tune growth in response to internal and external signals, and discuss the breakthroughs in understanding the mode of action of SAURs in relation to their complex regulation.

Flavonoid biosynthetic pathways in plants: Versatile targets for metabolic engineering

Plants, fungi, and microorganisms are equipped with biosynthesis machinery for producing thousands of secondary metabolites. These compounds have important functions in nature as a defence against predators or competitors as well as other ecological significances. The full utilization of these compounds for food, medicine, and other purposes requires a thorough understanding of their structures and the distinct biochemical pathways of their production in cellular systems. In this review, flavonoids as classical examples of secondary metabolites are employed to highlight recent advances in understanding how valuable compounds can be regulated at various levels. With extensive diversity in their chemistry and pharmacology, understanding the metabolic engineering of flavonoids now allows us to fine-tune the eliciting of their production, accumulation, and extraction from living systems. More specifically, recent advances in the shikimic acid and acetate biosynthetic pathways of flavonoids production from metabolic engineering point of view, from genes expression to multiple principles of regulation, are addressed. Specific examples of plants and microorganisms as the sources of flavonoids-based compounds with particular emphasis on therapeutic applications are also discussed.