The Chromosome-Level Reference Genome of Tea Tree Unveils Recent Bursts of Non-autonomous LTR Retrotransposons in Driving Genome Size Evolution

Evolution of the biochemistry underpinning purine alkaloid metabolism in plants

Purine alkaloids are naturally occurring nitrogenous methylated derivatives of purine nucleotide degradation products, having essential roles in medicine, food and various other aspects of our daily lives. They are generated through convergent evolution in different plant species. The pivotal reaction steps within the purine alkaloid metabolic pathways have been largely elucidated, and the convergent evolution of purine alkaloids has been substantiated through bioinformatic, biochemical and other research perspectives within S-adenosyl-ʟ-methionine-dependent N-methyltransferases. Currently, the biological and ecological roles of purine alkaloids, further refinement of the purine alkaloid metabolic pathways and the investigation of purine alkaloid adaptive evolutionary mechanisms continue to attract widespread research interest. The exploration of the purine alkaloid metabolic pathways also enhances our comprehension of the biochemical mechanism, providing insights into inter-species interactions and adaptive evolution and offering potential value in drug development and agricultural applications. Here, we review the progress of research in the distribution, metabolic pathway elucidation and regulation, evolutionary mechanism and ecological roles of purine alkaloids in plants. The opportunities and challenges involved in elucidating the biochemical basis and evolutionary mechanisms of the purine alkaloid metabolic pathways, as well as other research aspects, are also discussed. This article is part of the theme issue 'The evolution of plant meta-bolism'.

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.

Genome-wide epigenetic dynamics of tea leaves under mechanical wounding stress during oolong tea postharvest processing

Understanding the epigenetic responses to mechanical wounding stress during the postharvest processing of oolong tea provides insight into the reprogramming of the tea genome and its impact on tea quality. Here, we characterized the 5mC DNA methylation and chromatin accessibility landscapes of tea leaves subjected to mechanical wounding stress during the postharvest processing of oolong tea. Analysis of the differentially methylated regions and preferentially accessible promoters revealed many overrepresented TF-binding motifs, highlighting sets of TFs that are likely important for the quality of oolong tea. Within these sets, we constructed a chromatin accessibility-mediated gene regulatory network specific to mechanical wounding stress. In combination with the results of the TF-centred yeast one-hybrid assay, we identified potential binding sites of CsMYC2 and constructed a gene regulatory network centred on CsMYC2, clarifying the potential regulatory role of CsMYC2 in the postharvest processing of oolong tea. Interestingly, highly accessible chromatin and hypomethylated cytosine were found to coexist in the promoter region of the indole biosynthesis gene (tryptophan synthase β-subunit, CsTSB) under wounding stress, which indicates that these two important epigenetic regulatory mechanisms are jointly involved in regulating the synthesis of indole during the postharvest processing of oolong tea. These findings improve our understanding of the epigenetic regulatory mechanisms involved in quality formation during the postharvest processing of oolong tea.

Chromosome-Scale Genome Assembly and Characterization of Top-Quality Japanese Green Tea Cultivar 'Seimei'

Japanese green tea, an essential beverage in Japanese culture, is characterized by the initial steaming of freshly harvested leaves during production. This process efficiently inactivates endogenous enzymes such as polyphenol oxidases, resulting in the production of sencha, gyokuro and matcha that preserves the vibrant green color of young leaves. Although genome sequences of several tea cultivars and germplasms have been published, no reference genome sequences are available for Japanese green tea cultivars. Here, we constructed a reference genome sequence of the cultivar 'Seimei', which is used to produce high-quality Japanese green tea. Using the PacBio HiFi and Hi-C technologies for chromosome-scale genome assembly, we obtained 15 chromosome sequences with a total genome size of 3.1 Gb and an N50 of 214.9 Mb. By analyzing the genomic diversity of 23 Japanese tea cultivars and lines, including the leading green tea cultivars 'Yabukita' and 'Saemidori', it was revealed that several candidate genes could be related to the characteristics of Japanese green tea. The reference genome of 'Seimei' and information on genomic diversity of Japanese green tea cultivars should provide crucial information for effective breeding of such cultivars in the future.

Genomic analyses of agronomic traits in tea plants and related Camellia species

The genus Camellia contains three types of domesticates that meet various needs of ancient humans: the ornamental C. japonica, the edible oil-producing C. oleifera, and the beverage-purposed tea plant C. sinensis. The genomic drivers of the functional diversification of Camellia domesticates remain unknown. Here, we present the genomic variations of 625 Camellia accessions based on a new genome assembly of C. sinensis var. assamica ('YK10'), which consists of 15 pseudo-chromosomes with a total length of 3.35 Gb and a contig N50 of 816,948 bp. These accessions were mainly distributed in East Asia, South Asia, Southeast Asia, and Africa. We profiled the population and subpopulation structure in tea tree Camellia to find new evidence for the parallel domestication of C. sinensis var. assamica (CSA) and C. sinensis var. sinensis (CSS). We also identified candidate genes associated with traits differentiating CSA, CSS, oilseed Camellia, and ornamental Camellia cultivars. Our results provide a unique global view of the genetic diversification of Camellia domesticates and provide valuable resources for ongoing functional and molecular breeding research.

In-depth exploration of the genomic diversity in tea varieties based on a newly constructed pangenome of Camellia sinensis

Tea, one of the most widely consumed beverages globally, exhibits remarkable genomic diversity in its underlying flavour and health-related compounds. In this study, we present the construction and analysis of a tea pangenome comprising a total of 11 genomes, with a focus on three newly sequenced genomes comprising the purple-leaved assamica cultivar "Zijuan", the temperature-sensitive sinensis cultivar "Anjibaicha" and the wild accession "L618" whose assemblies exhibited excellent quality scores as they profited from latest sequencing technologies. Our analysis incorporates a detailed investigation of transposon complement across the tea pangenome, revealing shared patterns of transposon distribution among the studied genomes and improved transposon resolution with long read technologies, as shown by long terminal repeat (LTR) Assembly Index analysis. Furthermore, our study encompasses a gene-centric exploration of the pangenome, exploring the genomic landscape of the catechin pathway with our study, providing insights on copy number alterations and gene-centric variants, especially for Anthocyanidin synthases. We constructed a gene-centric pangenome by structurally and functionally annotating all available genomes using an identical pipeline, which both increased gene completeness and allowed for a high functional annotation rate. This improved and consistently annotated gene set will allow for a better comparison between tea genomes. We used this improved pangenome to capture the core and dispensable gene repertoire, elucidating the functional diversity present within the tea species. This pangenome resource might serve as a valuable resource for understanding the fundamental genetic basis of traits such as flavour, stress tolerance, and disease resistance, with implications for tea breeding programmes.

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.

Genome-wide identification and expression pattern analysis of WRKY transcription factors in response to biotic and abiotic stresses in tea plants (Camellia sinensis)

Plants would encounter various biotic and abiotic stresses during the growth and development. WRKY transcription factors (TFs) as plant-specific TFs, play an important role in responding to various adverse circumstances. Despite some advances were achieved in functional studies of WRKY TFs in tea plants, systematic analysis of the involvement of CsWRKY TFs when facing cold, salt, drought stresses and pathogen and insect attack was lacked. In present study, a total of 78 CsWRKY TFs were identified following the genomic and transcript databases. The expression patterns of CsWRKYs in various organs of tea plants and the expression profiles in response to biotic and abiotic stresses were investigated by examining representative RNA-seq data. Moreover, the effects of hormone treatments (SA and MeJA) on the transcription levels of WRKY TFs were also investigated. The phylogenetic tree of CsWRKY TFs from different species indicated the functional diversity of WRKY TFs was not closely related to their protein classification. Concurrently, CsWRKY70-2 TF was identified as a positive regulator in response to drought stress. This study provided solid and valuable information, helping us better understand the functional diversity of CsWRKY TFs, and laid the foundation for further research on the function of key WRKY genes in tea plants.

New insights into evolution and functional diversification of Camellia sinensis LRR-RLKs

Leucine-rich repeat receptor-like kinases (LRR-RLKs) represent the largest subgroup of receptor-like kinases (RLKs) in plants. While some LRR-RLK members play a role in regulating various plant growth processes related to morphogenesis, disease resistance, and stress response, the functions of most LRR-RLK genes remain unclear. In this study, we identified 397 LRR-RLK genes from the genome of Camellia sinensis and categorized them into 16 subfamilies. Approximately 62% of CsLRR-RLK genes are situated in regions resulting from segmental duplications, suggesting that the expansion of CsLRR-RLK genes is due to segmental duplications. Analysis of gene expression patterns revealed differential expression of CsLRR-RLK genes across different tissues and in response to stress. Furthermore, we demonstrated that CssEMS1 localizes to the cell membrane and can complement Arabidopsis ems1 mutant. This study is the initial in-depth evolutionary examination of LRR-RLKs in tea and provides a basis for future investigations into their functionality.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01458-1.

Depicting the genetic and metabolic panorama of chemical diversity in the tea plant

As a frequently consumed beverage worldwide, tea is rich in naturally important bioactive metabolites. Combining genetic, metabolomic and biochemical methodologies, here, we present a comprehensive study to dissect the chemical diversity in tea plant. A total of 2837 metabolites were identified at high-resolution with 1098 of them being structurally annotated and 63 of them were structurally identified. Metabolite-based genome-wide association mapping identified 6199 and 7823 metabolic quantitative trait loci (mQTL) for 971 and 1254 compounds in young leaves (YL) and the third leaves (TL), respectively. The major mQTL (i.e., P < 1.05 × 10-5, and phenotypic variation explained (PVE) > 25%) were further interrogated. Through extensive annotation of the tea metabolome as well as network-based analysis, this study broadens the understanding of tea metabolism and lays a solid foundation for revealing the natural variations in the chemical composition of the tea plant. Interestingly, we found that galloylations, rather than hydroxylations or glycosylations, were the largest class of conversions within the tea metabolome. The prevalence of galloylations in tea is unusual, as hydroxylations and glycosylations are typically the most prominent conversions of plant specialized metabolism. The biosynthetic pathway of flavonoids, which are one of the most featured metabolites in tea plant, was further refined with the identified metabolites. And we demonstrated the further mining and interpretation of our GWAS results by verifying two identified mQTL (including functional candidate genes CsUGTa, CsUGTb, and CsCCoAOMT) and completing the flavonoid biosynthetic pathway of the tea plant.

The miR166 targets CsHDZ3 genes to negatively regulate drought tolerance in tea plant (Camellia sinensis)

Drought is the stressor with a significant adverse impact on the yield stability of tea plants. HD-ZIP III transcription factors (TFs) play important regulatory roles in plant growth, development, and stress responses. However, whether and how HD-ZIP III TFs are involved in drought response and tolerance in tea plants remains unclear. Here, we identified seven HD-ZIP III genes (CsHDZ3-1 to CsHDZ3-7) in tea plant genome. The evolutionary analysis demonstrated that CsHDZ3 members were subjected to purify selection. Subcellular localization analysis revealed that all seven CsHDZ3s located in the nucleus. Yeast self-activation and dual-luciferase reporter assays demonstrated that CsHDZ3-1 to CsHDZ3-4 have trans-activation ability whereas CsHDZ3-5 to CsHDZ3-7 served as transcriptional inhibitors. The qRT-PCR assay showed that all seven CsHDZ3 genes could respond to simulated natural drought stress and polyethylene glycol treatment. Further assays verified that all CsHDZ3 genes can be cleaved by csn-miR166. Overexpression of csn-miR166 inhibited the expression of seven CsHDZ3 genes and weakened drought tolerance of tea leaves. In contrast, suppression of csn-miR166 promoted the expression of seven CsHDZ3 genes and enhanced drought tolerance of tea leaves. These findings established the foundation for further understanding the mechanism of CsHDZ3-miR166 modules' participation in drought responses and tolerance.

Facile General Injectable Gelatin/Metal/Tea Polyphenol Double Nanonetworks Remodel Wound Microenvironment and Accelerate Healing

Treating the most widespread complication of diabetes: diabetic wounds poses a significant clinical obstacle due to the intricate nature of wound healing in individuals with diabetes. Here a novel approach is proposed using easily applicable injectable gelatin/metal/tea polyphenol double nanonetworks, which effectively remodel the wound microenvironment and accelerates the healing process. The gelatin(Gel) crosslink with metal ions (Zr4+ ) through the amino acids, imparting advantageous mechanical properties like self-healing, injectability, and adhesion. The nanonetwork's biological functions are further enhanced by incorporating the tea polyphenol metal nanonetwork through in situ doping of the epigallocatechin gallate (EGCG) with great antibacterial, self-healing, antioxidant, and anticancer capabilities. The in vitro and in vivo tests show that this double nanonetworks hydrogel exhibits faster cell migration and favorable anti-inflammatory and antioxidant properties and can greatly reshape the microenvironment of diabetic wounds and accelerate the wound healing rate. In addition, this hydrogel is completely degraded after subcutaneous injection for 7 days, with nondetectable cytotoxicity in H&E staining of major mice organs and the serum level of liver function indicators. Considering the above-mentioned merits of this hydrogel, it is believed that the injectable gelatin/metal/tea polyphenol double nanonetworks have broad biomedical potential, especially in diabetic wound repair and tissue engineering.

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.

Multiomics Analysis Reveals the Genetic Basis of Volatile Terpenoid Formation in Oolong Tea

Oolong tea has gained great popularity in China due to its pleasant floral and fruity aromas. Although numerous studies have investigated the aroma differences across various tea cultivars, the genetic mechanism is unclear. This study performed multiomics analysis of three varieties suitable for oolong tea and three others with different processing suitability. Our analysis revealed that oolong tea varieties contained higher levels of cadinane sesquiterpenoids. PanTFBS was developed to identify variants of transcription factor binding sites (TFBSs). We found that the CsDCS gene had two TFBS variants in the promoter sequence and a single nucleotide polymorphism (SNP) in the coding sequence. Integrating data on genetic variations, gene expression, and protein-binding sites indicated that CsDCS might be a pivotal gene involved in the biosynthesis of cadinane sesquiterpenoids. These findings advance our understanding of the genetic factors involved in the aroma formation of oolong tea and offer insights into the enhancement of tea aroma.

Gene mining and genomics-assisted breeding empowered by the pangenome of tea plant Camellia sinensis

Tea is one of the world's oldest crops and is cultivated to produce beverages with various flavours. Despite advances in sequencing technologies, the genetic mechanisms underlying key agronomic traits of tea remain unclear. In this study, we present a high-quality pangenome of 22 elite cultivars, representing broad genetic diversity in the species. Our analysis reveals that a recent long terminal repeat burst contributed nearly 20% of gene copies, introducing functional genetic variants that affect phenotypes such as leaf colour. Our graphical pangenome improves the efficiency of genome-wide association studies and allows the identification of key genes controlling bud flush timing. We also identified strong correlations between allelic variants and flavour-related chemistries. These findings deepen our understanding of the genetic basis of tea quality and provide valuable genomic resources to facilitate its genomics-assisted breeding.

CRISPR/Cas9 Editing Sites Identification and Multi-Elements Association Analysis in Camellia sinensis

CRISPR/Cas9 is an efficient genome-editing tool, and the identification of editing sites and potential influences in the Camellia sinensis genome have not been investigated. In this study, bioinformatics methods were used to characterise the Camellia sinensis genome including editing sites, simple sequence repeats (SSRs), G-quadruplexes (GQ), gene density, and their relationships. A total of 248,134,838 potential editing sites were identified in the genome, and five PAM types, AGG, TGG, CGG, GGG, and NGG, were observed, of which 66,665,912 were found to be specific, and they were present in all structural elements of the genes. The characteristic region of high GC content, GQ density, and PAM density in contrast to low gene density and SSR density was identified in the chromosomes in the joint analysis, and it was associated with secondary metabolites and amino acid biosynthesis pathways. CRISPR/Cas9, as a technology to drive crop improvement, with the identified editing sites and effector elements, provides valuable tools for functional studies and molecular breeding in Camellia sinensis.

Plant exosome nanovesicles (PENs): green delivery platforms

Natural plants have been attracting increasing attention in biomedical research due to their numerous benefits. Plant exosome-derived vesicles, some of the plant's components, are small nanoscale vesicles secreted by plant cells. These vesicles are rich in bioactive substances and play significant roles in intercellular communication, information transfer, and maintaining homeostasis in organisms. They also hold promise for treating diseases, and their vesicular structures make them suitable carriers for drug delivery, with large-scale production feasible. Therefore, this paper aims to provide an overview of nanovesicles from different plant sources and their extraction methods. We also outline the biological activities of nanovesicles, including their anti-inflammatory, anti-viral, and anti-tumor properties, and systematically introduce their applications in drug delivery. These applications include transdermal delivery, targeted drug delivery, gene delivery, and their potential use in the modern food industry. This review provides new ideas and methods for future research on plant exosomes, including their empowerment by artificial intelligence and gene editing, as well as their potential application in the biomedicine, food, and agriculture industries.

The complete chloroplast genome and phylogenomic analysis of Camellia sinensis var. sinensis cultivar 'Liupao', a landrace from Guangxi, China

Liupao tea is one of the well-known Chinese tea brands and a famous local specialty in Wuzhou, Guangxi, China. However, the genetic background and phylogenetic relationship of the native resource plants of Liupao tea need study, especially at the genomic level. In this study, we reported the complete chloroplast (cp) genome sequence of Camellia sinensis var. sinensis cultivar 'Liupao' (LP, Liupao tea population) and inferred its phylogenetic relationship to other tea plant variants or cultivars. The cp genome had a total length of 157,097 bp and the overall GC content was 37.3%. The cp genome contained one LSC region (86,641 bp) and one SSC region (18,276 bp), which were separated by two IR regions (26,090 bp, respectively). Moreover, the cp genomes were composed of 130 genes, including 86 protein-coding genes, 36 tRNA genes, and eight rRNA genes. The phylogenetic analysis showed that LP was closely related to C. sinensis var. pabilimba cv. 'Lingyunbaihao'. This study will provide useful information for further investigating the genetic background, evolution, and breeding of LP as well as other tea cultivars and varieties.

Characterization of two O-methyltransferases involved in the biosynthesis of O-methylated catechins in tea plant

Tea is known for having a high catechin content, with the main component being (-)-epigallocatechin gallate (EGCG), which has significant bioactivities, including potential anti-cancer and anti-inflammatory activity. The poor intestinal stability and permeability of EGCG, however, undermine these health-improving benefits. O-methylated EGCG derivatives, found in a few tea cultivars in low levels, have attracted considerable interest due to their increased bioavailability. Here, we identify two O-methyltransferases from tea plant: CsFAOMT1 that has a specific O-methyltransferase activity on the 3''-position of EGCG to generate EGCG3''Me, and CsFAOMT2 that predominantly catalyzes the formation of EGCG4″Me. In different tea tissues and germplasms, the transcript levels of CsFAOMT1 and CsFAOMT2 are strongly correlated with the amounts of EGCG3''Me and EGCG4''Me, respectively. Furthermore, the crystal structures of CsFAOMT1 and CsFAOMT2 reveal the key residues necessary for 3''- and 4''-O-methylation. These findings may provide guidance for the future development of tea cultivars with high O-methylated catechin content.

Structural and Functional Analysis of the MADS-Box Genes Reveals Their Functions in Cold Stress Responses and Flower Development in Tea Plant (Camellia sinensis)

MADS-box genes comprise a large family of transcription factors that play crucial roles in all aspects of plant growth and development. However, no detailed information on the evolutionary relationship and functional characterization of MADS-box genes is currently available for some representative lineages, such as the Camellia plant. In this study, 136 MADS-box genes were detected from a reference genome of the tea plant (Camellia sinensis) by employing a 569 bp HMM (Hidden Markov Model) developed using nucleotide sequencing including 73 type I and 63 type II genes. An additional twenty-seven genes were identified, with five MIKC-type genes. Truncated and/or inaccurate gene models were manually verified and curated to improve their functional characterization. Subsequently, phylogenetic relationships, chromosome locations, conserved motifs, gene structures, and gene expression profiles were systematically investigated. Tea plant MIKC genes were divided into all 14 major eudicot subfamilies, and no gene was found in Mβ. The expansion of MADS-box genes in the tea plant was mainly contributed by WGD/fragment and tandem duplications. The expression profiles of tea plant MADS-box genes in different tissues and seasons were analyzed, revealing widespread evolutionary conservation and genetic redundancy. The expression profiles linked to cold stress treatments suggested the wide involvement of MADS-box genes from the tea plant in response to low temperatures. Moreover, a floral 'ABCE' model was proposed in the tea plant and proved to be both conserved and ancient. Our analyses offer a detailed overview of MADS-box genes in the tea plant, allowing us to hypothesize the potential functions of unknown genes and providing a foundation for further functional characterizations.

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.

Genome-wide characterization of long terminal repeat retrotransposons provides insights into trait evolution of four cucurbit species

Cucurbits are a diverse plant family that includes economically important crops, such as cucumber, watermelon, melon, and pumpkin. Knowledge of the roles that long terminal repeat retrotransposons (LTR-RTs) have played in diversification of cucurbit species is limited; to add to understanding of the roles of LTR-RTs, we assessed their distributions in four cucurbit species. We identified 381, 578, 1086, and 623 intact LTR-RTs in cucumber (Cucumis sativus L. var. sativus cv. Chinese Long), watermelon (Citrullus lanatus subsp. vulgaris cv. 97103), melon (Cucumis melo cv. DHL92), and Cucurbita (Cucurbita moschata var. Rifu), respectively. Among these LTR-RTs, the Ale clade of the Copia superfamily was the most abundant in all the four cucurbit species. Insertion time and copy number analysis revealed that an LTR-RT burst occurred approximately 2 million years ago in cucumber, watermelon, melon, and Cucurbita, and may have contributed to their genome size variation. Phylogenetic and nucleotide polymorphism analyses suggested that most LTR-RTs were formed after species diversification. Analysis of gene insertions by LTR-RTs revealed that the most frequent insertions were of Ale and Tekay and that genes related to dietary fiber synthesis were the most commonly affected by LTR-RTs in Cucurbita. These results increase our understanding of LTR-RTs and their roles in genome evolution and trait characterization in cucurbits.

Insights into the convergent evolution of fructan biosynthesis in angiosperms from the highly characteristic chicory genome

Fructans in angiosperms play essential roles in physiological functions and environmental adaptations. As a major source of industrial fructans (especially inulin-type), chicory (Cichorium intybus L.) is a model species for studying fructan biosynthesis. However, the genes underlying this process and their evolutionary history in angiosperms remain elusive. We combined multiple sequencing technologies to assemble and annotate the chicory genome and scan its (epi)genomic features, such as genomic components, DNA methylation, and three-dimensional (3D) structure. We also performed a comparative genomics analysis to uncover the associations between key traits and gene families. We achieved a nearly complete chicory genome assembly and found that continuous bursts of a few highly active retrotransposon families largely shaped the (epi)genomic characteristics. The highly methylated genome with its unique 3D structure potentially influences critical biological processes. Our comprehensive comparative genomics analysis deciphered the genetic basis for the rich sesquiterpene content in chicory and indicated that the fructan-accumulating trait resulted from convergent evolution in angiosperms due to shifts in critical sites of fructan-active enzymes. The highly characterized chicory genome provides insight into Asteraceae evolution and fructan biosynthesis in angiosperms.

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.

A reference-grade genome assembly for Astragalus mongholicus and insights into the biosynthesis and high accumulation of triterpenoids and flavonoids in its roots

Astragalus membranaceus var. mongholicus (AMM), a member of the Leguminosae, is one of the most important medicinal plants worldwide. The dried roots of AMM have a wide range of pharmacological effects and are a traditional Chinese medicine. Here, we report the first chromosome-level reference genome of AMM, comprising nine pseudochromosomes with a total size of 1.47 Gb and 27 868 protein-encoding genes. Comparative genomic analysis reveals that AMM has not experienced an independent whole-genome duplication (WGD) event after the WGD event shared by the Papilionoideae species. Analysis of long terminal repeat retrotransposons suggests a recent burst of these elements at approximately 0.13 million years ago, which may explain the large size of the AMM genome. Multiple gene families involved in the biosynthesis of triterpenoids and flavonoids were expanded, and our data indicate that tandem duplication has been the main driver for expansion of these families. Among the expanded families, the phenylalanine ammonia-lyase gene family was primarily expressed in the roots of AMM, suggesting their roles in the biosynthesis of phenylpropanoid compounds. The functional versatility of 2,3-oxidosqualene cyclase genes in cluster III may play a critical role in the diversification of triterpenoids in AMM. Our findings provide novel insights into triterpenoid and flavonoid biosynthesis and can facilitate future research on the genetics and medical applications of AMM.

Molecular regulation of immunity in tea plants

Tea, which is mainly produced using the young leaves and buds of tea plants (Camellia sinensis (L.) O. Kuntze), is one of the most common non-alcoholic beverages consumed in the world. The standard of tea mostly depends on the variety and quality of tea plants, which generally grow in subtropical areas, where the warm and humid conditions are also conducive to the occurrence of diseases. In fighting against pathogens, plants rely on their sophisticated innate immune systems which has been extensively studied in model plants. Many components involved in pathogen associated molecular patterns (PAMPs) triggered immunity (PTI) and effector triggered immunity (ETI) have been found. Nevertheless, the molecular regulating network against pathogens (e.g., Pseudopestalotiopsis sp., Colletotrichum sp. and Exobasidium vexans) causing widespread disease (such as grey blight disease, anthracnose, and blister blight) in tea plants is still unclear. With the recent release of the genome data of tea plants, numerous genes involved in tea plant immunity have been identified, and the molecular mechanisms behind tea plant immunity is being studied. Therefore, the recent achievements in identifying and cloning functional genes/gene families, in finding crucial components of tea immunity signaling pathways, and in understanding the role of secondary metabolites have been summarized and the opportunities and challenges in the future studies of tea immunity are highlighted in this review.

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.

Transcriptome-based variations effectively untangling the intraspecific relationships and selection signals in Xinyang Maojian tea population

As one of the world's top three popular non-alcoholic beverages, tea is economically and culturally valuable. Xinyang Maojian, this elegant green tea, is one of the top ten famous tea in China and has gained prominence for thousands of years. However, the cultivation history of Xinyang Maojian tea population and selection signals of differentiation from the other major variety Camellia sinensis var. assamica (CSA) remain unclear. We newly generated 94 Camellia sinensis (C. sinensis) transcriptomes including 59 samples in the Xinyang area and 35 samples collected from 13 other major tea planting provinces in China. Comparing the very low resolution of phylogeny inferred from 1785 low-copy nuclear genes with 94 C. sinensis samples, we successfully resolved the phylogeny of C. sinensis samples by 99,115 high-quality SNPs from the coding region. The sources of tea planted in the Xinyang area were extensive and complex. Specifically, Shihe District and Gushi County were the two earliest tea planting areas in Xinyang, reflecting a long history of tea planting. Furthermore, we identified numerous selection sweeps during the differentiation of CSA and CSS and these positive selection genes are involved in many aspects such as regulation of secondary metabolites synthesis, amino acid metabolism, photosynthesis, etc. Numerous specific selective sweeps of modern cultivars were annotated with functions in various different aspects, indicating the CSS and CSA populations possibly underwent independent specific domestication processes. Our study indicated that transcriptome-based SNP-calling is an efficient and cost-effective method in untangling intraspecific phylogenetic relationships. This study provides a significant understanding of the cultivation history of the famous Chinese tea Xinyang Maojian and unravels the genetic basis of physiological and ecological differences between the two major tea subspecies.

Genome-Wide Identification of AMT2-Type Ammonium Transporters Reveal That CsAMT2.2 and CsAMT2.3 Potentially Regulate NH4+ Absorption among Three Different Cultivars of Camellia sinensis

Ammonium (NH₄⁺), as a major inorganic source of nitrogen (N) for tea plant growth, is transported and distributed across membranes by the proteins of ammonium transporters (AMTs). However, the AMT2-type AMTs from tea plants remain poorly understood. In this study, five CsAMT2 subfamily genes were identified in tea plant genomes, and their full-length coding sequences (CDS) were isolated from roots. Then, a NH₄⁺ uptake kinetic comparison of Fudingdabaicha (FD), Huangdan (HD), and Maoxie (MX) showed that FD was a high N efficiency (HNE) cultivar that had a wide range of adaptability to NH₄⁺, HD was a high N efficiency under high N conditions (HNEH) cultivar, in which it was easy to obtain higher yield in a high N environment, and MX was a high N efficiency under low N conditions (HNEL) cultivar, which had a higher affinity for NH₄⁺ than the other two. Tissue-specific expression analysis suggested that CsAMT2.2 and CsAMT2.3 were highly expressed in the roots, indicating that these two members may be unique in the CsAMT2 subfamily. This is further supported by our findings from the temporal expression profiles in the roots among these three different N adaptation cultivars. Expression levels of CsAMT2.2 and CsAMT2.3 in FD and HD were upregulated by a short time (2 h) under high NH₄⁺ treatment, while under low NH₄⁺ treatment, CsAMT2.2 and CsAMT2.3 were highly expressed at 0 h and 2 h in the HNEL-type cultivar-MX. Furthermore, the functional analysis illustrated that CsAMT2.2 and CsAMT2.3 could make a functional complementation of NH₄⁺-defective mutant yeast cells at low NH₄⁺ levels, and the transport efficiency of CsAMT2.3 was higher than that of CsAMT2.2. Thus, we concluded that CsAMT2.2 and CsAMT2.3 might play roles in controlling the NH₄⁺ uptake from the soil to the roots. These results will further the understanding of the NH₄⁺ signal networks of AMT2-type proteins in tea plants.

Improving Plant Photosynthesis through Light-Harvesting Upconversion Nanoparticles

Nanotechnology is considered as an emerging effective means to augment plant photosynthesis. However, there is still a lot of work to be done in this field. Here, we applied the upconversion nanoparticles (UCNPs) on lettuce leaves and found that the UCNPs were able to transport into the lettuce body and colocalize with the chloroplasts. It was proved that UCNPs could harvest the near-infrared light of sunlight and increase the electron transfer rate in the photosynthesis process, thus increasing the photosynthesis rate. The gene expression analysis showed that more than 90% of gene expression in photosynthesis was upregulated. After spraying the UCNP solution on the leaves of lettuce and placing the lettuce under sunlight for 1 week, the wet/dry weight of the leaves increased by 53.33% and 45.71%, respectively. This nanoengineering of light-harvesting UCNPs may have great potential for applications in agriculture.

TeaPVs: a comprehensive genomic variation database for tea plant (Camellia sinensis)

Genome variation not only plays an important role in plant phenotypic modeling and adaptive evolution, but also enhances population genetic diversity and regulates gene expression. The tea tree (Camellia sinensis) has a large genome (~ 3.0 Gb), making the identification of genome-wide variants time-consuming and expensive. With the continuous publication of a large number of different types of population sequencing data, there is a lack of an open platform to integrate these data and identify variants in the tea plant genome.To integrate the genetic variation confidence in the tea plant population genome, 238 whole-genome resequencing, 213 transcriptome sequencing, and 96 hybrid F1 individuals with a total of more than 20 Tb were collected for mutation site identification. Based on these variations information, we constructed the first tea tree variation web service database TeaPVs ( http://47.106.184.91:8025/ and http://liushang.top:8025/ ). It supports users to search all SNP, Indel, SV mutations and SSR/Polymorphic SSR sequences by location or gene ID. Furthermore, the website also provides the functions of gene expression search of different transcriptome, sequence blast, sequence extraction of CDS and mutation loci, etc.The features of the TeaPVs database make it a comprehensive tea plant genetic variation bioinformatics platform for researchers, and will also be helpful for revealing new functional mutations in the tea plant genome and molecular marker-assisted breeding.

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.

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.

United States tea: A synopsis of ongoing tea research and solutions to United States tea production issues

Tea is a steeped beverage made from the leaves of Camellia sinensis. Globally, this healthy, caffeine-containing drink is one of the most widely consumed beverages. At least 50 countries produce tea and most of the production information and tea research is derived from international sources. Here, we discuss information related to tea production, genetics, and chemistry as well as production issues that affect or are likely to affect emerging tea production and research in the United States. With this review, we relay current knowledge on tea production, threats to tea production, and solutions to production problems to inform this emerging market in the United States.

Identification of gene family members and a key structural variation reveal important roles of OVATE genes in regulating tea (Camellia sinensis) leaf development

OVATE genes are a new class of transcriptional repressors with important regulatory roles in plant growth and development. Through genome-wide analysis of the OVATE gene family of tea plants, 26 and 13 family members were identified in cultivated and ancient tea plants, respectively. Syntenic results showed that OVATE gene family in cultivated tea plants may have experienced a special expansion event. Based on phylogenetic tree analysis, all OVATE genes were divided into four groups, and the third group had the largest number, reaching 16. Transcriptome data from different organs and populations indicated that many OVATE family members were highly expressed in young shoots and leaves, and their expression levels gradually decreased as tea leaves developed. Finally, the expression trends of the six key candidate genes were verified by RT-qPCR, which were consistent with the transcriptome results, indicating that the ovate gene family plays an important role in regulating the process of tea leaf development. In addition, we identified a key structural variation with a length of 184 bp, and the population genotyping showed that it was closely related to the area of tea leaves. Our research provides an important clue for further exploring the function of ovate gene family in tea plants and the development mechanism of tea leaves.

Transposons and non-coding regions drive the intrafamily differences of genome size in insects

Genome size (GS) can vary considerably between phylogenetically close species, but the landscape of GS changes in insects remain largely unclear. To better understand the specific evolutionary factors that determine GS in insects, we examined flow cytometry-based published GS data from 1,326 insect species, spanning 700 genera, 155 families, and 21 orders. Model fitting showed that GS generally followed an Ornstein-Uhlenbeck adaptive evolutionary model in Insecta overall. Ancestral reconstruction indicated a likely GS of 1,069 Mb, suggesting that most insect clades appeared to undergo massive genome expansions or contractions. Quantification of genomic components in 56 species from nine families in four insect orders revealed that the proliferation of transposable elements contributed to high variation in GS between close species, such as within Coleoptera. This study sheds lights on the pattern of GS variation in insects and provides a better understanding of insect GS evolution.

Whole-genome resequencing reveals the origin of tea in Lincang

Phylogeographic, population genetics and diversity analysis are crucial for local tea resource conservation and breeding programs. Lincang in Yunnan has been known as the possible place of domestication for tea worldwide, yet, its genetic makeup and unique Lincang origin are little understood. Here, we reported a large-scale whole-genome resequencing based population genomic analysis in eight main tea-producing areas of Lincang in Yunnan (1,350 accessions), and the first comprehensive map of tea genome variation in Lincang was constructed. Based on the population structure, tea sample in Lincang was divided into three subgroups, and inferred Xigui and Nahan Tea Mountain in Linxiang, Baiying Mountain Ancient Tea Garden in Yun, and Jinxiu Village of Xiaowan Town in Fengqing, which belong to the birthplace of the three subgroups, were all likely to be the origin center of Lincang tea. Meanwhile, the history population sizes analysis show that similar evolutionary patterns were observed for the three subgroups of Lincang. It also was observed that the hybrid among eight areas of Lincang was noticeable, resulting in insignificant genetic differentiation between geographical populations and low genetic diversity. The findings of this study clarified the genetic make-up and evolutionary traits of the local population of tea, which gave some insight into the development of Lincang tea.

The genome and gene editing system of sea barleygrass provide a novel platform for cereal domestication and stress tolerance studies

The tribe Triticeae provides important staple cereal crops and contains elite wild species with wide genetic diversity and high tolerance to abiotic stresses. Sea barleygrass (Hordeum marinum Huds.), a wild Triticeae species, thrives in saline marshlands and is well known for its high tolerance to salinity and waterlogging. Here, a 3.82-Gb high-quality reference genome of sea barleygrass is assembled de novo, with 3.69 Gb (96.8%) of its sequences anchored onto seven chromosomes. In total, 41 045 high-confidence (HC) genes are annotated by homology, de novo prediction, and transcriptome analysis. Phylogenetics, non-synonymous/synonymous mutation ratios (Ka/Ks), and transcriptomic and functional analyses provide genetic evidence for the divergence in morphology and salt tolerance among sea barleygrass, barley, and wheat. The large variation in post-domestication genes (e.g. IPA1 and MOC1) may cause interspecies differences in plant morphology. The extremely high salt tolerance of sea barleygrass is mainly attributed to low Na⁺ uptake and root-to-shoot translocation, which are mainly controlled by SOS1, HKT, and NHX transporters. Agrobacterium-mediated transformation and CRISPR/Cas9-mediated gene editing systems were developed for sea barleygrass to promote its utilization for exploration and functional studies of hub genes and for the genetic improvement of cereal crops.

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.

Genome-wide identification of PME gene family and expression of candidate genes associated with aluminum tolerance in tea plant (Camellia sinensis)

BACKGROUND

The major aluminum (Al) detoxication mechanism of tea plant (Camellia sinensis), as an Al hyperaccumulator plant, is the fixation of almost 70% of Al in the cell walls. Pectin is the primary constituent of cell walls, a degree of methylation of pectin polysaccharides regulated by the pectin methylesterase (PME) genes can greatly affect the Al binding capacity. The knowledge on PME gene family in tea plant is still poor.

RESULTS

We identified 66 (CsPME1-CsPME66) PME genes from C. sinensis genome. We studied their protein characterization, conserved motifs, gene structure, systematic evolution and gene expression under Al treatments, to establish a basis for in-depth research on the function of PMEs in tea plant. Gene structures analysis revealed that the majority of PME genes had 2-4 exons. Phylogenetic results pointed out that the PME genes from the same species displayed comparatively high sequence consistency and genetic similarity. Selective pressure investigation suggested that the Ka/Ks value for homologous genes of PME family was less than one. The expression of CsPMEs under three Al concentration treatments was tissue specific, eight PME genes in leaves and 15 in roots displayed a trend similar to of the Al contents and PME activities under Al concentration treatments, indicating that the degree of pectin de-esterification regulated by PME was crucial for Al tolerance of tea plant.

CONCLUSIONS

Sixty-six CsPME genes were identified for the first time in tea plant. The genome-wide identification, classification, evolutionary and transcription analyses of the PME gene family provided a new direction for further research on the function of PME gene in Al tolerance of tea plant.

Genome-Wide Investigation of the MiR166 Family Provides New Insights into Its Involvement in the Drought Stress Responses of Tea Plants (Camellia sinensis (L.) O. Kuntze)

MicroRNA166 (miR166) is a highly conserved plant miRNA that plays a crucial role in plant growth and the resistance to various abiotic stresses. However, the miR166s in tea (Camellia sinensis (L.) O. Kuntze) have not been comprehensively identified and analyzed. This study identified 30 mature miR166s and twelve pre-miR166s in tea plants. An evolutionary analysis revealed that csn-miR166s originating from the 3′ arm of their precursors were more conserved than the csn-miR166s derived from the 5′ arm of their precursors. The twelve pre-miR166s in tea were divided into two groups, with csn-MIR166 Scaffold364-2 separated from the other precursors. The Mfold-based predictions indicated that the twelve csn-MIR166s formed typical and stable structures comprising a stem-loop hairpin, with minimum free energy ranging from −110.90 to −71.80 kcal/mol. An analysis of the CsMIR166 promoters detected diverse cis-acting elements, including those related to light responses, biosynthesis and metabolism, abiotic stress defenses, and hormone responses. There was no one-to-one relationship between the csn-miR166s and their targets, but most csn-miR166s targeted HD-Zip III genes. Physiological characterization of tea plants under drought stress showed that leaf water content proportionally decreased with the aggravation of drought stress. In contrast, tea leaves’ malondialdehyde (MDA) content proportionally increased. Moreover, the cleavage site of the ATHB-15-like transcript was identified according to a modified 5′ RNA ligase-mediated rapid amplification of cDNA ends. The RT-qPCR data indicated that the transcription of nine csn-miR166s was negatively correlated with their target gene.

The Chromosome-Scale Reference Genome of Macadamia tetraphylla Provides Insights Into Fatty Acid Biosynthesis

Macadamia is an evergreen tree belonging to the Proteaceae family. The two commercial macadamia species, Macadamia integrifolia and M. tetraphylla, are highly prized for their edible kernels. The M. integrifolia genome was recently sequenced, but the genome of M. tetraphylla has to date not been published, which limits the study of biological research and breeding in this species. This study reports a high-quality genome sequence of M. tetraphylla based on the Oxford Nanopore Technologies technology and high-throughput chromosome conformation capture techniques (Hi-C). An assembly of 750.87 Mb with 51.11 Mb N50 length was generated, close to the 740 and 758 Mb size estimates by flow cytometry and k-mer analysis, respectively. Genome annotation indicated that 61.42% of the genome is composed of repetitive sequences and 34.95% is composed of long terminal repeat retrotransposons. Up to 31,571 protein-coding genes were predicted, of which 92.59% were functionally annotated. The average gene length was 6,055 bp. Comparative genome analysis revealed that the gene families associated with defense response, lipid transport, steroid biosynthesis, triglyceride lipase activity, and fatty acid metabolism are expanded in the M. tetraphylla genome. The distribution of fourfold synonymous third-codon transversion showed a recent whole-genome duplication event in M. tetraphylla. Genomic and transcriptomic analysis identified 187 genes encoding 33 crucial oil biosynthesis enzymes, depicting a comprehensive map of macadamia lipid biosynthesis. Besides, the 55 identified WRKY genes exhibited preferential expression in root as compared to that in other tissues. The genome sequence of M. tetraphylla provides novel insights for breeding novel varieties and genetic improvement of agronomic traits.

From Genome Sequencing to CRISPR-Based Genome Editing for Climate-Resilient Forest Trees

Due to the economic and ecological importance of forest trees, modern breeding and genetic manipulation of forest trees have become increasingly prevalent. The CRISPR-based technology provides a versatile, powerful, and widely accepted tool for analyzing gene function and precise genetic modification in virtually any species but remains largely unexplored in forest species. Rapidly accumulating genetic and genomic resources for forest trees enabled the identification of numerous genes and biological processes that are associated with important traits such as wood quality, drought, or pest resistance, facilitating the selection of suitable gene editing targets. Here, we introduce and discuss the latest progress, opportunities, and challenges of genome sequencing and editing for improving forest sustainability.

Combining novel technologies with interdisciplinary basic research to enhance horticultural crops

Horticultural crops mainly include fruits, vegetables, ornamental trees and flowers, and tea trees (Melaleuca alternifolia). They produce a variety of nutrients for the daily human diet in addition to the nutrition provided by staple crops, and some of them additionally possess ornamental and medicinal features. As such, horticultural crops make unique and important contributions to both food security and a colorful lifestyle. Under the current climate change scenario, the growing population and limited arable land means that agriculture, and especially horticulture, has been facing unprecedented challenges to meet the diverse demands of human daily life. Breeding horticultural crops with high quality and adaptability and establishing an effective system that combines cultivation, post-harvest handling, and sales becomes increasingly imperative for horticultural production. This review discusses characteristic and recent research highlights in horticultural crops, focusing on the breeding of quality traits and the mechanisms that underpin them. It additionally addresses challenges and potential solutions in horticultural production and post-harvest practices. Finally, we provide a prospective as to how emerging technologies can be implemented alongside interdisciplinary basic research to enhance our understanding and exploitation of horticultural crops.

Chromosome-Level Genome Assembly Provides New Insights into Genome Evolution and Tuberous Root Formation of Potentilla anserina

Potentilla anserina is a perennial stoloniferous plant with edible tuberous roots in Rosaceae, served as important food and medicine sources for Tibetans in the Qinghai-Tibetan Plateau (QTP), China, over thousands of years. However, a lack of genome information hindered the genetic study. Here, we presented a chromosome-level genome assembly using single-molecule long-read sequencing, and the Hi-C technique. The assembled genome was 454.28 Mb, containing 14 chromosomes, with contig N50 of 2.14 Mb. A total of 46,495 protein-coding genes, 169.74 Mb repeat regions, and 31.76 Kb non-coding RNA were predicted. P. anserina diverged from Potentilla micrantha ∼28.52 million years ago (Mya). Furthermore, P. anserina underwent a recent tetraploidization ∼6.4 Mya. The species-specific genes were enriched in Starch and sucrose metabolism and Galactose metabolism pathways. We identified the sub-genome structures of P. anserina, with A sub-genome was larger than B sub-genome and closer to P. micrantha phylogenetically. Despite lacking significant genome-wide expression dominance, the A sub-genome had higher homoeologous gene expression in shoot apical meristem, flower and tuberous root. The resistance genes was contracted in P. anserina genome. Key genes involved in starch biosynthesis were expanded and highly expressed in tuberous roots, which probably drives the tuber formation. The genomics and transcriptomics data generated in this study advance our understanding of the genomic landscape of P. anserina, and will accelerate genetic studies and breeding programs.

Identification and expression analysis of PUB genes in tea plant exposed to anthracnose pathogen and drought stresses

The plant U-box (PUB) gene family, one of the major ubiquitin ligase families in plants, plays important roles in multiple cellular processes including environmental stress responses and resistance. The function of U-box genes has been well characterized in Arabidopsis and other plants. However, little is known about the tea plant (Camellia sinensis) PUB genes. Here, 89 U-box proteins were identified from the chromosome-scale referenced genome of tea plant. According to the domain organization and phylogenetic analysis, the tea plant PUB family were classified into ten classes, named Class I to X, respectively. Using previously released stress-related RNA-seq data in tea plant, we identified 34 stress-inducible CsPUB genes. Specifically, eight CsPUB genes were expressed differentially under both anthracnose pathogen and drought stresses. Moreover, six of the eight CsPUBs were upregulated in response to these two stresses. Expression profiling performed by qRT-PCR was consistent with the RNA-seq analysis, and stress-related cis-acting elements were identified in the promoter regions of the six upregulated CsPUB genes. These results strongly implied the putative functions of U-box ligase genes in response to biotic and abiotic stresses in tea plant.

Gapless indica rice genome reveals synergistic contributions of active transposable elements and segmental duplications to rice genome evolution

The ultimate goal of genome assembly is a high-accuracy gapless genome. Here, we report a new assembly pipeline that is used to produce a gapless genome for the indica rice cultivar Minghui 63. The resulting 397.71-Mb final assembly is composed of 12 contigs with a contig N50 size of 31.93 Mb. Each chromosome is represented by a single contig and the genomic sequences of all chromosomes are gapless. Quality evaluation of this gapless genome assembly showed that gene regions in our assembly have the highest completeness compared with the other 15 reported high-quality rice genomes. Further comparison with the japonica rice genome revealed that the gapless indica genome assembly contains more transposable elements (TEs) and segmental duplications (SDs), the latter of which produce many duplicated genes that can affect agronomic traits through dose effect or sub-/neo-functionalization. The insertion of TEs can also affect the expression of duplicated genes, which may drive the evolution of these genes. Furthermore, we found the expansion of nucleotide-binding site with leucine-rich repeat disease-resistance genes and cis-zeatin-O-glucosyltransferase growth-related genes in SDs in the gapless indica genome assembly, suggesting that SDs contribute to the adaptive evolution of rice disease resistance and developmental processes. Collectively, our findings suggest that active TEs and SDs synergistically contribute to rice genome evolution.