Camellia oil extracted from Camellia seeds is rich in unsaturated fatty acids (UFAs) and secondary metabolites beneficial to human health. However, no oil-tea tree genome has yet been published, which is a major obstacle to investigating the heredity improvement of oil-tea trees. Here, using both Illumina and PicBio sequencing technologies, we present the first chromosome-level genome sequence of the oil-tea tree species Camellia chekiangoleosa Hu. (CCH). The assembled genome consists of 15 pseudochromosomes with a genome size of 2.73 Gb and a scaffold N50 of 185.30 Mb. At least 2.16 Gb of the genome assembly consists of repetitive sequences, and the rest involves a high-confidence set of 64 608 protein-coding gene models. Comparative genomic analysis revealed that the CCH genome underwent a whole-genome duplication (WGD) event shared across the Camellia genus at ~57.48 MYA and a γ-WGT event shared across all core eudicot plants at ~120 MYA. Gene family clustering revealed that the genes involved in terpenoid biosynthesis have undergone rapid expansion. Furthermore, we determined the expression patterns of oleic acid accumulation- and terpenoid biosynthesis-associated genes in six tissues. We found that these genes tend to be highly expressed in leaves, pericarp tissues, roots, and seeds. The first chromosome-level genome of oil-tea trees will provide valuable resources for determining Camellia evolution and utilizing the germplasm of this taxon.

Camellia chekiangoleosa has a higher oleic acid content and a shorter reproductive cycle than typical oil tea plants. It was intensively sampled over six C. chekiangoleosa seed development stages. The content of fatty acids determined by GC showed that the accumulation of fatty acids gradually increased from the S1 to S5 stages, and the maximum concentration was reached in S5. Then, fatty acids declined slightly in S6. The main fatty acid component showed the same accumulation trend as the total fatty acids, except linolenic acid, which remained at a low level throughout seed developmental stages. Changes in the expression of fatty acid accumulation-related genes were monitored using second-generation and SMRT full-length transcriptome sequencing. Finally, 18.92 G accurate and reliable data were obtained. Differential expression analysis and weighted coexpression analysis revealed two "gene modules" significantly associated with oleic acid and linoleic acid contents, and the high expression of ENR, KAS I, and KAS II, which accumulate substrates for oleic acid synthesis, was thought to be responsible for the rapid accumulation of fatty acids in the early stage. The rapid increase in fatty acids in the second stage may be closely related to the synergy between the high expression of SAD and low expression of FAD2. In addition, many transcription factors, such as ERF, GRAS, GRF, MADS, MYB and WRKY, may be involved in the fatty acid synthesis. Our data provide a rich resource for further studies on the regulation of fatty acid synthesis in C. chekiangoleosa.

Shell-color polymorphism is a common phenomenon in several mollusk species and has been associated with thermal capacity, developmental stability, shell strength, and immunity. Shell-color polymorphism has been related to the differential expression of genes in several signal transduction pathways; however, the functions of micro-RNAs (miRNAs) in shell-color formation remain unclear. In the present study, we compared high-quality, small-RNA transcriptomes in three strains of the Manila clam Ruditapes philippinarum with specific shell-color patterns, artificially selected for six generations. Totals of 114 known and 208 novel miRNAs were identified by high-throughput sequencing, of which nine known and one novel miRNA were verified by stem-loop quantitative real time-polymerase chain reaction. Predicted miRNA targets were subjected to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. miR-137 and miR-216b and the Hedgehog signaling pathway and Wnt signaling pathway were identified as being potentially involved in pigment formation and regulation in R. philippinarum. These results may help to clarify the role of miRNAs in shell coloration and shed light on the mechanisms regulating color formation in bivalve shells.

Red-flowered oil-tea camellia (ROC) is an important woody oil species growing in the south, and its oil has high nutritional value. There are four main species of ROC in China, namely, Camellia chekiangoleosa (CCH), Camellia polyodonta (CPO), Camellia semiserrata (CSE) and Camellia reticulata (CRE). Reports on the comprehensive comparative analysis of ROC are limited. This study investigated the fruit characteristics and nutritional components of four ROC fruits, and the results showed that ROC had high oil content with levels of 39.13%–58.84%, especially the CCH fruit, which reached 53.6–58.84%. The contents of lipid concomitants of ROC oil were also substantial, including β-amyrin (0.87 mg/g–1.41 mg/g), squalene (0.43 mg/g–0.69 mg/g), β-sitosterin (0.47 mg/g–0.63 mg/g) and α-tocopherol (177.52 μg/g–352.27 μg/g). Moreover, the transverse diameter(TD)/longitudinal diameter (LD) of fruits showed a significant positive correlation with the oil content, and ROC fruits with thinner peels seemed to have better oil quality, which is similar to the result of the oil quality evaluation obtained by the gray correlation coefficient evaluation method. Four ROC oils were evaluated using the gray correlation coefficient method based on 11 indicators related to the nutritional value of ROC. CCH oil had the highest score of 0.8365, and YS-2 (a clone of CCH) was further evaluated as the best CCH oil. Finally, the results of heatmap analysis showed that triglycerides could be used as a characteristic substance to distinguish CCH oil from the other three ROC oils. The PLSDA (Partial least squares regression analysis) model and VIP (Variable important in projection) values further showed that P/S/O, P/O/O, P/L/L, P/L/Ln, S/S/O, S/O/O and P/S/S (these all represent abbreviations for fatty acids) could be used as characteristic differential triglycerides among the four ROC oils. This study provides a convenient way for planters to assess the nutritional quality of seed oil depending on fruit morphology and a potential way to distinguish between various ROC oils.

MADS-box genes encode transcription factors that affect plant growth and development. Camellia chekiangoleosa is an oil tree species with ornamental value, but there have been few molecular biological studies on the developmental regulation of this species. To explore their possible role in C. chekiangoleosa and lay a foundation for subsequent research, 89 MADS-box genes were identified across the whole genome of C. chekiangoleosa for the first time. These genes were present on all the chromosomes and were found to have expanded by tandem duplication and fragment duplication. Based on the results of a phylogenetic analysis, the 89 MADS-box genes could be divided into either type I (38) or type II (51). Both the number and proportion of the type II genes were significantly greater than those of Camellia sinensis and Arabidopsis thaliana, indicating that C. chekiangoleosa type II genes experienced a higher duplication rate or a lower loss rate. The results of both a sequence alignment and a conserved motif analysis suggest that the type II genes are more conserved, meaning that they may have originated and differentiated earlier than the type I genes did. At the same time, the presence of extra-long amino acid sequences may be an important feature of C. chekiangoleosa. Gene structure analysis revealed the number of introns of MADS-box genes: twenty-one type I genes had no introns, and 13 type I genes contained only 1~2 introns. The type II genes have far more introns and longer introns than the type I genes do. Some MIKC^C genes have super large introns (≥15 kb), which are rare in other species. The super large introns of these MIKC^C genes may indicate richer gene expression. Moreover, the results of a qPCR expression analysis of the roots, flowers, leaves and seeds of C. chekiangoleosa showed that the MADS-box genes were expressed in all those tissues. Overall, compared with that of the type I genes, the expression of the type II genes was significantly higher. The CchMADS31 and CchMADS58 genes (type II) were highly expressed specifically in the flowers, which may in turn regulate the size of the flower meristem and petals. CchMADS55 was expressed specifically in the seeds, which might affect seed development. This study provides additional information for the functional characterization of the MADS-box gene family and lays an important foundation for in-depth study of related genes, such as those involved in the development of the reproductive organs of C. chekiangoleosa.