Genetic analysis of the early bud flush trait of tea plants (Camellia sinensis) in the cultivar ‘Emei Wenchun’ and its open-pollinated offspring

Abstract

The timing of bud flush (TBF) in the spring is one of the most important agronomic traits of tea plants (Camellia sinensis). In this study, we designed an open-pollination breeding program using ‘Emei Wenchun’ (EW, a clonal tea cultivar with extra-early TBF) as a female parent. A half-sib population (n = 388) was selected for genotyping using specific-locus amplified fragment sequencing. The results enabled the identification of paternity for 294 (75.8%) of the offspring, including 11 (2.8%) from EW selfing and 217 (55.9%) assigned to a common father, ‘Chuanmu 217’ (CM). The putative EW × CM full-sib population was used to construct a linkage map. The map has 4244 markers distributed in 15 linkage groups, with an average marker distance of 0.34 cM. A high degree of collinearity between the linkage map and physical map was observed. Sprouting index, a trait closely related to TBF, was recorded for the offspring population in 2020 and 2021. The trait had moderate variation, with coefficients of variation of 18.5 and 17.6% in 2020 and 2021, respectively. Quantitative trait locus (QTL) mapping that was performed using the linkage map identified two major QTLs and three minor QTLs related to the sprouting index. These QTLs are distributed on Chr3, Chr4, Chr5, Chr9, and Chr14 of the reference genome. A total of 1960 predicted genes were found within the confidence intervals of QTLs, and 22 key candidate genes that underlie these QTLs were preliminarily screened. These results are important for breeding and understanding the genetic base of the TBF trait of tea plants.

References

[1] Drew L. The growth of tea. Nature. 2019; 566: S2-2.

[2] Zhang X, Chen S, Shi L et al. Haplotype-resolved genome assembly provides insights into evolutionary history of the tea plant Camellia sinensis . Nat Genet. 2021; 53: 1250-9.

[3] Xia E, Tong W, Hou Y et al. The reference genome of tea plant and resequencing of 81 diverse accessions provide insights into its genome evolution and adaptation. Mol Plant. 2020; 13: 1013-26.

[4] Wang X, Feng H, Chang Y et al. Population sequencing enhances understanding of tea plant evolution. Nat Commun. 2020; 11: 4447.

[5] Yu X, Xiao J, Chen S et al. Metabolite signatures of diverse Camellia sinensis tea populations. Nat Commun. 2020; 11: 1829.

[6] Dong C, Li F, Yang T et al. Theanine transporters identified in tea plants (Camellia sinensis L.) . Plant J. 2020; 101: 57-70.

[7] Fuchinoue Y . Analysis of self-incompatibility alleles of major varieties of tea. Jpn Agric Res Q. 1979; 13: 43-8.

[8] Yang Y, Liang Y . Chinese Clone Tea Plant Cultivars . Shanghai: Shanghai Scientific and Technical Publisher; 2014.

[9] Wang X, Wang L, Hao X et al. Tea genetics and breeding. Progress during the 13th five-year plan period and development direction in the 14th five-year plan period. China Tea. 2021; 43: 50-7.

[10] Wang RJ, Gao XF, Yang J et al. Genome-wide association study to identify favorable SNP allelic variations and candidate genes that control the timing of spring bud flush of tea (Camellia sinensis) using SLAF-seq . J Agric Food Chem. 2019; 67: 10380-91.

[11] Wang X, Hao X, Ma C et al. Identification of differential gene expression profiles between winter dormant and sprouting axillary buds in tea plant (Camellia sinensis) by suppression subtractive hybridization . Tree Genet Genomes. 2014; 10: 1149-59.

[12] Hao X, Yang Y, Yue C et al. Comprehensive transcriptome analyses reveal differential gene expression profiles of Camellia sinensis axillary buds at para-, endo-, ecodormancy, and bud flush stages. Front Plant Sci. 2017; 8: 553.

[13] Tan L, Wang L, Zhou B et al. Comparative transcriptional analysis reveled genes related to short winter-dormancy regulation in Camellia sinensis . Plant Growth Regul. 2020; 92: 401-15.

[14] Tan LQ, Wang LY, Xu LY et al. SSR-based genetic mapping and QTL analysis for timing of spring bud flush, young shoot color, and mature leaf size in tea plant (Camellia sinensis) . Tree Genet Genomes. 2016; 12: 52.

[15] Liu Q-L, Li X-X, Xu Y-H et al. Analysis of inbreeding depression in the tea plant (Camellia sinensis) based on the self-compatible cultivar ‘Ziyan’ . Sci Hortic. 2021; 284: 110120.

[16] Stevens EL, Heckenberg G, Roberson ED et al. Inference of relationships in population data using identity-by-descent and identity-by-state. PLoS Genet. 2011; 7: e1002287.

[17] Miller-Crews I, Matz MV, Hofmann HA . A 2b-RAD parentage analysis pipeline for complex and mixed DNA samples. Forensic Sci Int Genet. 2021; 55: 102590.

[18] Kalinowski ST, Taper ML, Marshall TC . Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007; 16: 1099-106.

[19] Xu L-Y, Wang L-Y, Wei K et al. High-density SNP linkage map construction and QTL mapping for flavonoid-related traits in a tea plant (Camellia sinensis) using 2b-RAD sequencing . BMC Genomics. 2018; 19: 955.

[20] Luo Y-P, Tang M, Cai W-z et al. Study on the optimum machine-plucking period for high quality tea. J Tea Sci. 2008; 28: 9-13.

[21] van Ooijen JW . MapQTL6, software for the mapping of quantitative trait loci in experimental populations . Wageningen: Kyazma BV.

[22] Howard NP, Peace C, Silverstein KA et al. The use of shared haplotype length information for pedigree reconstruction in asexually propagated outbreeding crops, demonstrated for apple and sweet cherry. Hortic Res. 2021; 8: 1-13.

[23] Zhao X, Liu Z . Study on crossing-breeding in tea - 1. Seed setting rate and fruit traits of different hybrid combinations. Tea. 1981; 3: 13-7.

[24] Ma JQ, Ming-Zhe Y, Chun-Lei M et al. Construction of a SSR-based genetic map and identification of QTLs for catechins content in tea plant (Camellia sinensis) . PLoS One. 2014; 9: 1-11.

[25] Koech RK, Malebe PM, Nyarukowa C et al. Functional annotation of putative QTL associated with black tea quality and drought tolerance traits. Sci Rep. 2019; 9: 1-11.

[26] Huang R, Wang J-Y, Yao M-Z et al. Quantitative trait loci mapping for free amino acid content using an albino population and SNP markers provides insight into the genetic improvement of tea plants. Hortic Res. 2022; 9: 1-13.

[27] Cooke JE, Eriksson ME, Junttila O . The dynamic nature of bud dormancy in trees: environmental control and molecular mechanisms. Plant Cell Environ. 2012; 35: 1707-28.

[28] Rushton DL, Tripathi P, Rabara RC et al. WRKY transcription factors: key components in abscisic acid signalling. Plant Biotechnol J. 2012; 10: 2-11.

[29] Yordanov YS, Ma C, Strauss SH et al. EARLY BUD-BREAK 1 (EBB1) is a regulator of release from seasonal dormancy in poplar trees. Proc Natl Acad Sci USA 2014; 111: 10001-6.

[30] Kim SY, Warpeha KM, Huber SC . The brassinosteroid receptor kinase, BRI1, plays a role in seed germination and the release of dormancy by cold stratification. J Plant Physiol. 2019; 241: 153031.

[31] Hao X, Hu T, Bo W et al. Gene characterization and expression analysis reveal the importance of auxin signaling in bud dormancy regulation in tea plant. J Plant Growth Regul. 2019; 38: 225-40.

[32] Prudencio ÁS, Werner O, Martínez-García PJ et al. DNA methylation analysis of dormancy release in almond (Prunus dulcis) flower buds using epi-genotyping by sequencing . Int J Mol Sci. 2018; 19: 3542.

[33] Chen J, Liu J, Jiang J et al. F-box protein CFK1 interacts with and degrades de novo DNA methyltransferase in Arabidopsis . New Phytol. 2021; 229: 3303-17.

[34] Liu H-Y, Yu X, Cui D-Y et al. The role of water channel proteins and nitric oxide signaling in rice seed germination. Cell Res. 2007; 17: 638-49.

[35] Yooyongwech S, Horigane AK, Yoshida M et al. Changes in aquaporin gene expression and magnetic resonance imaging of water status in peach tree flower buds during dormancy. Physiol Plant. 2008; 134: 522-33.

[36] Zhao Y-X, Zhao Z, Chen C-S et al. Characterization of self-incompatibility and expression profiles of CsMCU2 related to pollination in different varieties of tea plants. Sci Hortic. 2022; 293: 110693.

[37] Wachira FN, Kamunya S . Pseudo-self-incompatibility in some tea clones (Camellia sinensis (L.) O. Kuntze) . J Hortic Sci Biotechnol. 2005; 80: 716-20.

[38] Zhang C-C, Wang L-Y, Wei K et al. Transcriptome analysis reveals self-incompatibility in the tea plant (Camellia sinensis) might be under gametophytic control . BMC Genomics. 2016; 17: 289.

[39] Wei C, Yang H, Wang S et al. Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality . Proc Natl Acad Sci USA 2018; 115: E4151-8.

[40] Sun X, Liu D, Zhang X et al. SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS One. 2013; 8: e58700.

[41] Li H, Durbin R . Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25: 1754-60.

[42] McKenna A, Hanna M, Banks E et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20: 1297-303.

[43] Danecek P, Auton A, Abecasis G et al. The variant call format and VCFtools. Bioinformatics. 2011; 27: 2156-8.

[44] Chang CC, Chow CC, Tellier LC et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015; 4: 7.

[45] Alexander DH, Novembre J, Lange K . Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009; 19: 1655-64.

[46] Price AL, Patterson NJ, Plenge RM et al. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006; 38: 904-9.

[47] Liu D, Ma C, Hong W et al. Construction and analysis of high-density linkage map using high-throughput sequencing data. PLoS One. 2014; 9: e98855.

[48] Lei X, Wang Y, Zhou Y et al. TeaPGDB: tea plant genome database. Beverage Plant Res. 2021; 1: 5.

[49] Kim D, Langmead B, Salzberg SL . HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015; 12: 357-60.

[50] Pertea M, Pertea GM, Antonescu CM et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015; 33: 290-5.

[51] Love MI, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15: 1-21.