PacBio Iso-Seq and RNA-Seq Reveal Gene Expression Divergence in Early-Flowering Gynoecious and Androecious Ramie Mutants

Abstract

Ramie (Boehmeria nivea L.) is an economically valuable bast fiber crop, yet a high-quality full-length transcriptome reference for sex-divergent lines remains lacking. To fill this gap, we combined PacBio Iso-Seq with Illumina RNA-Seq to construct a full-length transcriptome for the gynoecious mutant GBN09 and compared transcriptomes of four tissues (leaves, bast fibers, stems, and flowers) between GBN09 and the androecious mutant GBN10. The analysis generated 63,082 high-quality isoforms and annotated 30,885 genes. Regulatory element mining revealed 565 transcription factors (predominantly IAA, AP2/ERF, and WRKY families), 16,617 long non-coding RNAs (lncRNAs), and 31,759 simple sequence repeats (SSRs). Transcriptomic comparisons identified 8414 differentially expressed genes (DEGs), with the most profound transcriptional shifts in floral and young stem tissues. Enrichment analysis linked these DEGs primarily to linoleic acid metabolism, phytohormone signaling cascades, and secondary metabolite biosynthesis, highlighting their roles in sexual differentiation. Furthermore, we detected 274 sex-biased differential alternative splicing (DAS) events, many of which mapped to cell wall biogenesis and core auxin signaling components such as AFB2 and NPY2. These results indicate that post-transcriptional modifications provide an additional regulatory layer during ramie floral sex differentiation. Validated by qRT-PCR, this dataset establishes a robust molecular resource and identifies candidate genes to elucidate sex differentiation and support molecular breeding in ramie.

Keywords

Ramie / Sex differentiation / PacBio Iso-Seq / RNA-Seq

References

[1] Huang C, Zhou J, Jie Y, Xing H, Zhong Y, Yu W, et al. A Ramie bZIP Transcription Factor BnbZIP2 Is Involved in Drought, Salt, and Heavy Metal Stress Response. DNA Cell Biol. 2016, 35, 776-786. DOI: 10.1089/dna.2016.3251

[2] Fu H, Wang C, Cui G, She W, Zhao L. Ramie Yield Estimation Based on UAV RGB Images. Sensors 2021, 21, 669. DOI: 10.3390/s21020669

[3] Li F, Zeng Z, Huang R, Wang Y, Liu T. Identification of proteins associated with bast fiber growth of ramie by differential proteomic analysis. BMC Genom. 2021, 22, 865. DOI: 10.1186/s12864-021-08195-9

[4] Zeng Z, Li F, Huang R, Wang Y, Liu T. Phosphoproteome analysis reveals an extensive phosphorylation of proteins associated with bast fiber growth in ramie. BMC Plant Biol. 2021, 21, 473. DOI: 10.1186/s12870-021-03252-7

[5] Kandimalla R, Kalita S, Choudhury B, Devi D, Kalita D, Kalita K, et al. Fiber from ramie plant (Boehmeria nivea): A novel suture biomaterial. Mater. Sci. Eng. C 2016, 62, 816-822. DOI: 10.1016/j.msec.2016.02.040

[6] Liu T, Zhu S, Tang Q, Tang S. Genome-wide transcriptomic profiling of ramie (Boehmeria nivea L. Gaud) in response to cadmium stress. Gene 2015, 558, 131-137. DOI: 10.1016/j.gene.2014.12.057

[7] Mu L, Cai M, Wang Z, Liu J, Liu T, Wanapat M, et al. Assessment of ramie leaf (Boehmeria nivea L. Gaud) as an animal feed supplement in P.R. China. Trop. Anim. Health Prod. 2020, 52, 115-121. DOI: 10.1007/s11250-019-01997-w

[8] Zhang H, He Y, Yu L, Tang S, Tan Z, Zhou C, et al. The potential of ramie as forage for ruminants: Impacts on growth, digestion, ruminal fermentation, carcass characteristics and meat quality of goats. Anim. Sci. J. 2019, 90, 481-492. DOI: 10.1111/asj.13165

[9] Bjornson M, Kajala K, Zipfel C, Ding P. Low-cost and High-throughput RNA-Seq Library Preparation for Illumina Sequencing from Plant Tissue. Bio-Protocol 2020, 10, e3799. DOI: 10.21769/BioProtoc.3799

[10] Liu YH, Wang L, Yu L. The principle and application of the single-molecule real-time sequencing technology. Hereditas 2015, 37, 259-268. DOI: 10.16288/j.yczz.14-323

[11] Li W, Liu Z, Feng H, Yang J, Li C. Characterization of the Gene Expression Profile Response to Drought Stress in Populus ussuriensis Using PacBio SMRT and Illumina Sequencing. Int. J. Mol. Sci. 2022, 23, 3840. DOI: 10.3390/ijms23073840

[12] Savarese M, Qureshi T, Torella A, Laine P, Giugliano T, Jonson PH, et al. Identification and Characterization of Splicing Defects by Single-Molecule Real-Time Sequencing Technology (PacBio). J. Neuromuscul. Dis. 2020, 7, 477-481. DOI: 10.3233/JND-200523

[13] Chen J, Pei Z, Dai L, Wang B, Liu L, An X, et al. Transcriptome profiling using pyrosequencing shows genes associated with bast fiber development in ramie (Boehmeria nivea L.) . BMC Genom. 2014, 15, 919. DOI: 10.1186/1471-2164-15-919

[14] She W, Zhu S, Jie Y, Xing H, Cui G. Expression Profiling of Cadmium Response Genes in Ramie (Boehmeria nivea L.) Root. Bull. Environ. Contam. Toxicol. 2015, 94, 453-459. DOI: 10.1007/s00128-015-1502-z

[15] Shang L, Zhou Y, Wen S, Wang K, Li Y, Zhang M, et al. Construction of heat stress regulation networks based on Illumina and SMRT sequencing data in potato. Front. Plant Sci. 2023, 14, 1271084. DOI: 10.3389/fpls.2023.1271084

[16] Jian H, Sun H, Liu R, Zhang W, Shang L, Wang J, et al. Construction of drought stress regulation networks in potato based on SMRT and RNA sequencing data. BMC Plant Biol. 2022, 22, 381. DOI: 10.1186/s12870-022-03758-8

[17] Yang F, Lv G. Characterization of the gene expression profile response to drought stress in Haloxylon using PacBio single-molecule real-time and Illumina sequencing. Front. Plant Sci. 2022, 13, 981029. DOI: 10.3389/fpls.2022.981029

[18] Zhu Z, Chen H, Xie K, Liu C, Li L, Liu L, et al. Characterization of Drought-Responsive Transcriptome During Seed Germination in Adzuki Bean (Vigna angularis L.) by PacBio SMRT and Illumina Sequencing. Front. Genet. 2020, 11, 996. DOI: 10.3389/fgene.2020.00996

[19] Wei M, Wang B, Li C, Li X, He C, Li Y. Integrated PacBio SMRT and Illumina sequencing uncovers transcriptional and physiological responses to drought stress in whole-plant Nitraria tangutorum. Front. Genet. 2024, 15, 1474259. DOI: 10.3389/fgene.2024.1474259

[20] Xiong Y, Xiang X, Xiao C, Zhang N, Cheng H, Rao S, et al. Illumina RNA and SMRT Sequencing Reveals the Mechanism of Uptake and Transformation of Selenium Nanoparticles in Soybean Seedlings. Plants 2023, 12, 789. DOI: 10.3390/plants12040789

[21] Rhoads A, Au KF. PacBio Sequencing and its Applications. Genom. Proteom. Bioinform. 2015, 13, 278-289. DOI: 10.1016/j.gpb.2015.08.002

[22] Au KF, Underwood JG, Lee L, Wong WH. Improving PacBio Long Read Accuracy by Short Read Alignment. PLoS ONE 2012, 7, e46679. DOI: 10.1371/journal.pone.0046679

[23] Li W, Godzik A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658-1659. DOI: 10.1093/bioinformatics/btl158

[24] Liaquat F, Munis MFH, Arif S, Haroon U, Shi J, Saqib S, et al. PacBio Single-Molecule Long-Read Sequencing Reveals Genes Tolerating Manganese Stress in Schima superba Saplings. Front. Genet. 2021, 12, 635043. DOI: 10.3389/fgene.2021.635043

[25] Zheng Y, Jiao C, Sun H, Rosli HG, Pombo MA, Zhang P, et al. iTAK: A Program for Genome-wide Prediction and Classification of Plant Transcription Factors, Transcriptional Regulators, and Protein Kinases. Mol. Plant 2016, 9, 1667-1670. DOI: 10.1016/j.molp.2016.09.014

[26] Sun L, Luo H, Bu D, Zhao G, Yu K, Zhang C, et al. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013, 41, e166. DOI: 10.1093/nar/gkt646

[27] Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, et al. CPC: Assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007, 35, W345-W349. DOI: 10.1093/nar/gkm391

[28] Zhao N, Zhou E, Miao Y, Xue D, Wang Y, Wang K, et al. High-quality faba bean reference transcripts generated using PacBio and Illumina RNA-Seq data. Sci. Data 2024, 11, 359. DOI: 10.1038/s41597-024-03204-4

[29] Li A, Zhang J, Zhou Z. PLEK: A tool for predicting long non-coding RNAs and messenger RNAs based on an improved k-mer scheme. BMC Bioinform. 2014, 15, 311. DOI: 10.1186/1471-2105-15-311

[30] Li W, Fu Y, Lv W, Zhao S, Feng H, Shao L, et al. Characterization of the early gene expression profile in Populus ussuriensis under cold stress using PacBio SMRT sequencing integrated with RNA-Seq reads. Tree Physiol. 2022, 42, 646-663. DOI: 10.1093/treephys/tpab130

[31] Weng X, Zhou X, Xie S, Gu J, Wang ZY. Identification of cassava alternative splicing-related genes and functional characterization of MeSCL30 involvement in drought stress. Plant Physiol. Biochem. 2021, 160, 130-142. DOI: 10.1016/j.plaphy.2021.01.016

[32] Foissac S, Sammeth M. ASTALAVISTA: Dynamic and flexible analysis of alternative splicing events in custom gene datasets. Nucleic Acids Res. 2007, 35, W297-W299. DOI: 10.1093/nar/gkm311

[33] Robinson MD, McCarthy DJ, Smyth GK. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139-140. DOI: 10.1093/bioinformatics/btp616

[34] Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674-3676. DOI: 10.1093/bioinformatics/bti610

[35] Kanehisa M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27-30. DOI: 10.1093/nar/28.1.27

[36] Trincado JL, Entizne JC, Hysenaj G, Singh B, Skalic M, Elliott DJ, et al. SUPPA2: Fast, accurate, and uncertainty-aware differential splicing analysis across multiple conditions. Genome Biol. 2018, 19, 40. DOI: 10.1186/s13059-018-1417-1

[37] Figueiredo MRAD, Strader LC. Intrinsic and extrinsic regulators of Aux/IAA protein degradation dynamics. Trends Biochem. Sci. 2022, 47, 865-874. DOI: 10.1016/j.tibs.2022.06.004

[38] Chen Q, Li N, Cui X, Ge F. AP2/ERF transcription factors regulate the biosynthesis of terpenoids, phenolics, and alkaloids in plants. Hortic. Res. 2026, 13, uhaf280. DOI: 10.1093/hr/uhaf280

[39] Song H, Duan Z, Zhang J. WRKY transcription factors modulate flowering time and response to environmental changes. Plant Physiol. Biochem. 2024, 210, 108630. DOI: 10.1016/j.plaphy.2024.108630

[40] Pannell JR. Plant Sex Determination. Curr. Biol. 2017, 27, R191-R197. DOI: 10.1016/j.cub.2017.01.052

[41] Andreuzza S. Y Keep Your X? Insights into the Genetic Basis of Plant Sex Chromosome Evolution. Plant Cell 2020, 32, 1769-1770. DOI: 10.1105/tpc.20.00265

[42] Li N, Meng Z, Tao M, Wang Y, Zhang Y, Li S, et al. Comparative transcriptome analysis of male and female flowers in Spinacia oleracea L. BMC Genom. 2020, 21, 850. DOI: 10.1186/s12864-020-07277-4

[43] Zhang H, Li W, Zhang X, Pan R, Tang M, Peng Y, et al. Integrating transcriptome and metabolome analyses to characterize flower development in Salix variegata franch., a typical dioecious plant. Sci. Rep. 2025, 15, 8010. DOI: 10.1038/s41598-025-91317-0

[44] Zhang K, Cao Y, Guo X, Kong F, Sun H, Jing T, et al. Comparative transcriptome analysis of differentially expressed genes and pathways in male and female flowers of Fraxinus mandshurica. PLoS ONE 2024, 19, e0308013. DOI: 10.1371/journal.pone.0308013