Dof1 and Dof2 transcription factors are associated with expression of multiple genes involved in carbon metabolism in maize

Investigating the cis-Regulatory Basis of C3 and C4 Photosynthesis in Grasses at Single-Cell Resolution

While considerable knowledge exists about the enzymes pivotal for C4 photosynthesis, much less is known about the cis-regulation important for specifying their expression in distinct cell types. Here, we use single-cell-indexed ATAC-seq to identify cell-type-specific accessible chromatin regions (ACRs) associated with C4 enzymes for five different grass species. This study spans four C4 species, covering three distinct photosynthetic subtypes: Zea mays and Sorghum bicolor (NADP-ME), Panicum miliaceum (NAD-ME), Urochloa fusca (PEPCK), along with the C3 outgroup Oryza sativa. We studied the cis-regulatory landscape of enzymes essential across all C4 species and those unique to C4 subtypes, measuring cell-type-specific biases for C4 enzymes using chromatin accessibility data. Integrating these data with phylogenetics revealed diverse co-option of gene family members between species, showcasing the various paths of C4 evolution. Besides promoter proximal ACRs, we found that, on average, C4 genes have two to three distal cell-type-specific ACRs, highlighting the complexity and divergent nature of C4 evolution. Examining the evolutionary history of these cell-type-specific ACRs revealed a spectrum of conserved and novel ACRs, even among closely related species, indicating ongoing evolution of cis-regulation at these C4 loci. This study illuminates the dynamic and complex nature of CRE evolution in C4 photosynthesis, particularly highlighting the intricate cis-regulatory evolution of key loci. Our findings offer a valuable resource for future investigations, potentially aiding in the optimization of C3 crop performance under changing climatic conditions.

Transcription factor ZmDof22 enhances drought tolerance by regulating stomatal movement and antioxidant enzymes activities in maize (Zea mays L.)

KEY MESSAGE

The Dof22 gene encoding a deoxyribonucleic acid binding with one finger in maize, which is associated with its drought tolerance. The identification of drought stress regulatory genes is essential for the genetic improvement of maize yield. Deoxyribonucleic acid binding with one finger (Dof), a plant-specific transcription factor family, is involved in signal transduction, morphogenesis, and environmental stress responses. In present study, by weighted correlation network analysis (WGCNA) and gene co-expression network analysis, 15 putative Dof genes were identified from maize that respond to drought and rewatering. A real-time fluorescence quantitative PCR showed that these 15 genes were strongly induced by drought and ABA treatment, and among them ZmDof22 was highly induced by drought and ABA treatment. Its expression level increased by nearly 200 times after drought stress and more than 50 times after ABA treatment. After the normal conditions were restored, the expression levels were nearly 100 times and 40 times of those before treatment, respectively. The Gal4-LexA/UAS system and transcriptional activation analysis indicate that ZmDof22 is a transcriptional activator regulating drought tolerance and recovery ability in maize. Further, overexpressed transgenic and mutant plants of ZmDof22 by CRISPR/Cas9, indicates that the ZmDof22, improves maize drought tolerance by promoting stomatal closure, reduces water loss, and enhances antioxidant enzyme activity by participating in the ABA pathways. Taken together, our findings laid a foundation for further functional studies of the ZmDof gene family and provided insights into the role of the ZmDof22 regulatory network in controlling drought tolerance and recovery ability of maize.

The regulatory pathway of transcription factor MYB36 from Trichoderma asperellum Tas653 resistant to poplar leaf blight pathogen Alternaria alternata Aal004

In fungi, MYB transcription factors (TFs) mainly regulate growth, development, and resistance to stress. However, as major disease-resistance TFs, they have rarely been studied in biocontrol fungi. In this study, MYB36 of Trichoderma asperellum Tas653 (Ta) was shown to respond strongly to the stress caused by Alternaria alternata Aa1004. Compared with wild-type Ta (Ta-Wt), the inhibition rate of the MYB36 knockout strain (Ta-Kn) on Aa1004 decreased by 11.06%; the superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities decreased by 82.15 U/g, 0.19 OD470/min/g, and 1631.2 μmol/min/g, respectively. The MYB36 overexpression strain (Ta-Oe) not only enhanced hyperparasitism on Aa1004, caused its hyphae to swell, deform, or even rupture, but also reduced the incidence rate of poplar leaf blight. MYB36 regulates downstream (TFs, detoxification genes, defense genes, and other antifungal-related genes by binding to the cis-acting elements "ACAT" and "ATCG". Zinc finger TFs, as the main antifungal TFs, account for 90% of the total TFs, and Zn37.5 (23.24-) and Zn83.7 (23.18-fold) showed the greatest expression difference when regulated directly by MYB36. The detoxification genes mainly comprised 11 major major facilitator superfamily (MFS) genes, among which MYB36 directly increased the expression levels of three genes by more than 2-3.44-fold. The defense genes mainly encoded cytochrome P450 (P450) and hydrolases. e.g., P45061.3 (2-10.95-), P45060.2 (2-7.07-), and Hyd44.6 (2-2.30-fold). This study revealed the molecular mechanism of MYB36 regulation of the resistance of T. asperellum to A. alternata and provides theoretical guidance for the biocontrol of poplar leaf blight and the anti-disease mechanism of biocontrol fungi.

Genome-wide identification and characterization of flowering genes in Citrus sinensis (L.) Osbeck: a comparison among C. Medica L., C. Reticulata Blanco, C. Grandis (L.) Osbeck and C. Clementina

BACKGROUND

Flowering plays an important role in completing the reproductive cycle of plants and obtaining next generation of plants. In case of citrus, it may take more than a year to achieve progeny. Therefore, in order to fasten the breeding processes, the juvenility period needs to be reduced. The juvenility in plants is regulated by set of various flowering genes. The citrus fruit and leaves possess various medicinal properties and are subjected to intensive breeding programs to produce hybrids with improved quality traits. In order to break juvenility in Citrus, it is important to study the role of flowering genes. The present study involved identification of genes regulating flowering in Citrus sinensis L. Osbeck via homology based approach. The structural and functional characterization of these genes would help in targeting genome editing techniques to induce mutations in these genes for producing desirable results.

RESULTS

A total of 43 genes were identified which were located on all the 9 chromosomes of citrus. The in-silico analysis was performed to determine the genetic structure, conserved motifs, cis-regulatory elements (CREs) and phylogenetic relationship of the genes. A total of 10 CREs responsible for flowering were detected in 33 genes and 8 conserved motifs were identified in all the genes. The protein structure, protein-protein interaction network and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed to study the functioning of these genes which revealed the involvement of flowering proteins in circadian rhythm pathways. The gene ontology (GO) and gene function analysis was performed to functionally annotate the genes. The structure of the genes and proteins were also compared among other Citrus species to study the evolutionary relationship among them. The expression study revealed the expression of flowering genes in floral buds and ovaries. The qRT-PCR analysis revealed that the flowering genes were highly expressed in bud stage, fully grown flower and early stage of fruit development.

CONCLUSIONS

The findings suggested that the flowering genes were highly conserved in citrus species. The qRT-PCR analysis revealed the tissue specific expression of flowering genes (CsFT, CsCO, CsSOC, CsAP, CsSEP and CsLFY) which would help in easy detection and targeting of genes through various forward and reverse genetic approaches.

Evolution of plant cell-type-specific cis -regulatory elements

Cis -regulatory elements (CREs) are critical in regulating gene expression, and yet our understanding of CRE evolution remains a challenge. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa , integrating data from 104,029 nuclei representing 128 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species ( Zea mays, Sorghum bicolor, Panicum miliaceum , and Urochloa fusca ). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1 epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially critical silencer CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of cell-type-specific CRE evolution in plants.

Isolation and Functional Characterization of a Constitutive Promoter in Upland Cotton (Gossypium hirsutum L.)

Multiple cis-acting elements are present in promoter sequences that play critical regulatory roles in gene transcription and expression. In this study, we isolated the cotton FDH (Fiddlehead) gene promoter (pGhFDH) using a real-time reverse transcription-PCR (qRT-PCR) expression analysis and performed a cis-acting elements prediction analysis. The plant expression vector pGhFDH::GUS was constructed using the Gateway approach and was used for the genetic transformation of Arabidopsis and upland cotton plants to obtain transgenic lines. Histochemical staining and a β-glucuronidase (GUS) activity assay showed that the GUS protein was detected in the roots, stems, leaves, inflorescences, and pods of transgenic Arabidopsis thaliana lines. Notably, high GUS activity was observed in different tissues. In the transgenic lines, high GUS activity was detected in different tissues such as leaves, stalks, buds, petals, androecium, endosperm, and fibers, where the pGhFDH-driven GUS expression levels were 3-10-fold higher compared to those under the CaMV 35S promoter at 10-30 days post-anthesis (DPA) during fiber development. The results indicate that pGhFDH can be used as an endogenous constitutive promoter to drive the expression of target genes in various cotton tissues to facilitate functional genomic studies and accelerate cotton molecular breeding.

Compartmentation of photosynthesis gene expression in C4 maize depends on time of day

Compared with the ancestral C3 state, C4 photosynthesis occurs at higher rates with improved water and nitrogen use efficiencies. In both C3 and C4 plants, rates of photosynthesis increase with light intensity and are maximal around midday. We determined that in the absence of light or temperature fluctuations, photosynthesis in maize (Zea mays) peaks in the middle of the subjective photoperiod. To investigate the molecular processes associated with these temporal changes, we performed RNA sequencing of maize mesophyll and bundle sheath strands over a 24-h time course. Preferential expression of C4 cycle genes in these cell types was strongest between 6 and 10 h after dawn when rates of photosynthesis were highest. For the bundle sheath, DNA motif enrichment and gene coexpression analyses suggested members of the DNA binding with one finger (DOF) and MADS (MINICHROMOSOME MAINTENANCE FACTOR 1/AGAMOUS/DEFICIENS/Serum Response Factor)-domain transcription factor families mediate diurnal fluctuations in C4 gene expression, while trans-activation assays in planta confirmed their ability to activate promoter fragments from bundle sheath expressed genes. The work thus identifies transcriptional regulators and peaks in cell-specific C4 gene expression coincident with maximum rates of photosynthesis in the maize leaf at midday.

Identification and Functional Characterization of Two Major Loci Associated with Resistance against Brown Planthoppers (Nilaparvata lugens (Stål)) Derived from Oryza nivara

The brown planthopper (BPH) is a highly destructive pest of rice, causing significant economic losses in various regions of South and Southeast Asia. Researchers have made promising strides in developing resistance against BPH in rice. Introgression line RPBio4918-230S, derived from Oryza nivara, has shown consistent resistance to BPH at both the seedling and adult stages of rice plants. Segregation analysis has revealed that this resistance is governed by two recessive loci, known as bph39(t) and bph40(t), contributing to 21% and 22% of the phenotypic variance, respectively. We later mapped the genes using a backcross population derived from a cross between Swarna and RPBio4918-230S. We identified specific marker loci, namely RM8213, RM5953, and R4M17, on chromosome 4, flanking the bph39(t) and bph40(t) loci. Furthermore, quantitative expression analysis of candidate genes situated between the RM8213 and R4M17 markers was conducted. It was observed that eight genes exhibited up-regulation in RPBio4918-230S and down-regulation in Swarna after BPH infestation. One gene of particular interest, a serine/threonine-protein kinase receptor (STPKR), showed significant up-regulation in RPBio4918-230S. In-depth sequencing of the susceptible and resistant alleles of STPKR from Swarna and RPBio4918-230S, respectively, revealed numerous single nucleotide polymorphisms (SNPs) and insertion-deletion (InDel) mutations, both in the coding and regulatory regions of the gene. Notably, six of these mutations resulted in amino acid substitutions in the coding region of STPKR (R5K, I38L, S120N, T319A, T320S, and F348S) when compared to Swarna and the reference sequence of Nipponbare. Further validation of these mutations in a set of highly resistant and susceptible backcross inbred lines confirmed the candidacy of the STPKR gene with respect to BPH resistance controlled by bph39(t) and bph40(t). Functional markers specific for STPKR have been developed and validated and can be used for accelerated transfer of the resistant locus to elite rice cultivars.

Identification and promoter analysis of a GA-stimulated transcript 1 gene from Jatropha curcas

KEY MESSAGE

Overexpression of JcGAST1 promotes plant growth but inhibits pistil development. The pyrimidine box and CGTCA motif of the JcGAST1 promoter were responsible for the GA and MeJA responses. Members of the gibberellic acid-stimulated Arabidopsis (GASA) gene family play roles in plant growth and development, particularly in flower induction and seed development. However, there is still relatively limited knowledge of GASA genes in Jatropha curcas. Herein, we identified a GASA family gene from Jatropha curcas, namely, JcGAST1, which encodes a protein containing a conserved GASA domain. Sequence alignment showed that the JcGAST1 protein shares 76% sequence identity and 80% sequence similarity with SlGAST1. JcGAST1 had higher expression and protein levels in the female flowers than in the male flowers. Overexpression of JcGAST1 in tobacco promotes plant growth but inhibits pistil development. JcGAST1 expression was upregulated by GA and downregulated by MeJA. Promoter analysis indicated that the pyrimidine box and CGTCA motif were the GA- and MeJA-responsive elements of the JcGAST1 promoter. Using a Y1H screen, six transcription factors were found to interact with the pyrimidine box, and three transcription factors were found to interact with the CGTCA motif. Overall, the results of this study improve our understanding of the JcGAST1 gene and provide useful information for further studies.

Wheat DOF transcription factors TaSAD and WPBF regulate glutenin gene expression in cooperation with SPA

Grain storage proteins (GSPs) quantity and composition determine the end-use value of wheat flour. GSPs consists of low-molecular-weight glutenins (LMW-GS), high-molecular-weight glutenins (HMW-GS) and gliadins. GSP gene expression is controlled by a complex network of DNA-protein and protein-protein interactions, which coordinate the tissue-specific protein expression during grain development. The regulatory network has been most extensively studied in barley, particularly the two transcription factors (TFs) of the DNA binding with One Finger (DOF) family, barley Prolamin-box Binding Factor (BPBF) and Scutellum and Aleurone-expressed DOF (SAD). They activate hordein synthesis by binding to the Prolamin box, a motif in the hordein promoter. The BPBF ortholog previously identified in wheat, WPBF, has a transcriptional activity in expression of some GSP genes. Here, the wheat ortholog of SAD, named TaSAD, was identified. The binding of TaSAD to GSP gene promoter sequences in vitro and its transcriptional activity in vivo were investigated. In electrophoretic mobility shift assays, recombinant TaSAD and WPBF proteins bound to cis-motifs like those located on HMW-GS and LMW-GS gene promoters known to bind DOF TFs. We showed by transient expression assays in wheat endosperms that TaSAD and WPBF activate GSP gene expression. Moreover, co-bombardment of Storage Protein Activator (SPA) with WPBF or TaSAD had an additive effect on the expression of GSP genes, possibly through conserved cooperative protein-protein interactions.

Physiological and transcriptome analysis of Dendrobium officinale under low nitrogen stress

Nitrogen (N) is the main nutrient of plants, and low nitrogen usually affects plant growth and crop yield. The traditional Chinese herbal medicine Dendrobium officinale Kimura et. Migo is a typical low nitrogen-tolerant plant, and its mechanism in response to low nitrogen stress has not previously been reported. In this study, physiological measurements and RNA-Seq analysis were used to analyse the physiological changes and molecular responses of D. officinale under different nitrogen concentrations. The results showed that under low nitrogen levels, the growth, photosynthesis and superoxide dismutase activity were found to be significantly inhibited, while the activities of peroxidase and catalase, the content of polysaccharides and flavonoids significantly increased. Differentially expressed genes (DEGs) analysis showed that nitrogen and carbon metabolisms, transcriptional regulation, antioxidative stress, secondary metabolite synthesis and signal transduction all made a big difference in low nitrogen stress. Therefore, copious polysaccharide accumulation, efficient assimilation and recycling of nitrogen, as well as rich antioxidant components play critical roles. This study is helpful for understanding the response mechanism of D. officinale to low nitrogen levels, which might provide good guidance for practical production of high quality D. officinale .

Nitrogen-dependent binding of the transcription factor PBF1 contributes to the balance of protein and carbohydrate storage in maize endosperm

Fluctuations in nitrogen (N) availability influence protein and starch levels in maize (Zea mays) seeds, yet the underlying mechanism is not well understood. Here, we report that N limitation impacted the expression of many key genes in N and carbon (C) metabolism in the developing endosperm of maize. Notably, the promoter regions of those genes were enriched for P-box sequences, the binding motif of the transcription factor prolamin-box binding factor 1 (PBF1). Loss of PBF1 altered accumulation of starch and proteins in endosperm. Under different N conditions, PBF1 protein levels remained stable but PBF1 bound different sets of target genes, especially genes related to the biosynthesis and accumulation of N and C storage products. Upon N-starvation, the absence of PBF1 from the promoters of some zein genes coincided with their reduced expression, suggesting that PBF1 promotes zein accumulation in the endosperm. In addition, PBF1 repressed the expression of sugary1 (Su1) and starch branching enzyme 2b (Sbe2b) under normal N supply, suggesting that, under N-deficiency, PBF1 redirects the flow of C skeletons for zein toward the formation of C compounds. Overall, our study demonstrates that PBF1 modulates C and N metabolism during endosperm development in an N-dependent manner.

Allele-specific expression and chromatin accessibility contribute to heterosis in tea plants (Camellia sinensis)

Heterosis is extensively used to improve crop productivity, yet its allelic and chromatin regulation remains unclear. Based on our resolved genomes of the maternal TGY and paternal HD, we analyzed the contribution of allele-specific expression (ASE) and chromatin accessibility of JGY and HGY, the artificial hybrids of oolong tea with the largest cultivated area in China. The ASE genes (ASEGs) of tea hybrids with maternal-biased were mainly related to the energy and terpenoid metabolism pathways, whereas the ASEGs with paternal-biased tend to be enriched in glutathione metabolism, and these parental bias of hybrids may coordinate and lead to the acquisition of heterosis in more biological pathways. ATAC-seq results showed that hybrids have significantly higher accessible chromatin regions (ACRs) compared with their parents, which may confer broader and stronger transcriptional activity of genes in hybrids. The number of ACRs with significantly increased accessibility in hybrids was much greater than decreased, and the associated alleles were also affected by differential ACRs across different parents, suggesting enhanced positive chromatin regulation and potential genetic effects in hybrids. Core ASEGs of terpene and purine alkaloid metabolism pathways with significant positive heterosis have greater chromatin accessibility in hybrids, and were potentially regulated by several members of the MYB, DOF and TRB families. The binding motif of CsMYB85 in the promoter ACR of the rate-limiting enzyme CsDXS was verified by DAP-seq. These results suggest that higher numbers and more accessible ACRs in hybrids contribute to the regulation of ASEGs, thereby affecting the formation of heterotic metabolites.

Stress-induced higher vein density in the C3-C4 intermediate Moricandia suffruticosa under drought and heat stress

The evolution of C4 photosynthesis involved multiple anatomical and physiological modifications, yet our knowledge of the genetic regulation involved remains elusive. In this study, systematic analyses were conducted comparing the C3-C4 intermediate Moricandia suffruticosa and its C3 relative Brassica napus (rapeseed). We found that the leaves of M. suffruticosa had significantly higher vein density than those of B. napus, and the vein density was further increased in M. suffruticosa under drought and heat stress. Moreover, the bundle sheath distance, as the mean distance from the outer wall of one bundle sheath to the outer wall of an adjacent one, decreased and the number of centripetal chloroplasts in bundle sheath cells was found to be altered in M. suffruticosa leaves under drought and heat treatments. These results suggest that abiotic stress can induce a change in an intermediate C3-C4 anatomy towards a C4-like anatomy in land plants. By integrating drought and heat factors, co-expression network and comparative transcriptome analyses between M. suffruticosa and B. napus revealed that inducible auxin signaling regulated vascular development, and autophagy-related vesicle trafficking processes were associated with this stress-induced anatomical change. Overexpressing three candidate genes, MsERF02, MsSCL01, and MsDOF01, increased leaf vein density and/or enhanced photosynthetic assimilation and drought adaptability in the transgenic lines. The findings of this study may improve our understanding of the genetic regulation and evolution of C4 anatomy.

Evidence for Light and Tissue Specific Regulation of Genes Involved in Fructan Metabolism in Agave tequilana

Plant Glycoside Hydrolase Family 32 (PGHF32) contains the fructosyltransferases and fructan exohydrolase enzymes responsible for fructan metabolism, in addition to closely related vacuolar and cell wall acid invertases. Agave species produce complex and dynamic fructan molecules (agavins) requiring 4 different fructosyltransferase activities (1-SST, 1-FFT, 6G-FFT and 6-SFT) for their synthesis. Combined analysis of RNAseq and genome data for A. tequilana led to the characterization of the genes encoding 3 fructosyltransferases for this species and support the hypothesis that no separate 6-SFT type enzyme exists in A. tequilana, suggesting that at least one of the fructosyltransferases identified may have multiple enzymatic activities. Structures for PGHF32 genes varied for A. tequilana and between other plant species but were conserved for different enzyme types within a species. The observed patterns are consistent with the formation of distinct gene structures by intron loss. Promoter analysis of the PGHF32 genes identified abundant putative regulatory motifs for light regulation and tissue-specific expression, and these regulatory mechanisms were confirmed experimentally for leaf tissue. Motifs for phytohormone response, carbohydrate metabolism and dehydration responses were also uncovered. Based on the regulatory motifs, full-length cDNAs for MYB, GATA, DOF and GBF transcription factors were identified and their phylogenetic distribution determined by comparison with other plant species. In silico expression analysis for the selected transcription factors revealed both tissue-specific and developmental patterns of expression, allowing candidates to be identified for detailed analysis of the regulation of fructan metabolism in A. tequilana at the molecular level.

Selection and functional identification of Dof genes expressed in response to nitrogen in Populus simonii × Populus nigra

In plants, Dof transcription factors are involved in regulating the expression of a series of genes related to N uptake and utilization. Therefore, the present study investigated how DNA-binding with one finger (Dof) genes are expressed in response to nitrogen (N) form and concentration to clarify the role of Dof genes and their functions in promoting N assimilation and utilization in poplar. The basic characteristics and expression patterns of Dof genes in poplar were analyzed by the use of bioinformatics methods. Dof genes expressed in response to N were screened, after which the related genes were cloned and transformed into Arabidopsis thaliana; the physiological indexes and the expression of related genes were subsequently determined. The function of Dof genes was then verified in Arabidopsis thaliana plants grown in the presence of different N forms and concentrations. Forty-four Dof genes were identified, most of which were expressed in the roots and young leaves, and some of the Dof genes were expressed under ammonia- and nitrate-N treatments. Three genes related to N induction were cloned, their proteins were found to localize in the nucleus, and PnDof30 was successfully transformed into Arabidopsis thaliana for functional verification. On comparing Arabidopsis thaliana with WT Arabidopsis thaliana plants, Arabidopsis thaliana plants overexpressing the Dof gene grew better under low N levels; the contents of soluble proteins and chlorophyll significantly increased, while the soluble sugar content significantly decreased. The expressions of several AMT, NRT, and GS genes were upregulated, while the expressions of several others were downregulated, and the expression of PEPC and PK genes significantly increased. In addition, the activity of PEPC, PK, GS, and NR enzymes significantly increased. The results showed that overexpression of PnDof30 significantly increased the level of carbon and N metabolism and improved the growth of transgenic Arabidopsis thaliana plants under low-N conditions. The study revealed the biological significance of poplar Dof transcription factors in N response and regulation of related downstream gene expression and provided some meaningful clues to explain the huge difference between poplar and Arabidopsis thaliana transformed by exogenous Dof gene, which could promote the comprehensive understanding of the molecular mechanism of efficient N uptake and utilization in trees.

The Original Form of C4-Photosynthetic Phosphoenolpyruvate Carboxylase Is Retained in Pooids but Lost in Rice

Poaceae is the most prominent monocot family that contains the primary cereal crops wheat, rice, and maize. These cereal species exhibit physiological diversity, such as different photosynthetic systems and environmental stress tolerance. Phosphoenolpyruvate carboxylase (PEPC) in Poaceae is encoded by a small multigene family and plays a central role in C4-photosynthesis and dicarboxylic acid metabolism. Here, to better understand the molecular basis of the cereal species diversity, we analyzed the PEPC gene family in wheat together with other grass species. We could designate seven plant-type and one bacterial-type grass PEPC groups, ppc1a, ppc1b, ppc2a, ppc2b, ppc3, ppc4, ppcC4, and ppc-b, respectively, among which ppc1b is an uncharacterized type of PEPC. Evolutionary inference revealed that these PEPCs were derived from five types of ancient PEPCs (ppc1, ppc2, ppc3, ppc4, and ppc-b) in three chromosomal blocks of the ancestral Poaceae genome. C4-photosynthetic PEPC (ppcC4 ) had evolved from ppc1b, which seemed to be arisen by a chromosomal duplication event. We observed that ppc1b was lost in many Oryza species but preserved in Pooideae after natural selection. In silico analysis of cereal RNA-Seq data highlighted the preferential expression of ppc1b in upper ground organs, selective up-regulation of ppc1b under osmotic stress conditions, and nitrogen response of ppc1b. Characterization of wheat ppc1b showed high levels of gene expression in young leaves, transcriptional responses under nitrogen and abiotic stress, and the presence of a Dof1 binding site, similar to ppcC4 in maize. Our results indicate the evolving status of Poaceae PEPCs and suggest the functional association of ppc1-derivatives with adaptation to environmental changes.

Network and epigenetic characterization of subsets of genes specifically expressed in maize bundle sheath cells

Bundle sheath (BS) cells exhibit dramatically structural differences and functional variations at physiological, biochemical and epigenetic levels as compared to mesophyll (M) cells in maize. The regulatory mechanisms controlling functional divergences between M and BS have been extensively investigated. However, BS cell-related regulatory networks are still not completely characterized. To address this, we herein conducted WGCNA-related co-expression assays using bulk M and BS cell RNA-seq data sets and identified a module containing 384 genes highly expressed in BS cells (including 20 hub TFs) instead of M cells. According to the hub TF centered regulatory network, we found that Dof22 and Dof30 might act as key regulators in the regulation of expression of BS-specific genes, and several MYB TFs exhibited a high collaboration with Dof TFs. By comparing the enrichment levels of histone modifications, we found that genes in the aforementioned module were more enriched with histone acetylation as compared to other BS-enriched DEGs with similar expression levels. Moreover, we found that a subset of genes functioning in photosynthesis, protein auto processing and enzymatic activities were significantly enriched with broad H3K4me3. Thus, our study provides evidence showing that regulatory network and histone modifications may play vital roles in the regulation of a subset of genes with important functions in BS cells.

Wheat glutamine synthetase TaGSr-4B is a candidate gene for a QTL of thousand grain weight on chromosome 4B

Glutamine synthetase TaGSr-4B is a candidate gene for a QTL of thousand grain weight on 4B, and the gene marker is ready for wheat breeding. A QTL for thousand grain weight (TGW) in wheat was previously mapped on chromosome 4B in a DH population of Westonia × Kauz. For identifying the candidate genes of the QTL, wheat 90 K SNP array was used to saturate the existing linkage map, and four field trials plus one glasshouse experiment over five locations were conducted to refine the QTL. Three nitrogen levels were applied to two of those field trials, resulting in a TGW phenotype data set from nine environments. A robust TGW QTL cluster including 773 genes was detected in six environments with the highest LOD value of 13.4. Based on differentiate gene expression within the QTL cluster in an RNAseq data of Westonia and Kauz during grain filling, a glutamine synthesis gene (GS: TaGSr-4B) was selected as a potential candidate gene for the QTL. A SNP on the promoter region between Westonia and Kauz was used to develop a cleaved amplified polymorphic marker for TaGSr-4B gene mapping and QTL reanalysing. As results, TGW QTL appeared in seven environments, and in four out of seven environments, the TGW QTL were localized on the TaGSr-4B locus and showed significant contributions to the phenotype. Based on the marker, two allele groups of Westonia and Kauz formed showed significant differences on TGW in eight environments. In agreement with the roles of GS genes on nitrogen and carbon remobilizations, TaGSr-4B is likely the candidate gene of the TGW QTL on 4B and the TaGSr-4B gene marker is ready for wheat breeding.

Characterization of sucrose binding protein as a seed-specific promoter in transgenic tobacco Nicotiana tabacum L

Seed-specific expression using appropriate promoters is a recommended strategy for the efficiently producing valuable metabolites in transgenic plants. In the present study, we investigated the sequence of sucrose binding protein (SBP) as a seed-specific promoter to find the cis-acting elements specific to gene expression in seeds. The 1860 bp SBP sequence was analyzed using Plant Care and PLACE databases to find cis-acting elements, which resulted in a finding of 22 cis-acting elements required for seed expression. In addition, we have discovered cis- acting elements that are indirectly involved in triacylglycerol synthesis (GATABOX, DOFCOREZM, CACGTGMOTIF). The seed specificity of SBP was analyzed by generating a stable transgenic tobacco plant harboring β-glucuronidase (GUS) reporter gene under the control of the SBP promoter. Histochemical analysis of these transgenic tobacco plants indicated decreasing GUS activity in the leaves during the vegetative stage. However, the mature seeds of transgenic plants showed GUS activity. Moreover, the SBP promoter function in the seed oil content was evaluated by the expression of DGAT1. The expression analysis of DGAT1 in SBP-DGAT1 transgenic tobacco seeds using quantitative real-time PCR revealed a 7.8-fold increase in DGAT1 than in non-transgenic plants. Moreover, oil content increased up to 2.19 times more than in non-transgenic plants. And the oil content of the SBP-DGAT1 transgenic tobacco leaves did not change compared to the control plant. Therefore, we suggested that the SBP promoter could be used as a seed-specific promoter for targeted expression of desired genes in the metabolite engineering of oilseed crops.

Genome-wide analysis of AAAG and ACGT cis-elements in Arabidopsis thaliana reveals their involvement with genes downregulated under jasmonic acid response in an orientation independent manner

Cis-regulatory elements are regions of noncoding DNA that regulate the transcription of neighboring genes. The study of cis-element architecture that functions in transcription regulation are essential. AAAG and ACGT are a class of cis-regulatory elements, known to interact with Dof and bZIP transcription factors respectively, and are known to regulate the expression of auxin response, gibberellin response, floral development, light response, seed storage proteins genes, biotic and abiotic stress genes in plants. Analysis of the frequency of occurrence of AAAG and ACGT motifs from varying spacer lengths (0–30 base pair) between these 2 motifs in both possible orientations—AAAG _(N) ACGT and ACGT _(N) AAAG, in the promoters and genome of Arabidopsis thaliana which indicated preferred orientation of AAAG _(N) ACGT over ACGT _(N) AAAG across the genome and in promoters. Further, microarray analysis revealed the involvement of these motifs in the genes downregulated under jasmonic acid response in an orientation-independent manner. These results were further confirmed by the transient expression studies with promoter-reporter cassettes carrying AAAG and ACGT motifs in both orientations. Furthermore, cluster analysis on genes with AAAG _(N) ACGT and ACGT _(N) AAAG motifs orientations revealed clusters of genes to be involved in ABA signaling, transcriptional regulation, DNA binding, and metal ion binding. These findings can be utilized in designing synthetic promoters for the development of stress-tolerant transgenic plants and also provides an insight into the roles of these motifs in transcriptional regulation.

Single-Cell Transcriptome and Network Analyses Unveil Key Transcription Factors Regulating Mesophyll Cell Development in Maize

Background: Maize mesophyll (M) cells play important roles in various biological processes such as photosynthesis II and secondary metabolism. Functional differentiation occurs during M-cell development, but the underlying mechanisms for regulating M-cell development are largely unknown. Results: We conducted single-cell RNA sequencing (scRNA-seq) to profile transcripts in maize leaves. We then identified coregulated modules by analyzing the resulting pseudo-time-series data through gene regulatory network analyses. WRKY, ERF, NAC, MYB and Heat stress transcription factor (HSF) families were highly expressed in the early stage, whereas CONSTANS (CO)-like (COL) and ERF families were highly expressed in the late stage of M-cell development. Construction of regulatory networks revealed that these transcript factor (TF) families, especially HSF and COL, were the major players in the early and later stages of M-cell development, respectively. Integration of scRNA expression matrix with TF ChIP-seq and Hi-C further revealed regulatory interactions between these TFs and their targets. HSF1 and COL8 were primarily expressed in the leaf bases and tips, respectively, and their targets were validated with protoplast-based ChIP-qPCR, with the binding sites of HSF1 being experimentally confirmed. Conclusions: Our study provides evidence that several TF families, with the involvement of epigenetic regulation, play vital roles in the regulation of M-cell development in maize.

Chromatin and regulatory differentiation between bundle sheath and mesophyll cells in maize

C₄ plants partition photosynthesis enzymes between the bundle sheath (BS) and the mesophyll (M) cells for the better delivery of CO₂ to RuBisCO and to reduce photorespiration. To better understand how C₄ photosynthesis is regulated at the transcriptional level, we performed RNA-seq, ATAC-seq, ChIP-seq and Bisulfite-seq (BS-seq) on BS and M cells isolated from maize leaves. By integrating differentially expressed genes with chromatin features, we found that chromatin accessibility coordinates with epigenetic features, especially H3K27me3 modification and CHH methylation, to regulate cell type-preferentially enriched gene expression. Not only the chromatin-accessible regions (ACRs) proximal to the genes (pACRs) but also the distal ACRs (dACRs) are determinants of cell type-preferentially enriched expression. We further identified cell type-preferentially enriched motifs, e.g. AAAG for BS cells and TGACC/T for M cells, and determined their corresponding transcription factors: DOFs and WRKYs. The complex interaction between cis and trans factors in the preferential expression of C₄ genes was also observed. Interestingly, cell type-preferentially enriched gene expression can be fine-tuned by the coordination of multiple chromatin features. Such coordination may be critical in ensuring the cell type-specific function of key C₄ genes. Based on the observed cell type-preferentially enriched expression pattern and coordinated chromatin features, we predicted a set of functionally unknown genes, e.g. Zm00001d042050 and Zm00001d040659, to be potential key C₄ genes. Our findings provide deep insight into the architectures associated with C₄ gene expression and could serve as a valuable resource to further identify the regulatory mechanisms present in C₄ species.

Genome-wide analysis of Dof transcription factors and their response to cold stress in rice (Oryza sativa L.)

Background

Rice (Oryza sativa L.) is a food crop for humans worldwide. However, temperature has an effect during the vegetative and reproductive stages. In high-latitude regions where rice is cultivated, cold stress is a major cause of yield loss and plant death. Research has identified a group of plant-specific transcription factors, DNA binding with one zinc fingers (DOFs), with a diverse range of functions, including stress signaling and stress response during plant growth. The aim of this study was to identify Dof genes in two rice subspecies, indica and japonica, and screen for Dof genes that may be involved in cold tolerance during plant growth.

Results

A total of 30 rice Dofs (OsDofs) were identified using bioinformatics and genome-wide analyses and phylogenetically analyzed. The 30 OsDOFs were classified into six subfamilies, and 24 motifs were identified based on protein sequence alignment. The chromosome locations of OsDofs were determined and nine gene duplication events were identified. A joint phylogenetic analysis was performed on DOF protein sequences obtained from four monocotyledon species to examine the evolutionary relationship of DOF proteins. Expression profiling of OsDofs from two japonica cultivars (Longdao5, which is cold-tolerant, and Longjing11, which is cold-sensitive) revealed that OsDof1 and OsDof19 are cold-inducible genes. We examined the seed setting rates in OsDof1- and OsDof19-overexpression and RNAi lines and found that OsDof1 showed a response to cold stress.

Conclusions

Our investigation identified OsDof1 as a potential target for genetic breeding of rice with enhanced cold tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08104-0.

PsEND1 Is a Key Player in Pea Pollen Development Through the Modulation of Redox Homeostasis

Redox homeostasis has been linked to proper anther and pollen development. Accordingly, plant cells have developed several Reactive Oxygen Species (ROS)-scavenging mechanisms to maintain the redox balance. Hemopexins constitute one of these mechanisms preventing heme-associated oxidative stress in animals, fungi, and plants. Pisum sativum ENDOTHECIUM 1 (PsEND1) is a pea anther-specific gene that encodes a protein containing four hemopexin domains. We report the functional characterization of PsEND1 and the identification in its promoter region of cis-regulatory elements that are essential for the specific expression in anthers. PsEND1 promoter deletion analysis revealed that a putative CArG-like regulatory motif is necessary to confer promoter activity in developing anthers. Our data suggest that PsEND1 might be a hemopexin regulated by a MADS-box protein. PsEND1 gene silencing in pea, and its overexpression in heterologous systems, result in similar defects in the anthers consisting of precocious tapetum degradation and the impairment of pollen development. Such alterations were associated to the production of superoxide anion and altered activity of ROS-scavenging enzymes. Our findings demonstrate that PsEND1 is essential for pollen development by modulating ROS levels during the differentiation of the anther tissues surrounding the microsporocytes.

MP3RNA-seq: Massively parallel 3' end RNA sequencing for high-throughput gene expression profiling and genotyping

Transcriptome deep sequencing (RNA-seq) has become a routine method for global gene expression profiling. However, its application to large-scale experiments remains limited by cost and labor constraints. Here we describe a massively parallel 3' end RNA-seq (MP3RNA-seq) method that introduces unique sample barcodes during reverse transcription to permit sample pooling immediately following this initial step. MP3RNA-seq allows for handling of hundreds of samples in a single experiment, at a cost of about $6 per sample for library construction and sequencing. MP3RNA-seq is effective for not only high-throughput gene expression profiling, but also genotyping. To demonstrate its utility, we applied MP3RNA-seq to 477 double haploid lines of maize. We identified 19,429 genes expressed in at least 50% of the lines and 35,836 high-quality single nucleotide polymorphisms for genotyping analysis. Armed with these data, we performed expression and agronomic trait quantitative trait locus (QTL) mapping and identified 25,797 expression QTLs for 15,335 genes and 21 QTLs for plant height, ear height, and relative ear height. We conclude that MP3RNA-seq is highly reproducible, accurate, and sensitive for high-throughput gene expression profiling and genotyping, and should be generally applicable to most eukaryotic species.

Genome-Wide Analysis and Expression Profiles of the Dof Family in Cleistogenes songorica under Temperature, Salt and ABA Treatment

The DNA-binding with one zinc finger (Dof) family of plant-specific transcription factors has a variety of important functions in gene transcriptional regulation, development, and stress responses. However, the structure and expression patterns of Dof family have not been identified in Cleistogenes songorica, which is an important xerophytic and perennial gramineous grass in desert grassland. In this study, 50 Dof genes were identified in C. songorica and could be classified into four groups. According to genome-wide analysis, 46 of 50 Dof genes were located on 20 chromosomes, and the gene structure and conserved protein motif of these proteins were analyzed. In addition, phylogenetic analysis of Dof genes in C. songorica, Arabidopsis thaliana, Oryza sativa, and Brachypodium distachyon estimated the evolutionary relationships, and these genes were grouped into seven clusters. Moreover, the expression profiles of these Dof genes in C. songorica were analyzed in response to high/low temperature, salinity, and ABA treatments. These results will provide valuable information for future studies on gene classification, cloning, and functional characterization of this family in C. songorica.

Low Light/Darkness as Stressors of Multifactor-Induced Senescence in Rice Plants

Leaf senescence, as an integral part of the final development stage for plants, primarily remobilizes nutrients from the sources to the sinks in response to different stressors. The premature senescence of leaves is a critical challenge that causes significant economic losses in terms of crop yields. Although low light causes losses of up to 50% and affects rice yield and quality, its regulatory mechanisms remain poorly elucidated. Darkness-mediated premature leaf senescence is a well-studied stressor. It initiates the expression of senescence-associated genes (SAGs), which have been implicated in chlorophyll breakdown and degradation. The molecular and biochemical regulatory mechanisms of premature leaf senescence show significant levels of redundant biomass in complex pathways. Thus, clarifying the regulatory mechanisms of low-light/dark-induced senescence may be conducive to developing strategies for rice crop improvement. This review describes the recent molecular regulatory mechanisms associated with low-light response and dark-induced senescence (DIS), and their effects on plastid signaling and photosynthesis-mediated processes, chloroplast and protein degradation, as well as hormonal and transcriptional regulation in rice.

Genome-Wide In Silico Identification and Comparative Analysis of Dof Gene Family in Brassica napus

DNA binding with one finger (DOF) proteins are plant-specific transcription factors that play roles in diverse plant functions. However, little is known about the DOF protein repertoire of the allopolyploid crop, Brassica napus. This in silico study identified 117 Brassica napus Dof genes (BnaDofs) and classified them into nine groups (A, B1, B2, C1, C2.1, C2.2, C3, D1, and D2), based on phylogenetic analysis. Most members belonging to a particular group displayed conserved gene structural organisation and protein motif distribution. Evolutionary analysis exemplified that the divergence of the Brassica genus from Arabidopsis, the whole-genome triplication event, and the hybridisation of Brassica oleracea and Brassica rapa to form B. napus, followed by gene loss and rearrangements, led to the expansion and divergence of the Dof transcription factor (TF) gene family in B. napus. So far, this is the largest number of Dof genes reported in a single eudicot species. Functional annotation of BnaDof proteins, cis-element analysis of their promoters, and transcriptomic analysis suggested potential roles in organ development, the transition from the vegetative to the reproductive stage, light responsiveness, phytohormone responsiveness, as well as potential regulatory roles in abiotic stress. Overall, our results provide a comprehensive understanding of the molecular structure, evolution, and possible functional roles of Dof genes in plant development and abiotic stress response.

Cloning and functional characterization of seed-specific LEC1A promoter from peanut (Arachis hypogaea L.)

LEAFY COTYLEDON1 (LEC1) is a HAP3 subunit of CCAAT-binding transcription factor, which controls several aspects of embryo and postembryo development, including embryo morphogenesis, storage reserve accumulation and skotomorphogenesis. Herein, using the method of chromosomal walking, a 2707bp upstream sequence from the ATG initiation codon site of AhLEC1A which is a homolog of Arabidopsis LEC1 was isolated in peanut. Its transcriptional start site confirmed by 5’ RACE was located at 82 nt from 5’ upstream of ATG. The bioinformatics analysis revealed that there existed many tissue-specific elements and light responsive motifs in its promoter. To identify the functional region of the AhLEC1A promoter, seven plant expression vectors expressing the GUS (β-glucuronidase) gene, driven by 5’ terminal series deleted fragments of AhLEC1A promoter, were constructed and transformed into Arabidopsis. Results of GUS histochemical staining showed that the regulatory region containing 82bp of 5’ UTR and 2228bp promoter could facilitate GUS to express preferentially in the embryos at different development periods of Arabidopsis. Taken together, it was inferred that the expression of AhLEC1A during seed development of peanut might be controlled positively by several seed-specific regulatory elements, as well as negatively by some other regulatory elements inhibiting its expression in other organs. Moreover, the GUS expression pattern of transgenic seedlings in darkness and in light was relevant to the light-responsive elements scattered in AhLEC1A promoter segment, implying that these light-responsive elements harbored in the AhLEC1A promoter regulate skotomorphogenesis of peanut seeds, and AhLEC1A expression was inhibited after the germinated seedlings were transferred from darkness to light.

Recent Advances in Carbon and Nitrogen Metabolism in C3 Plants

C and N are the most important essential elements constituting organic compounds in plants. The shoots and roots depend on each other by exchanging C and N through the xylem and phloem transport systems. Complex mechanisms regulate C and N metabolism to optimize plant growth, agricultural crop production, and maintenance of the agroecosystem. In this paper, we cover the recent advances in understanding C and N metabolism, regulation, and transport in plants, as well as their underlying molecular mechanisms. Special emphasis is given to the mechanisms of starch metabolism in plastids and the changes in responses to environmental stress that were previously overlooked, since these changes provide an essential store of C that fuels plant metabolism and growth. We present general insights into the system biology approaches that have expanded our understanding of core biological questions related to C and N metabolism. Finally, this review synthesizes recent advances in our understanding of the trade-off concept that links C and N status to the plant's response to microorganisms.

Co-expression network analysis of protein phosphatase 2A (PP2A) genes with stress-responsive genes in Arabidopsis thaliana reveals 13 key regulators

Abiotic and biotic stresses adversely affect plant growth and development and eventually result in less yield and threaten food security worldwide. In plants, several studies have been carried out to understand molecular responses to abiotic and biotic stresses. However, the complete circuitry of stress-responsive genes that plants utilise in response to those environmental stresses are still unknown. The protein phosphatase 2A (PP2A) gene has been known to have a crucial role in abiotic and biotic stresses; but how it regulates the stress response in plants is still not known completely. In this study, we constructed gene co-expression networks of PP2A genes with stress-responsive gene datasets from cold, drought, heat, osmotic, genotoxic, salt, and wounding stresses to unveil their relationships with the PP2A under different conditions of stress. The graph analysis identified 13 hub genes and several influential genes based on closeness centrality score (CCS). Our findings also revealed the count of unique genes present in different settings of stresses and subunits. We also formed clusters of influential genes based on the stress, CCS, and co-expression value. Analysis of cis-regulatory elements (CREs), recurring in promoters of these genes was also performed. Our study has led to the identification of 16 conserved CREs.

Genome-Wide Identification and Expression Analysis of the Dof Gene Family in Medicago sativa L. Under Various Abiotic Stresses

The Dof transcription factor is a plant-specific transcriptional regulator that plays important roles in plant development and acts as a mediator in plant external stress responses. However, Dofs have previously been identified in several plants but not in alfalfa (Medicago sativa L.), one of the most widely cultivated forage legumes. In the present study, a total of 40 MsDof genes were identified, and the phylogenetic reconstruction, classification, conserved motifs, and expression patterns under abscisic acid (ABA), cold, heat, drought and salt stresses of these Dof genes were comprehensively analyzed. The Dof genes family in alfalfa could be classified into eight classes. Gene ontology (GO) and tissue-specific analysis indicated that most MsDof genes may be involved in biological functions during plant growth. Moreover, the expression profiles and quantitative real-time PCR analysis indicated that eight candidate abiotic tolerance genes were induced in response to four abiotic stresses. This study identified the possibility of abiotic tolerance candidate genes playing various roles in stress resistance at the whole genome level, which would provide new information on the Dof family in alfalfa.

Genome-wide identification of the maize 2OGD superfamily genes and their response to Fusarium verticillioides and Fusarium graminearum

In maize, eat rot and stalk rot caused by Fusarium verticillioides and Fusarium graminearum lead to contamination of moldy grains to produce mycotoxins. Identification of resistance genes against these pathogens for maize breeding is an effective way for disease control. Several 2-oxoglutarate-dependent dioxygenase (2OGD) proteins have been found to confer resistance to different pathogens in diverse plant species. However, little is known about the 2OGD superfamily in maize. Here, we identified 103 putative 2OGD genes in maize from a genome-wide analysis, and divided them into three classes - DOXA, DOXB, and DOXC. We further comprehensively investigated their gene structure, chromosome distribution, phylogenetic tree, gene-function enrichment, and expression profiles among different tissues. The genes encoding three 2OGD proteins, ACO, F3H, and NCS involved in ethylene biosynthesis, flavonoids biosynthesis, and alkaloids biosynthesis pathways, respectively, were identified to be induced by F. verticillioides and F. graminearum. The promoters of the three genes contain the binding sites for the transcription factor ZmDOF and ZmHSF, which are also induced by the two pathogens. The results imply that the three 2OGDs and the two transcription factors might be involved in the resistance to the two pathogens. This study provided a comprehensive understanding of the 2OGD superfamily in maize and laid the foundation for the further functional analysis of their roles in maize resistance to eat rot and stalk rot.

Genome-wide identification and evolution of Dof transcription factor family in cultivated and ancestral cotton species

The DNA-binding with one finger (Dof) proteins are transcription factors involved in many biological processes in plants. To predict the evolutionary pattern, a genome-wide in-silico analysis of Dof TFs family in diploid (Gossypium arboreum and Gossypium raimondii) and allotetraploid (Gossypium hirsutum and Gossypium barbadense) cotton species were carried out. In G. arboreum, we have identified 58 non-redundant genes encoding Dof proteins renamed as GaDof (G. arboreum Dof), 55 Dof genes were identified in G. raimondii (GrDof), 89 were predicted ffrom G. hirsutum (GhDof) and the highest, 110 Dof genes were identified in G. barbadense (GbDof). The phylogenetic analysis, physical location, gene structure, conserved domain analyses were also investigated for G. arboreum, G. raimondii, and G. hirsutum. The gene expression pattern in G. hirsutum, at different growth stages, revealing the probable involvement of some GhDof genes in growth and development. These genes may improve seed germination and growth in cotton.

CDF transcription factors: plant regulators to deal with extreme environmental conditions

In terrestrial environments, water and nutrient availabilities and temperature conditions are highly variable, and especially in extreme environments limit survival, growth, and reproduction of plants. To sustain growth and maintain cell integrity under unfavourable environmental conditions, plants have developed a variety of biochemical and physiological mechanisms, orchestrated by a large set of stress-responsive genes and a complex network of transcription factors. Recently, cycling DOF factors (CDFs), a group of plant-specific transcription factors (TFs), were identified as components of the transcriptional regulatory networks involved in the control of abiotic stress responses. The majority of the members of this TF family are activated in response to a wide range of adverse environmental conditions in different plant species. CDFs regulate different aspects of plant growth and development such as photoperiodic flowering-time control and root and shoot growth. While most of the functional characterization of CDFs has been reported in Arabidopsis, recent data suggest that their diverse roles extend to other plant species. In this review, we integrate information related to structure and functions of CDFs in plants, with special emphasis on their role in plant responses to adverse environmental conditions.

Heavy Metal Stress-Associated Proteins in Rice and Arabidopsis: Genome-Wide Identification, Phylogenetics, Duplication, and Expression Profiles Analysis

Heavy metal exposure is a serious environmental stress in plants. However, plants have evolved several strategies to improve their heavy metal tolerance. Heavy metal-associated proteins (HMPs) participate in heavy metal detoxification. Here, we identified 46 and 55 HMPs in rice and Arabidopsis, respectively, and named them OsHMP 1–46 and AtHMP 1–55 according to their chromosomal locations. The HMPs from both plants were divided into six clades based on the characteristics of their heavy metal-associated domains (HMA). The HMP gene structures and motifs varied greatly among the different classifications. The HMPs had high collinearity and were segmentally duplicated. A cis-element analysis revealed that the HMPs may be regulated by different transcription factors. An expression profile analysis disclosed that only eight OsHMPs were constitutive in rice tissues. Of these, the expression of OsHMP37 was far higher than that of the other seven genes while OsHMP28 was expressed exclusively in the roots. For Arabidopsis, nine AtHMPs presented with very high transcript levels in all organs. Most of the selected OsHMPs were differentially expressed in various tissues under different heavy metal stresses. Only OsHMP09, OsHMP18, and OsHMP22 showed higher expression levels in all tissues under different heavy metal stresses. In contrast, most of the selected AtHMPs had nearly constant expression levels in different tissues under various heavy metal stresses. The AtHMP20, AtHMP23, AtHMP25, AtHMP31, AtHMP35, AtHMP46 expression levels under different heavy metal stresses were higher in the leaves and roots. The foregoing discoveries elucidated HMP evolution in monocotyledonous and dicotyledonous plants and may helpful functionally characterize HMPs in the future.

Regulation and Evolution of C4 Photosynthesis

C₄ photosynthesis evolved multiple times independently from ancestral C₃ photosynthesis in a broad range of flowering land plant families and in both monocots and dicots. The evolution of C₄ photosynthesis entails the recruitment of enzyme activities that are not involved in photosynthetic carbon fixation in C₃ plants to photosynthesis. This requires a different regulation of gene expression as well as a different regulation of enzyme activities in comparison to the C₃ context. Further, C₄ photosynthesis relies on a distinct leaf anatomy that differs from that of C₃, requiring a differential regulation of leaf development in C₄. We summarize recent progress in the understanding of C₄-specific features in evolution and metabolic regulation in the context of C₄ photosynthesis.

Soybean Mutants Lacking Abundant Seed Storage Proteins Are Impaired in Mobilization of Storage Reserves and Germination

Spontaneous and radiation-induced mutants of soybean, despite loss of abundant seed proteins, have been reported to grow and reproduce normally without any apparent physiological abnormalities. Here, we report the development and characterization of a soybean line (BSH-2) that lacks several abundant seed storage proteins. One-dimensional and high-resolution two-dimensional gel electrophoresis revealed the absence of the α' and α subunits of β-conglycinin and G1, G2, G4, and G5 glycinin in the newly developed mutant line (BSH-2). Like our earlier developed soybean mutant line (BSH-3), the seeds of BSH-2 also accumulated high levels of free amino acids as compared with wild-type DN47 seeds. An examination of the germination rates revealed that both BSH-2 and BSH-3 had significantly lower germination rates compared with the parent line DN47. Two-dimensional gel electrophoresis analysis demonstrated that these mutants had slower rates of mobilization of seed storage proteins. The delayed mobilization of storage proteins in BSH-2 and BSH-3 seeds was also correlated with a delayed induction of proteolytic activity in the mutants when compared to DN47. Similarly, qRT-PCR analysis revealed distinct expression pattern of genes involved in proteolytic pathway in the mutants when compared to DN47. Transmission electron microscopy observation of soybean seeds at two germination stages revealed striking differences in the breakdown of protein storage vacuoles and lipid bodies in the mutants. Our study demonstrates that BSH-2 and BSH-3 are compromised in mobilization of storage reserves and the absence of abundant storage proteins may affect the seed germination efficiency and post-germinative seedling establishment.

Genome-wide analysis of wheat DNA-binding with one finger (Dof) transcription factor genes: evolutionary characteristics and diverse abiotic stress responses

Background

DNA binding with one finger (Dof) transcription factors play important roles in plant growth and abiotic stress responses. Although genome-wide identification and analysis of the DOF transcription factor family has been reported in other species, no relevant studies have emerged in wheat. The aim of this study was to investigate the evolutionary and functional characteristics associated with plant growth and abiotic stress responses by genome-wide analysis of the wheat Dof transcription factor gene family.

Results

Using the recently released wheat genome database (IWGSC RefSeq v1.0), we identified 96 wheat Dof gene family members, which were phylogenetically clustered into five distinct subfamilies. Gene duplication analysis revealed a broad and heterogeneous distribution of TaDofs on the chromosome groups 1 to 7, and obvious tandem duplication genes were present on chromosomes 2 and 3.Members of the same gene subfamily had similar exon-intron structures, while members of different subfamilies had obvious differences. Functional divergence analysis indicated that type-II functional divergence played a major role in the differentiation of the TaDof gene family. Positive selection analysis revealed that the Dof gene family experienced different degrees of positive selection pressure during the process of evolution, and five significant positive selection sites (30A, 31 T, 33A, 102G and 104S) were identified. Additionally, nine groups of coevolving amino acid sites, which may play a key role in maintaining the structural and functional stability of Dof proteins, were identified. The results from the RNA-seq data and qRT-PCR analysis revealed that TaDof genes exhibited obvious expression preference or specificity in different organs and developmental stages, as well as in diverse abiotic stress responses. Most TaDof genes were significantly upregulated by heat, PEG and heavy metal stresses.

Conclusions

The genome-wide analysis and identification of wheat DOF transcription factor family and the discovery of important amino acid sites are expected to provide new insights into the structure, evolution and function of the plant Dof gene family.

Allele specific expression of Dof genes responding to hormones and abiotic stresses in sugarcane

Dof transcription factors plant-specific and associates with growth and development in plants. We conducted comprehensive and systematic analyses of Dof transcription factors in sugarcane, and identified 29 SsDof transcription factors in sugarcane genome. Those SsDof genes were divided into five groups, with similar gene structures and conserved motifs within the same groups. Segmental duplications are predominant in the evolution of Dof in sugarcane. Cis-element analysis suggested that the functions of SsDofs were involved in growth and development, hormones and abiotic stresses responses in sugarcane. Expression patterns indicated that SsDof7, SsDof23 and SsDof24 had a comparatively high expression in all detected tissues, indicating these genes are crucial in sugarcane growth and development. Moreover, we examined the transcription levels of SsDofs under four plant hormone treatments, SsDof7-3 and SsDof7-4 were down-regulated after ABA treatment, while SsDof7-1 and SsDof7-2 were induced after the same treatment, indicating different alleles may play different roles in response to plant hormones. We also analyzed SsDofs' expression profiling under four abiotic stresses, SsDof5 and SsDof28 significantly responded to these four stresses, indicating they are associate with abiotic stresses responses. Collectively, our results yielded allele specific expression of Dof genes responding to hormones and abiotic stresses in sugarcane, and their cis-elements could be crucial for sugarcane improvement.

Mechanisms Regulating the Dynamics of Photosynthesis Under Abiotic Stresses

Photosynthesis sustains plant life on earth and is indispensable for plant growth and development. Factors such as unfavorable environmental conditions, stress regulatory networks, and plant biochemical processes limits the photosynthetic efficiency of plants and thereby threaten food security worldwide. Although numerous physiological approaches have been used to assess the performance of key photosynthetic components and their stress responses, though, these approaches are not extensive enough and do not favor strategic improvement of photosynthesis under abiotic stresses. The decline in photosynthetic capacity of plants due to these stresses is directly associated with reduction in yield. Therefore, a detailed information of the plant responses and better understanding of the photosynthetic machinery could help in developing new crop plants with higher yield even under stressed environments. Interestingly, cracking of signaling and metabolic pathways, identification of some key regulatory elements, characterization of potential genes, and phytohormone responses to abiotic factors have advanced our knowledge related to photosynthesis. However, our understanding of dynamic modulation of photosynthesis under dramatically fluctuating natural environments remains limited. Here, we provide a detailed overview of the research conducted on photosynthesis to date, and highlight the abiotic stress factors (heat, salinity, drought, high light, and heavy metal) that limit the performance of the photosynthetic machinery. Further, we reviewed the role of transcription factor genes and various enzymes involved in the process of photosynthesis under abiotic stresses. Finally, we discussed the recent progress in the field of biodegradable compounds, such as chitosan and humic acid, and the effect of melatonin (bio-stimulant) on photosynthetic activity. Based on our gathered researched data set, the logical concept of photosynthetic regulation under abiotic stresses along with improvement strategies will expand and surely accelerate the development of stress tolerance mechanisms, wider adaptability, higher survival rate, and yield potential of plant species.

Characterization of a vacuolar sucrose transporter, HbSUT5, from Hevea brasiliensis: involvement in latex production through regulation of intracellular sucrose transport in the bark and laticifers

BACKGROUND

Sucrose (Suc), as the precursor molecule for rubber biosynthesis in Hevea brasiliensis, is transported via phloem-mediated long-distance transport from leaves to laticifers in trunk bark, where latex (cytoplasm of laticifers) is tapped for rubber. In our previous report, six Suc transporter (SUT) genes have been cloned in Hevea tree, among which HbSUT3 is verified to play an active role in Suc loading to the laticifers. In this study, another latex-abundant SUT isoform, HbSUT5, with expressions only inferior to HbSUT3 was characterized especially for its roles in latex production.

RESULTS

Both phylogenetic analysis and subcellular localization identify HbSUT5 as a tonoplast-localized SUT protein under the SUT4-clade (=type III). Suc uptake assay in baker's yeast reveals HbSUT5 to be a typical Suc-H⁺ symporter, but its high affinity for Suc (Km = 2.03 mM at pH 5.5) and the similar efficiency in transporting both Suc and maltose making it a peculiar SUT under the SUT4-clade. At the transcript level, HbSUT5 is abundantly and preferentially expressed in Hevea barks. The transcripts of HbSUT5 are conspicuously decreased both in Hevea latex and bark by two yield-stimulating treatments of tapping and ethephon, the patterns of which are contrary to HbSUT3. Under the ethephon treatment, the Suc level in latex cytosol decreases significantly, but that in latex lutoids (polydispersed vacuoles) changes little, suggesting a role of the decreased HbSUT5 expression in Suc compartmentalization in the lutoids and thus enhancing the Suc sink strength in laticifers.

CONCLUSIONS

Our findings provide insights into the roles of a vacuolar sucrose transporter, HbSUT5, in Suc exchange between lutoids and cytosol in rubber-producing laticifers.

Genome-wide identification and transcriptional profiling of small heat shock protein gene family under diverse abiotic stress conditions in Sorghum bicolor (L.)

The small heat shock proteins (sHsps/Hsp20s) are the molecular chaperones that maintain proper folding, trafficking and disaggregation of proteins under diverse abiotic stress conditions. In the present investigation, a genome-wide scan revealed the presence of a total of 47 sHsps in Sorghum bicolor (SbsHsps), distributed across 10 subfamilies, the major subfamily being P (plastid) group with 17 genes. Chromosomes 1 and 3 appear as the hot spot regions for SbsHsps, and majority of them were found acidic, hydrophilic, unstable and intron less. Interestingly, promoter analysis indicated that they are associated with both biotic and abiotic stresses, as well as plant development. Sorghum sHsps exhibited 15 paralogous and 20 orthologous duplications. Expression analysis of 15 genes selected from different subfamilies showed high transcript levels in roots and leaves implying that they are likely to participate in the developmental processes. SbsHsp genes were highly induced by diverse abiotic stresses inferring their critical role in mediating the environmental stress responses. Gene expression data revealed that SbsHsp-02 is a candidate gene expressed in all the tissues under varied stress conditions tested. Our results contribute to the understanding of the complexity of SbsHsp genes and help to analyse them further for functional validation.

Identification and Functional Characterization of a Soybean (Glycine max) Thioesterase that Acts on Intermediates of Fatty Acid Biosynthesis

(1) Background: Plants possess many acyl-acyl carrier protein (acyl-ACP) thioesterases (TEs) with unique specificity. One such TE is methylketone synthase 2 (MKS2), an enzyme with a single-hotdog-fold structure found in several tomato species that hydrolyzes 3-ketoacyl-ACPs to give free 3-ketoacids. (2) Methods: In this study, we identified and characterized a tomato MKS2 homolog gene, namely, GmMKS2, in the genome of soybean (Glycine max). (3) Results: GmMKS2 underwent alternative splicing to produce three alternative transcripts, but only one encodes a protein with thioesterase activity when recombinantly expressed in Escherichia coli. Heterologous expression of the main transcript of GmMKS2, GmMKS2-X2, in E. coli generated various types of fatty acids, including 3-ketoacids-with 3-ketotetradecenoic acid (14:1) being the most abundant-cis-Δ5-dodecanoic acid, and 3-hydroxyacids, suggesting that GmMKS2 acts as an acyl-ACP thioesterase. In plants, the GmMKS2-X2 transcript level was found to be higher in the roots compared to other examined organs. In silico analysis revealed that there is a substantial enrichment of putative cis-regulatory elements related to disease-resistance responses and abiotic stress responses in the promoter of this gene. (4) Conclusions: GmMKS2 showed broad substrate specificities toward a wide range of acyl-ACPs that varied in terms of chain length, oxidation state, and saturation degree. Our results suggest that GmMKS2 might have a stress-related physiological function in G. max.

Rice DNA-Binding One Zinc Finger 24 (OsDOF24) Delays Leaf Senescence in a Jasmonate-Mediated Pathway

Leaf senescence is the final stage of leaf development and in cereal crops, the timing of senescence relative to grain filling has major effects on agronomic traits such as yield. Although many genetic factors are involved in the regulation of leaf senescence in cereals, the key regulators remain to be determined. Plant transcription factors with a conserved DOF (DNA-binding one zinc finger) domain play roles in multiple physiological processes. Here, we show a novel function for OsDOF24 as a repressor of leaf senescence in rice (Oryza sativa). In wild-type leaves, OsDOF24 expression rapidly decreased during natural senescence (NS) and dark-induced senescence (DIS). The gain-of-function mutant osdof24-D, which contains an enhancer-trap T-DNA in the OsDOF24 promoter, exhibited delayed leaf yellowing during NS and DIS. Transgenic plants overexpressing OsDOF24 showed the same phenotype during DIS. Reverse-transcription quantitative real-time PCR analysis revealed that senescence-associated genes (Osl85, Osl57 and OsNAP) and chlorophyll degradation genes (NYC1, NYC3 and SGR) were downregulated in the osdof24-D mutant during dark incubation. Among the phytohormones, only methyl jasmonate induced OsDOF24 expression. Furthermore, the reduced expression of jasmonate biosynthesis-related genes (OsLOX2, OsLOX8, OsHI-LOX, OsAOS1 and OsAOS2) in osdof24-D decreased endogenous jasmonate levels, resulting in delayed leaf senescence under DIS conditions. Yeast one-hybrid assays showed that OsDOF24 binds to the promoter region of OsAOS1. Taken together, our results demonstrate that OsDOF24 suppresses the induction of leaf senescence during vegetative growth by deactivating jasmonate biosynthetic pathways.

Exploring regulatory networks in plants: transcription factors of starch metabolism

Biological networks are complex (non-linear), redundant (cyclic) and compartmentalized at the subcellular level. Rational manipulation of plant metabolism may have failed due to inherent difficulties of a comprehensive understanding of regulatory loops. We first need to identify key factors controlling the regulatory loops of primary metabolism. The paradigms of plant networks are revised in order to highlight the differences between metabolic and transcriptional networks. Comparison between animal and plant transcription factors (TFs) reveal some important differences. Plant transcriptional networks function at a lower hierarchy compared to animal regulatory networks. Plant genomes contain more TFs than animal genomes, but plant proteins are smaller and have less domains as animal proteins which are often multifunctional. We briefly summarize mutant analysis and co-expression results pinpointing some TFs regulating starch enzymes in plants. Detailed information is provided about biochemical reactions, TFs and cis regulatory motifs involved in sucrose-starch metabolism, in both source and sink tissues. Examples about coordinated responses to hormones and environmental cues in different tissues and species are listed. Further advancements require combined data from single-cell transcriptomic and metabolomic approaches. Cell fractionation and subcellular inspection may provide valuable insights. We propose that shuffling of promoter elements might be a promising strategy to improve in the near future starch content, crop yield or food quality.

Banana (Musa acuminata) transcriptome profiling in response to rhizobacteria: Bacillus amyloliquefaciens Bs006 and Pseudomonas fluorescens Ps006

Background

Banana is one of the most important crops in tropical and sub-tropical regions. To meet the demands of international markets, banana plantations require high amounts of chemical fertilizers which translate into high farming costs and are hazardous to the environment when used excessively. Beneficial free-living soil bacteria that colonize the rhizosphere are known as plant growth-promoting rhizobacteria (PGPR). PGPR affect plant growth in direct or indirect ways and hold great promise for sustainable agriculture.

Results

PGPR of the genera Bacillus and Pseudomonas in banana cv. Williams were evaluated. These plants were produced through in vitro culture and inoculated individually with two rhizobacteria, Bacillus amyloliquefaciens strain Bs006 and Pseudomonas fluorescens strain Ps006. Control plants without microbial inoculum were also evaluated. These plants were kept in a controlled climate growth room with conditions required to favor plant-microorganism interactions. These interactions were evaluated at 1-, 48- and 96-h using transcriptome sequencing after inoculation to establish differentially expressed genes (DEGs) in plants elicited by the interaction with the two rhizobacteria. Additionally, droplet digital PCR was performed at 1, 48, 96 h, and also at 15 and 30 days to validate the expression patterns of selected DEGs. The banana cv. Williams transcriptome reported differential expression in a large number of genes of which 22 were experimentally validated. Genes validated experimentally correspond to growth promotion and regulation of specific functions (flowering, photosynthesis, glucose catabolism and root growth) as well as plant defense genes. This study focused on the analysis of 18 genes involved in growth promotion, defense and response to biotic or abiotic stress.

Conclusions

Differences in banana gene expression profiles in response to the rhizobacteria evaluated here (Bacillus amyloliquefaciens Bs006 and Pseudomonas fluorescens Ps006) are influenced by separate bacterial colonization processes and levels that stimulate distinct groups of genes at various points in time.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5763-5) contains supplementary material, which is available to authorized users.