Comparative transcriptomics and weighted gene co-expression correlation network analysis (WGCNA) reveal potential regulation mechanism of carotenoid accumulation in Chrysanthemum × morifolium

Exploring the molecular basis of efficient feed utilization in low residual feed intake slow-growing ducks based on breast muscle transcriptome

Residual feed intake (RFI) has recently gained attention as a key indicator of feed efficiency in poultry. In this study, 800 slow-growing ducks with similar initial body weights were reared in an experimental facility until they were culled at 42 d of age. Thirty high RFI (HRFI) and 30 low RFI (LRFI) birds were selected to evaluate their growth performance, carcass characteristics, and muscle development. Transcriptome and weighted gene co-expression correlation network analyses of pectoral muscles were conducted on six LRFI and six HRFI ducks. The results revealed that selecting for LRFI significantly reduced feed consumption (P < 0.05) and improved feed efficiency without affecting the growth performance, slaughter rate, or meat quality of ducks (P > 0.05). Moreover, compared with HRFI ducks, LRFI ducks had a lower pectoral muscle fat content (P < 0.05), larger muscle fiber diameter and area (P < 0.05), and lower muscle fiber density (P < 0.05). There were significant differences in gene expression between LRFI and HRFI ducks, with 102 upregulated and 258 downregulated genes, which were enriched in the PPAR signaling pathway, adipocytokine signaling pathway, actin cytoskeleton regulation, ECM-receptor interaction, and focal adhesion. The expression of genes associated with fat and energy metabolism, including ACSL6, PCK1, APOC3, HMGCS2, PRKAG3, and G6PC1, was downregulated in LRFI ducks, and weighted gene co-expression correlation network analysis identified PRKAG3 as a hub gene. Our findings indicate that reduced mitochondrial energy metabolism may contribute to the RFI of slow-growing ducks, with PRKAG3 playing a pivotal role in this biological process. These findings provide novel insights into the molecular changes underlying RFI variation in slow-growing ducks.

Unveiling Salt Tolerance Mechanisms and Hub Genes in Alfalfa (Medicago sativa L.) Through Transcriptomic and WGCNA Analysis

Alfalfa (Medicago sativa L.), an important forage crop with high nutritional value and good palatability, plays a vital role in the development of animal husbandry in China. In Northeast China, there are vast areas of saline-alkali land that remain undeveloped. Given that alfalfa is a highly adaptable forage crop, exploring its salt tolerance at the molecular transcriptional level and identifying salt-tolerant genes has great significance for breeding salt-resistant alfalfa varieties. This also provides valuable genetic resources for better utilization of saline-alkali land. In this study, we conducted two rounds of screening on 41 alfalfa varieties and identified WL168 as a salt-sensitive variety and Longmu801 as a salt-tolerant variety. After 7 days of 300 mM salt stress, both varieties showed a decreasing trend in plant height, fresh weight, and dry weight over time, but Longmu801 demonstrated better water retention ability compared to WL168. Chlorophyll content also declined, but chlorophyll a and total chlorophyll levels in Longmu801 were higher than in WL168. Hydrogen peroxide and malondialdehyde levels increased overall, but Longmu801 had significantly lower levels than WL168 under prolonged stress. Both varieties showed increasing trends in soluble sugars, proline, and antioxidant enzymes (SOD, POD, CAT), with Longmu801 significantly outperforming WL168. This suggests that the two varieties share similar growth and physiological response mechanisms, with their differences primarily arising from variations in indicator levels. In the above, comparisons between varieties were conducted based on the relative values of the indicators in relation to their controls. Transcriptomic analysis revealed that under salt stress, Longmu801 had 16,485 differentially expressed genes (DEGs) relative to its control, while WL168 had 18,726 DEGs compared to its control. Among these, 2164 DEGs shared the same expression trend, with GO functions enriched in response to oxidative stress, nucleus, plasma membrane, and others. The KEGG pathways were enriched in phenylpropanoid biosynthesis, protein processing in the endoplasmic reticulum, starch and sucrose metabolism, and others. This suggests that alfalfa's transcriptional response mechanism to salt stress involves these pathways. Additionally, the variety-specific DEGs were also enriched in the same KEGG pathways and GO functions, indicating that the differences between the two varieties stem from their unique stress-responsive DEGs, while their overall mechanisms for coping with stress remain similar. To further identify salt stress-related genes, this study conducted WGCNA analysis using 32,683 genes and physiological indicators. Six modules closely related to physiological traits were identified, and the top five genes ranked by degree in each module were selected as hub genes. Further analysis of these hub genes identified five genes directly related to salt stress: Msa085011, Msa0605650, Msa0397400, Msa1258740, and Msa0958830. Mantel test analysis revealed that these genes showed strong correlations with physiological indicators. This study will provide important insights for breeding salt-tolerant alfalfa varieties.

Integrated transcriptome and endogenous hormone analyses reveal the factors affecting the yield of Camellia oleifera

Camellia oleifera is an important woody oil tree in China, in which the flowers and fruits appear during the same period. The endogenous hormone changes and transcription expression levels in different parts of the flower tissue (sepals, petals, stamens, and pistils), flower buds, leaves, and seeds of Changlin 23 high-yield (H), Changlin low-yield (L), and control (CK) C. oleifera groups were studied. The abscisic acid (ABA) content in the petals and stamens in the L group was significantly higher than that in the H and CK groups, which may be related to flower and fruit drops. The high N6-isopentenyladenine (iP) and indole acetic acid (IAA) contents in the flower buds may be associated with a high yield. Comparative transcriptome analysis showed that the protein phosphatase 2C (PP2C), jasmonate-zim-domain protein (JAZ), and WRKY-related differentially expressed genes (DEGs) may play an important role in determining leaf color. Gene set enrichment analysis (GSEA) comparison showed that jasmonic acid (JA) and cytokinin play an important role in determining the pistil of the H group. In this study, endogenous hormone and transcriptome analyses were carried out to identify the factors influencing the large yield difference in C. oleifera in the same year, which provides a theoretical basis for C. oleifera in the future.

Integrated metabolomic and transcriptomic analysis reveals biosynthesis mechanism of flavone and caffeoylquinic acid in chrysanthemum

BACKGROUND

Chrysanthemum morifolium 'HangBaiJu', a popular medicinal and edible plant, exerts its biological activities primarily through the presence of flavones and caffeoylquinic acids (CQAs). However, the regulatory mechanism of flavone and CQA biosynthesis in the chrysanthemum capitulum remains unclear.

RESULTS

In this study, the content of flavones and CQAs during the development of chrysanthemum capitulum was determined by HPLC, revealing an accumulation pattern with higher levels at S1 and S2 and a gradual decrease at S3 to S5. Transcriptomic analysis revealed that CmPAL1/2, CmCHS1/2, CmFNS, CmHQT, and CmHCT were key structural genes in flavones and CQAs biosynthesis. Furthermore, weighted gene co-expression correlation network analysis (WGCNA), k-means clustering, correlation analysis and protein interaction prediction were carried out in this study to identify transcription factors (TFs) associated with flavone and CQA biosynthesis, including MYB, bHLH, AP2/ERF, and MADS-box families. The TFs CmERF/PTI6 and CmCMD77 were proposed to act as upstream regulators of CmMYB3 and CmbHLH143, while CmMYB3 and CmbHLH143 might form a complex to directly regulate the structural genes CmPAL1/2, CmCHS1/2, CmFNS, CmHQT, and CmHCT, thereby controlling flavone and CQA biosynthesis.

CONCLUSIONS

Overall, these findings provide initial insights into the TF regulatory network underlying flavones and CQAs accumulation in the chrysanthemum capitulum, which laid a theoretical foundation for the quality improvement of C. morifolium 'HangBaiJu' and the high-quality development of the industry.

Co-regulatory network analysis of the main secondary metabolite (SM) biosynthesis in Crocus sativus L

Saffron (Crocus sativus L.) is being embraced as the most important medicinal plant and the commercial source of saffron spice. Despite the beneficial economic and medicinal properties of saffron, the regulatory mechanism of the correlation of TFs and genes related to the biosynthesis of the apocarotenoids pathway is less obvious. Realizing these regulatory hierarchies of gene expression networks related to secondary metabolites production events is the main challenge owing to the complex and extensive interactions between the genetic behaviors. Recently, high throughput expression data have been highly feasible for constructing co-regulation networks to reveal the regulated processes and identifying novel candidate hub genes in response to complex processes of the biosynthesis of secondary metabolites. Herein, we performed Weighted Gene Co-expression Network Analysis (WGCNA), a systems biology method, to identify 11 regulated modules and hub TFs related to secondary metabolites. Three specialized modules were found in the apocarotenoids pathway. Several hub TFs were identified in notable modules, including MADS, C2H2, ERF, bZIP, HD-ZIP, and zinc finger protein MYB and HB, which were potentially associated with apocarotenoid biosynthesis. Furthermore, the expression levels of six hub TFs and six co-regulated genes of apocarotenoids were validated with RT-qPCR. The results confirmed that hub TFs specially MADS, C2H2, and ERF had a high correlation (P < 0.05) and a positive effect on genes under their control in apocarotenoid biosynthesis (CCD2, GLT2, and ADH) among different C. sativus ecotypes in which the metabolite contents were assayed. Promoter analysis of the co-expressed genes of the modules involved in apocarotenoids biosynthesis pathway suggested that not only are the genes co-expressed, but also share common regulatory motifs specially related to hub TFs of each module and that they may describe their common regulation. The result can be used to engineer valuable secondary metabolites of C. sativus by manipulating the hub regulatory TFs.

Analysis of huanglongbing-associated RNA-seq data reveals disturbances in biological processes within Citrus spp. triggered by Candidatus Liberibacter asiaticus infection

Introduction

Huanglongbing (HLB), a disease that's ubiquitous worldwide, wreaks havoc on the citrus industry. The primary culprit of HLB is the gram-negative bacterium Candidatus Liberibacter asiaticus (CLas) that infects the phloem, but its damaging mechanism is yet to be fully understood.

Methods and results

In this study, a multitude of tools including weighted correlation network analysis (WGCNA), protein-protein interaction (PPI) network analysis and gene expression profiling are employed to unravel the intricacies of its pathogenesis. The investigation pinpoints various central genes, such as the ethylene-responsive transcription factor 9 (ERF9) and thioredoxin reductase 1 (TrxR1), that are associated with CLas invasion and resultant disturbances in numerous biological operations. Additionally, the study uncovers a range of responses through the detection of differential expressed genes (DEGs) across different experiments. The discovery of core DEGs leads to the identification of pivotal genes such as the sieve element occlusion (SEO) and the wall-associated receptor kinase-like 15 (WAKL15). PPI network analysis highlights potential vital proteins, while GO and KEGG pathway enrichment analysis illustrate a significant impact on multiple defensive and metabolic pathways. Gene set enrichment analysis (GSEA) indicates significant alterations in biological processes such as leaf senescence and response to biotic stimuli.

Discussion

This all-encompassing approach extends valuable understanding into the pathogenesis of CLas, potentially aiding future research and therapeutic strategies for HLB.

Comparative analysis of differential gene expression reveals novel insights into the heteroblastic foliage functional traits of Pinus massoniana seedlings

Pinus massoniana needles, rich in medicinal polysaccharides and flavonoids, undergo heteroblastic foliage, transitioning from primary needles (PN) to secondary needles (SN) during growth, resulting in altered functional traits. Despite its significance, the molecular regulatory mechanisms governing these traits remain unclear. This study employs Iso-Seq and RNA-Seq analyses to explore differentially expressed genes (DEGs) associated with functional traits throughout the main growth season of heteroblastic foliage. Co-expression network analysis identified 34 hub genes and 17 key transcription factors (TFs) influencing light-harvesting antenna, photosystem I and II, crucial in photosynthesis regulation. Additionally, 14 genes involved in polysaccharide metabolism pathways, synthesizing sucrose, glucose, UDP sugars, and xylan, along with four genes in flavonoid biosynthesis pathways, regulating p-coumaroyl-CoA, quercetin, galangin, and myricetin production, exhibited differential expression between PN and SN. Further analysis unveils a highly interconnected network among these genes, forming a pivotal cascade of TFs and DEGs. Therefore, heteroblastic changes significantly impact needle functional traits, potentially affecting the pharmacological properties of PN and SN. Thus, these genomic insights into understanding the molecular-level differences of heteroblastic foliage, thereby establishing a foundation for advancements in the pharmaceutical industry related to needle-derived products.

Transcriptomic and functional analyzes reveal that the brassinosteroid insensitive 1 receptor (OsBRI1) regulates cold tolerance in rice

Brassinosteroids (BR) play crucial roles in plant development and abiotic stresses in plants. Exogenous application of BR can significantly enhance cold tolerance in rice. However, the regulatory relationship between cold tolerance and the BR signaling pathway in rice remains largely unknown. Here, we characterized functions of the BR receptor OsBRI1 in response to cold tolerance in rice using its loss-of-function mutant (d61-1). Our results showed that mutant d61-1 was less tolerant to cold stress than wild-type (WT). Besides, d61-1 had lower levels than WT for some physiological parameters, including catalase activity (CAT), superoxide dismutase activity (SOD), peroxidase activity (POD), peroxidase activity (PRO), soluble protein, and soluble sugar content, while malondialdehyde content (MDA) and relative electrical conductivity (REC) levels in d61-1 were higher than those in WT plants. These results indicated that the loss of OsBRI1 function resulted in decreased cold tolerance in rice. In addition, we performed RNA sequencing (RNA-seq) of WT and d61-1 mutant under cold stress. Numerous common and unique differentially expressed genes (DEGs) with up- and down-regulation were observed in WT and d61-1 mutant. Some DEGs were expressed to various degrees, even opposite, between CK1 vs. T1 (WT) and CK2 vs. T2 (d61-1). Among these specific DEGs, some typical genes are involved in plant tolerance to cold stress. Through weighted correlation network analysis (WGCNA), 50 hub genes were screened in the turquoise and blue module. Many genes were involved in cold stress and plant hormone, such as Os01g0279800 (BRI1-associated receptor kinase 1 precursor), Os10g0513200 (Dwarf and tiller-enhancing 1, DTE1), Os02g0706400 (MYB-related transcription factor, OsRL3), etc. Differential expression levels of some genes were verified in WT and d61-1 under cold stress using qRT-PCR. These valuable findings and gene resources will be critical for understanding the regulatory relationships between cold stress tolerance and the BR signaling pathways in rice.

Transcriptomic and Phenotypic Analyses Reveal the Molecular Mechanism of Dwarfing in Tetraploid Robinia pseudoacacia L

Polyploid breeding techniques aid in the cultivation of new forestry cultivars, thus expanding the suite of strategies for the improvement of arboreal traits and innovation within the field of forestry. Compared to diploid Robinia pseudoacacia L. (black locust) 'D26-5①' (2×), its dwarfed homologous tetraploid 'D26-5②' (4×) variety has better application prospects in garden vegetation guardrails and urban landscape. However, the molecular mechanism of the generation and growth of this dwarf variety is still unclear. Here, plant growth and development as well as histological differences between the diploid and its autotetraploid were investigated. Levels of endogenous hormones at three different developmental stages (20, 40, and 70 days) of 2× and homologous 4× tissue culture plantlets were assessed, and it was found that the brassinosteroid (BR) contents of the former were significantly higher than the latter. Transcriptome sequencing data analysis of 2× and homologous 4× showed that differentially expressed genes (DEGs) were significantly enriched in plant hormone synthesis and signal transduction, sugar and starch metabolism, and the plant circadian rhythm pathway, which are closely related to plant growth and development. Therefore, these biological pathways may be important regulatory pathways leading to dwarfism and slow growth in tetraploids. Additionally, utilizing weighted gene coexpression network analysis (WGCNA), we identified three crucial differentially expressed genes (DEGs)-PRR5, CYP450, and SPA1-that potentially underlie the observed ploidy variation. This study provides a new reference for the molecular mechanism of dwarfism in dwarfed autotetraploid black locusts. Collectively, our results of metabolite analysis and comparative transcriptomics confirm that plant hormone signaling and the circadian rhythm pathway result in dwarfism in black locusts.

Involvement of Transcription Factors and Regulatory Proteins in the Regulation of Carotenoid Accumulation in Plants and Algae

Carotenoids are essential for photosynthesis and photoprotection in photosynthetic organisms, which are widely used in food coloring, feed additives, nutraceuticals, cosmetics, and pharmaceuticals. Carotenoid biofortification in crop plants or algae has been considered as a sustainable strategy to improve human nutrition and health. However, the regulatory mechanisms of carotenoid accumulation are still not systematic and particularly scarce in algae. This article focuses on the regulatory mechanisms of carotenoid accumulation in plants and algae through regulatory factors (transcription factors and regulatory proteins), demonstrating the complexity of homeostasis regulation of carotenoids, mainly including transcriptional regulation as the primary mechanism, subsequent post-translational regulation, and cross-linking with other metabolic processes. Different organs of plants and different plant/algal species usually have specific regulatory mechanisms for the biosynthesis, storage, and degradation of carotenoids in response to the environmental and developmental signals. In plants and algae, regulators such as MYB, bHLH, MADS, bZIP, AP2/ERF, WRKY, and orange proteins can be involved in the regulation of carotenoid metabolism. And many more regulators, regulatory networks, and mechanisms need to be explored. Our paper will provide a basis for multitarget or multipathway engineering for carotenoid biofortification in plants and algae.

Comprehensive transcriptome and WGCNA analysis reveals the potential function of anthocyanins in low-temperature resistance of a red flower mutant tobacco

The anthocyanin is a protective substance in various plants, and plays important roles in resisting to low-temperature. Here, we explored transcriptome analysis of pink flower (as CK) and the natural mutant red flower (as research objects) under low-temperature conditions, and aimed to reveal the potential functions of anthocyanins and anthocyanin-related regulatory factors in resistance to low-temperature. Our results showed that most of the differentially expressed genes (DEGs) encoding key enzymes in the late stage of anthocyanin metabolism in the mutant were significantly up-regulated. Meanwhile, several genes significantly differentially expressed in CK or mutant were obtained by classification and analysis of transcription factors (TFs), phytohormones and osmoregulators. Additionally, WGCNA was carried out to mine hub genes resistanted to low-temperature stress in flavonoid pathway. Finally, one UFGT family gene, three MYB and one bHLH were obtained as the future hub genes of this study. Combined with the above information, we concluded that the ability of the red flower mutant to grow and develop normally at low-temperatures was the result of a combination of flavonoids and cold resistance genes.

Integration of comparative transcriptomics and WGCNA characterizes the regulation of anthocyanin biosynthesis in mung bean (Vigna radiata L.)

Mung bean is a dual-use crop widely cultivated in Southeast Asia as a food and medicine resource. The development of new functional mung bean varieties demands identifying new genes regulating anthocyanidin synthesis and investigating their molecular mechanism. In this study, we used high-throughput sequencing technology to generate transcriptome sequence of leaves, petioles, and hypocotyls for investigating the anthocyanins accumulation in common mung bean variety as well as anthocyanidin rich mung bean variety, and to elucidate their molecular mechanisms. 29 kinds of anthocyanin compounds were identified. Most of the anthocyanin components contents were significantly higher in ZL23 compare with AL12. Transcriptome analysis suggested that a total of 93 structural genes encoding the anthocyanin biosynthetic pathway and 273 regulatory genes encoding the ternary complex of MYB-bHLH-WD40 were identified, of which 26 and 78 were differentially expressed in the two varieties. Weighted gene co-expression network analysis revealed that VrMYB3 and VrMYB90 might have enhanced mung bean anthocyanin content by inducing the expression of structural genes such as PAL, 4CL, F3'5'H, LDOX, and F3'H, which was consistent with qRT-PCR results. These findings are envisaged to provide a reference for studying the molecular mechanism of anthocyanin accumulation in mung beans.

Phytochemicals, therapeutic benefits and applications of chrysanthemum flower: A review

Chrysanthemum is a flowering plant belonging to a genus of the dicotyledonous herbaceous annual flowering plant of the Asteraceae (Compositae) family. It is a perpetual flowering plant, mostly cultivated for medicinal purposes; generally, used in popular drinks due to its aroma and flavor. It is primarily cultivated in China, Japan, Europe, and United States. These flowers were extensively used in various healthcare systems and for treating various diseases. Chrysanthemum flowers are rich in phenolic compounds and exhibit strong properties including antioxidant, antimicrobial, anti-inflammatory, anticancer, anti-allergic, anti-obesity, immune regulation, hepatoprotective, and nephroprotective activities. The main aim of the present review was to investigate the nutritional profile, phytochemistry, and biological activities of flowers of different Chrysanthemum species. Also, a critical discussion of the diverse metabolites or bioactive constituents of the Chrysanthemum flowers is highlighted in the present review. Moreover, the flower extracts of Chrysanthemum have been assessed to possess a rich phytochemical profile, including compounds such as cyanidin-3-O-(6″-O-malonyl) glucoside, delphinidin 3-O-(6" -O-malonyl) glucoside-3', rutin, quercetin, isorhamnetin, rutinoside, and others. These profiles exhibit potential health benefits, leading to their utilization in the production of supplementary food products and pharmaceutical drugs within the industry. However, more comprehensive research studies/investigations are still needed to further discover the potential benefits for human and animal utilization.

Integration of transcriptomic and metabolomic analysis unveils the response mechanism of sugar metabolism in Cyclocarya paliurus seedlings subjected to PEG-induced drought stress

Cyclocarya paliurus (Batal.) Iljinskaja is a multiple function tree species used for functional food and valued timber production. Carbohydrates, especially water-soluble carbohydrates, play an important role in osmotic protection, signal transduction and carbon storage. Under the circumstance of global climate change the abiotic stress would restrict the development of C. paliurus plantation, whereas there is few knowledge on the regulatory mechanisms of sugar metabolism under drought stress in C. paliurus. To investigate the drought response of C. paliurus at molecular level, we conducted an integrated analysis of transcriptomic and metabolomic of C. paliurus at three PEG-induced drought stress levels (0%: control; 15%: moderate drought; 25%: severe drought) in short term. Both moderate and severe drought treatments activated the chemical defense with lowering relative water content, and enhancing the contents of soluble protein, proline and malondialdehyde in the leaves. Meanwhile, alterations in the expression of differentially expressed genes and carbohydrate metabolism profiles were observed among the treatments. Weighted gene co-expression network analysis (WGCNA) showed 3 key modules, 8 structural genes (such as genes encoding beta-fructofuranosidase (INV), sucrose synthase (SUS), raffinose synthase (RS)) and 14 regulatory transcription factors were closely linked to sugar metabolism. Our results provided the foundation to understand the response mechanism of sugar metabolism in C. paliurus under drought stress, and would drive progress in breeding of drought-tolerant varieties and plantation development of the species.

Weighted gene co-expression network analysis of nitrogen (N)-responsive genes and the putative role of G-quadruplexes in N use efficiency (NUE) in rice

Rice is an important target to improve crop nitrogen (N) use efficiency (NUE), and the identification and shortlisting of the candidate genes are still in progress. We analyzed data from 16 published N-responsive transcriptomes/microarrays to identify, eight datasets that contained the maximum number of 3020 common genes, referred to as N-responsive genes. These include different classes of transcription factors, transporters, miRNA targets, kinases and events of post-translational modifications. A Weighted gene co-expression network analysis (WGCNA) with all the 3020 N-responsive genes revealed 15 co-expression modules and their annotated biological roles. Protein-protein interaction network analysis of the main module revealed the hub genes and their functional annotation revealed their involvement in the ubiquitin process. Further, the occurrences of G-quadruplex sequences were examined, which are known to play important roles in epigenetic regulation but are hitherto unknown in N-response/NUE. Out of the 3020 N-responsive genes studied, 2298 contained G-quadruplex sequences. We compared these N-responsive genes containing G-quadruplex sequences with the 3601 genes we previously identified as NUE-related (for being both N-responsive and yield-associated). This analysis revealed 389 (17%) NUE-related genes containing G-quadruplex sequences. These genes may be involved in the epigenetic regulation of NUE, while the rest of the 83% (1811) genes may regulate NUE through genetic mechanisms and/or other epigenetic means besides G-quadruplexes. A few potentially important genes/processes identified as associated with NUE were experimentally validated in a pair of rice genotypes contrasting for NUE. The results from the WGCNA and G4 sequence analysis of N-responsive genes helped identify and shortlist six genes as candidates to improve NUE. Further, the hitherto unavailable segregation of genetic and epigenetic gene targets could aid in informed interventions through genetic and epigenetic means of crop improvement.

Floral Development Stage-Specific Transcriptomic Analysis Reveals the Formation Mechanism of Different Shapes of Ray Florets in Chrysanthemum

The formation mechanism of different ray floret shapes of chrysanthemum (Chrysanthemum × morifolium) remains elusive due to its complex genetic background. C. vestitum, with the basic ray floret shapes of the flat, spoon, and tubular types, is considered a model material for studying ray floret morphogenesis. In this study, the flat and tubular type lines of C. vestitum at specific stages were used to investigate the key genes that regulate morphological differences in ray florets. We found that the expression levels of genes related to auxin synthesis, transport, and response were generally higher in the tubular type than in the flat type. CvARF3 was highly expressed in the flat type, while CvARF5 and CvARF6 were highly expressed in the tubular type. Additionally, the transcription levels of Class B and E genes closely related to petal development, including CvPI, CvAP3, Cvdefh21, CvSEP3, and CvCDM77, were expressed at higher levels in the tubular type than the flat type. Based on the results, it is proposed that auxin plays a key role in the development of ray florets, and auxin-related genes, especially CvARFs, may be key genes to control the morphological difference of ray florets. Simultaneously, MADS-box genes are involved in the co-regulation of ray floret morphogenesis. The results provide novel insights into the molecular mechanism of different petal type formation and lay a theoretical foundation for the directional breeding of petal type in chrysanthemums.

Identification of Key Modules and Candidate Genes for Powdery Mildew Resistance of Wheat-Agropyron cristatum Translocation Line WAT-2020-17-6 by WGCNA

As one of the serious diseases of wheat, powdery mildew (Blumeria graminis f. sp. tritici) is a long-term threat to wheat production. Therefore, it is of great significance to explore new powdery mildew-resistant genes for breeding. The wild relative species of wheat provide gene resources for resistance to powdery mildew breeding. Agropyron cristatum (2n = 4x = 28, genomes PPPP) is an important wild relative of wheat, carrying excellent genes for high yield, disease resistance, and stress resistance, which can be used for wheat improvement. To understand the molecular mechanism of powdery mildew resistance in the wheat-A. cristatum translocation line WAT2020-17-6, transcriptome sequencing was performed, and the resistance genes were analyzed by weighted gene co-expression network analysis (WGCNA). In the results, 42,845 differentially expressed genes were identified and divided into 18 modules, of which six modules were highly correlated with powdery mildew resistance. Gene ontology (GO) enrichment analysis showed that the six interested modules related to powdery mildew resistance were significantly enriched in N-methyltransferase activity, autophagy, mRNA splicing via spliceosome, chloroplast envelope, and AMP binding. The candidate hub genes of the interested modules were further identified, and their regulatory relationships were analyzed based on co-expression data. The temporal expression pattern of the 12 hub genes was verified within 96 h after powdery mildew inoculation by RT-PCR assay. In this study, we preliminarily explained the resistance mechanism of the wheat-A. cristatum translocation lines and obtained the hub candidate genes, which laid a foundation in the exploration of resistance genes in A. cristatum for powdery mildew-resistant breeding in wheat.

Transcriptomic analysis reveals the formation mechanism of anemone-type flower in chrysanthemum

BACKGROUND

The ray and disc florets on the chrysanthemum capitulum are morphologically diverse and have remarkably abundant variant types, resulting in a rich variety of flower types. An anemone shape with pigmented and elongated disk florets is an important trait in flower shape breeding of chrysanthemums. The regulatory mechanism of their anemone-type disc floret formation was not clear, thus limiting the directional breeding of chrysanthemum flower types. In this study, we used morphological observation, transcriptomic analysis, and gene expression to investigate the morphogenetic processes and regulatory mechanisms of anemone-type chrysanthemum.

RESULT

Scanning electron microscopy (SEM) observation showed that morphological differences between non-anemone-type disc florets and anemone-type disc florets occurred mainly during the petal elongation period. The anemone-type disc florets elongated rapidly in the later stages of development. Longitudinal paraffin section analysis revealed that the anemone-type disc florets were formed by a great number of cells in the middle layer of the petals with vigorous division. We investigated the differentially expressed genes (DEGs) using ray and disc florets of two chrysanthemum cultivars, 082 and 068, for RNA-Seq and their expression patterns of non-anemone-type and anemone-type disc florets. The result suggested that the CYCLOIDEA2 (CYC2s), MADS-box genes, and phytohormone signal-related genes appeared significantly different in both types of disc florets and might have important effects on the formation of anemone-type disc florets. In addition, it is noteworthy that the auxin and jasmonate signaling pathways might play a vital role in developing anemone-type disc florets.

CONCLUSIONS

Based on our findings, we propose a regulatory network for forming non-anemone-type and anemone-type disc florets. The results of this study lead the way to further clarify the mechanism of the anemone-type chrysanthemum formation and lay the foundation for the directive breeding of chrysanthemum petal types.

High-quality Japanese flounder genome aids in identifying stress-related genes using gene coexpression network

The Japanese flounder is one of the most economically important marine flatfish. However, due to the increased frequency of extreme weather events and high-density industrial farming, an increasing number of environmental stresses have become severe threats to the healthy development of the Japanese flounder culture industry. Herein, we produced a high-quality chromosome-scale Japanese flounder genome using PacBio Circular Consensus Sequencing technologies. The assembled Japanese flounder genome spanned 588.22 Mb with a contig N50 size of 24.35 Mb. In total, 105.89 Mb of repetitive sequences and 22,565 protein-coding genes were identified by genome annotation. In addition, 67 candidate genes responding to distinct stresses were identified by gene coexpression network analysis based on 16 published stress-related RNA-seq datasets encompassing 198 samples. A high-quality chromosome-scale Japanese flounder genome and candidate stress-related gene set will not only serve as key resources for genomics studies and further research on the underlying stress responsive molecular mechanisms in Japanese flounder but will also advance the progress of genetic improvement and comprehensive stress-resistant molecular breeding of Japanese flounder.

Measurement(s)	genome assembly
Technology Type(s)	PacBio RS II

Putative Transcription Factor Genes Associated with Regulation of Carotenoid Biosynthesis in Chili Pepper Fruits Revealed by RNA-Seq Coexpression Analysis

During the ripening process, the pericarp of chili pepper (Capsicum spp.) fruits accumulates large amounts of carotenoids. Although the carotenoid biosynthesis pathway in the Capsicum genus has been widely studied from different perspectives, the transcriptional regulation of genes encoding carotenoid biosynthetic enzymes has not been elucidated in this fruit. We analyzed RNA-Seq transcriptomic data from the fruits of 12 accessions of Capsicum annuum during the growth, development, and ripening processes using the R package named Salsa. We performed coexpression analyses between the standardized expression of genes encoding carotenoid biosynthetic enzymes (target genes (TGs)) and the genes of all expressed transcription factors (TFs). Additionally, we analyzed the promoter region of each biosynthetic gene to identify putative binding sequences for each selected TF candidate. We selected 83 TFs as putative regulators of the carotenogenic structural genes. From them, putative binding sites in the promoters of the carotenoid-biosynthesis-related structural genes were found for only 54 TFs. These results could guide the search for transcription factors involved in the regulation of the carotenogenic pathway in chili pepper fruits and might facilitate the collection of corresponding experimental evidence to corroborate their participation in the regulation of this biosynthetic pathway in Capsicum spp.

Gene co-expression network analysis reveals the positive impact of endocytosis and mitochondria-related genes over nitrogen metabolism in Saccharomyces cerevisiae

Nitrogen metabolism is essential for most cellular activities. Therefore, a deep understanding of its regulatory mechanisms is necessary for the efficient utilization of nitrogen sources for Saccharomyces cerevisiae. In this study, a gene co-expression network was constructed for S. cerevisiae S288C with different nitrogen sources. From this, a key gene co-expression module related to nitrogen source preference utilization was obtained, and 10 hub genes centrally located in the co-expression network were identified. Functional studies verified that the endocytosis-related genes CAP1 and END3 significantly increased the utilization of multiple non-preferred amino acids and reduced the accumulation of the harmful nitrogen metabolite precursor urea by regulating amino acid transporters and TOR pathway. The mitochondria-related gene ATP12, MRPL22, MRP1 and NAM9 significantly increased the utilization of multiple non-preferred amino acids and reduced accumulation of the urea by coordinately regulating nitrogen catabolism repression, Ssy1p-Ptr3p-Ssy5p signaling sensor system, amino acid transporters, TOR pathway and urea metabolism-related pathways. Furthermore, these data revealed the potential positive effects of endocytosis and mitochondrial ribosomes protein translation on nitrogen source preference. This study provides new analytical perspectives for complex regulatory networks involving nitrogen metabolism in S. cerevisiae.

Proteomics and Co-expression Network Analysis Reveal the Importance of Hub Proteins and Metabolic Pathways in Nicotine Synthesis and Accumulation in Tobacco (Nicotiana tabacum L.)

Nicotine is a unique alkaloid present in tobacco that is widely used in cigarettes and in the agricultural, chemical, and pharmaceutical industries. However, the research on nicotine is mostly limited to its synthesis pathways, and only a few studies have explored the effects of other metabolic pathways on nicotine precursors. Regulating the nicotine content in tobacco can greatly promoting the application of nicotine in other fields. In this study, we performed global data-independent acquisition proteomics analysis of four tobacco varieties. Of the four varieties, one had high nicotine content and three had a low nicotine content. A total of 31,259 distinct peptides and 6,018 proteins across two samples were identified. A total of 45 differentially expressed proteins (DEPs) co-existed in the three comparison groups and were mainly involved in the transport and metallic processes of the substances. Most DEPs were enriched in the biosynthesis of secondary metals, glutathione metabolism, carbon metabolism, and glycolysis/gluconeogenesis. In addition, the weighted gene co-expression network analysis identified an expression module closely related to the nicotine content (Brown, r = 0.74, P = 0.006). Gene Ontology annotation and Kyoto Encyclopaedia of Genes and Genomes enrichment analysis showed that the module proteins were mainly involved in the synthesis and metabolism of nicotine precursors such as arginine, ornithine aspartate, proline, and glutathione. The increased levels of these precursors lead to the synthesis and accumulation of nicotine in plants. More importantly, these proteins regulate nicotine synthesis by affecting the formation of putrescine, which is the core intermediate product in nicotine anabolism. Our results provide a reference for tobacco variety selection with a suitable nicotine content and regulation of the nicotine content. Additionally, the results highlight the importance of other precursor metabolism in nicotine synthesis.

Biosynthetic pathways of triterpenoids and strategies to improve their Biosynthetic Efficiency

"Triterpenoids" can be considered natural products derived from the cyclization of squalene, yielding 3-deoxytriterpenes (hydrocarbons) or 3-hydroxytriterpenes. Triterpenoids are metabolites of these two classes of triterpenes, produced by the functionalization of their carbon skeleton. They can be categorized into different groups based on their structural formula/design. Triterpenoids are an important group of compounds that are widely used in the fields of pharmacology, food, and industrial biotechnology. However, inadequate synthetic methods and insufficient knowledge of the biosynthesis of triterpenoids, such as their structure, enzymatic activity, and the methods used to produce pure and active triterpenoids, are key problems that limit the production of these active metabolites. Here, we summarize the derivatives, pharmaceutical properties, and biosynthetic pathways of triterpenoids and review the enzymes involved in their biosynthetic pathway. Furthermore, we concluded the screening methods, identified the genes involved in the pathways, and highlighted the appropriate strategies used to enhance their biosynthetic production to facilitate the commercial process of triterpenoids through the synthetic biology method.

Functional analysis of the ScAG and ScAGL11 MADS-box transcription factors for anthocyanin biosynthesis and bicolour pattern formation in Senecio cruentus ray florets

Cineraria (Senecio cruentus) is an ornamental plant with pure colour and bicolour cultivars, widely used for landscaping. Anthocyanin biosynthesis influences coloration patterns in cineraria. However, how anthocyanins accumulate and distribute in cineraria is poorly understood. This study investigated the molecular mechanisms underlying anthocyanin biosynthesis and bicolour formation in cineraria using pure colour and bicolour cultivars. Transcriptome and gene expression analysis showed that five genes, ScCHS2, ScF3H1, ScDFR3, ScANS, and ScbHLH17, were inhibited in the white cultivar and colourless regions of bicolour cultivars. In contrast, two MADS-box genes, ScAG and ScAGL11, showed significantly higher expression in the colourless regions of bicolour cultivars. ScAG and ScAGL11 were localized in the nucleus and co-expressed with the bicolour trait. Further functional analysis verified that ScAG inhibits anthocyanin accumulation in tobacco (Nicotiana tabacum). However, virus-induced gene silencing (VIGS) experiments showed that silencing of ScAG and ScAGL11 increases anthocyanin content in cineraria leaves. Similar results were observed when ScAG and ScAGL11 were silenced in the cineraria capitulum, accompanied by the smaller size of the colourless region, specifically in the ScAG/ScAGL11-silenced plants. The expression of ScCHS2, ScDFR3, and ScF3H1 increased in silenced cineraria leaves and capitulum. Furthermore, yeast two-hybrid and bimolecular fluorescence complementation experiments demonstrated that ScAG interacts with ScAGL11. Moreover, ScAG directly inhibited the transcription of ScF3H1 while ScAGL11 inhibited ScDFR3 expression by binding to their promoters separately. The findings reported herein indicate that ScAG and ScAGL11 negatively regulate anthocyanin biosynthesis in cineraria ray florets, and their differential expression in ray florets influences the bicolour pattern appearance.

Identification of key gene networks controlling anthocyanin biosynthesis in peach flower

Flavonoids, particularly anthocyanin is the main pigment that determined the red color of peach flowers, and help the plant to attract pollinators, protect the reproductive organs of flower from photo-oxidative effects of light and various non-communicable diseases. Through weightage gene coexpression network analysis (WGCNA) we identified a network of 15 hub genes that co-expressed throughout peach flower development including 5 genes coded for the key enzymes (CHI, F3'H, DFR, LAR and UFGT) of flavonoid biosynthetic pathway and 1 gene Prupe.1G111700 identified as R2R3 family transcription factor MYB108. Over expression of PpMYB108 significantly increased anthocyanin biosynthesis in Tobacco flowers. Moreover, the expression correlation between PpMYB108 and PpDFR, suggests that PpMYB108 play the role of transcriptional activator for PpDFR. This was further supported by a 6 bp insertion of MYB biding site in the core promoter region of PpDFR in red flower. The positive interaction of PpMYB108 with PpDFR promoter from red flower was confirmed in yeast one hybrid assay. These findings may be helpful in peach breeding programs as well as in identifying anthocyanin related genes in other species.

The transcription factor complex CmAP3-CmPI-CmUIF1 modulates carotenoid metabolism by directly regulating carotenogenic gene CmCCD4a-2 in chrysanthemum

Carotenoids are one of the most important pigments for the coloring in many plants, fruits and flowers. Recently, significant progress has been made in carotenoid metabolism. However, the specific understanding on transcriptional regulation controlling the expression of carotenoid metabolic genes remains extremely limited. Anemone-type chrysanthemum, as a special group of chrysanthemum cultivars, contain elongated disc florets in capitulum, which usually appear in different colors compared with the ray florets since accumulating distinct content of carotenoids. In this study, the carotenoid composition and content of the ray and disc florets of an anemone-type chrysanthemum cultivar 'Dong Li Fen Gui' were analyzed by high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) and the key structural gene CmCCD4a-2, of which differential expression resulted in the distinct content of carotenoids accumulated in these two types of florets, was identified. Then the promoter sequence of CmCCD4a-2 was used as bait to screen a chrysanthemum flower cDNA library and two transcription factors, CmAP3 and CmUIF1 were identified. Y2H, BiFC and Y3H experiments demonstrated that these two TFs were connected by CmPI to form CmAP3-CmPI-CmUIF1 TF complex. This TF complex regulated carotenoid metabolism through activating the expression of CmCCD4a-2 directly. Furthermore, a large number of target genes regulated directly by the CmAP3-CmPI-CmUIF1 TF complex, including carotenoid biosynthetic genes, flavonoid biosynthetic genes and flower development-related genes, were identified by DNA-affinity purification sequencing (DAP-seq), which indicated that the CmAP3-CmPI-CmUIF1 TF complex might participate in multiple processes. These findings expand our knowledge for the transcriptional regulation of carotenoid metabolism in plants and will be helpful to manipulating carotenoid accumulation in chrysanthemum.

Integrative Analysis of Metabolome and Transcriptome Reveals the Mechanism of Color Formation in Liriope spicata Fruit

Liriope spicata is an important ornamental ground cover plant, with a fruit color that turns from green to black during the development and ripening stages. However, the material basis and regulatory mechanism of the color variation remains unclear. In this study, a total of 31 anthocyanins and 2 flavonols were identified from the skin of L. spicata fruit via integrative analysis on the metabolome and transcriptome of three developmental stages. The pigments of black/mature fruits are composed of five common anthocyanin compounds, of which Peonidin 3–O–rutinoside and Delphinidin 3–O–glucoside are the most differential metabolites for color conversion. Using dual-omics joint analysis, the mechanism of color formation was obtained as follows. The expression of structural genes including 4CL, F3H, F3′H, F3′5′H and UFGT were activated due to the upregulation of transcription factor genes MYB and bHLH. As a result, a large amount of precursor substances for the synthesis of flavonoids accumulated. After glycosylation, stable pigments were generated which promoted the accumulation of anthocyanins and the formation of black skin.

Transcriptomic and Chemical Analyses Reveal the Hub Regulators of Flower Color Variation from Camellia japonica Bud Sport

Camellia japonica is a woody ornamental plant with multiple flower color variations caused by bud sport; however, the molecular mechanism remains unclear. Here, chemical and transcriptomic analyses of C. japonica were performed with white, pink, red, and dark red flowers caused by bud sport. Seven anthocyanins were detected in these samples, except in C. japonica ’YuDan’ (white petals). The total anthocyanin content of C. japonica ’JinBiHuiHuang’ was the highest, and cyanidin 3-O-β-glucoside (Cy3G) was the main anthocyanin affecting the redness of petals. Furthermore, the ratio of Cy3G and cyanidin-3-O-(6-O-(E)-p-coumaroyl)-B-glucoside) was significantly correlated with the red petal phenotype. In total, 5673 genes were identified as differentially expressed genes (DEGs). The potential co-expression modules related to anthocyanin accumulation were established, which featured transcription factors, anthocyanin biosynthesis, and plant hormone signal transduction. Thirteen structural genes in the anthocyanin biosynthetic pathway were identified as DEGs, most of them were upregulated with deepening of flower redness. An integrated promoter and cluster analysis suggested that CjMYB62, CjMYB52, and CjGATA may play important roles in anthocyanin accumulation. These results provide insight and candidate genes for the transcriptional mechanism responsible for the bud sport phenotype.

From gene to biomolecular networks: a review of evidences for understanding complex biological function in plants

Although at the infancy stage, biomolecular network biology is a comprehensive approach to understand complex biological function in plants. Recent advancements in the accumulation of multi-omics data coupled with computational approach have accelerated our current understanding of the complexities of gene function at the system level. Biomolecular networks such as protein-protein interaction, co-expression and gene regulatory networks have extensively been used to decipher the intricacies of transcriptional reprogramming of hundreds of genes and their regulatory interaction in response to various environmental perturbations mainly in the model plant Arabidopsis. This review describes recent applications of network-based approaches to understand the biological functions in plants and focuses on the challenges and opportunities to harness the full potential of the approach.

Optimizing weighted gene co-expression network analysis with a multi-threaded calculation of the topological overlap matrix

Biomolecular networks are often assumed to be scale-free hierarchical networks. The weighted gene co-expression network analysis (WGCNA) treats gene co-expression networks as undirected scale-free hierarchical weighted networks. The WGCNA R software package uses an Adjacency Matrix to store a network, next calculates the topological overlap matrix (TOM), and then identifies the modules (sub-networks), where each module is assumed to be associated with a certain biological function. The most time-consuming step of WGCNA is to calculate TOM from the Adjacency Matrix in a single thread. In this paper, the single-threaded algorithm of the TOM has been changed into a multi-threaded algorithm (the parameters are the default values of WGCNA). In the multi-threaded algorithm, Rcpp was used to make R call a C++ function, and then C++ used OpenMP to start multiple threads to calculate TOM from the Adjacency Matrix. On shared-memory MultiProcessor systems, the calculation time decreases as the number of CPU cores increases. The algorithm of this paper can promote the application of WGCNA on large data sets, and help other research fields to identify sub-networks in undirected scale-free hierarchical weighted networks. The source codes and usage are available at https://github.com/do-somethings-haha/multi-threaded_calculate_unsigned_TOM_from_unsigned_or_signed_Adjacency_Matrix_of_WGCNA.

Comparison of chrysanthemum flowers grown under hydroponic and soil-based systems: yield and transcriptome analysis

BACKGROUND

Flowers of Chrysanthemum × morifolium Ramat. are used as tea in traditional Chinese cuisine. However, with increasing population and urbanization, water and land availability have become limiting for chrysanthemum tea production. Hydroponic culture enables effective, rapid nutrient exchange, while requiring no soil and less water than soil cultivation. Hydroponic culture can reduce pesticide residues in food and improve the quantity or size of fruits, flowers, and leaves, and the levels of active compounds important for nutrition and health. To date, studies to improve the yield and active compounds of chrysanthemum have focused on soil culture. Moreover, the molecular effects of hydroponic and soil culture on chrysanthemum tea development remain understudied.

RESULTS

Here, we studied the effects of soil and hydroponic culture on yield and total flavonoid and chlorogenic acid contents in chrysanthemum flowers (C. morifolium 'wuyuanhuang'). Yield and the total flavonoids and chlorogenic acid contents of chrysanthemum flowers were higher in the hydroponic culture system than in the soil system. Transcriptome profiling using RNA-seq revealed 3858 differentially expressed genes (DEGs) between chrysanthemum flowers grown in soil and hydroponic conditions. Gene Ontology (GO) enrichment annotation revealed that these differentially transcribed genes are mainly involved in "cytoplasmic part", "biosynthetic process", "organic substance biosynthetic process", "cell wall organization or biogenesis" and other processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed enrichment in "metabolic pathways", "biosynthesis of secondary metabolites", "ribosome", "carbon metabolism", "plant hormone signal transduction" and other metabolic processes. In functional annotations, pathways related to yield and formation of the main active compounds included phytohormone signaling, secondary metabolism, and cell wall metabolism. Enrichment analysis of transcription factors also showed that under the hydroponic system, bHLH, MYB, NAC, and ERF protein families were involved in metabolic pathways, biosynthesis of secondary metabolites, and plant hormone signal transduction.

CONCLUSIONS

Hydroponic culture is a simple and effective way to cultivate chrysanthemum for tea production. A transcriptome analysis of chrysanthemum flowers grown in soil and hydroponic conditions. The large number of DEGs identified confirmed the difference of the regulatory machinery under two culture system.

Comparison of the transcriptomic responses of two Chrysanthemum morifolium cultivars to low light

BACKGROUND

Low light is a primary regulator of chrysanthemum growth. Our aim was to analyse the different transcriptomic responses of two Chrysanthemum morifolium cultivars to low light.

METHODS AND RESULTS

We conducted a transcriptomic analysis of leaf samples from the 'Nannonggongfen' and 'Nannongxuefeng' chrysanthemum cultivars following a 5-day exposure to optimal light (70%, control [CK]) or low-light (20%, LL) conditions. Gene Ontology (GO) classification of upregulated genes revealed these genes to be associated with 11 cellular components, 9 molecular functions, and 15 biological processes, with the majority being localized to the chloroplast, highlighting the role of chloroplast proteins as regulators of shading tolerance. Downregulated genes were associated with 11 cellular components, 8 molecular functions, and 16 biological processes. Heat map analyses suggested that basic helix-loop-helix domain genes and elongation factors were markedly downregulated in 'Nannongxuefeng' leaves, consistent with the maintenance of normal stem length, whereas no comparable changes were observed in 'Nanonggongfen' leaves. Subsequent qPCR analyses revealed that phytochrome-interacting factors and dormancy-associated genes were significantly upregulated under LL conditions relative to CK conditions, while succinate dehydrogenase 1, elongated hypocotyls 5, and auxin-responsive gene of were significantly downregulated under LL conditions.

CONCLUSIONS

These findings suggest that LL plants were significantly lower than those of the CK plants. Low-light tolerant chrysanthemum cultivars may maintain reduced indole-3-acetic acid (IAA) and elongation factor expression as a means of preventing the onset of shade-avoidance symptoms.

Transcriptome-Based WGCNA Analysis Reveals Regulated Metabolite Fluxes between Floral Color and Scent in Narcissus tazetta Flower

A link between the scent and color of Narcissus tazetta flowers can be anticipated due to their biochemical origin, as well as their similar biological role. Despite the obvious aesthetic and ecological significance of these colorful and fragrant components of the flowers and the molecular profiles of their pigments, fragrant formation has addressed in some cases. However, the regulatory mechanism of the correlation of fragrant components and color patterns is less clear. We simultaneously used one way to address how floral color and fragrant formation in different tissues are generated during the development of an individual plant by transcriptome-based weighted gene co-expression network analysis (WGCNA). A spatiotemporal pattern variation of flavonols/carotenoids/chlorophyll pigmentation and benzenoid/phenylpropanoid/ monoterpene fragrant components between the tepal and corona in the flower tissues of Narcissus tazetta, was exhibited. Several candidate transcription factors: MYB12, MYB1, AP2-ERF, bZIP, NAC, MYB, C2C2, C2H2 and GRAS are shown to be associated with metabolite flux, the phenylpropanoid pathway to the production of flavonols/anthocyanin, as well as related to one branch of the phenylpropanoid pathway to the benzenoid/phenylpropanoid component in the tepal and the metabolite flux between the monoterpene and carotenoids biosynthesis pathway in coronas. It indicates that potential competition exists between floral pigment and floral fragrance during Narcissus tazetta individual plant development and evolutionary development.

Comparative Transcriptome and Weighted Gene Co-expression Network Analysis Identify Key Transcription Factors of Rosa chinensis 'Old Blush' After Exposure to a Gradual Drought Stress Followed by Recovery

Rose is one of the most fundamental ornamental crops, but its yield and quality are highly limited by drought. The key transcription factors (TFs) and co-expression networks during rose’s response to drought stress and recovery after drought stress are still limited. In this study, the transcriptomes of leaves of 2-year-old cutting seedlings of Rosa chinensis ‘Old Blush’ from three continuous droughted stages (30, 60, 90 days after full watering) and rewatering were analyzed using RNA sequencing. Weighted gene co-expression network analysis (WGCNA) was used to construct a co-expression network, which was associated with the physiological traits of drought response to discovering the hub TFs involved in drought response. More than 45 million high-quality clean reads were generated from the sample and used for comparison with the rose reference genome. A total of 46433 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that drought stress caused significant changes in signal transduction, plant hormones including ABA, auxin, brassinosteroid (BR), cytokinin, ethylene (ET), jasmonic acid (JA) and salicylic acid (SA), primary and secondary metabolism, and a certain degree of recovery after rewatering. Gene co-expression analysis identified 18 modules, in which four modules showed a high degree of correlation with physiological traits. In addition, 42 TFs including members of NACs, WRKYs, MYBs, AP2/ERFs, ARFs, and bHLHs with high connectivity in navajowhite1 and blue modules were screened. This study provides the transcriptome sequencing report of R. chinensis ‘Old Blush’ during drought stress and rewatering process. The study also identifies the response of candidate TFs to drought stress, providing guidelines for improving the drought tolerance of the rose through molecular breeding in the future.

Elucidation of the Regulatory Network of Flavonoid Biosynthesis by Profiling the Metabolome and Transcriptome in Tartary Buckwheat

The characteristics of flavonoid metabolism in different Tartary buckwheat (TB) tissues and the related gene regulation network are still unclear at present. One hundred forty-seven flavonoids were identified from six TB tissues using the ultra performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS) method. The roadmap of the rutin synthesis pathway was revealed. Through transcriptomic analysis it was revealed that the differentially expressed genes (DEGs) are mainly enriched in the "Phenylpropanoid biosynthesis" pathway. Fifty-two DEGs involved in the "flavonol synthesis" pathway were identified. The weighted gene correlation network analysis revealed four co-expression network modules correlated with six flavonol metabolites. Eventually, 74 genes revealed from MEblue and MElightsteelblue modules were potentially related to flavonol synthesis. Of them, 7 MYB transcript factors had been verified to regulate flavonoid synthesis. Furthermore, overexpressed FtMYB31 enhanced the rutin content in vivo. The present findings provide a dynamic flavonoid metabolism profile and co-expression network related to rutin synthesis and are thus valuable in understanding the molecular mechanisms of rutin synthesis in TB.

Identification of Chlorophyll Metabolism- and Photosynthesis-Related Genes Regulating Green Flower Color in Chrysanthemum by Integrative Transcriptome and Weighted Correlation Network Analyses

Green chrysanthemums are difficult to breed but have high commercial value. The molecular basis for the green petal color in chrysanthemum is not fully understood. This was investigated in the present study by RNA sequencing analysis of white and green ray florets collected at three stages of flower development from the F₁ progeny of the cross between Chrysanthemum × morifolium “Lüdingdang” with green-petaled flowers and Chrysanthemum vistitum with white-petaled flowers. The chlorophyll content was higher and chloroplast degradation was slower in green pools than in white pools at each developmental stage. Transcriptome analysis revealed that genes that were differentially expressed between the two pools were enriched in pathways related to chlorophyll metabolism and photosynthesis. We identified the transcription factor genes CmCOLa, CmCOLb, CmERF, and CmbHLH as regulators of the green flower color in chrysanthemum by differential expression analysis and weighted gene co-expression network analysis. These findings can guide future efforts to improve the color palette of chrysanthemum flowers through genetic engineering.

Transcriptomic and chemical analyses to identify candidate genes involved in color variation of sainfoin flowers

BACKGROUND

Sainfoin (Onobrychis viciifolia Scop) is not only a high-quality legume forage, but also a nectar-producing plant. Therefore, the flower color of sainfoin is an important agronomic trait, but the factors affecting its flower phenotype are still unclear. To gain insights into the regulatory networks associated with metabolic pathways of coloration compounds (flavonoids or anthocyanins) and identify the key genes, we conducted a comprehensive analysis of the phenotype, metabolome and transcriptome of WF and AF of sainfoin.

RESULTS

Delphinidin, petunidin and malvidin derivatives were the main anthocyanin compounds in the AF of sainfoin. These substances were not detected in the WF of sainfoin. The transcriptomes of WF and AF in sainfoin at the S1 and S3 stages were obtained using the Illumina HiSeq4000 platform. Overall, 10,166 (4273 upregulated and 5893 downregulated) and 15,334 (8174 upregulated and 7160 downregulated) DEGs were identified in flowers at S1 and S3 stages, respectively (WF-VS-AF). KEGG pathway annotations showed that 6396 unigenes were annotated to 120 pathways and contained 866 DEGs at S1 stages, and 6396 unigenes were annotated to 131 pathways and included 1546 DEGs at the S3 stage. Nine DEGs belonging to the "flavonoid biosynthesis"and "phenylpropanoid biosynthesis" pathways involved in flower color formation were identified and verified by RT-qPCR analyses. Among these DEGs, 4CL3, FLS, ANS, CHS, DFR and CHI2 exhibited downregulated expression, and F3H exhibited upregulated expression in the WF compared to the AF, resulting in a decrease in anthocyanin synthesis and the formation of WF in sainfoin.

CONCLUSIONS

This study is the first to use transcriptome technology to study the mechanism of white flower formation in sainfoin. Our transcriptome data will be a great enrichment of the genetic information for sainfoin. In addition, the data presented herein will provide valuable molecular information for genetic breeding and provide insight into the future study of flower color polymorphisms in sainfoin.

Isolation and Functional Analysis of Genes Involved in Polyacylated Anthocyanin Biosynthesis in Blue Senecio cruentus

Polyacylated anthocyanins with multiple glycosyl and aromatic acyl groups tend to make flowers display bright and stable blue colours. However, there are few studies on the isolation and functional characterization of genes involved in the polyacylated anthocyanin biosynthesis mechanism, which limits the molecular breeding of truly blue flowers. Senecio cruentus is an important potted ornamental plant, and its blue flowers contain 3',7-polyacylated delphinidin-type anthocyanins that are not reported in any other plants, suggesting that it harbours abundant gene resources for the molecular breeding of blue flowers. In this study, using high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) analysis of blue, carmine and white colours of cineraria cultivars "Venezia" (named VeB, VeC, and VeW, respectively), we found that 3',7-polyacylated anthocyanin, cinerarin, was the main pigment component that determined the blue colour of ray florets of cineraria. Based on the transcriptome sequencing and differential gene expression (DEG) analysis combined with RT- and qRT-PCR, we found two genes encoding uridine diphosphate glycosyltransferase, named ScUGT1 and ScUGT4; two genes encoding acyl-glucoside-dependent glucosyltransferases which belong to glycoside hydrolase family 1 (GH1), named ScAGGT11 and ScAGGT12; one gene encoding serine carboxypeptidase-like acyltransferase ScSCPL2; and two MYB transcriptional factor genes ScMYB2 and ScMYB4, that were specifically highly expressed in the ray florets of VeB, which indicated that these genes may be involved in cinerarin biosynthesis. The function of ScSCPL2 was analysed by virus-induced gene silencing (VIGS) in cineraria leaves combined with HPLC-MS/MS. ScSCPL2 mainly participated in the 3' and 7-position acylation of cinerarin. These results will provide new insight into the molecular basis of the polyacylated anthocyanin biosynthesis mechanism in higher plants and are of great significance for blue flower molecular breeding of ornamental plants.

Gene co-expression network analysis to identify critical modules and candidate genes of drought-resistance in wheat

AIM

To establish a gene co-expression network for identifying principal modules and hub genes that are associated with drought resistance mechanisms, analyzing their mechanisms, and exploring candidate genes.

METHODS AND FINDINGS

42 data sets including PRJNA380841 and PRJNA369686 were used to construct the co-expression network through weighted gene co-expression network analysis (WGCNA). A total of 1,896,897,901 (284.30 Gb) clean reads and 35,021 differentially expressed genes (DEGs) were obtained from 42 samples. Functional enrichment analysis indicated that photosynthesis, DNA replication, glycolysis/gluconeogenesis, starch and sucrose metabolism, arginine and proline metabolism, and cell cycle were significantly influenced by drought stress. Furthermore, the DEGs with similar expression patterns, detected by K-means clustering, were grouped into 29 clusters. Genes involved in the modules, such as dark turquoise, yellow, and brown, were found to be appreciably linked with drought resistance. Twelve central, greatly correlated genes in stage-specific modules were subsequently confirmed and validated at the transcription levels, including TraesCS7D01G417600.1 (PP2C), TraesCS5B01G565300.1 (ERF), TraesCS4A01G068200.1 (HSP), TraesCS2D01G033200.1 (HSP90), TraesCS6B01G425300.1 (RBD), TraesCS7A01G499200.1 (P450), TraesCS4A01G118400.1 (MYB), TraesCS2B01G415500.1 (STK), TraesCS1A01G129300.1 (MYB), TraesCS2D01G326900.1 (ALDH), TraesCS3D01G227400.1 (WRKY), and TraesCS3B01G144800.1 (GT).

CONCLUSIONS

Analyzing the response of wheat to drought stress during different growth stages, we have detected three modules and 12 hub genes that are associated with drought resistance mechanisms, and five of those genes are newly identified for drought resistance. The references provided by these modules will promote the understanding of the drought-resistance mechanism. In addition, the candidate genes can be used as a basis of transgenic or molecular marker-assisted selection for improving the drought resistance and increasing the yields of wheat.