DNA-binding specificities of plant transcription factors and their potential to define target genes

Exploiting Genomic Features to Improve the Prediction of Transcription Factor-Binding Sites in Plants

The identification of transcription factor (TF) target genes is central in biology. A popular approach is based on the location by pattern matching of potential cis-regulatory elements (CREs). During the last few years, tools integrating next-generation sequencing data have been developed to improve the performance of pattern matching. However, such tools have not yet been comprehensively evaluated in plants. Hence, we developed a new streamlined method aiming at predicting CREs and target genes of plant TFs in specific organs or conditions. Our approach implements a supervised machine learning strategy, which allows decision rule models to be learnt using TF ChIP-chip/seq experimental data. Different layers of genomic features were integrated in predictive models: the position on the gene, the DNA sequence conservation, the chromatin state and various CRE footprints. Among the tested features, the chromatin features were crucial for improving the accuracy of the method. Furthermore, we evaluated the transferability of predictive models across TFs, organs and species. Finally, we validated our method by correctly inferring the target genes of key TFs controlling metabolite biosynthesis at the organ level in Arabidopsis. We developed a tool-Wimtrap-to reproduce our approach in plant species and conditions/organs for which ChIP-chip/seq data are available. Wimtrap is a user-friendly R package that supports an R Shiny web interface and is provided with pre-built models that can be used to quickly get predictions of CREs and TF gene targets in different organs or conditions in Arabidopsis thaliana, Solanum lycopersicum, Oryza sativa and Zea mays.

Insights into Four NAC Transcription Factors Involved in Grain Development and in Response to Moderate Heat in the Triticeae Tribe

NAC (NAM (no apical meristem)−ATAF (Arabidopsis transcription activation factor)−CUC (cup-shaped cotyledons)) are among the largest transcription factor families in plants, involved in a plethora of physiological mechanisms. This study focused on four NAC genes previously identified in bread wheat as specifically grain-expressed which could be considered as candidate genes for yield improvement under climate changes. Using in silico analyses, the Triticum aestivum “Grain-NAC” (TaGNAC) orthologs in 14 cereal species were identified. A conserved protein motif was identified only in Triticeae. The expression of TaGNAC and einkorn TmGNAC was studied in response to moderate heat stress during grain development and showed a similar expression pattern that is accelerated during cell division stages under heat stress. A conserved structure was found in the promoter of the Triticeae GNAC orthologs, which is absent in the other Poaceae species. A specific model of promoter structure in Triticeae was proposed, based on the presence of key cis-elements involved in the regulation of seed development, hormonal regulation and response to biotic and abiotic stresses. In conclusion, GNAC genes could play a central role in the regulation of grain development in the Triticeae tribe, particularly in the accumulation of storage proteins, as well as in response to heat stress and could be used as candidate genes for breeding.

Rice NAC17 transcription factor enhances drought tolerance by modulating lignin accumulation

Land plants have developed a comprehensive system to cope with the drought stress, and it is operated by intricate signaling networks, including transcriptional regulation. Herein, we identified the function of OsNAC17, a member of NAC (NAM, ATAF, and CUC2) transcription factor family, in drought tolerance. OsNAC17 is localized to the nucleus, and its expression was significantly induced under drought conditions. A transactivation assay in yeast revealed that the OsNAC17 is a transcriptional activator, harboring an activation domain in the C-terminal region. Overexpressing (OsNAC17^OX) transgenic plants showed drought-tolerant, and knock-out (OsNAC17^KO) plants exhibited drought susceptible phenotype compared to non-transgenic plants. Further investigation revealed that OsNAC17 positively regulates several lignin biosynthetic genes and promotes lignin accumulation in leaves and roots. Together, our results show that OsNAC17 contributes to drought tolerance through lignin biosynthesis in rice.

Genome-Wide Survey and Functional Verification of the NAC Transcription Factor Family in Wild Emmer Wheat

The NAC transcription factor (TF) family is one of the largest TF families in plants, which has been widely reported in rice, maize and common wheat. However, the significance of the NAC TF family in wild emmer wheat (Triticum turgidum ssp. dicoccoides) is not yet well understood. In this study, a genome-wide investigation of NAC genes was conducted in the wild emmer genome and 249 NAC family members (TdNACs) were identified. The results showed that all of these genes contained NAM/NAC-conserved domains and most of them were predicted to be located on the nucleus. Phylogenetic analysis showed that these 249 TdNACs can be classified into seven clades, which are likely to be involved in the regulation of grain protein content, starch synthesis and response to biotic and abiotic stresses. Expression pattern analysis revealed that TdNACs were highly expressed in different wheat tissues such as grain, root, leaves and shoots. We found that TdNAC8470 was phylogenetically close to NAC genes that regulate either grain protein or starch accumulation. Overexpression of TdNAC8470 in rice showed increased grain starch concentration but decreased grain Fe, Zn and Mn contents compared with wild-type plants. Protein interaction analysis indicated that TdNAC8470 might interact with granule-bound starch synthase 1 (TdGBSS1) to regulate grain starch accumulation. Our work provides a comprehensive understanding of the NAC TFs family in wild emmer wheat and establishes the way for future functional analysis and genetic improvement of increasing grain starch content in wheat.

Ectopic Expression of Arabidopsis thaliana zDof1.3 in Tomato (Solanum lycopersicum L.) Is Associated with Improved Greenhouse Productivity and Enhanced Carbon and Nitrogen Use

A large collection of transgenic tomato lines, each ectopically expressing a different Arabidopsis thaliana transcription factor, was screened for variants with alterations in leaf starch. Such lines may be affected in carbon partitioning, and in allocation to the sinks. We focused on ‘L4080’, which harbored an A. thaliana zDof (DNA-binding one zinc finger) isoform 1.3 (AtzDof1.3) gene, and which had a 2−4-fold higher starch-to-sucrose ratio in source leaves over the diel (p < 0.05). Our aim was to determine whether there were associated effects on productivity. L4080 plants were altered in nitrogen (N) and carbon (C) metabolism. The N-to-C ratio was higher in six-week-old L4080, and when treated with 1/10 N, L4080 growth was less inhibited compared to the wild-type and this was accompanied by faster root elongation (p < 0.05). The six-week-old L4080 acquired 42% more dry matter at 720 ppm CO2, compared to ambient CO2 (p < 0.05), while the wild-type (WT) remained unchanged. GC-MS-TOF data showed that L4080 source leaves were enriched in amino acids compared to the WT, and at 49 DPA, fruit had 25% greater mass, higher sucrose, and increased yield (25%; p < 0.05) compared to the WT. An Affymetrix cDNA array analysis suggested that only 0.39% of the 9000 cDNAs were altered by 1.5-fold (p < 0.01) in L4080 source leaves. 14C-labeling of fruit disks identified potential differences in 14-DPA fruit metabolism suggesting that post-transcriptional regulation was important. We conclude that AtzDof1.3 and the germplasm derived therefrom, should be investigated for their ‘climate-change adaptive’ potential.

BES1/BZR1 Family Transcription Factors Regulate Plant Development via Brassinosteroid-Dependent and Independent Pathways

The BES1/BZR1 family is a plant-specific small group of transcription factors possessing a non-canonical bHLH domain. Genetic and biochemical analyses within the last two decades have demonstrated that members of this family are key transcription factors in regulating the expression of brassinosteroid (BR) response genes. Several recent genetic and evolutionary studies, however, have clearly indicated that the BES1/BZR1 family transcription factors also function in regulating several aspects of plant development via BR-independent pathways, suggesting they are not BR specific. In this review, we summarize our current understanding of this family of transcription factors, the mechanisms regulating their activities, DNA binding motifs, and target genes. We selectively discuss a number of their biological functions via BR-dependent and particularly independent pathways, which were recently revealed by loss-of-function genetic analyses. We also highlight a few possible future directions.

The sea-island cotton GbTCP4 transcription factor positively regulates drought and salt stress responses

TCP transcription factors play important regulatory roles in plant growth and development; however, their function in response to salt and drought stress in sea-island cotton (Gossypium barbadense) is unknown. Here, GbTCP4 expression was induced by abscisic acid (ABA), drought, and NaCl treatments. Under drought stress, compared to wild-type (WT) Arabidopsis, transgenic GbTCP4-overexpressing Arabidopsis showed increased seed germination rate, root length and survival rate; additionally, it was ABA-insensitive at the germination stage but ABA-sensitive at the seedling stage, showing reduced stomatal opening and ABA enrichment. Under salt stress, compared to WT Arabidopsis, transgenic GbTCP4-overexpressing Arabidopsis showed greater root length, survival rate, and SPAD value and lower malondialdehyde (MDA) content. Conversely, under drought or salt stress, virus-induced gene-silenced GbTCP4 cotton showed decreased root length, area and volume and increased MDA content and sensitivity to drought and salt stress compared with control cotton. RNA-seq and quantitative real-time PCR analyses showed that GbTCP4 affected the transcription levels of genes across multiple abiotic stress-related metabolic pathways. Furthermore, GbTCP4 activated the transcription of GbUVR8 and GbbHLH130 by binding to their promoters. These results suggest that GbTCP4 positively regulates drought and salt stress responses and is a suitable candidate gene for improving plant drought and salt tolerance.

Plant-DTI: Extending the landscape of TF protein and DNA interaction in plants by a machine learning-based approach

As a sessile organism, plants hold elaborate transcriptional regulatory systems that allow them to adapt to variable surrounding environments. Current understanding of plant regulatory mechanisms is greatly constrained by limited knowledge of transcription factor (TF)-DNA interactions. To mitigate this problem, a Plant-DTI predictor (Plant DBD-TFBS Interaction) was developed here as the first machine-learning model that covered the largest experimental datasets of 30 plant TF families, including 7 plant-specific DNA binding domain (DBD) types, and their transcription factor binding sites (TFBSs). Plant-DTI introduced a novel TFBS feature construction, called TFBS base-preference, which enhanced the specificity of TFBS to DBD types. The proposed model showed better predictive performance with the TFBS base-preference than the simple binary representation. Plant-DTI was validated with 22 independent ChIP-seq datasets. It accurately predicted the measured DBD-TFBS pairs along with their TFBS motifs, and effectively predicted interactions of other TFs containing similar DBD types. Comparing to the existing state-of-art methods, Plant-DTI prediction showed a figure of merit in sensitivity and specificity with respect to the position weight matrix (PWM) and TSPTFBS methods. Finally, the proposed Plant-DTI model helped to fill the knowledge gap in the regulatory mechanisms of the cassava sucrose synthase 1 gene (MeSUS1). Plant-DTI predicted MeERF72 as a regulator of MeSUS1 in consistence with the yeast one-hybrid (Y1H) experiment. Taken together, Plant-DTI would help facilitate the prediction of TF-TFBS and TF-target gene (TG) interactions, thereby accelerating the study of transcriptional regulatory systems in plant species.

SUPPRESSOR OF PHYTOCHROME B-4 #3 reduces the expression of PIF-activated genes and increases expression of growth repressors to regulate hypocotyl elongation in short days

SUPPRESSOR OF PHYTOCHROME B-4 #3 (SOB3) is a member of the AT-HOOK MOTIF CONTAINING NUCLEAR LOCALIZED (AHL) family of transcription factors that are involved in light-mediated growth in Arabidopsis thaliana, affecting processes such as hypocotyl elongation. The majority of the research on the AHLs has been conducted in continuous light. However, there are unique molecular events that promote growth in short days (SD) compared to constant light conditions. Therefore, we investigated how AHLs affect hypocotyl elongation in SD. Firstly, we observed that AHLs inhibit hypocotyl growth in SD, similar to their effect in constant light. Next, we identified AHL-regulated genes in SD-grown seedlings by performing RNA-seq in two sob3 mutants at different time points. Our transcriptomic data indicate that PHYTOCHROME INTERACTING FACTORS (PIFs) 4, 5, 7, and 8 along with PIF-target genes are repressed by SOB3 and/or other AHLs. We also identified PIF target genes that are repressed and have not been previously described as AHL-regulated, including PRE1, PIL1, HFR1, CDF5, and XTR7. Interestingly, our RNA-seq data also suggest that AHLs activate the expression of growth repressors to control hypocotyl elongation, such as HY5 and IAA17. Notably, many growth-regulating and other genes identified from the RNA-seq experiment were differentially regulated between these two sob3 mutants at the time points tested. Surprisingly, our ChIP-seq data suggest that SOB3 mostly binds to similar genes throughout the day. Collectively, these data suggest that AHLs affect gene expression in a time point-specific manner irrespective of changes in binding to DNA throughout SD.

Zinc Finger-Homeodomain and Mini Zinc Finger proteins are key players in plant growth and responses to environmental stresses

The ZINC FINGER-HOMEODOMAIN (ZHD) protein family is a plant-specific family of transcription factors containing two conserved motifs: a non-canonical C5H3 zinc finger domain (ZF) and a DNA-binding homeodomain (HD). The MINI ZINC FINGER (MIF) proteins belong to this family, but were possibly derived from the ZHDs by losing the HD. Information regarding the function of ZHD and MIF proteins is scarce. However, different studies have shown that ZHD/MIF proteins play important roles not only in plant growth and development, but also in response to environmental stresses, including drought and pathogen attack. Here we review recent advances relative to ZHD/MIF functions in multiple species, to provide new insights into the diverse roles of these transcription factors in plants. Their mechanism of action in relation to their ability to interact with other proteins and DNA is also discussed. We then propose directions for future studies to understand better their important roles and pinpoint strategies for potential applications in crop improvement.

Transcription-coupled and epigenome-encoded mechanisms direct H3K4 methylation

Mono-, di-, and trimethylation of histone H3 lysine 4 (H3K4me1/2/3) are associated with transcription, yet it remains controversial whether H3K4me1/2/3 promote or result from transcription. Our previous characterizations of Arabidopsis H3K4 demethylases suggest roles for H3K4me1 in transcription. However, the control of H3K4me1 remains unexplored in Arabidopsis, in which no methyltransferase for H3K4me1 has been identified. Here, we identify three Arabidopsis methyltransferases that direct H3K4me1. Analyses of their genome-wide localization using ChIP-seq and machine learning reveal that one of the enzymes cooperates with the transcription machinery, while the other two are associated with specific histone modifications and DNA sequences. Importantly, these two types of localization patterns are also found for the other H3K4 methyltransferases in Arabidopsis and mice. These results suggest that H3K4me1/2/3 are established and maintained via interplay with transcription as well as inputs from other chromatin features, presumably enabling elaborate gene control.

Genome‑wide identification, phylogenetic and expression pattern analysis of GATA family genes in foxtail millet (Setaria italica)

Background

Transcription factors (TFs) play important roles in plants. Among the major TFs, GATA plays a crucial role in plant development, growth, and stress responses. However, there have been few studies on the GATA gene family in foxtail millet (Setaria italica). The release of the foxtail millet reference genome presents an opportunity for the genome-wide characterization of these GATA genes.

Results

In this study, we identified 28 GATA genes in foxtail millet distributed on seven chromosomes. According to the classification method of GATA members in Arabidopsis, SiGATA was divided into four subfamilies, namely subfamilies I, II, III, and IV. Structural analysis of the SiGATA genes showed that subfamily III had more introns than other subfamilies, and a large number of cis-acting elements were abundant in the promoter region of the SiGATA genes. Three tandem duplications and five segmental duplications were found among SiGATA genes. Tissue-specific results showed that the SiGATA genes were mainly expressed in foxtail millet leaves, followed by peels and seeds. Many genes were significantly induced under the eight abiotic stresses, such as SiGATA10, SiGATA16, SiGATA18, and SiGATA25, which deserve further attention.

Conclusions

Collectively, these findings will be helpful for further in-depth studies of the biological function of SiGATA, and will provide a reference for the future molecular breeding of foxtail millet.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08786-0.

TDTHub, a web server tool for the analysis of transcription factor binding sites in plants

Transcriptional regulation underlies most developmental programs and physiological responses to environmental changes in plants. Transcription factors (TFs) play a key role in the regulation of gene expression by binding specifically to short DNA sequences in the regulatory regions of genes: the TF binding sites (TFBSs). In recent years, several bioinformatic tools have been developed to detect TFBSs in candidate genes, either by de novo prediction or by directly mapping experimentally known TFBSs. However, most of these tools contain information for only a few species or require multi-step procedures, and are not always intuitive for non-experienced researchers. Here we present TFBS-Discovery Tool Hub (TDTHub), a web server for quick and intuitive studies of transcriptional regulation in plants. TDTHub uses pre-computed TFBSs in 40 plant species and allows the choice of two mapping algorithms, providing a higher versatility. Besides the main TFBS enrichment tool, TDTHub includes additional tools to assist in the analysis and visualization of data. In order to demonstrate the effectiveness of TDTHub, we analyzed the transcriptional regulation of the anthocyanin biosynthesis pathway. We also analyzed the transcriptional cascades in response to jasmonate and wounding in Arabidopsis and tomato (Solanum lycopersicum), respectively. In these studies, TDTHub helped to verify the most relevant TF nodes and to propose new ones with a prominent role in these pathways. TDTHub is available at http://acrab.cnb.csic.es/TDTHub/, and it will be periodically upgraded and expanded for new species and gene annotations.

Deciphering the molecular basis of tissue-specific gene expression in plants: Can synthetic biology help?

Gene expression differences between distinct cell types are orchestrated by specific sets of transcription factors and epigenetic regulators acting upon the genome. In plants, the mechanisms underlying tissue-specific gene activity remain largely unexplored. Although transcriptional and epigenetic profiling of individual organs, tissues, and more recently, of single cells can easily detect the molecular signatures of different biological samples, how these unique cell identities are established at the mechanistic level is only beginning to be decoded. Computational methods, including machine learning, used in combination with experimental approaches, enable the identification and validation of candidate cis-regulatory elements driving cell-specific expression. Synthetic biology shows great promise not only as a means of testing candidate DNA motifs but also for establishing the general rules of nature driving promoter architecture and for the rational design of genetic circuits in research and agriculture to confer tissue-specific expression to genes or molecular pathways of interest.

Genome-wide identification and characterization of GATA family genes in wheat

BACKGROUND

Transcription factors GATAs were a member of zinc finger protein, which could bind DNA regulatory regions to control expression of target genes, thus influencing plant growth and development either in normal condition or environmental stresses. Recently, GATA genes have been found and functionally characterized in a number of plant species. However, little information of GATA genes were annotated in wheat.

RESULTS

In the current study, 79 GATA genes were identified in wheat, which were unevenly located on 21 chromosomes. According to the analysis of phylogenetic tree and functional domain structures, TaGATAs were classified into four subfamilies (I, II, III, and IV), consist of 35, 21, 12, and 11 genes, respectively. Meanwhile, the amino acids of 79 TaGATAs exhibited apparent difference in four subfamilies according to GATA domains comparison, gene structures and conserved motif analysis. We then analyze the gene duplication and synteny between the genomes of wheat and Arabidopsis, rice and barley, which provided insights into evolutionary characteristics. In addition, expression patterns of TaGATAs were analyzed, and they showed obvious difference in diverse tissues and abiotic stresses.

CONCLUSION

In general, these results provide useful information for future TaGATA gene function analysis, and it helps to better understand molecular breeding and stress response in wheat.

Computational analysis of potential candidate genes involved in the cold stress response of ten Rosaceae members

Background

Plant species from Rosaceae family are economically important. One of the major environmental factors impacting those species is cold stress. Although several Rosaceae plant genomes have recently been sequenced, there have been very few research conducted on cold upregulated genes and their promoter binding sites. In this study, we used computational approaches to identify and analyse potential cold stress response genes across ten Rosaceae family members.

Results

Cold stress upregulated gene data from apple and strawberry were used to identify syntelogs in other Rosaceae species. Gene duplication analysis was carried out to better understand the distribution of these syntelog genes in different Rosaceae members. A total of 11,145 popular abiotic stress transcription factor-binding sites were identified in the upstream region of these potential cold-responsive genes, which were subsequently categorised into distinct transcription factor (TF) classes. MYB classes of transcription factor binding site (TFBS) were abundant, followed by bHLH, WRKY, and AP2/ERF. TFBS patterns in the promoter regions were compared among these species and gene families, found to be quite different even amongst functionally related syntelogs. A case study on important cold stress responsive transcription factor family, AP2/ERF showed less conservation in TFBS patterns in the promoter regions. This indicates that syntelogs from the same group may be comparable at the gene level but not at the level of cis-regulatory elements. Therefore, for such genes from the same family, different repertoire of TFs could be recruited for regulation and expression. Duplication events must have played a significant role in the similarity of TFBS patterns amongst few syntelogs of closely related species.

Conclusions

Our study overall suggests that, despite being from the same gene family, different combinations of TFs may play a role in their regulation and expression. The findings of this study will provide information about potential genes involved in the cold stress response, which will aid future functional research of these gene families involved in many important biological processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08751-x.

Auxin response factors are keys to the many auxin doors

In plants, most developmental programs depend on the action of auxin. The best described model of the auxin signaling pathway, which explains most, but not all, of the auxin transcriptional responses, relies on a de-repression mechanism. The auxin/indole-3-acetic acid repressors (Aux/IAAs) interact with the auxin response factors (ARFs), the transcription factors of the auxin signaling pathway, leading to repression of the ARF-controlled genes. Auxin induces Aux/IAA degradation, releases ARFs and activates transcription. However, this elegant model is not suitable for all ARFs. Indeed, in Arabidopsis, which has 22 ARFs, only five of them fit into the model since they are the ones able to interact with Aux/IAAs. The remaining 17 have a limited capacity to interact with the repressors, and their mechanisms of action are still unclear. The differential interactions between ARF and Aux/IAA proteins constitute one of many examples of the biochemical and structural diversification of ARFs that affect their action and therefore affect auxin transcriptional responses. A deeper understanding of the structural properties of ARFs is fundamental to obtaining a better explanation of the action of auxin in plants.

Overexpression of a novel heat-inducible ethylene-responsive factor gene LlERF110 from Lilium longiflorum decreases thermotolerance

AP2/ERF (APETALA2/ethylene-responsive factor) family transcription factors are involved in various plant-specific processes, especially in plant development and response to abiotic stress. However, their roles in thermotolerance are still largely unknown. In the current study, we identified a heat-inducible ERF member LlERF110 from Lilium longiflorum that was rapidly induced by high temperature. Its protein was localized in the nucleus, and transcriptional activation activity was observed in yeast and plant cells. In addition, LlERF110 was able to bind to GCC- and CGG-elements, but not to DRE-elements. Overexpression of LlERF110 conferred delayed bolting and bushy phenotype, with decreased thermotolerance accompanied by a disrupted ROS (reactive oxygen species) homeostasis in transgenic plants. The accumulation of LlERF110 may activate certain repressors related to heat stress response (HSR) and indirectly damage the normal expression of heat stress (HS)-protective genes such as AtHSFA2, which consequently leads to reduced thermotolerance. Our results implied that LlERF110 might function as a heat-inducible gene but may hinder the establishment of thermotolerance.

MaMYB4 is a negative regulator and a substrate of RING-type E3 ligases MaBRG2/3 in controlling banana fruit ripening

Fruit ripening is a complex developmental process, which is modulated by both transcriptional and post-translational events. Control of fruit ripening is important in maintaining moderate quality traits and minimizing postharvest deterioration. In this study, we discovered that the transcription factor MaMYB4 acts as a negative regulator of fruit ripening in banana. The protein levels of MaMYB4 decreased gradually with banana fruit ripening, paralleling ethylene production, and decline in firmness. DNA affinity purification sequencing combined with RNA-sequencing analyses showed that MaMYB4 preferentially binds to the promoters of various ripening-associated genes including ethylene biosynthetic and cell wall modifying genes. Furthermore, ectopic expression of MaMYB4 in tomato delayed tomato fruit ripening, which was accompanied by downregulation of ethylene biosynthetic and cell wall modifying genes. Importantly, two RING finger E3 ligases MaBRG2/3, whose protein accumulation increased progressively with fruit ripening, were found to interact with and ubiquitinate MaMYB4, contributing to decreased accumulation of MaMYB4 during fruit ripening. Transient overexpression of MaMYB4 and MaBRG2/3 in banana fruit ripening delayed or promoted fruit ripening by inhibiting or stimulating ethylene biosynthesis, respectively. Taken together, we demonstrate that MaMYB4 negatively modulates banana fruit ripening, and that MaMYB4 abundance could be regulated by protein ubiquitination, thus providing insights into the role of MaMYB4 in controlling fruit ripening at both transcriptional and post-translational levels.

Genome-wide clonal variability in European pear "Rocha" using high-throughput sequencing

Pears (Pyrus) are one of the most economically important fruits worldwide. The Pyrus genus is characterized by a high degree of genetic variability between species and interspecific hybrids, and several studies have been performed to assess this variability for both cultivated and wild accessions. These studies have mostly been limited by the resolving power of traditional molecular markers, although in the recent past the availability of reference genome sequences or SNP arrays for pear have enhanced the capability of high-resolution genomics studies. These tools can also be applied to better understand the intra-varietal (or clonal) variability in pear. Here we report the first high resolution genomics analysis of a pear clonal population using whole genome sequencing (WGS). Results showed unique signatures for the accumulation of mutations and transposable element insertions in each clone, which are likely related to their history of propagation and cultivation. The nucleotide diversity remained low in the clonal collection with the exception of few genomic windows, suggesting that balancing selection may be occurring. These windows included mainly genes related to plant fertility. Regions with higher mutational load were partially associated with transcription factors, probably reflecting the distinctive phenotypes in the collection. The annotation of variants also revealed the theoretical disruption of relevant genes in pear. Taken together, the results from this study show that pear clones accumulate mutations differently, and that those mutations can play a role on pear phenotypes, meaning that the study of pear clonal populations can be relevant in genetic studies, mainly when comparing with traditional association studies.

Genome-Wide Identification and Expression Analysis of Homeodomain Leucine Zipper Subfamily IV (HD-ZIP IV) Gene Family in Cannabis sativa L

The plant-specific homeodomain zipper family (HD-ZIP) of transcription factors plays central roles in regulating plant development and environmental resistance. HD-ZIP transcription factors IV (HDZ IV) have been involved primarily in the regulation of epidermal structure development, such as stomata and trichomes. In our study, we identified nine HDZ IV-encoding genes in Cannabis sativa L. by conducting a computational analysis of cannabis genome resources. Our analysis suggests that these genes putatively encode proteins that have all the conserved domains of HDZ IV transcription factors. The phylogenetic analysis of HDZ IV gene family members of cannabis, rice (Oryza sativa), and Arabidopsis further implies that they might have followed distinct evolutionary paths after divergence from a common ancestor. All the identified cannabis HDZ IV gene promoter sequences have multiple regulation motifs, such as light- and hormone-responsive elements. Furthermore, experimental evidence shows that different HDZ IV genes have different expression patterns in root, stem, leaf, and flower tissues. Four genes were primarily expressed in flowers, and the expression of CsHDG5 (XP_030501222.1) was also correlated with flower maturity. Fifty-nine genes were predicted as targets of HDZ IV transcription factors. Some of these genes play central roles in pathogen response, flower development, and brassinosteroid signaling. A subcellular localization assay indicated that one gene of this family is localized in the Arabidopsis protoplast nucleus. Taken together, our work lays fundamental groundwork to illuminate the function of cannabis HDZ IV genes and their possible future uses in increasing cannabis trichome morphogenesis and secondary metabolite production.

Transcriptional competition shapes proteotoxic ER stress resolution

Through dynamic activities of conserved master transcription factors (mTFs), the unfolded protein response (UPR) relieves proteostasis imbalance of the endoplasmic reticulum (ER), a condition known as ER stress^1,2. Because dysregulated UPR is lethal, the competence for fate changes of the UPR mTFs must be tightly controlled^3,4. However, the molecular mechanisms underlying regulatory dynamics of mTFs remain largely elusive. Here, we identified the abscisic acid-related regulator G-class bZIP TF2 (GBF2) and the cis-regulatory element G-box as regulatory components of the plant UPR led by the mTFs, bZIP28 and bZIP60. We demonstrate that, by competing with the mTFs at G-box, GBF2 represses UPR gene expression. Conversely, a gbf2 null mutation enhances UPR gene expression and suppresses the lethality of a bzip28 bzip60 mutant in unresolved ER stress. By demonstrating that GBF2 functions as a transcriptional repressor of the UPR, we address the long-standing challenge of identifying shared signalling components for a better understanding of the dynamic nature and complexity of stress biology. Furthermore, our results identify a new layer of UPR gene regulation hinged upon an antagonistic mTFs-GFB2 competition for proteostasis and cell fate determination.

A chromosome-level reference genome of Ensete glaucum gives insight into diversity and chromosomal and repetitive sequence evolution in the Musaceae

BACKGROUND

Ensete glaucum (2n = 2x = 18) is a giant herbaceous monocotyledonous plant in the small Musaceae family along with banana (Musa). A high-quality reference genome sequence assembly of E. glaucum is a resource for functional and evolutionary studies of Ensete, Musaceae, and the Zingiberales.

FINDINGS

Using Oxford Nanopore Technologies, chromosome conformation capture (Hi-C), Illumina and RNA survey sequence, supported by molecular cytogenetics, we report a high-quality 481.5 Mb genome assembly with 9 pseudo-chromosomes and 36,836 genes. A total of 55% of the genome is composed of repetitive sequences with predominantly LTR-retroelements (37%) and DNA transposons (7%). The single 5S ribosomal DNA locus had an exceptionally long monomer length of 1,056 bp, more than twice that of the monomers at multiple loci in Musa. A tandemly repeated satellite (1.1% of the genome, with no similar sequence in Musa) was present around all centromeres, together with a few copies of a long interspersed nuclear element (LINE) retroelement. The assembly enabled us to characterize in detail the chromosomal rearrangements occurring between E. glaucum and the x = 11 species of Musa. One E. glaucum chromosome has the same gene content as Musa acuminata, while others show multiple, complex, but clearly defined evolutionary rearrangements in the change between x= 9 and 11.

CONCLUSIONS

The advance towards a Musaceae pangenome including E. glaucum, tolerant of extreme environments, makes a complete set of gene alleles, copy number variation, and a reference for structural variation available for crop breeding and understanding environmental responses. The chromosome-scale genome assembly shows the nature of chromosomal fusion and translocation events during speciation, and features of rapid repetitive DNA change in terms of copy number, sequence, and genomic location, critical to understanding its role in diversity and evolution.

OsOSCA1.1 Mediates Hyperosmolality and Salt Stress Sensing in Oryza sativa

OSCA (reduced hyperosmolality-induced [Ca²⁺]_i increase) is a family of mechanosensitive calcium-permeable channels that play a role in osmosensing and stomatal immunity in plants. Oryza sativa has 11 OsOSCA genes; some of these were shown to complement hyperosmolality-induced [Ca²⁺]_cyt increases (OICI_cyt), salt stress-induced [Ca²⁺]_cyt increases (SICI_cyt), and the associated growth phenotype in the Arabidopsis thaliana mutant osca1. However, their biological functions in rice remain unclear. In this paper, we found that OsOSCA1.1 mediates OICI_cyt and SICI_cyt in rice roots, which are critical for stomatal closure, plant survival, and gene expression in shoots, in response to hyperosmolality and the salt stress treatment of roots. Compared with wild-type (Zhonghua11, ZH11) plants, OICI_cyt and SICI_cyt were abolished in the roots of 10-day-old ososca1.1 seedlings, in response to treatment with 250 mM of sorbitol and 100 mM of NaCl, respectively. Moreover, hyperosmolality- and salt stress-induced stomatal closure were also disrupted in a 30-day-old ososca1.1 mutant, resulting in lower stomatal resistance and survival rates than that in ZH11. However, overexpression of OsOSCA1.1 in ososca1.1 complemented stomatal movement and survival, in response to hyperosmolality and salt stress. The transcriptomic analysis further revealed the following three types of OsOSCA1.1-regulated genes in the shoots: 2416 sorbitol-responsive, 2349 NaCl-responsive and 1844 common osmotic stress-responsive genes after treated with 250 mM of sorbitol and 125 mM NaCl of in 30-day-old rice roots for 24 h. The Gene Ontology enrichment analysis showed that these OsOSCA1.1-regulated genes were relatively enriched in transcription regulation, hormone response, and phosphorylation terms of the biological processes category, which is consistent with the Cis-regulatory elements ABRE, ARE, MYB and MYC binding motifs that were overrepresented in 2000-bp promoter regions of these OsOSCA1.1-regulated genes. These results indicate that OsOSCA-mediated calcium signaling specifically regulates gene expression, in response to drought and salt stress in rice.

Transcriptomic Data Meta-Analysis Sheds Light on High Light Response in Arabidopsis thaliana L

The availability and intensity of sunlight are among the major factors of growth, development and metabolism in plants. However, excessive illumination disrupts the electronic balance of photosystems and leads to the accumulation of reactive oxygen species in chloroplasts, further mediating several regulatory mechanisms at the subcellular, genetic, and molecular levels. We carried out a comprehensive bioinformatic analysis that aimed to identify genetic systems and candidate transcription factors involved in the response to high light stress in Arabidopsis thaliana L. using resources GEO NCBI, string-db, ShinyGO, STREME, and Tomtom, as well as programs metaRE, CisCross, and Cytoscape. Through the meta-analysis of five transcriptomic experiments, we selected a set of 1151 differentially expressed genes, including 453 genes that compose the gene network. Ten significantly enriched regulatory motifs for TFs families ZF-HD, HB, C2H2, NAC, BZR, and ARID were found in the promoter regions of differentially expressed genes. In addition, we predicted families of transcription factors associated with the duration of exposure (RAV, HSF), intensity of high light treatment (MYB, REM), and the direction of gene expression change (HSF, S1Fa-like). We predicted genetic components systems involved in a high light response and their expression changes, potential transcriptional regulators, and associated processes.

Plant circadian clock control of Medicago truncatula nodulation via regulation of nodule cysteine-rich peptides

Legumes house nitrogen-fixing endosymbiotic rhizobia in specialized polyploid cells within root nodules, which undergo tightly regulated metabolic activity. By carrying out expression analysis of transcripts over time in Medicago truncatula nodules, we found that the circadian clock enables coordinated control of metabolic and regulatory processes linked to nitrogen fixation. This involves the circadian clock-associated transcription factor LATE ELONGATED HYPOCOTYL (LHY), with lhy mutants being affected in nodulation. Rhythmic transcripts in root nodules include a subset of nodule-specific cysteine-rich peptides (NCRs) that have the LHY-bound conserved evening element in their promoters. Until now, studies have suggested that NCRs act to regulate bacteroid differentiation and keep the rhizobial population in check. However, these conclusions came from the study of a few members of this very large gene family that has complex diversified spatio-temporal expression. We suggest that rhythmic expression of NCRs may be important for temporal coordination of bacterial activity with the rhythms of the plant host, in order to ensure optimal symbiosis.

BBX16 mediates the repression of seedling photomorphogenesis downstream of the GUN1/GLK1 module during retrograde signalling

Plastid-to-nucleus retrograde signalling (RS) initiated by dysfunctional chloroplasts impact photomorphogenic development. We have previously shown that the transcription factor GLK1 acts downstream of the RS regulator GUN1 in photodamaging conditions to regulate not only the well established expression of photosynthesis-associated nuclear genes (PhANGs) but also to regulate seedling morphogenesis. Specifically, the GUN1/GLK1 module inhibits the light-induced phytochrome-interacting factor (PIF)-repressed transcriptional network to suppress cotyledon development when chloroplast integrity is compromised, modulating the area exposed to potentially damaging high light. However, how the GUN1/GLK1 module inhibits photomorphogenesis upon chloroplast damage remained undefined. Here, we report the identification of BBX16 as a novel direct target of GLK1. BBX16 is induced and promotes photomorphogenesis in moderate light and is repressed via GUN1/GLK1 after chloroplast damage. Additionally, we showed that BBX16 represents a regulatory branching point downstream of GUN1/GLK1 in the regulation of PhANG expression and seedling development upon RS activation. The gun1 phenotype in lincomycin and the gun1-like phenotype of GLK1OX are markedly suppressed in gun1bbx16 and GLK1OXbbx16. This study identified BBX16 as the first member of the BBX family involved in RS, and defines a molecular bifurcation mechanism operated by GLK1/BBX16 to optimise seedling de-etiolation, and to ensure photoprotection in unfavourable light conditions.

The SEEL motif and members of the MYB-related REVEILLE transcription factor family are important for the expression of LORELEI in the synergid cells of the Arabidopsis female gametophyte

Synergid cells in the micropylar end of the female gametophyte are required for critical cell-cell signaling interactions between the pollen tube and the ovule that precede double fertilization and seed formation in flowering plants. LORELEI (LRE) encodes a putative GPI-anchored protein that is expressed primarily in the synergid cells, and together with FERONIA, a receptor-like kinase, it controls pollen tube reception by the receptive synergid cell. Still, how LRE expression is controlled in synergid cells remains poorly characterized. We identified candidate cis-regulatory elements enriched in LRE and other synergid cell-expressed genes. One of the candidate motifs ('TAATATCT') in the LRE promoter was an uncharacterized variant of the Evening Element motif that we named as the Short Evening Element-like (SEEL) motif. Deletion or point mutations in the SEEL motif of the LRE promoter resulted in decreased reporter expression in synergid cells, demonstrating that the SEEL motif is important for expression of LRE in synergid cells. Additionally, we found that LRE expression is decreased in the loss of function mutants of REVEILLE (RVE) transcription factors, which are clock genes known to bind the SEEL and other closely related motifs. We propose that RVE transcription factors regulate LRE expression in synergid cells by binding to the SEEL motif in the LRE promoter. Identification of cis-regulatory elements and transcription factors involved in the expression of LRE will serve as a foundation to characterize the gene regulatory networks in synergid cells.

Cis-regulatory sequences in plants: Their importance, discovery, and future challenges

The identification and characterization of cis-regulatory DNA sequences and how they function to coordinate responses to developmental and environmental cues is of paramount importance to plant biology. Key to these regulatory processes are cis-regulatory modules (CRMs), which include enhancers and silencers. Despite the extraordinary advances in high-quality sequence assemblies and genome annotations, the identification and understanding of CRMs, and how they regulate gene expression, lag significantly behind. This is especially true for their distinguishing characteristics and activity states. Here, we review the current knowledge on CRMs and breakthrough technologies enabling identification, characterization, and validation of CRMs; we compare the genomic distributions of CRMs with respect to their target genes between different plant species, and discuss the role of transposable elements harboring CRMs in the evolution of gene expression. This is an exciting time to study cis-regulomes in plants; however, significant existing challenges need to be overcome to fully understand and appreciate the role of CRMs in plant biology and in crop improvement.

Genome-wide identification of C2H2-type zinc finger gene family members and their expression during abiotic stress responses in orchardgrass (Dactylis glomerata)

The C2H2-type zinc finger protein (ZFP) family is one of the largest transcription factor families in the plant kingdom and its members are involved in plant growth, development, and stress responses. As an economically valuable perennial graminaceous forage crop, orchardgrass (Dactylis glomerata) is an important feedstuff resource owing to its high yield and quality. In this study, 125 C2H2-type ZFPs in orchardgrass (Dg-ZFPs) were identified and further classified by phylogenetic analysis. The members with similar gene structures were generally clustered into the same groups, with proteins containing the conserved QALGGH motif being concentrated in groups VIII and IX. Gene ontology and miRNA target analyses indicated that Dg-ZFPs likely perform diverse biological functions through their gene interactions. The RNA-seq data revealed differentially expressed genes across tissues and development phases, suggesting that some Dg-ZFPs might participate in growth and development regulation. Abiotic stress responses of Dg-ZFP genes were verified by qPCR and Saccharomyces cerevisiae transformation, revealing that Dg-ZFP125 could enhance the tolerance of yeasts to osmotic and salt stresses. Our study performed a novel systematic analysis of Dg-ZFPs in orchardgrass, providing a reference for this gene family in other grasses and revealing new insights for enhancing gene utilization.

Plant transcription factors - being in the right place with the right company

Transcriptional regulation underlies many of the growth and developmental processes that shape plants as well as their adaptation to their environment. Key to transcriptional control are transcription factors, DNA-binding proteins that serve two essential functions: to find the appropriate DNA contact sites in their target genes; and to recruit other proteins to execute transcriptional transactions. In recent years, protein structural, genomic, bioinformatic, and proteomic analyses have led to new insights into how these central functions are regulated. Here, we review new findings relating to plant transcription factor function and to their role in shaping transcription in the context of chromatin.

How do plants defend themselves against pathogens-Biochemical mechanisms and genetic interventions

In agro-ecosystem, plant pathogens hamper food quality, crop yield, and global food security. Manipulation of naturally occurring defense mechanisms in host plants is an effective and sustainable approach for plant disease management. Various natural compounds, ranging from cell wall components to metabolic enzymes have been reported to protect plants from infection by pathogens and hence provide specific resistance to hosts against pathogens, termed as induced resistance. It involves various biochemical components, that play an important role in molecular and cellular signaling events occurring either before (elicitation) or after pathogen infection. The induction of reactive oxygen species, activation of defensive machinery of plants comprising of enzymatic and non-enzymatic antioxidative components, secondary metabolites, pathogenesis-related protein expression (e.g. chitinases and glucanases), phytoalexin production, modification in cell wall composition, melatonin production, carotenoids accumulation, and altered activity of polyamines are major induced changes in host plants during pathogen infection. Hence, the altered concentration of biochemical components in host plants restricts disease development. Such biochemical or metabolic markers can be harnessed for the development of "pathogen-proof" plants. Effective utilization of the key metabolites-based metabolic markers can pave the path for candidate gene identification. This present review discusses the valuable information for understanding the biochemical response mechanism of plants to cope with pathogens and genomics-metabolomics-based sustainable development of pathogen proof cultivars along with knowledge gaps and future perspectives to enhance sustainable agricultural production.

Pseudomonas syringae activates ZAT18 to inhibit salicylic acid accumulation by repressing EDS1 transcription for bacterial infection

Phytopathogens can manipulate plant hormone signaling to counteract immune responses; however, the underlying mechanism is mostly unclear. Here, we report that Pseudomonas syringae pv tomato (Pst) DC3000 induces expression of C2H2 zinc finger transcription factor ZAT18 in a jasmonic acid (JA)-signaling-dependent manner. Biochemical assays further confirmed that ZAT18 is a direct target of MYC2, which is a very important regulator in JA signaling. CRISPR/Cas9-generated zat18-cr mutants exhibited enhanced resistance to Pst DC3000, while overexpression of ZAT18 resulted in impaired disease resistance. Genetic characterization of ZAT18 mutants demonstrated that ZAT18 represses defense responses by inhibiting the accumulation of the key plant immune signaling molecule salicylic acid (SA), which is dependent on its EAR motif. ZAT18 exerted this inhibitory effect by directly repressing the transcription of Enhanced Disease Susceptibility 1 (EDS1), which is the key signaling component of pathogen-induced SA accumulation. Overexpression of ZAT18 resulted in decreased SA content, while loss of function of ZAT18 showed enhanced SA accumulation upon pathogen infection. Furthermore, enhanced resistance and SA content in zat18-cr mutants was abolished by the mutation in EDS1. Our data indicate that pathogens induce ZAT18 expression to repress the transcription of EDS1, further antagonising SA accumulation for bacterial infection.

Prediction of conserved and variable heat and cold stress response in maize using cis-regulatory information

Changes in gene expression are important for responses to abiotic stress. Transcriptome profiling of heat- or cold-stressed maize genotypes identifies many changes in transcript abundance. We used comparisons of expression responses in multiple genotypes to identify alleles with variable responses to heat or cold stress and to distinguish examples of cis- or trans-regulatory variation for stress-responsive expression changes. We used motifs enriched near the transcription start sites (TSSs) for thermal stress-responsive genes to develop predictive models of gene expression responses. Prediction accuracies can be improved by focusing only on motifs within unmethylated regions near the TSS and vary for genes with different dynamic responses to stress. Models trained on expression responses in a single genotype and promoter sequences provided lower performance when applied to other genotypes but this could be improved by using models trained on data from all three genotypes tested. The analysis of genes with cis-regulatory variation provides evidence for structural variants that result in presence/absence of transcription factor binding sites in creating variable responses. This study provides insights into cis-regulatory motifs for heat- and cold-responsive gene expression and defines a framework for developing models to predict expression responses across multiple genotypes.

BUB1B (BUB1 Mitotic Checkpoint Serine/Threonine Kinase B) Promotes Lung Adenocarcinoma by Interacting with Zinc Finger Protein ZNF143 and Regulating Glycolysis

Lung adenocarcinoma (LUAD) is one of the most common causes of cancer death in men. BUB1B (BUB1 mitotic checkpoint serine/threonine kinase B) has been reported to contribute to the initiation and development of several cancers. Here, we aimed to explore the potential role of BUB1B in LUAD. We found BUB1B was upregulated in LUAD, suggesting its potential role as a biomarker for LUAD diagnosis. Significantly, LUAD patients with high BUB1B expression had a shorter survival time than those with low BUB1B expression. Knocking-out BUB1B resulted in suppression of cell proliferation, migration, and invasion in vitro, and inhibition of tumor growth in the xenograft experiment. Further analysis revealed that BUB1B regulates glycolysis in LUAD and interacting with ZNF143 in LUAD cells. The interaction was demonstrated by silencing ZNF143, which led to a decrease in proliferation, migration, and invasion in LUAD cells, whereas overexpressing BUB1B had the opposite effects. Our study suggested that the ZNF143/BUB1B axis plays a pivotal role in LUAD progression, which might be a potential target for LUAD management.

Genome-Wide Identification of APETALA2/ETHYLENE RESPONSIVE FACTOR Transcription Factors in Cucurbita moschata and Their Involvement in Ethylene Response

APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF), a plant-specific transcription factor (TF) family, plays an essential role in the growth and development of plants, and in their response to biotic and abiotic stresses. However, information on AP2/ERF in Cucurbita moschata (pumpkin), an edible and medicinal vegetable used worldwide, is scarce. A total of 212 AP2/ERF genes were identified in the C. moschata genome (CmoAP2/ERFs). Based on phylogenetic analysis, they were divided into four groups-28 AP2s, 92 ERFs, 86 dehydration-responsive element-binding (DREB) factors, and 6 ABI3/VPs (RAV). The 212 AP2/ERF genes were unevenly distributed on the 20 chromosomes of C. moschata. The results of structural analysis showed the absence of introns on 132 CmoAP2/ERFs. Four pairs of tandem duplication and 155 pairs of segmental duplication events were identified, which indicated that segmental duplications might be the main reason for the expansion of the CmoAP2/ERF family. The analysis of cis-regulatory elements (CREs) showed that most of the CmoAP2/ERFs contained hormone response elements (ABREs, EREs) in their promoters, suggesting that AP2/ERFs could contribute to the processes regulated by ethylene and abscisic acid. By comparing the transcriptome of ethephon-treated and control plants, we found that 16 CmoAP2/ERFs were significantly upregulated after ethephon treatment. Furthermore, we determined the expression patterns of these genes at different developmental stages of female and male flowers. This study provides insights into the identification, classification, physicochemical property, phylogenetic analysis, chromosomal location, gene structure, motif identification, and CRE prediction of the AP2/ERF superfamily in C. moschata. Sixteen CmoAP2/ERF genes were identified as ethylene-inducible genes. The results of this study will be valuable for understanding the roles of CmoAP2/ERFs in ethylene response and should provide a foundation for elucidating the function of AP2/ERF TFs in C. moschata.

A Cys2His2 Zinc Finger Transcription Factor BpSZA1 Positively Modulates Salt Stress in Betula platyphylla

Zinc finger proteins (ZFPs) are widely involved in plant growth and abiotic stress responses, however, few of these proteins have been functionally characterized in tree species. In this study, we cloned and characterized the BpSZA1 gene encoding a C2H2-type ZFP from Betula platyphylla. BpSZA1 is a transcription factor localized in the nucleus, with a transcription activation domain located at the N-terminus. BpSZA1 was predominantly expressed in stems and was induced by salt. We generated transgenic birch lines displaying overexpression (OE) or RNAi silencing (Ri) of BpSZA1 and exposed these along with wild-type birch seedlings to salinity. Phenotypic and physiological parameters such as superoxide dismutase, peroxisome, H₂O₂ content, proline content, water loss rate, and malondialdehyde content were examined. Overexpression of BpSZA1 in birch conferred increased salt tolerance. Chromatin immunoprecipitation-qPCR and RNA-seq showed that BpSZA1 binds to the GAGA-motif in the promoter of downstream target genes including BpAPX1, BpAPX2, BpCAT, and Bp6PGDH to activate their transcription. BpSZA1 also participates in abscisic acid (ABA) biosynthesis, proline biosynthesis, and the ABA/jasmonic acid pathway to enhance the salt stress of B. platyphylla.

The interface of central metabolism with hormone signaling in plants

Amongst the myriad of metabolites produced by plants, primary metabolites and hormones play crucial housekeeping roles in the cell and are essential for proper plant growth and development. While the biosynthetic pathways of primary metabolism are well characterized, those of hormones are yet to be completely defined. Central metabolism provides precursors for hormone biosynthesis and the regulation and function of primary metabolites and hormones are tightly entwined. The combination of reverse genetics and technological advances in our ability to evaluate the levels of the molecular entities of the cell (transcripts, proteins and metabolites) has led to considerable improvements in our understanding of both the regulatory interaction between primary metabolites and hormones and its coordination in response to different conditions. Here, we provide an overview of the interaction of primary and hormone metabolism at the metabolic and signaling levels, as well as a perspective regarding the tools that can be used to tackle our current knowledge gaps at the signaling level.

NIN-like protein7 and PROTEOLYSIS6 functional interaction enhances tolerance to sucrose, ABA, and submergence

Nitrate (NO3) assimilation and signaling regulate plant growth through the relevant function of the transcription factor NIN-like Protein7 (NLP7). NO3 is also the main source for plants to produce nitric oxide (NO), which regulates growth and stress responses. NO-mediated regulation requires efficient sensing via the PROTEOLYSIS6 (PRT6)-mediated proteasome-triggered degradation of group VII of ethylene response transcription factors through the Cys/Arg N-degron pathway. The convergence of NO3 signaling and N-degron proteolysis on NO-mediated regulation remains largely unknown. Here, we investigated the functional interaction between NLP7 and PRT6 using Arabidopsis (Arabidopsis thaliana) double prt6 nlp7 mutant plants as well as complementation lines overexpressing NLP7 in different mutant genetic backgrounds. prt6 nlp7 mutant plants displayed several potentiated prt6 characteristic phenotypes, including slower vegetative growth, increased NO content, and diminished tolerance to abiotic stresses such as high-sucrose concentration, abscisic acid, and hypoxia-reoxygenation. Although NLP7 has an N-terminus that could be targeted by the N-degron proteolytic pathway, it was not a PRT6 substrate. The potential PRT6- and NO-regulated nucleocytoplasmic translocation of NLP7, which is likely modulated by posttranslational modifications, is proposed to act as a regulatory loop to control NO homeostasis and action.

The BTB/POZ domain protein GmBTB/POZ promotes the ubiquitination and degradation of the soybean AP2/ERF-like transcription factor GmAP2 to regulate the defense response to Phytophthora sojae

Phytophthora root and stem rot in soybean (Glycine max) is a destructive disease worldwide, and hence improving crop resistance to the causal pathogen, P. sojae, is a major target for breeders. However, it remains largely unclear how the pathogen regulates the various affected signaling pathways in the host, which consist of complex networks including key transcription factors and their targets. We have previously demonstrated that GmBTB/POZ enhances soybean resistance to P. sojae and the associated defense response. Here, we demonstrate that GmBTB/POZ interacts with the transcription factor GmAP2 and promotes its ubiquitination. GmAP2-RNAi transgenic soybean hairy roots exhibited enhanced resistance to P. sojae, whereas roots overexpressing GmAP2 showed hypersensitivity. GmWRKY33 was identified as a target of GmAP2, which represses its expression by directly binding to the promoter. GmWRKY33 acts as a positive regulator in the response of soybean to P. sojae. Overexpression of GmBTB/POZ released the GmAP2-regulated suppression of GmWRKY33 in hairy roots overexpressing GmAP2 and increased their resistance to P. sojae. Taken together, our results indicate that GmBTB/POZ-GmAP2 modulation of the P. sojae resistance response forms a novel regulatory mechanism, which putatively regulates the downstream target gene GmWRKY33 in soybean.

Characterization of Chromatin Accessibility and Gene Expression upon Cold Stress Reveals that the RAV1 Transcription Factor Functions in Cold Response in Vitis Amurensis

Cold tolerance is regulated by a variety of transcription factors (TFs) and their target genes. Except for the well-characterized C-repeat binding factors (CBFs)-dependent transcriptional cascade, the mechanisms of cold tolerance mediated by other transcriptional regulatory networks are still largely unknown. Here, we used the assay for transposase-accessible chromatin with sequencing (ATAC-seq) and RNA-seq to identify cold responsive TFs in Vitis amurensis, a grape species with high cold hardiness. Nine TFs, including CBF4, RAV1 and ERF104, were identified after cold treatment. Weighted gene co-expression network analysis (WGCNA) and gene ontology (GO) analysis revealed that these TFs may regulate cold response through different pathways. As a prime candidate TF, overexpression of VaRAV1 in grape cells improved its cold tolerance. The transgenic cells exhibited low electrolyte leakage and malondialdehyde content and high peroxidase activity. Moreover, the TF gene TCP8 and a gene involving in homogalacturonan biosynthesis were found to be regulated by VaRAV1, suggesting that the contribution of VaRAV1 to cold tolerance may be achieved by enhancing the stability of cell membrane and regulating the expression of target genes involved in plant cell wall composition. Our work provides novel insights into plant response to cold stress and demonstrates the utility of ATAC-seq and RNA-seq for the rapid identification of TFs in response to cold stress in grapevine. VaRAV1 may play an important role in adaption to cold stress.

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.

Transcriptional regulation of flavonoid biosynthesis in Artemisia annua by AaYABBY5

Artemisia annua is a medicinal plant rich in terpenes and flavonoids with useful biological activities such as antioxidant, anticancer, and antimalarial activities. The transcriptional regulation of flavonoid biosynthesis in A. annua has not been well-studied. In this study, we identified a YABBY family transcription factor, AaYABBY5, as a positive regulator of anthocyanin and total flavonoid contents in A. annua. AaYABBY5 was selected based on its similar expression pattern to the phenylalanine ammonia lyase (PAL), chalcone synthase (CHS), chalcone isomerase (CHI), and flavonol synthase (FLS) genes. A transient dual-luciferase assay in Nicotiana bethamiana with the AaYABBY5 effector showed a significant increase in the activity of the downstream LUC gene, with reporters AaPAL, AaCHS, AaCHI, and AaUFGT. The yeast one-hybrid system further confirmed the direct activation of these promoters by AaYABBY5. Gene expression analysis of stably transformed AaYABBY5 overexpression, AaYABBY5 antisense, and control plants revealed a significant increase in the expression of AaPAL, AaCHS, AaCHI, AaFLS, AaFSII, AaLDOX, and AaUFGT in AaYABBY5 overexpression plants. Moreover, their total flavonoid content and anthocyanin content were also found to increase. AaYABBY5 antisense plants showed a significant decrease in the expression of flavonoid biosynthetic genes, as well as a decrease in anthocyanin and total flavonoid contents. In addition, phenotypic analysis revealed deep purple-pigmented stems, an increase in the leaf lamina size, and higher trichome densities in AaYABBY5 overexpression plants. Together, these data proved that AaYABBY5 is a positive regulator of flavonoid biosynthesis in A. annua. Our study provides candidate transcription factors for the improvement of flavonoid concentrations in A. annua and can be further extended to elucidate its mechanism of regulating trichome development.

Identification and expression analysis of candidate genes related to seed dormancy and germination in the wheat GATA family

GATA transcription factors have been reported to function in plant growth and development and during various biotic/abiotic stresses in Arabidopsis and rice. However, the functions of wheat GATAs, particularly in the regulation of seed dormancy and germination, remain unclear. Here, we identified 78 TaGATAs in wheat and divided them into five subfamilies. Sixty-four paralogous pairs and 52 orthologous pairs were obtained, and Ka/Ks ratios showed that the TaGATAs had undergone strong purifying election during the evolutionary process. Triplet analysis indicated that a high homologue retention rate could explain the large number of TaGATAs in wheat. Gene structure analysis revealed that most members of the same subfamily had similar structures, and subcellular localization prediction indicated that most TaGATAs were located in the nucleus. Gene ontology annotation results showed that most TaGATAs had molecular functions in DNA and zinc binding, and promoter analysis suggested that they may play important roles in growth, development, and biotic/abiotic stress response. We combined three microarray datasets with qRT-PCR expression data from wheat varieties of contrasting dormancy and pre-harvest sprouting resistance levels during imbibition in order to identify ten candidate genes (TaGATA17/-25/-34/-37/-40/-46/-48/-51/-72/-73) that may be involved in the regulation of seed dormancy and germination in wheat. These findings provide valuable information for further dissection of TaGATA functions in the regulation of seed dormancy and germination, thereby enabling the improvement of wheat pre-harvest sprouting resistance by gene pyramiding.

Transcriptome-Based Identification and Functional Characterization of NAC Transcription Factors Responsive to Drought Stress in Capsicum annuum L

Capsicum annuum L. is one of the most cultivated Solanaceae species, and in the open field, water limitation leading to drought stress affects its fruit quality, fruit setting, fruit size and ultimately yield. We identified stage-specific and a common core set of differentially expressed genes, following RNA-seq transcriptome analyses of a breeding line subjected to acute drought stress followed by recovery (rewatering), at three stages of plant development. Among them, two NAC transcription factor (TF) genes, i.e., CaNAC072 and CaNAC104, were always upregulated after drought stress and downregulated after recovery. The two TF proteins were observed to be localized in the nucleus following their transient expression in Nicotiana benthamiana leaves. The expression of the two NACs was also induced by NaCl, polyethylene glycol (PEG) and abscisic acid (ABA) treatments, suggesting that CaNAC072 is an early, while CaNAC104 is a late abiotic stress-responsive gene. Virus-induced gene silencing (VIGS) of CaNAC104 did not affect the pepper plantlet’s tolerance to drought stress, while VIGS of CaNAC072 increased drought tolerance. Heterologous expression of CaNAC072 in Arabidopsis thaliana as well as in plants mutated for its homolog ANAC072 did not increase drought stress tolerance. This highlights a different role of the two NAC homologs in the two species. Here, we discuss the complex role of NACs as transcriptional switches in the response to drought stress in bell pepper.

DNA features beyond the transcription factor binding site specify target recognition by plant MYC2-related bHLH proteins

Transcription factors (TFs) regulate gene expression by binding to cis-regulatory sequences in the promoters of target genes. Recent research is helping to decipher in part the cis-regulatory code in eukaryotes, including plants, but it is not yet fully understood how paralogous TFs select their targets. Here we addressed this question by studying several proteins of the basic helix-loop-helix (bHLH) family of plant TFs, all of which recognize the same DNA motif. We focused on the MYC-related group of bHLHs, that redundantly regulate the jasmonate (JA) signaling pathway, and we observed a high correspondence between DNA-binding profiles in vitro and MYC function in vivo. We demonstrated that A/T-rich modules flanking the MYC-binding motif, conserved from bryophytes to higher plants, are essential for TF recognition. We observed particular DNA-shape features associated with A/T modules, indicating that the DNA shape may contribute to MYC DNA binding. We extended this analysis to 20 additional bHLHs and observed correspondence between in vitro binding and protein function, but it could not be attributed to A/T modules as in MYCs. We conclude that different bHLHs may have their own codes for DNA binding and specific selection of targets that, at least in the case of MYCs, depend on the TF-DNA interplay.

Underwater CAM photosynthesis elucidated by Isoetes genome

To conserve water in arid environments, numerous plant lineages have independently evolved Crassulacean Acid Metabolism (CAM). Interestingly, Isoetes, an aquatic lycophyte, can also perform CAM as an adaptation to low CO2 availability underwater. However, little is known about the evolution of CAM in aquatic plants and the lack of genomic data has hindered comparison between aquatic and terrestrial CAM. Here, we investigate underwater CAM in Isoetes taiwanensis by generating a high-quality genome assembly and RNA-seq time course. Despite broad similarities between CAM in Isoetes and terrestrial angiosperms, we identify several key differences. Notably, Isoetes may have recruited the lesser-known 'bacterial-type' PEPC, along with the 'plant-type' exclusively used in other CAM and C4 plants for carboxylation of PEP. Furthermore, we find that circadian control of key CAM pathway genes has diverged considerably in Isoetes relative to flowering plants. This suggests the existence of more evolutionary paths to CAM than previously recognized.

Genome-wide identification, phylogeny and expression analysis of AP2/ERF transcription factors family in sweet potato

Background

In recent years, much attention has been given to AP2/ERF transcription factors because they play indispensable roles in many biological processes, such as plant development and biotic and abiotic stress responses. Although AP2/ERFs have been thoroughly characterised in many plant species, the knowledge about this family in the sweet potato, which is a vital edible and medicinal crop, is still limited. In this study, a comprehensive genome-wide investigation was conducted to characterise the AP2/ERF gene family in the sweet potato.

Results

Here, 198 IbAP2/ERF transcription factors were obtained. Phylogenetic analysis classified the members of the IbAP2/ERF family into three groups, namely, ERF (172 members), AP2 (21 members) and RAV (5 members), which was consistent with the analysis of gene structure and conserved protein domains. The evolutionary characteristics of these IbAP2/ERF genes were systematically investigated by analysing chromosome location, conserved protein motifs and gene duplication events, indicating that the expansion of the IbAP2/ERF gene family may have been caused by tandem duplication. Furthermore, the analysis of cis-acting elements in IbAP2/ERF gene promoters implied that these genes may play crucial roles in plant growth, development and stress responses. Additionally, the available RNA-seq data and quantitative real-time PCR (qRT-PCR) were used to investigate the expression patterns of IbAP2/ERF genes during sweet potato root development as well as under multiple forms of abiotic stress, and we identified several developmental stage-specific and stress-responsive IbAP2/ERF genes. Furthermore, g59127 was differentially expressed under various stress conditions and was identified as a nuclear protein, which was in line with predicted subcellular localization results.

Conclusions

This study originally revealed the characteristics of the IbAP2/ERF superfamily and provides valuable resources for further evolutionary and functional investigations of IbAP2/ERF genes in the sweet potato.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08043-w.