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

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.

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.

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.

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.

•

The most comprehensive variation pattern of insects genome size (GS) to date

•

GS evolution of the Insecta was reflected by the adaptive Ornstein–Uhlenbeck model

•

Ancestral state of insect GS was estimated to be ∼1 Gb

•

Intrafamily GS variations were driven by recent transpositions and non-coding regions

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.

An ancient whole-genome duplication event and its contribution to flavor compounds in the tea plant (Camellia sinensis)

Tea, coffee, and cocoa are the three most popular nonalcoholic beverages in the world and have extremely high economic and cultural value. The genomes of four tea plant varieties have recently been sequenced, but there is some debate regarding the characterization of a whole-genome duplication (WGD) event in tea plants. Whether the WGD in the tea plant is shared with other plants in order Ericales and how it contributed to tea plant evolution remained unanswered. Here we re-analyzed the tea plant genome and provided evidence that tea experienced only WGD event after the core-eudicot whole-genome triplication (WGT) event. This WGD was shared by the Polemonioids-Primuloids-Core Ericales (PPC) sections, encompassing at least 17 families in the order Ericales. In addition, our study identified eight pairs of duplicated genes in the catechins biosynthesis pathway, four pairs of duplicated genes in the theanine biosynthesis pathway, and one pair of genes in the caffeine biosynthesis pathway, which were expanded and retained following this WGD. Nearly all these gene pairs were expressed in tea plants, implying the contribution of the WGD. This study shows that in addition to the role of the recent tandem gene duplication in the accumulation of tea flavor-related genes, the WGD may have been another main factor driving the evolution of tea flavor.

A consensus and saturated genetic map provides insight into genome anchoring, synteny of Solanaceae and leaf- and fruit-related QTLs in wolfberry (Lycium Linn.)

BACKGROUND

Lycium Linn. (Solanaceae) is a genus of economically important plants producing fruits and leaves with high nutritional value and medicinal benefits. However, genetic analysis of this plant and molecular breeding for quality improvement are limited by the lack of sufficient molecular markers.

RESULTS

In this study, two parental strains, 'Ningqi No. 1' (Lycium barbarum L.) and 'Yunnan Gouqi' (Lycium yunnanense Kuang et A.M. Lu), and 200 F₁ hybrid individuals were resequenced for genetic analysis. In total, 8,507 well-selected SNPs were developed, and a high-density genetic map (NY map) was constructed with a total genetic distance of 2,122.24 cM. A consensus genetic map was established by integrating the NY map and a previously published genetic map (NC map) containing 15,240 SNPs, with a total genetic distance of 3,058.19 cM and an average map distance of 0.21 cM. The 12 pseudochromosomes of the Lycium reference genome were anchored using this consensus genetic map, with an anchoring rate of 64.3%. Moreover, weak collinearities between the consensus map and the pepper, potato, and tomato genomes were observed. Twenty-five stable QTLs were identified for leaf- and fruit-related phenotypes, including fruit weight, fruit longitude, leaf length, the fruit index, and the leaf index; these stable QTLs were mapped to four different linkage groups, with LOD scores ranging from 2.51 to 19.37 and amounts of phenotypic variance explained from 6.2% to 51.9%. Finally, 82 out of 188 predicted genes underlying stable QTLs for fruit-related traits were differentially expressed according to RNA-seq analysis.

CONCLUSIONS

A chromosome-level assembly can provide a foundation for further functional genomics research for wolfberry. The genomic regions of these stably expressed QTLs could be used as targets for further fine mapping and development of molecular markers for marker-assisted selection (MAS). The present study provided valuable information on saturated SNP markers and reliable QTLs for map-based cloning of functional genes related to yield and morphological traits in Lycium spp.

The Welwitschia genome reveals a unique biology underpinning extreme longevity in deserts

The gymnosperm Welwitschia mirabilis belongs to the ancient, enigmatic gnetophyte lineage. It is a unique desert plant with extreme longevity and two ever-elongating leaves. We present a chromosome-level assembly of its genome (6.8 Gb/1 C) together with methylome and transcriptome data to explore its astonishing biology. We also present a refined, high-quality assembly of Gnetum montanum to enhance our understanding of gnetophyte genome evolution. The Welwitschia genome has been shaped by a lineage-specific ancient, whole genome duplication (~86 million years ago) and more recently (1-2 million years) by bursts of retrotransposon activity. High levels of cytosine methylation (particularly at CHH motifs) are associated with retrotransposons, whilst long-term deamination has resulted in an exceptionally GC-poor genome. Changes in copy number and/or expression of gene families and transcription factors (e.g. R2R3MYB, SAUR) controlling cell growth, differentiation and metabolism underpin the plant's longevity and tolerance to temperature, nutrient and water stress.

Comparative Multi-Omics of Tender Shoots from a Novel Evergrowing Tea Cultivar Provide Insight into the Winter Adaptation Mechanism

Tea (Camellia sinensis [L.] O. Kuntze) tree is a perennial plant in which winter dormancy is an important biological adaptation to environmental changes. We discovered and reported a novel tea tree cultivar that can generate tender shoots in winter several years ago, but the molecular mechanism for this unique phenotype remains unknown . Here, we conducted comparative transcriptomics, proteomics and metabolomics along with phytohormone quantitation between the winter and spring tender shoots to investigate the physiological basis and putative regulatory mechanisms of its evergrowing character during winter. Our multi-omics study has led to the following findings. Gibberellin (GA) levels and key enzymes for GA biosynthesis and the signal transduction pathway were increased in the winter shoots, causing the ABA/GA content ratio to decrease, which might play a key regulatory role in maintaining normal growth during winter. The abundance of proteins, genes and metabolites involved in energy metabolism was all increased in winter shoots, indicating that energy is critical for continuous growth under the relatively weak-light and low-temperature environment. Abiotic resistance-related proteins and free amino acids were also increased in abundance in the winter shoots, which possibly represents an adaptation response to winter conditions. These results allowed us to hypothesize a novel molecular mechanism of adaptation for this unique tender shoot evergrowing in winter.

Genetic basis of high aroma and stress tolerance in the oolong tea cultivar genome

Tea plants (Camellia sinensis) are commercially cultivated in >60 countries, and their fresh leaves are processed into tea, which is the most widely consumed beverage in the world. Although several chromosome-level tea plant genomes have been published, they collapsed the two haplotypes and ignored a large number of allelic variations that may underlie important biological functions in this species. Here, we present a phased chromosome-scale assembly for an elite oolong tea cultivar, "Huangdan", that is well known for its high levels of aroma. Based on the two sets of haplotype genome data, we identified numerous genetic variations and a substantial proportion of allelic imbalance related to important traits, including aroma- and stress-related alleles. Comparative genomics revealed extensive structural variations as well as expansion of some gene families, such as terpene synthases (TPSs), that likely contribute to the high-aroma characteristics of the backbone parent, underlying the molecular basis for the biosynthesis of aroma-related chemicals in oolong tea. Our results uncovered the genetic basis of special features of this oolong tea cultivar, providing fundamental genomic resources to study evolution and domestication for the economically important tea crop.

Chromatin accessibility and translational landscapes of tea plants under chilling stress

Plants have evolved regulatory mechanisms at multiple levels to regulate gene expression in order to improve their cold adaptability. However, limited information is available regarding the stress response at the chromatin and translational levels. Here, we characterize the chromatin accessibility, transcriptional, and translational landscapes of tea plants in vivo under chilling stress for the first time. Chilling stress significantly affected both the transcription and translation levels as well as the translation efficiency of tea plants. A total of 3010 genes that underwent rapid and independent translation under chilling stress were observed, and they were significantly enriched in the photosynthesis-antenna protein and phenylpropanoid biosynthesis pathways. A set of genes that were significantly responsive to cold at the transcription and translation levels, including four (+)-neomenthol dehydrogenases (MNDs) and two (E)-nerolidol synthases (NESs) arranged in tandem on the chromosomes, were also found. We detected potential upstream open reading frames (uORFs) on 3082 genes and found that tea plants may inhibit the overall expression of genes by enhancing the translation of uORFs under chilling stress. In addition, we identified distal transposase hypersensitive sites (THSs) and proximal THSs and constructed a transcriptional regulatory network for tea plants under chilling stress. We also identified 13 high-confidence transcription factors (TFs) that may play a crucial role in cold regulation. These results provide valuable information regarding the potential transcriptional regulatory network in plants and help to clarify how plants exhibit flexible responses to chilling stress.

Effects of Long-Term Nitrogen Fertilization on the Formation of Metabolites Related to Tea Quality in Subtropical China

As a main agronomic intervention in tea cultivation, nitrogen (N) application is useful to improve tea yield and quality. However, the effects of N application on the formation of tea quality-related metabolites have not been fully studied, especially in long-term field trials. In this study, a 10-year field experiment was conducted to investigate the effect of long-term N application treatments on tea quality-related metabolites, their precursors, and related gene expression. Long-term N application up-regulated the expression of key genes for chlorophyll synthesis and promoted its synthesis, thus increasing tea yield. It also significantly increased the contents of total free amino acids, especially l-theanine, in fresh tea leaves, while decreasing the catechin content, which is conducive to enhancing tea liquor freshness. However, long-term N application significantly reduced the contents of benzyl alcohol and 2-phenylethanol in fresh tea leaves, and also reduced (E)-nerolidol and indole in withered leaves, which were not conducive to the formation of floral and fruity aroma compounds. In general, an appropriate amount of N fertilizer (225 kg/hm²) balanced tea yield and quality. These results not only provide essential information on how N application affects tea quality, but also provide detailed experimental data for field fertilization.

The chromosome-level reference genome assembly for Panax notoginseng and insights into ginsenoside biosynthesis

Panax notoginseng, a perennial herb of the genus Panax in the family Araliaceae, has played an important role in clinical treatment in China for thousands of years because of its extensive pharmacological effects. Here, we report a high-quality reference genome of P. notoginseng, with a genome size up to 2.66 Gb and a contig N50 of 1.12 Mb, produced with third-generation PacBio sequencing technology. This is the first chromosome-level genome assembly for the genus Panax. Through genome evolution analysis, we explored phylogenetic and whole-genome duplication events and examined their impact on saponin biosynthesis. We performed a detailed transcriptional analysis of P. notoginseng and explored gene-level mechanisms that regulate the formation of characteristic tubercles. Next, we studied the biosynthesis and regulation of saponins at temporal and spatial levels. We combined multi-omics data to identify genes that encode key enzymes in the P. notoginseng terpenoid biosynthetic pathway. Finally, we identified five glycosyltransferase genes whose products catalyzed the formation of different ginsenosides in P. notoginseng. The genetic information obtained in this study provides a resource for further exploration of the growth characteristics, cultivation, breeding, and saponin biosynthesis of P. notoginseng.