The genome of oil-Camellia and population genomics analysis provide insights into seed oil domestication

Genome-wide identification and molecular evolution of Dof gene family in Camellia oleifera

DNA binding with one finger(Dof) gene family is a class of transcription factors which play an important role on plant growth and development. Genome-wide identification results indicated that there were 45 Dof genes(ColDof) in C.oleifera genome. All 45 ColDof proteins were non-transmembrane and non-secretory proteins. Phosphorylation site analysis showed that biological function of ColDof proteins were mainly realized by phosphorylation at serine (Ser) site. The secondary structure of 44 ColDof proteins was dominated by random coil, and only one ColDof protein was dominated by α-helix. ColDof genes' promoter region contained a variety of cis-acting elements, including light responsive regulators, gibberellin responsive regulators, abscisic acid responsive regulators, auxin responsive regulators and drought induction responsive regulators. The SSR sites analysis showed that the proportion of single nucleotide repeats and the frequency of A/T in ColDof genes were the largest. Non-coding RNA analysis showed that 45 ColDof genes contained 232 miRNAs. Transcription factor binding sites of ColDof genes showed that ColDof genes had 5793 ERF binding sites, 4381 Dof binding sites, 2206 MYB binding sites, 3702 BCR-BPC binding sites. ColDof9, ColDof39 and ColDof44 were expected to have the most TFBSs. The collinearity analysis showed that there were 40 colinear locis between ColDof proteins and AtDof proteins. Phylogenetic analysis showed that ColDof gene family was most closely related to that of Camellia sinensis var. sinensis cv.Biyun and Camellia lanceoleosa. Protein-protein interaction analysis showed that ColDof34, ColDof20, ColDof28, ColDof35, ColDof42 and ColDof26 had the most protein interactions. The transcriptome analysis of C. oleifera seeds showed that 21 ColDof genes were involved in the growth and development process of C. oleifera seeds, and were expressed in 221 C. oleifera varieties. The results of qRT-PCR experiments treated with different concentrations NaCl and PEG6000 solutions indicated that ColDof1, ColDof2, ColDof14 and ColDof36 not only had significant molecular mechanisms for salt stress tolerance, but also significant molecular functions for drought stress tolerance in C. oleifera. The results of this study provide a reference for further understanding of the function of ColDof genes in C.oleifera.

Genomic variation of 363 diverse tea accessions unveils the genetic diversity, domestication, and structural variations associated with tea adaptation

Domestication has shaped the population structure and agronomic traits of tea plants, yet the complexity of tea population structure and genetic variation that determines these traits remains unclear. We here investigated the resequencing data of 363 diverse tea accessions collected extensively from almost all tea distributions and found that the population structure of tea plants was divided into eight subgroups, which were basically consistent with their geographical distributions. The genetic diversity of tea plants in China decreased from southwest to east as latitude increased. Results also indicated that Camellia sinensis var. assamica (CSA) illustrated divergent selection signatures with Camellia sinensis var. sinensis (CSS). The domesticated genes of CSA were mainly involved in leaf development, flavonoid and alkaloid biosynthesis, while the domesticated genes in CSS mainly participated in amino acid metabolism, aroma compounds biosynthesis, and cold stress. Comparative population genomics further identified ~730 Mb novel sequences, generating 6,058 full-length protein-encoding genes, significantly expanding the gene pool of tea plants. We also discovered 217,376 large-scale structural variations and 56,583 presence and absence variations (PAVs) across diverse tea accessions, some of which were associated with tea quality and stress resistance. Functional experiments demonstrated that two PAV genes (CSS0049975 and CSS0006599) were likely to drive trait diversification in cold tolerance between CSA and CSS tea plants. The overall findings not only revealed the genetic diversity and domestication of tea plants, but also underscored the vital role of structural variations in the diversification of tea plant traits.

Sugar import mediated by sugar transporters and cell wall invertases for seed development in Camellia oleifera

Seed development and yield depend on the transport and supply of sugar. However, an insufficient supply of nutrients from maternal tissues to embryos results in seed abortion and yield reduction in Camellia oleifera. In this study, we systematically examined the route and regulatory mechanisms of sugar import into developing C. oleifera seeds using a combination of histological observations, transcriptome profiling, and functional analysis. Labelling with the tracer carboxyfluorescein revealed a symplasmic route in the integument and an apoplasmic route for postphloem transport at the maternal-filial interface. Enzymatic activity and histological observation showed that at early stages [180-220 days after pollination (DAP)] of embryo differentiation, the high hexose/sucrose ratio was primarily mediated by acid invertases, and the micropylar endosperm/suspensor provides a channel for sugar import. Through Camellia genomic profiling, we identified three plasma membrane-localized proteins including CoSWEET1b, CoSWEET15, and CoSUT2 and one tonoplast-localized protein CoSWEET2a in seeds and verified their ability to transport various sugars via transformation in yeast mutants and calli. In situ hybridization and profiling of glycometabolism-related enzymes further demonstrated that CoSWEET15 functions as a micropylar endosperm-specific gene, together with the cell wall acid invertase CoCWIN9, to support early embryo development, while CoSWEET1b, CoSWEET2a, and CoSUT2 function at transfer cells and chalazal nucellus coupled with CoCWIN9 and CoCWIN11 responsible for sugar entry in bulk into the filial tissue. Collectively, our findings provide the first comprehensive evidence of the molecular regulation of sugar import into and within C. oleifera seeds and provide a new target for manipulating seed development.

The complex hexaploid oil-Camellia genome traces back its phylogenomic history and multi-omics analysis of Camellia oil biosynthesis

Oil-Camellia (Camellia oleifera), belonging to the Theaceae family Camellia, is an important woody edible oil tree species. The Camellia oil in its mature seed kernels, mainly consists of more than 90% unsaturated fatty acids, tea polyphenols, flavonoids, squalene and other active substances, which is one of the best quality edible vegetable oils in the world. However, genetic research and molecular breeding on oil-Camellia are challenging due to its complex genetic background. Here, we successfully report a chromosome-scale genome assembly for a hexaploid oil-Camellia cultivar Changlin40. This assembly contains 8.80 Gb genomic sequences with scaffold N50 of 180.0 Mb and 45 pseudochromosomes comprising 15 homologous groups with three members each, which contain 135 868 genes with an average length of 3936 bp. Referring to the diploid genome, intragenomic and intergenomic comparisons of synteny indicate homologous chromosomal similarity and changes. Moreover, comparative and evolutionary analyses reveal three rounds of whole-genome duplication (WGD) events, as well as the possible diversification of hexaploid Changlin40 with diploid occurred approximately 9.06 million years ago (MYA). Furthermore, through the combination of genomics, transcriptomics and metabolomics approaches, a complex regulatory network was constructed and allows to identify potential key structural genes (SAD, FAD2 and FAD3) and transcription factors (AP2 and C2H2) that regulate the metabolism of Camellia oil, especially for unsaturated fatty acids biosynthesis. Overall, the genomic resource generated from this study has great potential to accelerate the research for the molecular biology and genetic improvement of hexaploid oil-Camellia, as well as to understand polyploid genome evolution.

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.

De novo assembling a high-quality genome sequence of Amur grape (Vitis amurensis Rupr.) gives insight into Vitis divergence and sex determination

To date, there has been no high-quality sequence for genomes of the East Asian grape species, hindering biological and breeding efforts to improve grape cultivars. This study presents ~522 Mb of the Vitis amurensis (Va) genome sequence containing 27 635 coding genes. Phylogenetic analysis indicated that Vitis riparia (Vr) may have first split from the other two species, Va and Vitis vinifera (Vv). Divergent numbers of duplicated genes reserved among grapes suggests that the core eudicot-common hexaploidy (ECH) and the subsequent genome instability still play a non-negligible role in species divergence and biological innovation. Prominent accumulation of sequence variants might have improved cold resistance in Va, resulting in a more robust network of regulatory cold resistance genes, explaining why it is extremely cold-tolerant compared with Vv and Vr. In contrast, Va has preserved many fewer nucleotide binding site (NBS) disease resistance genes than the other grapes. Notably, multi-omics analysis identified one trans-cinnamate 4-monooxygenase gene positively correlated to the resveratrol accumulated during Va berry development. A selective sweep analysis revealed a hypothetical Va sex-determination region (SDR). Besides, a PPR-containing protein-coding gene in the hypothetical SDR may be related to sex determination in Va. The content and arrangement order of genes in the putative SDR of female Va were similar to those of female Vv. However, the putative SDR of female Va has lost one flavin-containing monooxygenase (FMO) gene and contains one extra protein-coding gene uncharacterized so far. These findings will improve the understanding of Vitis biology and contribute to the improvement of grape breeding.

Comparative Genome-Wide Identification of the Fatty Acid Desaturase Gene Family in Tea and Oil Tea

Camellia oil is valuable as an edible oil and serves as a base material for a range of high-value products. Camellia plants of significant economic importance, such as Camellia sinensis and Camellia oleifera, have been classified into sect. Thea and sect. Oleifera, respectively. Fatty acid desaturases play a crucial role in catalyzing the formation of double bonds at specific positions of fatty acid chains, leading to the production of unsaturated fatty acids and contributing to lipid synthesis. Comparative genomics results have revealed that expanded gene families in oil tea are enriched in functions related to lipid, fatty acid, and seed processes. To explore the function of the FAD gene family, a total of 82 FAD genes were identified in tea and oil tea. Transcriptome data showed the differential expression of the FAD gene family in mature seeds of tea tree and oil tea tree. Furthermore, the structural analysis and clustering of FAD proteins provided insights for the further exploration of the function of the FAD gene family and its role in lipid synthesis. Overall, these findings shed light on the role of the FAD gene family in Camellia plants and their involvement in lipid metabolism, as well as provide a reference for understanding their function in oil synthesis.

Genome assembly of the rare and endangered Grantham's camellia, Camellia granthamiana

Grantham's camellia (Camellia granthamiana Sealy) is a rare and endangered tea species discovered in Hong Kong in 1955 and endemic to southern China. Despite its high conservation value, the genomic resources of C. granthamiana are limited. Here, we present a chromosome-scale draft genome of the tetraploid C. granthamiana (2n = 4x = 60), combining PacBio long-read sequencing and Omni-C data. The assembled genome size is ∼2.4 Gb, with most sequences anchored to 15 pseudochromosomes resembling a monoploid genome. The genome has high contiguity, with a scaffold N50 of 139.7 Mb, and high completeness (97.8% BUSCO score). Our gene model prediction resulted in 68,032 protein-coding genes (BUSCO score of 90.9%). We annotated 1.65 Gb of repeat content (68.48% of the genome). Our Grantham's camellia genome assembly is a valuable resource for investigating Grantham's camellia's biology, ecology, and phylogenomic relationships with other Camellia species, and provides a foundation for further conservation measures.

Genome-wide identification and comparative analysis of YABBY transcription factors in oil tea and tea tree

The plant-specific transcription factor gene family, YABBY, plays an important role in plant development and stress response. Although YABBY genes have been identified in numerous species, a comprehensive characterization of YABBYs in tea tree and oil tea has been lacking. In this study, ten and three YABBY genes were identified in Camellia sinensis and C. oleifera, respectively. YABBY proteins could be divided into five subfamilies. Most YABBY genes in the same clade had similar structures and conserved motifs. Protein evolutionary analysis revealed that FIL/YAB3 displayed high conservation in all positions, followed by INO, YAB2, YAB5, and CRC. Specific site analysis suggested that the YABBY family was polyphyletic during the evolution. Compared to C. oleifera, two segmentally duplicated gene pairs were formed in C. sinensis during recent WGD events generated 30.69 and 45.08 Mya, respectively. Cis-acting element indicated that most YABBY genes contain box4, ARE, and MYB elements. A total of 120 SSR loci were found within CsYABBYs, consisting of six types, while 48 SSR loci were identified within CoYABBY, consisting of three types. Transcriptome results revealed that CRC and INO clades were specifically expressed in floral organs. The expression of CsYABBY10 and CsYABBY5 was significantly up-regulated under drought and salt treatments, respectively, as confirmed by qRT-PCR. CoYABBY genes were more susceptible to salt stress, as CoYABBY3 increased by about 15-fold. Furthermore, functional differentiation may have occurred in duplicated genes. These discoveries provide important information for further research on YABBYs in tea tree and oil tea.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03940-9.

The Idesia polycarpa genome provides insights into its evolution and oil biosynthesis

The deciduous tree Idesia polycarpa can provide premium edible oil with high polyunsaturated fatty acid contents. Here, we generate its high-quality reference genome, which is ∼1.21 Gb, comprising 21 pseudochromosomes and 42,086 protein-coding genes. Phylogenetic and genomic synteny analyses show that it diverged with Populus trichocarpa about 16.28 million years ago. Notably, most fatty acid biosynthesis genes are not only increased in number in its genome but are also highly expressed in the fruits. Moreover, we identify, through genome-wide association analysis and RNA sequencing, the I. polycarpa SUGAR TRANSPORTER 5 (IpSTP5) gene as a positive regulator of high oil accumulation in the fruits. Silencing of IpSTP5 by virus-induced gene silencing causes a significant reduction of oil content in the fruits, suggesting it has the potential to be used as a molecular marker to breed the high-oil-content cultivars. Our results collectively lay the foundation for breeding the elite cultivars of I. polycarpa.

Genomics insights into flowering and floral pattern formation: regional duplication and seasonal pattern of gene expression in Camellia

BACKGROUND

The formation and domestication of ornamental traits are influenced by various aspects, such as the recognition of esthetic values and cultural traditions. Camellia japonica is widely appreciated and domesticated around the world mainly due to its rich variations in ornamental traits. Ornamental camellias have a diverse range of resources, including different bud variations from Camellia spp. as well as inter- and intra- specific hybridization. Despite research on the formation of ornamental traits, a basic understanding of their genetics and genomics is still lacking.

RESULTS

Here, we report the chromosomal-level reference genome of C. japonica through combining multiple DNA-sequencing technologies and obtain a high-density genetic linkage map of 4255 markers by sequencing 98 interspecific F1 hybrids between C. japonica and C. chekiangoleosa. We identify two whole-genome duplication events in C. japonica: one is a shared ancient γ event, and the other is revealed to be specific to genus Camellia. Based on the micro-collinearity analysis, we find large-scale segmental duplication of chromosome 8, resulting to two copies of the AGAMOUS loci, which may play a key role in the domestication of floral shapes. To explore the regulatory mechanisms of seasonal flowering, we have analyzed year-round gene expression patterns of C. japonica and C. azalea-a sister plant of continuous flowering that has been widely used for cross breeding. Through comparative analyses of gene co-expression networks and annual gene expression patterns, we show that annual expression rhythms of some important regulators of seasonal growth and development, including GIGANTEA and CONSTANS of the photoperiod pathway, have been disrupted in C. azalea. Furthermore, we reveal that the distinctive expression patterns of FLOWERING LOCUS T can be correlated with the seasonal activities of flowering and flushing. We demonstrate that the regulatory module involved in GIGANTEA, CONSTANS, and FLOWERING LOCUS T is central to achieve seasonality.

CONCLUSIONS

Through the genomic and comparative genomics characterizations of ornamental Camellia spp., we propose that duplication of chromosomal segments as well as the establishment of gene expression patterns has played a key role in the formation of ornamental traits (e.g., flower shape, flowering time). This work provides a valuable genomic platform for understanding the molecular basis of ornamental traits.

Single-cell RNA-seq reveals a link of ovule abortion and sugar transport in Camellia oleifera

Camellia oleifera is the most important woody oil crop in China. Seed number per fruit is an important yield trait in C. oleifera. Ovule abortion is generally observed in C. oleifera and significantly decreases the seed number per fruit. However, the mechanisms of ovule abortion remain poorly understood at present. Single-cell RNA sequencing (scRNA-seq) was performed using mature ovaries of two C. oleifera varieties with different ovule abortion rates (OARs). In total, 20,526 high-quality cells were obtained, and 18 putative cell clusters were identified. Six cell types including female gametophyte, protoxylem, protophloem, procambium, epidermis, and parenchyma cells were identified from three main tissue types of ovule, placenta, and pericarp inner layer. A comparative analysis on scRNA-seq data between high- and low-OAR varieties demonstrated that the overall expression of CoSWEET and CoCWINV in procambium cells, and CoSTP in the integument was significantly upregulated in the low-OAR variety. Both the infertile ovule before pollination and the abortion ovule producing after compatible pollination might be attributed to selective abortion caused by low sugar levels in the apoplast around procambium cells and a low capability of hexose uptake in the integument. Here, the first single-cell transcriptional landscape is reported in woody crop ovaries. Our investigation demonstrates that ovule abortion may be related to sugar transport in placenta and ovules and sheds light on further deciphering the mechanism of regulating sugar transport and the improvement of seed yield in C. oleifera.

Decoding the complete organelle genomic architecture of Stewartia gemmata: an early-diverging species in Theaceae

BACKGROUND

Theaceae, comprising 300 + species, holds significance in biodiversity, economics, and culture, notably including the globally consumed tea plant. Stewartia gemmata, a species of the earliest diverging tribe Stewartieae, is critical to offer insights into Theaceae's origin and evolutionary history.

RESULT

We sequenced the complete organelle genomes of Stewartia gemmata using short/long reads sequencing technologies. The chloroplast genome (158,406 bp) exhibited a quadripartite structure including the large single-copy region (LSC), a small single-copy region (SSC), and a pair of inverted repeat regions (IRs); 114 genes encoded 80 proteins, 30 tRNAs, and four rRNAs. The mitochondrial genome (681,203 bp) exhibited alternative conformations alongside a monocyclic structure: 61 genes encoding 38 proteins, 20 tRNAs, three rRNAs, and RNA editing-impacting genes, including ATP6, RPL16, COX2, NAD4L, NAD5, NAD7, and RPS1. Comparative analyses revealed frequent recombination events and apparent rRNA gene gains and losses in the mitochondrial genome of Theaceae. In organelle genomes, the protein-coding genes exhibited a strong A/U bias at codon endings; ENC-GC3 analysis implies selection-driven codon bias. Transposable elements might facilitate interorganelle sequence transfer. Phylogenetic analysis confirmed Stewartieae's early divergence within Theaceae, shedding light on organelle genome characteristics and evolution in Theaceae.

CONCLUSIONS

We studied the detailed characterization of organelle genomes, including genome structure, composition, and repeated sequences, along with the identification of lateral gene transfer (LGT) events and complexities. The discovery of a large number of repetitive sequences and simple sequence repeats (SSRs) has led to new insights into molecular phylogenetic markers. Decoding the Stewartia gemmata organellar genome provides valuable genomic resources for further studies in tea plant phylogenomics and evolutionary biology.

TPIA2: an updated tea plant information archive for Camellia genomics

The genus Camellia consists of about 200 species, which include many economically important species widely used for making tea, ornamental flowers and edible oil. Here, we present an updated tea plant information archive for Camellia genomics (TPIA2; http://tpia.teaplants.cn) by integrating more novel large-scale genomic, transcriptomic, metabolic and genetic variation datasets as well as a variety of useful tools. Specifically, TPIA2 hosts all currently available and well assembled 10 Camellia genomes and their comprehensive annotations from three major sections of Camellia. A collection of 15 million SNPs and 950 950 small indels from large-scale genome resequencing of 350 diverse tea accessions were newly incorporated, followed by the implementation of a novel 'Variation' module to facilitate data retrieval and analysis of the functionally annotated variome. Moreover, 116 Camellia transcriptomes were newly assembled and added, leading to a significant extension of expression profiles of Camellia genes to 13 developmental stages and eight abiotic/biotic treatments. An updated 'Expression' function has also been implemented to provide a comprehensive gene expression atlas for Camellia. Two novel analytic tools (e.g. Gene ID Convert and Population Genetic Analysis) were specifically designed to facilitate the data exchange and population genomics in Camellia. Collectively, TPIA2 provides diverse updated valuable genomic resources and powerful functions, and will continue to be an important gateway for functional genomics and population genetic studies in Camellia.

A chromosome-level genome assembly provides insights into Cornus wilsoniana evolution, oil biosynthesis, and floral bud development

Cornus wilsoniana W. is a woody oil plant with high oil content and strong hypolipidemic effects, making it a valuable species for medicinal, landscaping, and ecological purposes in China. To advance genetic research on this species, we employed PacBio together with Hi-C data to create a draft genome assembly for C. wilsoniana. Based on an 11-chromosome anchored chromosome-level assembly, the estimated genome size was determined to be 843.51 Mb. The N50 contig size and N50 scaffold size were calculated to be 4.49 and 78.00 Mb, respectively. Furthermore, 30 474 protein-coding genes were annotated. Comparative genomics analysis revealed that C. wilsoniana diverged from its closest species ~12.46 million years ago (Mya). Furthermore, the divergence between Cornaceae and Nyssaceae occurred >62.22 Mya. We also found evidence of whole-genome duplication events and whole-genome triplication γ, occurring at ~44.90 and 115.86 Mya. We further inferred the origins of chromosomes, which sheds light on the complex evolutionary history of the karyotype of C. wilsoniana. Through transcriptional and metabolic analysis, we identified two FAD2 homologous genes that may play a crucial role in controlling the oleic to linoleic acid ratio. We further investigated the correlation between metabolites and genes and identified 33 MADS-TF homologous genes that may affect flower morphology in C. wilsoniana. Overall, this study lays the groundwork for future research aimed at identifying the genetic basis of crucial traits in C. wilsoniana.

Oilomics: An important branch of foodomics dealing with oil science and technology

Oil is one of three nutritious elements. The application of omics techniques in the field of oil science and technology is attracted increasing attention. Oilomics, which emerged as an important branch of foodomics, has been widely used in various aspects of oil science and technology. However, there are currently no articles systematically reviewing the application of oilomics. This paper aims to provide a critical overview of the advantages and value of oilomics technology compared to traditional techniques in various aspects of oil science and technology, including oil nutrition, oil processing, oil quality, safety, and traceability. Moreover, this article intends to review major issues in oilomics and give a comprehensive, critical overview of the current state of the art, future challenges and trends in oilomics, with a view to promoting the optimal application and development of oilomics technology in oil science and technology.

Extensive natural Agrobacterium-induced transformation in the genus Camellia

MAIN CONCLUSION

The genus Camellia underwent extensive natural transformation by Agrobacterium. Over a period of 15 million years, at least 12 different inserts accumulated in 72 investigated Camellia species. Like a wide variety of other wild and cultivated plants, Camellia species carry cellular T-DNA sequences (cT-DNAs) in their nuclear genomes, resulting from natural Agrobacterium-mediated transformation. Short and long DNA sequencing reads of 435 accessions belonging to 72 Camellia species (representing 12 out of 14 sections) were investigated for the occurrence of cT-DNA insertions. In all, 12 different cT-DNAs were recovered, either completely or partially, called CaTA to CaTL. Divergence analysis of internal cT-DNA repeats revealed that the insertion events span a period from 0.075 to 15 Mio years ago, and yielded an average transformation frequency of one event per 1.25 Mio years. The two oldest inserts, CaTA and CaTD, have been modified by spontaneous deletions and inversions, and by insertion of various plant sequences. In those cases where enough accessions were available (C. japonica, C. oleifera, C. chekiangoleosa, C. sasanqua and C. pitardii), the younger cT-DNA inserts showed a patchy distribution among different accessions of each species, indicating that they are not genetically fixed. It could be shown that Camellia breeding has led to intersectional transfer of cT-DNAs. Altogether, the cT-DNAs cover 374 kb, and carry 47 open reading frames (ORFs). Two Camellia cT-DNA genes, CaTH-orf358 and CaTK-orf8, represent new types of T-DNA genes. With its large number of cT-DNA sequences, the genus Camellia constitutes an interesting model for the study of natural Agrobacterium transformants.

Ectopic expression of Camellia oleifera Abel. gibberellin 20-oxidase gene increased plant height and promoted secondary cell walls deposition in Arabidopsis

MAIN CONCLUSION

Ectopic expression of Camellia oleifera Abel. gibberellin 20-oxidase 1 caused a taller phenotype, promoted secondary cell wall deposition, leaf enlargement, and early flowering, and reduced chlorophyll and anthocyanin accumulation and seed enlargement phenotype in Arabidopsis. Plant height and secondary cell wall (SCW) deposition are important plant traits. Gibberellins (GAs) play important roles in regulating plant height and SCWs deposition. Gibberellin 20-oxidase (GA20ox) is an important enzyme involved in GA biosynthesis. In the present study, we identified a GA synthesis gene in Camellia oleifera. The total length of the CoGA20ox1 gene sequence was 1146 bp, encoding 381 amino acids. Transgenic plants with CoGA20ox1 had a taller phenotype; a seed enlargement phenotype; promoted SCWs deposition, leaf enlargement, and early flowering; and reduced chlorophyll and anthocyanin accumulation. Genetic analysis showed that the mutant ga20ox1-3 Arabidopsis partially rescued the phenotype of CoGA20ox1 overexpression plants. The results showed that CoGA20ox1 participates in the growth and development of C. oleifera. The morphological changes in CoGA20ox1 overexpressed plants provide a theoretical basis for further exploration of GA biosynthesis and analysis of the molecular mechanism in C. oleifera.

Application of Multi-Perspectives in Tea Breeding and the Main Directions

Tea plants are an economically important crop and conducting research on tea breeding contributes to enhancing the yield and quality of tea leaves as well as breeding traits that satisfy the requirements of the public. This study reviews the current status of tea plants germplasm resources and their utilization, which has provided genetic material for the application of multi-omics, including genomics and transcriptomics in breeding. Various molecular markers for breeding were designed based on multi-omics, and available approaches in the direction of high yield, quality and resistance in tea plants breeding are proposed. Additionally, future breeding of tea plants based on single-cellomics, pangenomics, plant-microbe interactions and epigenetics are proposed and provided as references. This study aims to provide inspiration and guidance for advancing the development of genetic breeding in tea plants, as well as providing implications for breeding research in other crops.

Development of a protoplast isolation system for functional gene expression and characterization using petals of Camellia Oleifera

Protoplasts preparation and purification have been frequently used in plant genetics and breeding studies, whereas application of protoplasts in woody plants is still in its infancy. Although transient gene expression using purified protoplasts is well-documented and widely used in model plants and agriculture crops, no instance of either stable transformation or transient gene expression in the woody plant Camellia Oleifera has as of yet been reported. Here, we developed a protoplast preparation and purification method using C. oleifera petals by optimizing osmotic condition with D-mannitol and polysaccharide-degrading enzyme concentrations for petal cell wall digestion, to reach a high efficiency of protoplast productivity and viability. The achieved protoplasts yield was approximately 1.42 × 107 cells per gram of petal material and the viability of protoplasts was up to 89%. In addition, we explored influencing factors of protoplast transformation, including concentrations of PEG4000 and plasmid DNA. The transformation efficiency of 81% could be reached under the optimized condition. This protoplast isolation and transient expression system were deployed to further identify the functional regulation of C. oleifera related genes and the subcellular distribution of their encoded products. In summary, the protoplast isolation and transient expression system we established using oil-tea tree petals is an efficient, versatile and time-saving system, being suitable for gene function characterization and molecular mechanism analysis.

Whole-genome resequencing of the native sheep provides insights into the microevolution and identifies genes associated with reproduction traits

BACKGROUND

Sheep genomes undergo numerous genes losses, gains and mutation that generates genome variability among breeds of the same species after long time natural and artificial selection. However, the microevolution of native sheep in northwest China remains elusive. Our aim was to compare the genomes and relevant reproductive traits of four sheep breeds from different climatic environments, to unveil the selection challenges that this species cope with, and the microevolutionary differences in sheep genomes. Here, we resequenced the genomes of 4 representative sheep breeds in northwest China, including Kazakh sheep and Duolang sheep of native breeds, and Hu sheep and Suffolk sheep of exotic breeds with different reproductive characteristics.

RESULTS

We found that these four breeds had a similar expansion experience from ~ 10,000 to 1,000,000 years ago. In the past 10,000 years, the selection intensity of the four breeds was inconsistent, resulting in differences in reproductive traits. We explored the sheep variome and selection signatures by F_ST and θπ. The genomic regions containing genes associated with different reproductive traits that may be potential targets for breeding and selection were detected. Furthermore, non-synonymous mutations in a set of plausible candidate genes and significant differences in their allele frequency distributions across breeds with different reproductive characteristics were found. We identified PAK1, CYP19A1 and PER1 as a likely causal gene for seasonal reproduction in native sheep through qPCR, Western blot and ELISA analyses. Also, the haplotype frequencies of 3 tested gene regions related to reproduction were significantly different among four sheep breeds.

CONCLUSIONS

Our results provide insights into the microevolution of native sheep and valuable genomic information for identifying genes associated with important reproductive traits in sheep.

Integrative Metabolic and Transcriptomic Profiling in Camellia oleifera and Camellia meiocarpa Uncover Potential Mechanisms That Govern Triacylglycerol Degradation during Seed Desiccation

Camellia seed oil is a top-end quality of cooking oil in China. The oil quality and quantity are formed during seed maturation and desiccation. So far, it remains largely unresolved whether lipid degradation occurs and contributes to Camellia oil traits. In this study, three different Camellia germplasms, C. oleifera cv. Min 43 (M43), C. meiocarpa var. Qingguo (QG), and C. meiocarpa cv Hongguo (HG) were selected, their seed oil contents and compositions were quantified across different stages of seed desiccation. We found that at the late stage of desiccation, M43 and QG lost a significant portion of seed oil, while such an event was not observed in HG. To explore the molecular bases for the oil loss In M43, the transcriptomic profiling of M43 and HG was performed at the early and the late seed desiccation, respectively, and differentially expressed genes (DEGs) from the lipid metabolic pathway were identified and analyzed. Our data demonstrated that different Camellia species have diverse mechanisms to regulate seed oil accumulation and degradation, and that triacylglycerol-to-terpenoid conversion could account for the oil loss in M43 during late seed desiccation.

Metabolomic profiles and health-promoting functions of Camellia drupifera mature-seeds were revealed relate to their geographical origins using comparative metabolomic analysis and network pharmacology approach

To insight into the chemical components and their health-promoting function of Camellia drupifera mature-seeds (CMS) in Hainan and Liangguang, UPLC-MS/MS- and HS-SPME/GC-MS-based metabolomic analyses and network pharmacology approaches were combined preformed to Camellia drupifera mature-seeds samples (CMSSs). Totally, 1057 metabolites were identified, of which 76 and 99 metabolites were annotated as key active ingredients in Traditional Chinese Medicines and the active pharmaceutical ingredients for seven human disease-resistance, respectively. Comparative analysis revealed different metabolomic profiles of CMSSs from Hainan and Liangguang. KEGG annotation and enrichment analysis showed secondary metabolic pathways, especially "flavone and flavonol biosynthesis", were played important roles. Finally, 22 metabolites that only detected in CMSSs from Hainan or Liangguang were explored as potential indicators to separate CMS from Hainan out of Liangguang. Our findings enhanced the understanding of chemical compositions of CMS and provided valuable information for the healthy development of oil-tea Camellia industry in Hainan.

Genome-wide identification of the WRKY gene family in Camellia oleifera and expression analysis under phosphorus deficiency

Camellia oleifera Abel. is an economically important woody edible-oil species that is mainly cultivated in hilly areas of South China. The phosphorus (P) deficiency in the acidic soils poses severe challenges for the growth and productivity of C. oleifera. WRKY transcription factors (TFs) have been proven to play important roles in biological processes and plant responses to various biotic/abiotic stresses, including P deficiency tolerance. In this study, 89 WRKY proteins with conserved domain were identified from the C. oleifera diploid genome and divided into three groups, with group II further classified into five subgroups based on the phylogenetic relationships. WRKY variants and mutations were detected in the gene structure and conserved motifs of CoWRKYs. Segmental duplication events were considered as the primary driver in the expanding process of WRKY gene family in C. oleifera. Based on transcriptomic analysis of two C. oleifera varieties characterized with different P deficiency tolerances, 32 CoWRKY genes exhibited divergent expression patterns in response to P deficiency stress. qRT-PCR analysis demonstrated that CoWRKY11, -14, -20, -29 and -56 had higher positive impact on P-efficient CL40 variety compared with P-inefficient CL3 variety. Similar expression trends of these CoWRKY genes were further observed under P deficiency with longer treatment period of 120d. The result indicated the expression sensitivity of CoWRKYs on the P-efficient variety and the C. oleifera cultivar specificity on the P deficiency tolerance. Tissue expression difference showed CoWRKYs may play a crucial role in the transportation and recycling P in leaves by affecting diverse metabolic pathways. The available evidences in the study conclusively shed light on the evolution of the CoWRKY genes in C. oleifera genome and provided a valuable resource for further investigation of functional characterization of WRKY genes involved to enhance the P deficiency tolerance in C. oleifera.

Transcriptome Analysis of the Late-Acting Self-Incompatibility Associated with RNase T2 Family in Camellia oleifera

The Camellia oil tree (Camellia oleifera Abel.) is an important nonwood forest species in China, and the majority of its cultivars are late-acting self-incompatibility (LSI) types. Although several studies have examined the mechanism of LSI, the process is quite complicated and unclear. In this study, pollen tube growth and fruit setting of two Camellia oil tree cultivars Huashuo (HS) and Huajin (HJ) were investigated after non and self-pollination, and transcriptomic analysis of the ovaries was performed 48 h after self-pollination to identify the potential genes implicated in the LSI of Camellia oil trees. The results showed that the fruit set of HS was significantly higher than that of HJ after self-pollination. Transcriptomic analysis revealed that plant hormone signal transduction, the phosphatidylinositol signaling system, ATP-binding cassette (ABC) transporters, reactive oxygen species (ROS) metabolism, and Ca2+ signaling were mainly contributed in the LSI of reaction of Camellia oil tree. Moreover, nine RNase T2 genes were identified from the transcriptome analysis, which also showed that CoRNase7 participated in the self-incompatibility reaction in HS. Based on phylogenetic analysis, CoRNase6 was closely related to S-RNase from coffee, and CoRNase7 and CoRNase8 were closely related to S-RNase from Camellia sinensis. The 9 RNase T2 genes successfully produced proteins in prokaryotes. Subcellular localization indicated that CoRNase1 and CoRNase5 were cytoplasmic proteins, while CoRNase7 was a plasma membrane protein. These results screened the main metabolic pathways closely related to LSI in Camellia oil tree, and SI signal transduction might be regulated by a large molecular regulatory network. The discovery of T2 RNases provided evidence that Camellia oil tree might be under RNase-based gametophytic self-incompatibility.

Resequencing of Rosa rugosa accessions revealed the history of population dynamics, breed origin, and domestication pathways

BACKGROUND

Rosa rugosa is a shrub that originated in China and has economic and ecological value. However, during the development of R. rugosa, the genetic background was chaotic, and the genetic structure among different wild populations was unclear, as well as wild and cultivated accessions. Here, we report whole-genome resequencing of wild and cultivated R. rugosa accessions.

RESULTS

A total of 19,041,284 SNPs were identified in 188 R. rugosa accessions and 3 R. chinensis accessions by resequencing. Population genetic analysis revealed that cultivated and wild groups were separated very early. All R. rugosa accessions were divided into 8 categories based on genetic structure: (1) Weihai, Yantai, and Liaoning category, (2) Jilin category, and (3) Hammonasset category (above three are wild); (4) traditional varieties, (5) hybrids between R. rugosa and R. chinensis, (6) Zizhi Rose, (7) Kushui Rose, (8) hybrids between R. rugosa and R. multiflora. We found that the heterozygosity and genetic diversity of wild accessions were generally lower than those of cultivated individuals. The genes that were selected during cultivation were identified, and it was found that these genes were mainly related to environmental adaptation and growth.

CONCLUSIONS

The Jilin population was the oldest population and later migrated to Liaoning and then migrated to Yantai and Weihai by sea regression in the Bohai Basin. The Hammonasset naturalized population probably originated from the Jilin population and then experienced separate differentiation. The long-term asexual reproduction pattern of R. rugosa decreased genetic diversity in the wild population. During R. rugosa cultivation, the ancestors of the Jilin population were involved in breeding traditional varieties, after which almost no wild individuals were engaged in breeding. However, in recent decades, cross breeding of R. rugosa started the utilization of wild germplasms. In comparison, some other species play important roles in variety formation. Few genes related to economic traits were selected, suggesting no directional domestication in the R. rugosa cultivation process.

Transcriptome and Anatomical Comparisons Reveal the Effects of Methyl Jasmonate on the Seed Development of Camellia oleifera

Seed is a major storage organ that determines the yield and quality of Camellia oleifera (C. oleifera). Methyl jasmonate (MeJA) is a signaling molecule involved in plant growth and development. However, the role of MeJA in the development of C. oleifera seeds remains a mystery. This study demonstrated that the larger seeds induced by MeJA resulted from more cell numbers and a larger cell area in the outer seed coat and embryo at the cellular level. At the molecular level, MeJA could regulate the expression of factors in the known signaling pathways of seed size control as well as cell proliferation and expansion, resulting in larger seeds. Furthermore, the accumulation of oil and unsaturated fatty acids due to MeJA-inducement was attributed to the increased expression of fatty acid biosynthesis-related genes but reduced expression of fatty acid degradation-related genes. CoMYC2, a key regulator in jasmonate signaling, was considered a potential hub regulator which directly interacted with three hub genes (CoCDKB2-3, CoCYCB2-3, and CoXTH9) related to the seed size and two hub genes (CoACC1 and CoFAD2-3) related to oil accumulation and fatty acid biosynthesis by binding to their promoters. These findings provide an excellent target for the improvement of the yield and quality in C. oleifera.

The lineage-specific evolution of the oleosin family in Theaceae

Oleosins play essential roles in stabilization of lipid droplets (LDs) and seed oil production. However, evolution of this gene family has not been reported in Theaceae, a large plant family that contains many important tea and oil tea species. In this study, a total of 65 oleosin genes were identified in nine genome-sequenced Theaceae species. Among these genomes, the gene number of oleosin showed significant difference, with Camellia sinensis var. sinensis cv. Shuchazao and Camellia lanceoleosa displayed more oleosin numbers than other species. Phylogenetic analyses revealed that Theaceae oleosin genes were classified into three clades (U, SL, SH) respectively. Proteins within the same clade had similar gene structure and motif composition. Segmental duplication was the primary driving force for the evolution of oleosin genes in Shuchazao (SCZ), Huangdan (HD), C.lanceoleosa (Cla), and wild tea (DASZ). Synteny analysis showed that most oleosin genes displayed inter-species synteny among tea and oil tea species. Expression analysis demonstrated that oleosin genes were specifically expressed in seed and kernel of Huangdan (HD) and C.lanceoleosa. Moreover, expression divergence was observed in paralogous pairs and ∼1-2 oleosin genes in each clade have become activate. This study leads to a comprehensive understanding of evolution of oleosin family in Theaceae, and provides a rich resource to further address the functions of oleosin in tea and oil tea species.

Co-regulatory effects of hormone and mRNA-miRNA module on flower bud formation of Camellia oleifera

Few flower buds in a high-yield year are the main factors restricting the yield of Camellia oleifera in the next year. However, there are no relevant reports on the regulation mechanism of flower bud formation. In this study, hormones, mRNAs, and miRNAs were tested during flower bud formation in MY3 ("Min Yu 3," with stable yield in different years) and QY2 ("Qian Yu 2," with less flower bud formation in a high-yield year) cultivars. The results showed that except for IAA, the hormone contents of GA₃, ABA, tZ, JA, and SA in the buds were higher than those in the fruit, and the contents of all hormones in the buds were higher than those in the adjacent tissues. This excluded the effect of hormones produced from the fruit on flower bud formation. The difference in hormones showed that 21-30 April was the critical period for flower bud formation in C. oleifera; the JA content in MY3 was higher than that in QY2, but a lower concentration of GA₃ contributed to the formation of the C. oleifera flower bud. JA and GA₃ might have different effects on flower bud formation. Comprehensive analysis of the RNA-seq data showed that differentially expressed genes were notably enriched in hormone signal transduction and the circadian system. Flower bud formation in MY3 was induced through the plant hormone receptor TIR1 (transport inhibitor response 1) of the IAA signaling pathway, the miR535-GID1c module of the GA signaling pathway, and the miR395-JAZ module of the JA signaling pathway. In addition, the expression of core clock components GI (GIGANTEA) and CO (CONSTANS) in MY3 increased 2.3-fold and 1.8-fold over that in QY2, respectively, indicating that the circadian system also played a role in promoting flower bud formation in MY3. Finally, the hormone signaling pathway and circadian system transmitted flowering signals to the floral meristem characteristic genes LFY (LEAFY) and AP1 (APETALA 1) via FT (FLOWERING LOCUS T) and SOC1 (SUPPRESSOR OF OVEREXPRESSION OF CO 1) to regulate flower bud formation. These data will provide the basis for understanding the mechanism of flower bud alternate formation and formulating high yield regulation measures for C. oleifera.

Identification of Camellia oleifera WRKY transcription factor genes and functional characterization of CoWRKY78

Camellia oleifera Abel is a highly valued woody edible oil tree, which is endemic to China. It has great economic value because C. oleifera seed oil contains a high proportion of polyunsaturated fatty acids. C. oleifera anthracnose caused by Colletotrichum fructicola, poses a serious threat to C. oleifera growth and yield and causes the benefit of the C. oleifera industry to suffer directly. The WRKY transcription factor family members have been widely characterized as vital regulators in plant response to pathogen infection. Until now, the number, type and biological function of C. oleifera WRKY genes are remains unknown. Here, we identified 90 C. oleifera WRKY members, which were distributed across 15 chromosomes. C. oleifera WRKY gene expansion was mainly attributed to segmental duplication. We performed transcriptomic analyses to verify the expression patterns of CoWRKYs between anthracnose-resistant and -susceptible cultivars of C. oleifera. These results demonstrated that multiple candidate CoWRKYs can be induced by anthracnose and provide useful clues for their functional studies. CoWRKY78, an anthracnose-induced WRKY gene, was isolated from C. oleifera. It was significantly down-regulated in anthracnose-resistant cultivars. Overexpression of CoWRKY78 in tobacco markedly reduced resistance to anthracnose than WT plants, as evidenced by more cell death, higher malonaldehyde content and reactive oxygen species (ROS), but lower activities of superoxide dismutase (SOD), peroxidase (POD), as well as phenylalanine ammonia-lyase (PAL). Furthermore, the expression of multiple stress-related genes, which are associated with ROS-homeostasis (NtSOD and NtPOD), pathogen challenge (NtPAL), and pathogen defense (NtPR1, NtNPR1, and NtPDF1.2) were altered in the CoWRKY78-overexpressing plants. These findings increase our understanding of the CoWRKY genes and lay the foundation for the exploration of anthracnose resistance mechanisms and expedite the breeding of anthracnose-resistant C. oleifera cultivars.

Understanding the Origin and Evolution of Tea (Camellia sinensis [L.]): Genomic Advances in Tea

Tea, which is processed by the tender shoots or leaves of tea plant (Camellia sinensis), is one of the most popular nonalcoholic beverages in the world and has numerous health benefits for humans. Along with new progress in biotechnologies, the refined chromosome-scale reference tea genomes have been achieved, which facilitates great promise for the understanding of fundamental genomic architecture and evolution of the tea plants. Here, we summarize recent achievements in genome sequencing in tea plants and review the new progress in origin and evolution of tea plants by population sequencing analysis. Understanding the genomic characterization of tea plants is import to improve tea quality and accelerate breeding in tea plants.

Genetic differentiation and genetic structure of mixed-ploidy Camellia hainanica populations

Camellia hainanica, which is common in China's Hainan Province, is an important woody olive tree species. Due to many years of geographic isolation, C. hainanica has not received the attention it deserves, which limits the exploitation of germplasm resources. Therefore, it is necessary to study population genetic characteristics for further utilization and conservation of C. hainanica. In this study, 96 individuals in six wild Camellia hainanica populations were used for ploidy analysis of the chromosome number, and the genetic diversity and population structure were investigated using 12 pairs of SSR primers. The results show complex ploidy differentiation in C. hainanica species. The ploidy of wild C. hainanica includes tetraploid, pentaploid, hexaploid, heptaploid, octoploid and decaploid species. Genetic analysis shows that genetic diversity and genetic differentiation among populations are low. Populations can be divided into two clusters based on their genetic structure, which matches their geographic location. Finally, to further maintain the genetic diversity of C. hainanica, ex-situ cultivation and in-situ management measures should be considered to protect it in the future.

Novel insight into anthocyanin metabolism and molecular characterization of its key regulators in Camellia sasanqua

Flower color is a trait that affects the ornamental value of a plant. Camellia sasanqua is a horticultural plant with rich flower color, but little is known about the regulatory mechanism of color diversity in this plant. Here, the anthocyanin profile of 20 C. sasanqua cultivars revealed and quantified 11 anthocyanin derivatives (five delphinidin-based and six cyanidin-based anthocyanins) for the first time. Cyanidin-3-O-(6-O-(E)-p-coumaroyl)-glucoside was the main contributor to flower base color, and the accumulation of cyanidin and delphinidin derivatives differed in the petals. To further explore the molecular mechanism of color divergence, a transcriptome analysis was performed using C. sasanqua cultivars 'YingYueYe', 'WanXia', 'XueYueHua', and'XiaoMeiGui'. The co-expression network related to differences in delphinidin and cyanidin derivatives accumulation was identified. Eleven candidate genes encoding key enzymes (e.g., F3H, F3'H, and ANS) were involved in anthocyanin biosynthesis. Moreover, 27 transcription factors were screened as regulators of the two types of accumulating anthocyanins. The association was suggested by correlation analysis between the expression levels of the candidate genes and the different camellia cultivars. We concluded that cyanidin and delphinidin derivatives are the major drivers of color diversity in C. sasanqua. This finding provides valuable resources for the study of flower color in C. sasanqua and lays a foundation for genetic modification of anthocyanin biosynthesis.

Gene coexpression analysis reveals key pathways and hub genes related to late-acting self-incompatibility in Camellia oleifera

INTRODUCTION

Self-incompatibility (SI) is an important strategy for plants to maintain abundant variation to enhance their adaptability to the environment. Camellia oleifera is one of the most important woody oil plants and is widely cultivated in China. Late acting self-incompatibility (LSI) in C. oleifera results in a relatively poor fruit yield in the natural state, and understanding of the LSI mechanism remains limited.

METHODS

To better understand the molecular expression and gene coexpression network in the LSI reaction in C. oleifera, we conducted self- and cross-pollination experiments at two different flower bud developmental stages (3-4 d before flowering and 1 d before flowering), and cytological observation, fruit setting rate (FSR) investigation and RNA-Seq analysis were performed to investigate the mechanism of the male -female interaction and identify hub genes responsible for the LSI in C. oleifera.

RESULTS

Based on the 21 ovary transcriptomes, a total of 7669 DEGs were identified after filtering out low-expression genes. Weighted gene coexpression network analysis (WGCNA) divided the DEGs into 15 modules. Genes in the blue module (1163 genes) were positively correlated with FSR, and genes in the pink module (339 genes) were negatively correlated with FSR. KEGG analysis indicated that flavonoid biosynthesis, plant MAPK signaling pathways, ubiquitin-mediated proteolysis, and plant-pathogen interaction were the crucial pathways for the LSI reaction. Fifty four transcription factors (TFs) were obtained in the two key modules, and WRKY and MYB were potentially involved in the LSI reaction in C. oleifera. Network establishment indicated that genes encoding G-type lectin S-receptor-like serine (lecRLK), isoflavone 3'-hydroxylase-like (CYP81Q32), cytochrome P450 87A3-like (CYP87A3), and probable calcium-binding protein (CML41) were the hub genes that positively responded to the LSI reaction. The other DEGs inside the two modules, including protein RALF-like 10 (RALF), F-box and pectin acetylesterase (MTERF5), might also play vital roles in the LSI reaction in C. oleifera.

DISCUSSION

Overall, our study provides a meaningful resource for gene network studies of the LSI reaction process and subsequent analyses of pollen-pistil interactions and TF roles in the LSI reaction, and it also provides new insights for exploring the mechanisms of the LSI response.

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.

BSR and Full-Length Transcriptome Approaches Identified Candidate Genes for High Seed Ratio in Camellia vietnamensis

(1) Background: C. vietnamensis is very suitable for growth in the low hilly areas of southern subtropical regions. Under appropriate conditions, the oil yield of C. vietnamensis can reach 1125 kg/ha (the existing varieties can reach 750 kg/ha). Moreover, the fruit of C. vietnamensis is large and the pericarp is thick (>5 cm). Therefore, a high seed ratio has become the main target economic trait for the breeding of C. vietnamensis. (2) Methods: A half-sibling population of C. vietnamensis plants with a combination of high and low seed ratios was constructed by crossing a C. vietnamensis female parent. Bulked segregant RNA analysis and full-length transcriptome sequencing were performed to determine the molecular mechanisms underlying a high seed ratio. (3) Results: Seed ratio is a complex quantitative trait with a normal distribution, which is significantly associated with four other traits of fruit (seed weight, seed number, fruit diameter, and pericarp thickness). Two candidate regions related to high seed ratio (HSR) were predicted. One spanned 140.8−148.4 Mb of chromosome 2 and was associated with 97 seed-yield-related candidate genes ranging in length from 278 to 16,628 bp. The other spanned 35.3−37.3 Mb on chromosome 15 and was associated with 38 genes ranging in length from 221 to 16,928 bp. Using the full-length transcript as a template, a total of 115 candidate transcripts were obtained, and 78 transcripts were predicted to be functionally annotated. The DEGs from two set pairs of cDNA sequencing bulks were enriched to cytochrome p450 CYP76F14 (KOG0156; GO:0055114, HSR4, HSR7), the gibberellin phytohormone pathway (GO:0016787, HSR5), the calcium signaling pathway (GO:0005509, HSR6), the polyubiquitin-PPAR signaling pathway (GO:0005515, HSR2, HSR3), and several main transcription factors (bZIP transcription factor, HSR1) in C. vietnamensis.

Genetic Diversity Analysis and Core Germplasm Collection Construction of Camellia oleifera Based on Fruit Phenotype and SSR Data

Many Camellia oleifera germplasm resources were collected from Guizhou Province, but the fruit morphological variation and genetic diversity of C. oleifera germplasm resources remain unclear. The genetic diversity of C. oleifera germplasms resources in Guizhou was studied based on fruit traits and simple sequence repeat (SSR) molecular markers to build a core collection. This paper aims to provide a scientific basis for the collection, management, development, and utilization of C. oleifera resources in Guizhou province. The variation coefficients among and within varieties of seven fruit phenotypic traits of C. oleifera ranged from 11.79% to 61.76% and from 8.15% to 42.31%, respectively, showing rich phenotypic variation. Furthermore, 12 SSR markers were used to analyze the genetic diversity. These primers generated 214 polymorphic bands, and the average number was 17.833. The average number of effective alleles (Ne), Shannon's information index (I), observed heterozygosity (Ho), expected heterozygosity (He), polymorphic information content (PIC), and major allele frequency (MAF) were 8.999, 2.301, 0.965, 0.50, 0.836, and 0.238, respectively. The results showed that 12 SSR markers had high polymorphism, and the genetic diversity of 167 C. oleifera germplasm resources was high. Based on SSR molecular marker information and fruit traits clustering, 167 C. oleifera germplasm resources were divided into three groups. When constructing core collections based on fruit traits and molecular marker information, the PowerCore-25 of core collections greatly preserves fruit traits and improves genetic diversity. This paper can provide a reference for the genetic diversity and fruit traits variation of C. camellia germplasm resources in Guizhou Province. It is significant for establishing a core collection, thus promoting germplasm innovation and the development of the oil tea industry in Guizhou.

Overexpression of dihydroflavonol 4-reductase (CoDFR) boosts flavonoid production involved in the anthracnose resistance

The outbreak of anthracnose caused by Colletotrichum spp. represents a devastating epidemic that severely affects oil tea (Camellia oleifera) production in China. However, the unknown resistance mechanism to anthracnose in C. oleifera has impeded the progress of breeding disease-resistant varieties. In this study, we investigated the physiological responses of resistant and susceptible lines during C. gloeosporioides infection. Our results showed that the accumulation of malondialdehyde (MDA), catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD) in both disease-resistant and susceptible lines increased by C. gloeosporioides infection. Also, disease-resistant lines exhibited lower MDA, but higher POD, SOD, and CAT activities compared to susceptible lines. The accumulation of flavonoids in both resistant and susceptible C. oleifera leaves increased following C. gloeosporioides infection, and the increase was greater in resistant lines. Further, we identified and functionally characterized the dihydroflavonol 4-reductase (CoDFR) from the resistant C. oleifera line. We showed that the full-length coding sequence (CDS) of CoDFR is 1044 bp encoding 347 amino acids. The overexpression of CoDFR in tobacco altered the expression of flavonoid biosynthetic genes, resulting in an increased flavonoid content in leaves. CoDFR transgenic tobacco plants exhibited increased anthracnose resistance. Furthermore, the transgenic plants had higher salicylic acid content. These findings offer potential insights into the pivotal role of CoDFR involved in flavonoid-mediated defense mechanisms during anthracnose invasion in resistant C. oleifera.

Genome Survey and SSR Analysis of Camellia nitidissima Chi (Theaceae)

Camellia nitidissima Chi (CNC), a species of golden Camellia, is well known as "the queen of camellias." It is an ornamental, medicinal, and edible plant grown in China. In this study, we conducted a genome survey sequencing analysis and simple sequence repeat (SSR) identification of CNC using the Illumina sequencing platform. The 21-mer analysis predicted its genome size to be 2,778.82 Mb, with heterozygosity and repetition rates of 1.42% and 65.27%, respectively. The CNC genome sequences were assembled into 9,399,197 scaffolds, covering ∼2,910 Mb and an N50 of 869 base pair. Its genomic characteristics were found to be similar to those of Camellia oleifera. In addition, 1,940,616 SSRs were identified from the genome data, including mono-(61.85%), di-(28.71%), tri-(6.51%), tetra-(1.85%), penta-(0.57%), and hexanucleotide motifs (0.51%). We believe these data will provide a useful foundation for the development of novel molecular markers for CNC as well as for further whole-genome sequencing of CNC.

Transcriptome Analysis Reveals Putative Induction of Floral Initiation by Old Leaves in Tea-Oil Tree (Camellia oleifera ‘changlin53’)

Floral initiation is a major phase change in the spermatophyte, where developmental programs switch from vegetative growth to reproductive growth. It is a key phase of flowering in tea-oil trees that can affect flowering time and yield, but very little is known about the molecular mechanism of floral initiation in tea-oil trees. A 12-year-old Camellia oleifera (cultivar ‘changlin53’) was the source of experimental materials in the current study. Scanning electron microscopy was used to identify the key stage of floral initiation, and transcriptome analysis was used to reveal the transcriptional regulatory network in old leaves involved in floral initiation. We mined 5 DEGs related to energy and 55 DEGs related to plant hormone signal transduction, and we found floral initiation induction required a high level of energy metabolism, and the phytohormones signals in the old leaves regulate floral initiation, which occurred at stage I and II. Twenty-seven rhythm-related DEGs and 107 genes associated with flowering were also identified, and the circadian rhythm interacted with photoperiod pathways to induce floral initiation. Unigene0017292 (PSEUDO-RESPONSE REGULATOR), Unigene0046809 (LATE ELONGATED HYPOCOTYL), Unigene0009932 (GIGANTEA), Unigene0001842 (CONSTANS), and Unigene0084708 (FLOWER LOCUS T) were the key genes in the circadian rhythm-photoperiod regulatory network. In conjunction with morphological observations and transcriptomic analysis, we concluded that the induction of floral initiation by old leaves in C. oleifera ‘changlin53’ mainly occurred during stages I and II, floral initiation was completed during stage III, and rhythm–photoperiod interactions may be the source of the main signals in floral initiation induced by old leaves.

A High-Quality Genome Assembly of the Mitochondrial Genome of the Oil-Tea Tree Camellia gigantocarpa (Theaceae)

Camellia gigantocarpa is one of the oil-tea trees whose seeds can be used to extract high-quality vegetable oil. To date, there are no data on the mitochondrial genome of the oil-tea tree, in contrast to the tea-tree C. sinensis, which belongs to the same genus. In this paper, we present the first complete mitochondrial genomes of C. gigantocarpa obtained using PacBio Hi-Fi (high-fidelity) and Hi-C sequencing technologies to anchor the 970,410 bp genome assembly into a single sequence. A set of 44 protein-coding genes, 22 non-coding genes, 746 simple sequence repeats (SSRs), and more than 201 kb of repetitive sequences were annotated in the genome assembly. The high percentage of repetitive sequences in the mitochondrial genome of C. gigantocarpa (20.81%) and C.sinensis (22.15%, tea tree) compared to Arabidopsis thaliana (4.96%) significantly increased the mitogenome size in the genus Camellia. The comparison of the mitochondrial genomes between C. gigantocarpa and C. sinensis revealed genes exhibit high variance in gene order and low substitution rate within the genus Camellia. Information on the mitochondrial genome provides a better understanding of the structure and evolution of the genome in Camellia and may contribute to further study of the after-ripening process of oil-tea trees.

Floral organ transcriptome in Camellia sasanqua provided insight into stamen petaloid

BACKGROUND

The cultivated Camellia sasanqua forms a divergent double flower pattern, and the stamen petaloid is a vital factor in the phenomenon. However, the regulation mechanism remains largely unclear.

RESULTS

Here, a comprehensive comparative transcriptome analysis of the wild-type, "semi-double", "peony double", and "rose double" was performed. The cluster analysis of global gene expression level showed petal and stamen difficulty separable in double flower. The crucial pathway and genes related to double flower patterns regulation were identified by pairwise comparisons and weighted gene coexpression network (WGCNA). Divergent genes expression, such as AUX1 and AHP, are involved in plant hormone signaling and photosynthesis, and secondary metabolites play an important role. Notably, the diversity of a petal-specific model exhibits a similar molecular signature to the stamen, containing extensin protein and PSBO1, supporting the stamen petaloid point. Moreover, the expansion of class A gene activity influenced the double flower formation, showing that the key function of gene expression was probably demolished.

CONCLUSIONS

Overall, this work confirmed the ABCE model and provided new insights for elucidating the molecular signature of double formation.

Chromosome-scale genome assembly of Camellia sinensis combined with multi-omics provides insights into its responses to infestation with green leafhoppers

The tea plant (Camellia sinensis) is an important economic crop, which is becoming increasingly popular worldwide, and is now planted in more than 50 countries. Tea green leafhopper is one of the major pests in tea plantations, which can significantly reduce the yield and quality of tea during the growth of plant. In this study, we report a genome assembly for DuyunMaojian tea plants using a combination of Oxford Nanopore Technology PromethION^™ with high-throughput chromosome conformation capture technology and used multi-omics to study how the tea plant responds to infestation with tea green leafhoppers. The final genome was 3.08 Gb. A total of 2.97 Gb of the genome was mapped to 15 pseudo-chromosomes, and 2.79 Gb of them could confirm the order and direction. The contig N50, scaffold N50 and GC content were 723.7 kb, 207.72 Mb and 38.54%, respectively. There were 2.67 Gb (86.77%) repetitive sequences, 34,896 protein-coding genes, 104 miRNAs, 261 rRNA, 669 tRNA, and 6,502 pseudogenes. A comparative genomics analysis showed that DuyunMaojian was the most closely related to Shuchazao and Yunkang 10, followed by DASZ and tea-oil tree. The multi-omics results indicated that phenylpropanoid biosynthesis, α-linolenic acid metabolism, flavonoid biosynthesis and 50 differentially expressed genes, particularly peroxidase, played important roles in response to infestation with tea green leafhoppers (Empoasca vitis Göthe). This study on the tea tree is highly significant for its role in illustrating the evolution of its genome and discovering how the tea plant responds to infestation with tea green leafhoppers will contribute to a theoretical foundation to breed tea plants resistant to insects that will ultimately result in an increase in the yield and quality of tea.

Analysis of Camellia oleifera transcriptome reveals key pathways and hub genes involved during different photoperiods

BACKGROUND

Camellia oleifera Abel. (C. oleifera) is an important traditional woody species in China that produces edible oil. However, the current literature lacks a proper understanding of C. oleifera's ability to adapt to different photoperiods.

RESULTS

Our results indicate that the photoperiod can significantly impact flowering time in C. oleifera. We grew a total of nine samples under the short day condition (SD), middle day condition (MD) and long day condition (LD). Transcriptome analysis yielded 66.94 Gb of high-quality clean reads, with an average of over 6.73 Gb of reads for per sample. Following assembly, a total of 120,080 transcripts were obtained and 94,979 unigenes annotated. A total of 3475 differentially expressed genes (DEGs) were identified between the SD_MD, SD_LD, and MD_LD gene sets. Moreover, WGCNA identified ten gene modules. Genes in pink module (92 genes) were positively correlated with SD, and negatively correlated with both MD and LD. Genes in the magenta module (42 genes) were positively correlated with MD and negatively correlated with both LD and SD. Finally, genes in the yellow module (1758 genes) were positively correlated with both SD and MD, but negatively correlated with LD. KEGG enrichment analysis revealed that genes in the pink, magenta, and yellow modules were involved in flavonoid biosynthesis, amino sugar and nucleotide sugar metabolism and circadian rhythm pathways. Additionally, eight hub genes (GI, AP2, WRKY65, SCR, SHR, PHR1, ERF106, and SCL3) were obtained through network analysis. The hub genes had high connectivity with other photoperiod-sensitive DEGs. The expression levels of hub genes were verified by qRT-PCR analysis.

CONCLUSION

An increase in light duration promotes earlier flowering of C. oleifera. Flavonoid biosynthesis, amino sugar and nucleotide sugar metabolism, and circadian rhythm pathways may function in the photoperiodic flowering pathway of C. oleifera. We also identified eight hub genes that may play a role in this pathway. Ultimately, this work contributes to our understanding of the photoperiodic flowering pathway of C. oleifera and further informs molecular breeding programs on the plant's photoperiodic sensitivity.

Identification of long non-coding RNAs and microRNAs involved in anther development in the tropical Camellia oleifera

BACKGROUND

Explored the molecular science of anther development is important for improving productivity and overall yield of crops. Although the role of regulatory RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), in regulating anther development has been established, their identities and functions in Camellia oleifera, an important industrial crop, have yet not been clearly explored. Here, we report the identification and characterization of genes, lncRNAs and miRNAs during three stages of the tropical C. oleifera anther development by single-molecule real-time sequencing, RNA sequencing and small RNA sequencing, respectively.

RESULTS

These stages, viz. the pollen mother cells stage, tetrad stage and uninucleate pollen stage, were identified by analyzing paraffin sections of floral buds during rapid expansion periods. A total of 18,393 transcripts, 414 putative lncRNAs and 372 miRNAs were identified, of which 5,324 genes, 115 lncRNAs, and 44 miRNAs were differentially accumulated across three developmental stages. Of these, 44 and 92 genes were predicted be regulated by 37 and 30 differentially accumulated lncRNAs and miRNAs, respectively. Additionally, 42 differentially accumulated lncRNAs were predicted as targets of 27 miRNAs. Gene ontology enrichment indicated that potential target genes of lncRNAs were enriched in photosystem II, regulation of autophagy and carbohydrate phosphatase activity, which are essential for anther development. Functional annotation of genes targeted by miRNAs indicated that they are relevant to transcription and metabolic processes that play important roles in microspore development. An interaction network was built with 2 lncRNAs, 6 miRNAs and 10 mRNAs. Among these, miR396 and miR156 family were up-regulated, while their targets, genes (GROWTH REGULATING FACTORS and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes) and lncRNAs, were down-regulated. Further, the trans-regulated targets of these lncRNAs, like wall-associated kinase2 and phosphomannose isomerase1, are involved in pollen wall formation during anther development.

CONCLUSIONS

This study unravels lncRNAs, miRNAs and miRNA-lncRNA-mRNA networks involved in development of anthers of the tropical C. oleifera lays a theoretical foundation for further elucidation of regulatory roles of lncRNAs and miRNAs in anther development.

Phylogenomics Resolves the Phylogeny of Theaceae by Using Low-Copy and Multi-Copy Nuclear Gene Makers and Uncovers a Fast Radiation Event Contributing to Tea Plants Diversity

Simple Summary

The Theaceae includes more than 300 species of great morphological diversity and has immense economic, cultural, and ornamental values. However, the evolutionary history of this family remains elusive. We integrated 91 genomes and transcriptome datasets of Theaceae and successfully resolved the phylogeny of Theaceae including relatives of cultivated tea plants from both extensive low-copy and multi-copy nuclear gene markers. Bayes-based molecular dating revealed that the ancestor of the tea family originated slightly earlier than the K-Pg boundary (Mass extinction events including the extinction of dinosaurs) with early diversification of three tribes associated with the Early Eocene Climatic Optimum. Further speciation analysis suggested a sole significant diversification shift rate in the common ancestor of Camellia associated with the Mid-Miocene Climatic Optimum. Collectively, polyploidy events, and key morphological innovation characters, such as pericarp with seed coat hardening, could possibly contribute to the Theaceae species diversity.

Abstract

Tea is one of the three most popular nonalcoholic beverages globally and has extremely high economic and cultural value. Currently, the classification, taxonomy, and evolutionary history of the tea family are largely elusive, including phylogeny, divergence, speciation, and diversity. For understanding the evolutionary history and dynamics of species diversity in Theaceae, a robust phylogenetic framework based on 1785 low-copy and 79,103 multi-copy nuclear genes from 91 tea plant genomes and transcriptome datasets had been reconstructed. Our results maximumly supported that the tribes Stewartieae and Gordonieae are successive sister groups to the tribe Theeae from both coalescent and super matrix ML tree analyses. Moreover, in the most evolved tribe, Theeae, the monophyletic genera Pyrenaria, Apterosperma, and Polyspora are the successive sister groups of Camellia. We also yield a well-resolved relationship of Camellia, which contains the vast majority of Theaceae species richness. Molecular dating suggests that Theaceae originated in the late L-Cretaceous, with subsequent early radiation under the Early Eocene Climatic Optimal (EECO) for the three tribes. A diversification rate shift was detected in the common ancestors of Camellia with subsequent acceleration in speciation rate under the climate optimum in the early Miocene. These results provide a phylogenetic framework and new insights into factors that likely have contributed to the survival of Theaceae, especially a successful radiation event of genus Camellia members to subtropic/tropic regions. These novel findings will facilitate the efficient conservation and utilization of germplasm resources for breeding cultivated tea and oil-tea. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Theaceae.

Applications of Chinese Camellia oleifera and its By-Products: A Review

Camellia oleifera is a woody oil tree species unique to China that has been cultivated and used in China for more than 2,300 years. Most biological research on C. oleifera in recent years has focused on the development of new varieties and breeding. Novel genomic information has been generated for C. oleifera, including a high-quality reference genome at the chromosome level. Camellia seeds are used to process high-quality edible oil; they are also often used in medicine, health foods, and daily chemical products and have shown promise for the treatment and prevention of diseases. C. oleifera by-products, such as camellia seed cake, saponin, and fruit shell are widely used in the daily chemical, dyeing, papermaking, chemical fibre, textile, and pesticide industries. C. oleifera shell can also be used to prepare activated carbon electrodes, which have high electrochemical performance when used as the negative electrode of lithium-ion batteries. C. oleifera is an economically valuable plant with diverse uses, and accelerating the utilization of its by-products will greatly enhance its industrial value.

Critical metabolic pathways and SAD/FADs, WRI1s, and DGATs cooperate for high-oleic acid oil production in developing oil tea (Camellia oleifera) seeds

Oil tea trees produce high-quality edible oils with desirably high oleic acid (18:1) and low linoleic (18:2) and linolenic (18:3) fatty acid (FA) levels, but limited understanding of tea oil biosynthesis and regulation has become a significant obstacle for the breeding of high-yield and -quality oil tea varieties. By integrating metabolite and transcriptome analyses of developing oil tea seeds, we dissected the critical metabolic pathways, including glycolysis, fatty acid, and triacylglycerol (TAG) biosynthesis, as well as genes essential for tea seed oil production. Two plastidic stearoyl-acyl carrier protein desaturases (CoSAD1 and 2) and two endoplasmic reticulum-localized FA desaturases (CoFAD2 and 3) were functionally characterized as responsible for high 18:1 and low 18:2 and 18:3 proportions in tea oils. Two diacylglycerol O-acyltransferases (CoDGAT1 and 2) that may prefer to synthesize 18:1-TAG were functionally characterized and might be also important for high 18:1-TAG production. The highly expressed CoWRI1a and b were identified and characterized as activators of glycolysis and regulators of directing source carbon flux into FA biosynthesis in developing oil tea seeds. The upregulated CoSADs with downregulated CoFAD2 and CoFAD3 at the late seed developmental stages mainly accounted for high 18:1 levels. Two CoDGATs might be responsible for assembling TAGs with oleoyl acyl chains, whilst two CoWRI1s regulated carbons from parental sources, partitioning into oil production in oil tea embryo sinks. This study provides a deep understanding of the biosynthesis of tea seed oils and information on genes that may be used as molecular markers to breed oil tea varieties with higher oil yield and quality.

Genetic relationships and low diversity among the tea-oil Camellia species in Sect. Oleifera, a bulk woody oil crop in China

Tea-oil Camellia is one of the four woody oil crops in the world and has high ecological, economic and medicinal values. However, there are great differences in the classification and merging of tea-oil Camellia Sect. Oleifera species, which brings difficulties to the innovative utilization and production of tea-oil Camellia resources. Here, ISSR, SRAP and chloroplast sequence markers were analyzed in 18 populations of tea-oil Camellia Sect. Oleifera species to explore their phylogenetic relationships and genetic diversity. The results showed that their genetic diversity were low, with mean H and π values of 0.16 and 0.00140, respectively. There was high among-population genetic differentiation, with ISSR and SRAP markers showing an Fst of 0.38 and a high Nm of 1.77 and cpDNA markers showing an Fst of 0.65 and a low Nm of 0.27. The C. gauchowensis, C. vietnamensis and Hainan Island populations formed a single group, showing the closest relationships, and supported being the same species for them with the unifying name C. drupifera and classifying the resources on Hainan Island as C. drupifera. The tea-oil Camellia resources of Hainan Island should be classified as a special ecological type or variety of C. drupifera. However, cpDNA marker-based STRUCTURE analysis showed that the genetic components of the C. osmantha population formed an independent, homozygous cluster; hence, C. osmantha should be a new species in Sect. Oleifera. The C. oleifera var. monosperma and C. oleifera populations clustered into two distinct clades, and the C. oleifera var. monosperma populations formed an independent cluster, accounting for more than 99.00% of its genetic composition; however, the C. oleifera populations contained multiple different cluster components, indicating that C. oleifera var. monosperma significantly differs from C. oleifera and should be considered the independent species C. meiocarpa. Haplotype analysis revealed no rapid expansion in the tested populations, and the haplotypes of C. oleifera, C. meiocarpa and C. osmantha evolved from those of C. drupifera. Our results support the phylogenetic classification of Camellia subgenera, which is highly significant for breeding and production in tea-oil Camellia.