Integrated RNA-seq and small RNA-seq profiling reveals tissue-associated transcriptional programs and candidate miRNA–mRNA relationships across tea floral tissues

Abstract

Introduction Tea flowers comprise anatomically and functionally specialized tissues, yet the transcriptional and post-transcriptional programs that distinguish these tissues remain poorly characterized. Methods In this study, flower buds of Camellia sinensis cv. ‘Fudingdabaicha’ at the bud-white stage were dissected into seven tissues: ovary, outer green petal, inner white petal, filament, sepal, anther, and style. RNA-seq and small RNA-seq libraries were prepared from the same 21 biological samples, enabling direct comparison and integration of mRNA and miRNA expression profiles. Results Across the 21 pairwise tissue comparisons, 33,333 nonredundant differentially expressed genes were identified, revealing substantial variation in transcriptomic divergence among tissues. The outer green and inner white petals had the most similar transcriptomic profiles, and the genes that distinguished them were enriched mainly in plastid-associated and photosynthetic processes. By contrast, the anther and style showed the greatest transcriptomic divergence, involving metabolic, structural, and signaling functions. Weighted gene co-expression network analysis identified 24 modules, including three representative modules whose eigengenes showed positive correlations with the anther, sepal, and style traits. These modules were enriched in processes related to pollen germination and polarized growth, cell wall remodeling and organ growth, and RNA metabolism and chromatin-associated developmental regulation, respectively. RT-qPCR analysis of six representative genes showed tissue expression patterns broadly concordant with the RNA-seq data. Small RNA profiling identified 215 miRNAs, including 54 homologs of known plant miRNA families and 161 putative novel miRNAs; 175 miRNAs were differentially expressed among tissues. Intersecting psRNATarget predictions with candidate hub genes from the three representative modules yielded 45 candidate miRNA–mRNA relationships involving 39 miRNAs and 28 candidate hub genes. Six representative predicted pairs showed significant negative correlations across the matched samples. Discussion These findings characterize tissue-dependent mRNA and miRNA expression programs in tea flowers and prioritize candidate genes and miRNA–mRNA pairs for functional studies of floral tissue differentiation and reproductive development.