Comparative Transcriptomics Analysis of Brassica napus L. during Seed Maturation Reveals Dynamic Changes in Gene Expression between Embryos and Seed Coats and Distinct Expression Profiles of Acyl-CoA-Binding Proteins for Lipid Accumulation

Interaction of Soybean (Glycine max (L.) Merr.) Class II ACBPs with MPK2 and SAPK2 Kinases: New Insights into the Regulatory Mechanisms of Plant ACBPs

Plant acyl-CoA-binding proteins (ACBPs) function in plant development and stress responses, with some ACBPs interacting with protein partners. This study tested the interaction between two Class II GmACBPs (Glycine max ACBPs) and seven kinases, using yeast two-hybrid (Y2H) assays and bimolecular fluorescence complementation (BiFC). The results revealed that both GmACBP3.1 and GmACBP4.1 interact with two soybean kinases, a mitogen-activated protein kinase MPK2, and a serine/threonine-protein kinase SAPK2, highlighting the significance of the ankyrin-repeat (ANK) domain in facilitating protein-protein interactions. Moreover, an in vitro kinase assay and subsequent Phos-tag SDS-PAGE determined that GmMPK2 and GmSAPK2 possess the ability to phosphorylate Class II GmACBPs. Additionally, the kinase-specific phosphosites for Class II GmACBPs were predicted using databases. The HDOCK server was also utilized to predict the binding models of Class II GmACBPs with these two kinases, and the results indicated that the affected residues were located in the ANK region of Class II GmACBPs in both docking models, aligning with the findings of the Y2H and BiFC experiments. This is the first report describing the interaction between Class II GmACBPs and kinases, suggesting that Class II GmACBPs have potential as phospho-proteins that impact signaling pathways.

Ribosome footprint profiling enables elucidating the systemic regulation of fatty acid accumulation in Acer truncatum

BACKGROUND

The accumulation of fatty acids in plants covers a wide range of functions in plant physiology and thereby affects adaptations and characteristics of species. As the famous woody oilseed crop, Acer truncatum accumulates unsaturated fatty acids and could serve as the model to understand the regulation and trait formation in oil-accumulation crops. Here, we performed Ribosome footprint profiling combing with a multi-omics strategy towards vital time points during seed development, and finally constructed systematic profiling from transcription to proteomes. Additionally, we characterized the small open reading frames (ORFs) and revealed that the translational efficiencies of focused genes were highly influenced by their sequence features.

RESULTS

The comprehensive multi-omics analysis of lipid metabolism was conducted in A. truncatum. We applied the Ribo-seq and RNA-seq techniques, and the analyses of transcriptional and translational profiles of seeds collected at 85 and 115 DAF were compared. Key members of biosynthesis-related structural genes (LACS, FAD2, FAD3, and KCS) were characterized fully. More meaningfully, the regulators (MYB, ABI, bZIP, and Dof) were identified and revealed to affect lipid biosynthesis via post-translational regulations. The translational features results showed that translation efficiency tended to be lower for the genes with a translated uORF than for the genes with a non-translated uORF. They provide new insights into the global mechanisms underlying the developmental regulation of lipid metabolism.

CONCLUSIONS

We performed Ribosome footprint profiling combing with a multi-omics strategy in A. truncatum seed development, which provides an example of the use of Ribosome footprint profiling in deciphering the complex regulation network and will be useful for elucidating the metabolism of A. truncatum seed oil and the regulatory mechanisms.

Organo-mineral complexes alter bacterial composition and induce carbon and nitrogen cycling in the rhizosphere

It is widely thought that organo-mineral complexes (OMCs) stabilize organic matter via mineral adsorption. Recent studies have demonstrated that root exudates can activate OMCs, but the influence of OMCs on plant rhizosphere, which is among the most active areas for microbes, has not been thoroughly researched. In this study, a pot experiment using Brassica napus was conducted to investigate the effects of OMCs on plant rhizosphere. The result showed that OMC addition significantly promoted the growth of B. napus compared to the prevalent fertilization (PF, chemical fertilizer + chicken compost) treatment. Specifically, OMC addition increased the relative abundance (RA) of nitrogen-fixing bacteria and the bacterial α-diversity, and the operational taxonomic unit (OTU) group with RA > 0.5% in the OMC-treated rhizosphere was the result of a deterministic assembly process with homogeneous selection. Gene abundance related to nitrogen cycling and the soil chemical analysis demonstrated that the OMC-altered bacterial community induced nitrogen fixation and converted nitrate to ammonium. The upregulated carbon sequestration pathway genes and the increased soil microbial biomass carbon (23.68%) demonstrated that the bacterial-induced carbon storage in the rhizosphere was activated. This study shows that the addition of OMCs can influence the biogeochemical carbon and nitrogen cycling via regulating microorganisms in the rhizosphere. The findings provide fresh insights into the effects of OMCs on the biogeochemical cycling of important elements and suggest a promising strategy for improving soil productivity.

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.

Roles of acyl-CoA-binding proteins in plant reproduction

Acyl-CoA-binding proteins (ACBPs) constitute a well-conserved family of proteins in eukaryotes that are important in stress responses and development. Past studies have shown that ACBPs are involved in maintaining, transporting and protecting acyl-CoA esters during lipid biosynthesis in plants, mammals, and yeast. ACBPs show differential expression and various binding affinities for acyl-CoA esters. Hence, ACBPs can play a crucial part in maintaining lipid homeostasis. This review summarizes the functions of ACBPs during the stages of reproduction in plants and other organisms. A comprehensive understanding on the roles of ACBPs during plant reproduction may lead to opportunities in crop improvement in agriculture.

Transgenic manipulation of triacylglycerol biosynthetic enzymes in B. napus alters lipid-associated gene expression and lipid metabolism

Oilseed rape (Brassica napus) is an important crop that is cultivated for the oil (mainly triacylglycerol; TAG) it produces in its seeds. TAG synthesis is controlled mainly by key enzymes in the Kennedy pathway, such as glycerol 3-phosphate acyltransferase (GPAT), lysophosphatidate acyltransferase (LPAT) and diacylglycerol acyltransferase (DGAT) but can also be produced from phosphoglycerides such as phosphatidylcholine (PC) by the activity of the enzyme phospholipid: diacylglycerol acyltransferase (PDAT). To evaluate the potential for these enzymes to alter oil yields or composition, we analysed transgenic B. napus lines which overexpressed GPAT, LPAT or PDAT using heterologous transgenes from Arabidopsis and Nasturtium and examined lipid profiles and changes in gene expression in these lines compared to WT. Distinct changes in PC and TAG abundance and spatial distribution in embryonic tissues were observed in some of the transgenic lines, together with altered expression of genes involved generally in acyl-lipid metabolism. Overall our results show that up-regulation of these key enzymes differentially affects lipid composition and distribution as well as lipid-associated gene expression, providing important information which could be used to improve crop properties by metabolic engineering.

Genome-Wide Analysis of Myeloblastosis-Related Genes in Brassica napus L. and Positive Modulation of Osmotic Tolerance by BnMRD107

Myeloblastosis (MYB)-related transcription factors comprise a large subfamily of the MYB family. They play significant roles in plant development and in stress responses. However, MYB-related proteins have not been comprehensively investigated in rapeseed (Brassica napus L.). In the present study, a genome-wide analysis of MYB-related transcription factors was performed in rapeseed. We identified 251 Brassica napus MYB (BnMYB)-related members, which were divided phylogenetically into five clades. Evolutionary analysis suggested that whole genome duplication and segmental duplication events have played a significant role in the expansion of BnMYB-related gene family. Selective pressure of BnMYB-related genes was estimated using the Ka/Ks ratio, which indicated that BnMYB-related genes underwent strong purifying selection during evolution. In silico analysis showed that various development-associated, phytohormone-responsive, and stress-related cis-acting regulatory elements were enriched in the promoter regions of BnMYB-related genes. Furthermore, MYB-related genes with tissue or organ-specific, stress-responsive expression patterns were identified in B. napus based on temporospatial and abiotic stress expression profiles. Among the stress-responsive MYB-related genes, BnMRD107 was strongly induced by drought stress, and was therefore selected for functional study. Rapeseed seedlings overexpressing BnMRD107 showed improved resistance to osmotic stress. Our findings not only lay a foundation for further functional characterization of BnMYB-related genes, but also provide valuable clues to determine candidate genes for future genetic improvement of B. napus.

In silico Analysis of Acyl-CoA-Binding Protein Expression in Soybean

Plant acyl-CoA-binding proteins (ACBPs) form a highly conserved protein family that binds to acyl-CoA esters as well as other lipid and protein interactors to function in developmental and stress responses. This protein family had been extensively studied in non-leguminous species such as Arabidopsis thaliana (thale cress), Oryza sativa (rice), and Brassica napus (oilseed rape). However, the characterization of soybean (Glycine max) ACBPs, designated GmACBPs, has remained unreported although this legume is a globally important crop cultivated for its high oil and protein content, and plays a significant role in the food and chemical industries. In this study, 11 members of the GmACBP family from four classes, comprising Class I (small), Class II (ankyrin repeats), Class III (large), and Class IV (kelch motif), were identified. For each class, more than one copy occurred and their domain architecture including the acyl-CoA-binding domain was compared with Arabidopsis and rice. The expression profile, tertiary structure and subcellular localization of each GmACBP were predicted, and the similarities and differences between GmACBPs and other plant ACBPs were deduced. A potential role for some Class III GmACBPs in nodulation, not previously encountered in non-leguminous ACBPs, has emerged. Interestingly, the sole member of Class III ACBP in each of non-leguminous Arabidopsis and rice had been previously identified in plant-pathogen interactions. As plant ACBPs are known to play important roles in development and responses to abiotic and biotic stresses, the in silico expression profiles on GmACBPs, gathered from data mining of RNA-sequencing and microarray analyses, will lay the foundation for future studies in their applications in biotechnology.

RICE ACYL-COA-BINDING PROTEIN6 Affects Acyl-CoA Homeostasis and Growth in Rice

BACKGROUNDS

Acyl-coenzyme A (CoA) esters are important intermediates in lipid metabolism with regulatory properties. Acyl-CoA-binding proteins bind and transport acyl-CoAs to fulfill these functions. RICE ACYL-COA-BINDING PROTEIN6 (OsACBP6) is currently the only one peroxisome-localized plant ACBP that has been proposed to be involved in β-oxidation in transgenic Arabidopsis. The role of the peroxisomal ACBP (OsACBP6) in rice (Oryza sativa) was investigated.

RESULTS

Here, we report on the function of OsACBP6 in rice. The osacbp6 mutant showed diminished growth with reduction in root meristem activity and leaf growth. Acyl-CoA profiling and lipidomic analysis revealed an increase in acyl-CoA content and a slight triacylglycerol accumulation caused by the loss of OsACBP6. Comparative transcriptomic analysis discerned the biological processes arising from the loss of OsACBP6. Reduced response to oxidative stress was represented by a decline in gene expression of a group of peroxidases and peroxidase activities. An elevation in hydrogen peroxide was observed in both roots and shoots/leaves of osacbp6. Taken together, loss of OsACBP6 not only resulted in a disruption of the acyl-CoA homeostasis but also peroxidase-dependent reactive oxygen species (ROS) homeostasis. In contrast, osacbp6-complemented transgenic rice displayed similar phenotype to the wild type rice, supporting a role for OsACBP6 in the maintenance of the acyl-CoA pool and ROS homeostasis. Furthermore, quantification of plant hormones supported the findings observed in the transcriptome and an increase in jasmonic acid level occurred in osacbp6.

CONCLUSIONS

In summary, OsACBP6 appears to be required for the efficient utilization of acyl-CoAs. Disruption of OsACBP6 compromises growth and led to provoked defense response, suggesting a correlation of enhanced acyl-CoAs content with defense responses.

Characterization and function of a sunflower (Helianthus annuus L.) Class II acyl-CoA-binding protein

Acyl-CoA-binding proteins (ACBP) bind to long-chain acyl-CoA esters and phospholipids, enhancing the activity of different acyltransferases in animals and plants. Nevertheless, the role of these proteins in the synthesis of triacylglycerols (TAGs) remains unclear. Here, we cloned a cDNA encoding HaACBP1, a Class II ACBP from sunflower (Helianthus annuus), one of the world's most important oilseed crop plants. Transcriptome analysis of this gene revealed strong expression in developing seeds from 16 to 30 days after flowering. The recombinant protein (rHaACBP1) was expressed in Escherichia coli and purified to be studied by in vitro isothermal titration calorimetry and for phospholipid binding. Its high affinity for saturated palmitoyl-CoA (16:0-CoA; K_D 0.11 μM) and stearoyl-CoA (18:0-CoA; K_D 0.13 μM) esters suggests that rHaACBP1 could act in acyl-CoA transfer pathways that involve saturated acyl derivatives. Furthermore, rHaACBP1 also binds to both oleoyl-CoA (18:1-CoA; K_D 6.4 μM) and linoleoyl-CoA (18:2-CoA; K_D 21.4 μM) esters, the main acyl-CoA substrates used to synthesise the TAGs that accumulate in sunflower seeds. Interestingly, rHaACBP1 also appears to bind to different species of phosphatidylcholines (dioleoyl-PC and dilinoleoyl-PC), glycerolipids that are also involved in TAG synthesis, and while it interacts with dioleoyl-PA, this is less prominent than its binding to the PC derivative. Expression of rHaACBP in yeast alters its fatty acid composition, as well as the composition and size of the host acyl-CoA pool. These results suggest that HaACBP1 may potentially fulfil a role in the transport and trafficking of acyl-CoAs during sunflower seed development.

Global Transcriptome and Correlation Analysis Reveal Cultivar-Specific Molecular Signatures Associated with Fruit Development and Fatty Acid Determination in Camellia oleifera Abel

Background

Oil-tea Camellia is a very important edible oil plant widely distributed in southern China. Tea oil extracted from the oil-tea Camellia seeds is beneficial to health and is considered as a health edible oil. We attempt to identify genes related to fatty acid biosynthesis in an oil-tea Camellia seed kernel, generated a comprehensive transcriptome analysis of the seed kernel at different developmental stages, and explore optimal picking time of fruit. Material and Methods. A gas chromatography-mass spectrometer was used to detect the content of various fatty acids in samples. Transcriptome analysis was performed to detect gene dynamics and corresponding functions.

Results

Multiple phenotypic data were counted in detail, including the oil content, oleic acid content, linoleic acid content, linolenic acid content, fruit weight, fruit height, fruit diameter, single seed weight, seed length, and seed width in different developmental stages, which indicate that a majority of indicators increased with the development of oil-tea Camellia. The transcriptomics was conducted to perform a comprehensive and system-level view on dynamic gene expression networks for different developmental stages. Short Time-series Expression Miner (STEM) analysis of XL106 (the 6 time points) and XL210 (8 time points) was performed to screen related fatty acid (FA) gene set, from which 1041 candidate genes related to FA were selected in XL106 and 202 related genes were screened in XL210 based on GO and KEGG enrichment. Then, candidate genes and trait dataset were combined to conduct correlation analysis, and 10 genes were found to be strongly connected with several key traits.

Conclusions

The multiple phenotypic data revealed the dynamic law of changes during the picking stage. Transcriptomic analysis identified a large number of potential key regulatory factors that can control the oil content of dried kernels, oleic acid, linoleic acid, linolenic acid, fresh seed rate, and kernel-to-seed ratio, thereby providing a new insight into the molecular networks underlying the picking stage of oil-tea Camellia, which provides a theoretical basis for the optimal fruit picking point.

De novo transcriptome assembly for the five major organs of Zanthoxylum armatum and the identification of genes involved in terpenoid compound and fatty acid metabolism

Background

Zanthoxylum armatum (Z. armatum) is a highly economically important tree that presents a special numbing taste. However, the underlying regulatory mechanism of the numbing taste remains poorly understood. Thus, the elucidation of the key genes associated with numbing taste biosynthesis pathways is critical for providing genetic information on Z. armatumand the breeding of high-quality germplasms of this species.

Results

Here, de novo transcriptome assembly was performed for the five major organs of Z. armatum, including the roots, stems, leaf buds, mature leaves and fruits. A total of 111,318 unigenes were generated with an average length of 1014 bp. Additionally, a large number of SSRs were obtained to improve our understanding of the phylogeny and genetics of Z. armatum. The organ-specific unigenes of the five major samples were screened and annotated via GO and KEGG enrichment analysis. A total of 53 and 34 unigenes that were exclusively upregulated in fruit samples were identified as candidate unigenes for terpenoid biosynthesis or fatty acid biosynthesis, elongation and degradation pathways, respectively. Moreover, 40 days after fertilization (Fr4 stage) could be an important period for the accumulation of terpenoid compounds during the fruit development and maturation of Z. armatum. The Fr4 stage could be a key point at which the first few steps of the fatty acid biosynthesis process are promoted, and the catalysis of subsequent reactions could be significantly induced at 62 days after fertilization (Fr6 stage).

Conclusions

The present study realized de novo transcriptome assembly for the five major organs of Z. armatum. To the best of our knowledge, this study provides the first comprehensive analysis revealing the genes underlying the special numbing taste of Z. armatum. The assembled transcriptome profiles expand the available genetic information on this species and will contribute to gene functional studies, which will aid in the engineering of high-quality cultivars of Z. armatum.

Big data from small tissues: extraction of high-quality RNA for RNA-sequencing from different oilseed Brassica seed tissues during seed development

BACKGROUND

Obtaining high-quality RNA for gene expression analyses from different seed tissues is challenging due to the presence of various contaminants, such as polyphenols, polysaccharides and lipids which interfere with RNA extraction methods. At present, the available protocols for extracting RNA from seeds require high amounts of tissue and are mainly focused on extracting RNA from whole seeds. However, extracting RNA at the tissue level enables more detailed studies regarding tissue specific transcriptomes during seed development.

RESULTS

Seeds from heart stage embryo to mature developmental stages of Brassica napus and B. oleracea were sampled for isolation of the embryo, endosperm and seed coat tissues. Ovules and ovary wall tissue were also collected from pre-fertilized buds. Subsequent to testing several RNA extraction methods, modifications applied to E.Z.N.A. Plant RNA and Picopure RNA Isolation kit extraction methods resulted in RNA with high yield and quality. Furthermore, the use of polyvinylpolypyrrolidone for seed coats and endosperm at green stages resulted in high-quality RNA. As a result of the introduced modifications to established RNA extraction methods, the RNA from all the above-mentioned tissues presented clear 28S and 18S bands and high RIN values, ranging from 7.0 to 10.0. The protocols reported in this study are not only suitable for different and challenging seed tissue types, but also enable the extraction of high-quality RNA using only 2 to 3 mg of starting tissue.

CONCLUSIONS

Here, we present efficient, reproducible and reliable high-quality RNA extraction methods for diverse oilseed Brassica spp reproductive tissue types including pre-fertilization and developing seed tissues for diploid and polyploid species. The high-quality RNA obtained is suitable for RNA-Sequencing and subsequent gene expression analysis.