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.

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.

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.

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.

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.

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.

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.

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.

Theanine transporters identified in tea plants (Camellia sinensis L.).

Theanine, a unique non-proteinogenic amino acid, is an important component of tea, as it confers the umami taste and relaxation effect of tea as a beverage. Theanine is primarily synthesized in tea roots and is subsequently transported to young shoots, which are harvested for tea production. Currently, the mechanism for theanine transport in the tea plant remains unknown. Here, by screening a yeast mutant library, followed by functional analyses, we identified the glutamine permease, GNP1 as a specific transporter for theanine in yeast. Although there is no GNP1 homolog in the tea plant, we assessed the theanine transport ability of nine tea plant amino acid permease (AAP) family members, with six exhibiting transport activity. We further determined that CsAAP1, CsAAP2, CsAAP4, CsAAP5, CsAAP6, and CsAAP8 exhibited moderate theanine affinities and transport was H+ -dependent. The tissue-specific expression of these six CsAAPs in leaves, vascular tissues, and the root suggested their broad roles in theanine loading and unloading from the vascular system, and in targeting to sink tissues. Furthermore, expression of these CsAAPs was shown to be seasonally regulated, coincident with theanine transport within the tea plant. Finally, CsAAP1 expression in the root was highly correlated with root-to-bud transport of theanine, in seven tea plant cultivars. Taken together, these findings support the hypothesis that members of the CsAAP family transport theanine and participate in its root-to-shoot delivery in the tea plant.

Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality.

Tea, one of the world's most important beverage crops, provides numerous secondary metabolites that account for its rich taste and health benefits. Here we present a high-quality sequence of the genome of tea, Camellia sinensis var. sinensis (CSS), using both Illumina and PacBio sequencing technologies. At least 64% of the 3.1-Gb genome assembly consists of repetitive sequences, and the rest yields 33,932 high-confidence predictions of encoded proteins. Divergence between two major lineages, CSS and Camellia sinensis var. assamica (CSA), is calculated to ∼0.38 to 1.54 million years ago (Mya). Analysis of genic collinearity reveals that the tea genome is the product of two rounds of whole-genome duplications (WGDs) that occurred ∼30 to 40 and ∼90 to 100 Mya. We provide evidence that these WGD events, and subsequent paralogous duplications, had major impacts on the copy numbers of secondary metabolite genes, particularly genes critical to producing three key quality compounds: catechins, theanine, and caffeine. Analyses of transcriptome and phytochemistry data show that amplification and transcriptional divergence of genes encoding a large acyltransferase family and leucoanthocyanidin reductases are associated with the characteristic young leaf accumulation of monomeric galloylated catechins in tea, while functional divergence of a single member of the glutamine synthetase gene family yielded theanine synthetase. This genome sequence will facilitate understanding of tea genome evolution and tea metabolite pathways, and will promote germplasm utilization for breeding improved tea varieties.

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.

Optimized sequencing depth and de novo assembler for deeply reconstructing the transcriptome of the tea plant, an economically important plant species.

BACKGROUND: Tea is the oldest and among the world's most popular non-alcoholic beverages, which has important economic, health and cultural values. Tea is commonly produced from the leaves of tea plants (Camellia sinensis), which belong to the genus Camellia of family Theaceae. In the last decade, many studies have generated the transcriptomes of tea plants at different developmental stages or under abiotic and/or biotic stresses to investigate the genetic basis of secondary metabolites that determine tea quality. However, these results exhibited large differences, particularly in the total number of reconstructed transcripts and the quality of the assembled transcriptomes. These differences largely result from limited knowledge regarding the optimized sequencing depth and assembler for transcriptome assembly of structurally complex plant species genomes. RESULTS: We employed different amounts of RNA-sequencing data, ranging from 4 to 84 Gb, to assemble the tea plant transcriptome using five well-known and representative transcript assemblers. Although the total number of assembled transcripts increased with increasing sequencing data, the proportion of unassembled transcripts became saturated as revealed by plant BUSCO datasets. Among the five representative assemblers, the Bridger package shows the best performance in both assembly completeness and accuracy as evaluated by the BUSCO datasets and genome alignment. In addition, we showed that Bridger and BinPacker harbored the shortest runtimes followed by SOAPdenovo and Trans-ABySS. CONCLUSIONS: The present study compares the performance of five representative transcript assemblers and investigates the key factors that affect the assembly quality of the transcriptome of the tea plants. This study will be of significance in helping the tea research community obtain better sequencing and assembly of tea plant transcriptomes under conditions of interest and may thus help to answer major biological questions currently facing the tea industry.

Tea Plant Information Archive: a comprehensive genomics and bioinformatics platform for tea plant.

Tea is the world's widely consumed nonalcohol beverage with essential economic and health benefits. Confronted with the increasing large-scale omics-data set particularly the genome sequence released in tea plant, the construction of a comprehensive knowledgebase is urgently needed to facilitate the utilization of these data sets towards molecular breeding. We hereby present the first integrative and specially designed web-accessible database, Tea Plant Information Archive (TPIA; http://tpia.teaplant.org). The current release of TPIA employs the comprehensively annotated tea plant genome as framework and incorporates with abundant well-organized transcriptomes, gene expressions (across species, tissues and stresses), orthologs and characteristic metabolites determining tea quality. It also hosts massive transcription factors, polymorphic simple sequence repeats, single nucleotide polymorphisms, correlations, manually curated functional genes and globally collected germplasm information. A variety of versatile analytic tools (e.g. JBrowse, blast, enrichment analysis, etc.) are established helping users to perform further comparative, evolutionary and functional analysis. We show a case application of TPIA that provides novel and interesting insights into the phytochemical content variation of section Thea of genus Camellia under a well-resolved phylogenetic framework. The constructed knowledgebase of tea plant will serve as a central gateway for global tea community to better understand the tea plant biology that largely benefits the whole tea industry.

Characterization and Alternative Splicing Profiles of the Lipoxygenase Gene Family in Tea Plant (Camellia sinensis).

Oxylipins, including jasmonic acid (JA) and volatiles, are important for signaling in plants, and these are formed by the lipoxygenase (LOX) enzyme family. There is a large gap in understanding of the underlying molecular basis of their roles in tea plants. Here, we identified 11 CsLOX genes from the tea plant (Camellia sinensis), and characterized their phylogeny, gene structure and protein features into three subclasses. We then examined their enzymatic activities, LOX expression and alternative splicing of transcripts during development and in response to abiotic or biotic stresses in tea plants. In vitro expressed protein assays showed that the CsLOX2, 3 and 9 enzymatically function to produce 9/13-HPOT, 13-HPOT and 9-HPOT, respectively. CsLOX2 and CsLOX9 green fluorescent protein (GFP) fusion proteins localized to chloroplasts and the cytoplasm, respectively. RNA sequencing, quantitative reverse transcription-PCR and Northern blot analysis suggested that CsLOX5, 6 and 9 were predominantly expressed in seeds, flowers and roots, respectively. CsLOX2, 3, 4, 6 and 7 were up-regulated after attack by the insect Ectropis oblique, while CsLOX1 was induced after infection with the pathogen Glomerella cingulata. CsLOX3, 7 and 10 were up-regulated by JA but not ABA or salicylic acid. Long-term cold stress down-regulated CsLOX expression while a short duration of cold induced the expression of CsLOX1, 6 and 7. Alternatively spliced transcripts of six CsLOX genes were dynamically regulated through time and varied in relative abundances under the investigated stresses; we propose a mechanism of competing or compensating regulation between isoforms. This study improves our understanding of evolution of LOXs and regulation of their diverse functions in plants.

Circular RNA architecture and differentiation during leaf bud to young leaf development in tea (Camellia sinensis).

MAIN CONCLUSION: Circular RNA (circRNA) discovery, expression patterns and experimental validation in developing tea leaves indicates its correlation with circRNA-parental genes and potential roles in ceRNA interaction network. Circular RNAs (circRNAs) have recently emerged as a novel class of abundant endogenous stable RNAs produced by circularization with regulatory potential. However, identification of circRNAs in plants, especially in non-model plants with large genomes, is challenging. In this study, we undertook a systematic identification of circRNAs from different stage tissues of tea plant (Camellia sinensis) leaf development using rRNA-depleted circular RNA-seq. By combining two state-of-the-art detecting tools, we characterized 3174 circRNAs, of which 342 were shared by each approach, and thus considered high-confidence circRNAs. A few predicted circRNAs were randomly chosen, and 20 out of 24 were experimental confirmed by PCR and Sanger sequencing. Similar in other plants, tissue-specific expression was also observed for many C. sinensis circRNAs. In addition, we found that circRNA abundances were positively correlated with the mRNA transcript abundances of their parental genes. qRT-PCR validated the differential expression patterns of circRNAs between leaf bud and young leaf, which also indicated the low expression abundance of circRNAs compared to the standard mRNAs from the parental genes. We predicted the circRNA-microRNA interaction networks, and 54 of the differentially expressed circRNAs were found to have potential tea plant miRNA binding sites. The gene sets encoding circRNAs were significantly enriched in chloroplasts related GO terms and photosynthesis/metabolites biosynthesis related KEGG pathways, suggesting the candidate roles of circRNAs in photosynthetic machinery and metabolites biosynthesis during leaf development.

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.

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.

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.