Molecular interaction between PHO2 and GIGANTEA reveals a new crosstalk between flowering time and phosphate homeostasis in Oryza sativa

Functional divergence of FTL9 and FTL10 in flowering control in rice

BACKGROUND

Floral transition in cereals is a critical phenomenon influenced by exogenous and endogenous signals, determining crop yield and reproduction. Flowering Locus T-like (FT-like) genes encode a mobile florigen, the main signaling molecule for flowering.

RESULTS

In this study, we characterized two FT-like genes, FTL9 and FTL10, to study their functional diversity in flowering control in rice. We compared independent mutant lines of ftl10 with WT and observed negligible differences in the flowering phenotype, or agronomic traits implying potentially redundant roles of FTL10 loss-of-function in flowering control in rice. Nevertheless, we found that overexpression of FTL10, but not FTL9, substantially accelerated flowering, indicating the flowering-promoting role of FTL10 and the divergent functions between FTL9 and FTL10 in flowering. Besides flowering, additive agronomic roles were observed for FTL10-OE regulating the number of effective panicles per plant, the number of primary branches per panicle, and spikelets per panicle without regulating seed size. Mechanistically, our Y2H and BiFC analyses demonstrate that FTL10, in contrast to FTL9, can interact with FD1 and GF14c, forming a flowering activation complex and thereby regulating flowering.

CONCLUSION

Altogether, our results elucidate the regulatory roles of FTL9 and FTL10 in flowering control, unveiling the molecular basis of functional divergence between FTL10 and FTL9, which provides mechanistic insights into shaping the dynamics of flowering time regulation in rice.

Transcriptome analysis with different leaf blades identifies the phloem-specific phosphate transporter OsPHO1;3 required for phosphate homeostasis in rice

Phosphate (Pi) is essential for plant growth and development. One strategy to improve Pi use efficiency is to enhance Pi remobilization among leaves. Using transcriptome analysis with first (top) and fourth (down) leaf blades from rice (Oryza sativa) in Pi-sufficient and deficient conditions, we identified 1384 genes differentially expressed among these leaf blades. These genes were involved in physiological processes, metabolism, transport, and photosynthesis. Moreover, we identified the Pi efflux transporter gene, OsPHO1;3, responding to Pi-supplied conditions among these leaf blades. OsPHO1;3 is highly expressed in companion cells of phloem, but not xylem, in leaf blades and induced by Pi starvation. Mutation of OsPHO1;3 led to Pi accumulation in second to fourth leaves under Pi-sufficient conditions, but enhanced Pi levels in first leaves under Pi-deficient conditions. These Pi accumulations in leaves of Ospho1;3 mutants resulted from induction of OsPHT1;2 and OsPHT1;8 in root and reduction of Pi remobilization in leaf blades, revealed by the decreased Pi in phloem of leaves. Importantly, lack of OsPHO1;3 caused growth defects under a range of Pi-supplied conditions. These results demonstrate that Pi remobilization is essential for Pi homeostasis and plant growth irrespective of Pi-supplied conditions, and OsPHO1;3 plays an essential role in Pi remobilization for normal plant growth.

Transcriptomic Analysis Reveals the Response Mechanisms of Bell Pepper (Capsicum annuum) to Phosphorus Deficiency

Phosphorus (P) is an important nutritional element needed by plants. Roots obtain P as inorganic phosphate (Pi), mostly in H2PO-4 form. It is vital for plants to have a sufficient supply of Pi since it participates in important processes like photosynthesis, energy transfer, and protein activation, among others. The physicochemical properties and the organic material usually make Pi bioavailability in soil low, causing crops and undomesticated plants to experience variations in accessibility or even a persistent phosphate limitation. In this study, transcriptome data from pepper roots under low-Pi stress was analyzed in order to identify Pi starvation-responsive genes and their relationship with metabolic pathways and functions. Transcriptome data were obtained from pepper roots with Pi deficiency by RNASeq and analyzed with bioinformatic tools. A total of 97 differentially expressed genes (DEGs) were identified; Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment revealed that metabolic pathways, such as porphyrin and chlorophyll metabolism, were down-regulated, and galactose and fatty acid metabolism were up-regulated. The results indicate that bell pepper follows diverse processes related to low Pi tolerance regulation, such as the remobilization of internal Pi, alternative metabolic pathways to generate energy, and regulators of root development.

Rice domestication-associated transcription factor PROSTRATE GROWTH 1 controls plant and panicle architecture by regulating the expression of LAZY 1 and OsGIGANTEA, respectively

Plant architecture and panicle architecture are two critical agronomic traits that greatly affect the yield of rice (Oryza sativa). PROSTRATE GROWTH 1 (PROG1) encodes a key C2H2-type zinc-finger transcription factor and has pleiotropic effects on the regulation of both plant and panicle architecture, thereby influencing the grain yield. However, the molecular mechanisms through which PROG1 controls plant and panicle architecture remain unclear. In this study, we showed that PROG1 directly binds the LAZY 1 (LA1) promoter and acts as a repressor to inhibit LA1 expression. Conversely, LA1 acts as a repressor of PROG1 by directly binding to the PROG1 promoter. These two genes play antagonistic roles in shaping plant architecture by regulating both tiller angle and tiller number. Interestingly, our data showed that PROG1 controls panicle architecture through direct binding to the intragenic regulatory regions of OsGIGANTEA (OsGI) and subsequent activation of its expression. Collectively, we have identified two crucial targets of PROG1, LA1 and OsGI, shedding light on the molecular mechanisms underlying plant and panicle architecture control by PROG1. Our study provides valuable insights into the regulation of key domestication-related traits in rice and identifies potential targets for future high-yield rice breeding.

INDETERMINATE1 autonomously regulates phosphate homeostasis upstream of the miR399-ZmPHO2 signaling module in maize

The macronutrient phosphorus is essential for plant growth and development. Plants have evolved multiple strategies to increase the efficiency of phosphate (Pi) acquisition to protect themselves from Pi starvation. However, the crosstalk between Pi homeostasis and plant development remains to be explored. Here, we report that overexpressing microRNA399 (miR399) in maize (Zea mays) is associated with premature senescence after pollination. Knockout of ZmPHO2 (Phosphate 2), a miR399 target, resulted in a similar premature senescence phenotype. Strikingly, we discovered that INDETERMINATE1 (ID1), a floral transition regulator, inhibits the transcription of ZmMIR399 genes by directly binding to their promoters, alleviating the repression of ZmPHO2 by miR399 and ultimately contributing to the maintenance of Pi homeostasis in maize. Unlike ZmMIR399 genes, whose expression is induced by Pi deficiency, ID1 expression was independent of the external inorganic orthophosphate status, indicating that ID1 is an autonomous regulator of Pi homeostasis. Furthermore, we show that ZmPHO2 was under selection during maize domestication and cultivation, resulting in a more sensitive response to Pi starvation in temperate maize than in tropical maize. Our study reveals a direct functional link between Pi-deprivation sensing by the miR399-ZmPHO2 regulatory module and plant developmental regulation by ID1.

Genome-wide identification and characterization of mungbean CIRCADIAN CLOCK ASSOCIATED 1 like genes reveals an important role of VrCCA1L26 in flowering time regulation

Background

CIRCADIAN CLOCK ASSOCIATED 1 like (CCA1L) proteins are important components that participate in plant growth and development, and now have been characterized in multiple plant species. However, information on mungbean CCA1L genes is limited.

Results

In this study, we identified 27 VrCCA1L genes from the mungbean genome. VrCCA1L genes were unevenly distributed on 10 of the 11 chromosomes and showed one tandem and two interchromosomal duplication events. Two distinct kinds of conserved MYB domains, MYB 1 and MYB 2, were found, and the conserved SHAQK(Y/F) F sequence was found at the C terminus of each MYB 2 domain. The VrCCA1Ls displayed a variety of exon-intron organizations, and 24 distinct motifs were found among these genes. Based on phylogenetic analysis, VrCCA1L proteins were classified into five groups; group I contained the most members, with 11 VrCCA1Ls. VrCCA1L promoters contained different types and numbers of cis-acting elements, and VrCCA1Ls showed different expression levels in different tissues. The VrCCA1Ls also displayed distinct expression patterns under different photoperiod conditions throughout the day in leaves. VrCCA1L26 shared greatest homology to Arabidopsis CCA1 and LATE ELONGATED HYPOCOTYL (LHY). It delayed the flowering time in Arabidopsis by affecting the expression levels of CONSTANS (CO), FLOWERING LOCUS T (FT), and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1).

Conclusion

We identified and characterized 27 VrCCA1L genes from mungbean genome, and investigated their spatio-temporal expression patterns. Further analysis revealed that VrCCA1L26 delayed flowering time in transgenic Arabidopsis plants. Our results provide useful information for further functional characterization of the VrCCA1L genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08620-7.

Histone ubiquitination controls organ size in cotton (Gossypium hirsutum)

Ubiquitination plays a vital role in modifying protein activity and destiny. Ub-conjugating enzyme E2 is one of the enzymes that participates in this precise process. There are at least 169 E2 proteins in the allotetraploid cotton (Gossypium hirsutum), but their function remains unknown. Here we identify an E2 gene GhUBC2L and show its positive role in cell proliferation and expansion. Complete knock-down of GhUBC2L in cotton resulted in retarded growth and reduced organ size. Conversely, overexpression of GhUBC2L promoted cotton growth, generating enlarged organs in size. Monoubiquitination of H2A and H2B was strongly impaired in GhUBC2L-suppressed cotton but slightly enhanced in GhUBC2L-overexpressed plant. GhUbox8, a U-box type E3 ligase protein, was found to interact with GhUBC2L both in vivo and in vitro, indicating their synergistical function in protein ubiquitination. Furthermore, GhUbox8 was shown to interact with a series of histone proteins, including histone H2A and H2B, indicating its potential monoubiquitination on H2A and H2B. Expression of genes relating to cell cycle and organ development were altered when the expression of GhUBC2L was changed. Our results show that GhUBC2L modulates histone monoubiquitination synergistically with GhUbox8 to regulate the expression of genes involved in organ development and cell cycle, thus controlling organ size in cotton. This research provides new insights into the role of protein ubiquitination in organ size control. Histone monoubiquitination plays an important role in plant development. Here, we identified an E2 enzyme GhUBC2L that modulates histone monoubiquitination synergistically with an E3 ligase GhUbox8 to mediate organ size control in cotton.

Heterosis Derived From Nonadditive Effects of the BnFLC Homologs Coordinates Early Flowering and High Yield in Rapeseed (Brassica napus L.)

Early flowering facilitates crops to adapt multiple cropping systems or growing regions with a short frost-free season; however, it usually brings an obvious yield loss. In this study, we identified that the three genes, namely, BnFLC.A2, BnFLC.C2, and BnFLC.A3b, are the major determinants for the flowering time (FT) variation of two elite rapeseed (Brassica napus L.) accessions, i.e., 616A and R11. The early-flowering alleles (i.e., Bnflc.a2 and Bnflc.c2) and late-flowering allele (i.e., BnFLC.A3b) from R11 were introgressed into the recipient parent 616A through a breeding strategy of marker-assisted backcross, giving rise to eight homozygous near-isogenic lines (NILs) associated with these three loci and 19 NIL hybrids produced by the mutual crossing of these NILs. Phenotypic investigations showed that NILs displayed significant variations in both FT and plant yield (PY). Notably, genetic analysis indicated that BnFLC.A2, BnFLC.C2, and BnFLC.A3b have additive effects of 1.446, 1.365, and 1.361 g on PY, respectively, while their dominant effects reached 3.504, 2.991, and 3.284 g, respectively, indicating that the yield loss caused by early flowering can be successfully compensated by exploring the heterosis of FT genes in the hybrid NILs. Moreover, we further validated that the heterosis of FT genes in PY was also effective in non-NIL hybrids. The results demonstrate that the exploration of the potential heterosis underlying the FT genes can coordinate early flowering (maturation) and high yield in rapeseed (B. napus L.), providing an effective strategy for early flowering breeding in crops.

Insights Into MicroRNA-Mediated Regulation of Flowering Time in Cotton Through Small RNA Sequencing

The timing of flowering is a key determinant for plant reproductive. It has been demonstrated that microRNAs (miRNAs) play an important role in transition from the vegetative to reproductive stage in cotton; however, knowledge remains limited about the regulatory role of miRNAs involved in flowering time regulation in cotton. To elucidate the molecular basis of miRNAs in response to flowering time in cotton, we performed high-throughput small RNA sequencing at the fifth true leaf stage. We identified 56 and 43 miRNAs that were significantly up- and downregulated in two elite early flowering cultivars (EFC) compared with two late flowering cultivars (LFC), respectively. The miRNA targets by RNA sequencing analysis showed that GhSPL4 in SBP transcription factor family targeted by GhmiR156 was significantly upregulated in EFCs. Co-expression regulatory network analysis (WGCNA) revealed that GhSOC1, GhAP1, GhFD, GhCOL3, and GhAGL16 act as node genes in the auxin- and gibberellin-mediated flowering time regulatory networks in cotton. Therefore, elucidation of miRNA-mediated flowering time regulatory network will contribute to our understanding of molecular mechanisms underlying flowering time in cotton.

tae-miR399-UBC24 Module Enhances Freezing Tolerance in Winter Wheat via a CBF Signaling Pathway

Although the regulation of Pi homeostasis by miR399 has been studied in various plant species, its underlying molecular mechanism in response to freezing stress is still poorly understood. In this work, we found that the expression of tae-miR399 and its target gene TaUBC24 in the tillering nodes of the strong cold-resistant winter wheat cultivar Dongnongdongmai1 (Dn1) was not only significantly altered after severe winters but also responsive to short-term freezing stress. TaUBC24 physically interacted with TaICE1. Enhanced freezing tolerance was observed for tae-miR399-overexpressing Arabidopsis lines. Under freezing stress, overexpression of tae-miR399 ultimately decreased the expression of AtUBC24, inhibiting the degradation of AtICE1, which increased the expression of genes involved in the CBF signaling pathway and starch metabolism and promoted the activities of antioxidant enzymes. These results will improve our understanding of the molecular mechanism through which the miR399-UBC24 module plays a cardinal role in regulating plant freezing stress tolerance through mediation of downstream pathways.

Characterization and Comparative Analysis of RWP-RK Proteins from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea

RWP-RK proteins are important factors involved in nitrate response and gametophyte development in plants, and the functions of RWP-RK proteins have been analyzed in many species. However, the characterization of peanut RWP-RK proteins is limited. In this study, we identified 16, 19, and 32 RWP-RK members from Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively, and investigated their evolution relationships. The RWP-RK proteins were classified into two groups, RWP-RK domain proteins and NODULE-INCEPTION-like proteins. Chromosomal distributions, gene structures, and conserved motifs of RWP-RK genes were compared among wild and cultivated peanuts. In addition, we identified 12 orthologous gene pairs from the two wild peanut species, 13 from A. duranensis and A. hypogaea, and 13 from A. ipaensis and A. hypogaea. One, one, and seventeen duplicated gene pairs were identified within the A. duranensis, A. ipaensis, and A. hypogaea genomes, respectively. Moreover, different numbers of cis-acting elements in the RWP-RK promoters were found in wild and cultivated species (87 in A. duranensis, 89 in A. ipaensis, and 92 in A. hypogaea), and as a result, many RWP-RK genes showed distinct expression patterns in different tissues. Our study will provide useful information for further functional and evolutionary analysis of the RWP-RK genes.

Exogenous Cry1Ab/c Protein Recruits Different Endogenous Proteins for Its Function in Plant Growth and Development

Current risk assessments of transgenic crops do not take into consideration whether exogenous proteins interact with endogenous proteins and thereby induce unintended effects in the crops. Therefore, the unintended effects through protein interactions in insect-resistant transgenic rice merit investigation. Here, a yeast two-hybrid assay was used to evaluate interactions between Bacillus thuringiensis (Bt) protein-derived Cry1Ab/c insect resistance rice Huahui-1 and the endogenous proteins of its parental rice Minghui-63. The authenticity of the strongest interactions of Cry1Ab/c and 14 endogenous proteins involved in photosynthesis and stress resistance, which may be primarily responsible for the significant phenotypic differences between transgenic Huahui-1 and parental Minghui-63, were then analyzed and validated by subcellular co-localization, bimolecular fluorescence complementation and co-immunoprecipitation. As the exogenous full-length Cry1Ab/c protein was found to have self-activating activity, we cleaved it - into three segments based on its three domains, and these were screened for interaction with host proteins using the yeast two-hybrid assay. Sixty endogenous proteins related to the regulation of photosynthesis, stress tolerance, and substance metabolism were found to interact with the Cry1Ab/c protein. The results of bimolecular fluorescence complementation and co-immunoprecipitation verified the interactions between the full-length Cry1Ab/c protein and 12 endogenous proteins involved in photosynthesis 23KD, G, PSBP, Rubisco, Trx, THF1 and stress resistance CAMTAs, DAHP, E3s, HKMTs, KIN13A, FREE1. We used a combination of yeast two-hybrid, bimolecular fluorescence complementation, and co-immunoprecipitation to identify Cry1Ab/c interacting with rice proteins that seem to be associated with the observed unintended effects on photosynthesis and stress resistance between Huahui-1 and Minghui-63 rice plants, and analyze the possible interaction mechanisms by comparing differences in cell localization and interaction sites between these interactions. The results herein provide a molecular analytical system to qualify and quantify the interactions between exogenous proteins and the endogenous proteins of the recipient crop. It could help elucidate both the positive and negative effects of creating transgenic plants and predict their potential risks as well as net crop quality and yield.

B-box Proteins in Arachis duranensis: Genome-Wide Characterization and Expression Profiles Analysis

B-box (BBX) proteins are important factors involved in plant growth and developmental regulation, and they have been identified in many species. However, information on the characteristics and transcription patterns of BBX genes in wild peanut are limited. In this study, we identified and characterized 24 BBX genes from a wild peanut, Arachis duranensis. Many characteristics were analyzed, including chromosomal locations, phylogenetic relationships, and gene structures. Arachis duranensis B-box (AdBBX) proteins were grouped into five classes based on the diversity of their conserved domains: I (3 genes), II (4 genes), III (4 genes), IV (9 genes), and V (4 genes). Fifteen distinct motifs were found in the 24 AdBBX proteins. Duplication analysis revealed the presence of two interchromosomal duplicated gene pairs, from group II and IV. In addition, 95 kinds of cis-acting elements were found in the genes’ promoter regions, 53 of which received putative functional predictions. The numbers and types of cis-acting elements varied among different AdBBX promoters, and, as a result, AdBBX genes exhibited distinct expression patterns in different tissues. Transcriptional profiling combined with synteny analysis suggests that AdBBX8 may be a key factor involved in flowering time regulation. Our study will provide essential information for further functional investigation of AdBBX genes.

Molecular and transcriptional characterization of phosphatidyl ethanolamine-binding proteins in wild peanuts Arachis duranensis and Arachis ipaensis

BACKGROUND

Phosphatidyl ethanolamine-binding proteins (PEBPs) are involved in the regulation of plant architecture and flowering time. The functions of PEBP genes have been studied in many plant species. However, little is known about the characteristics and expression profiles of PEBP genes in wild peanut species, Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanuts.

RESULTS

In this study, genome-wide identification methods were used to identify and characterize a total of 32 peanut PEBP genes, 16 from each of the two wild peanut species, A. duranensis and A. ipaensis. These PEBP genes were classified into 3 groups (TERMINAL FLOWER1-like, FLOWERING LOCUS T-like, and MOTHER OF FT AND TFL1-like) based on their phylogenetic relationships. The gene structures, motifs, and chromosomal locations for each of these PEBPs were analyzed. In addition, 4 interchromosomal duplications and 1 tandem duplication were identified in A. duranensis, and 2 interchromosomal paralogs and 1 tandem paralog were identified in A. ipaensis. Ninety-five different cis-acting elements were identified in the PEBP gene promoter regions and most genes had different numbers and types of cis-elements. As a result, the transcription patterns of these PEBP genes varied in different tissues and under long day and short day conditions during different growth phases, indicating the functional diversities of PEBPs in different tissues and their potential functions in plant photoperiod dependent developmental pathways. Moreover, our analysis revealed that AraduF950M/AraduWY2NX in A. duranensis, and Araip344D4/Araip4V81G in A. ipaensis are good candidates for regulating plant architecture, and that Aradu80YRY, AraduYY72S, and AraduEHZ9Y in A. duranensis and AraipVEP8T in A. ipaensis may be key factors regulating flowering time.

CONCLUSION

Sixteen PEBP genes were identified and characterized from each of the two diploid wild peanut genomes, A. duranensis and A. ipaensis. Genetic characterization and spatio-temporal expression analysis support their importance in plant growth and development. These findings further our understanding of PEBP gene functions in plant species.

Regulation of Ubiquitination Is Central to the Phosphate Starvation Response

As sessile organisms, plants have developed numerous strategies to overcome the limiting availability of the essential nutrient phosphate in nature. Recent studies reveal that post-translational modification (PTM) by ubiquitination is an important and central regulation mechanism in the plant phosphate starvation response (PSR). Ubiquitination precisely modulates the stability and trafficking of proteins in response to the heterogeneous phosphate supplement. Induction of autophagy provides novel insights into the molecular mechanisms under phosphate starvation. In this review, we present and discuss novel findings on the regulation of diverse PSRs through ubiquitination. Resolving these regulation mechanisms will pave the way to improve phosphate acquisition and utilization efficiency in crops.

Conserved MicroRNA Act Boldly During Sprout Development and Quality Formation in Pingyang Tezaocha (Camellia sinensis)

Tea tree [Camellia sinensis (L.) O. Kuntze] is an important leaf (sometimes tender stem)-using commercial plant with many medicinal uses. The development of newly sprouts would directly affect the yield and quality of tea product, especially significant for Pingyang Tezaocha (PYTZ) which takes up a large percent in the early spring tea market. MicroRNA (miRNA), particularly the conserved miRNAs, often position in the center of subtle and complex gene regulatory systems, precisely control the biological processes together with other factors in a spatio-temporal pattern. Here, quality-determined metabolites catechins, theanine and caffeine in PYTZ sprouts including buds (sBud), different development stages of leaves (sL1, sL2) and stems (sS1, sS2) were quantified. A total of 15 miRNA libraries of the same tissue with three repetitions for each were constructed to explore vital miRNAs during the biological processes of development and quality formation. We analyzed the whole miRNA profiles during the sprout development and defined conserved miRNA families in the tea plant. The differentially expressed miRNAs related to the expression profiles buds, leaves, and stems development stages were described. Twenty one miRNAs and eight miRNA-TF pairs that most likely to participate in regulating development, and at least two miRNA-TF-metabolite triplets that participate in both development and quality formation had been filtered. Our results indicated that conserved miRNA act boldly during important biological processes, they are (i) more likely to be linked with morphological function in primary metabolism during sprout development, and (ii) hold an important position in secondary metabolism during quality formation in tea plant, also (iii) coordinate with transcription factors in forming networks of complex multicellular organism regulation.

The FvPHR1 transcription factor control phosphate homeostasis by transcriptionally regulating miR399a in woodland strawberry

Plants have evolved phosphate (Pi) starvation response to adapt the low-Pi environment. The regulation of adaptive responses to phosphorus deficiency by the PHR1-miR399-PHO2 module has been well studied in Arabidopsis thaliana but not in strawberry. Transcription factor PHR1 as the central regulator in the Pi starvation signaling has been revealed in a few plant species. However, the function of PHR1 homologues in strawberry is still unknown. In this study, a total of 13 MYB-CC genes were identified in the woodland strawberry (Fragaria vesca) genome and the FvPHR1 gene was characterized. FvPHR1 contains MYB domain and coiled-coil (CC) domain and is localized in the nucleus. FvPHR1 also exhibits trans-activation ability. Furthermore, the P content in leaves of FvPHR1-overexpressing woodland strawberries was significantly increased by 1.38-fold to 1.78-fold compared with that in the wild type. FvPHR1 was also demonstrated to directly bind to the FvMIR399a promoter and positively regulate the expression of FvmiR399a in woodland strawberry. These results showed that PHR1-miR399 module is involved in the regulation of phosphate-signaling pathway and phosphate homeostasis in woodland strawberry.

Molecular Characterization and Expression Profile Analysis of Heat Shock Transcription Factors in Mungbean

Heat shock transcription factors (Hsfs) are essential elements in plant signal transduction pathways that mediate gene expression in response to various abiotic stresses. Mungbean (Vigna radiata) is an important crop worldwide. The emergence of a genome database now allows for functional analysis of mungbean genes. In this study, we dissect the mungbean Hsfs using genome-wide identification and expression profiles. We characterized a total of 24 VrHsf genes and classified them into three groups (A, B, and C) based on their phylogeny and conserved domain structures. All VrHsf genes exhibit highly conserved exon-intron organization, with two exons and one intron. In addition, all VrHsf proteins contain 16 distinct motifs. Chromosome location analysis revealed that VrHsf genes are located on 8 of the 11 mungbean chromosomes, and that seven duplicated gene pairs had formed among them. Moreover, transcription patterns of VrHsf genes varied in different tissues, indicating their different roles in plant growth and development. We identified multiple stress related cis-elements in VrHsf promoter regions 2 kb upstream of the translation initiation codons, and the expression of most VrHsf genes was altered under different stress conditions, suggesting their potential functions in stress resistance pathways. These molecular characterization and expression profile analyses of VrHsf genes provide essential information for further function investigation.

Reference gene selection for qRT-PCR assays in Stellera chamaejasme subjected to abiotic stresses and hormone treatments based on transcriptome datasets

Background

Stellera chamaejasme Linn, an important poisonous plant of the China grassland, is toxic to humans and livestock. The rapid expansion of S. chamaejasme has greatly damaged the grassland ecology and, consequently, seriously endangered the development of animal husbandry. To draft efficient prevention and control measures, it has become more urgent to carry out research on its adaptive and expansion mechanisms in different unfavorable habitats at the genetic level. Quantitative real-time polymerase chain reaction (qRT-PCR) is a widely used technique for studying gene expression at the transcript level; however, qRT-PCR requires reference genes (RGs) as endogenous controls for data normalization and only through appropriate RG selection and qRT-PCR can we guarantee the reliability and robustness of expression studies and RNA-seq data analysis. Unfortunately, little research on the selection of RGs for gene expression data normalization in S. chamaejasme has been reported.

Method

In this study, 10 candidate RGs namely, 18S, 60S, CYP, GAPCP1, GAPDH2, EF1B, MDH, SAND, TUA1, and TUA6, were singled out from the transcriptome database of S. chamaejasme, and their expression stability under three abiotic stresses (drought, cold, and salt) and three hormone treatments (abscisic acid, ABA; gibberellin, GA; ethephon, ETH) were estimated with the programs geNorm, NormFinder, and BestKeeper.

Result

Our results showed that GAPCP1 and EF1B were the best combination for the three abiotic stresses, whereas TUA6 and SAND, TUA1 and CYP, GAPDH2 and 60S were the best choices for ABA, GA, and ETH treatment, respectively. Moreover, GAPCP1 and 60S were assessed to be the best combination for all samples, and 18S was the least stable RG for use as an internal control in all of the experimental subsets. The expression patterns of two target genes (P5CS2 and GI) further verified that the RGs that we selected were suitable for gene expression normalization.

Discussion

This work is the first attempt to comprehensively estimate the stability of RGs in S. chamaejasme. Our results provide suitable RGs for high-precision normalization in qRT-PCR analysis, thereby making it more convenient to analyze gene expression under these experimental conditions.

Root Cell-Specific Regulators of Phosphate-Dependent Growth

Cellular specialization in abiotic stress responses is an important regulatory feature driving plant acclimation. Our in silico approach of iterative coexpression, interaction, and enrichment analyses predicted root cell-specific regulators of phosphate starvation response networks in Arabidopsis (Arabidopsis thaliana). This included three uncharacterized genes termed Phosphate starvation-induced gene interacting Root Cell Enriched (PRCE1, PRCE2, and PRCE3). Root cell-specific enrichment of 12 candidates was confirmed in promoter-GFP lines. T-DNA insertion lines of 11 genes showed changes in phosphate status and growth responses to phosphate availability compared with the wild type. Some mutants (cbl1, cipk2, prce3, and wdd1) displayed strong biomass gain irrespective of phosphate supply, while others (cipk14, mfs1, prce1, prce2, and s6k2) were able to sustain growth under low phosphate supply better than the wild type. Notably, root or shoot phosphate accumulation did not strictly correlate with organ growth. Mutant response patterns markedly differed from those of master regulators of phosphate homeostasis, PHOSPHATE STARVATION RESPONSE1 (PHR1) and PHOSPHATE2 (PHO2), demonstrating that negative growth responses in the latter can be overcome when cell-specific regulators are targeted. RNA sequencing analysis highlighted the transcriptomic plasticity in these mutants and revealed PHR1-dependent and -independent regulatory circuits with gene coexpression profiles that were highly correlated to the quantified physiological traits. The results demonstrate how in silico prediction of cell-specific, stress-responsive genes uncovers key regulators and how their manipulation can have positive impacts on plant growth under abiotic stress.