Camellia oleifera is a major tree species for producing edible oil. Its seed oil is well known for the high level of oleic acids; however, little is known regarding the molecular mechanism of lipid biosynthesis in C. oleifera. Here, we measured the oil contents and fatty acid (FA) compositions at four developmental stages and investigated the global gene expression profiles through transcriptomics sequencing. We identified differentially-expressed genes (DEGs) among the developmental stages and found that the distribution of numbers of DEGs was associated with the accumulation pattern of seed oil. Gene Ontology (GO) enrichment analysis revealed some critical biological processes related to oil accumulation, including lipid metabolism and phosphatidylcholine metabolism. Furthermore, we investigated the expression patterns of lipid biosynthesis genes. We showed that most of the genes were identified with single or multiple copies, and some had correlated profiles along oil accumulation. We proposed that the higher levels of stearoyl-ACP desaturases (SADs) coupled with lower activities of fatty acid desaturase 2 (FAD2) might be responsive to the boost of oleic acid at the late stage of C. oleifera seeds’ development. This work presents a comprehensive transcriptomics study of C. oleifera seeds and uncovers valuable DEGs that are associated with the seed oil accumulation.

Camellia oleifera (Ca. oleifera) is a woody tree species cultivated for the production of edible oil from its seed. The growth and yield of tea-oil trees are severely affected by anthracnose (caused by Colletotrichum gloeosporioides). In this study, the transcriptomic and metabolomic analyses were performed to detect the key transcripts and metabolites associated with differences in the susceptibility between anthracnose-resistant (ChangLin150) and susceptible (ChangLin102) varieties of Ca. oleifera. In total, 5001 differentially expressed genes (DEGs) were obtained, of which 479 DEGs were common between the susceptible and resistant varieties and further analyzed. KEGG enrichment analysis showed that these DEGs were significantly enriched in tyrosine metabolism, phenylpropanoid biosynthesis, flavonoid biosynthesis and isoquinoline alkaloid biosynthesis pathways. Furthermore, 68 differentially accumulated metabolites (DAMs) were detected, including flavonoids, such as epicatechin, phenethyl caffeate and procyanidin B2. Comparison of the DEGs and DAMs revealed that epicatechin, procyanidin B2 and arachidonic acid (peroxide free) are potentially important. The expression patterns of genes involved in flavonoid biosynthesis were confirmed by qRT-PCR. These results suggested that flavonoid biosynthesis might play an important role in the fight against anthracnose. This study provides valuable molecular information about the response of Ca. oleifera to Co. gloeosporioides infection and will aid the selection of resistant varieties using marker-assisted breeding.

With the increasing power of DNA sequencing, the genomics-based approach is becoming a promising resolution to dissect the molecular mechanism of domestication of complex traits in trees. Genus Camellia possesses rich resources with a substantial value for producing beverage, ornaments, edible oil and more. Currently, a vast number of genetic and genomic research studies in Camellia plants have emerged and provided an unprecedented opportunity to expedite the molecular breeding program. In this paper, we summarize the recent advances of gene expression and genomic resources in Camellia species and focus on identifying genes related to key economic traits such as flower and fruit development and stress tolerances. We investigate the genetic alterations and genomic impacts under different selection programs in closely related species. We discuss future directions of integrating large-scale population and quantitative genetics and multiple omics to identify key candidates to accelerate the breeding process. We propose that future work of exploiting the genomic data can provide insights related to the targets of domestication during breeding and the evolution of natural trait adaptations in genus Camellia.

Camellia is the largest genus in the family Theaceae. Due to phenotypic diversity, frequent hybridization, and polyploidization, an understanding of the phylogenetic relationships between Camellia species remains challenging. Comparative chloroplast (cp) genomics provides an informative resource for phylogenetic analyses of Camellia. In this study, 12 chloroplast genome sequences from nine Camellia species were determined using Illumina sequencing technology via de novo assembly. The cp genome sizes ranged from 156,545 to 157,021 bp and were organized into quadripartite regions with the typical angiosperm cp genomes. Each genome harbored 87 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes in the same order and orientation. Differences in long and short sequence repeats, SNPs, and InDels were detected across the 12 cp genomes. Combining with the complete cp sequences of seven other species in the genus Camellia, a total of nine intergenic sequence divergent hotspots and 14 protein-coding genes with high sequence polymorphism were identified. These hotspots, especially the InDel (~400 bp) located in atpH-atpI region, had sufficient potential to be used as barcode markers for further phylogenetic analysis and species identification. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These cp genomes could facilitate the development of new molecular markers, accurate species identification, and investigations of the phylogenomic relationships of the genus Camellia.