Comparative Transcriptome Analysis of Gene Expression and Regulatory Characteristics Associated with Different Vernalization Periods in Brassica rapa

Identification of NPF Family Genes in Brassica rapa Reveal Their Potential Functions in Pollen Development and Response to Low Nitrate Stress

Nitrate Transporter 1/Peptide Transporter Family (NPF) genes encode membrane transporters involved in the transport of diverse substrates. However, little is known about the diversity and functions of NPFs in Brassica rapa. In this study, 85 NPFs were identified in B. rapa (BrNPFs) which comprised eight subfamilies. Gene structure and conserved motif analysis suggested that BrNFPs were conserved throughout the genus. Stress and hormone-responsive cis-acting elements and transcription factor binding sites were identified in BrNPF promoters. Syntenic analysis suggested that tandem duplication contributed to the expansion of BrNPFs in B. rapa. Transcriptomic profiling analysis indicated that BrNPF2.6, BrNPF2.15, BrNPF7.6, and BrNPF8.9 were expressed in fertile floral buds, suggesting important roles in pollen development. Thirty-nine BrNPFs were responsive to low nitrate availability in shoots or roots. BrNPF2.10, BrNPF2.19, BrNPF2.3, BrNPF5.12, BrNPF5.16, BrNPF5.8, and BrNPF6.3 were only up-regulated in roots under low nitrate conditions, indicating that they play positive roles in nitrate absorption. Furthermore, many genes were identified in contrasting genotypes that responded to vernalization and clubroot disease. Our results increase understanding of BrNPFs as candidate genes for genetic improvement studies of B. rapa to promote low nitrate availability tolerance and for generating sterile male lines based on gene editing methods.

Identification of the Genes Encoding B3 Domain-Containing Proteins Related to Vernalization of Beta vulgaris

Vernalization is the process of exposure to low temperatures, which is crucial for the transition from vegetative to reproductive growth of plants. In this study, the global landscape vernalization-related mRNAs and long noncoding RNAs (lncRNAs) were identified in Beta vulgaris. A total of 22,159 differentially expressed mRNAs and 4418 differentially expressed lncRNAs were uncovered between the vernalized and nonvernalized samples. Various regulatory proteins, such as zinc finger CCCH domain-containing proteins, F-box proteins, flowering-time-related proteins FY and FPA, PHD finger protein EHD3 and B3 domain proteins were identified. Intriguingly, a novel vernalization-related lncRNA-mRNA target-gene co-expression regulatory network and the candidate vernalization genes, VRN1, VRN1-like, VAL1 and VAL2, encoding B3 domain-containing proteins were also unveiled. The results of this study pave the way for further illumination of the molecular mechanisms underlying the vernalization of B. vulgaris.

Transcriptome analysis reveals anthocyanin regulation in Chinese cabbage (Brassica rapa L.) at low temperatures

Chinese cabbage that prefers cold conditions is also affected by low-temperature stress, such as the accumulation of leaf anthocyanins. Research on anthocyanin biosynthesis and regulation mechanisms has made great progress. However, research on anthocyanin accumulation for resistance to biological and non-biological stress is still lacking. To study the relationship between anthocyanin accumulation of Chinese cabbage and resistance under low-temperature conditions, RNA sequencing (RNA-seq) was performed on Chinese cabbage ‘Xiao Baojian’ grown at a low temperature for four time periods and at a control temperature for five time periods. In Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, 7954 differentially expressed genes (DEGs) were enriched, of which 587 DEGs belonged to "biosynthesis of other secondary metabolites." Gene temporal expression patterns were used to discover enriched genes related to phenylpropanoid biosynthesis; flavonoid biosynthesis and anthocyanin biosynthesis pathways were found in cluster 1. The interaction networks were constructed, and hub genes were selected, showing that flavonoid biosynthesis pathway genes (DFR, ANS, F3H, FLS1, CHS1, CHS3, and TT8) and defense mechanisms-related genes (DFR, SNL6, and TKPR1) interact with each other. Anthocyanin biosynthesis DEGs in Chinese cabbage were evaluated under low-temperature conditions to map the relevant pathways, and expression maps of transcription factors in the flavonoid pathway were created at various periods. Low temperature upregulated the expression of genes related to anthocyanin biosynthesis. Taken together, our results provide further analysis of the relationship between plant anthocyanin synthesis and stress resistance and may also provide further insights for the future development of high-quality color and cold-tolerant Chinese cabbage germplasm resources.

Whole-transcriptome sequencing reveals a vernalization-related ceRNA regulatory network in chinese cabbage (Brassica campestris L. ssp. pekinensis)

BACKGROUND

The transition from vegetative growth to reproductive growth involves various pathways. Vernalization is a crucial process for floral organ formation and regulation of flowering time that is widely utilized in plant breeding. In this study, we aimed to identify the global landscape of mRNAs, microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) related to vernalization in Chinese cabbage. These data were then used to construct a competitive endogenous RNA (ceRNA) network that provides valuable information to better understand the vernalization response.

RESULTS

In this study, seeds sampled from the Chinese cabbage doubled haploid (DH) line 'FT' with or without vernalization treatment were used for whole-transcriptome sequencing. A total of 2702 differentially expressed (DE) mRNAs, 151 DE lncRNAs, 16 DE circRNAs, and 233 DE miRNAs were identified in the vernalization-treated seeds. Various transcription factors, such as WRKY, MYB, NAC, bHLH, MADS-box, zinc finger protein CONSTANS-like gene, and B3 domain protein, and regulatory proteins that play important roles in the vernalization pathway were identified. Additionally, we constructed a vernalization-related ceRNA-miRNA-target gene network and obtained 199 pairs of ceRNA relationships, including 108 DEmiRNA‒DEmRNA, 67 DEmiRNA‒DElncRNA, and 12 DEmiRNA‒DEcircRNA interactions, in Chinese cabbage. Furthermore, several important vernalization-related genes and their interacting lncRNAs, circRNAs, and miRNAs, which are involved in the regulation of flowering time, floral organ formation, bolting, and flowering, were identified.

CONCLUSIONS

Our results reveal the potential mRNA and non-coding RNAs involved in vernalization, providing a foundation for further studies on the molecular mechanisms underlying vernalization in Chinese cabbage.

Germinating seed can sense low temperature for the floral transition and vernalization of winter rapeseed (Brassica rapa)

The hybrid production of winter rapeseed is limited by the difficult vernalization processes. Thus, floral regulation of winter rapeseed parental lines cannot be executed through selection of sowing time during hybrid production. Therefore, in this study, strong winter rapeseed was used as the material to analyse the floral transition mechanism of germinating seed vernalization. Results demonstrated that germinating seeds could sense low temperatures and complete vernalization following a low temperature treatment for 56.5 d with a 100 % vernalization rate. The regression equation between vernalization rate (y) and vernalization treatment days (x) was determined as y = 0.019x - 0.0765 (R² = 0.8529). When the vernalization treatment time was prolonged, the vernalization rate and fruiting ability increased rapidly, and variations were observed in the membrane lipid oxidation and physiological characteristics. Furthermore, at the prolonged treatment time of 10-50 d, the salicylic acid (SA) content continued to decrease, with values significantly lower than those of the control. SA content is significantly positively correlated with the level of BrFLC transcription and a significantly negatively correlated with the vernalization rate of germinating seeds. Moreover, the expressions of genes associated with SA biosynthesis, SA signal transduction, the flowering key negative regulators were suppressed and that of positive regulators were promoted during vernalization. These results suggest that SA as a floral repressor is involved in the regulation of the vernalization process of winter rapeseed germination seeds. In addition, SA may be related to the counting dosage of vernalization.

Gene co-expression network analysis of the heat-responsive core transcriptome identifies hub genes in Brassica rapa

Gene co-expression network analysis of the heat-responsive core transcriptome in two contrasting Brassica rapa accessions reveals the main metabolic pathways, key modules and hub genes, are involved in long-term heat stress. Brassica rapa is a widely cultivated and economically important vegetable in Asia. High temperature is a common stress that severely impacts leaf head formation in B. rapa, resulting in reduced quality and production. The purpose of this study was thus to identify candidate heat tolerance genes by comparative transcriptome analysis of two contrasting B. rapa accessions in response to long-term heat stress. Two B. rapa accessions, '268' and '334', which showed significant differences in heat tolerance, were used for RNA sequencing analysis. We identified a total of 11,055 and 8921 differentially expressed genes (DEGs) in '268' and '334', respectively. Functional enrichment analyses of all of the identified DEGs, together with the genes identified from weighted gene co-expression network analyses (WGCNA), revealed that the autophagy pathway, glutathione metabolism, and ribosome biogenesis in eukaryotes were significantly up-regulated, whereas photosynthesis was down-regulated, in the heat resistance of B. rapa '268'. Furthermore, when B. rapa '334' was subjected to long-term high-temperature stress, heat stress caused significant changes in the expression of certain functional genes linked to protein processing in the endoplasmic reticulum and plant hormone signal transduction pathways. Autophagy-related genes might have been induced by persistent heat stress and remained high during recovery. Several hub genes like HSP17.6, HSP17.6B, HSP70-8, CLPB1, PAP1, PYR1, ADC2, and GSTF11 were discussed in this study, which may be potential candidates for further analyses of the response to long-term heat stress. These results should help elucidate the molecular mechanisms of heat stress adaptation in B. rapa.

Transcriptome architecture reveals genetic networks of bolting regulation in spinach

BACKGROUND

Bolting refers to the early flowering stem production on agricultural and horticultural crops before harvesting. Indeed, bolting is an event induced by the coordinated effects of various environmental factors and endogenous genetic components, which cause a large reduction in the quality and productivity of vegetable crops like spinach. However, little is known about the signaling pathways and molecular functions involved in bolting mechanisms in spinach. The genetic information regarding the transition from vegetative growth to the reproductive stage in spinach would represent an advantage to regulate bolting time and improvement of resistant cultivars to minimize performance loss.

RESULTS

To investigate the key genes and their genetic networks controlling spinach bolting, we performed RNA-seq analysis on early bolting accession Kashan and late-bolting accession Viroflay at both vegetative and reproductive stages and found a significant number of differentially expressed genes (DEGs) ranging from 195 to 1230 in different comparisons. These genes were mainly associated with the signaling pathways of vernalization, photoperiod/circadian clock, gibberellin, autonomous, and aging pathways. Gene ontology analysis uncovered terms associated with carbohydrate metabolism, and detailed analysis of expression patterns for genes of Fructose-1, 6-bisphosphate aldolase, TREHALOSE-6-PHOSPHATE SYNTHASE 1, FLOWERING PROMOTING FACTOR 1, EARLY FLOWERING, GIGANTEA, and MADS-box proteins revealed their potential roles in the initiating or delaying of bolting.

CONCLUSION

This study is the first report on identifying bolting and flowering-related genes based on transcriptome sequencing in spinach, which provides insight into bolting control and can be useful for molecular breeding programs and further study in the regulation of the genetic mechanisms related to bolting in other vegetable crops.

Gene co-expression network analysis reveals key pathways and hub genes in Chinese cabbage (Brassica rapa L.) during vernalization

BACKGROUND

Vernalization is a type of low temperature stress used to promote rapid bolting and flowering in plants. Although rapid bolting and flowering promote the reproduction of Chinese cabbages (Brassica rapa L. ssp. pekinensis), this process causes their commercial value to decline. Clarifying the mechanisms of vernalization is essential for its further application. We performed RNA sequencing of gradient-vernalization in order to explore the reasons for the different bolting process of two Chinese cabbage accessions during vernalization.

RESULTS

There was considerable variation in gene expression between different-bolting Chinese cabbage accessions during vernalization. Comparative transcriptome analysis and weighted gene co-expression network analysis (WGCNA) were performed for different-bolting Chinese cabbage during different vernalization periods. The biological function analysis and hub gene annotation of highly relevant modules revealed that shoot system morphogenesis and polysaccharide and sugar metabolism caused early-bolting 'XBJ' to bolt and flower faster; chitin, ABA and ethylene-activated signaling pathways were enriched in late-bolting 'JWW'; and leaf senescence and carbohydrate metabolism enrichment were found in the two Chinese cabbage-related modules, indicating that these pathways may be related to bolting and flowering. The high connectivity of hub genes regulated vernalization, including MTHFR2, CPRD49, AAP8, endoglucanase 10, BXLs, GATLs, and WRKYs. Additionally, five genes related to flower development, BBX32 (binds to the FT promoter), SUS1 (increases FT expression), TSF (the closest homologue of FT), PAO and NAC029 (plays a role in leaf senescence), were expressed in the two Chinese cabbage accessions.

CONCLUSION

The present work provides a comprehensive overview of vernalization-related gene networks in two different-bolting Chinese cabbages during vernalization. In addition, the candidate pathways and hub genes related to vernalization identified here will serve as a reference for breeders in the regulation of Chinese cabbage production.

Comparative transcriptome analysis revealed differential gene expression in multiple signaling pathways at flowering in polyploid Brassica rapa

BACKGROUND

Polyploidy is widespread in angiosperms and has a significant impact on plant evolution, diversity, and breeding program. However, the changes in the flower development regulatory mechanism in autotetraploid plants remains relatively limited. In this study, RNA-seq analysis was used to investigate changes in signaling pathways at flowering in autotetraploid Brassica rapa.

RESULTS

The study findings showed that the key genes such as CO, CRY2, and FT which promotes floral formation were down-regulated, whereas floral transition genes FPF1 and FD were up-regulated in autotetraploid B. rapa. The data also demonstrated that the positive regulators GA1 and ELA1 in the gibberellin's biosynthesis pathway were negatively regulated by polyploidy in B. rapa. Furthermore, transcriptional factors (TFs) associated with flower development were significantly differentially expressed including the up-regulated CIB1 and AGL18, and the down-regulated AGL15 genes, and by working together such genes affected the expression of the down-stream flowering regulator FLOWERING LOCUS T in polyploid B. rapa. Compared with that in diploids autotetrapoid plants consist of differential expression within the signaling transduction pathway, with 13 TIFY gens up-regulated and 17 genes related to auxin pathway down-regulated.

CONCLUSION

Therefore, polyploidy is more likely to integrate multiple signaling pathways to influence flowering in B. rapa after polyploidization. In general, the present results shed new light on our global understanding of flowering regulation in polyploid plants during breeding program.