Transcriptome Sequencing and Metabolome Analysis Reveal Genes Involved in Pigmentation of Green-Colored Cotton Fibers

Hub Gene Mining and Co-Expression Network Construction of Low-Temperature Response in Maize of Seedling by WGCNA

Weighted gene co-expression network analysis (WGCNA) is a research method in systematic biology. It is widely used to identify gene modules related to target traits in multi-sample transcriptome data. In order to further explore the molecular mechanism of maize response to low-temperature stress at the seedling stage, B144 (cold stress tolerant) and Q319 (cold stress sensitive) provided by the Maize Research Institute of Heilongjiang Academy of Agricultural Sciences were used as experimental materials, and both inbred lines were treated with 5 °C for 0 h, 12 h, and 24 h, with the untreated material as a control. Eighteen leaf samples were used for transcriptome sequencing, with three biological replicates. Based on the above transcriptome data, co-expression networks of weighted genes associated with low-temperature-tolerance traits were constructed by WGCNA. Twelve gene modules significantly related to low-temperature tolerance at the seedling stage were obtained, and a number of hub genes involved in low-temperature stress regulation pathways were discovered from the four modules with the highest correlation with target traits. These results provide clues for further study on the molecular genetic mechanisms of low-temperature tolerance in maize at the seedling stage.

Comparative Transcriptome Analysis Revealed Key Genes Regulating Gossypol Synthesis in Tetraploid Cultivated Cotton

Tetraploid cultivated cotton (Gossypium spp.) produces cottonseeds rich in protein and oil. Gossypol and related terpenoids, stored in the pigment glands of cottonseeds, are toxic to human beings and monogastric animals. However, a comprehensive understanding of the genetic basis of gossypol and gland formation is still lacking. We performed a comprehensive transcriptome analysis of four glanded versus two glandless tetraploid cultivars distributed in Gossypium hirsutum and Gossypium barbadense. A weighted gene co-expression network analysis (WGCNA) based on 431 common differentially expressed genes (DEGs) uncovered a candidate module that was strongly associated with the reduction in or disappearance of gossypol and pigment glands. Further, the co-expression network helped us to focus on 29 hub genes, which played key roles in the regulation of related genes in the candidate module. The present study contributes to our understanding of the genetic basis of gossypol and gland formation and serves as a rich potential source for breeding cotton cultivars with gossypol-rich plants and gossypol-free cottonseed, which is beneficial for improving food safety, environmental protection, and economic gains of tetraploid cultivated cotton.

Cloning and functional analysis of GhDFR1, a key gene of flavonoid synthesis pathway in naturally colored cotton

BACKGROUND

The naturally colored brown cotton fiber is the most widely used environmentally friendly textile material, which primarily contains proanthocyanidins and their derivatives. Many structural genes in the flavonoid synthesis pathway are known to improve the genetic resources of naturally colored cotton. Among them, DFR is a crucial late enzyme to synthesis both anthocyanins and proanthocyanidins in the plant flavonoid pathway.

METHODS

The protein sequences of GhDFRs were analyzed using bioinformatic tools. The expression levels of GhDFRs in various tissues and organs of upland cotton Zongxu1 (ZX1), were analyzed by quantitative real-time PCR, and the expression pattern of GhDFR1 during fiber development of white cotton and brown cotton was analyzed further. The function of GhDFR1 in NCC ZX1 was preliminarily analyzed by virus induced gene silencing (VIGS) technology.

RESULTS

Bioinformatic analysis revealed that GhDFRs sequences in upland cotton genome were extremely conserved. Furthermore, evolutionary tree analysis revealed that the functions of GhDFR1 and GhDFR2, and GhDFR3 and GhDFR4, presented different and shared some similarities. Our study showed GhDFR1 and GhDFR2 were specifically expressed in fibers, while GhDFR3 and GhDFR4 were specifically expressed in petals. GhDFR1 was exclusively expressed in brown cotton fiber at various stages of development and progressively increased with the growth of fiber, but the trend of expression in white cotton was quite the opposite. We silenced GhDFR1 expression in brown cotton fiber using VIGS technology, and observed the VIGS-interference plants. After reducing the expression level of GhDFR1, the period for significant GhDFR1 expression in the developing fibers changed, reducing the content of anthocyanins, and lightening the color of mature cotton fibers.

CONCLUSION

GhDFR1 was preferentially expressed in brown cotton during fiber development. The timing of GhDFR1 expression for flavonoid synthesis altered, resulting in anthocyanin contents reduced and the fiber color of the GhDFR1i lines lightened. These findings showed the role of GhDFR1 in fiber coloration of NCC and provided a new candidate for NCC genetic improvement.

Identification and Functional Analysis of the Promoter of a Leucoanthocyanidin Reductase Gene from Gossypium hirsutum

Leucoanthocyanidin reductase (LAR) is the critical enzyme in the synthesis pathway of proanthocyanidins, which are the primary pigments in brown cotton fibers. Our previous study has revealed significant differences in the expression levels of GhLAR1 between white and brown cotton fibers at 10 DPA. In this work, the expression pattern of the GhLAR1 gene was further studied, and the promoter of GhLAR1 (1780 bp) was isolated and characterized. Bioinformatic analysis indicated that GhLAR1 promoter contained many known light response elements and several defenses related to transcriptional factor-binding boxes, which may partially explain the response of the GhLAR1 to temperature, NaCl, and PEG treatments. Furthermore, GhLAR1 was preferentially and strongly expressed in fibers and flowers of cotton, and the expression levels in all tested tissues (especially fibers) of brown cotton were significantly higher than those in white cotton. Consistent with the expression analysis, the GhLAR1 promoter mainly drove GUS expression in epidermal trichomes and floral organs.

Comprehensive genomic identification and expression analysis 4CL gene family in apple

To clarify the structural characteristics, phylogeny, biological function and regulation of 4-coumarate-CoAligase (4CL) in anthocyanin synthesis, the 4CL gene family members in apples were identified and bioinformatic analysis was performed. qRT-PCR was used to analyze the expression levels of this gene family members in different apple varieties, and the role of the 4CL gene in apple anthocyanin synthesis was preliminaries clarified, which provided a certain theoretical basis for the regulatory network of apple anthocyanin synthesis. The results showed that a total of 69 members of the 4CL gene family were identified in the apple (Malus domestica Brokh.), encoding amino acids ranging from 97 to 2310 with theoretical isoelectric points ranging from 5.28 to 9.84. The 69 4CL family members were distributed on 17 chromosomes in the apple, among which chromosome 17 had the largest distribution (9 members), followed by chromosome 9 (7 members), chromosomes 16 and 14 (6 members each), and chromosomes 15 and 13 (5 members each). The subcellular localization prediction showed that apple 4CL gene family members were mainly expressed in cytoplasm, chloroplast, nucleus and cell membrane, with a small amount of expression in mitochondria, vacuoles, peroxisomes, cytoskeleton, golgi and cell matrix, but not in endoplasmic reticulum. The secondary structures are mainly α-helices and irregular coils. Microarray expression profile analysis showed that the expression levels of each member in apple were related to fruit variety and tissue structure, and the expression levels were mainly higher in fruit, flower and leaf. Real-time PCR analysis showed that the expression level of each member was directly proportional to the degree of fruit coloring and anthocyanin accumulation. The expression levels of Md4CL10 and Md4CL23 in 'Astar' (G4) apple fruit skin with the highest anthocyanin content were 516, 20 and 2 times higher than those in 'Chengji NO.1' (G1), 'Golden Delicious' (G2) and 'Ruixue' (G3), respectively.

Transcriptome and Metabolome Profiling Unveil Pigment Formation Variations in Brown Cotton Lines (Gossypium hirsutum L.)

Naturally brown colored cotton (NBCC) is becoming increasingly popular due to its natural properties of coloration. However, poor fiber quality and color fading are key issues that are hindering the cultivation of naturally colored cotton. In this study, based on transcriptome and metabolome of 18 days post-anthesis (DPA), we compared the variations of pigment formation in two brown cotton fibers (DCF and LCF), with white cotton fiber (WCF) belonging to a near-isogenic line. A transcriptome study revealed a total of 15,785 differentially expressed genes significantly enriched in the flavonoid biosynthesis pathway. Furthermore, for flavonoid biosynthesis-related genes, such as flavonoid 3′5′-hydroxylase (F3′5′H), anthocyanidin synthase (ANS), anthocyanidin reductase (ANR), chalcone synthase (CHS), dihydroflavonol 4-reductase (DFR), and chalcone isomerase (CHI), their expressions significantly increased in LCF compared with DCF and WCF. Moreover, transcription factors MYB and bHLH were significantly expressed in LCF and DCF. Most flavonoid-related metabolites (myricetin naringenin, catechin, epicatechin-epiafzelechin, and epigallocatechin) were found to be more highly up-regulated in LCF and DCF than WCF. These findings reveal the regulatory mechanism controlling different brown pigmentation in cotton fibers and elucidate the need for the proper selection of high-quality brown cotton fiber breeding lines for promising fiber quality and durable brown color pigmentation.

Lipidomic and transcriptomic profiles of glycerophospholipid metabolism during Hemerocallis citrina Baroni flowering

BACKGROUND

Hemerocallis citrina Baroni (daylily) is a horticultural ornamental plant and vegetable with various applications as a raw material in traditional Chinese medicine and as a flavouring agent. Daylily contains many functional substances and is rich in lecithin, which is mostly composed of glycerophospholipids. To study the comprehensive dynamic changes in glycerophospholipid during daylily flowering and the underlying signalling mechanisms, we performed comprehensive, time-resolved lipidomic and transcriptomic analyses of 'Datong Huanghua 6' daylily.

RESULTS

Labelling with PKH67 fluorescent antibodies clearly and effectively helped visualise lipid changes in daylily, while relative conductivity and malonaldehyde content detection revealed that the early stages of flowering were controllable processes; however, differences became non-significant after 18 h, indicating cellular damage. In addition, phospholipase D (PLD) and lipoxygenase (LOX) activities increased throughout the flowering process, suggesting that lipid hydrolysis and oxidation had intensified. Lipidomics identified 558 lipids that changed during flowering, with the most different lipids found 12 h before and 12 h after flowering. Transcriptome analysis identified 13 key functional genes and enzymes in the glycerophospholipid metabolic pathway. The two-way orthogonal partial least squares analysis showed that diacylglycerol diphosphate phosphatase correlated strongly and positively with phosphatidic acid (PA)(22:0/18:2), PA(34:2), PA(34:4), and diacylglycerol(18:2/21:0) but negatively with phospholipase C. In addition, ethanolamine phosphotransferase gene and phospholipid-N-methyltransferase gene correlated positively with phosphatidylethanolamine (PE)(16:0/18:2), PE(16:0/18:3), PE(33:2), and lysophosphatidylcholine (16:0) but negatively with PE(34:1).

CONCLUSIONS

Overall, this study elucidated changes in the glycerophospholipid metabolism pathway during the daylily flowering process, as well as characteristic genes, thus providing a basis for future studies of glycerophospholipids and signal transduction in daylilies.

Metabolomics-centered mining of plant metabolic diversity and function: Past decade and future perspectives

Plants are natural experts in organic synthesis, being able to generate large numbers of specific metabolites with widely varying structures that help them adapt to variable survival challenges. Metabolomics is a research discipline that integrates the capabilities of several types of research including analytical chemistry, statistics, and biochemistry. Its ongoing development provides strategies for gaining a systematic understanding of quantitative changes in the levels of metabolites. Metabolomics is usually performed by targeting either a specific cell, a specific tissue, or the entire organism. Considerable advances in science and technology over the last three decades have propelled us into the era of multi-omics, in which metabolomics, despite at an earlier developmental stage than genomics, transcriptomics, and proteomics, offers the distinct advantage of studying the cellular entities that have the greatest influence on end phenotype. Here, we summarize the state of the art of metabolite detection and identification, and illustrate these techniques with four case study applications: (i) comparing metabolite composition within and between species, (ii) assessing spatio-temporal metabolic changes during plant development, (iii) mining characteristic metabolites of plants in different ecological environments and upon exposure to various stresses, and (iv) assessing the performance of metabolomics as a means of functional gene identification , metabolic pathway elucidation, and metabolomics-assisted breeding through analyzing plant populations with diverse genetic variations. In addition, we highlight the prominent contributions of joint analyses of plant metabolomics and other omics datasets, including those from genomics, transcriptomics, proteomics, epigenomics, phenomics, microbiomes, and ion-omics studies. Finally, we discuss future directions and challenges exploiting metabolomics-centered approaches in understanding plant metabolic diversity.

Transcriptome and metabolome analyses reveal the key genes related to grain size of big grain mutant in Tartary Buckwheat (Fagopyrum tartaricum)

Grain size with high heritability and stability is an important selection target during Tartary buckwheat breeding. However, the mechanisms that regulate Tartary buckwheat grain development are unknown. We generated transcriptome and metabolome sequencing from 10 and 15 days past anthesis (DPA) grains of big grain mutant (bg1) and WT, and identified 4108 differentially expressed genes (DEGs) including 93 significantly up-regulated differential genes and 85 significantly down-regulated genes in both stages, simultaneously. Meanwhile, we identified DEGs involved in ubiquitin-proteasome pathway, HAI-KU (IKU) pathway, mitogen-activated protein kinase (MAPK) signaling pathway, plant hormone (auxin, brassinosteroids and cytokinins) transduction pathway and five transcription factor families, including APETALA (AP2), GROWTH-REGULATING FACTORS (GRF), AUXIN RESPONSE FACTOR (ARF), WRKY and MYB. Weighted gene co-expression network analysis (WGCNA) was performed and obtained 9 core DEGs. Conjoint analyses of transcriptome and metabolome sequencing screened out 394 DEGs. Using a combined comprehensive analysis, we identified 24 potential candidate genes that encode E3 ubiquitin-protein ligase HIP1, EMBRYO-DEFECTIVE (EMB) protein, receptor-like protein kinase FERONIA (FER), kinesin-4 protein SRG1, and so on, which may be associated with the big-grain mutant bg1. Finally, a quantitative real-time Polymerase Chain Reaction (qRT-PCR) assay was conducted to validate the identified DEGs. Our results provide additional knowledge for identification and functions of causal candidate genes responsible for the variation in grain size and will be an invaluable resource for the genetic dissection of Tartary buckwheat high-yield molecular breeding.

Integrated Transcriptome and Metabolome Analysis to Identify Sugarcane Gene Defense against Fall Armyworm (Spodoptera frugiperda) Herbivory

Sugarcane is the most important sugar crop, contributing ≥80% to total sugar production around the world. Spodoptera frugiperda is one of the main pests of sugarcane, potentially causing severe yield and sugar loss. The identification of key defense factors against S. frugiperda herbivory can provide targets for improving sugarcane resistance to insect pests by molecular breeding. In this work, we used one of the main sugarcane pests, S. frugiperda, as the tested insect to attack sugarcane. Integrated transcriptome and metabolomic analyses were performed to explore the changes in gene expression and metabolic processes that occurred in sugarcane leaf after continuous herbivory by S. frugiperda larvae for 72 h. The transcriptome analysis demonstrated that sugarcane pest herbivory enhanced several herbivory-induced responses, including carbohydrate metabolism, secondary metabolites and amino acid metabolism, plant hormone signaling transduction, pathogen responses, and transcription factors. Further metabolome analysis verified the inducement of specific metabolites of amino acids and secondary metabolites by insect herbivory. Finally, association analysis of the transcriptome and metabolome by the Pearson correlation coefficient method brought into focus the target defense genes against insect herbivory in sugarcane. These genes include amidase and lipoxygenase in amino acid metabolism, peroxidase in phenylpropanoid biosynthesis, and pathogenesis-related protein 1 in plant hormone signal transduction. A putative regulatory model was proposed to illustrate the sugarcane defense mechanism against insect attack. This work will accelerate the dissection of the mechanism underlying insect herbivory in sugarcane and provide targets for improving sugarcane variety resistance to insect herbivory by molecular breeding.

Comparative Metabolome and Transcriptome Analysis of Anthocyanin Biosynthesis in White and Pink Petals of Cotton (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) is one of the important fiber crops. Cotton flowers usually appear white (or cream-colored) without colored spots at the petal base, and turn pink on the next day after flowering. In this study, using a mutant showing pink petals with crimson spots at their base, we conducted comparative metabolome and transcriptome analyses to investigate the molecular mechanism of coloration in cotton flowers. Metabolic profiling showed that cyanidin-3-O-glucoside and glycosidic derivatives of pelargonidins and peonidins are the main pigments responsible for the coloration of the pink petals of the mutant. A total of 2443 genes differentially expressed (DEGs) between the white and pink petals were identified by RNA-sequencing. Many DEGs are structural genes and regulatory genes of the anthocyanin biosynthesis pathway. Among them, MYB21, UGT88F3, GSTF12, and VPS32.3 showed significant association with the accumulation of cyanidin-3-O-glucoside in the pink petals. Taken together, our study preliminarily revealed the metabolites responsible for the pink petals and the key genes regulating the biosynthesis and accumulation of anthocyanins in the pink petals. The results provide new insights into the biochemical and molecular mechanism underlying anthocyanin biosynthesis in upland cotton.

Series-temporal transcriptome profiling of cotton reveals the response mechanism of phosphatidylinositol signaling system in the early stage of drought stress

Plants are sessile organisms suffering severe environmental conditions. Drought stress is one of the major environmental issues that affect plant growth and productivity. Although complex regulatory gene networks of plants under drought stress have been analyzed extensively, the response mechanism in the early stage of drought stress is still rarely mentioned. Here, we performed transcriptome analyses on cotton samples treated for a short time (10 min, 30 min, 60 min, 180 min) using 10% PEG, which is used to simulate drought stress. The analysis of differently expressed genes (DEGs) showed that the number of DEGs in roots was obviously more than that in stems and leaves at the four time points and maintained >2000 FDEGs (DEGs appearing for the first time) from 10 min, indicating that root tissues of plants respond to drought stress quickly and continuously strongly. Gene ontology (GO) analysis showed that DEGs in roots were mainly enriched in protein modification and microtubule-based process. DEGs were found significantly enriched in phosphatidylinositol signaling system at 10 min through Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, implying the great importance of phosphatidylinositol signal in the early stage of drought stress. What was more, two co-expression modules, which were significantly positively correlated with drought stress, were found by Weighted Gene Co-expression Network Analysis (WGCNA). From one of the co-expression modules, we identified a hub-gene Gohir.A07G058200, which is annotated as "phosphatidylinositol 3- and 4-kinase" in phosphatidylinositol signaling system, and found this gene may interact with auxin-responsive protein. This result suggested that Gohir.A07G058200 may be involved in the crosstalk of phosphatidylinositol signal and auxin signal in the early stage of drought stress. In summary, through transcriptome sequencing, we found that phosphatidylinositol signaling system is an important signal transduction pathway in early stage in response to drought stress, and it may interact with auxin signal transduction through phosphatidylinositol 3- and 4-kinase.

Function deficiency of GhOMT1 causes anthocyanidins over-accumulation and diversifies fibre colours in cotton (Gossypium hirsutum)

Naturally coloured cotton (NCC) fibres need little or no dyeing process in textile industry to low-carbon emission and are environment-friendly. Proanthocyanidins (PAs) and their derivatives were considered as the main components causing fibre coloration and made NCCs very popular and healthy, but the monotonous fibre colours greatly limit the wide application of NCCs. Here a G. hirsutum empurpled mutant (HS2) caused by T-DNA insertion is found to enhance the anthocyanidins biosynthesis and accumulate anthocyanidins in the whole plant. HPLC and LC/MS-ESI analysis confirmed the anthocyanidins methylation and peonidin, petunidin and malvidin formation are blocked. The deficiency of GhOMT1 in HS2 was associated with the activation of the anthocyanidin biosynthesis and the altered components of anthocyanidins. The transcripts of key genes in anthocyanidin biosynthesis pathway are significantly up-regulated in HS2, while transcripts of the genes for transport and decoration were at similar levels as in WT. To investigate the potential mechanism of GhOMT1 deficiency in cotton fibre coloration, HS2 mutant was crossed with NCCs. Surprisingly, offsprings of HS2 and NCCs enhanced PAs biosynthesis and increased PAs levels in their fibres from the accumulated anthocyanidins through up-regulated GhANR and GhLAR. As expected, multiple novel lines with improved fibre colours including orange red and navy blue were produced in their generations. Based on this work, a new strategy for breeding diversified NCCs was brought out by promoting PA biosynthesis. This work will help shed light on mechanisms of PA biosynthesis and bring out potential molecular breeding strategy to increase PA levels in NCCs.

Revealing Genetic Differences in Fiber Elongation between the Offspring of Sea Island Cotton and Upland Cotton Backcross Populations Based on Transcriptome and Weighted Gene Coexpression Networks

Fiber length is an important indicator of cotton fiber quality, and the time and rate of cotton fiber cell elongation are key factors in determining the fiber length of mature cotton. To gain insight into the differences in fiber elongation mechanisms in the offspring of backcross populations of Sea Island cotton Xinhai 16 and land cotton Line 9, we selected two groups with significant differences in fiber length (long-fiber group L and short-fiber group S) at different fiber development stages 0, 5, 10 and 15 days post-anthesis (DPA) for transcriptome comparison. A total of 171.74 Gb of clean data was obtained by RNA-seq, and eight genes were randomly selected for qPCR validation. Data analysis identified 6055 differentially expressed genes (DEGs) between two groups of fibers, L and S, in four developmental periods, and gene ontology (GO) term analysis revealed that these DEGs were associated mainly with microtubule driving, reactive oxygen species, plant cell wall biosynthesis, and glycosyl compound hydrolase activity. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that plant hormone signaling, mitogen-activated protein kinase (MAPK) signaling, and starch and sucrose metabolism pathways were associated with fiber elongation. Subsequently, a sustained upregulation expression pattern, profile 19, was identified and analyzed using short time-series expression miner (STEM). An analysis of the weighted gene coexpression network module uncovered 21 genes closely related to fiber development, mainly involved in functions such as cell wall relaxation, microtubule formation, and cytoskeletal structure of the cell wall. This study helps to enhance the understanding of the Sea Island–Upland backcross population and identifies key genes for cotton fiber development, and these findings will provide a basis for future research on the molecular mechanisms of fiber length formation in cotton populations.

Identification of Candidate Genes for Lint Percentage and Fiber Quality Through QTL Mapping and Transcriptome Analysis in an Allotetraploid Interspecific Cotton CSSLs Population

Upland cotton (Gossypium hirsutum) has long been an important fiber crop, but the narrow genetic diversity of modern G. hirsutum limits the potential for simultaneous improvement of yield and fiber quality. It is an effective approach to broaden the genetic base of G. hirsutum through introgression of novel alleles from G. barbadense with excellent fiber quality. In the present study, an interspecific chromosome segment substitution lines (CSSLs) population was established using G. barbadense cultivar Pima S-7 as the donor parent and G. hirsutum cultivar CCRI35 as the recipient parent. A total of 105 quantitative trait loci (QTL), including 85 QTL for fiber quality and 20 QTL for lint percentage (LP), were identified based on phenotypic data collected from four environments. Among these QTL, 25 stable QTL were detected in two or more environments, including four for LP, eleven for fiber length (FL), three for fiber strength (FS), six for fiber micronaire (FM), and one for fiber elongation (FE). Eleven QTL clusters were observed on nine chromosomes, of which seven QTL clusters harbored stable QTL. Moreover, eleven major QTL for fiber quality were verified through analysis of introgressed segments of the eight superior lines with the best comprehensive phenotypes. A total of 586 putative candidate genes were identified for 25 stable QTL associated with lint percentage and fiber quality through transcriptome analysis. Furthermore, three candidate genes for FL, GH_A08G1681 (GhSCPL40), GH_A12G2328 (GhPBL19), and GH_D02G0370 (GhHSP22.7), and one candidate gene for FM, GH_D05G1346 (GhAPG), were identified through RNA-Seq and qRT-PCR analysis. These results lay the foundation for understanding the molecular regulatory mechanism of fiber development and provide valuable information for marker-assisted selection (MAS) in cotton breeding.

Transcriptome Analysis Reveals Differences in Anthocyanin Accumulation in Cotton (Gossypium hirsutum L.) Induced by Red and Blue Light

Many factors, including illumination, affect anthocyanin biosynthesis and accumulation in plants. light quality is the key factor affecting the process of photoinduced anthocyanin biosynthesis and accumulation. We observed that the red color of the Upland cotton accession Huiyuan with the R1 mutation turned to normal green color under light-emitting diodes (LEDs), which inspired us to investigate the effect of red and blue lights on the biosynthesis and accumulation of anthocyanins. We found that both red and blue lights elevated accumulation of anthocyanins. Comparative transcriptomic analyses, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and GSEA, revealed that genes differentially expressed under different light conditions were enriched with the pathways of circadian rhythm, phenylpropanoid biosynthesis, anthocyanin biosynthesis, and flavone and flavonol biosynthesis. Not surprisingly, all the major structural genes related to biosynthesis of anthocyanins, including the key regulatory MYB transcription factor (GhPAP1D) and anthocyanin transporter (GhGSTF12), were induced by red or blue light treatment. However, LARs and MATEs related to biosynthesis of proanthocyanidins were more significantly up-regulated by red light radiation than by blue light radiation. Vice versa, the accumulation of anthocyanins under red light was not as high as that under blue light. In addition, we demonstrated a potential role of GhHY5, a key regulator in plant circadian rhythms, in regulation of anthocyanin accumulation, which could be achieved via interaction with GhPAP1D. Together, these results indicate different effect of red and blue lights on biosynthesis and accumulation of anthocyanins and a potential module including GhHY5 and GhPAP1D in regulation of anthocyanin accumulation in cotton. These results also suggest that the substrates responsible the synthesis of anthocyanins under blue light is diverted to biosynthesis of proanthocyanidin under red light.

Rice metabolic regulatory network spanning the entire life cycle

As one of the most important crops in the world, rice (Oryza sativa) is a model plant for metabolome research. Although many studies have focused on the analysis of specific tissues, the changes in metabolite abundance across the entire life cycle have not yet been determined. In this study, combining both targeted and nontargeted metabolite profiling methods, a total of 825 annotated metabolites were quantified in rice samples from different tissues covering the entire life cycle. The contents of metabolites in different tissues of rice were significantly different, with various metabolites accumulating in the plumule and radicle during seed germination. Combining these data with transcriptome data obtained from the same time period, we constructed the Rice Metabolic Regulation Network. The metabolites and co-expressed genes were further divided into 12 clusters according to their accumulation patterns, with members within each cluster displaying a uniform and clear pattern of abundance across development. Using this dataset, we established a comprehensive metabolic profile of the rice life cycle and used two independent strategies to identify novel transcription factors-namely the use of known regulatory genes as bait to screen for new networks underlying lignin metabolism and the unbiased identification of new glycerophospholipid metabolism regulators on the basis of tissue specificity. This study thus demonstrates how guilt-by-association analysis of metabolome and transcriptome data spanning the entire life cycle in cereal crops provides novel resources and tools to aid in understanding the mechanisms underlying important agronomic traits.

Combined transcriptome and metabolome integrated analysis of Acer mandshuricum to reveal candidate genes involved in anthocyanin accumulation

The red color formation of Acer mandshuricum leaves is caused by the accumulation of anthocyanins primarily, but the molecular mechanism researches which underlie anthocyanin biosynthesis in A. mandshuricum were still lacking. Therefore, we combined the transcriptome and metabolome and analyzed the regulatory mechanism and accumulation pattern of anthocyanins in three different leaf color states. In our results, 26 anthocyanins were identified. Notably, the metabolite cyanidin 3-O-glucoside was found that significantly correlated with the color formation, was the predominant metabolite in anthocyanin biosynthesis of A. mandshuricum. By the way, two key structural genes ANS (Cluster-20561.86285) and BZ1 (Cluster-20561.99238) in anthocyanidin biosynthesis pathway were significantly up-regulated in RL, suggesting that they might enhance accumulation of cyanidin 3-O-glucoside which is their downstream metabolite, and contributed the red formation of A. mandshuricum leaves. Additionally, most TFs (e.g., MYBs, bZIPs and bHLHs) were detected differentially expressed in three leaf color stages that could participate in anthocyanin accumulation. This study sheds light on the anthocyanin molecular regulation of anthocyanidin biosynthesis and accumulation underlying the different leaf color change periods in A. mandshuricum, and it could provide basic theory and new insight for the leaf color related genetic improvement of A. mandshuricum.

Breeding Next-Generation Naturally Colored Cotton

Use of naturally colored cotton (NCC) in the textile industry is limited due to its inferior fiber quality. Recent advances in understanding pigmentation in NCC fiber and the genetic relationship between pigmentation and fiber development may overcome this hindrance and can guide developing high fiber quality and yield NCC.

GhGSTF12, a glutathione S-transferase gene, is essential for anthocyanin accumulation in cotton (Gossypium hirsutum L.)

Anthocyanins are flavonoid pigments providing plants a range of colors from red, pink, orange to blue. Anthocyanins are synthesized in the cytosol but accumulate predominantly in the vacuoles through vacuolar sequestration involving glutathione S-transferases (GSTs) and multidrug and toxic compound extrusion (MATE) and the ATP binding cassette (ABC) transporters. However, little is known about anthocyanin-related GSTs in Upland cotton (Gossypium hirsutum L.). In this study, we performed genome-wide identification of GST genes in Upland cotton and identified GST genes functioning in accumulation of anthocyanins. We demonstrated that GhGSTF12 was able to complement the defective leaf color phenotypes of the Arabidopsis tt19 mutant caused by mutation in a GSTF gene. Virus-induced silencing of GhGSTF12 in the red leaf cultivar turned its red color to green and transient overexpression of GhGSTF12 accelerated anthocyanin accumulation in the red leaf cultivar but not in the green leaf cultivar. Collectively, GhGSTF12 may be involved in transport of anthocyanins from cytosol to vacuoles in cotton. These results also demonstrated a conserved function of plant GSTF genes in anthocyanin accumulation and provide a candidate gene for manipulating pigmentation in cotton tissues.

The Cotton Lignin Biosynthetic Gene Gh4CL30 Regulates Lignification and Phenolic Content and Contributes to Verticillium Wilt Resistance

Verticillium wilt is a vascular disease causing tremendous damage to cotton production worldwide. However, our knowledge of the mechanisms of cotton resistance or susceptibility to this disease is very limited. In this study, we compared the defense transcriptomes of cotton (Gossypium hirsutum) cultivars Shidalukang 1 (Verticillium dahliae resistant, HR) and Junmian 1 (V. dahliae susceptible, HS) before and after V. dahliae infection, identified hub genes of the network associated with responses to V. dahliae infection, and functionally characterized one of the hub genes involved in biosynthesis of lignin and phenolics. We identified 6,831 differentially expressed genes (DEGs) between the basal transcriptomes of HR and HS; 3,685 and 3,239 of these DEGs were induced in HR and HS, respectively, at different time points after V. dahliae infection. KEGG pathway analysis indicated that DEGs were enriched for genes involved in lignin biosynthesis. In all, 23 hub genes were identified based on a weighted gene coexpression network analysis of the 6,831 DEGs and their expression profiles at different time points after V. dahliae infection. Knockdown of Gh4CL30, one of the hub genes related to the lignin biosynthesis pathway, by virus-induced gene silencing, led to a decreased content of flavonoids, lignin, and S monomer but an increased content of G monomer, G/S lignin monomer, caffeic acid, and ferulic acid, and enhanced cotton resistance to V. dahliae. These results suggest that Gh4CL30 is a key gene modulating the outputs of different branches of the lignin biosynthesis pathway, and provide new insights into cotton resistance to V. dahliae.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Comparative transcriptome analysis of leaves during early stages of chilling stress in two different chilling-tolerant brown-fiber cotton cultivars

Chilling stress generates significant inhibition of normal growth and development of cotton plants and lead to severe reduction of fiber quality and yield. Currently, little is known for the molecular mechanism of brown-fiber cotton (BFC) to respond to chilling stress. Herein, RNA-sequencing (RNA-seq)-based comparative analysis of leaves under 4°C treatment in two different-tolerant BFC cultivars, chilling-sensitive (CS) XC20 and chilling-tolerant (CT) Z1612, was performed to investigate the response mechanism. A total of 72650 unigenes were identified with eight commonly used databases. Venn diagram analysis identified 1194 differentially expressed genes (DEGs) with significant up-regulation in all comparison groups. Furthermore, enrichment analyses of COG and KEGG, as well as qRT-PCR validation, indicated that 279 genes were discovered as up-regulated DEGs (UDEGs) with constant significant increased expression in CT cultivar Z1612 groups at the dimensions of both each comparison group and treatment time, locating in the enriched pathways of signal transduction, protein and carbohydrate metabolism, and cell component. Moreover, the comprehensive analyses of gene expression, physiological index and intracellular metabolite detections, and ascorbate antioxidative metabolism measurement validated the functional contributions of these identified candidate genes and pathways to chilling stress. Together, this study for the first time report the candidate key genes and metabolic pathways responding to chilling stress in BFC, and provide the effective reference for understanding the regulatory mechanism of low temperature adaptation in cotton.

Transcriptome Analysis and RNA Interference Reveal GhGDH2 Regulating Cotton Resistance to Verticillium Wilt by JA and SA Signaling Pathways

Verticillium wilt, caused by Verticillium dahliae, is one of the most damaging and widespread soil-borne cotton diseases. The molecular mechanisms underlying the cotton defense against V. dahliae remain largely elusive. Here, we compared the transcriptional differences between Upland cotton cultivars: one highly resistant (HR; Shidalukang 1) and one highly susceptible (HS; Junmian 1). This was done at multiple time points after V. dahliae inoculation, which identified 2010 and 1275 differentially expressed genes (DEGs) in HR and HS, respectively. Plant hormone signal transduction-related genes were enriched in HR, whereas genes related to lignin biosynthesis were enriched in both HR and HS. Weighted gene co-expression network analysis (WGCNA) using the 2868 non-redundant genes differentially expressed between the V. dahliae infected and uninfected samples in HR or HS identified 10 different gene network modules and 22 hub genes with a potential role in regulating cotton defense against V. dahliae infection. GhGDH2, encoding glutamate dehydrogenase (GDH), was selected for functional characterization. Suppressing the expression level of GhGDH2 by virus-induced gene silencing (VIGS) in HS led to inhibition of the salicylic acid (SA) biosynthesis/signaling pathways and activation of the jasmonic acid (JA) biosynthesis/signaling pathways, which resulted in an increase of 42.1% JA content and a reduction of 78.9% SA content in cotton roots, and consequently enhanced V. dahliae resistance. Our finding provides new insights on the molecular mechanisms of cotton resistance to V. dahliae infection and candidate genes for breeding V. dahliae resistance cotton cultivars by genetic modification.

Quantitative metabolome and transcriptome analysis reveals complex regulatory pathway underlying photoinduced fiber color formation in cotton

As an important plant single cell model and textile application materials, poorly known about fiber color formation in cotton, which is sensitively regulated by environmental signals. Our studies underline the importance of photo signal on sensitive fiber color formation and characterize fiber color early initiation (15 DPA) and late accumulated metabolites (45 DPA) in different lighting condition. The results revealed 236 differential metabolites between control and shading, of which phenylpropanoids metabolites accounted for 20%, including uncharacterized novel metabolites and pathways. Furthermore, the early initiation specific genes respond to the absence of light are highly correlated with phenylpropanoid metabolites related to pigmentation. The current study reveals the complex pathways involving early initiation regulation and late metabolic pathways. In addition, the collection composed of uncharacterized photoinduced metabolites and early initiation signaling/regulatory genes were identified, which are important resources for understanding fiber color formation. This report provides new insight into molecular regulatory and biochemical basis underlying photoinduced fiber color formation in cotton.

Understanding the role of phytohormones in cotton fiber development through omic approaches; recent advances and future directions

Cotton is among the most important fiber crops for the textile-based industry, thanks to its cellulose-rich mature fibers. The fiber initiation and elongation are one of the best models for deciphering mechanisms of single-cell differentiation and growth, that also target of fiber development programs. During the last couple of decades, high yielding omics approaches (genomics, transcriptomics, and proteomics), have helped in the identification of several genes and gene products involved in fiber development along with functional relationship to phytohormones. For example, MYB transcription factor family and Sus gene family have been evidenced by controlling cotton fiber initiation. Most importantly, the biosynthesis, responses, and transporting of phytohormones is documented to participate in the initiation of cotton fibers. Herein, in this review, the reliable genetic evidence by manipulating the above genes in cotton have been summarized to describe the relationships among key phytohormones, transcription factors, proteins, and downstream fiber growth-related genes such as Sus. The effect of other important factors such as ROS, fatty acid metabolism, and actin (globular multi-functional proteins) over fiber development has also been discussed. The challenges and deficiencies in the research of cotton fiber development have been mentioned along with a future perspective to discover new crucial genes using multiple omics analysis.