Bulked Segregant RNA-Seq Reveals Different Gene Expression Patterns and Mutant Genes Associated with the Zigzag Pattern of Tea Plants (Camellia sinensis).

The unique zigzag-patterned tea plant is a rare germplasm resource. However, the molecular mechanism behind the formation of zigzag stems remains unclear. To address this, a BC1 genetic population of tea plants with zigzag stems was studied using histological observation and bulked segregant RNA-seq. The analysis revealed 1494 differentially expressed genes (DEGs) between the upright and zigzag stem groups. These DEGs may regulate the transduction and biosynthesis of plant hormones, and the effects on the phenylpropane biosynthesis pathways may cause the accumulation of lignin. Tissue sections further supported this finding, showing differences in cell wall thickness between upright and curved stems, potentially due to lignin accumulation. Additionally, 262 single-nucleotide polymorphisms (SNPs) across 38 genes were identified as key SNPs, and 5 genes related to zigzag stems were identified through homologous gene function annotation. Mutations in these genes may impact auxin distribution and content, resulting in the asymmetric development of vascular bundles in curved stems. In summary, we identified the key genes associated with the tortuous phenotype by using BSR-seq on a BC1 population to minimize genetic background noise.

Deeply functional identification of TCS1 alleles provides efficient technical paths for low-caffeine breeding of tea plants.

Caffeine is an important functional component in tea, which has the effect of excitement and nerve stimulation, but excessive intake can cause insomnia and dysphoria. Therefore, the production of tea with low-caffeine content can meet the consumption needs of certain people. Here, in addition to the previous alleles of the tea caffeine synthase (TCS1) gene, a new allele (TCS1h) from tea germplasms was identified. Results of in vitro activity analysis showed that TCS1h had both theobromine synthase (TS) and caffeine synthase (CS) activities. Site-directed mutagenesis experiments of TCS1a, TCS1c, and TCS1h demonstrated that apart from the 225th amino acid residue, the 269th amino acid also determined the CS activity. GUS histochemical analysis and dual-luciferase assay indicated the low promoter activity of TCS1e and TCS1f. In parallel, insertion and deletion mutations in large fragments of alleles and experiments of site-directed mutagenesis identified a key cis-acting element (G-box). Furthermore, it was found that the contents of purine alkaloids were related to the expression of corresponding functional genes and alleles, and the absence or presence and level of gene expression determined the content of purine alkaloids in tea plants to a certain extent. In summary, we concluded TCS1 alleles into three types with different functions and proposed a strategy to effectively enhance low-caffeine tea germplasms in breeding practices. This research provided an applicable technical avenue for accelerating the cultivation of specific low-caffeine tea plants.

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.

Transcriptomic and Metabolomic Analyses Provide Insights Into an Aberrant Tissue of Tea Plant (Camellia sinensis).

Tea plant (Camellia sinensis (L.) O. Kuntze) is one of the most important economic crops with multiple mutants. Recently, we found a special tea germplasm that has an aberrant tissue on its branches. To figure out whether this aberrant tissue is associated with floral bud (FB) or dormant bud (DB), we performed tissue section, transcriptome sequencing, and metabolomic analysis of these tissues. Longitudinal sections indicated the aberrant tissue internal structure was more like a special bud (SB), but was similar to that of DB. Transcriptome data analysis showed that the number of heterozygous and homozygous SNPs was significantly different in the aberrant tissue compared with FB and DB. Further, by aligning the unmapped sequences of the aberrant tissue to the Non-Redundant Protein Sequences (NR) database, we observed that 36.13% of unmapped sequences were insect sequences, which suggested that the aberrant tissue might be a variation of dormant bud tissue influenced by the interaction of tea plants and insects or pathogens. Metabolomic analysis showed that the differentially expressed metabolites (DEMs) between the aberrant tissue and DB were significantly enriched in the metabolic pathways of biosynthesis of plant hormones and biosynthesis of phenylpropanoids. Subsequently, we analyzed the differentially expressed genes (DEGs) in the above mentioned two tissues, and the results indicated that photosynthetic capacity in the aberrant tissue was reduced, whereas the ethylene, salicylic acid and jasmonic acid signaling pathways were activated. We speculated that exogenous infection induced programmed cell death (PCD) and increased the lignin content in dormant buds of tea plants, leading to the formation of this aberrant tissue. This study advanced our understanding of the interaction between plants and insects or pathogens, providing important clues about biotic stress factors and key genes that lead to mutations and formation of the aberrant tissue.

Identifying Conserved Functional Gene Modules Underlying the Dynamic Regulation of Tea Plant Development and Secondary Metabolism.

Tea plants adjust development and metabolism by integrating environmental and endogenous signals in complex but poorly defined gene networks. Here, we present an integrative analysis framework for the identification of conserved modules controlling important agronomic traits using a comprehensive collection of RNA-seq datasets in Camellia plants including 189 samples. In total, 212 secondary metabolism-, 182 stress response-, and 182 tissue development-related coexpressed modules were revealed. Functional modules (e.g., drought response, theobromine biosynthesis, and new shoot development-related modules) and potential regulators that were highly conserved across diverse genetic backgrounds and/or environmental conditions were then identified by cross-experiment comparisons and consensus clustering. Moreover, we investigate the preservation of gene networks between Camellia sinensis and other Camellia species. This revealed that the coexpression patterns of several recently evolved modules related to secondary metabolism and environmental adaptation were rewired and showed higher connectivity in tea plants. These conserved modules are excellent candidates for modeling the core mechanism of tea plant development and secondary metabolism and should serve as a great resource for hypothesis generation and tea quality improvement.

The chromosome-scale genome reveals the evolution and diversification after the recent tetraploidization event in tea plant.

Tea is one of the most popular nonalcoholic beverages due to its characteristic secondary metabolites with numerous health benefits. Although two draft genomes of tea plant (Camellia sinensis) have been published recently, the lack of chromosome-scale assembly hampers the understanding of the fundamental genomic architecture of tea plant and potential improvement. Here, we performed a genome-wide chromosome conformation capture technique (Hi-C) to obtain a chromosome-scale assembly based on the draft genome of C. sinensis var. sinensis and successfully ordered 2984.7 Mb (94.7%) scaffolds into 15 chromosomes. The scaffold N50 of the improved genome was 218.1 Mb, ~157-fold higher than that of the draft genome. Collinearity comparison of genome sequences and two genetic maps validated the high contiguity and accuracy of the chromosome-scale assembly. We clarified that only one Camellia recent tetraploidization event (CRT, 58.9-61.7 million years ago (Mya)) occurred after the core-eudicot common hexaploidization event (146.6-152.7 Mya). Meanwhile, 9243 genes (28.6%) occurred in tandem duplication, and most of these expanded after the CRT event. These gene duplicates increased functionally divergent genes that play important roles in tea-specific biosynthesis or stress response. Sixty-four catechin- and caffeine-related quantitative trait loci (QTLs) were anchored to chromosome assembly. Of these, two catechin-related QTL hotspots were derived from the CRT event, which illustrated that polyploidy has played a dramatic role in the diversification of tea germplasms. The availability of a chromosome-scale genome of tea plant holds great promise for the understanding of genome evolution and the discovery of novel genes contributing to agronomically beneficial traits in future breeding programs.

Repressed Gene Expression of Photosynthetic Antenna Proteins Associated with Yellow Leaf Variation as Revealed by Bulked Segregant RNA-seq in Tea Plant Camellia sinensis.

The young leaves and shoots of albino tea cultivars are usually characterized as having a yellow or pale color, high amino acid, and low catechin. Increasing attention has been paid to albino tea cultivars in recent years because their tea generally shows high umami and reduced astringency. However, the genetic mechanism of yellow-leaf variation in albino tea cultivar has not been elucidated clearly. In this study, bulked segregant RNA-seq (BSR-seq) was performed on bulked yellow- and green-leaf hybrid progenies from a leaf color variation population. A total of 359 and 1134 differentially expressed genes (DEGs) were identified in the yellow and green hybrid bulked groups (Yf vs Gf) and parent plants (Yp vs Gp), respectively. The significantly smaller number of DEGs in Yf versus Gf than in Yp versus Gp indicated that individual differences could be reduced within the same hybrid progeny. Analysis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes revealed that the photosynthetic antenna protein was most significantly enriched in either the bulked groups or their parents. Interaction was found among light-harvesting chlorophyll a/b -binding proteins (LHC), heat shock proteins (HSPs), and enzymes involved in cuticle formation. Combined with the transcriptomic expression profile, results showed that the repressed genes encoding LHC were closely linked to aberrant chloroplast development in yellow-leaf tea plants. Furthermore, the photoprotection and light stress response possessed by genes involved in HSP protein interaction and cuticle formation were discussed. The expression profile of DEGs was verified via quantitative real-time PCR analysis of the bulked samples and other F1 individuals. In summary, using BSR-seq on a hybrid population eliminated certain disturbing effects of genetic background and individual discrepancy, thereby helping this study to intensively focus on the key genes controlling leaf color variation in yellow-leaf tea plants.

Identification and distribution of a single nucleotide polymorphism responsible for the catechin content in tea plants.

Catechins are the predominant products in tea plants and have essential functions for both plants and humans. Several genes encoding the enzymes regulating catechin biosynthesis have been identified, and the identification of single nucleotide polymorphisms (SNPs) resulting in nonsynonymous mutations within these genes can be used to establish a functional link to catechin content. Therefore, the transcriptomes of two parents and four filial offspring were sequenced using next-generation sequencing technology and aligned to the reference genome to enable SNP mining. Subsequently, 176 tea plant accessions were genotyped based on candidate SNPs using kompetitive allele-specific polymerase chain reaction (KASP). The catechin contents of these samples were characterized by high-performance liquid chromatography (HPLC), and analysis of variance (ANOVA) was subsequently performed to determine the relationship between genotypes and catechin content. As a result of these efforts, a SNP within the chalcone synthase (CHS) gene was shown to be functionally associated with catechin content. Furthermore, the geographical and interspecific distribution of this SNP was investigated. Collectively, these results will contribute to the early evaluation of tea plants and serve as a rapid tool for accelerating targeted efforts in tea breeding.

Gene Coexpression Networks Reveal Key Drivers of Flavonoid Variation in Eleven Tea Cultivars (Camellia sinensis).

Following the recent completion of the draft genome sequence of the tea plant, high-throughput decoding of gene function, especially for those involved in complex secondary metabolic pathways, has become a major challenge. Here, we profiled the metabolome and transcriptome of 11 tea cultivars, and then illustrated a weighted gene coexpression network analysis (WGCNA)-based system biological strategy to interpret metabolomic flux, predict gene functions, and mine key regulators involved in the flavonoid biosynthesis pathway. We constructed a multilayered regulatory network, which integrated the gene coexpression relationship with the microRNA target and promoter cis-regulatory element information. This allowed us to reveal new uncharacterized TFs (e.g., MADSs, WRKYs, and SBPs) and microRNAs (including 17 conserved and 15 novel microRNAs) that are potentially implicated in different steps of the catechin biosynthesis. Furthermore, we applied metabolic-signature-based association method to capture additional key regulators involved in catechin pathway. This provides important clues for the functional characterization of five SCPL1A acyltransferase family members, which might be implicated in the production balance of anthocyanins, galloylated catechins, and proanthocyanins. Application of an "omics"-based system biology strategy should facilitate germplasm utilization and provide valuable resources for tea quality improvement.

Gossypium barbadense genome sequence provides insight into the evolution of extra-long staple fiber and specialized metabolites.

Of the two cultivated species of allopolyploid cotton, Gossypium barbadense produces extra-long fibers for the production of superior textiles. We sequenced its genome (AD)2 and performed a comparative analysis. We identified three bursts of retrotransposons from 20 million years ago (Mya) and a genome-wide uneven pseudogenization peak at 11-20 Mya, which likely contributed to genomic divergences. Among the 2,483 genes preferentially expressed in fiber, a cell elongation regulator, PRE1, is strikingly At biased and fiber specific, echoing the A-genome origin of spinnable fiber. The expansion of the PRE members implies a genetic factor that underlies fiber elongation. Mature cotton fiber consists of nearly pure cellulose. G. barbadense and G. hirsutum contain 29 and 30 cellulose synthase (CesA) genes, respectively; whereas most of these genes (>25) are expressed in fiber, genes for secondary cell wall biosynthesis exhibited a delayed and higher degree of up-regulation in G. barbadense compared with G. hirsutum, conferring an extended elongation stage and highly active secondary wall deposition during extra-long fiber development. The rapid diversification of sesquiterpene synthase genes in the gossypol pathway exemplifies the chemical diversity of lineage-specific secondary metabolites. The G. barbadense genome advances our understanding of allopolyploidy, which will help improve cotton fiber quality.