Microarray analysis of gene expression in seeds of Brassica napus planted in Nanjing (altitude: 8.9 m), Xining (altitude: 2261.2 m) and Lhasa (altitude: 3658 m) with different oil content

Increasing oil content in Brassica oilseed species

Brassica oilseed species are the third most important in the world, providing approximately 15 % of the total vegetable oils. Three species (Brassica rapa, B. juncea, B. napus) dominate with B. napus being the most common in Canada, China and Europe. Originally, B. napus was a crop producing seed with high erucic acid content, which still persists today, to some extent, and is used for industrial purposes. In contrast, cultivars which produce seed used for food and feed are low erucic acid cultivars which also have reduced glucosinolate content. Because of the limit to agricultural land, recent efforts have been made to increase productivity of oil crops, including Brassica oilseed species. In this article, we have detailed research in this regard. We have covered modern genetic, genomic and metabolic control analysis approaches to identifying potential targets for the manipulation of seed oil content. Details of work on the use of quantitative trait loci, genome-wide association and comparative functional genomics to highlight factors influencing seed oil accumulation are given and functional proteins which can affect this process are discussed. In summary, a wide variety of inputs are proving useful for the improvement of Brassica oilseed species, as major sources of global vegetable oil.

Effects of Low Nighttime Temperature on Fatty Acid Content in Developing Seeds from Brassica napus L. Based on RNA-Seq and Metabolome

Brassica napus L. is a vital plant oil resource worldwide. The fatty acid biosynthesis and oil accumulation in its seeds are controlled by several genetic and environmental factors, including daytime and nighttime temperatures. We analyzed changes in oleic and erucic acid content in two double haploid (DH) lines, DH0729, a weakly temperature-sensitive line, and DH0815, a strongly temperature-sensitive line, derived from B. napus plants grown at different altitudes (1600, 1800, 2000, 2200, and 2400 m a.s.l., 28.85° N, 112.35° E) and nighttime temperatures (20/18, 20/16, 20/13 and 20/10 °C, daytime/nighttime temperature). Based on medium- and long-chain fatty acid metabolites, the total oleic acid content 35 and 43 days after flowering was significantly lower in low nighttime temperature (LNT, 20/13 °C) plants than in high nighttime temperature (HNT, 20/18 °C) plants (HNT: 58-62%; LNT: 49-54%; an average decrease of 9%), and the total erucic acid content was significantly lower in HNT than in LNT plants (HNT: 1-2%; LNT: 8-13%; an average increase of 10%). An RNA-seq analysis showed that the expression levels of SAD (LOC106366808), ECR (LOC106396280), KCS (LOC106419344), KAR (LOC106367337), HB1(LOC106430193), and DOF5 (LOC111211868) in STSL seeds increased under LNT conditions. In STSL seeds, a base mutation in the cis-acting element involved in low-temperature responsiveness (LTR), the HB1 and KCS promoter caused loss of sensitivity to low temperatures, whereas that of the KCS promoter caused increased sensitivity to low temperatures.

Temporal and Spatial Variation of PM2.5 in Xining, Northeast of the Qinghai–Xizang (Tibet) Plateau

PM2.5 was sampled from January 2017 to May 2018 at an urban, suburban, industrial, and rural sites in Xining. The annual mean of PM2.5 was highest at the urban site and lowest at the rural site, with an average of 51.5 ± 48.9 and 26.4 ± 17.8 μg·m−3, respectively. The average PM2.5 concentration of the industrial and suburban sites was 42.8 ± 27.4 and 37.2 ± 23.7 μg·m−3, respectively. All sites except for the rural had concentrations above the ambient air quality standards of China (GB3095-2012). The highest concentration of PM2.5 at all sites was observed in winter, followed by spring, autumn, and summer. The concentration of major constituents showed statistically significant seasonal and spatial variation. The highest concentrations of organic carbon (OC), elemental carbon (EC), water-soluble organic carbon (WSOC), and water-soluble inorganic ions (WSIIs) were found at the urban site in winter. The average concentration of F− was higher than that in many studies, especially at the industrial site where the annual average concentration of F− was 1.5 ± 1.7 μg·m−3. The range of sulfur oxidation ratio (SOR) was 0.1–0.18 and nitrogen oxidation ratio (NOR) was 0.02–0.1 in Xining. The higher SO42−/NO3− indicates that coal combustion has greater impact than vehicle emissions. The results of the potential source contribution function (PSCF) suggest that air mass from middle- and large-scale transport from the western areas of Xining have contributed to the higher level of PM2.5. On the basis of the positive matrix factorization (PMF) model, it was found that aerosols from salt lakes and dust were the main sources of PM2.5 in Xining, accounting for 26.3% of aerosol total mass. During the sandstorms, the concentration of PM2.5 increased sharply, and the concentrations of Na+, Ca2+ and Mg2+ were 1.13–2.70, 1.68–4.41, and 1.15–5.12 times higher, respectively, than annual average concentration, implying that aerosols were mainly from dust and the largest saltwater lake, Qinghai Lake, and many other salt lakes in the province of Qinghai. Time-of-flight secondary ion mass spectrometry (ToF-SIMS) was utilized to study the surface components of PM2.5 and F− was found to be increasingly distributed from the surface to inside the particles. We determined that the extremely high PM2.5 concentration appears to be due to an episode of heavy pollution resulting from the combination of sandstorms and the burning of fireworks.

Diurnal biomarkers reveal key photosynthetic genes associated with increased oil palm yield

To investigate limiters of photosynthate assimilation in the carbon-source limited crop, oil palm (Elaeis guineensis Jacq.), we measured differential metabolite, gene expression and the gas exchange in leaves in an open field for palms with distinct mesocarp oil content. We observed higher concentrations of glucose 1-phosphate, glucose 6-phosphate, sucrose 6-phosphate, and sucrose in high-oil content palms with the greatest difference being at 11:00 (p-value ≤0.05) immediately after the period of low morning light intensity. Three important photosynthetic genes were identified using differentially expressed gene analysis (DEGs) and were found to be significantly enriched through Gene Ontology (GO) and pathway enrichment: chlorophyll a-b binding protein (CAB-13), photosystem I (PSI), and Ferredoxin-NADP reductase (FNR), particularly for sampling points at non-peak light (11:00 and 19:00), ranging from 3.3-fold (PSI) and 5.6-fold (FNR) to 10.3-fold (CAB-13). Subsequent gas exchange measurements further supported increased carbon assimilation through higher level of internal CO₂ concentration (Ci), stomatal conductance (g_s) and transpiration rate (E) in high-oil content palms. The selection for higher expression of key photosynthesis genes together with CO₂ assimilation under low light is likely to be important for crop improvement, in particular at full maturity and under high density planting regimes where light competition exists between palms.

Comparative RNA-Seq Analysis of High- and Low-Oil Yellow Horn During Embryonic Development

Yellow horn (Xanthoceras sorbifolium Bunge) is an endemic oil-rich shrub that has been widely cultivated in northern China for bioactive oil production. However, little is known regarding the molecular mechanisms that contribute to oil content in yellow horn. Herein, we measured the oil contents of high- and low-oil yellow horn embryo tissues at four developmental stages and investigated the global gene expression profiles through RNA-seq. The results found that at 40, 54, 68, and 81 days after anthesis, a total of 762, 664, 599, and 124 genes, respectively, were significantly differentially expressed between the high- and low-oil lines. Gene ontology (GO) enrichment analysis revealed some critical GO terms related to oil accumulation, including acyl-[acyl-carrier-protein] desaturase activity, pyruvate kinase activity, acetyl-CoA carboxylase activity, and seed oil body biogenesis. The identified differentially expressed genes also included several transcription factors, such as, AP2-EREBP family members, B3 domain proteins and C2C2-Dof proteins. Several genes involved in fatty acid (FA) biosynthesis, glycolysis/gluconeogenesis, and pyruvate metabolism were also up-regulated in the high-oil line at different developmental stages. Our findings indicate that the higher oil accumulation in high-oil yellow horn could be mostly driven by increased FA biosynthesis and carbon supply, i.e. a source effect.

Differential gene expression at different stages of mesocarp development in high- and low-yielding oil palm

BACKGROUND

The oil yield trait of oil palm is expected to involve multiple genes, environmental influences and interactions. Many of the underlying mechanisms that contribute to oil yield are still poorly understood. In this study, we used a microarray approach to study the gene expression profiles of mesocarp tissue at different developmental stages, comparing genetically related high- and low- oil yielding palms to identify genes that contributed to the higher oil-yielding palm and might contribute to the wider genetic improvement of oil palm breeding populations.

RESULTS

A total of 3412 (2001 annotated) gene candidates were found to be significantly differentially expressed between high- and low-yielding palms at at least one of the different stages of mesocarp development evaluated. Gene Ontologies (GO) enrichment analysis identified 28 significantly enriched GO terms, including regulation of transcription, fatty acid biosynthesis and metabolic processes. These differentially expressed genes comprise several transcription factors, such as, bHLH, Dof zinc finger proteins and MADS box proteins. Several genes involved in glycolysis, TCA, and fatty acid biosynthesis pathways were also found up-regulated in high-yielding oil palm, among them; pyruvate dehydrogenase E1 component Subunit Beta (PDH), ATP-citrate lyase, β- ketoacyl-ACP synthases I (KAS I), β- ketoacyl-ACP synthases III (KAS III) and ketoacyl-ACP reductase (KAR). Sucrose metabolism-related genes such as Invertase, Sucrose Synthase 2 and Sucrose Phosphatase 2 were found to be down-regulated in high-yielding oil palms, compared to the lower yield palms.

CONCLUSIONS

Our findings indicate that a higher carbon flux (channeled through down-regulation of the Sucrose Synthase 2 pathway) was being utilized by up-regulated genes involved in glycolysis, TCA and fatty acid biosynthesis leading to enhanced oil production in the high-yielding oil palm. These findings are an important stepping stone to understand the processes that lead to production of high-yielding oil palms and have implications for breeding to maximize oil production.

Differential gene expression at different stages of mesocarp development in high- and low-yielding oil palm

Background

The oil yield trait of oil palm is expected to involve multiple genes, environmental influences and interactions. Many of the underlying mechanisms that contribute to oil yield are still poorly understood. In this study, we used a microarray approach to study the gene expression profiles of mesocarp tissue at different developmental stages, comparing genetically related high- and low- oil yielding palms to identify genes that contributed to the higher oil-yielding palm and might contribute to the wider genetic improvement of oil palm breeding populations.

Results

A total of 3412 (2001 annotated) gene candidates were found to be significantly differentially expressed between high- and low-yielding palms at at least one of the different stages of mesocarp development evaluated. Gene Ontologies (GO) enrichment analysis identified 28 significantly enriched GO terms, including regulation of transcription, fatty acid biosynthesis and metabolic processes. These differentially expressed genes comprise several transcription factors, such as, bHLH, Dof zinc finger proteins and MADS box proteins. Several genes involved in glycolysis, TCA, and fatty acid biosynthesis pathways were also found up-regulated in high-yielding oil palm, among them; pyruvate dehydrogenase E1 component Subunit Beta (PDH), ATP-citrate lyase, β- ketoacyl-ACP synthases I (KAS I), β- ketoacyl-ACP synthases III (KAS III) and ketoacyl-ACP reductase (KAR). Sucrose metabolism-related genes such as Invertase, Sucrose Synthase 2 and Sucrose Phosphatase 2 were found to be down-regulated in high-yielding oil palms, compared to the lower yield palms.

Conclusions

Our findings indicate that a higher carbon flux (channeled through down-regulation of the Sucrose Synthase 2 pathway) was being utilized by up-regulated genes involved in glycolysis, TCA and fatty acid biosynthesis leading to enhanced oil production in the high-yielding oil palm. These findings are an important stepping stone to understand the processes that lead to production of high-yielding oil palms and have implications for breeding to maximize oil production.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3855-7) contains supplementary material, which is available to authorized users.

Microarray analysis of differentially expressed genes between Brassica napus strains with high- and low-oleic acid contents

An increase in oleic acid (C18:1) content is a desirable trait. Despite the critical roles of the two desaturases, FAD2 and FAD3, in the control of fatty acid desaturation, a dispute remains over whether inactivation of their genes alone is sufficient enough to generate the high-oleic trait. To address this question, we employed microarray technology to investigate the difference in gene expression profile between two different Brassica napus strains with high-C18:1 (71.71%) and low-C18:1 (55.6%) contents, respectively. Our study revealed 562 differentially expressed genes, of which 194 genes were up-regulated and 368 down-regulated. Based on the Gene Ontology classification, these genes were classified into 23 functional categories. Three of the up-regulated genes represent B. napus homologs of Arabidopsis genes encoding a cytosolic isoform of pyruvate kinase (AT3G55810), Δ9 acyl-lipid desaturase (AT1G06080, ADS1) and fatty acyl-ACP thioesterase B (AT1G08510), respectively. Conversely, the homologs of two Arabidopsis sequences encoding Δ9 acyl-lipid desaturase (AT2G31360, ADS2) and FAD3 desaturase (AT2G29980) were down-regulated in the high-oleic acid strain. Furthermore, 60 differentially expressed genes were classified as associated with relevant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Collectively, our results suggest that expressing the high-oleic acid trait may require a coordinated regulation of diverse regulatory and metabolic gene networks in addition to inactivation of the FAD2 and FAD3 genes in the oilseed. A set of the differentially expressed genes identified in this study will facilitate our efforts to tap the germplasms with the potential to express the high-oleic acid trait.

Small RNA profiling in two Brassica napus cultivars identifies microRNAs with oil production- and development-correlated expression and new small RNA classes

MicroRNAs (miRNAs) and small interfering RNAs are important regulators of plant development and seed formation, yet their population and abundance in the oil crop Brassica napus are still not well understood, especially at different developmental stages and among cultivars with varied seed oil contents. Here, we systematically analyzed the small RNA expression profiles of Brassica napus seeds at early embryonic developmental stages in high-oil-content and low-oil-content B. napus cultivars, both cultured in two environments. A total of 50 conserved miRNAs and 9 new miRNAs were identified, together with some new miRNA targets. Expression analysis revealed some miRNAs with varied expression levels in different seed oil content cultivars or at different embryonic developmental stages. A large number of 23-nucleotide small RNAs with specific nucleotide composition preferences were also identified, which may present new classes of functional small RNAs.

Enhanced seed oil production in canola by conditional expression of Brassica napus LEAFY COTYLEDON1 and LEC1-LIKE in developing seeds

The seed oil content in oilseed crops is a major selection trait to breeders. In Arabidopsis (Arabidopsis thaliana), LEAFY COTYLEDON1 (LEC1) and LEC1-LIKE (L1L) are key regulators of fatty acid biosynthesis. Overexpression of AtLEC1 and its orthologs in canola (Brassica napus), BnLEC1 and BnL1L, causes an increased fatty acid level in transgenic Arabidopsis plants, which, however, also show severe developmental abnormalities. Here, we use truncated napin A promoters, which retain the seed-specific expression pattern but with a reduced expression level, to drive the expression of BnLEC1 and BnL1L in transgenic canola. Conditional expression of BnLEC1 and BnL1L increases the seed oil content by 2% to 20% and has no detrimental effects on major agronomic traits. In the transgenic canola, expression of a subset of genes involved in fatty acid biosynthesis and glycolysis is up-regulated in developing seeds. Moreover, the BnLEC1 transgene enhances the expression of several genes involved in Suc synthesis and transport in developing seeds and the silique wall. Consistently, the accumulation of Suc and Fru is increased in developing seeds of the transgenic rapeseed, suggesting the increased carbon flux to fatty acid biosynthesis. These results demonstrate that BnLEC1 and BnL1L are reliable targets for genetic improvement of rapeseed in seed oil production.

Characterization and expression profile analysis of a sucrose synthase gene from common bean (Phaseolus vulgaris L.) during seed development

A full-length cDNA encoding common bean (Phaseolus vulgaris L.) sucrose synthase (designated as Pv_BAT93 Sus), which catalyses the synthesis and cleavage of sucrose, was isolated from seeds at 15 days after pollination (DAP) by rapid amplification of cDNA ends (RACE). The full-length cDNA of Pv_BAT93 Sus had a 2,418 bp open reading frame (ORF) encoding a protein of 806 amino acid residues. Sequence comparison analysis showed that Pv_BAT93 Sus was very similar to several members of the sucrose synthase family of other plant species. Tissue expression pattern analysis showed that Pv_BAT93 Sus was expressed in leaves, flowers, stems, roots, cotyledons, and particularly during seed development. Expression studies using in situ hybridization revealed altered spatial and temporal patterns of Sus expression in the EMS mutant relative to wild-type and confirmed Sus expression in common bean developing seeds. The expression and accumulation of Sus mRNA was clearly shown in several tissues, such as the suspensor and embryo, but also in the transfer cells and endothelium. The results highlight the diverse roles that Sus might play during seed development in common bean.

Isolation and characterization of a novel BcMF14 gene from Brassica campestris ssp. chinensis

A putative RALF (rapid alkalinization factor)-like gene (GenBank accession number EF523517), named BcMF14, was isolated from Chinese cabbage (Brassica campestris L. ssp. chinensis Makino, syn. B. rapa ssp. chinensis) by rapid amplification of cDNA ends (RACE) based on a cDNA-AFLP differential fragment exclusively expressed in fertile line. Semi-quantitative reverse transcription-polymerase chain reaction (RT-PCR) discovered that BcMF14 was prominently expressed in stage four and five flower buds of fertile line, no expression in vegetative structures or in sterility line. Detailed RT-PCR illuminated its strong expression in stamens. Successful suppression of BcMF14 gene expression greatly reduced the normal pollen grains. The frequency of abnormal pollen grains was 48.95% in the mutant with many shriveled pollen grains with irregular shape and some larger ones with deep hollows along the germination ditch. Pollen germination was stopped because of the severely twisted pollen tubes. These results demonstrate a potential role of the BcMF14 gene in the development of male gametogenesis in Chinese cabbage.

Cloning, biochemical characterisation, tissue localisation and possible post-translational regulatory mechanism of the cytosolic phosphoglucose isomerase from developing sunflower seeds

Lipid biosynthesis in developing sunflower (Helianthus annuus L.) seeds requires reducing power. One of the main sources of cellular NADPH is the oxidative pentose phosphate pathway (OPPP), generated from the oxidation of glucose-6-phosphate. This glycolytic intermediate, which can be imported to the plastid and enter in the OPPP, is the substrate and product of cytosolic phosphoglucose isomerase (cPGI, EC 5.3.1.9). In this report, we describe the cloning of a full-length cDNA encoding cPGI from developing sunflower seeds. The sequence was predicted to code for a protein of 566 residues characterised by the presence of two sugar isomerase domains. This cDNA was heterologously expressed in Escherichia coli as a His-tagged protein. The recombinant protein was purified using immobilised metal ion affinity chromatography and biochemically characterised. The enzyme had a specific activity of 1,436 micromol min(-1) mg(-1) and 1,011 micromol min(-1) mg(-1) protein when the reaction was initiated with glucose-6-phosphate and fructose-6-phosphate, respectively. Activity was not affected by erythrose-4-phosphate, but was inhibited by 6-P gluconate and glyceraldehyde-3-phosphate. A polyclonal immune serum was raised against the purified enzyme, allowing the study of protein levels during the period of active lipid synthesis in seeds. These results were compared with PGI activity profiles and mRNA expression levels obtained from Q-PCR studies. Our results point to the existence of a possible post-translational regulatory mechanism during seed development. Immunolocalisation of the protein in seed tissues further indicated that cPGI is highly expressed in the procambial ring.

Analysis of the differential expression of the genes related to Brassica napus seed development

To screen the genes related to Brassica napus seed development at pattern formation and maturation stages, the suppression subtractive cDNA libraries of B. napus cultivar Zhongshuang 6 were constructed with its embryos at 10 days after flowering (10 DAF) and 30 days after flowering (30 DAF) through suppression subtractive hybridization (SSH) technology. The positive clones were screened by PCR and dot blot hybridization, and then sequenced. High quality expressed sequence tags (ESTs) were used for COG functional classification with COGNITOR software, as well as analysis and annotation with BLAST software. Tissue-specific detection of five genes screened was performed by RT-PCR in root, stem, leaf, flower, bud, pod, and embryo tissues. The insert size ranged from 100 to 900 bp, with an average size of about 500 bp. According to COG functional classification database, the differentially expressed genes mainly involved in lipid and amino acid metabolism, signal transduction, post-transcriptional modification, and etc. The results from RT-PCR detection of the five differentially expressed genes indicated that genes 2-96 and 2-352 presented embryo-specific expression, gene 1-385 expressed in parts of tissues, and genes 1-71 and 1-682 expressed in all tissues. Two genes were found to be involved in seed development, lipid and protein metabolisms, two genes may be involved in signal transduction, one gene could not match with the homologous sequences known to date and was likely a new gene. These results are helpful for future gene cloning and their functional analysis.

A GA-insensitive dwarf mutant of Brassica napus L. correlated with mutation in pyrimidine box in the promoter of GID1

A dwarf mutant from Brassica napus, namely NDF-1, which was derived from a high doubled haploid (DH) line '3529'(Brassica napus L.) of which seeds were jointly treated with chemical inducers and fast neutron bombardment, was revealed that dwarfism is under the control of a major gene(designated as ndf1) with a mainly additive effect and non-significant dominance effect. The germination and hypocotyls elongation response of dwarf mutants after exogenous GA and uniconazol application showed NDF-1 was a gibberellin insensitive dwarf. We cloned the Brassica napus GID1 gene, named BnGID1, and found it was the ortholog of AtGID1a. The sequence blasting of the BnGID1 genes from NDF-1 and wild type showed there was no mutant in the gene. But the quantitative RT-PCR analysis of GID1 EST pointed out the mutation was caused by the low-level expression of BnGID1 gene. After sequenced the BnGID1 gene's upstream, we found three bases mutated in the pyrimidine box (P-box) of the BnGID1 promoter, which is linkage with the dwarf mutant.