Haplotype-resolved genome assembly for tetraploid Chinese cherry (Prunus pseudocerasus) offers insights into fruit firmness.

Chinese cherry (Prunus pseudocerasus) holds considerable importance as one of the primary stone fruit crops in China. However, artificially improving its traits and genetic analysis are challenging due to lack of high-quality genomic resources, which mainly result from difficulties associated with resolving its tetraploid and highly heterozygous genome. Herein, we assembled a chromosome-level, haplotype-resolved genome of the cultivar 'Zhuji Duanbing', comprising 993.69 Mb assembled into 32 pseudochromosomes using PacBio HiFi, Oxford Nanopore, and Hi-C. Intra-haplotype comparative analyses revealed extensive intra-genomic sequence and expression consistency. Phylogenetic and comparative genomic analyses demonstrated that P. pseudocerasus was a stable autotetraploid species, closely related to wild P. pusilliflora, with the two diverging ~18.34 million years ago. Similar to other Prunus species, P. pseudocerasus underwent a common whole-genome duplication event that occurred ~139.96 million years ago. Because of its low fruit firmness, P. pseudocerasus is unsuitable for long-distance transportation, thereby restricting its rapid development throughout China. At the ripe fruit stage, P. pseudocerasus cv. 'Zhuji Duanbing' was significantly less firm than P. avium cv. 'Heizhenzhu'. The difference in firmness is attributed to the degree of alteration in pectin, cellulose, and hemicellulose contents. In addition, comparative transcriptomic analyses identified GalAK-like and Stv1, two genes involved in pectin biosynthesis, which potentially caused the difference in firmness between 'Zhuji Duanbing' and 'Heizhenzhu'. Transient transformations of PpsGalAK-like and PpsStv1 increase protopectin content and thereby enhance fruit firmness. Our study lays a solid foundation for functional genomic studies and the enhancement of important horticultural traits in Chinese cherries.

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.

Interfamily Grafted Hybrids Vitis vinifera/Schisandra chinensis Resulted in Transcriptomic, Phenotypic, and Metabolic Changes.

Long-distance transfer of genetic material and metabolites between rootstock and scions is well documented in homo-grafted hybrids but has rarely been reported in genetically-distant grafts where the rootstock and scion belong to different families. In this study, we grafted Vitis vinifera scions onto Schisandra chinensis stocks and obtained 20 vegetative hybrids, Vitis vinifera/Schisandra chinensis (Vs). After 25 years of growth, we found that the phenotypes of the leaves, internodes, and fruits of the Vs hybrids above the graft union resembled an intermediate phenotype between V. vinifera and S. chinensis, and the new traits were stable when propagated vegetatively. We further analyzed genetic differences between Vv plants and Vs hybrids using high-throughput sequencing, while metabolomes were analyzed by liquid chromatography-mass spectrometry (LC-MS). We found a total of 2113 differentially expressed genes (DEGs). GO annotation and KEGG pathway enrichment analysis showed that these DEGs enriched mainly in oxidation-reduction and metabolic processes. Seventy-nine differentially expressed miRNAs (DEMs) containing 27 known miRNAs and 52 novel miRNAs were identified. A degradation analysis detected 840 target genes corresponding to 252 miRNAs, of which 12 DEMs and their corresponding target gene expression levels were mostly negatively correlated. Furthermore, 1188 differential metabolic compounds were identified. In particular, in Vs hybrids, the abundance of the metabolites schizandrin and gomisin as the main medicinal ingredients in S. chinensis were down-regulated and up-regulated, respectively. Our data demonstrated the effects of interfamily grafts on the phenotype, transcript profile and metabolites of the scion, and also provided new insight into the genetic, phenotypic, and metabolic plasticity associated with genetically distant grafted hybrids.

LUX ARRHYTHMO links CBF pathway and jasmonic acid metabolism to regulate cold tolerance of tea plants.

Cold stress declines the quality and yield of tea, yet the molecular basis underlying cold tolerance of tea plants (Camellia sinensis) remains largely unknown. Here, we identified a circadian rhythm component LUX ARRHYTHMO (LUX) that potentially regulates cold tolerance of tea plants through a genome-wide association study and transcriptomic analysis. The expression of CsLUX phased with sunrise and sunset and was strongly induced by cold stress. Genetic assays indicated that CsLUX is a positive regulator of freezing tolerance in tea plants. CsLUX was directly activated by CsCBF1 and repressed the expression level of CsLOX2, which regulates the cold tolerance of tea plants through dynamically modulating jasmonic acid content. Furthermore, we showed that the CsLUX-CsJAZ1 complex attenuated the physical interaction of CsJAZ1 with CsICE1, liberating CsICE1 with transcriptional activities to withstand cold stress. Notably, a single-nucleotide variation of C-to-A in the coding region of CsLUX was functionally validated as the potential elite haplotype for cold response, which provided valuable molecular markers for future cold resistance breeding in tea plants.

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.

Genetic insights into the dissolution of dioecy in diploid persimmon Diospyros oleifera Cheng.

BACKGROUND: Dioecy, a sexual system of single-sexual (gynoecious/androecious) individuals, is rare in flowering plants. This rarity may be a result of the frequent transition from dioecy into systems with co-sexual individuals. RESULTS: In this study, co-sexual expression (monoecy and hermaphroditic development), previously thought to be polyploid-specific in Diospyros species, was identified in the diploid D. oleifeara historically. We characterized potential genetic mechanisms that underlie the dissolution of dioecy to monoecy and andro(gyno)monoecy, based on multiscale genome-wide investigations of 150 accessions of Diospyros oleifera. We found all co-sexual plants, including monoecious and andro(gyno)monoecious individuals, possessed the male determinant gene OGI, implying the presence of genetic factors controlling gynoecia development in genetically male D. oleifera. Importantly, discrepancies in the OGI/MeGI module were found in diploid monoecious D. oleifera compared with polyploid monoecious D. kaki, including no Kali insertion on the promoter of OGI, no different abundance of smRNAs targeting MeGI (a counterpart of OGI), and no different expression of MeGI between female and male floral buds. On the contrary, in both single- and co-sexual plants, female function was expressed in the presence of a genome-wide decrease in methylation levels, along with sexually distinct regulatory networks of smRNAs and their targets. Furthermore, a genome-wide association study (GWAS) identified a genomic region and a DUF247 gene cluster strongly associated with the monoecious phenotype and several regions that may contribute to andromonoecy. CONCLUSIONS: Collectively, our findings demonstrate stable breakdown of the dioecious system in D. oleifera, presumably also a result of genomic features of the Y-linked region.

Chromosome-Level Genome Assembly and Population Genomic Analyses Reveal Geographic Variation and Population Genetic Structure of Prunus tenella.

Prunus tenella is a rare and precious relict plant in China. It is an important genetic resource for almond improvement and an indispensable material in ecological protection and landscaping. However, the research into molecular breeding and genetic evolution has been severely restricted due to the lack of genome information. In this investigation, we created a chromosome-level genomic pattern of P. tenella, 231 Mb in length with a contig N50 of 18.1 Mb by Hi-C techniques and high-accuracy PacBio HiFi sequencing. The present assembly predicted 32,088 protein-coding genes, and an examination of the genome assembly indicated that 94.7% among all assembled transcripts were alignable to the genome assembly; most (97.24%) were functionally annotated. By phylogenomic genome comparison, we found that P. tenella is an ancient group that diverged approximately 13.4 million years ago (mya) from 13 additional closely related species and about 6.5 Mya from the cultivated almond. Collinearity analysis revealed that P. tenella is highly syntenic and has high sequence conservation with almond and peach. However, this species also exhibits many presence/absence variants. Moreover, a large inversion at the 7588 kb position of chromosome 5 was observed, which may have a significant association with phenotypic traits. Lastly, population genetic structure analysis in eight different populations indicated a high genetic differentiation among the natural distribution of P. tenella. This high-quality genome assembly provides critical clues and comprehensive information for the systematic evolution, genetic characteristics, and functional gene research of P. tenella. Moreover, it provides a valuable genomic resource for in-depth study in protection, developing, and utilizing P. tenella germplasm resources.

Transcriptome and Metabolome Reveal Distinct Sugar Accumulation Pattern between PCNA and PCA Mature Persimmon Fruit.

Persimmon (Diospyros kaki) fruit have significant variation between pollination-constant non-astringent (PCNA) and pollination-constant astringent (PCA) persimmons. The astringency type affects not only the soluble tannin concentration but also the accumulation of individual sugars. Thus, we comprehensively investigate the gene expression and metabolite profiles of individual sugars to resolve the formation of flavor differences in PCNA and PCA persimmon fruit. The results showed that soluble sugar, starch content, sucrose synthase, and sucrose invertase were significantly different between PCNA and PCA persimmon fruit. The sucrose and starch metabolism pathway was considerably enriched, and six sugar metabolites involving this pathway were significantly differentially accumulated. In addition, the expression patterns of diferentially expressed genes (such as bglX, eglC, Cel, TPS, SUS, and TREH genes) were significantly correlated with the content of deferentially accumulated metabolites (such as starch, sucrose, and trehalose) in the sucrose and starch metabolism pathway. These results indicated that the sucrose and starch metabolism pathway maintained a central position of sugar metabolism between PCNA and PCA persimmon fruit. Our results provide a theoretical basis for exploring functional genes related to sugar metabolism and provide useful resources for future studies on the flavor differences between PCNA and PCA persimmon fruit.

Allele-aware chromosome-level genome assembly of the autohexaploid Diospyros kaki Thunb.

Artificially improving persimmon (Diospyros kaki Thunb.), one of the most important fruit trees, remains challenging owing to the lack of reference genomes. In this study, we generated an allele-aware chromosome-level genome assembly for the autohexaploid persimmon 'Xiaoguotianshi' (Chinese-PCNA type) using PacBio CCS and Hi-C technology. The final assembly contained 4.52 Gb, with a contig N50 value of 5.28 Mb and scaffold N50 value of 44.01 Mb, of which 4.06 Gb (89.87%) of the assembly were anchored onto 90 chromosome-level pseudomolecules comprising 15 homologous groups with 6 allelic chromosomes in each. A total of 153,288 protein-coding genes were predicted, of which 98.60% were functionally annotated. Repetitive sequences accounted for 64.02% of the genome; and 110,480 rRNAs, 12,297 tRNAs, 1,483 miRNAs, and 3,510 snRNA genes were also identified. This genome assembly fills the knowledge gap in the autohexaploid persimmon genome, which is conducive in the study on the regulatory mechanisms underlying the major economically advantageous traits of persimmons and promoting breeding programs.

Integrated physiological, metabolite and proteomic analysis reveal the glyphosate stress response mechanism in tea plant (Camellia sinensis).

Glyphosate residues can tremendously impact the physiological mechanisms of tea plants, thus threatening tea security and human health. Herein, integrated physiological, metabolite, and proteomic analyses were performed to reveal the glyphosate stress response mechanism in tea plant. After exposure to glyphosate (≥1.25 kg ae/ha), the leaf ultrastructure was damaged, and chlorophyll content and relative fluorescence intensity decreased significantly. The characteristic metabolites catechins and theanine decreased significantly, and the 18 volatile compounds content varied significantly under glyphosate treatments. Subsequently, tandem mass tags (TMT)-based quantitative proteomics was employed to identify the differentially expressed proteins (DEPs) and to validate their biological functions at the proteome level. A total of 6287 proteins were identified and 326 DEPs were screened. These DEPs were mainly catalytic, binding, transporter and antioxidant active proteins, involved in photosynthesis and chlorophyll biosynthesis, phenylpropanoid and flavonoid biosynthesis, sugar and energy metabolism, amino acid metabolism, and stress/defense/detoxification pathway, etc. A total of 22 DEPs were validated by parallel reaction monitoring (PRM), demonstrating that the protein abundances were consistent between TMT and PRM data. These findings contribute to our understanding of the damage of glyphosate to tea leaves and molecular mechanism underlying the response of tea plants to glyphosate.

Regulatory mechanism of MeGI on sexuality in Diospyros oleifera.

Dioecy system is an important strategy for maintaining genetic diversity. The transcription factor MeGI, contributes to dioecy by promoting gynoecium development in Diospyros lotus and D. kaki. However, the function of MeGI in D. oleifera has not been identified. In this study, we confirmed that MeGI, cloned from D. oleifera, repressed the androecium development in Arabidopsis thaliana. Subsequently, chromatin immunoprecipitation-sequencing (ChIP-seq), DNA affinity purification-sequencing (DAP-seq), and RNA-seq were used to uncover the gene expression response to MeGI. The results showed that the genes upregulated and downregulated in response to MeGI were mainly enriched in the circadian rhythm-related and flavonoid biosynthetic pathways, respectively. Additionally, the WRKY DNA-binding protein 28 (WRKY28) gene, which was detected by ChIP-seq, DAP-seq, and RNA-seq, was emphasized. WRKY28 has been reported to inhibit salicylic acid (SA) biosynthesis and was upregulated in MeGI-overexpressing A. thaliana flowers, suggesting that MeGI represses the SA level by increasing the expression level of WRKY28. This was confirmed that SA level was lower in D. oleifera female floral buds than male. Overall, our findings indicate that the MeGI mediates its sex control function in D. oleifera mainly by regulating genes in the circadian rhythm, SA biosynthetic, and flavonoid biosynthetic pathways.

Evolutionary Landscape of Tea Circular RNAs and Its Contribution to Chilling Tolerance of Tea Plant.

Chilling stress threatens the yield and distribution pattern of global crops, including the tea plant (Camellia sinensis), one of the most important cash crops around the world. Circular RNA (circRNA) plays roles in regulating plant growth and biotic/abiotic stress responses. Understanding the evolutionary characteristics of circRNA and its feedbacks to chilling stress in the tea plant will help to elucidate the vital roles of circRNAs. In the current report, we systematically identified 2702 high-confidence circRNAs under chilling stress in the tea plant, and interestingly found that the generation of tea plant circRNAs was associated with the length of their flanking introns. Repetitive sequences annotation and DNA methylation analysis revealed that the longer flanking introns of circRNAs present more repetitive sequences and higher methylation levels, which suggested that repeat-elements-mediated DNA methylation might promote the circRNAs biogenesis in the tea plant. We further detected 250 differentially expressed circRNAs under chilling stress, which were functionally enriched in GO terms related to cold/stress responses. Constructing a circRNA-miRNA-mRNA interaction network discovered 139 differentially expressed circRNAs harboring potential miRNA binding sites, which further identified 14 circRNAs that might contribute to tea plant chilling responses. We further characterized a key circRNA, CSS-circFAB1, which was significantly induced under chilling stress. FISH and silencing experiments revealed that CSS-circFAB1 was potentially involved in chilling tolerance of the tea plant. Our study emphasizes the importance of circRNA and its preliminary role against low-temperature stress, providing new insights for tea plant cold tolerance breeding.

Integrating metabolite and transcriptome analysis revealed the different mechanisms of characteristic compound biosynthesis and transcriptional regulation in tea flowers.

The flowers of tea plants (Camellia sinensis), as well as tea leaves, contain abundant secondary metabolites and are big potential resources for the extraction of bioactive compounds or preparation of functional foods. However, little is known about the biosynthesis and transcriptional regulation mechanisms of those metabolites in tea flowers, such as terpenoid, flavonol, catechins, caffeine, and theanine. This study finely integrated target and nontarget metabolism analyses to explore the metabolic feature of developing tea flowers. Tea flowers accumulated more abundant terpenoid compounds than young leaves. The transcriptome data of developing flowers and leaves showed that a higher expression level of later genes of terpenoid biosynthesis pathway, such as Terpene synthases gene family, in tea flowers was the candidate reason of the more abundant terpenoid compounds than in tea leaves. Differently, even though flavonol and catechin profiling between tea flowers and leaves was similar, the gene family members of flavonoid biosynthesis were selectively expressed by tea flowers and tea leaves. Transcriptome and phylogenetic analyses indicated that the regulatory mechanism of flavonol biosynthesis was perhaps different between tea flowers and leaves. However, the regulatory mechanism of catechin biosynthesis was perhaps similar between tea flowers and leaves. This study not only provides a global vision of metabolism and transcriptome in tea flowers but also uncovered the different mechanisms of biosynthesis and transcriptional regulation of those important compounds.

Effects of Different Chemicals on Sexual Regulation in Persimmon (Diospyros kaki Thunb.) Flowers.

Research on crop sexuality is important for establishing systems for germplasm innovation and cultivating improved varieties. In this study, androecious persimmon trees were treated with various concentrations of ethrel (100, 500, and 1,000 mg/L) and zeatin (1, 5, and 10 mg/L) to investigate the morphological, physiological, and molecular characteristics of persimmon. Ethrel at 1,000 mg/L and zeatin at 10 mg/L both significantly reduced the stamen length and pollen grain diameter in androecious trees. Ethrel treatment also led to reduced stamen development with degenerated cellular contents; zeatin treatment promoted the development of arrested pistils via maintaining relatively normal mitochondrial morphology. Both treatments altered carbohydrate, amino acid, and endogenous phytohormone contents, as well as genes associated with hormone production and floral organ development. Thereafter, we explored the combined effects of four chemicals, including ethrel and zeatin, as well as zebularine and 5-azacytidine, both of which are DNA methylation inhibitors, on androecious persimmon flower development. Morphological comparisons showed that stamen length, pollen viability, and pollen grain diameter were significantly inhibited after combined treatment. Large numbers of genes involving in carbohydrate metabolic, mitogen-activated protein kinase (MAPK) signaling, and ribosome pathways, and metabolites including uridine monophosphate (UMP) and cyclamic acid were identified in response to the treatment, indicating complex regulatory mechanisms. An association analysis of transcriptomic and metabolomic data indicated that ribosomal genes have distinct effects on UMP and cyclamic acid metabolites, explaining how male floral buds of androecious persimmon trees respond to these exogenous chemicals. These findings extend the knowledge concerning sexual differentiation in persimmon; they also provide a theoretical basis for molecular breeding, high-yield cultivation, and quality improvement in persimmon.

Comparative transcriptomic analysis unveils the deep phylogeny and secondary metabolite evolution of 116 Camellia plants.

Camellia plants include more than 200 species of great diversity and immense economic, ornamental, and cultural values. We sequenced the transcriptomes of 116 Camellia plants from almost all sections of the genus Camellia. We constructed a pan-transcriptome of Camellia plants with 89 394 gene families and then resolved the phylogeny of genus Camellia based on 405 high-quality low-copy core genes. Most of the inferred relationships are well supported by multiple nuclear gene trees and morphological traits. We provide strong evidence that Camellia plants shared a recent whole genome duplication event, followed by large expansions of transcription factor families associated with stress resistance and secondary metabolism. Secondary metabolites, particularly those associated with tea quality such as catechins and caffeine, were preferentially heavily accumulated in the Camellia plants from section Thea. We thoroughly examined the expression patterns of hundreds of genes associated with tea quality, and found that some of them exhibited significantly high expression and correlations with secondary metabolite accumulations in Thea species. We also released a web-accessible database for efficient retrieval of Camellia transcriptomes. The reported transcriptome sequences and obtained novel findings will facilitate the efficient conservation and utilization of Camellia germplasm towards a breeding program for cultivated tea, camellia, and oil-tea plants.

Metabolite Profiling and Transcriptome Analysis Revealed the Conserved Transcriptional Regulation Mechanism of Caffeine Biosynthesis in Tea and Coffee Plants.

Caffeine is a characteristic bioactive compound in tea and coffee plants, which is synthesized and accumulated extensively in leaves and seeds. However, little is known about the regulatory mechanism of caffeine synthesis in plants. This study compared the caffeine metabolite between tea and coffee plants. We found that tea leaves contained significantly higher caffeine than coffee leaves, which is perhaps due to more members of N-methyltransferase (NMT) genes as well as higher expression levels in tea plants. Substantial numbers of transcription factors were predicted to be involved in caffeine biosynthesis regulation, combining weighted gene co-expression network analysis and the cis-element of NMT promoter analysis in tea and coffee plants. Furthermore, analysis of the transcription factors from the caffeine-related modules suggested that the regulatory mechanism of caffeine biosynthesis was probably partly conservative in tea and coffee plants. This study provides an essential resource for the regulatory mechanism of caffeine biosynthesis in plants.

Divergent DNA methylation contributes to duplicated gene evolution and chilling response in tea plants.

The tea plant (Camellia sinensis) is a thermophilic cash crop and contains a highly duplicated and repeat-rich genome. It is still unclear how DNA methylation regulates the evolution of duplicated genes and chilling stress in tea plants. We therefore generated a single-base-resolution DNA methylation map of tea plants under chilling stress. We found that, compared with other plants, the tea plant genome is highly methylated in all three sequence contexts, including CG, CHG and CHH (where H = A, T, or C), which is further proven to be correlated with its repeat content and genome size. We show that DNA methylation in the gene body negatively regulates the gene expression of tea plants, whereas non-CG methylation in the flanking region enables a positive regulation of gene expression. We demonstrate that transposable element-mediated methylation dynamics significantly drives the expression divergence of duplicated genes in tea plants. The DNA methylation and expression divergence of duplicated genes in the tea plant increases with evolutionary age and selective pressure. Moreover, we detect thousands of differentially methylated genes, some of which are functionally associated with chilling stress. We also experimentally reveal that DNA methyltransferase genes of tea plants are significantly downregulated, whereas demethylase genes are upregulated at the initial stage of chilling stress, which is in line with the significant loss of DNA methylation of three well-known cold-responsive genes at their promoter and gene body regions. Overall, our findings underscore the importance of DNA methylation regulation and offer new insights into duplicated gene evolution and chilling tolerance in tea plants.

The Reference Genome of Tea Plant and Resequencing of 81 Diverse Accessions Provide Insights into Its Genome Evolution and Adaptation.

Tea plant is an important economic crop, which is used to produce the world's oldest and most widely consumed tea beverages. Here, we present a high-quality reference genome assembly of the tea plant (Camellia sinensis var. sinensis) consisting of 15 pseudo-chromosomes. LTR retrotransposons (LTR-RTs) account for 70.38% of the genome, and we present evidence that LTR-RTs play critical roles in genome size expansion and the transcriptional diversification of tea plant genes through preferential insertion in promoter regions and introns. Genes, particularly those coding for terpene biosynthesis proteins, associated with tea aroma and stress resistance were significantly amplified through recent tandem duplications and exist as gene clusters in tea plant genome. Phylogenetic analysis of the sequences of 81 tea plant accessions with diverse origins revealed three well-differentiated tea plant populations, supporting the proposition for the southwest origin of the Chinese cultivated tea plant and its later spread to western Asia through introduction. Domestication and modern breeding left significant signatures on hundreds of genes in the tea plant genome, particularly those associated with tea quality and stress resistance. The genomic sequences of the reported reference and resequenced tea plant accessions provide valuable resources for future functional genomics study and molecular breeding of improved cultivars of tea plants.

Transcriptomic analyses reveal a systemic defense role of the uninfested adjacent leaf in tea plant (Camellia sinensis) attacked by tea geometrids (Ectropis obliqua).

To get a more detailed understanding of the interaction between tea plant (Camellia sinensis) and tea geometrids (Ectropis obliqua), transcriptomic profile in undamaged adjacent leaf (TGL) of tea geometrids fed local leaves (LL) was investigated for the first time. Here, approximately 245 million clean reads contained 39.39 Gb of sequence data were obtained from TGL. Further analysis revealed that systemic response was induced in TGL after tea geometrids feeding on LL, although the defense response was weaker than that in LL. The differentially expressed genes (DEGs) identification analysis showed little overlap of DEGs between TGL and LL. Comparative transcriptome analysis suggested that JA signal regulated resistant pathway was induced in LL; whereas primary metabolism pathway was activated in TGL in response to tea geometrids feeding. This study reveals a novel resistance mechanism of TGL to tea geometrids feeding.