Genome-level diversification of eight ancient tea populations in the Guizhou and Yunnan regions identifies candidate genes for core agronomic traits

Analysis of population structure and genetic diversity of Camellia tachangensis in Guizhou based on SNP markers

BACKGROUND

Camellia tachangensis F. C. Zhang is a five-compartment species in the ovary of tea group plants, which represents the original germline of early differentiation of some tea group plants.

METHODS AND RESULTS

In this study, we analyzed single-nucleotide polymorphisms (SNPs) at the genome level, constructed a phylogenetic tree, analyzed the genetic diversity, and further investigated the population structure of 100 C. tachangensis accessions using the genotyping-by-sequencing (GBS) method. A total of 91,959 high-quality SNPs were obtained. Population structure analysis showed that the 100 C. tachangensis accessions clustered into three groups: YQ-1 (Village Group), YQ-2 (Forest Group) and YQ-3 (Transition Group), which was further consistent with the results of phylogenetic analysis and principal component analyses (PCA). In addition, a comparative analysis of the genetic diversity among the three populations (Forest, Village, and Transition Groups) detected the highest genetic diversity in the Transition Group and the highest differentiation between Forest and Village Groups.

CONCLUSIONS

C. tachangensis plants growing in the forest had different genetic backgrounds from those growing in villages. This study provides a basis for the effective protection and utilization of C. tachangensis populations and lays a foundation for future C. tachangensis breeding.

CsPAT1, a GRAS transcription factor, promotes lignin accumulation by antagonistic interacting with CsWRKY13 in tea plants

Lignin is an important component of plant cell walls and plays crucial roles in the essential agronomic traits of tea quality and tenderness. However, the molecular mechanisms underlying the regulation of lignin biosynthesis in tea plants remain unclear. CsWRKY13 acts as a negative regulator of lignin biosynthesis in tea plants. In this study, we identified a GRAS transcription factor, phytochrome A signal transduction 1 (CsPAT1), that interacts with CsWRKY13. Silencing CsPAT1 expression in tea plants and heterologous overexpression in Arabidopsis demonstrated that CsPAT1 positively regulates lignin accumulation. Further investigation revealed that CsWRKY13 directly binds to the promoters of CsPAL and CsC4H and suppresses transcription of CsPAL and CsC4H. CsPAT1 indirectly affects the promoter activities of CsPAL and CsC4H by interacting with CsWRKY13, thereby facilitating lignin biosynthesis in tea plants. Compared with the expression of CsWRKY13 alone, the co-expression of CsPAT1 and CsWRKY13 in Oryza sativa significantly increased lignin biosynthesis. Conversely, compared with the expression of CsPAT1 alone, the co-expression of CsPAT1 and CsWRKY13 in O. sativa significantly reduced lignin accumulation. These results demonstrated the antagonistic regulation of the lignin biosynthesis pathway by CsPAT1 and CsWRKY13. These findings improve our understanding of lignin biosynthesis mechanisms in tea plants and provide insights into the role of the GRAS transcription factor family in lignin accumulation.

CsEXL3 regulate mechanical harvest-related droopy leaves under the transcriptional activation of CsBES1.2 in tea plant

Due to a labor shortage, the mechanical harvesting of tea plantations has become a focal point. However, mechanical harvest efficiency was hampered by droopy leaves, leading to a high rate of broken tea shoots and leaves. Here, we dissected the genetic structure of leaf droopiness in tea plants using genome-wide association studies (GWAS) on 146 accessions, combined with transcriptome from two accessions with contrasting droopy leaf phenotypes. A set of 16 quantitative trait loci (QTLs) containing 54 SNPs and 34 corresponding candidate genes associated with droopiness were then identified. Among these, CsEXL3 (EXORDIUM-LIKE 3) from Chromosome 1 emerged as a candidate gene. Further investigations revealed that silencing CsEXL3 in tea plants resulted in weaker vascular cell malformation and brassinosteroid-induced leaf droopiness. Additionally, brassinosteroid signal factor CsBES1.2 was proved to participate in CsEXL3-induced droopiness and vascular cell malformation via using the CsBES1.2-silencing tea plant. Notably, CsBES1.2 bound on the E-box of CsEXL3 promoter to transcriptionally activate CsEXL3 expression as CUT&TAG based ChIP-qPCR and ChIP-seq suggested in vivo as well as EMSA and Y1H indicated in vitro. Furthermore, CsEXL3 instead of CsBES1.2 decreased lignin content and the expressing levels of lignin biosynthesis genes. Overall, our findings suggest that CsEXL3 regulates droopy leaves, partially through the transcriptional activation of CsBES1.2, with the potential to improve mechanical harvest efficiency in tea plantations.

Genetic analyses of ancient tea trees provide insights into the breeding history and dissemination of Chinese Assam tea (Camellia sinensis var. assamica)

Chinese Assam tea (Camellia sinensis var. assamica) is an important tea crop with a long history of cultivation in Yunnan, China. Despite its potential value as a genetic resource, its genetic diversity and domestication/breeding history remain unclear. To address this issue, we genotyped 469 ancient tea plant trees representing 26 C. sinensis var. assamica populations, plus two of its wild relatives (six and three populations of C. taliensis and C. crassicolumna, respectively) using 16 nuclear microsatellite loci. Results showed that Chinese Assam tea has a relatively high, but comparatively lower gene diversity (HS = 0.638) than the wild relative C. crassicolumna (HS = 0.658). Clustering in STRUCTURE indicated that Chinese Assam tea and its two wild relatives formed distinct genetic groups, with considerable interspecific introgression. The Chinese Assam tea accessions clustered into three gene pools, corresponding well with their geographic distribution. However, NewHybrids analysis indicated that 68.48% of ancient Chinese Assam tea plants from Xishuangbanna were genetic intermediates between the Puer and Lincang gene pools. In addition, 10% of the ancient Chinese Assam tea individuals were found to be hybrids between Chinese Assam tea and C. taliensis. Our results suggest that Chinese Assam tea was domesticated separately in three gene pools (Puer, Lincang and Xishuangbanna) in the Mekong River valley and that the hybrids were subsequently selected during the domestication process. Although the domestication history of Chinese Assam tea in southwestern Yunnan remains complex, our results will help to identify valuable genetic resources that may be useful in future tea breeding programs.

Genome-Wide Association Studies of Biluochun Tea Plant Populations in Dongting Mountain and Comprehensive Identification of Candidate Genes Associated with Core Agronomic Traits by Four Analysis Models

The elite germplasm resources are key to the beautiful appearance and pleasant flavor of Biluochun tea. We collected and measured the agronomic traits of 95 tea plants to reveal the trait diversity and breeding value of Biluochun tea plant populations. The results revealed that the agronomic traits of Biluochun tea plant populations were diverse and had high breeding value. Additionally, we resequenced these tea plant populations to reveal genetic diversity, population structure, and selection pressure. The Biluochun tea plant populations contained two groups and were least affected by natural selection based on the results of population structure and selection pressure. More importantly, four non-synonymous single nucleotide polymorphisms (nsSNPs) and candidate genes associated with (-)-gallocatechin gallate (GCG), (-)-gallocatechin (GC), and caffeine (CAF) were detected using at least two GWAS models. The results will promote the development and application of molecular markers and the utilization of elite germplasm from Biluochun populations.

The transcription factor CsS40 negatively regulates TCS1 expression and caffeine biosynthesis in connection to leaf senescence in Camellia sinensis

Caffeine is considered as one of the most important bioactive components in the popular plant beverages tea, cacao, and coffee, but as a wide-spread plant secondary metabolite its biosynthetic regulation at transcription level remains largely unclear. Here, we report a novel transcription factor Camellia sinensis Senescnece 40 (CsS40) as a caffeine biosynthesis regulator, which was discovered during screening a yeast expression library constructed from tea leaf cDNAs for activation of tea caffeine synthase (TCS1) promoter. Besides multiple hits of the non-self-activation CsS40 clones that bound to and activated TCS1 promoter in yeast-one-hybrid assays, a split-luciferase complementation assay demonstrated that CsS40 acts as a transcription factor to activate the CsTCS1 gene and EMSA assay also demonstrated that CsS40 bound to the TCS1 gene promoter. Consistently, immunofluorescence data indicated that CsS40-GFP fusion was localized in the nuclei of tobacco epidermal cells. The expression pattern of CsS40 in 'Fuding Dabai' developing leaves was opposite to that of TCS1; and knockdown and overexpression of CsS40 in tea leaf calli significantly increased and decreased TCS1 expression levels, respectively. The expression levels of CsS40 were also negatively correlated to caffeine accumulation in developing leaves and transgenic calli of 'Fuding Dabai'. Furthermore, overexpression of CsS40 reduced the accumulation of xanthine and hypoxanthine in tobacco plants, meanwhile, increased their susceptibility to aging. CsS40 expression in tea leaves was also induced by senescence-promoting hormones and environmental factors. Taken together, we showed that a novel senescence-related factor CsS40 negatively regulates TCS1 and represses caffeine accumulation in tea cultivar 'Fuding Dabai'. The study provides new insights into caffeine biosynthesis regulation by a plant-specific senescence regulator in tea plants in connection to leaf senescence and hormone signaling.

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.

Molecular markers in tea plant (Camellia sinensis): Applications to evolution, genetic identification, and molecular breeding

Tea plants have a long cultivation history in the world, and the beverage (tea) made from its leaves is well known in the world. Due to the characteristics of self-incompatibility, long-term natural and artificial hybridization, tea plants have a very complex genetic background, which make the classification of tea plants unclear. Molecular marker, one type of genetic markers, has the advantages of stable inheritance, large amount of information, and high reliability. The development of molecular marker has facilitated the understanding of complex tea germplasm resources. So far, molecular markers had played important roles in the study of the origin and evolution, the preservation and identification of tea germplasms, and the excellent cultivars breeding of tea plants. However, the information is scattered, making it difficult to understand the advance of molecular markers in tea plants. In this paper, we summarized the development process and types of molecular markers in tea plants. In addition, the application advance of these molecular markers in tea plants was reviewed. Perspectives of molecular markers in tea plants were also systematically provided and discussed. The elaboration of molecular markers in this paper should help us to renew understanding of its application in tea plants.

Genome-wide association study of leaf-related traits in tea plant in Guizhou based on genotyping-by-sequencing

BACKGROUND

Studying the genetic characteristics of tea plant (Camellia spp.) leaf traits is essential for improving yield and quality through breeding and selection. Guizhou Plateau, an important part of the original center of tea plants, has rich genetic resources. However, few studies have explored the associations between tea plant leaf traits and single nucleotide polymorphism (SNP) markers in Guizhou.

RESULTS

In this study, we used the genotyping-by-sequencing (GBS) method to identify 100,829 SNP markers from 338 accessions of tea germplasm in Guizhou Plateau, a region with rich genetic resources. We assessed population structure based on high-quality SNPs, constructed phylogenetic relationships, and performed genome-wide association studies (GWASs). Four inferred pure groups (G-I, G-II, G-III, and G-IV) and one inferred admixture group (G-V), were identified by a population structure analysis, and verified by principal component analyses and phylogenetic analyses. Through GWAS, we identified six candidate genes associated with four leaf traits, including mature leaf size, texture, color and shape. Specifically, two candidate genes, located on chromosomes 1 and 9, were significantly associated with mature leaf size, while two genes, located on chromosomes 8 and 11, were significantly associated with mature leaf texture. Additionally, two candidate genes, located on chromosomes 1 and 2 were identified as being associated with mature leaf color and mature leaf shape, respectively. We verified the expression level of two candidate genes was verified using reverse transcription quantitative polymerase chain reaction (RT-qPCR) and designed a derived cleaved amplified polymorphism (dCAPS) marker that co-segregated with mature leaf size, which could be used for marker-assisted selection (MAS) breeding in Camellia sinensis.

CONCLUSIONS

In the present study, by using GWAS approaches with the 338 tea accessions population in Guizhou, we revealed a list of SNPs markers and candidate genes that were significantly associated with four leaf traits. This work provides theoretical and practical basis for the genetic breeding of related traits in tea plant leaves.

Geochemical characteristics of heavy metals of bedrock, soil, and tea in a metamorphic rock area of Guizhou Province, China

Through field surveys, sampling, and laboratory experiments, the content, enrichment, and migration characteristics of heavy metals in the rock-soil-tea system in the Guizhou metamorphic rock area were analyzed herein. The results show that the As and Hg contents in metamorphic sandstone in the study area are higher than those in slate, while the Mn, Zn, Pb, and Ni contents in slate are higher. The content of heavy metals in the sandstone of the Jialu formation of the Xiajiang group in Neoproterozoic is the highest, and that of the slate of the Fanzhao formation is the lowest. The content of As in the rock samples in the study area is higher than that in the Epicontinental crust, while Cd is found to be deficient. The heavy metals at some soil points exceed the standard, but this does not affect the growth and development of tea. The content of heavy metals in most tea is lower than the safety limit standards, and the THQ value is less than 1, indicating that the heavy metals in tea in Guizhou metamorphic rock areas pose no concern for human health. In the metamorphic rock distribution areas of Guizhou, tea does not enrich As, Cr, Hg, and Pb in soil, but enriches Mn, Cu, Ni, and Zn to varying degrees, with the enrichment of Mn being the strongest.

Topics and trends in fresh tea (Camellia sinensis) leaf research: A comprehensive bibliometric study

Tea plant (Camellia sinensis) is a widely cultivated cash crop and tea is a favorite functional food in the world. Fresh tea leaves (FTLs) play a critical role in bridging the two fields closely related to tea cultivation and tea processing, those are, tea plant biology and tea biochemistry. To provide a comprehensive overview of the development stages, authorship collaboration, research topics, and hotspots and their temporal evolution trends in the field of FTLs research, we conducted a bibliometric analysis, based on 971 publications on FTLs-related research published during 2001-2021 from Web of Science Core Collection. CiteSpace, R package Bibliometrix, and VOSviewer were employed in this research. The results revealed that the development history can be roughly divided into three stages, namely initial stage, slow development stage and rapid development stage. Journal of Agricultural & Food Chemistry published most articles in this field, while Frontiers in Plant Science held the highest total citations and h-index. The most influential country, institution, and author in this field was identified as China, the Chinese Academy of Agricultural Sciences, and Xiaochun Wan, respectively. FTLs-related research can be categorized into three main topics: the regulation mechanism of key genes, the metabolism and features of essential compounds, and tea plants' growth and stress responses. The most concerning hotspots are the application of advanced technologies, essential metabolites, leaf color variants, and effective cultivation treatments. There has been a shift from basic biochemical and enzymatic studies to studies of molecular mechanisms that depend on multi-omics technologies. We also discussed the future development in this field. This study provides a comprehensive summary of the research field, making it easier for researchers to be informed about its development history, status, and trends.

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.

Evaluating the productivity of ancient Pu'er tea trees (Camellia sinensis var. assamica): a multivariate modeling approach

BACKGROUND

The demand for productive economic plant resources is increasing with the continued growth of the human population. Ancient Pu'er tea trees [Camellia sinensis var. assamica (J. W. Mast.) Kitam.] are an important ecological resource with high economic value and large interests. The study intends to explore and evaluate critical drivers affecting the species' productivity, then builds formulas and indexes to make predicting the productivity of such valuable plant resources possible and applicable.

RESULTS

Our analysis identified the ideal values of the seven most important environmental variables and their relative contribution (shown in parentheses) to the distribution of ancient Pu'er tea trees: annual precipitation, ca. 1245 mm (28.73%); min temperature of coldest month, ca. 4.2 °C (18.25%); precipitation of driest quarter, ca. 47.5 mm (14.45%); isothermality, 49.9% to 50.4% (14.11%); precipitation seasonality, ca. 89.2 (6.77%); temperature seasonality, ca. 391 (4.46%); and solar radiation, 12,250 to 13,250 kJ m^-2 day^-1 (3.28%). Productivity was indicated by the total value (viz. fresh leaf harvested multiplied by unit price) of each tree. Environmental suitability, tree growth, and management positively affected productivity; regression weights were 0.325, 0.982, and 0.075, respectively. The degree of productivity was classified as follows: > 0.8, "highly productive"; 0.5-0.8, "productive"; 0.3-0.5, "poorly productive"; and < 0.3, "unproductive". Overall, 53% of the samples were categorized as "poorly productive" or "unproductive"; thus, the management of these regions require attention.

CONCLUSIONS

This model improves the accuracy of the predictions of ancient Pu'er tea tree productivity and will aid future analyses of distribution shifts under climate change, as well as the identification of areas suitable for Pu'er tea tree plantations. Our modeling framework provides insights that facilitate the interpretation of abstract concepts and could be applied to other economically valuable plant resources.

Genetic analyses reveal regional structure and demographic expansion of the predominant tea pest Empoasca onukii (Hemiptera: Cicadellidae) in China

BACKGROUND

The tea green leafhopper, Empoasca onukii Matsuda, is the most destructive insect pest of tea plantations in East Asia. Despite its economic importance and previous studies on this species, it remains unclear as to how this small-sized pest can have such wide range.

RESULTS

By sequencing three mitochondrial genes and 17 microsatellite loci, we revealed the regional structure and demographic expansion of 59 E. onukii populations in China. Bayesian analysis of population genetic structure (BAPS) on microsatellites identified four genetic groups with spatial discontinuities, while analysis on mitochondrial genes inferred five nested and differentiated clusters. Both the Mantel test and the generalized linear model indicated a significant pattern of isolation by geographic distance in E. onukii populations. Based on the approximate Bayesian computation approach, E. onukii was found to have originated from southwestern China and expanded northward and eastward. While MIGRATE-N and Bayesian stochastic search variable selection (BSSVS) procedure in BEAST confirmed the possible eastward and northward dispersal from Yunnan, they also detected more gene flow from the derived populations in central and southeastern China.

CONCLUSION

Our results suggest that the current distribution and structure of E. onukii is complicatedly influenced by human activities of cultivation, wide dissemination of tea in ancient China as well as recent transportation of tea seedlings for establishing new tea plantations. Insights into genetic differentiation and demographic expansion patterns from this study provide an important basis for the development of area-wide management of the E. onukii populations. © 2022 Society of Chemical Industry.

The Future Potential Distribution and Sustainable Management of Ancient Pu’er Tea Trees (Camellia sinensis var. assamica (J. W. Mast.) Kitam.)

Ancient Pu’er tea trees (Camellia sinensis var. assamica (J. W. Mast.) Kitam.) are an important ecological resource with high economic value. Knowledge of the environmental variables shaping the original distribution and the effects of climate change on the future potential distribution of these trees, as well as the identification of sustainable management approaches, is essential for ensuring their future health and production. Here, we used 28 current environmental variables and the future climate data to model the suitable areas for ancient Pu’er tea trees. We also compared the health of these ancient trees in areas under different local management strategies. The results suggested the general distribution is likely to remain stable, but there are environmentally suitable areas outside its current habitats. To achieve more sustainable management, the main areas in which the management of poorly-managed trees can be improved include learning from managers of well-managed trees and following the common technical management regulations stipulated by the local government. The suitable value ranges for environmental factors, potentially suitable areas under climate change, and assessment of management approaches will aid the future cultivation and transplantation of ancient Pu’er tea trees. The methodology includes management-level analysis and provides practical insights that could be applied to regions outside the most suitable areas identified.

The reference genome of camellia chekiangoleosa provides insights into camellia evolution and tea oil biosynthesis

Camellia oil extracted from Camellia seeds is rich in unsaturated fatty acids (UFAs) and secondary metabolites beneficial to human health. However, no oil-tea tree genome has yet been published, which is a major obstacle to investigating the heredity improvement of oil-tea trees. Here, using both Illumina and PicBio sequencing technologies, we present the first chromosome-level genome sequence of the oil-tea tree species Camellia chekiangoleosa Hu. (CCH). The assembled genome consists of 15 pseudochromosomes with a genome size of 2.73 Gb and a scaffold N50 of 185.30 Mb. At least 2.16 Gb of the genome assembly consists of repetitive sequences, and the rest involves a high-confidence set of 64 608 protein-coding gene models. Comparative genomic analysis revealed that the CCH genome underwent a whole-genome duplication (WGD) event shared across the Camellia genus at ~57.48 MYA and a γ-WGT event shared across all core eudicot plants at ~120 MYA. Gene family clustering revealed that the genes involved in terpenoid biosynthesis have undergone rapid expansion. Furthermore, we determined the expression patterns of oleic acid accumulation- and terpenoid biosynthesis-associated genes in six tissues. We found that these genes tend to be highly expressed in leaves, pericarp tissues, roots, and seeds. The first chromosome-level genome of oil-tea trees will provide valuable resources for determining Camellia evolution and utilizing the germplasm of this taxon.