Global transcriptome and coexpression network analyses reveal cultivar-specific molecular signatures associated with seed development and seed size/weight determination in chickpea

Insights into Bactericera cockerelli and Candidatus Liberibacter solanacearum interaction: a tissue-specific transcriptomic approach

The tomato-potato psyllid, Bactericera cockerelli (Šulc), belonging to the Hemiptera order, is an insect pest of solanaceous crops and vectors a fastidious bacterium, Candidatus Liberibacter solanacearum (CLso), the presumptive causal agent of zebra chip and vein greening diseases in potatoes and tomatoes, respectively. The genome of B. cockerelli has been sequenced recently, providing new avenues to elucidate mechanistic insights into pathogenesis in vegetable crops. In this study, we performed RNA-sequencing of the critical psyllid organs (salivary glands and ovaries) involved in CLso pathology and transmission to host plants. Transcriptome analysis revealed differentially expressed genes and organ-specific enrichment of gene ontology (GO) terms related to metabolic processes, response to stress/stimulus, phagocytosis, proteolysis, endocytosis, and provided candidate genes encoding transcription factors (TFs). To examine gene regulatory networks across the psyllid organs under CLso(-) and CLso(+) conditions, we performed weighted gene co-expression network analysis (WGCNA), and unique modules differentiating the psyllid organs were identified. A comparative GO analysis of the unique gene modules revealed functional terms enriched in response to stress, gene regulation, and cell division processes in the ovaries. In contrast, respiration, transport, and neuronal transmission-related GO terms were enriched in the salivary glands. Altogether, this study reveals new insights into tissue-specific expression of the psyllid organs in the absence or presence of CLso bacterium. This knowledge can be leveraged to develop new pest and disease management strategies by delineating the regulatory networks involved in the psyllid-CLso interaction.

Phenotypic and transcriptomic responses of diverse rice accessions to transient heat stress during early grain development

Heat stress (HS) occurring during the grain-filling period has a detrimental effect on grain yield and quality in rice (Oryza sativa). The development of heat-resilient cultivars could partly solve this issue if tolerant alleles can be identified and incorporated into the germplasm. In this study, we posit that some of the phenotypic variations for heat resilience during grain development could be due to variations in gene expression among accessions. To test this, we characterized the HS response of 10 diverse rice accessions from three major sub-populations using physiological and transcriptome analyses. At a single-grain level, grain width and grain thickness emerged as the most heat-sensitive traits. During a transient HS, IND-3 was categorized as highly sensitive, while five accessions exhibited moderate heat sensitivity, and four accessions were tolerant. Only a core set of 29.4% of the differentially expressed genes was common to the three rice sub-populations. Heat-tolerant accession TEJ-5 uniquely triggered an unfolded protein response (UPR) under HS, as evident from the induction of OsbZIP50 and downstream UPR genes. OsbZIP58, a gene that positively regulates grain filling, was more highly induced by HS in IND-2 despite its moderate heat sensitivity. Collectively, our analysis suggests that both unique gene expression responses and variation in the level of responses for a given pathway distinguish diverse accessions. Only some of these responses are associated with single-grain phenotypes in a manner consistent with the known roles of these genes and pathways.

Genome-wide mapping of DNase I hypersensitive sites revealed differential chromatin accessibility and regulatory DNA elements under drought stress in rice cultivars

Drought stress (DS) is one of the major constraints limiting yield in crop plants including rice. Gene regulation under DS is largely governed by accessibility of the transcription factors (TFs) to their cognate cis-regulatory elements (CREs). In this study, we used DNase I hypersensitive assays followed by sequencing to identify the accessible chromatin regions under DS in a drought-sensitive (IR64) and a drought-tolerant (N22) rice cultivar. Our results indicated that DNase I hypersensitive sites (DHSs) were highly enriched at transcription start sites (TSSs) and numerous DHSs were detected in the promoter regions. DHSs were concurrent with epigenetic marks and the genes harboring DHSs in their TSS and promoter regions were highly expressed. In addition, DS induced changes in DHSs (∆DHSs) in TSS and promoter regions were positively correlated with upregulation of several genes involved in drought/abiotic stress response, those encoding TFs and located within drought-associated quantitative trait loci, much preferentially in the drought-tolerant cultivar. The CREs representing the binding sites of TFs involved in DS response were detected within the ∆DHSs, suggesting differential accessibility of TFs to their cognate sites under DS in different rice cultivars, which may be further deployed for enhancing drought tolerance in rice.

Domestication has altered gene expression and secondary metabolites in pea seed coat

The mature seed in legumes consists of an embryo and seed coat. In contrast to knowledge about the embryo, we know relatively little about the seed coat. We analyzed the gene expression during seed development using a panel of cultivated and wild pea genotypes. Gene co-expression analysis identified gene modules related to seed development, dormancy, and domestication. Oxidoreductase genes were found to be important components of developmental and domestication processes. Proteomic and metabolomic analysis revealed that domestication favored proteins involved in photosynthesis and protein metabolism at the expense of seed defense. Seed coats of wild peas were rich in cell wall-bound metabolites and the protective compounds predominated in their seed coats. Altogether, we have shown that domestication altered pea seed development and modified (mostly reduced) the transcripts along with the protein and metabolite composition of the seed coat, especially the content of the compounds involved in defense. We investigated dynamic profiles of selected identified phenolic and flavonoid metabolites across seed development. These compounds usually deteriorated the palatability and processing of the seeds. Our findings further provide resources to study secondary metabolism and strategies for improving the quality of legume seeds which comprise an important part of the human protein diet.

Genome-wide transcriptome analysis reveals key regulatory networks and genes involved in the determination of seed hardness in vegetable soybean

Seed hardness is an important quality trait of vegetable soybean. To determine the factors underlying seed hardness, two landraces with contrasting seed hardness, Niumaohuang (low seed hardness) and Pixiansilicao (high seed hardness), were selected from 216 soybean accessions originating from 26 provinces in China. The contents of the main components in vegetable soybean seeds such as water, soluble sugar, starch, protein and oil were measured, and transcriptome analyses performed during five stages of seed developmental. Transcriptome analysis indicates that during the middle and late stages of seed development, a large number of genes involved in the synthesis or degradation of starch, storage protein, and fatty acids were differentially expressed, leading to differences in the accumulation of stored substances during seed maturation among Niumaohuang and Pixiansilicao. The activity of cell proliferation and the formation of cell walls in the middle and late stages of seed development may also affect the hardness of seeds to a certain extent. In addition, weighted gene co-expression network analysis (WGCNA) was undertaken to identify co-expressed gene modules and hub genes that regulate seed hardness. Overexpression of a candidate seed hardness regulatory hub gene, GmSWEET2, resulted in increased seed hardness. In this study, the important role of GmSWEET2 in regulating the hardness of vegetable soybean seeds was verified and numerous potential key regulators controlling seed hardness and the proportion of seed components were identified, laying the groundwork for improving the texture of vegetable soybean.

Increased DNA methylation of the splicing regulator SR45 suppresses seed abortion in litchi

The gene regulatory networks that govern seed development are complex, yet very little is known about the genes and processes that are controlled by DNA methylation. Here, we performed single-base resolution DNA methylome analysis and found that CHH methylation increased significantly throughout seed development in litchi. Based on the association analysis of differentially methylated regions and weighted gene co-expression network analysis (WGCNA), 46 genes were identified as essential DNA methylation-regulated candidate genes involved in litchi seed development, including LcSR45, a homolog of the serine/arginine-rich (SR) splicing regulator SR45. LcSR45 is predominately expressed in the funicle, embryo, and seed integument, and displayed increased CHH methylation in the promoter during seed development. Notably, silencing of LcSR45 in a seed-aborted litchi cultivar significantly improved normal seed development, whereas the ectopic expression of LcSR45 in Arabidopsis caused seed abortion. Furthermore, LcSR45-dependent alternative splicing events were found to regulate genes involved in seed development. Together, our findings demonstrate that LcSR45 is hypermethylated, and plays a detrimental role in litchi seed development, indicating a global increase in DNA methylation at this stage.

Integrated Comparative Transcriptome and Weighted Gene Co-Expression Network Analysis Provide Valuable Insights into the Mechanisms of Pinhead Initiation in Chinese Caterpillar Mushroom Ophiocordyceps sinensis (Ascomycota)

The initiation and formation of the "pinhead" is the key node in growth process of Ophiocordyceps sinensis (Chinese Cordyceps). The research on the mechanism of changes in this growth stage is the basis for realizing the industrialization of its artificial cultivation. Clarifying the mechanisms of pinhead initiation is essential for its further application. Here, we performed a comprehensive transcriptome analysis of pinhead initiation process in O. sinensis. Comparative transcriptome analysis revealed remarkable variation in gene expression and enriched pathways at different pinhead initiation stages. Gene co-expression network analysis by WGCNA identified 4 modules highly relevant to different pinhead initiation stages, and 23 hub genes. The biological function analysis and hub gene annotation of these identified modules demonstrated that transmembrane transport and nucleotide excision repair were the topmost enriched in pre-pinhead initiation stage, carbohydrate metabolism and protein glycosylation were specially enriched in pinhead initiation stage, nucleotide binding and DNA metabolic process were over-represented after pinhead stage. These key regulators are mainly involved in carbohydrate metabolism, synthesis of proteins and nucleic acids. This work excavated the candidate pathways and hub genes related to the pinhead initiation stage, which will serve as a reference for realizing the industrialization of artificial cultivation in O. sinensis.

Exploring the role of transcriptional and post-transcriptional processes in mRNA co-expression

Co-expression of two or more genes at the single-cell level is usually associated with functional co-regulation. While mRNA co-expression-measured as the correlation in mRNA levels-can be influenced by both transcriptional and post-transcriptional events, transcriptional regulation is typically considered dominant. We review and connect the literature describing transcriptional and post-transcriptional regulation of co-expression. To enhance our understanding, we integrate four datasets spanning single-cell gene expression data, single-cell promoter activity data and individual transcript half-lives. Confirming expectations, we find that positive co-expression necessitates promoter coordination and similar mRNA half-lives. Surprisingly, negative co-expression is favored by differences in mRNA half-lives, contrary to initial predictions from stochastic simulations. Notably, this association manifests specifically within clusters of genes. We further observe a striking compensation between promoter coordination and mRNA half-lives, which additional stochastic simulations suggest might give rise to the observed co-expression patterns. These findings raise intriguing questions about the functional advantages conferred by this compensation between distal kinetic steps.

Transcriptome and Metabolome Analyses Reveal Sugar and Acid Accumulation during Apricot Fruit Development

The apricot (Prunus armeniaca L.) is a fruit that belongs to the Rosaceae family; it has a unique flavor and is of important economic and nutritional value. The composition and content of soluble sugars and organic acids in fruit are key factors in determining the flavor quality. However, the molecular mechanism of sugar and acid accumulation in apricots remains unclear. We measured sucrose, fructose, glucose, sorbitol, starch, malate, citric acid, titratable acid, and pH, and investigated the transcriptome profiles of three apricots (the high-sugar cultivar 'Shushanggan', common-sugar cultivar 'Sungold', and low-sugar cultivar 'F43') at three distinct developmental phases. The findings indicated that 'Shushanggan' accumulates a greater amount of sucrose, glucose, fructose, and sorbitol, and less citric acid and titratable acid, resulting in a better flavor; 'Sungold' mainly accumulates more sucrose and less citric acid and starch for the second flavor; and 'F43' mainly accumulates more titratable acid, citric acid, and starch for a lesser degree of sweetness. We investigated the DEGs associated with the starch and sucrose metabolism pathways, citrate cycle pathway, glycolysis pathway, and a handful of sugar transporter proteins, which were considered to be important regulators of sugar and acid accumulation. Additionally, an analysis of the co-expression network of weighted genes unveiled a robust correlation between the brown module and sucrose, glucose, and fructose, with VIP being identified as a hub gene that interacted with four sugar transporter proteins (SLC35B3, SLC32A, SLC2A8, and SLC2A13), as well as three structural genes for sugar and acid metabolism (MUR3, E3.2.1.67, and CSLD). Furthermore, we found some lncRNAs and miRNAs that regulate these genes. Our findings provide clues to the functional genes related to sugar metabolism, and lay the foundation for the selection and cultivation of high-sugar apricots in the future.

Transcriptome Dynamics of Brassica juncea Leaves in Response to Omnivorous Beet Armyworm (Spodoptera exigua, Hübner)

Cruciferous plants manufacture glucosinolates (GSLs) as special and important defense compounds against insects. However, how insect feeding induces glucosinolates in Brassica to mediate insect resistance, and how plants regulate the strength of anti-insect defense response during insect feeding, remains unclear. Here, mustard (Brassica juncea), a widely cultivated Brassica plant, and beet armyworm (Spodoptera exigua), an economically important polyphagous pest of many crops, were used to analyze the changes in GSLs and transcriptome of Brassica during insect feeding, thereby revealing the plant-insect interaction in Brassica plants. The results showed that the content of GSLs began to significantly increase after 48 h of herbivory by S. exigua, with sinigrin as the main component. Transcriptome analysis showed that a total of 8940 DEGs were identified in mustard challenged with beet armyworm larvae. The functional enrichment results revealed that the pathways related to the biosynthesis of glucosinolate and jasmonic acid were significantly enriched by upregulated DEGs, suggesting that mustard might provide a defense against herbivory by inducing JA biosynthesis and then promoting GSL accumulation. Surprisingly, genes regulating JA catabolism and inactivation were also activated, and both JA signaling repressors (JAZs and JAMs) and activators (MYCs and NACs) were upregulated during herbivory. Taken together, our results indicate that the accumulation of GSLs regulated by JA signaling, and the regulation of active and inactive JA compound conversion, as well as the activation of JA signaling repressors and activators, collectively control the anti-insect defense response and avoid over-stunted growth in mustard during insect feeding.

Periodic changes in chain lengths distribution parameters of wheat starch during endosperm development

This study analyzed the molecular structure of developing wheat endosperm starch at different stages after anthesis (DAA) using chain length distribution analysis by size exclusion chromatography (SEC) and fluorophore-assisted carbohydrate electrophoresis. Our results revealed periodic changes in the content of both amylose and amylopectin fractions. Specifically, the content of amylose chains with a degree of polymerization (DP) > 100 significantly decreased from 5 to 10 DAA (28% to 21%) and from 15 to 20 DAA (29% to 26%), but increased between 10 and 15 DAA (21% to 29%) and 20 to 25 DAA (30.0% to 33%). Conversely, the content of short amylopectin chains with DP ≤ 32 showed the opposite trend. Interestingly, mRNA expression levels of key starch biosynthesis genes did not exhibit periodic changes. These findings contribute to our understanding of starch biosynthesis and provide important insights for the development of starch-based products.

Functional role of microRNA in the regulation of biotic and abiotic stress in agronomic plants

The increasing demand for food is the result of an increasing population. It is crucial to enhance crop yield for sustainable production. Recently, microRNAs (miRNAs) have gained importance because of their involvement in crop productivity by regulating gene transcription in numerous biological processes, such as growth, development and abiotic and biotic stresses. miRNAs are small, non-coding RNA involved in numerous other biological functions in a plant that range from genomic integrity, metabolism, growth, and development to environmental stress response, which collectively influence the agronomic traits of the crop species. Additionally, miRNA families associated with various agronomic properties are conserved across diverse plant species. The miRNA adaptive responses enhance the plants to survive environmental stresses, such as drought, salinity, cold, and heat conditions, as well as biotic stresses, such as pathogens and insect pests. Thus, understanding the detailed mechanism of the potential response of miRNAs during stress response is necessary to promote the agronomic traits of crops. In this review, we updated the details of the functional aspects of miRNAs as potential regulators of various stress-related responses in agronomic plants.

Global Transcriptome and Co-Expression Network Analyses Revealed Hub Genes Controlling Seed Size/Weight and/or Oil Content in Peanut

Cultivated peanut (Arachis hypogaea L.) is an important economic and oilseed crop worldwide, providing high-quality edible oil and high protein content. Seed size/weight and oil content are two important determinants of yield and quality in peanut breeding. To identify key regulators controlling these two traits, two peanut cultivars with contrasting phenotypes were compared to each other, one having a larger seed size and higher oil content (Zhonghua16, ZH16 for short), while the second cultivar had smaller-sized seeds and lower oil content (Zhonghua6, ZH6). Whole transcriptome analyses were performed on these two cultivars at four stages of seed development. The results showed that ~40% of the expressed genes were stage-specific in each cultivar during seed development, especially at the early stage of development. In addition, we identified a total of 5356 differentially expressed genes (DEGs) between ZH16 and ZH6 across four development stages. Weighted gene co-expression network analysis (WGCNA) based on DEGs revealed multiple hub genes with potential roles in seed size/weight and/or oil content. These hub genes were mainly involved in transcription factors (TFs), phytohormones, the ubiquitin-proteasome pathway, and fatty acid synthesis. Overall, the candidate genes and co-expression networks detected in this study could be a valuable resource for genetic breeding to improve seed yield and quality traits in peanut.

Laser capture microdissection-based spatiotemporal transcriptomes uncover regulatory networks during seed abortion in seedless Ponkan (Citrus reticulata)

Seed abortion is an important process in the formation of seedless characteristics in citrus fruits. However, the molecular regulatory mechanism underlying citrus seed abortion is poorly understood. Laser capture microdissection-based RNA-seq combined with Pacbio-seq was used to profile seed development in the Ponkan cultivars 'Huagan No. 4' (seedless Ponkan) (Citrus reticulata) and 'E'gan No. 1' (seeded Ponkan) (C. reticulata) in two types of seed tissue across three developmental stages. Through comparative transcriptome and dynamic phytohormone analyses, plant hormone signal, cell division and nutrient metabolism-related processes were revealed to play critical roles in the seed abortion of 'Huagan No. 4'. Moreover, several genes may play indispensable roles in seed abortion of 'Huagan No. 4', such as CrWRKY74, CrWRKY48 and CrMYB3R4. Overexpression of CrWRKY74 in Arabidopsis resulted in severe seed abortion. By analyzing the downstream regulatory network, we further determined that CrWRKY74 participated in seed abortion regulation by inducing abnormal programmed cell death. Of particular importance is that a preliminary model was proposed to depict the regulatory networks underlying seed abortion in citrus. The results of this study provide novel insights into the molecular mechanism across citrus seed development, and reveal the master role of CrWRKY74 in seed abortion of 'Huagan No. 4'.

Genetic Analysis for the Flag Leaf Heterosis of a Super-Hybrid Rice WFYT025 Combination Using RNA-Seq

The photosynthetic capacity of flag leaf plays a key role in grain yield in rice. Nevertheless, there are few studies on the heterosis of the rice flag leaf. Therefore, this study focuses on investigating the genetic basis of heterosis for flag leaf in the indica super hybrid rice combination WFYT025 in China using a high-throughput next-generation RNA-seq strategy. We analyzed the gene expression of flag leaf in different environments and different time periods between WFYT025 and its female parent. After obtaining the gene expression profile of the flag leaf, we further investigated the gene regulatory network. Weighted gene expression network analysis (WGCNA) was used to identify the co-expressed gene sets, and a total of 5000 highly expressed genes were divided into 24 co-expression groups. In CHT025, we found 13 WRKY family transcription factors in SDG_hps under the environment of early rice and 16 WRKY family genes in SDG_hps of under the environment of middle rice. We found nine identical transcription factors in the two stages. Except for five reported TFs, the other four TFs might play an important role in heterosis for grain number and photosynthesis. Transcription factors such as WRKY3, WRKY68, and WRKY77 were found in both environments. To eliminate the influence of the environment, we examined the metabolic pathway with the same SDG_hp (SSDG_hp) in two environments. There were 312 SSDG_hps in total. These SSDG_hps mainly focused on the phosphorus metallic process, phosphorylation, plasma membrane, etc. These results provide resources for studying heterosis during super hybrid rice flag leaf development.

Transcriptome analysis of ovules offers early developmental clues after fertilization in Cicer arietinum L

Chickpea (Cicer arietinum L.) seeds are valued for their nutritional scores and limited information on the molecular mechanisms of chickpea fertilization and seed development is available. In the current work, comparative transcriptome analysis was performed on two different stages of chickpea ovules (pre- and post-fertilization) to identify key regulatory transcripts. Two-staged transcriptome sequencing was generated and over 208 million reads were mapped to quantify transcript abundance during fertilization events. Mapping to the reference genome showed that the majority (92.88%) of high-quality Illumina reads were aligned to the chickpea genome. Reference-guided genome and transcriptome assembly yielded a total of 28,783 genes. Of these, 3399 genes were differentially expressed after the fertilization event. These involve upregulated genes including a protease-like secreted in CO(2) response (LOC101500970), amino acid permease 4-like (LOC101506539), and downregulated genes MYB-related protein 305-like (LOC101493897), receptor like protein 29 (LOC101491695). WGCNA analysis and pairwise comparison of datasets, successfully constructed four co-expression modules. Transcription factor families including bHLH, MYB, MYB-related, C2H2 zinc finger, ERF, WRKY and NAC transcription factor were also found to be activated after fertilization. Activation of these genes and transcription factors results in the accumulation of carbohydrates and proteins by enhancing their trafficking and biosynthesis. Total 17 differentially expressed genes, were randomly selected for qRT-PCR for validation of transcriptome analysis and showed statistically significant correlations with the transcriptome data. Our findings provide insights into the regulatory mechanisms underlying changes in fertilized chickpea ovules. This work may come closer to a comprehensive understanding of the mechanisms that initiate developmental events in chickpea seeds after fertilization.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03599-8.

Transcriptome and Anatomical Comparisons Reveal the Effects of Methyl Jasmonate on the Seed Development of Camellia oleifera

Seed is a major storage organ that determines the yield and quality of Camellia oleifera (C. oleifera). Methyl jasmonate (MeJA) is a signaling molecule involved in plant growth and development. However, the role of MeJA in the development of C. oleifera seeds remains a mystery. This study demonstrated that the larger seeds induced by MeJA resulted from more cell numbers and a larger cell area in the outer seed coat and embryo at the cellular level. At the molecular level, MeJA could regulate the expression of factors in the known signaling pathways of seed size control as well as cell proliferation and expansion, resulting in larger seeds. Furthermore, the accumulation of oil and unsaturated fatty acids due to MeJA-inducement was attributed to the increased expression of fatty acid biosynthesis-related genes but reduced expression of fatty acid degradation-related genes. CoMYC2, a key regulator in jasmonate signaling, was considered a potential hub regulator which directly interacted with three hub genes (CoCDKB2-3, CoCYCB2-3, and CoXTH9) related to the seed size and two hub genes (CoACC1 and CoFAD2-3) related to oil accumulation and fatty acid biosynthesis by binding to their promoters. These findings provide an excellent target for the improvement of the yield and quality in C. oleifera.

The developmental dynamics in cool season legumes with focus on chickpea

Chickpea is one of the most widely consumed grain legume world-wide. Advances in next-generation sequencing and genomics tools have led to genetic dissection and identification of potential candidate genes regulating agronomic traits in chickpea. However, the developmental particularities and its potential in reforming the yield and nutritional value remain largely unexplored. Studies in crops such as rice, maize, tomato and pea have highlighted the contribution of key regulator of developmental events in yield related traits. A comprehensive knowledge on the development aspects of a crop can pave way for new vistas to explore. Pea and Medicago are the close relatives of genus Cicer and the basic developmental events in these legumes are similar. However, there are some distinct developmental features in chickpea which hold potential for future crop improvement endeavours. The global chickpea germplasm encompasses wide range of diversities in terms of morphology at both vegetative and reproductive stages. There is an immediate need for understanding the genetic and molecular basis of this diversity and utilizing them for the yield contributing trait improvement. The review discusses some of the key developmental events which have potential in yield enhancement and the lessons which can be learnt from model legumes in this regard.

Comparative transcriptomics and co-expression networks reveal cultivar-specific molecular signatures associated with reproductive-stage cold stress in rice

The resistance of Huaidao5 results from the high constitutive expression of tolerance genes, while that of Huaidao9 is due to the cold-induced resistance in flag leaves and panicles. The regulation mechanism of rice seedlings' cold tolerance is relatively clear, and knowledge of its underlying mechanisms at the reproductive stage is limited. We performed differential expression and co-expression network analyses to transcriptomes from panicle and flag leaf tissues of a cold-tolerant cultivar (Huaidao5), and a sensitive cultivar (Huaidao9), under reproductive-stage cold stress. The results revealed that the expression levels of genes in stress-related pathways such as MAPK signaling pathway, diterpenoid biosynthesis, glutathione metabolism, plant-pathogen interaction and plant hormone signal transduction were constitutively highly expressed in Huaidao5, especially in panicles. Moreover, the Hudaidao5's panicle sample-specific (under cold) module contained some genes related to rice yield, such as GW5L, GGC2, SG1 and CTPS1. However, the resistance of Huaidao9 was derived from the induced resistance to cold in flag leaves and panicles. In the flag leaves, the responses included a series of stress response and signal transduction, while in the panicles nitrogen metabolism was severely affected, especially 66 endosperm-specific genes. Through integrating differential expression with co-expression networks, we predicted 161 candidate genes (79 cold-responsive genes common to both cultivars and 82 cold-tolerance genes associated with differences in cold tolerance between cultivars) potentially affecting cold response/tolerance, among which 85 (52.80%) were known to be cold-related genes. Moreover, 52 (65.82%) cold-responsive genes (e.g., TIFY11C, LSK1 and LPA) could be confirmed by previous transcriptome studies and 72 (87.80%) cold-tolerance genes (e.g., APX5, OsFbox17 and OsSTA109) were located within QTLs associated with cold tolerance. This study provides an efficient strategy for further discovery of mechanisms of cold tolerance in rice.

Transcriptome analysis reveals new insights in the starch biosynthesis of non-waxy and waxy broomcorn millet (Panicum miliaceum L.)

Broomcorn millet is a popular cereal with health benefits, and its grains are rich in starch. However, the differences in the pathway and key genes involved in starch biosynthesis of waxy and non-waxy broomcorn millet grain remain unclear. Therefore, the grain and starch physicochemical index and transcriptomic analyses of two genotypes of broomcorn millet were conducted at 3, 6, 9, 12, 15, 18, and 21 days after pollination. The phenotypic and physiological results indicated that the starch synthetic process of non-waxy and waxy broomcorn millet was significantly different. The amylose, amylopectin, and total starch contents of non-waxy broomcorn millet were 1.99, 4.74, and 6.73 mg/grain, while those of waxy broomcorn millet were 0.34, 5.94, and 6.28 mg/grain, respectively. The transcriptomic analysis revealed that 106 differentially expressed genes were identified, which were mainly enriched in the "amino sugar and nucleotide sugar metabolism", "pyruvate metabolism", "galactose metabolism", and "starch and sucrose metabolism" pathways. The WGCNA suggested that a total of 31 hub genes were correlated with starch biosynthesis. These findings provide a new approach to studying the starch synthesis in broomcorn millet.

The transcriptional regulatory network of hormones and genes under salt stress in tomato plants (Solanum lycopersicum L.)

Salt stress has become one of the main limiting factors affecting the normal growth and development of tomatoes as well as fruit quality and yields. To further reveal the regulatory relationships between tomato hormones under salt stress, the interaction between hormones and TF and the genome-wide gene interaction network were analyzed and constructed. After salt treatment, the levels of ABA, SA, and JA were significantly increased, the levels of GA were decreased, and IAA and tZ showed a trend of first increasing and then decreasing. The expression patterns of hormone biosynthesis and signal transduction related genes were analyzed based on RNA-seq analysis, the co-expression network of hormones and genome-wide co-expression networks were constructed using weighted gene co-expression network analysis (WGCNA). The expression patterns of specific transcription factors under salt stress were also systematically analyzed and identified 20 hormone-related candidate genes associated with salt stress. In conclusion, we first revealed the relationship between hormones and genes in tomatoes under salt stress based on hormone and transcriptome expression profiles and constructed a gene regulatory network. A transcriptional regulation model of tomato consisted of six types of hormones was also proposed. Our study provided valuable insights into the molecular mechanisms regulating salt tolerance in tomatoes.

The Integrated mRNA and miRNA Approach Reveals Potential Regulators of Flowering Time in Arundina graminifolia

Orchids are among the most precious flowers in the world. Regulation of flowering time is one of the most important targets to enhance their ornamental value. The beauty of Arundina graminifolia is its year-round flowering, although the molecular mechanism of this flowering ability remains masked. Therefore, we performed a comprehensive assessment to integrate transcriptome and miRNA sequencing to disentangle the genetic regulation of flowering in this valuable species. Clustering analyses provided a set of molecular regulators of floral transition and floral morphogenesis. We mined candidate floral homeotic genes, including FCA, FPA, GI, FT, FLC, AP2, SOC1, SVP, GI, TCP, and CO, which were targeted by a variety of miRNAs. MiR11091 targeted the highest number of genes, including candidate regulators of phase transition and hormonal control. The conserved miR156-miR172 pathway of floral time regulation was evident in our data, and we found important targets of these miRNAs in the transcriptome. Moreover, endogenous hormone levels were determined to decipher the hormonal control of floral buds in A. graminifolia. The qRT-PCR analysis of floral and hormonal integrators validated the transcriptome expression. Therefore, miRNA-mediated mining of candidate genes with hormonal regulation forms the basis for comprehending the complex regulatory network of perpetual flowering in precious orchids. The findings of this study can do a great deal to broaden the breeding programs for flowering time manipulation of orchids.

Weighted gene co-expression network analysis reveals key module and hub genes associated with the anthocyanin biosynthesis in maize pericarp

Anthocyanins are the visual pigments that present most of the colors in plants. Its biosynthesis requires the coordinated expression of structural genes and regulatory genes. Pericarps are the rich sources of anthocyanins in maize seeds. In the experiment, the transcriptomes of transparent and anthocyanins-enriched pericarps at 15, 20, and 25 DAP were obtained. The results output 110.007 million raw reads and 51407 genes' expression matrix. Using data filtration in R language, 2057 genes were eventually identified for weighted gene co-expression network analysis. The results showed that 2057 genes were classified into ten modules. The cyan module containing 183 genes was confirmed to be the key module with the highest correlation value of 0.98 to the anthocyanins trait. Among 183 genes, seven structural genes were mapped the flavonoid biosynthesis pathway, and a transcription factor Lc gene was annotated as an anthocyanin regulatory gene. Cluster heatmap and gene network analysis further demonstrated that Naringenin, 2-oxoglutarate 3-dioxygenase (Zm00001d001960), Dihydroflavonol 4-reductase (Zm00001d044122), Leucoanthocyanidin dioxygenase (Zm00001d014914), anthocyanin regulatory Lc gene (Zm00001d026147), and Chalcone synthase C2 (Zm00001d052673) participated in the anthocyanins biosynthesis. And the transcription factor anthocyanin regulatory Lc gene Zm00001d026147 may act on the genes Chalcone synthase C2 (Zm00001d052673) and Dihydroflavonol 4-reductase (Zm00001d044122). The yeast one-hybrid assays confirmed that the Lc protein could combine with the promoter region of C2 and directly regulate the anthocyanin biosynthesis in the pericarp. These results may provide a new sight to uncover the module and hub genes related to anthocyanins biosynthesis in plants.

An integrated transcriptome mapping the regulatory network of coding and long non-coding RNAs provides a genomics resource in chickpea

Large-scale transcriptome analysis can provide a systems-level understanding of biological processes. To accelerate functional genomic studies in chickpea, we perform a comprehensive transcriptome analysis to generate full-length transcriptome and expression atlas of protein-coding genes (PCGs) and long non-coding RNAs (lncRNAs) from 32 different tissues/organs via deep sequencing. The high-depth RNA-seq dataset reveal expression dynamics and tissue-specificity along with associated biological functions of PCGs and lncRNAs during development. The coexpression network analysis reveal modules associated with a particular tissue or a set of related tissues. The components of transcriptional regulatory networks (TRNs), including transcription factors, their cognate cis-regulatory motifs, and target PCGs/lncRNAs that determine developmental programs of different tissues/organs, are identified. Several candidate tissue-specific and abiotic stress-responsive transcripts associated with quantitative trait loci that determine important agronomic traits are also identified. These results provide an important resource to advance functional/translational genomic and genetic studies during chickpea development and environmental conditions.

A full-length transcriptome and expression atlas of protein-coding genes and long non-coding RNAs is generated in chickpea. Components of transcriptional regulatory networks and candidate tissue-specific transcripts associated with quantitative trait loci are identified.

Genome-wide transcriptome analysis reveals the regulatory network governing terpene trilactones biosynthesis in Ginkgo biloba

Ginkgo biloba L. is currently the only remaining gymnosperm of the Ginkgoaceae Ginkgo genus, and its history can be traced back to the Carboniferous 200 million years ago. Terpene trilactones (TTLs) are one of the main active ingredients in G. biloba, including ginkgolides and bilobalide. They have a good curative effect on cardiovascular and cerebrovascular diseases because of their special antagonistic effect on platelet-activating factors. Therefore, it is necessary to deeply mine genes related to TTLs and to analyze their transcriptional regulation mechanism, which will hold vitally important scientific and practical significance for quality improvement and regulation of G. biloba. In this study, we performed RNA-Seq on the root, stem, immature leaf, mature leaf, microstrobilus, ovulate strobilus, immature fruit and mature fruit of G. biloba. The TTL regulatory network of G. biloba in different organs was revealed by different transcriptomic analysis strategies. Weighted gene co-expression network analysis (WGCNA) revealed that the five modules were closely correlated with organs. The 12 transcription factors, 5 structural genes and 24 Cytochrome P450 (CYP450) were identified as candidate regulators for TTL accumulation by WGCNA and cytoscape visualization. Finally, 6 APETALA2/ethylene response factors, 2 CYP450s and bHLH were inferred to regulate the metabolism of TTLs by correlation analysis. This study is the comprehensive in authenticating transcription factors, structural genes and CYP450 involved in TTL biosynthesis, thereby providing molecular evidence for revealing the comprehensive regulatory network involved in TTL metabolism in G. biloba.

Genome-wide analysis suggests the potential role of lncRNAs during seed development and seed size/weight determination in chickpea

The integrated transcriptome data analyses suggested the plausible roles of lncRNAs during seed development in chickpea. The candidate lncRNAs associated with QTLs and those involved in miRNA-mediated seed size/weight determination in chickpea have been identified. Long non-coding RNAs (lncRNAs) are important regulators of various biological processes. Here, we identified lncRNAs at seven successive stages of seed development in small-seeded and large-seeded chickpea cultivars. In total, 4751 lncRNAs implicated in diverse biological processes were identified. Most of lncRNAs were conserved between the two cultivars, whereas only a few of them were conserved in other plants, suggesting their species-specificity. A large number of lncRNAs differentially expressed between the two chickpea cultivars associated with seed development-related processes were identified. The lncRNAs acting as precursors of miRNAs and those mimicking target protein-coding genes of miRNAs involved in seed size/weight determination, including HAIKU1, BIG SEEDS1, and SHB1, were also revealed. Further, lncRNAs located within seed size/weight associated quantitative trait loci were also detected. Overall, we present a comprehensive resource and identified candidate lncRNAs that may play important roles during seed development and seed size/weight determination in chickpea.

Comparative transcriptomics of Pinus massoniana organs provides insights on terpene biosynthesis regulation

Terpenoids are the most important natural products collected from conifer species. However, the molecular mechanisms and core factors underlying terpenoid biosynthesis in Pinus massoniana remain unclear. To clarify these mechanisms, this study aimed to identify potential genes that might participate in the terpenoid biosynthesis of P. massoniana. In this study, single molecule real-time (SMRT) sequencing and expression analysis were used to confirm the expression patterns of genes involved in the cones, immature needles, mature needles, immature branches, and mature branches of P. massoniana. A total of 31,331 lncRNAs and 71,240 mRNAs were identified from these organs, and the greatest number of differentially expressed genes (DEGs) was discovered between needles and branches. Weighted gene coexpression network analysis (WGCNA) classified all expressed genes into nine typical modules with 11 kinds of transcription factors (TFs), namely, AP2-ERF, ARF, AUX-IAA, C2H2, Dof, F-box, SBP, WRKY, bHLH, bZIP, and GRAS, and seven kinds of functional genes, namely, ABC transporter, cellulose synthase (CesA), leucine-rich repeats (LRR), cytochrome P450 (CYT P450), pathogenesis-related protein (PR), terpene synthase (TPS), and chlorophyllase enzyme. A molecular network was constructed for hub genes, TFs, and functional genes in three modules. The potential function of eight candidate genes, including PmbHLH2, PmERF1, PmRGA, PmGAI, PmbZIP1, PmLOB1, PmMADS1, and PmMYB1, was validated through correlation analysis between terpenoid contents and expression levels, subcellular localization, and transcriptional activation activity, which provides us with probable regulators of terpenoid biosynthesis in conifers.

Genetic variation in CaTIFY4b contributes to drought adaptation in chickpea

Chickpea production is vulnerable to drought stress. Identifying the genetic components underlying drought adaptation is crucial for enhancing chickpea productivity. Here, we present the fine mapping and characterization of 'QTL-hotspot', a genomic region controlling chickpea growth with positive consequences on crop production under drought. We report that a non-synonymous substitution in the transcription factor CaTIFY4b regulates seed weight and organ size in chickpea. Ectopic expression of CaTIFY4b in Medicago truncatula enhances root growth under water deficit. Our results suggest that allelic variation in 'QTL-hotspot' improves pre-anthesis water use, transpiration efficiency, root architecture and canopy development, enabling high-yield performance under terminal drought conditions. Gene expression analysis indicated that CaTIFY4b may regulate organ size under water deficit by modulating the expression of GRF-INTERACTING FACTOR1 (GIF1), a transcriptional co-activator of Growth-Regulating Factors. Taken together, our study offers new insights into the role of CaTIFY4b and on diverse physiological and molecular mechanisms underpinning chickpea growth and production under specific drought scenarios.

Transcriptomes across fertilization and seed development in the water lily Nymphaea thermarum (Nymphaeales): evidence for epigenetic patterning during reproduction

The first record of gene expression during seed development within the Nymphaeales provides evidence for a variety of biological processes, including dynamic epigenetic patterning during sexual reproduction in the water lily Nymphaea thermarum. Studies of gene expression during seed development have been performed for a growing collection of species from a phylogenetically broad sampling of flowering plants (angiosperms). However, angiosperm lineages whose origins predate the divergence of monocots and eudicots have been largely overlooked. In order to provide a new resource for understanding the early evolution of seed development in flowering plants, we sequenced transcriptomes of whole ovules and seeds from three key stages of reproductive development in the waterlily Nymphaea thermarum, an experimentally tractable member of the Nymphaeales. We first explore patterns of gene expression, beginning with mature ovules and continuing through fertilization into early- and mid-stages of seed development. We find patterns of gene expression that corroborate histological/morphological observations of seed development in this species, such as expression of genes involved in starch synthesis and transcription factors that have been associated with embryo and endosperm development in other species. We also find evidence for processes that were previously not known to be occurring during seed development in this species, such as epigenetic modification. We then examine the expression of genes associated with patterning DNA and histone methylation-processes that are essential for seed development in distantly related and structurally diverse monocots and eudicots. Around 89% of transcripts putatively homologous to DNA and histone methylation modifiers are expressed during seed development in N. thermarum, including homologs of genes known to pattern imprinting-related epigenetic modifications. Our results suggest that dynamic epigenetic patterning is a deeply conserved aspect of angiosperm seed development.

Analysis of transcriptome data and quantitative trait loci enables the identification of candidate genes responsible for fiber strength in Gossypium barbadense

Gossypium barbadense possesses a superior fiber quality because of its fiber length and strength. An in-depth analysis of the underlying genetic mechanism could aid in filling the gap in research regarding fiber strength and could provide helpful information for Gossypium barbadense breeding. Three quantitative trait loci related to fiber strength were identified from a Gossypium barbadense recombinant inbred line (PimaS-7 × 5917) for further analysis. RNA sequencing was performed in the fiber tissues of PimaS-7 × 5917 0–35 days postanthesis. Four specific modules closely related to the secondary wall-thickening stage were obtained using the weighted gene coexpression network analysis. In total, 55 genes were identified as differentially expressed from 4 specific modules. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes were used for enrichment analysis, and Gbar_D11G032910, Gbar_D08G020540, Gbar_D08G013370, Gbar_D11G033670, and Gbar_D11G029020 were found to regulate fiber strength by playing a role in the composition of structural constituents of cytoskeleton and microtubules during fiber development. Quantitative real-time PCR results confirmed the accuracy of the transcriptome data. This study provides a quick strategy for exploring candidate genes and provides new insights for improving fiber strength in cotton.

Dynamic transcriptome analysis suggests the key genes regulating seed development and filling in Tartary buckwheat (Fagopyrum tataricum Garetn.)

Tartary buckwheat is highly attractive for the richness of nutrients and quality, yet post-embryonic seed abortion greatly halts the yield. Seed development is crucial for determining grain yield, whereas the molecular basis and regulatory network of Tartary buckwheat seed development and filling is not well understood at present. Here, we assessed the transcriptional dynamics of filling stage Tartary buckwheat seeds at three developmental stages by RNA sequencing. Among the 4249 differentially expressed genes (DEGs), genes related to seed development were identified. Specifically, 88 phytohormone biosynthesis signaling genes, 309 TFs, and 16 expansin genes participating in cell enlargement, 37 structural genes involved in starch biosynthesis represented significant variation and were candidate key seed development genes. Cis-element enrichment analysis indicated that the promoters of differentially expressed expansin genes and starch biosynthesis genes are rich of hormone-responsive (ABA-, AUX-, ET-, and JA-), and seed growth-related (MYB, MYC and WRKY) binding sites. The expansin DEGs showed strong correlations with DEGs in phytohormone pathways and transcription factors (TFs). In total, phytohormone ABA, AUX, ET, BR and CTK, and related TFs could substantially regulate seed development in Tartary buckwheat through targeting downstream expansin genes and structural starch biosynthetic genes. This transcriptome data could provide a theoretical basis for improving yield of Tartary buckwheat.

Lysophosphatidic acid acyltransferase 2 and 5 commonly, but differently, promote seed oil accumulation in Brassica napus

Background

Increasing seed oil content (SOC) of Brassica napus has become one of the main plant breeding goals over the past decades. Lysophosphatidic acid acyltransferase (LPAT) performs an important molecular function by regulating the production of phosphatidic acid (PA), a key intermediate in the synthesis of membrane and storage lipids. However, the mechanism underlying the effect of LPAT on the SOC of B. napus remains unclear.

Results

In the present study, significant elevation of SOC was achieved by overexpressing BnLPAT2 and BnLPAT5 in B. napus. RNAi and CRISPR–Cas9 were also successfully used to knock down and knock out these two genes in B. napus where SOC significantly decreased. Meanwhile, we found an accumulation of lipid droplets and oil bodies in seeds of BnLPAT2 and BnLPAT5 overexpression lines, whereas an increase of sugar and protein in Bnlpat2 and Bnlpat5 mutant seeds. Sequential transcriptome analysis was further performed on the developing seeds of the BnLPAT2 and BnLPAT5 overexpression, knockdown, and knockout rapeseed lines. Most differentially expressed genes (DEGs) that were expressed in the middle and late stages of seed development were enriched in photosynthesis and lipid metabolism, respectively. The DEGs involved in fatty acid and lipid biosynthesis were active in the overexpression lines but were relatively inactive in the knockdown and knockout lines. Further analysis revealed that the biological pathways related to fatty acid/lipid anabolism and carbohydrate metabolism were specifically enriched in the BnLPAT2 overexpression lines.

Conclusions

BnLPAT2 and BnLPAT5 are essential for seed oil accumulation. BnLPAT2 preferentially promoted diacylglycerol synthesis to increase SOC, whereas BnLPAT5 tended to boost PA synthesis for membrane lipid generation. Taken together, BnLPAT2 and BnLPAT5 can jointly but differently promote seed oil accumulation in B. napus. This study provides new insights into the potential mechanisms governing the promotion of SOC by BnLPAT2 and BnLPAT5 in the seeds of B. napus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02182-2.

Transcriptome analysis reveals cell cycle-related transcripts as key determinants of varietal differences in seed size of Brassica juncea

Brassica juncea is an important oilseed crop, widely grown as a source of edible oil. Seed size is a pivotal agricultural trait in oilseed Brassicas. However, the regulatory mechanisms underlying seed size determination are poorly understood. To elucidate the transcriptional dynamics involved in the determination of seed size in B. juncea, we performed a comparative transcriptomic analysis using developing seeds of two varieties, small-seeded Early Heera2 (EH2) and bold-seeded Pusajaikisan (PJK), at three distinct stages (15, 30 and 45 days after pollination). We detected 112,550 transcripts, of which 27,186 and 19,522 were differentially expressed in the intra-variety comparisons and inter-variety comparisons, respectively. Functional analysis using pathway, gene ontology, and transcription factor enrichment revealed that cell cycle- and cell division-related transcripts stay upregulated during later stages of seed development in the bold-seeded variety but are downregulated at the same stage in the small-seeded variety, indicating that an extended period of cell proliferation in the later stages increased seed weight in PJK as compared to EH2. Further, k-means clustering and candidate genes-based analyses unravelled candidates for employing in seed size improvement of B. juncea. In addition, candidates involved in determining seed coat color, oil content, and other seed traits were also identified.

Transcriptome Profiling of a Common Mistletoe Species Parasitizing Four Typical Host Species in Urban Southwest China

Comparing gene expressions among parasitic plants infecting different host species can have significant implications for understanding host–parasite interactions. Taxillus nigrans is a common hemiparasitic species in Southwest China that parasitizes a variety of host species. However, a lack of nucleotide sequence data to date has hindered transcriptome-level research on T. nigrans. In this study, the transcriptomes of T. nigrans individuals parasitizing four typical host species (Broussonetia papyrifera (Bpap), a broad-leaved tree species; Cryptomeria fortunei (Cfor), a coniferous tree species; Cinnamomum septentrionale (Csep), an evergreen tree species; and Ginkgo biloba (Gbil), a deciduous-coniferous tree species) were sequenced, and the expression profiles and metabolic pathways were compared among hosts. A total of greater than 400 million reads were generated in nine cDNA libraries. These were de novo assembled into 293823 transcripts with an N50 value of 1790 bp. A large number of differentially expressed genes (DEGs) were identified when comparing T. nigrans individuals on different host species: Bpap vs. Cfor (1253 DEGs), Bpap vs. Csep (864), Bpap vs. Gbil (517), Cfor vs. Csep (259), Cfor vs. Gbil (95), and Csep vs. Gbil (40). Four hundred and fifteen unigenes were common to all six pairwise comparisons; these were primarily associated with Cytochrome P450 and environmental adaptation, as determined in a KEGG enrichment analysis. Unique unigenes were also identified, specific to Bpap vs. Cfor (808 unigenes), Bpap vs. Csep (329 unigenes), Bpap vs. Gbil (87 unigenes), Cfor vs. Csep (108 unigenes), Cfor vs. Gbil (32 unigenes), and Csep vs. Gbil comparisons (23 unigenes); partial unigenes were associated with the metabolism of terpenoids and polyketides regarding plant hormone signal transduction. Weighted gene co-expression network analysis (WGCNA) revealed four modules that were associated with the hosts. These results provide a foundation for further exploration of the detailed molecular mechanisms involved in plant parasitism.

Draft genome sequence of Indian mulberry (Morus indica) provides a resource for functional and translational genomics

Mulberry is an important crop plant for the sericulture industry. Here, we report high-quality genome sequence of a cultivated Indian mulberry (Morus indica cv K2) obtained by combining data from four different technologies, including Illumina, single-molecule real-time sequencing, chromosome conformation capture and optical mapping, with a gene completeness of 96.5%. Based on the genome sequence, we identified 49.2% of repetitive DNA and 27,435 high-confidence protein-coding genes with >90% of them supported by transcript evidence. A comparative analysis with other plant genomes identified 4.8% of species-specific genes in the M. indica genome. Transcriptome profiling revealed tissue-specific and differential expression across multiple accessions of ~4.7% and 2-5% of protein-coding genes, respectively, implicated in diverse biological processes. Whole genome resequencing of 21 accessions/species revealed ~2.5 million single nucleotide polymorphisms and ~ 0.2 million insertions/deletions. These data and results provide a comprehensive resource to accelerate the genomics research in mulberry for its improvement.

Transcriptome and Coexpression Network Analyses Provide In-Sights into the Molecular Mechanisms of Hydrogen Cyanide Synthesis during Seed Development in Common Vetch (Vicia sativa L.)

The common vetch (Vicia sativa L.) seed is an ideal plant-based protein food for humans, but its edible value is mainly limited by the presence of cyanogenic glycosides that hydrolyze to produce toxic hydrogen cyanide (HCN), and the genes that regulate HCN synthesis in common vetch are unknown. In this study, seeds from common vetch at 5, 10, 15, 20, 25, 30, and 35 days after anthesis were sampled, and the seven stages were further divided into five developmental stages, S1, S2, S3, S4, and S5, based on morphological and transcriptome analyses. A total of 16,403 differentially expressed genes were identified in the five developmental stages. The HCN contents of seeds in these five stages were determined by alkaline titration, and weighted gene coexpression network analysis was used to explain the molecular regulatory mechanism of HCN synthesis in common vetch seeds. Eighteen key regulatory genes for HCN synthesis were identified, including the VsGT2, VsGT17 and CYP71A genes, as well as the VsGT1 gene family. VsGT1, VsGT2, VsGT17 and CYP71A jointly promoted HCN synthesis, from 5 to 25 days after anthesis, with VsGT1-1, VsGT1-4, VsGT1-11 and VsGT1-14 playing major roles. The HCN synthesis was mainly regulated by VsGT1, from 25 to 35 days after anthesis. As the expression level of VsGT1 decreased, the HCN content no longer increased. In-depth elucidation of seed HCN synthesis lays the foundations for breeding common vetch with low HCN content.

Comparative Flower Transcriptome Network Analysis Reveals DEGs Involved in Chickpea Reproductive Success during Salinity

Salinity is increasingly becoming a significant problem for the most important yet intrinsically salt-sensitive grain legume chickpea. Chickpea is extremely sensitive to salinity during the reproductive phase. Therefore, it is essential to understand the molecular mechanisms by comparing the transcriptomic dynamics between the two contrasting genotypes in response to salt stress. Chickpea exhibits considerable genetic variation amongst improved cultivars, which show better yields in saline conditions but still need to be enhanced for sustainable crop production. Based on previous extensive multi-location physiological screening, two identified genotypes, JG11 (salt-tolerant) and ICCV2 (salt-sensitive), were subjected to salt stress to evaluate their phenological and transcriptional responses. RNA-Sequencing is a revolutionary tool that allows for comprehensive transcriptome profiling to identify genes and alleles associated with stress tolerance and sensitivity. After the first flowering, the whole flower from stress-tolerant and sensitive genotypes was collected. A total of ~300 million RNA-Seq reads were sequenced, resulting in 2022 differentially expressed genes (DEGs) in response to salt stress. Genes involved in flowering time such as FLOWERING LOCUS T (FT) and pollen development such as ABORTED MICROSPORES (AMS), rho-GTPase, and pollen-receptor kinase were significantly differentially regulated, suggesting their role in salt tolerance. In addition to this, we identify a suite of essential genes such as MYB proteins, MADS-box, and chloride ion channel genes, which are crucial regulators of transcriptional responses to salinity tolerance. The gene set enrichment analysis and functional annotation of these genes in flower development suggest that they can be potential candidates for chickpea crop improvement for salt tolerance.

Genetic insights into the regulatory pathways for continuous flowering in a unique orchid Arundina graminifolia

BACKGROUND

Manipulation of flowering time and frequency of blooming is key to enhancing the ornamental value of orchids. Arundina graminifolia is a unique orchid that flowers year round, although the molecular basis of this flowering pattern remains poorly understood.

RESULTS

We compared the A. graminifolia transcriptome across tissue types and floral developmental stages to elucidate important genetic regulators of flowering and hormones. Clustering analyses identified modules specific to floral transition and floral morphogenesis, providing a set of candidate regulators for the floral initiation and timing. Among candidate floral homeotic genes, the expression of two FT genes was positively correlated with flower development. Assessment of the endogenous hormone levels and qRT-PCR analysis of 32 pathway-responsive genes supported a role for the regulatory networks in floral bud control in A. graminifolia. Moreover, WGCNA showed that flowering control can be delineated by modules of coexpressed genes; especially, MEgreen presented group of genes specific to flowering.

CONCLUSIONS

Candidate gene selection coupled with hormonal regulators brings a robust source to understand the intricate molecular regulation of flowering in precious orchids.

Draft genome and transcriptome analyses of halophyte rice Oryza coarctata provide resources for salinity and submergence stress response factors

Oryza coarctata is a wild relative of rice that has adapted to diverse ecological environments, including high salinity and submergence. Thus, it can provide an important resource for discovering candidate genes/factors involved in tolerance to these stresses. Here, we report a draft genome assembly of 573 Mb comprised of 8877 scaffolds with N50 length of 205 kb. We predicted a total of 50,562 protein-coding genes, of which a significant fraction was found to be involved in secondary metabolite biosynthesis and hormone signal transduction pathways. Several salinity and submergence stress-responsive protein-coding and long noncoding RNAs involved in diverse biological processes were identified using RNA-sequencing data. Based on small RNA sequencing, we identified 168 unique miRNAs and 3219 target transcripts (coding and noncoding) involved in several biological processes, including abiotic stress responses. Further, whole genome bisulphite sequencing data analysis revealed at least 19%-48% methylcytosines in different sequence contexts and the influence of methylation status on gene expression. The genome assembly along with other datasets have been made publicly available at http://ccbb.jnu.ac.in/ory-coar. Altogether, we provide a comprehensive genomic resource for understanding the regulation of salinity and submergence stress responses and identification of candidate genes/factors involved for functional genomics studies.

Integrative Morphological, Physiological, Proteomics Analyses of Jujube Fruit Development Provide Insights Into Fruit Quality Domestication From Wild Jujube to Cultivated Jujube

Jujube (Ziziphus jujuba) was domesticated from wild jujube (Z. jujuba var. spinosa). Here, integrative physiological, metabolomic, and comparative proteomic analyses were performed to investigate the fruit expansion and fruit taste components in a jujube cultivar 'Junzao' and a wild jujube 'Qingjiansuanzao' with contrasting fruit size and taste. We revealed that the duration of cell division and expansion largely determined the final fruit size, while the intercellular space in the mesocarp dictated the ratio of mesocarp volume in mature fruits. The high levels of endogenous gibbereline3 (GA) and zeatin in the growing fruit of 'Junzao' were associated with their increased fruit expansion. Compared with 'Junzao,' wild jujube accumulated lower sugars and higher organic acids. Furthermore, several protein co-expression modules and important member proteins correlated with fruit expansion, sugar synthesis, and ascorbic acid metabolism were identified. Among them, GA20OX involved in GA biosynthesis was identified as a key protein regulating fruit expansion, whereas sucrose-6-phosphate synthase (SPS) and neutral invertase (NINV) were considered as key enzymes promoting sugar accumulation and as major factors regulating the ratio of sucrose to hexose in jujube fruits, respectively. Moreover, the increase of Nicotinamide adenine dinucleotide-Malate dehydrogenase (NAD-MDH) activity and protein abundance were associated with the malic acid accumulation, and the high accumulation of ascorbic acid in wild jujube was correlated with the elevated abundance of GalDH, ZjAPXs, and MDHAR1, which are involved in the ascorbic acid biosynthesis and recycling pathways. Overall, these results deepened the understanding of mechanisms regulating fruit expansion and sugar/acids metabolisms in jujube fruit.

Integrated Transcriptomic and Bioinformatics Analyses Reveal the Molecular Mechanisms for the Differences in Seed Oil and Starch Content Between Glycine max and Cicer arietinum

The seed oil and starch content of soybean are significantly different from that of chickpea. However, there are limited studies on its molecular mechanisms. To address this issue, we conducted integrated transcriptomic and bioinformatics analyses for species-specific genes and acyl-lipid-, starch-, and carbon metabolism-related genes. Among seven expressional patterns of soybean-specific genes, four were highly expressed at the middle- and late oil accumulation stages; these genes significantly enriched fatty acid synthesis and carbon metabolism, and along with common acetyl CoA carboxylase (ACCase) highly expressed at soybean middle seed development stage, common starch-degrading enzyme beta-amylase-5 (BAM5) was highly expressed at soybean early seed development stage and oil synthesis-related genes ACCase, KAS, KAR, ACP, and long-chain acyl-CoA synthetase (LACS) were co-expressed with WRI1, which may result in high seed oil content and low seed starch content in soybean. The common ADP-glucose pyrophosphorylase (AGPase) was highly expressed at chickpea middle seed development stage, along with more starch biosynthesis genes co-expressed with four-transcription-factor homologous genes in chickpea than in soybean, and the common WRI1 was not co-expressed with oil synthesis genes in chickpea, which may result in high seed starch content and low seed oil content in chickpea. The above results may be used to improve chickpea seed oil content in two ways. One is to edit CaWRI1 to co-express with oil synthesis-related genes, which may increase carbon metabolites flowing to oil synthesis, and another is to increase the expression levels of miRNA159 and miRNA319 to inhibit the expression of MYB33, which may downregulate starch synthesis-related genes, making more carbon metabolites flow into oil synthesis. Our study will provide a basis for future breeding efforts to increase the oil content of chickpea seeds.

Comparative Dynamic Transcriptome Reveals the Delayed Secondary-Cell-Wall Thickening Results in Altered Lint Percentage and Fiber Elongation in a Chromosomal Segment Substitution Line of Cotton (Gossypium hirsutum L.)

Lint percentage (LP) is an important yield component in cotton that is usually affected by initial fiber number and cell wall thickness. To explore how fiber cell wall development affects LP, phenotypic identification and dynamic transcriptome analysis were conducted using a single segment substitution line of chromosome 15 (SL15) that harbors a major quantitative trait locus (QTL) for LP. Compared to its recurrent parent LMY22, SL15 did not differ in initial fiber number, but the fiber cell wall thickness and single-fiber weight decreased significantly, altering LP. The comparative transcriptome profiles revealed that the secondary cell wall (SCW) development phase of SL15 was relatively delayed. Meanwhile, the expression of genes related to cell expansion decreased more slightly in SL15 with fiber development, resulting in relatively higher expression at SL15_25D than at LMY22_25D. SCW development-related genes, such as GhNACs and GhMYBs, in the putative NAC-MYB-CESA network differentially expressed at SL15_25D, along with the lower expression of CESA6, CSLC12, and CSLA2. The substituted chromosomal interval was further investigated, and found 6 of 146 candidate genes were differentially expressed in all four cell development periods including 10, 15, 20 and 25 DPA. Genetic variation and co-expression analysis showed that GH_D01G0052, GH_D01G0099, GH_D01G0100, and GH_D01G0140 may be important candidate genes associated with qLP-C15-1. Our results provide novel insights into cell wall development and its relationship with LP, which is beneficial for lint yield and fiber quality improvement.

Proso millet (Panicum miliaceum L.): A potential crop to meet demand scenario for sustainable saline agriculture

Proso millet (Panicum miliaceum L.) is resilient to abiotic stress, especially to land degradation caused by soil salinization. However, the mechanisms by which its roots adapt and tolerate salt stress are obscure. In this study, plants of a salt-sensitive cultivar (SS 212) and a salt-tolerant cultivar (ST 47) of proso millet were exposed to severe salt stress and subsequent re-watering. ST 47 exhibited greater salt tolerance than SS 212, as evidenced by higher increases in total root length (TRL), root surface area (RSA), root tip number (RTN). Moreover, microstructural analysis showed that relative to SS 212, the roots of ST 47 could maintain more intact internal structures and thicker cell walls under salt stress. Digital RNA sequence analysis revealed that ST 47 maintained better Na⁺/K⁺ balance to resist Na⁺ toxicity via a higher capability to restrict Na⁺ uptake, vacuolar Na⁺ sequestration, and Na⁺ exclusion. The mechanism for Na⁺ toxicity resistance in ST 47 involved promoting cell wall composition changes via efficient regulation of galactose metabolism and biosynthesis of cellulose and phenylpropanoids. Overall, this study provides valuable salt-tolerant cultivar resources and mechanisms for regulating salt tolerance, which could be applied for the rehabilitation of saline lands.

Genome-wide identification of seed storage protein gene regulators in wheat through coexpression analysis

Only a few transcriptional regulators of seed storage protein (SSP) genes have been identified in common wheat (Triticum aestivum L.). Coexpression analysis could be an efficient approach to characterize novel transcriptional regulators at the genome-scale considering the correlated expression between transcriptional regulators and target genes. As the A genome donor of common wheat, Triticum urartu is more suitable for coexpression analysis than common wheat considering the diploid genome and single gene copy. In this work, the transcriptome dynamics in endosperm of T. urartu throughout grain filling were revealed by RNA-Seq analysis. In the coexpression analysis, a total of 71 transcription factors (TFs) from 23 families were found to be coexpressed with SSP genes. Among these TFs, TuNAC77 enhanced the transcription of SSP genes by binding to cis-elements distributed in promoters. The homolog of TuNAC77 in common wheat, TaNAC77, shared an identical function, and the total SSPs were reduced by about 24% in common wheat when TaNAC77 was knocked down. This is the first genome-wide identification of transcriptional regulators of SSP genes in wheat, and the newly characterized transcriptional regulators will undoubtedly expand our knowledge of the transcriptional regulation of SSP synthesis.

Stage Specificity, the Dynamic Regulators and the Unique Orchid Arundina graminifolia

Orchids take years to reach flowering, but the unique bamboo orchid (Arundina graminifolia) achieves reproductive maturity in six months and then keeps on year round flowering. Therefore, studying different aspects of its growth, development and flowering is key to boost breeding programs for orchids. This study uses transcriptome tools to discuss genetic regulation in five stages of flower development and four tissue types. Stage specificity was focused to distinguish genes specifically expressed in different stages of flower development and tissue types. The top 10 highly expressed genes suggested unique regulatory patterns for each stage or tissue. The A. graminifolia sequences were blasted in Arabidopsis genome to validate stage specific genes and to predict important hormonal and cell regulators. Moreover, weighted gene co-expression network analysis (WGCNA) modules were ascertained to suggest highly influential hubs for early and late stages of flower development, leaf and root. Hormonal regulators were abundant in all data sets, such as auxin (LAX2, GH3.1 and SAUR41), cytokinin (LOG1), gibberellin (GASA3 and YAB4), abscisic acid (DPBF3) and sucrose (SWEET4 and SWEET13). Findings of this study, thus, give a fine sketch of genetic variability in Orchidaceae and broaden our understanding of orchid flower development and the involvement of multiple pathways.

Comparative transcriptome analysis reveals regulatory network and regulators associated with proanthocyanidin accumulation in persimmon

BACKGROUND

Proanthocyanidins (PAs) are important plant secondary metabolites that confer flavor, nutritional value, and resistance to pathogens. Persimmon is one of the PA richest crops. Mature fruits can be inedible because of the astringency caused by high PA levels and need to go through a de-astringency treatment before consumption. The molecular basis for PA accumulation is poorly known, particularly transcriptional regulators. We characterised three genotypes ('Luotiantianshi' (LT), 'Mopanshi' (MP), and 'Youhou' (YH)) with different PA accumulation patterns using an approach that combined PacBio full-length sequencing and Illumina-based RNA sequencing to build high-quality full-length transcriptomes. Additionally, we analysed transcriptome dynamics of the three genotypes (LT, MP, and YH) at four key fruit developmental stages.

RESULTS

A total of 96,463 transcripts were obtained. We identified 80,075 protein-coding sequences (CDSs), 71,137 simple sequence repeats (SSRs), and 27,845 long noncoding RNAs (lncRNAs). Pearson correlation coefficient (PCC), principal component analysis (PCA), and differentially expressed transcripts (DETs) analyses indicated that the four different developmental stages within a genotype exhibited similar transcriptome activities. A total of 2,164 transcripts specific to each fruit developmental stage were detected. The transcripts specific to early stages were attributed to phenylpropanoid and flavonoid biosynthesis. Co-expression network analyses revealed MEbrown and MEblue modules were strongly associated to PA accumulation. From these two modules, 20 hub TFs are potential regulators for PA accumulation. Among them, Cluster_78388 (SBP protein), Cluster_63454 (bZIP protein), and Cluster_66595 (MYB protein) appear to involve in the PA biosynthesis in Chinese genotypes.

CONCLUSIONS

This is the first high-quality reference transcriptome for commercial persimmon. Our work provides insights into the molecular pathways underlying PA accumulation and enhances our global understanding of transcriptome dynamics throughout fruit development.

Key Regulators of Sucrose Metabolism Identified through Comprehensive Comparative Transcriptome Analysis in Peanuts

Sucrose content is a crucial indicator of quality and flavor in peanut seed, and there is a lack of clarity on the molecular basis of sucrose metabolism in peanut seed. In this context, we performed a comprehensive comparative transcriptome study on the samples collected at seven seed development stages between a high-sucrose content variety (ICG 12625) and a low-sucrose content variety (Zhonghua 10). The transcriptome analysis identified a total of 8334 genes exhibiting significantly different abundances between the high- and low-sucrose varieties. We identified 28 differentially expressed genes (DEGs) involved in sucrose metabolism in peanut and 12 of these encoded sugars will eventually be exported transporters (SWEETs). The remaining 16 genes encoded enzymes, such as cell wall invertase (CWIN), vacuolar invertase (VIN), cytoplasmic invertase (CIN), cytosolic fructose-bisphosphate aldolase (FBA), cytosolic fructose-1,6-bisphosphate phosphatase (FBP), sucrose synthase (SUS), cytosolic phosphoglucose isomerase (PGI), hexokinase (HK), and sucrose-phosphate phosphatase (SPP). The weighted gene co-expression network analysis (WGCNA) identified seven genes encoding key enzymes (CIN, FBA, FBP, HK, and SPP), three SWEET genes, and 90 transcription factors (TFs) showing a high correlation with sucrose content. Furthermore, upon validation, six of these genes were successfully verified as exhibiting higher expression in high-sucrose recombinant inbred lines (RILs). Our study suggested the key roles of the high expression of SWEETs and enzymes in sucrose synthesis making the genotype ICG 12625 sucrose-rich. This study also provided insights into the molecular basis of sucrose metabolism during seed development and facilitated exploring key candidate genes and molecular breeding for sucrose content in peanuts.

Integrating transcriptome, co-expression and QTL-seq analysis reveals that primary root growth in maize is regulated via flavonoid biosynthesis and auxin signal transduction

The primary root is critical for early seedling growth and survival. To understand the molecular mechanisms governing primary root development, we performed a dynamic transcriptome analysis of two maize (Zea mays) inbred lines with contrasting primary root length at nine time points over a 12-day period. A total of 18 702 genes were differentially expressed between two lines or different time points. Gene enrichment, phytohormone content determination, and metabolomics analysis showed that auxin biosynthesis and signal transduction, as well as the phenylpropanoid and flavonoid biosynthesis pathways, were associated with root development. Co-expression network analysis revealed that eight modules were associated with lines/stages, as well as primary or lateral root length. In root-related modules, flavonoid metabolism accompanied by auxin biosynthesis and signal transduction constituted a complex gene regulatory network during primary root development. Two candidate genes (rootless concerning crown and seminal roots, rtcs and Zm00001d012781) involved in auxin signaling and flavonoid biosynthesis were identified by co-expression network analysis, QTL-seq and functional annotation. These results increase our understanding of the regulatory network controlling the development of primary and lateral root length, and provide a valuable genetic resource for improvement of root performance in maize.

Transcriptome profiling reveals the spatial-temporal dynamics of gene expression essential for soybean seed development

BACKGROUND

Seeds are the economic basis of oilseed crops, especially soybeans, the most widely cultivated oilseed crop worldwide. Seed development is accompanied by a multitude of diverse cellular processes, and revealing the underlying regulatory activities is critical for seed improvement.

RESULTS

In this study, we profiled the transcriptomes of developing seeds at 20, 25, 30, and 40 days after flowering (DAF), as these stages represent critical time points of seed development from early to full development. We identified a set of highly abundant genes and highlighted the importance of these genes in supporting nutrient accumulation and transcriptional regulation for seed development. We identified 8925 differentially expressed genes (DEGs) that exhibited temporal expression patterns over the course and expression specificities in distinct tissues, including seeds and nonseed tissues (roots, stems, and leaves). Genes specific to nonseed tissues might have tissue-associated roles, with relatively low transcript abundance in developing seeds, suggesting their spatially supportive roles in seed development. Coexpression network analysis identified several underexplored genes in soybeans that bridge tissue-specific gene modules.

CONCLUSIONS

Our study provides a global view of gene activities and biological processes critical for seed formation in soybeans and prioritizes a set of genes for further study. The results of this study help to elucidate the mechanism controlling seed development and storage reserves.

Transcriptome analysis reveals the mechanism associated with dynamic changes in fatty acid and phytosterol content in foxtail millet (Setaria italica) during seed development

Foxtail millet (Setaria italica) is an excellent source of beneficial natural fatty acids and phytosterols. However, the mechanisms underlying the dynamic changes of fatty acids and phytosterols during seed development are unknown. In this study, a comprehensive dynamic change analysis of the bioactive compounds during seed development was conducted in two cultivars with different crude fat content (high-fat, JG 35 [5.40%]; and low-fat, JG 39 [2.90%]). GC-FID/MS analysis showed that the proportion of unsaturated fatty acids (UFAs) were higher than the saturated fatty acids (SFAs). UFA content first increased, then decreased during seed development, while SFA content showed the opposite trend. Oil contents continuously increased with seed development, especially at the S2 stage. Phytosterol contents initially increased, then decreased with seed development. Transcriptome analysis revealed that 152 genes were associated with fatty acid metabolism and phytosterol biosynthesis, of which 46 and 62 were related to UFA and phytosterol biosynthesis, respectively. Furthermore, the key genes involved in fatty acid synthesis (ACCase and FATA/B), triacylglycerol biosynthesis (LACS, GPAT, and DGAT), and phytosterols synthesis (CAS1, STM1, EGR6, and DWF1) were overexpressed. This led to maximum UFA, oil, and phytosterol accumulation in JG 35 at the S2 stage. This study reveals the mechanism behind the dynamic changes of fatty acid and phytosterol contents in foxtail millet during seed development.