- Open Access
Transcriptome profiling of the fertile parent and sterile hybrid in tea plant flower buds
Hereditas volume 156, Article number: 12 (2019)
The tea plant is a crucial economic crop. The floral organ development consumes a large amount of nutrients, which affects the leaf yield. To understand the mechanism by which the tea plant produces sterile floral buds, we obtained a sterile tea plant by artificial hybridization. RNA-sequencing based transcriptome analysis was implemented in three samples to determine the differentially expressed genes (DEGs) related to flower development.
In this study, a total of 1991 DEGs were identified; 1057 genes were up-regulated and 934 genes were down-regulated in sterile hybrid floral buds. These were mainly distributed in the regulation of biological and metabolic processes. Significantly, auxin biosynthesis genes YUCCA, AUX1 and PIN were dramatically down-regulated, and ARF gene was up-regulated in the sterile hybrid floral buds, and flower development-related genes AP1, AP2 and SPL were changed. A total of 12 energy transfer-related genes were significantly decreased. Furthermore, the expression of 11 transcription factor genes was significantly different.
The transcriptome analysis suggested that the production of sterile floral buds is a complex bioprocess, and that low auxin-related gene levels result in the formation of sterile floral buds in the tea plant.
Sterility is a complex phenomenon in plants, the main factors of which include the accumulation of reactive oxygen species, abnormalities in energy metabolism, programmed cell death, imbalance in endogenous hormones and changes in environmental condition [1,2,3]. In Arabidopsis and cereal grains, floral organ degeneration and abiotic stresses result in sterility. Abortion or degeneration of developing stamens and pistil is the key mechanism used by plants to produce sterile flowers .
Among the sterility mechanisms of the plant, flower development is one of the main concerns. It requires early and later organ identity genes. PINFORMED (PIN1), PINOID (PID), YUCCA (YUC) and NAKED PINS IN YUC MUTANT (NPY) are necessary for flower development [5, 6]. Organ identity genes APETALA1 (AP1), APETALA2 (AP2), APETALA3 (AP3), PISTILLATA (PI) and AGAMOUS (AG) are essential for the ABC model in Arabidopsis . The mutation of these genes causes more petals, fewer stamens, fused floral organs, and valveless gynoecia. In addition, important genes for flower development include SPL, TPD1, AMS, DYT1, SHP, WUS [8,9,10,11,12,13]. Brassionsteroid and jasmonic acid play a positive role in promoting the formation of stamens and pollen. Gibberellin deficiencies are related to male sterility, and the formation of female flowers requires the presence of ethylene [14,15,16]. Moreover, flower development is regulated by the coordinated interaction of the transcription factor LEAFY and auxin . However, molecular data are limited for the tea plant; the molecular mechanism of flower development remains unknown.
The tea plant is a crucial cash crop widely distributed around the world. Tea leaves have been used to produce various tea beverages. The floral organ development results in the fall of the yield of tea leaves by consuming a large amount of nutrients. Cultivation of the sterile tea plant is key to increasing the yield. The tea plant is self-incompatible, thus we performed transcriptome sequencing and comparative analysis on three samples, including Foxiang2 (FBH), Fudingbaicha (MBH) and hybrid sterile flowers (ZDH). The aim was to analyze the differentially expressed genes between the fertile and sterile floral buds, and to identify their related bioprocesses and correlative factors. Our results will help to reveal important information on the mechanism of sterility in the tea plant.
Materials and methods
The plant materials were five-year-old tea plants (C. sinensis (L.) O. Kuntze) from the Tea Research Institute, Yunnan Academy of Agricultural Sciences called ‘Foxiang2’, used as FBH, ‘Fudingbaicha’, used as MBH, and a sterile hybrid, used as ZDH. Flower buds were separately stripped from three tea plants. All of the experiments were performed using three biological replicates. Flower buds were picked on October 16, 2016, after which they were frozen using liquid nitrogen and stored in a freezer at − 80 °C for subsequent mRNA analysis.
cDNA library construction and sequencing
The construction of cDNA libraries and transcriptome sequencing were completed by Beijing Novogene Technology (Beijing, China). Total RNA was used as the starting sample and was directly added to the 3′-terminal hydroxyl group and the uniquely structured complete phosphoryl groups at the 5′-terminal of the sRNA, followed by reverse transcription to synthesize cDNA. After polymerase chain reaction amplification, polyacrylamide gel electrophoresis was used to separate the target DNA fragments, and the gel was recovered, completing the cDNA library. The effective concentration of the library was > 2 nmol/L, and sequencing was performed using Illumina HiSeq 2000 after the library was certified. All of the experiments were performed using three replicates.
The original image data files obtained from the Illumina HiSeq 2000 were subjected to base calling analysis and converted into raw reads. Among the raw reads obtained from the sequencing, the low-quality reads with adaptors were processed to obtain clean reads. The clean reads were assembled separately, and TGICL was used to get the longest non-redundant unigenes. The transcriptome data were deposited to the NCBI SRA database (SRA accession: PRJNA503652).
Differentially expressed genes (DEGs) test
Differentially expressed genes (DEGs) analysis using DEGseq (three biological replicates per group). DESeq provides statistical routines for determining differential expression in digital gene expression data, using a model based on negative binomial distribution. The resulting P-values were adjusted using the Benjamin and Hochberg approach for controlling the false discovery rate. The genes with P-value < 0.05 found by DESeq were differentially expressed.
GO and KEGG analysis
DEGs were characterized according to Gene Ontology enrichment analysis. GO annotations were provided by the Blast2GO program. Then the GO classification graph was generated by the WEGO.
KEGG was used to analyze the biological process and unigenes annotation of pathway. The results were comparatively analyzed between the KEGG integrated database resource and our data.
Quantitative real-time PCR assays
Total RNA was isolated using TRlpure reagent (BioTeke, China) according to the manufacturer’s instructions. cDNA was synthesized from total RNA using a PrimeScript RT reagent kit (TaKaRa, Japan). The obtained cDNA was used as a template in SYBR green-based q-PCR (CFX-96, Bio-Rad, Hercules, CA, USA). GAPDH was used for normalization. The primers are shown in Table 1.
Morphological characteristics of sterile flower buds
The floral organ of the tea plant consists of a complete bisexual flower composed of a thalamus, calyxes, petals, stamens, and a pistil (Fig. 1a). In our study, the sterile flower buds were smaller than the male and female parent flower. Petals were improperly unfolding during the developmental process until flower buds abscission (Fig. 1a). The filaments were shorter than those of fertile flowers, and the anthers contained no pollen (Fig. 1b). In addition, sterile flowers contained two to four imperfect pistils, which were split into two to five smaller stigma (Fig. 1b).
Transcriptome profiling of the male/female parent and the hybrid bud
We used the Illumina HiSeq 2000 platform to create three cDNA libraries, FBH, MBH, and ZDH. As a result, a total of 64.2, 60.6 and 74.6 million clean reads were generated. The quality score (Q30) percentage was above 92%, with the GC content of each clean data above 43%. The ratio of mapped reads was 73.64, 75.27 and 71.43%, respectively (Table 2). Afterwards, Trinity was used in splicing for clean reads. A total of 268,289 transcripts were obtained, and the longest transcript for each gene was selected as the unigene, of which 173,248 were screened for differential expression analysis (DEGs).
DEGs and functional characterization
A total of 6395 DEGs were detected in FBH, MBH and ZDH. 1914 DEGs were screened out between FBH and MBH. Between ZDH and FBH, the number of DEGs was 5438, whereas it was 3208 between ZDH and MBH. 1991 DEGs were screened out in ZDH, the expression quantity was different with FBH and MBH (Fig. 2).
The unique DEGs were characterized using the GO and KEGG databases. The GO enrichment analysis categorized 1991 unigenes, which contained 47 subcategories of biological processes, cellular components, and molecular functions. For the biological process, “biological regulation”, “cellular process”, “metabolic process”, and “single organism process” were the most representative groups. In the subcategory of cellular component, “cell”, “cell part”, “macromolecular complex”, and “organelle” were the predominant groups. For the molecular function, “binding”, and “catalytic activity” were the most common. Moreover, “reproduction”, “reproductive process”, and “growth” were enriched in 24, 16, and 3 unigenes, respectively (Fig. 3).
The KEGG analysis revealed that 764 unigenes were mapped into KEGG, containing cellular processes, environmental information processes, generic information processes, metabolism, and organismal systems. The main KEGG pathways included “Carbon metabolism”, “Biosynthesis of amino acids”, “Starch and sucrose metabolism”, “Plant hormone signal transduction”, “Plant pathogen interaction”, and “Glycolysis/Gluconeogenesis” (Fig. 4).
Identification of auxin-related genes involved in DEGs
Auxin is essential for flower organ development . Among DEGs, the auxin biosynthesis gene YUCCA related to floral organ formation was down-regulated. The auxin flux-related PIN homolog gene related to gynoecium formation was down-regulated. In addition, we identified ten auxin response factors ARF (Cluster-23,036.113917, Cluster-23,036.14480, Cluster-23,036.14481, Cluster-23,036.85241, Cluster-23,036.27862, Cluster-23,036.29441, Cluster-23,036.70956, Cluster-23,036.66525, Cluster-23,036.54073 and Cluster-23,036.87364) related to flower maturation was up-regulated (Table 3).
Identification of flower development-related genes involved in DEGs
Abnormal flower development is the key factor responsible for the sterility of the tea plant. The formation of flowers is a key step in the plant life cycle, which is a complex process . Each stage of the process is regulated by flower development-related genes, specifically the ABC floral organ-identity genes . We identified the A-class genes AP1 and AP2 were up-regulated.In addition, flower development-related genes, such as floral organ formation SPL homolog genes (Cluster-23,036.89600, Cluster-23,036.17164, Cluster-23,036.96383, Cluster-23,036.48034 and Cluster-23,036.89141) were up-regulated (Table 4).
Identification of energy transfer-related genes in DEGs
Energy transfer is an important process of plant growth and development. It is the foundation of the plant’s life . We identified that some genes involved in energy transfer were differentially expressed, including ABC transporter B family member 1 (ABCB1), six transporting ATPase-related genes (ATPeF0D, ATPeF0O, ATPeF1B, ATPeV0A, ATPeV1C, and ATPeV1B), and five solute carrier-related genes (SLC2A8, SLC15A3, SLC25A11, SLC32A, and SLC35B3) (Table 5). A total of 12 genes were down-regulated. In particular, transporting ATPase subunit beta (ATPeF1B), solute carrier family 32 (SLC32A) and solute carrier family 35 (SLC35B3) were down-regulated 65.7, 60.9 and 54.9-fold, respectively, in ZDH compared with FBH. The expression levels of ATPeF1B, SLC32A and SLC35B3 were also down-regulated 96.8, 34.7 and 65.4-fold, respectively, in ZDH compared with MBH.
Identification of other transcription factor genes in DEGs
We also identified some transcription factor genes. Among these genes, the expression levels of transcription factor AS1 (AS1), transcription factor CPC (CPC), and nuclear transcription factor Y (NFYA) were up-regulated significantly. Furthermore, MADS-box transcription factor (K09264), two transcription factor bHLH (bHLH77 and bHLH79), three transcription factor MYB (GAMYB, MYBP and MYB21), transcription factor TGA6 (TGA6) and transcription factor GTE2 (GTE2) were down-regulated in ZDH (Table 6).
qPCR analysis of DEGs
To confirm the transcriptome sequencing results, six DEGs were selected for qPCR analysis. Among these genes, Cluster-23,036.57048 (AP2) and Cluster-23,036.14480 (ARF) genes were up-regulated, and Cluster-23,036.68881 (JAR1), Cluster-23,036.53622 (IAA7), Cluster-23,036.69380 (AUX1) and Cluster-23,036.74178 (ATL3) genes were down-regulated (Fig. 5). qRT-PCR results showed a consistent expression tendency compared with the RNA-Seq. The result further confirmed the reliability and accuracy of the transcriptome sequencing.
Here, we used transcriptome sequencing to explore candidate genes associated with sterile floral buds in the tea plant. The method has been applied to the study of the genome for the tea plant [22,23,24,25]. A total of 1991 DEGs were screened out from the comparison among the three cDNA libraries (Fig. 2). The complexity of DEGs function was demonstrated by Go and KEGG analysis (Figs. 3 and 4). Based on the analysis, 452 metabolism-related pathways were identified (Fig. 4), suggesting that metabolism may be crucial for the sterility mechanism in the tea plant.
In general, the flavin monooxygenase (YUCCA) is involved in the tryptophan-dependent pathway of auxin biosynthesis , mutants of YUCCA present aberrant flower phenotypes, short stamen filaments, and thus sterility . Our results show that the expression levels of YUCCA were down-regulated significantly (Table 3), the flower phenotype also appears as short stamen filaments (Fig. 1b). Thus, inferring that the low expression of YUCCA is a key factor affecting the biosynthesis of auxin leading to flower sterility. In addition, we also found that the expression of SPL homolog genes as an inhibiting factor of YUCCA was up-regulated in sterile flowers (Table 4). Previous studies showed that SPL represses YUCCA gene expression to regulate the development of lateral organs . Our results suggest that SPL regulates the homeostasis of auxin by inhibiting YUCCA, resulting in flower sterility.
ARF play pivotal roles in the growth of inflorescences, stamens, anthers, and pistils. Mutations in ARF6 and ARF8 caused male and female infertility . However, we found that the expression level of ten ARF homolog genes was up-regulated in sterile floral buds (Table 3). This is possibly due to the high expression of ARF inhibiting auxin signaling targets . Moreover, ARF acts as a positive or negative regulator by binding to the auxin response element TGTCTC , and it is possible that the Aux/IAA inhibitor is able to inhibit transcription through interaction with ARF .
Interestingly, we found that the A-class of ABC floral organ-identity genes AP1 and AP2 were up-regulated in sterile floral buds (Table 4). It has been reported that the initial expression of AP1 and AP2 was restricted to the first and second whorls, and was inhibited in the third and fourth whorls of flower development . Thus inferring that high expression of AP1 and AP2 represses flower development in first and second whorls. Detailed mechanisms will require further study.
As an important regulator of flower development, auxin is transported to each tissue through carriers and the AUX1/LAX . We found that the expression levels of AUX1/LAX (Table 3) and twelve energy transfer-related genes (Table 5) were down-regulated, suggesting that the transport of auxin may be impeded. This may also be an important factor in the sterile floral buds of the tea plant.
In summary, the present work provides four key factors for the development of sterile floral buds in the tea plant. We thus present a theoretical basis for further study of mechanisms by which sterile floral buds are produced.
Differentially expressed genes
Kyoto Encyclopedia of Genes and Genomes
Hybrid sterile flowers
Cheng Y, Dai X, Zhao Y. Auxin biosynthesis by the YUCCA flavin monooxygenases controls the formation of floral organs and vascular tissues in Arabidopsis. Genes Dev. 2006;20(13):1790–9.
Tang RS, Zheng JC, Jin ZQ, et al. Possible correlation between high temperature-induced floret sterility and endogenous levels of IAA, GAs and ABA in rice (Oryza sativa L.). Plant Growth Regul. 2008;54(1):37–43.
Wan C, Li S, Wen L, et al. Damage of oxidative stress on mitochondria during microspores development in Honglian CMS line of rice. Plant cell Rep. 2007;26(3):373–82.
Smith AR, Zhao D. Sterility Caused by Floral Organ Degeneration and Abiotic Stresses in Arabidopsis and Cereal Grains. Front Plant Sci. 2016;7:1503.
Cheng Y, Qin G, Dai X, et al. NPY genes and AGC kinases define two key steps in auxin-mediated organogenesis in Arabidopsis. Proc Natl Acad Sci. 2008;105(52):21017–21,022.
Okada K, Ueda J, Komaki MK, et al. Requirement of the auxin polar transport system in early stages of Arabidopsis floral bud formation. Plant Cell. 1991;3(7):677–84.
Ma H. Molecular genetic analyses of microsporogenesis and microgametogenesis in flowering plants. Annu Rev Plant Biol. 2005;56:393–434.
Ito T, Wellmer F, Yu H, et al. The homeotic protein AGAMOUS controls microsporogenesis by regulation of SPOROCYTELESS. Nat. 2004;430(6997):356–60.
Liljegren SJ, Ditta GS, Eshed Y, et al. SHATTERPROOF MADS-box genes control seed dispersal in Arabidopsis. Nat. 2000;404(6779):766–70.
Mayer KFX, Schoof H, Haecker A, et al. Role of WUSCHEL in regulating stem cell fate in the Arabidopsis shoot meristem. Cell. 1998;95(6):805–15.
Sorensen AM, Kröber S, Unte US, et al. The Arabidopsis ABORTED MICROSPORES (AMS) gene encodes a MYC class transcription factor. Plant J. 2003;33(2):413–23.
Yang SL, Jiang L, San Puah C, et al. Overexpression of TAPETUM DETERMINANT1 alters the cell fates in the Arabidopsis carpel and tapetum via genetic interaction with excess microsporocytes1/extra sporogenous cells. Plant Physiol. 2005;139(1):186–91.
Zhang W, Sun Y, Timofejeva L, et al. Regulation of Arabidopsis tapetum development and function by DYSFUNCTIONAL TAPETUM1 (DYT1) encoding a putative bHLH transcription factor. Dev. 2006;133(16):3085–95.
Cecchetti V, Altamura MM, Falasca G, et al. Auxin regulates Arabidopsis anther dehiscence, pollen maturation, and filament elongation. Plant Cell. 2008;20(7):1760–74.
Mandaokar A, Thines B, Shin B, et al. Transcriptional regulators of stamen development in Arabidopsis identified by transcriptional profiling. Plant J. 2006;46(6):984–1008.
Szekeres M, Németh K, Koncz-Kálmán Z, et al. Brassinosteroids rescue the deficiency of CYP90, a cytochrome P450, controlling cell elongation and de-etiolation in Arabidopsis. Cell. 1996;85(2):171–82.
Yamaguchi N, Wu MF, Winter CM, et al. LEAFY and polar auxin transport coordinately regulate Arabidopsis flower development. Plants. 2014;3(2):251–65.
Cecchetti V, Altamura M M, Falasca G, et al. Auxin regulates Arabidopsis anther dehiscence, pollen maturation, and filament elongation[J]. The Plant Cell, 2008;20(7):1760-74.
Grossniklaus U, Schneitz K. The molecular and genetic basis of ovule and megagametophyte development. Semin Cell Dev Biol. 1998;9(2):227–38.
Ainsworth C, Crossley S, Buchanan-Wollaston V, et al. Male and female flowers of the dioecious plant sorrel show different patterns of MADS box gene expression. Plant Cell. 1995;7(10):1583–98.
Baker LA, Habershon S. Photosynthetic pigment-protein complexes as highly connected networks: implications for robust energy transport. Proc Math Phys Eng Sci. 2017;473(2201):20170112.
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 Nat Acad Sci. 2018;115(18):E4151–8.
Li H, Huang W, Wang GL, et al. Transcriptomic analysis of the biosynthesis, recycling, and distribution of ascorbic acid during leaf development in tea plant (Camellia sinensis (L.) O. Kuntze). Sci Rep. 2017;7:46212.
Liu SC, Jin JQ, Ma JQ, et al. Transcriptomic analysis of tea plant responding to drought stress and recovery. PloS One. 2016;11(1):e0147306.
Tai Y, Wei C, Yang H, et al. Transcriptomic and phytochemical analysis of the biosynthesis of characteristic constituents in tea (Camellia sinensis) compared with oil tea (Camellia oleifera). BMC Plant Biol. 2015;15(1):190.
Abu-Zaitoon YM. Phylogenetic analysis of putative genes involved in the tryptophan-dependent pathway of auxin biosynthesis in rice. Appl Biochem Biotechnol. 2014;172(5):2480–95.
Salinas-Grenet H, Herrera-Vásquez A, Parra S, et al. Modulation of Auxin Levels in Pollen Grains Affects Stamen Development and Anther Dehiscence in Arabidopsis. Int J Mol Sci. 2018;19(9):2480.
Li LC, Qin GJ, Tsuge T, et al. SPOROCYTELESS modulates YUCCA expression to regulate the development of lateral organs in Arabidopsis. New Phytol. 2008;179(3):751–64.
Nagpal P, Ellis CM, Weber H, et al. Auxin response factors ARF6 and ARF8 promote jasmonic acid production and flower maturation. Dev. 2005;132(18):4107–18.
Hofmann NR. A Genome-wide Approach to Understanding a Non-Canonical ARF. Plant Cell. 2017;29(8):1798–9.
Gray WM, Kepinski S, Rouse D, et al. Auxin regulates SCF TIR1-dependent degradation of AUX/IAA proteins. Nat. 2001;414(6861):271.
Guilfoyle TJ. The PB1 domain in auxin response factor and Aux/IAA proteins: a versatile protein interaction module in the auxin response. Plant Cell. 2015;27(1):33–43.
Irish VF. The flowering of Arabidopsis flower development[J]. Plant J. 2010;61(6):1014–28.
Novoselova ES, Mironova VV, Omelyanchuk NA, et al. Mathematical modeling of auxin transport in protoxylem and protophloem of Arabidopsis thaliana root tips[J]. J Bioinform Comput Biol. 2013;11(01):1340010.
This work was supported by the National Natural Science Foundation of China (31460216), Personnel Training Program of Yunnan Province (2015HB105) and State Key Laboratory of Tea Plant Biology and Utilization (SKLTOF20150105, SKLTOF20170115).
This work was funded by NSFC31460216, 2015HB105, SKLTOF20150105, and SKLTOF20170115.
Availability of data and materials
We have provided detailed information about the materials and methods in our manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, L., Qu, H., Xia, L. et al. Transcriptome profiling of the fertile parent and sterile hybrid in tea plant flower buds. Hereditas 156, 12 (2019) doi:10.1186/s41065-019-0090-z
- Tea plant
- Sterile floral buds
- Differentially expressed genes