High throughput sequencing of small RNAs reveals dynamic microRNAs expression of lipid metabolism during Camellia oleifera and C. meiocarpa seed natural drying

The mechanism of Andrena camellia in digesting toxic sugars

Camellia oleifera is an economically and medicinally valuable oilseed crop. Honeybee, the most abundant pollinator, rarely visits C. oleifera because of the toxic sugars in the nectar and pollen. These toxic sugars cannot be fully digested by honeybees and inhibit the process of synthesizing trehalose in honeybees. C. oleifera exhibits self-incompatibility, and its pollination heavily depends on Andrena camellia. However, the mechanism by which A. camellia digests toxic sugars in C. oleifera nectar and pollen remains unknown. Consequently, we identified and validated four single-copy genes (α-N-acetyl galactosamine-like, galactokinase, galactose-1-phosphate uridyltransferase, and UDP-galactose-4'-epimerase, abbreviated as NAGA-like, GALK, GALT, and GALE) essential for detoxifying toxic sugars in vitro. Then, we cloned the four genes into Escherichia coli, and expressed enzyme successfully degraded the toxic sugars. The phylogeny suggests that the genes were conserved and functionally diverged among the evolution. These results provide novel insights into pollinator detoxification during co-evolution.

Chromosome-scale genome assembly of oil-tea tree Camellia crapnelliana

Camellia crapnelliana Tutch., belonging to the Theaceae family, is an excellent landscape tree species with high ornamental values. It is particularly an important woody oil-bearing plant species with high ecological, economic, and medicinal values. Here, we first report the chromosome-scale reference genome of C. crapnelliana with integrated technologies of SMRT, Hi-C and Illumina sequencing platforms. The genome assembly had a total length of ~2.94 Gb with contig N50 of ~67.5 Mb, and ~96.34% of contigs were assigned to 15 chromosomes. In total, we predicted 37,390 protein-coding genes, ~99.00% of which could be functionally annotated. The chromosome-scale genome of C. crapnelliana will become valuable resources for understanding the genetic basis of the fatty acid biosynthesis, and greatly facilitate the exploration and conservation of C. crapnelliana.

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

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

Genomic and genetic advances of oiltea-camellia (Camellia oleifera)

Oiltea-camellia (C. oleifera) is a widely cultivated woody oil crop in Southern China and Southeast Asia. The genome of oiltea-camellia was very complex and not well explored. Recently, genomes of three oiltea-camellia species were sequenced and assembled, multi-omic studies of oiltea-camellia were carried out and provided a better understanding of this important woody oil crop. In this review, we summarized the recent assembly of the reference genomes of oiltea-camellia, genes related to economic traits (flowering, photosynthesis, yield and oil component), disease resistance (anthracnose) and environmental stress tolerances (drought, cold, heat and nutrient deficiency). We also discussed future directions of integrating multiple omics for evaluating genetic resources and mining key genes of important traits, and the application of new molecular breeding and gene editing technologies to accelerate the breeding process of oiltea-camellia.

Untargeted metabolism approach reveals difference of varieties of bud and relation among characteristics of grafting seedlings in Camellia oleifera

Camellia oleifera is one of the essential wood oil trees in the world. C.oleifera was propagated by nurse seedling grafting. Since the scion of C.oleifera had a significant regulated effect on the properties of rootstock after grafting and impacted on the growth of the grafted seedlings, it was necessary to understand the characteristics of buds among varieties to cultivate high-quality grafted seedlings. The metabolome was thought to be a powerful tool for understanding connecting phenotype-genotype interactions, which has an important impact on plant growth and development. In this study, UPLC-MS was used to determine the metabolites of the apical buds of CL3, CL4, CL40, and CL53 spring shoots after 30 days of sprout and to measure the growth characteristics of roots and stems after grafting. Metabolomics analysis revealed 554 kinds of metabolites were significant differences among four varieties, and 29 metabolic pathways were identified to have significant changes (p< 0.05), including carboxylic acids and derivatives, fatty Acyls, organooxygen compounds, and prenol lipids metabolites. The metabolites appeared in all varieties, including phenethyl rutinoside in glycosyl compounds and hovenidulcioside A1 in terpene glycosides. Metabolite-metabolite correlations in varieties revealed more complex patterns in relation to bud and enabled the recognition of key metabolites (e.g., Glutamate, (±)Catechin, GA₅₂, ABA, and cs-Zeatin) affecting grafting and growth ability. Each variety has a unique metabolite type and correlation network relationship. Differentiated metabolites showed different growth trends for development after grafting. Many metabolites regulate the growth of scions in buds before grafting, which plays a crucial role in the growth of seedlings after grafting. It not only regulates the growth of roots but also affects the development of this stem. Finally, those results were associated with the genetic background of each cultivar, showing that metabolites could be potentially used as indicators for the genetic background, indicating that metabolites could potentially be used as indicators for seedling growth characteristics. Together, this study will enrich the theoretical basis of seedling growth and lay a foundation for further research on the molecular regulation mechanism interaction between rootstock and scion, rootstock growth, and the development of grafted seedlings after grafting.

Assessment of the Genetic Relationship and Population Structure in Oil-Tea Camellia Species Using Simple Sequence Repeat (SSR) Markers

Oil-tea camellia trees, the collective term for a class of economically valuable woody oil crops in China, have attracted extensive attention because of their rich nutritional and pharmaceutical value. This study aimed to analyze the genetic relationship and genetic diversity of oil-tea camellia species using polymorphic SSR markers. One-hundred and forty samples of five species were tested for genetic diversity using twenty-four SSR markers. In this study, a total of 385 alleles were identified using 24 SSR markers, and the average number of alleles per locus was 16.0417. The average Shannon's information index (I) was 0.1890, and the percentages of polymorphic loci (P) of oil-tea camellia trees were 7.79-79.48%, indicating that oil-tea camellia trees have low diversity. Analysis of molecular variance (AMOVA) showed that the majority of genetic variation (77%) was within populations, and a small fraction (23%) occurred among populations. Principal coordinate analysis (PCoA) results indicated that the first two principal axes explained 7.30% (PC1) and 6.68% (PC2) of the total variance, respectively. Both UPGMA and PCoA divided the 140 accessions into three groups. Camellia oleifera clustered into one class, Camellia vietnamensis and Camellia gauchowensis clustered into one class, and Camellia crapnelliana and Camellia chekiangoleosa clustered into another class. It could be speculated that the genetic relationship of C. vietnamensis and C. gauchowensis is quite close. SSR markers could reflect the genetic relationship among oil-tea camellia germplasm resources, and the results of this study could provide comprehensive information on the conservation, collection, and breeding of oil-tea camellia germplasms.

MicroRNA-mediated regulation of lipid metabolism in virus-infected Emiliania huxleyi

The interactions between Emiliania huxleyi and E. huxleyi virus (EhV) regulate marine carbon and sulfur biogeochemical cycles and play a prominent role in global climate change. As a large DNA virus, EhV has developed a novel "virocell metabolism" model to meet its high metabolic needs. Although it has been widely demonstrated that EhV infection can profoundly rewire lipid metabolism, the epigenetic regulatory mechanisms of lipid metabolism are still obscure. MicroRNAs (miRNAs) can regulate biological pathways by targeting hub genes in the metabolic processes. In this study, the transcriptome, lipidome, and miRNAome were applied to investigate the epigenetic regulation of lipid metabolism in E. huxleyi cells during a detailed time course of viral infection. Combined transcriptomic, lipidomic, and physiological experiments revealed reprogrammed lipid metabolism, along with mitochondrial dysfunction and calcium influx through the cell membrane. A total of 69 host miRNAs (including 1 known miRNA) and 7 viral miRNAs were identified, 27 of which were differentially expressed. Bioinformatic prediction revealed that miRNAs involved in the regulation of lipid metabolism and a dual-luciferase reporter assay suggested that phosphatidylinositol 3-kinase (PI3K) gene might be a target of ehx-miR5. Further qPCR and western blot analysis showed a significant negative correlation between the expression of ehx-miR5 and its target gene PI3K, along with the lower activity of its downstream components (p-Akt, p-TOR, SREBP), indicating that lipid metabolism might be regulated by ehx-miR5 through the PI3K-Akt-TOR signaling pathway. Our findings reveal several novel mechanisms of viral strategies to manipulate host lipid metabolism and provide evidence that ehx-miR5 negatively modulates the expression of PI3K and disturbs lipid metabolism in the interactions between E. huxleyi and EhV.

Comparison of the Main Metabolites in Different Maturation Stages of Camelliavietnamensis Huang Seeds

Camellia vietnamensis Huang is an important woody oil crop in China, which has attracted much attention because of its abundant nutritional components and pharmaceutical value. Its seeds undergo a complex series of physiological and biochemical changes during maturation, with consequent alterations in metabolites. In order to investigate the endogenous metabolism of C. vietnamensis on Hainan Island during seed development, in this study, ultra-high-performance liquid tandem chromatography coupled with quadrupole time-of-flight mass spectrometry (UHPLC/Q-TOF-MS) and multivariate statistical analysis (MSA) were used to analyze the differences in the chemical compounds of C. vietnamensis seeds among the four maturation stages. A total of 293 metabolites were identified from the methanol extract of the seeds of C. vietnamensis. Five metabolites, belonging to benzene and substituted derivatives, 5′-deoxyribonucleosides and linear 1,3-diarylpropanoids, were found in all three comparison groups, with consistently down-regulated trends. The Kyoto Encyclopedia of Genes and Genomes (KEGG) results showed that phloretin and 5′-methylthioadenosine were the differentially expressed metabolites when seeds were in the growth periods of S2 and S3, and indole and L-tryptophan were the differentially expressed metabolites when seeds were in the growth periods of S3 and S4. In addition, 34 flavonoid metabolites were detected, of which 4 were differentially expressed. It was indicated that flavonoids dynamically change during all the oil-tea camellia seed development stages. The findings provide data for the better understanding of endogenous metabolic pathways during C. vietnamensis seed development.

4D genetic networks reveal the genetic basis of metabolites and seed oil-related traits in 398 soybean RILs

Background

The yield and quality of soybean oil are determined by seed oil-related traits, and metabolites/lipids act as bridges between genes and traits. Although there are many studies on the mode of inheritance of metabolites or traits, studies on multi-dimensional genetic network (MDGN) are limited.

Results

In this study, six seed oil-related traits, 59 metabolites, and 107 lipids in 398 recombinant inbred lines, along with their candidate genes and miRNAs, were used to construct an MDGN in soybean. Around 175 quantitative trait loci (QTLs), 36 QTL-by-environment interactions, and 302 metabolic QTL clusters, 70 and 181 candidate genes, including 46 and 70 known homologs, were previously reported to be associated with the traits and metabolites, respectively. Gene regulatory networks were constructed using co-expression, protein–protein interaction, and transcription factor binding site and miRNA target predictions between candidate genes and 26 key miRNAs. Using modern statistical methods, 463 metabolite–lipid, 62 trait–metabolite, and 89 trait–lipid associations were found to be significant. Integrating these associations into the above networks, an MDGN was constructed, and 128 sub-networks were extracted. Among these sub-networks, the gene–trait or gene–metabolite relationships in 38 sub-networks were in agreement with previous studies, e.g., oleic acid (trait)–GmSEI–GmDGAT1a–triacylglycerol (16:0/18:2/18:3), gene and metabolite in each of 64 sub-networks were predicted to be in the same pathway, e.g., oleic acid (trait)–GmPHS–d-glucose, and others were new, e.g., triacylglycerol (16:0/18:1/18:2)–GmbZIP123–GmHD-ZIPIII-10–miR166s–oil content.

Conclusions

This study showed the advantages of MGDN in dissecting the genetic relationships between complex traits and metabolites. Using sub-networks in MGDN, 3D genetic sub-networks including pyruvate/threonine/citric acid revealed genetic relationships between carbohydrates, oil, and protein content, and 4D genetic sub-networks including PLDs revealed the relationships between oil-related traits and phospholipid metabolism likely influenced by the environment. This study will be helpful in soybean quality improvement and molecular biological research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02191-1.

Quality Evaluation of the Oil of Camellia spp

The oil of Camellia spp. has become a well-known high-quality edible oil because of its rich nutrition. It is of great significance to breed fine varieties of Camellia spp. for the sustainable growth of the Camellia spp. industry. This study mainly evaluated the quality and antioxidant capacity of the camellia seed from several sources. The fatty acid composition and main active components of 40 kinds of C. oleifera, C. vietnamensis, C. osmantha, and C. gigantocarpa seeds, and so on, from different regions, were tested using GC-MS and HPLC. The quality of different Camellia spp. germplasm resources was comprehensively evaluated using multiple indices. The unsaturated fatty acid content and the antioxidant capacity of C. vietnamensis from Hainan were higher than those of C. oleifera Abel. In addition, there were a few differences in the fatty acid compositions of Camellia spp. oil from different species. Correlation analysis confirmed that rutin, total saponin, total flavonoids, squalene, and vitamin E were strongly correlated to the antioxidant capacity of Camellia spp. In the comprehensive evaluation, the best quality and strongest antioxidant activity were found for Chengmai Dafeng (C. vietnamensis). These methods in the study were applied for the first time for the quality evaluation of the Camellia spp. species. This study provided new insights into the quality evaluation of the Camellia spp. species, thus facilitating further development of variety breeding along with quality evaluation.

Hormone analysis and candidate genes identification associated with seed size in Camellia oleifera

Camellia oleifera is an important woody oil species in China. Its seed oil has been widely used as a cooking oil. Seed size is a crucial factor influencing the yield of seed oil. In this study, the horizontal diameter, vertical diameter and volume of C. oleifera seeds showed a rapid growth tendency from 235 days after pollination (DAP) to 258 DAP but had a slight increase at seed maturity. During seed development, the expression of genes related to cell proliferation and expansion differ greatly. Auxin plays an important role in C. oleifera seeds; YUC4 and IAA17 were significantly downregulated. Weighted gene co-expression network analysis screened 21 hub transcription factors for C. oleifera seed horizontal diameter, vertical diameter and volume. Among them, SPL4 was significantly decreased and associated with all these three traits, while ABI4 and YAB1 were significantly increased and associated with horizontal diameter of C. oleifera seeds. Additionally, KLU significantly decreased (2040-fold). Collectively, our data advances the knowledge of factors related to seed size and provides a theoretical basis for improving the yield of C. oleifera seeds.

Sugar Metabolism and Transcriptome Analysis Reveal Key Sugar Transporters during Camellia oleifera Fruit Development

Camellia oleifera is a widely planted woody oil crop with economic significance because it does not occupy cultivated land. The sugar-derived acetyl-CoA is the basic building block in fatty acid synthesis and oil synthesis in C. oleifera fruit; however, sugar metabolism in this species is uncharacterized. Herein, the changes in sugar content and metabolic enzyme activity and the transcriptomic changes during C. oleifera fruit development were determined in four developmental stages (CR6: young fruit formation; CR7: expansion; CR9: oil transformation; CR10: ripening). CR7 was the key period of sugar metabolism since it had the highest amount of soluble sugar, sucrose, and glucose with a high expression of genes related to sugar transport (four sucrose transporters (SUTs) or and one SWEET-like gene, also known as a sugar, will eventually be exported transporters) and metabolism. The significant positive correlation between their expression and sucrose content suggests that they may be the key genes responsible for sucrose transport and content maintenance. Significantly differentially expressed genes enriched in the starch and sucrose metabolism pathway were observed in the CR6 versus CR10 stages according to KEGG annotation. The 26 enriched candidate genes related to sucrose metabolism provide a molecular basis for further sugar metabolism studies in C. oleifera fruit.

Integrative iTRAQ-based proteomic and transcriptomic analysis reveals the accumulation patterns of key metabolites associated with oil quality during seed ripening of Camellia oleifera

Camellia oleifera (C. oleifera) is one of the four major woody oil-bearing crops in the world and has relatively high ecological, economic, and medicinal value. Its seeds undergo a series of complex physiological and biochemical changes during ripening, which is mainly manifested as the accumulation and transformation of certain metabolites closely related to oil quality, especially flavonoids and fatty acids. To obtain new insights into the underlying molecular mechanisms, a parallel analysis of the transcriptome and proteome profiles of C. oleifera seeds at different maturity levels was conducted using RNA sequencing (RNA-seq) and isobaric tags for relative and absolute quantification (iTRAQ) complemented with gas chromatography-mass spectrometry (GC-MS) data. A total of 16,530 transcripts and 1228 proteins were recognized with significant differential abundances in pairwise comparisons of samples at various developmental stages. Among these, 317 were coexpressed with a poor correlation, and most were involved in metabolic processes, including fatty acid metabolism, α-linolenic acid metabolism, and glutathione metabolism. In addition, the content of total flavonoids decreased gradually with seed maturity, and the levels of fatty acids generally peaked at the fat accumulation stage; these results basically agreed with the regulation patterns of genes or proteins in the corresponding pathways. The expression levels of proteins annotated as upstream candidates of phenylalanine ammonia-lyase (PAL) and chalcone synthase (CHS) as well as their cognate transcripts were positively correlated with the variation in the flavonoid content, while shikimate O-hydroxycinnamoyltransferase (HCT)-encoding genes had the opposite pattern. The increase in the abundance of proteins and mRNAs corresponding to alcohol dehydrogenase (ADH) was associated with a reduction in linoleic acid synthesis. Using weighted gene coexpression network analysis (WGCNA), we further identified six unique modules related to flavonoid, oil, and fatty acid anabolism that contained hub genes or proteins similar to transcription factors (TFs), such as MADS intervening keratin-like and C-terminal (MIKC_MADS), type-B authentic response regulator (ARR-B), and basic helix-loop-helix (bHLH). Finally, based on the known metabolic pathways and WGCNA combined with the correlation analysis, five coexpressed transcripts and proteins composed of cinnamyl-alcohol dehydrogenases (CADs), caffeic acid 3-O-methyltransferase (COMT), flavonol synthase (FLS), and 4-coumarate: CoA ligase (4CL) were screened out. With this exploratory multiomics dataset, our results presented a dynamic picture regarding the maturation process of C. oleifera seeds on Hainan Island, not only revealing the temporal specific expression of key candidate genes and proteins but also providing a scientific basis for the genetic improvement of this tree species.

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.

Characterization and Cloning of Grape Circular RNAs Identified the Cold Resistance-Related Vv-circATS1

Circular RNAs (circRNAs) are widely distributed and play essential roles in a series of developmental processes, although none have been identified or characterized in grapevines (Vitis vinifera). In this study, we characterized the function of grape circRNA and uncovered thousands of putative back-splicing sites by global transcriptome analysis. Our results indicated that several reported circRNA prediction algorithms should be used simultaneously to obtain comprehensive and reliable circRNA predictions in plants. Furthermore, the length of introns flanking grape circRNAs was closely related to exon circularization. Although the longer introns flanking grape circRNAs appeared to circularize more efficiently, a 20- to 50-nt region seemed large enough to drive grape circRNA biogenesis. In addition, the endogenous introns flanking circularized exon(s) in conjunction with reverse complementary sequences could support the accurate and efficient circularization of various exons in grape, which constitutes a new tool for exploring the functional consequences caused by circRNA expression. Finally, we identified 475 differentially expressed circRNAs in grape leaves under cold stress. Overexpression of Vv-circATS1, a circRNA derived from glycerol-3-P acyltransferase, improved cold tolerance in Arabidopsis (Arabidopsis thaliana), while the linear RNA derived from the same sequence cannot. These results indicate the functional difference between circRNA and linear RNA, and provide new insight into plant abiotic stress resistance.

Small RNA sequencing reveals dynamic microRNA expression of important nutrient metabolism during development of Camellia oleifera fruit

To obtain insight into the function of miRNAs in the synthesis and storage of important nutrients during the development of Camellia oleifera fruit, Illumina sequencing of flower and fruit small-RNA was conducted. The results revealed that 797 miRNAs were significantly differentially expressed between flower and fruit samples of Camellia oleifera. Through integrated GO and KEGG function annotations, it was determined that the miRNA target genes were mainly involved in metabolic pathways, plant hormone signal transduction, fruit development, mitosis and regulation of biosynthetic processes. Carbohydrate accumulation genes were differentially regulated by miR156, miR390 and miR395 in the fruit growth and development process. MiR477 is the key miRNA functioning in regulation of genes and involved in fatty acid synthesis. Additionally, miR156 also has the function of regulating glycolysis and nutrient transformation genes.