Comparative transcriptome analysis reveals the regulatory networks of cytokinin in promoting the floral feminization in the oil plant Sapium sebiferum

Comparative RNA-Seq Analysis between Monoecious and Androecious Plants Reveals Regulatory Mechanisms Controlling Female Flowering in Cucurbita pepo

In the monoecious Cucurbita pepo, the transition to female flowering is the time at which the plant starts the production of female flowers after an initial male phase of development. Ethylene plays an essential role in this process since some ethylene deficient and ethylene-insensitive mutants are androecious and only produce male flowers. To gain insight into the molecular mechanisms regulating the specification and early development of female flowers, we have compared the transcriptomic changes occurring in the shoot apices of WT and androecious ethylene-insensitive etr1b mutant plants upon female flowering transition. There were 1160 female flowering-specific DEGs identified in WT plants upon female flowering, and 284 of them were found to be modulated by the ethylene-insensitive etr1b mutation. The function of these DEGs indicated that female flower specification depends on the adoption of a transcriptional program that includes previously identified sex-determining genes in the ethylene pathway, but also genes controlling the biosynthesis and signaling pathways of other phytohormones, and those encoding for many different transcription factors. The transcriptomic changes suggested that gibberellins play a negative role in female flowering, while ethylene, auxins, ABA and cytokinins are positive regulators. Transcription factors from 34 families, including NAC, ERF, bHLH, bZIP, MYB and C2H2/CH3, were found to be regulating female flowering in an ethylene-dependent or -independent manner. Our data open a new perspective of the molecular mechanisms that control the specification and development of female flowers in C. pepo.

Comparative transcriptome analysis reveals candidate genes related to the sex differentiation of Schisandra chinensis

Schisandra chinensis is a monoecious plant with unisex flowers. The fruit of S. chinensis is of high medical with economic value. The yield of S. chinensis fruit is related to the ratio of its female and male flowers. However, there is little research on its floral development and sex differentiation. To elucidate the possible mechanism for the sex differentiation of S. chinensis, we collected 18 samples of female and male flowers from three developmental stages and performed a comparative RNA-seq analysis aimed at identifying differentially expressed genes (DEGs) that may be related to sex differentiation. The results showed 936, 7179, and 6890 differentially expressed genes between female and male flowers at three developmental stages, respectively, and 466 candidate genes may play roles in sex differentiation. KEGG analysis showed genes involved in the flavonoid biosynthesis pathway and DNA replication pathway were essential for the development of female flowers. 51 MADS-box genes and 10 YABBY genes were identified in S. chinensis. The DEGs analysis indicated that MADS-box and YABBY genes were strongly related to the sex determination of S. chinensis. RT-qPCR confirmed the RNA-seq results of 20 differentially expressed genes, including three male-biased genes and 17 female-biased genes. A possible regulatory model of sex differentiation in S. chinensis was proposed according to our results. This study helps reveal the sex-differentiation mechanism of S. chinensis and lays the foundation for regulating the male-female ratio of S. chinensis in the future.

Dynamic cytological and transcriptomic analyses provide novel insights into the mechanisms of sex determination in Castanea henryi

Castanea henryi is a monoecious woody food tree species whose yield and industrialization potential are limited by its low female-to-male flower ratio. Here, the male flowers on the male inflorescence of C. henryi were converted to female flowers by triple applications of exogenous cytokinin (CK) (N-(2-chloro-4-pyridyl)-N'-phenylurea, CPPU). To study the role of exogenous CK in flower sex determination, cytological and transcriptomic analyses were performed on samples from the five stages after CK treatment. Cytological analysis showed that stage 3 (nine days after the last CK treatment) was the critical stage in the differential development of the pistil primordium and stamen primordium. On this basis, one key module and two modules with significant positive correlations with stage 3 were identified by weighted gene co-expression network analysis (WGCNA), combined with transcriptome data. The CK and GA biosynthesis- and signaling-related genes, three transcription factor (TF) families, and 11 floral organ identity genes were identified in the related modules. In particular, the TFs WRKY47, ERF021, and MYB4, and floral organ identity genes AGL11/15, DEF, and SEP1 with large differences are considered to be critical regulators of sex determination in C. henryi. Based on these results, a genetic regulatory network for exogenous CK in the sex determination of flowers in C. henryi is proposed. This study contributes to the understanding of the role of CK in the sex regulation of flowers and provides new insights into the regulatory network of sex determination in C. henryi.

Comparative transcriptomics analysis of developing peanut (Arachis hypogaea L.) pods reveals candidate genes affecting peanut seed size

Pod size is one of the most important agronomic features of peanuts, which directly affects peanut yield. Studies on the regulation mechanism underpinning pod size in cultivated peanuts remain hitherto limited compared to model plant systems. To better understand the molecular elements that underpin peanut pod development, we conducted a comprehensive analysis of chronological transcriptomics during pod development in four peanut accessions with similar genetic backgrounds, but varying pod sizes. Several plant transcription factors, phytohormones, and the mitogen-activated protein kinase (MAPK) signaling pathways were significantly enriched among differentially expressed genes (DEGs) at five consecutive developmental stages, revealing an eclectic range of candidate genes, including PNC, YUC, and IAA that regulate auxin synthesis and metabolism, CYCD and CYCU that regulate cell differentiation and proliferation, and GASA that regulates seed size and pod elongation via gibberellin pathway. It is plausible that MPK3 promotes integument cell division and regulates mitotic activity through phosphorylation, and the interactions between these genes form a network of molecular pathways that affect peanut pod size. Furthermore, two variant sites, GCP4 and RPPL1, were identified which are stable at the QTL interval for seed size attributes and function in plant cell tissue microtubule nucleation. These findings may facilitate the identification of candidate genes that regulate pod size and impart yield improvement in cultivated peanuts.

Integrated transcriptome and plant growth substance profiles to identify the regulatory factors involved in floral sex differentiation in Zanthoxylum armatum DC

Zanthoxylum armatum is a prominent plant for food industries. Its male flowers often occur in gynogenesis plants; however, the potential mechanism remains poorly understood. Herein, a total of 26 floral sex differentiation stages were observed to select four vital phases to reveal key factors by using RNA-seq, phytohormones and carbohydrates investigation. The results showed that a selective abortion of stamen or pistil primordia could result in the floral sex differentiation in Z. armatum. Carbohydrates might collaborate with cytokinin to effect the male floral differentiation, whereas female floral differentiation was involved in SA, GA₁, and ABA biosynthesis and signal transduction pathways. Meanwhile, these endogenous regulators associated with reproductive growth might be integrated into ABCDE model to regulate the floral organ differentiation in Z. armatum. Furthermore, the 21 crucial candidates were identified in co-expression network, which would contribute to uncovering their roles in floral sex differentiation of Z. armatum in further studies. To the best of our knowledge, this study was the first comprehensive investigation to link floral sex differentiation with multi-level endogenous regulatory factors in Z. armatum. It also provided new insights to explore the regulatory mechanism of floral sex differentiation, which would be benefited to cultivate high-yield varieties in Z. armatum.

Transcriptome Analysis Reveals the Role of GA3 in Regulating the Asynchronism of Floral Bud Differentiation and Development in Heterodichogamous Cyclocarya paliurus (Batal.) Iljinskaja

Cyclocarya paliurus is an important medical plant owing to the diverse bioactive compounds in its leaves. However, the heterodichogamy with female and male functions segregation within protandry (PA) or protogyny (PG) may greatly affect seed quality and its plantations for medicinal use. To speculate on the factor playing the dominant role in regulating heterodichogamy in C. paliurus, based on phenotypic observations, our study performed a multi comparison transcriptome analysis on female and male buds (PG and PA types) using RNA-seq. For the female and male bud comparisons, a total of 6753 differentially expressed genes (DEGs) were detected. In addition, functional analysis revealed that these DEGs were significantly enriched in floral development, hormone, and GA-related pathways. As the dominant hormones responsible for floral differentiation and development, gibberellins (GAs) in floral buds from PG and PA types were quantified using HPLC-MS. Among the tested GAs, GA₃ positively regulated the physiological differentiation (S0) and germination (S2) of floral buds. The dynamic changes of GA₃ content and floral morphological features were consistent with the expression levels of GA-related genes. Divergences of GA₃ contents at S0 triggered the asynchronism of physiological differentiation between male and female buds of intramorphs (PA-M vs. PA-F and PG-F vs. PG-M). A significant difference in GA₃ content enlarged this asynchronism at S2. Thus, we speculate that GA₃ plays the dominant role in the formation of heterodichogamy in C. paliurus. Meanwhile, the expression patterns of GA-related DEGs, including CPS, KO, GA20ox, GA2_OX, GID1, and DELLA genes, which play central roles in regulating flower development, coincided with heterodichogamous characteristics. These results support our speculations well, which should be further confirmed.

Karrikinolide alleviates salt stress in wheat by regulating the redox and K+/Na+ homeostasis

Karrikinolide (KAR₁), identified in biochars, has gained research attention because of its significant role in seed germination, seedling development, root development, and abiotic stresses. However, KAR₁ regulation of salt stress in wheat is elusive. This study investigated the physiological mechanism involved in KAR₁ alleviation of salt stress in wheat. The results showed KAR₁ boosted seed germination percentage under salinity stress via stimulating the relative expression of genes regulating gibberellins biosynthesis and decreasing the expression levels of abscisic acid biosynthesis and signaling genes. As seen in seed germination, exogenous supplementation of KAR₁ dramatically mitigated the salt stress also in wheat seedling, resulting in increased root and shoot growth as measured in biomass as compared to salt stress alone. Salt stress significantly induced the endogenous hydrogen peroxide and malondialdehyde levels, whereas KAR₁ strictly counterbalanced them. Under salt stress, KAR₁ supplementation showed significant induction in reduced glutathione (GSH) and reduction in oxidized glutathione (GSSG) content, which improved GSH/GSSG ratio in wheat seedlings. Exogenous supplementation of KAR₁ significantly promoted the activities of enzymatic antioxidants in wheat seedlings exposed to salt stress. KAR₁ induced the relative expression of genes regulating the biosynthesis of antioxidants in wheat seedlings under salinity. Moreover, KAR₁ induced the expression level of K⁺/Na⁺ homeostasis genes, reduced Na⁺ concentration, and induced K⁺ concentration in wheat seedling under salt stress. The results suggest that KAR₁ supplementation maintained the redox and K⁺/Na⁺ homeostasis in wheat seedling under salinity, which might be a crucial part of physiological mechanisms in KAR₁ induced tolerance to salt stress. In conclusion, we exposed the protective role of KAR₁ against salt stress in wheat.

Screening and identification of miRNAs related to sexual differentiation of strobili in Ginkgo biloba by integration analysis of small RNA, RNA, and degradome sequencing

BACKGROUND

Ginkgo biloba, a typical dioecious plant, is a traditional medicinal plant widely planted. However, it has a long juvenile period, which severely affected the breeding and cultivation of superior ginkgo varieties.

RESULTS

In order to clarify the complex mechanism of sexual differentiation in G. biloba strobili. Here, a total of 3293 miRNAs were identified in buds and strobili of G. biloba, including 1085 known miRNAs and 2208 novel miRNAs using the three sequencing approaches of transcriptome, small RNA, and degradome. Comparative transcriptome analysis screened 4346 and 7087 differentially expressed genes (DEGs) in male buds (MB) _vs_ female buds (FB) and microstrobilus (MS) _vs_ ovulate strobilus (OS), respectively. A total of 6032 target genes were predicted for differentially expressed miRNA. The combined analysis of both small RNA and transcriptome datasets identified 51 miRNA-mRNA interaction pairs that may be involved in the process of G. biloba strobili sexual differentiation, of which 15 pairs were verified in the analysis of degradome sequencing.

CONCLUSIONS

The comprehensive analysis of the small RNA, RNA and degradome sequencing data in this study provided candidate genes and clarified the regulatory mechanism of sexual differentiation of G. biloba strobili from multiple perspectives.

Gene expression trajectories during male and female reproductive development in balsam poplar (Populus balsamifera L.)

Plant reproductive development from the first appearance of reproductively committed axes through to floral maturation requires massive and rapid remarshalling of gene expression. In dioecious species such as poplar this is further complicated by divergent male and female developmental programs. We used seven time points in male and female balsam poplar (Populus balsamifera L.) buds and catkins representing the full annual flowering cycle, to elucidate the effects of time and sex on gene expression during reproductive development. Time (developmental stage) is dominant in patterning gene expression with the effect of sex nested within this. Here, we find (1) evidence for five successive waves of alterations to the chromatin landscape which may be important in setting the overall reproductive trajectory, regardless of sex. (2) Each individual developmental stage is further characterized by marked sex-differential gene expression. (3) Consistent sexually differentiated gene expression regardless of developmental stage reveal candidates for high-level regulators of sex and include the female-specific poplar ARR17 homologue. There is also consistent male-biased expression of the MADS-box genes PISTILLATA and APETALA3. Our work provides insights into expression trajectories shaping reproductive development, its potential underlying mechanisms, and sex-specific translation of the genome information into reproductive structures in balsam poplar.

Karrikin Improves Osmotic and Salt Stress Tolerance via the Regulation of the Redox Homeostasis in the Oil Plant Sapium sebiferum

Karrikins are reported to stimulate seed germination, regulate seedling growth, and increase the seedling vigor in abiotic stress conditions in plants. Nevertheless, how karrikins alleviate abiotic stress remains largely elusive. In this study, we found that karrikin (KAR1) could significantly alleviate both drought and salt stress in the important oil plant Sapium sebiferum. KAR1 supplementation in growth medium at a nanomolar (nM) concentration was enough to recover seed germination under salt and osmotic stress conditions. One nanomolar of KAR1 improved seedling biomass, increased the taproot length, and increased the number of lateral roots under abiotic stresses, suggesting that KAR1 is a potent alleviator of abiotic stresses in plants. Under abiotic stresses, KAR1-treated seedlings had a higher activity of the key antioxidative enzymes, such as superoxide dismutase, peroxidase, catalase, and ascorbate peroxidase, in comparison with the control, which leads to a lower level of hydrogen peroxide, malondialdehyde, and electrolyte leakage. Moreover, the metabolome analysis showed that KAR1 treatment significantly increased the level of organic acids and amino acids, which played important roles in redox homeostasis under stresses, suggesting that karrikins might alleviate abiotic stresses via the regulation of redox homeostasis. Under abiotic stresses, applications of karrikins did not increase the endogenous abscisic acid level but altered the expression of several ABA signaling genes, such as SNF1-RELATED PROTEIN KINASE2.3, SNF1-RELATED PROTEIN KINASE2.6, ABI3, and ABI5, suggesting potential interactions between karrikins and ABA signaling in the stress responses. Conclusively, we not only provided the physiological and molecular evidence to clarify the mechanism of karrikins in the regulation of stress adaptation in S. sebiferum but also showed the potential value of karrikins in agricultural practices, which will lay a foundation for further studies about the role of karrikins in abiotic stress alleviation in plants.

Flower-Specific Overproduction of Cytokinins Altered Flower Development and Sex Expression in the Perennial Woody Plant Jatropha curcas L

Jatropha curcas L. is monoecious with a low female-to-male ratio, which is one of the factors restricting its seed yield. Because the phytohormone cytokinins play an essential role in flower development, particularly pistil development, in this study, we elevated the cytokinin levels in J. curcas flowers through transgenic expression of a cytokinin biosynthetic gene (AtIPT4) from Arabidopsis under the control of a J. curcas orthologue of TOMATO MADS BOX GENE 6 (JcTM6) promoter that is predominantly active in flowers. As expected, the levels of six cytokinin species in the inflorescences were elevated, and flower development was modified without any alterations in vegetative growth. In the transgenic J. curcas plants, the flower number per inflorescence was significantly increased, and most flowers were pistil-predominantly bisexual, i.e., the flowers had a huge pistil surrounded with small stamens. Unfortunately, both the male and the bisexual flowers of transgenic J. curcas were infertile, which might have resulted from the continuously high expression of the transgene during flower development. However, the number and position of floral organs in the transgenic flowers were well defined, which suggested that the determinacy of the floral meristem was not affected. These results suggest that fine-tuning the endogenous cytokinins can increase the flower number and the female-to-male ratio in J. curcas.

Vitis flower types: from the wild to crop plants

Vitis vinifera can be divided into two subspecies, V. vinifera subsp. vinifera, one of the most important agricultural crops in the world, and its wild ancestor, V. vinifera subsp. sylvestris. Three flower types can be observed: hermaphrodite and female (on some varieties) in vinifera, and male or female flowers in sylvestris. It is assumed that the different flower types in the wild ancestor arose through specific floral patterns of organ abortion. A considerable amount of data about the diversity of sexual systems in grapevines has been collected over the past century. Several grapevine breeding studies led to the hypothesis that dioecy in vinifera is derived from a hermaphrodite ancestor and could be controlled by either, one or two linked genetic determinants following Mendelian inherence. More recently, experiments using molecular approaches suggested that these loci were located in a specific region of the chromosome 2 of vinifera. Based on the works published so far, its seems evident that a putative sex locus is present in chromosome 2. However, it is still not fully elucidated whether flower types are regulated by two linked loci or by one locus with three alleles. Nevertheless, several genes could contribute to sex determination in grapevine. This review presents the results from early studies, combined with the recent molecular approaches, which may contribute to the design of new experiments towards a better understanding of the sex inheritance in grapevine.

Overexpression of the Stress-Inducible SsMAX2 Promotes Drought and Salt Resistance via the Regulation of Redox Homeostasis in Arabidopsis

Recent studies have demonstrated that strigolactones (SLs) also participate in the regulation of stress adaptation; however, the regulatory mechanism remains elusive. In this study, the homolog of More Axillary Branches 2, which encodes a key component in SL signaling, in the perennial oil plant Sapium sebiferum was identified and functionally characterized in Arabidopsis. The results showed that the expression of SsMAX2 in S. sebiferum seedlings was stress-responsive, and SsMAX2 overexpression (OE) in Arabidopsis significantly promoted resistance to drought, osmotic, and salt stresses. Moreover, SsMAX2 OE lines exhibited decreased chlorophyll degradation, increased soluble sugar and proline accumulation, and lower water loss ratio in response to the stresses. Importantly, anthocyanin biosynthesis and the activities of several antioxidant enzymes, such as superoxide dismutase (SOD), peroxidase (POD), and ascorbate peroxidase (APX), were enhanced in the SsMAX2 OE lines, which further led to a significant reduction in hydrogen peroxide levels. Additionally, the SsMAX2 OE lines exhibited higher expression level of several abscisic acid (ABA) biosynthesis genes, suggesting potential interactions between SL and ABA in the regulation of stress adaptation. Overall, we provide physiological and biochemical evidence demonstrating the pivotal role of SsMAX2 in the regulation of osmotic, drought, and salt stress resistance and show that MAX2 can be a genetic target to improve stress tolerance.

Transcriptome Profiling of the Elongating Internode of Cotton (Gossypium hirsutum L.) Seedlings in Response to Mepiquat Chloride

The plant growth retardant mepiquat chloride (MC) has been extensively used to produce compact plant canopies and increase yield in cotton (Gossypium hirsutum L.). Previous studies showed that MC reduced plant height and internode length by inhibiting GA biosynthesis and cell elongation. However, whether there are other molecular mechanisms underlying MC-induced growth retardation has remained largely unknown. In the present study, we conducted histological, transcriptomic, and phytohormone analyses of the second elongating internodes of cotton seedlings treated with MC. Histological analysis revealed that the MC shortened the internodes through suppressing both cell division and cell elongation. Consistent with the observed phenotype, many genes related to cell growth were significantly downregulated by MC. Transcriptome profiling showed that the expression of genes related not only to GA, but also to auxin, brassinosteroid (BR), and ethylene metabolism and signaling was remarkably suppressed, whereas that of genes related to cytokinin (CK) and abscisic acid (ABA) metabolism was induced by MC. Consistent with the expression pattern, significant decrease of endogenous GA, auxin, and BR content, but an increase in CK content was observed after MC treatment. Most of these hormone related genes displayed opposite regulation pattern by exogenous GA₃ treatment compared to MC; simultaneous application of MC and GA₃ could alleviate the genes expression changes induced by MC treatment, indicating MC does not directly affect other plant hormones, but through the inhibition of the GA biosynthesis. In addition, the expression of genes related to secondary metabolism and many transcription factors (TFs) were differentially regulated by MC. In summary, we confirmed the important role of GA in MC-induced growth inhibition of cotton, and further found that other hormones were also involved in this process in a GA-dependent manner. This study provides novel insights into the molecular mechanism underlying the MC-mediated inhibition of internode elongation in cotton seedlings.