Transcriptome analysis of Brassica napus pod using RNA-Seq and identification of lipid-related candidate genes

Fine mapping and candidate gene analysis of a major QTL for oil content in the seed of Brassica napus

KEY MESSAGE

By integrating QTL fine mapping and transcriptomics, a candidate gene responsible for oil content in rapeseed was identified. The gene is anticipated to primarily function in photosynthesis and photosystem metabolism pathways. Brassica napus is one of the most important oil crops in the world, and enhancing seed oil content is an important goal in its genetic improvement. However, the underlying genetic basis for the important trait remains poorly understood in this crop. We previously identified a major locus, OILA5 responsible for seed oil content on chromosome A5 through genome-wide association study. To better understand the genetics of the QTL, we performed fine mapping of OILA5 with a double haploid population and a BC3F2 segregation population consisting of 6227 individuals. We narrowed down the QTL to an approximate 43 kb region with twelve annotated genes, flanked by markers ZDM389 and ZDM337. To unveil the potential candidate gene responsible for OILA5, we integrated fine mapping data with transcriptome profiling using high and low oil content near-isogenic lines. Among the candidate genes, BnaA05G0439400ZS was identified with high expression levels in both seed and silique tissues. This gene exhibited homology with AT3G09840 in Arabidopsis that was annotated as cell division cycle 48. We designed a site-specific marker based on resequencing data and confirmed its effectiveness in both natural and segregating populations. Our comprehensive results provide valuable genetic information not only enhancing our understanding of the genetic control of seed oil content but also novel germplasm for advancing high seed oil content breeding in B. napus and other oil crops.

Global transcriptome analysis reveals potential genes associated with genic male sterility of rapeseed (Brassica napus L.)

Rapeseed is the third leading source of edible oil in the world. Genic male sterility (GMS) lines provide crucial material for harnessing heterosis for rapeseed. GMS lines have been widely used successfully for rapeseed hybrid production. The physiological and molecular mechanism of pollen development in GMS lines of rapeseed (Brassica napus L.) need to be determined for the creation of hybrids and cultivation of new varieties. However, limited studies have focused on systematically mining genes that regulate the pollen development of GMS lines in B. napus. In the present study, to determine the stage at which pollen development begins to show abnormality in the GMS lines, we performed semi-thin section analysis of the anthers with five pollen development stages. The results indicated that the abnormal pollen development in DGMS lines might start at the meiotic stage, and abnormal pollen development in RGMS lines probably occurred before the tetrad stage. To investigate the critical genes and pathways involved in pollen development in GMS lines, we constructed and sequenced 24 transcriptome libraries for the flower buds from the fertile and sterile lines of two recessive GMS (RGMS) lines (6251AB and 6284AB) and two dominant GMS (DGMS) lines (4001AB and 4006AB). A total of 23,554 redundant DEGs with over two-fold change between sterile and fertile lines were obtained. A total of 346 DEGs were specifically related to DGMS, while 1,553 DEGs were specifically related to RGMS. A total of 1,545 DEGs were shared between DGMS and RGMS. And 253 transcription factors were found to be differentially expressed between the sterile and fertile lines of GMS. In addition, 6,099 DEGs possibly related to anther, pollen, and microspore development processes were identified. Many of these genes have been reported to be involved in anther and microspore developmental processes. Several DEGs were speculated to be key genes involved in the regulation of fertility. Three differentially expressed genes were randomly selected and their expression levels were verified by quantitative PCR (qRT-PCR). The results of qRT-PCR largely agreed with the transcriptome sequencing results. Our findings provide a global view of genes that are potentially involved in GMS occurrence. The expression profiles and function analysis of these DEGs were provided to expand our understanding of the complex molecular mechanism in pollen and sterility development in B. napus.

Mechanism and Utilization of Ogura Cytoplasmic Male Sterility in Cruciferae Crops

Hybrid production using lines with cytoplasmic male sterility (CMS) has become an important way to utilize heterosis in vegetables. Ogura CMS, with the advantages of complete pollen abortion, ease of transfer and a progeny sterility rate reaching 100%, is widely used in cruciferous crop breeding. The mapping, cloning, mechanism and application of Ogura CMS and fertility restorer genes in Brassica napus, Brassica rapa, Brassica oleracea and other cruciferous crops are reviewed herein, and the existing problems and future research directions in the application of Ogura CMS are discussed.

Spatio-temporal transcriptome profiling and subgenome analysis in Brassica napus

Brassica napus is an important oil crop and an allotetraploid species. However, the detailed analysis of gene function and homoeologous gene expression in all tissues at different developmental stages was not explored. In this study, we performed a global transcriptome analysis of 24 vegetative and reproductive tissues at six developmental stages (totally 111 tissues). These samples were clustered into eight groups. The gene functions of silique pericarp were similar to roots, stems and leaves. In particular, glucosinolate metabolic process was associated with root and silique pericarp. Genes involved in protein phosphorylation were often associated with stamen, anther and the early developmental stage of seeds. Transcription factor (TF) genes were more specific than structural genes. A total of 17 100 genes that were preferentially expressed in one tissue (tissue-preferred genes, TPGs), including 889 TFs (5.2%), were identified in the 24 tissues. Some TPGs were identified as hub genes in the co-expression network analysis, and some TPGs in different tissues were involved in different hormone pathways. About 67.0% of the homoeologs showed balanced expression, whereas biased expression of homoeologs was associated with structural divergence. In addition, the spatiotemporal expression of homoeologs was related to the presence of transposable elements (TEs) and regulatory elements (REs); more TEs and fewer REs in the promoters resulted in divergent expression in different tissues. This study provides a valuable transcriptional map for understanding the growth and development of B. napus, for identifying important genes for future crop improvement, and for exploring gene expression patterns in the B. napus.

Comparative Transcriptome Analysis of the Anthers from the Cytoplasmic Male-Sterile Pepper Line HZ1A and Its Maintainer Line HZ1B

Cytoplasmic male-sterility (CMS) is important for the utilization of crop heterosis and study of the molecular mechanisms involved in CMS could improve breeding programs. In the present study, anthers of the pepper CMS line HZ1A and its maintainer line HZ1B were collected from stages S1, S2, and S3 for transcriptome sequencing. A total of 47.95 million clean reads were obtained, and the reads were assembled into 31,603 unigenes. We obtained 42 (27 up-regulated and 15 down-regulated), 691 (346 up-regulated and 345 down-regulated), and 709 (281 up-regulated and 428 down-regulated) differentially expressed genes (DEGs) in stages S1, S2, and S3, respectively. Through Gene Ontology (GO) analysis, the DEGs were found to be composed of 46 functional groups. Two GO terms involved in photosynthesis, photosynthesis (GO:0015986) and photosystem I (GO:0009522), may be related to CMS. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, oxidative phosphorylation (ko00190) and phenylpropanoid biosynthesis (ko00940) were significantly enriched in the S1 and S2 stages, respectively. Through the analysis of 104 lipid metabolism-related DEGs, four significantly enriched KEGG pathways may help to regulate male sterility during anther development. The mitochondrial genes orf470 and atp6 were identified as candidate genes of male sterility for the CMS line HZ1A. Overall, the results will provide insights into the molecular mechanisms of pepper CMS.

Transcriptome Analysis of Seed Weight Plasticity in Brassica napus

A critical barrier to improving crop yield is the trade-off between seed weight (SW) and seed number (SN), which has been commonly reported in several crops, including Brassica napus. Despite the agronomic relevance of this issue, the molecular factors involved in the interaction between SW and SN are largely unknown in crops. In this work, we performed a detailed transcriptomic analysis of 48 seed samples obtained from two rapeseed spring genotypes subjected to different source-sink (S-S) ratios in order to examine the relationship between SW and SN under different field conditions. A multifactorial analysis of the RNA-seq data was used to identify a group of 1014 genes exclusively regulated by the S-S ratio. We found that a reduction in the S-S ratio during seed filling induces the expression of genes involved in sucrose transport, seed weight, and stress responses. Moreover, we identified five co-expression modules that are positively correlated with SW and negatively correlated with SN. Interestingly, one of these modules was significantly enriched in transcription factors (TFs). Furthermore, our network analysis predicted several NAC TFs as major hubs underlying SW and SN compensation. Taken together, our study provides novel insights into the molecular factors associated with the SW-SN relationship in rapeseed and identifies TFs as potential targets when improving crop yield.

Salt tolerance involved candidate genes in rice: an integrative meta-analysis approach

BACKGROUND

Salinity, as one of the main abiotic stresses, critically threatens growth and fertility of main food crops including rice in the world. To get insight into the molecular mechanisms by which tolerant genotypes responds to the salinity stress, we propose an integrative meta-analysis approach to find the key genes involved in salinity tolerance. Herein, a genome-wide meta-analysis, using microarray and RNA-seq data was conducted which resulted in the identification of differentially expressed genes (DEGs) under salinity stress at tolerant rice genotypes. DEGs were then confirmed by meta-QTL analysis and literature review.

RESULTS

A total of 3449 DEGs were detected in 46 meta-QTL positions, among which 1286, 86, 1729 and 348 DEGs were observed in root, shoot, seedling, and leaves tissues, respectively. Moreover, functional annotation of DEGs located in the meta-QTLs suggested some involved biological processes (e.g., ion transport, regulation of transcription, cell wall organization and modification as well as response to stress) and molecular function terms (e.g., transporter activity, transcription factor activity and oxidoreductase activity). Remarkably, 23 potential candidate genes were detected in Saltol and hotspot-regions overlying original QTLs for both yield components and ion homeostasis traits; among which, there were many unreported salinity-responsive genes. Some promising candidate genes were detected such as pectinesterase, peroxidase, transcription regulator, high-affinity potassium transporter, cell wall organization, protein serine/threonine phosphatase, and CBS domain cotaining protein.

CONCLUSIONS

The obtained results indicated that, the salt tolerant genotypes use qualified mechanisms particularly in sensing and signalling of the salt stress, regulation of transcription, ionic homeostasis, and Reactive Oxygen Species (ROS) scavenging in response to the salt stress.

Comparative analysis of sugarcane root transcriptome in response to the plant growth-promoting Burkholderia anthina MYSP113

The diazotrophic Burkholderia anthina MYSP113 is a vital plant growth-promoting bacteria and sugarcane root association. The present study based on a detailed analysis of sugarcane root transcriptome by using the HiSeq-Illumina platform in response to the strain MYSP113. The bacterium was initially isolated from the rhizosphere of sugarcane. To better understand biological, cellular, and molecular mechanisms, a de novo transcriptomic assembly of sugarcane root was performed. HiSeq-Illumina platformwas employed for the sequencing of an overall of 16 libraries at a 2×100 bp configuration. Differentially expressed genes (DEGs) analysis identified altered gene expression in 370 genes (total of 199 up-regulated genes and 171 down-regulated genes). Deciphering the huge datasets, concerning the functioning and production of biological systems, a high throughput genome sequencing analysis was attempted with Gene ontology functional analyses and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The report revealed a total of 148930 unigenes. 70414 (47.28%) of them were annotated successfully to Gene Ontology (GO) terms. 774 at 45 days, 4985 of 30 days and 15 days of 6846 terms were significantly regulated. GO analysis revealed that many genes involved in the metabolic, oxidation-reduction process and biological regulatory processes in response to strain MYSP113 and significantly enriched as compare to the control. Moreover, KEGG enriched results show that differentially expressed genes were classified into different pathway categories involved in various processes, such as nitrogen metabolism, plant hormone signal transduction, etc. The sample correlation analyses could help examine the similarity at the gene expression level. The reliability of the observed differential gene expression patterns was validated with quantitative real-time PCR (qRT-PCR). Additionally, plant enzymes activities such as peroxidase and superoxide dismutase were significantly increased in plant roots after the inoculation of strain MYSP113. The results of the research may help in understanding the plant growth-promoting rhizobacteria and plant interaction.

Protective Role of Leaf Variegation in Pittosporum tobira under Low Temperature: Insights into the Physio-Biochemical and Molecular Mechanisms

Leaf variegation has been demonstrated to have adaptive functions such as cold tolerance. Pittosporum tobira is an ornamental plant with natural leaf variegated cultivars grown in temperate regions. Herein, we investigated the role of leaf variegation in low temperature responses by comparing variegated “Variegatum” and non-variegated “Green Pittosporum” cultivars. We found that leaf variegation is associated with impaired chloroplast development in the yellow sector, reduced chlorophyll content, strong accumulation of carotenoids and high levels of ROS. However, the photosynthetic efficiency was not obviously impaired in the variegated leaves. Also, leaf variegation plays low temperature protective function since “Variegatum” displayed strong and efficient ROS-scavenging enzymatic systems to buffer cold (10 °C)-induced damages. Transcriptome analysis under cold conditions revealed 309 differentially expressed genes between both cultivars. Distinctly, the strong cold response observed in “Variegatum” was essentially attributed to the up-regulation of HSP70/90 genes involved in cellular homeostasis; up-regulation of POD genes responsible for cell detoxification and up-regulation of FAD2 genes and subsequent down-regulation of GDSL genes leading to high accumulation of polyunsaturated fatty acids for cell membrane fluidity. Overall, our results indicated that leaf variegation is associated with changes in physiological, biochemical and molecular components playing low temperature protective function in P. tobira.

Unravelling the MicroRNA-Mediated Gene Regulation in Developing Pongamia Seeds by High-Throughput Small RNA Profiling

Pongamia (Millettia pinnata syn. Pongamia pinnata) is a multipurpose biofuel tree which can withstand a variety of abiotic stresses. Commercial applications of Pongamia trees may substantially benefit from improvements in their oil-seed productivity, which is governed by complex regulatory mechanisms underlying seed development. MicroRNAs (miRNAs) are important molecular regulators of plant development, while relatively little is known about their roles in seed development, especially for woody plants. In this study, we identified 236 conserved miRNAs within 49 families and 143 novel miRNAs via deep sequencing of Pongamia seeds sampled at three developmental phases. For these miRNAs, 1327 target genes were computationally predicted. Furthermore, 115 differentially expressed miRNAs (DEmiRs) between successive developmental phases were sorted out. The DEmiR-targeted genes were preferentially enriched in the functional categories associated with DNA damage repair and photosynthesis. The combined analyses of expression profiles for DEmiRs and functional annotations for their target genes revealed the involvements of both conserved and novel miRNA-target modules in Pongamia seed development. Quantitative Real-Time PCR validated the expression changes of 15 DEmiRs as well as the opposite expression changes of six targets. These results provide valuable miRNA candidates for further functional characterization and breeding practice in Pongamia and other oilseed plants.

De novo sequencing of tree peony (Paeonia suffruticosa) transcriptome to identify critical genes involved in flowering and floral organ development

BACKGROUND

Tree peony (Paeonia suffruticosa Andrews) is a globally famous ornamental flower, with large and colorful flowers and abundant flower types. However, a relatively short and uniform flowering period hinders the applications and production of ornamental tree peony. Unfortunately, the molecular mechanism of regulating flowering time and floral organ development in tree peony has yet to be elucidated. Because of the absence of genomic information, 454-based transcriptome sequence technology for de novo transcriptomics was used to identify the critical flowering genes using re-blooming, non-re-blooming, and wild species of tree peonies.

RESULTS

A total of 29,275 unigenes were obtained from the bud transcriptome, with an N50 of 776 bp. The average length of unigenes was 677.18 bp, and the longest sequence was 5815 bp. Functional annotation showed that 22,823, 17,321, 13,312, 20,041, and 9940 unigenes were annotated by NCBI-NR, Swiss-Prot, COG, GO, and KEGG, respectively. Within the differentially expressed genes (DEGs) 64 flowering-related genes were identified and some important flowering genes were also characterized by bioinformatics methods, reverse transcript polymerase chain reaction (RT-PCR), and rapid-amplification of cDNA ends (RACE). Then, the putative genetic network of flowering induction pathways and a floral organ development model were put forward, according to the comparisons of DEGs in any two samples and expression levels of the important flowering genes in differentiated buds, buds from different developmental stages, and with GA or vernalization treated. In tree peony, five pathways (long day, vernalization, autonomous, age, and gibberellin) regulated flowering, and the floral organ development followed an ABCE model. Moreover, it was also found that the genes PsAP1, PsCOL1, PsCRY1, PsCRY2, PsFT, PsLFY, PsLHY, PsGI, PsSOC1, and PsVIN3 probably regulated re-blooming of tree peony.

CONCLUSION

This study provides a comprehensive report on the flowering-related genes in tree peony for the first time and investigated the expression levels of the critical flowering related genes in buds of different cultivars, developmental stages, differentiated primordium, and flower parts. These results could provide valuable insights into the molecular mechanisms of flowering time regulation and floral organ development.

Overexpression of Acyl-ACP Thioesterases, CpFatB4 and CpFatB5, Induce Distinct Gene Expression Reprogramming in Developing Seeds of Brassica napus

We examined the substrate preference of Cuphea paucipetala acyl-ACP thioesterases, CpFatB4 and CpFatB5, and gene expression changes associated with the modification of lipid composition in the seed, using Brassica napus transgenic plants overexpressing CpFatB4 or CpFatB5 under the control of a seed-specific promoter. CpFatB4 seeds contained a higher level of total saturated fatty acid (FA) content, with 4.3 times increase in 16:0 palmitic acid, whereas CpFatB5 seeds showed approximately 3% accumulation of 10:0 and 12:0 medium-chain FAs, and a small increase in other saturated FAs, resulting in higher levels of total saturated FAs. RNA-Seq analysis using entire developing pods at 8, 25, and 45 days after flowering (DAF) showed up-regulation of genes for β-ketoacyl-acyl carrier protein synthase I/II, stearoyl-ACP desaturase, oleate desaturase, and linoleate desaturase, which could increase unsaturated FAs and possibly compensate for the increase in 16:0 palmitic acid at 45 DAF in CpFatB4 transgenic plants. In CpFatB5 transgenic plants, many putative chloroplast- or mitochondria-encoded genes were identified as differentially expressed. Our results report comprehensive gene expression changes induced by alterations of seed FA composition and reveal potential targets for further genetic modifications.

Comparative transcriptomic analysis of high- and low-oil Camellia oleifera reveals a coordinated mechanism for the regulation of upstream and downstream multigenes for high oleic acid accumulation

Tea oil camellia (Camellia oleifera) is an important woody oil tree in southern China. However, little is known regarding the molecular mechanisms that contribute to high oleic acid accumulation in tea oil camellia. Herein, we measured the oil content and fatty acid compositions of high- and low-oil tea oil camellia seeds and investigated the global gene expression profiles by RNA-seq. The results showed that at the early, second and third seed developmental stages, a total of 64, 253, and 124 genes, respectively, were significantly differentially expressed between the high- and low-oil cultivars. Gene ontology (GO) enrichment analysis of the identified differentially expressed transcription factors (TFs; ABI3, FUS3, LEC1, WRI1, TTG2 and DOF4.6) revealed some critical GO terms associated with oil biosynthesis and fatty acid accumulation, including glycolysis, zinc ion binding, positive regulation of fatty acid biosynthetic process, triglyceride biosynthetic process, seed coat development, abscisic acid-mediated signaling pathway and embryo development. Comprehensive comparisons of transcriptomic profiles and expression analysis of multigenes based on qRT-PCR showed that coordinated high expression of the upstream genes HAD, EAR and KASI directly increased the relative levels of C16:0-ACP, which provided enough precursor resources for oleic acid biosynthesis. Continuous high expression of the SAD gene accelerated oleic acid synthesis and accumulation, and coordinated low expression of the downstream genes FAD2, FAD3, FAD7, FAD8 and FAE1 decreased the consumption of oleic acid for conversion. The coordinated regulation of these multigenes ensures the high accumulation of oleic acid in the seeds of tea oil camellia. Our data represent a comprehensive transcriptomic study of high- and low-oil tea oil camellia, not only increasing the number of sequences associated with lipid biosynthesis and fatty acid accumulation in public resource databases but also providing a scientific basis for genetic improvement of the oleic acid content in woody oil trees.

Dynamic Transcriptome Changes Related to Oil Accumulation in Developing Soybean Seeds

Soybean is one of the most important oil crops in the world. Revealing the molecular basis and exploring key candidate genes for seed oil synthesis has great significance for soybean improvement. In this study, we found that oil accumulation rates and gene expression levels changed dynamically during soybean seed development. The expression levels of genes in metabolic pathways such as carbon fixation, photosynthesis, glycolysis, and fatty acid biosynthesis were significantly up-regulated during the rapid accumulation of oil in developing soybean seeds. Through weighted correlation network analysis, we identified six co-expression modules associated with soybean seed oil content and the pink module was the most positively correlated (r = 0.83, p = 7 × 10-4) network. Through the integration of differential expression and co-expression analysis, we predicted 124 candidate genes potentially affecting soybean seed oil content, including seven genes in lipid metabolism pathway, two genes involved in glycolysis, one gene in sucrose metabolism, and 12 genes belonged to transcription factors as well as other categories. Among these, three genes (GmABI3b, GmNFYA and GmFAD2-1B) have been shown to control oil and fatty acid content in soybean seeds, and other newly identified candidate genes would broaden our knowledge to understand the molecular basis for oil accumulation in soybean seeds.

Comparative Transcriptome Analysis between Fertile and CMS Flower Buds in Wucai (Brassica campestris L.)

Background

Wucai (Brassica campestris L. ssp. chinensis var. rosularis Tsen) is a variant of nonheading Chinese cabbage (Brassica campestris L.), which is one of the major vegetables in China. Cytoplasmic male sterility (CMS) has been used for Wucai breeding in recent years. However, the underlying molecular mechanism of Wucai CMS remains unclear. In this study, the phenotypic and cytological features of Wucai CMS were observed by anatomical analysis, and a comparative transcriptome analysis was carried out to identify genes related to male sterility using Illumina RNA sequencing technology (RNA-Seq).

Results

Microscopic observation demonstrated that tapetum development was abnormal in the CMS line, which failed to produce fertile pollen. Bioinformatics analysis detected 4430 differentially expressed genes (DEGs) between the fertile and sterile flower buds. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to better understand the functions of these DEGs. Among the DEGs, 35 genes (53 DEGS) were implicated in anther and pollen development, and 11 genes were involved in pollen cell wall formation and modification; most of these showed downregulated expression in sterile buds. In addition, several genes related to tapetum development (A6, AMS, MS1, MYB39, and TSM1) and a few genes annotated to flowering (CO, AP3, VIN3, FLC, FT, and AGL) were detected and confirmed by qRT-PCR as being expressed at the meiosis, tetrad, and uninucleate microspore stages, thus implying possible roles in specifying or determining the fate and development of the tapetum, male gametophyte and stamen. Moreover, the top four largest transcription factor families (MYB, bHLH, NAC and WRKY) were analyzed, and most showed reduced expression in sterile buds. These differentially expressed transcription factors might result in abortion of pollen development in Wucai.

Conclusion

The present comparative transcriptome analysis suggested that many key genes and transcription factors involved in anther development show reduced gene expression patterns in the CMS line, which might contribute to male sterility in Wucai. This study provides valuable information for a better understanding of CMS molecular mechanisms and functional genome studies in Wucai.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5331-4) contains supplementary material, which is available to authorized users.

An integrated omics analysis reveals molecular mechanisms that are associated with differences in seed oil content between Glycine max and Brassica napus

BACKGROUND

Rapeseed (Brassica napus L.) and soybean (Glycine max L.) seeds are rich in both protein and oil, which are major sources of biofuels and nutrition. Although the difference in seed oil content between soybean (~ 20%) and rapeseed (~ 40%) exists, little is known about its underlying molecular mechanism.

RESULTS

An integrated omics analysis was performed in soybean, rapeseed, Arabidopsis (Arabidopsis thaliana L. Heynh), and sesame (Sesamum indicum L.), based on Arabidopsis acyl-lipid metabolism- and carbon metabolism-related genes. As a result, candidate genes and their transcription factors and microRNAs, along with phylogenetic analysis and co-expression network analysis of the PEPC gene family, were found to be largely associated with the difference between the two species. First, three soybean genes (Glyma.13G148600, Glyma.13G207900 and Glyma.12G122900) co-expressed with GmPEPC1 are specifically enriched during seed storage protein accumulation stages, while the expression of BnPEPC1 is putatively inhibited by bna-miR169, and two genes BnSTKA and BnCKII are co-expressed with BnPEPC1 and are specifically associated with plant circadian rhythm, which are related to seed oil biosynthesis. Then, in de novo fatty acid synthesis there are rapeseed-specific genes encoding subunits β-CT (BnaC05g37990D) and BCCP1 (BnaA03g06000D) of heterogeneous ACCase, which could interfere with synthesis rate, and β-CT is positively regulated by four transcription factors (BnaA01g37250D, BnaA02g26190D, BnaC01g01040D and BnaC07g21470D). In triglyceride synthesis, GmLPAAT2 is putatively inhibited by three miRNAs (gma-miR171, gma-miR1516 and gma-miR5775). Finally, in rapeseed there was evidence for the expansion of gene families, CALO, OBO and STERO, related to lipid storage, and the contraction of gene families, LOX, LAH and HSI2, related to oil degradation.

CONCLUSIONS

The molecular mechanisms associated with differences in seed oil content provide the basis for future breeding efforts to improve seed oil content.

Cloning and characterization of EgGDSL, a gene associated with oil content in oil palm

Oil palm (Elaeis guineensis, Jacq.) is a key tropical oil crop, which provides over one third of the global vegetable oil production, but few genes related to oil yield have been characterized. In this study, a GDSL esterase/lipase gene, which was significantly associated with oil content, was isolated from oil palm and designated as EgGDSL. Its functional characterization was carried out through ectopic expression in Arabidopsis ecotype Col-0. It was shown that expression of EgGDSL in Arabidopsis led to the increased total fatty acid content by 9.5% compared with the wild type. Further analysis of the fatty acid composition revealed that stearic acid (18:0) increased in the seeds of the transgenic lines, but the levels of linoleic acid (18:2) plus 11-eicosenoic acid drastically declined. Quantitative real-time PCR (qPCR) revealed that in oil palm, EgGDSL was highly expressed in mesocarp followed by leaf, and the expression level was very low in the root. The expression level of EgGDSL gene began to increase at two months after flowering (MAF) and reached its peak by four MAF, then declined rapidly, and reached its lowest level during the mature period (6 MAF). The EgGDSL gene was more highly expressed in oil palm trees with high oil content than that with low oil content, demonstrating that the transcription level of EgGDSL correlated with the amount of oil accumulation. The gene may be valuable for engineering fatty acid metabolism in crop improvement programmes and for marker-assisted breeding.

Temporal transcriptome profiling of developing seeds reveals a concerted gene regulation in relation to oil accumulation in Pongamia (Millettia pinnata)

BACKGROUND

Pongamia (Millettia pinnata syn. Pongamia pinnata), an oilseed legume species, is emerging as potential feedstock for sustainable biodiesel production. Breeding Pongamia for favorable traits in commercial application will rely on a comprehensive understanding of molecular mechanism regulating oil accumulation during its seed development. To date, only limited genomic or transcript sequences are available for Pongamia, while a temporal transcriptome profiling of developing seeds is still lacking in this species.

RESULTS

In this work, we conducted a time-series analysis of morphological and physiological characters, oil contents and compositions, as well as global gene expression profiles in developing Pongamia seeds. Firstly, three major developmental phases were characterized based on the combined evidences from embryonic shape, seed weight, seed moisture content, and seed color. Then, the gene expression levels at these three phases were quantified by RNA-Seq analyses with three biological replicates from each phase. Nearly 94% of unigenes were expressed at all three phases, whereas only less than 2% of unigenes were exclusively expressed at one of these phases. A total of 8881 differentially expressed genes (DEGs) were identified between phases. Furthermore, the qRT-PCR analyses for 10 DEGs involved in lipid metabolism demonstrated a good reliability of our RNA-Seq data in temporal gene expression profiling. We observed a dramatic increase in seed oil content from the embryogenesis phase to the early seed-filling phase, followed by a steady and moderate increase towards the maximum at the desiccation phase. We proposed that a highly active expression of most genes related to fatty acid (FA) and triacylglycerol (TAG) biosynthesis at the embryogenesis phase might trigger both the substantial oil accumulation and the membrane lipid synthesis for rapid cell proliferation at this phase, while a concerted reactivation of TAG synthesis-related genes at the desiccation phase might further promote storage lipid synthesis to achieve the maximum content of seed oils.

CONCLUSIONS

This study not only built a bridge between gene expression profiles and oil accumulation in developing seeds, but also laid a foundation for future attempts on genetic engineering of Pongamia varieties to acquire higher oil yield or improved oil properties for biofuel applications.

RNA-seq transcriptome analysis of the immature seeds of two Brassica napus lines with extremely different thousand-seed weight to identify the candidate genes related to seed weight

Brassica napus is an important oilseed crop worldwide. Although seed weight is the main determinant of seed yield, few studies have focused on the molecular mechanisms that regulate seed weight in B. napus. In this study, the immature seeds of G-42 and 7–9, two B. napus doubled haploid (DH) lines with extremely different thousand-seed weight (TSW), were selected for a transcriptome analysis to determine the regulatory mechanisms underlying seed weight at the whole gene expression level and to identify candidate genes related to seed weight. A total of 2,251 new genes and 2,205 differentially expressed genes (DEGs) were obtained via RNA-seq (RNA sequencing). Among these genes, 1,747 (77.61%) new genes and 2020 (91.61%) DEGs were successfully annotated. Of these DEGs, 1,118 were up-regulated and 1,087 were down-regulated in the large-seed line. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database analysis indicated that 15 DEGs were involved in ubiquitin-mediated proteolysis and proteasome pathways, which might participate in regulating seed weight. The Gene Ontology (GO) database indicated that 222 DEGs were associated with the biological process or molecular function categories related to seed weight, such as cell division, cell size and cell cycle regulation, seed development, nutrient reservoir activity, and proteasome-mediated ubiquitin-dependent protein catabolic processes. Moreover, 50 DEGs encoding key enzymes or proteins were identified that likely participate in regulating seed weight. A DEG (GSBRNA2T00037121001) identified by the transcriptome analysis was also previously identified in a quantitative trait locus (QTL) region for seed weight via SLAF-seq (Specific Locus Amplified Fragment sequencing). Finally, the expression of 10 DEGs with putative roles in seed weight and the expression of the DEG GSBRNA2T00037121001 were confirmed by a quantitative real-time reverse transcription PCR (qRT-PCR) analysis, and the results were consistent with the RNA sequencing data. This work has provided new insights on the molecular mechanisms underlying seed weight-related biosynthesis and has laid a solid foundation for further improvements to the seed yield of oil crops.

De Novo Transcriptome Sequencing and the Hypothetical Cold Response Mode of Saussurea involucrata in Extreme Cold Environments

Saussurea involucrata grows in high mountain areas covered by snow throughout the year. The temperature of this habitat can change drastically in one day. To gain a better understanding of the cold response signaling pathways and molecular metabolic reactions involved in cold stress tolerance, genome-wide transcriptional analyses were performed using RNA-Seq technologies. A total of 199,758 transcripts were assembled, producing 138,540 unigenes with 46.8 Gb clean data. Overall, 184,416 (92.32%) transcripts were successfully annotated. The 365 transcription factors identified (292 unigenes) belonged to 49 transcription factor families associated with cold stress responses. A total of 343 transcripts on the signal transduction (132 upregulated and 212 downregulated in at least any one of the conditions) were strongly affected by cold temperature, such as the CBL-interacting serine/threonine-protein kinase (CIPKs), receptor-like protein kinases, and protein kinases. The circadian rhythm pathway was activated by cold adaptation, which was necessary to endure the severe temperature changes within a day. There were 346 differentially expressed genes (DEGs) related to transport, of which 138 were upregulated and 22 were downregulated in at least any one of the conditions. Under cold stress conditions, transcriptional regulation, molecular transport, and signal transduction were involved in the adaptation to low temperature in S. involucrata. These findings contribute to our understanding of the adaptation of plants to harsh environments and the survival traits of S. involucrata. In addition, the present study provides insight into the molecular mechanisms of chilling and freezing tolerance.

Integration of omics approaches to understand oil/protein content during seed development in oilseed crops

Oilseed crops, especially soybean (Glycine max) and canola/rapeseed (Brassica napus), produce seeds that are rich in both proteins and oils and that are major sources of energy and nutrition worldwide. Most of the nutritional content in the seed is accumulated in the embryo during the seed filling stages of seed development. Understanding the metabolic pathways that are active during seed filling and how they are regulated are essential prerequisites to crop improvement. In this review, we summarize various omics studies of soybean and canola/rapeseed during seed filling, with emphasis on oil and protein traits, to gain a systems-level understanding of seed development. Currently, most (80-85%) of the soybean and rapeseed reference genomes have been sequenced (950 and 850 megabases, respectively). Parallel to these efforts, extensive omics datasets from different seed filling stages have become available. Transcriptome and proteome studies have detected preponderance of starch metabolism and glycolysis enzymes to be the possible cause of higher oil in B. napus compared to other crops. Small RNAome studies performed during the seed filling stages have revealed miRNA-mediated regulation of transcription factors, with the suggestion that this interaction could be responsible for transitioning the seeds from embryogenesis to maturation. In addition, progress made in dissecting the regulation of de novo fatty acid synthesis and protein storage pathways is described. Advances in high-throughput omics and comprehensive tissue-specific analyses make this an exciting time to attempt knowledge-driven investigation of complex regulatory pathways.

An integrated analysis of QTL mapping and RNA sequencing provides further insights and promising candidates for pod number variation in rapeseed (Brassica napus L.)

BACKGROUND

As the most important yield component in rapeseed (Brassica napus L.), pod number is determined by a series of successive growth and development processes. Pod number shows extensive variation in rapeseed natural germplasm, which is valuable for genetic improvement. However, the genetic and especially the molecular mechanism for this kind of variation are poorly understood. In this study, we conducted QTL mapping and RNA sequencing, respectively, using the BnaZNRIL population and its two parental cultivars Zhongshuang11 and No.73290 which showed significant difference in pod number, primarily due to the difference in floral organ number.

RESULT

A total of eight QTLs for pod number were identified using BnaZNRIL population with a high-density SNP linkage map, each was distributed on seven linkage groups and explained 5.8-11.9% of phenotypic variance. Then, they were integrated with those previously detected in BnaZNF₂ population (deriving from same parents) and resulted in 15 consensus-QTLs. Of which, seven QTLs were identical to other studies, whereas the other eight should be novel. RNA sequencing of the shoot apical meristem (SAM) at the formation stage of floral bud primordia identified 9135 genes that were differentially expressed between the two parents. Gene ontology (GO) analysis showed that the top two enriched groups were S-assimilation, providing an essential nutrient for the synthesis of diverse metabolites, and polyamine metabolism, serving as second messengers that play an essential role in flowering genes initiation. KEGG analysis showed that the top three overrepresented pathways were carbohydrate (707 genes), amino acid (390 genes) and lipid metabolisms (322 genes). In silico mapping showed that 647 DEGs were located within the confidence intervals of 15 consensus QTLs. Based on annotations of Arabidopsis homologs corresponding to DEGs, nine genes related to meristem growth and development were considered as promising candidates for six QTLs.

CONCLUSION

In this study, we discovered the first repeatable major QTL for pod number in rapeseed. In addition, RNA sequencing was performed for SAM in rapeseed, which provides new insights into the determination of floral organ number. Furthermore, the integration of DEGs and QTLs identified promising candidates for further gene cloning and mechanism study.

Genome-wide transcriptomic analysis uncovers the molecular basis underlying early flowering and apetalous characteristic in Brassica napus L

Floral transition and petal onset, as two main aspects of flower development, are crucial to rapeseed evolutionary success and yield formation. Currently, very little is known regarding the genetic architecture that regulates flowering time and petal morphogenesis in Brassica napus. In the present study, a genome-wide transcriptomic analysis was performed with an absolutely apetalous and early flowering line, APL01, and a normally petalled line, PL01, using high-throughput RNA sequencing. In total, 13,205 differential expressed genes were detected, of which 6111 genes were significantly down-regulated, while 7094 genes were significantly up-regulated in the young inflorescences of APL01 compared with PL01. The expression levels of a vast number of genes involved in protein biosynthesis were altered in response to the early flowering and apetalous character. Based on the putative rapeseed flowering genes, an early flowering network, mainly comprised of vernalization and photoperiod pathways, was built. Additionally, 36 putative upstream genes possibly governing the apetalous character of line APL01 were identified, and six genes potentially regulating petal origination were obtained by combining with three petal-related quantitative trait loci. These findings will facilitate understanding of the molecular mechanisms underlying floral transition and petal initiation in B. napus.

Time-Series Analyses of Transcriptomes and Proteomes Reveal Molecular Networks Underlying Oil Accumulation in Canola

Understanding the regulation of lipid metabolism is vital for genetic engineering of canola (Brassica napus L.) to increase oil yield or modify oil composition. We conducted time-series analyses of transcriptomes and proteomes to uncover the molecular networks associated with oil accumulation and dynamic changes in these networks in canola. The expression levels of genes and proteins were measured at 2, 4, 6, and 8 weeks after pollination (WAP). Our results show that the biosynthesis of fatty acids is a dominant cellular process from 2 to 6 WAP, while the degradation mainly happens after 6 WAP. We found that genes in almost every node of fatty acid synthesis pathway were significantly up-regulated during oil accumulation. Moreover, significant expression changes of two genes, acetyl-CoA carboxylase and acyl-ACP desaturase, were detected on both transcriptomic and proteomic levels. We confirmed the temporal expression patterns revealed by the transcriptomic analyses using quantitative real-time PCR experiments. The gene set association analysis show that the biosynthesis of fatty acids and unsaturated fatty acids are the most significant biological processes from 2-4 WAP and 4-6 WAP, respectively, which is consistent with the results of time-series analyses. These results not only provide insight into the mechanisms underlying lipid metabolism, but also reveal novel candidate genes that are worth further investigation for their values in the genetic engineering of canola.