Benchmarking long-read RNA sequencing for de novo transcriptome assembly in non-model plant species: insights from Moricandia arvensis
- 1Universidad de Granada
Abstract
De novo transcriptome assembly is the standard approach for constructing a reference transcriptome in non-model plants that lack a high-quality genome, yet short-read assemblies struggle to resolve full-length isoforms. Long-read Iso-Seq (PacBio) captures full-length transcripts directly, but its use as a primary reference and the choice of downstream assembly pipeline remains poorly benchmarked. Here, we construct genome-free Iso-Seq reference transcriptomes for two organs, flower and leaf, of the non-model species Moricandia arvensis (L.) DC. (Brassicaceae), and systematically compare pipeline strategies combining Iso-Seq clustering, CD-HIT redundancy reduction, and Cogent graph-based reconstruction, benchmarked by BUSCO completeness, RSEM short-read mapping, and TransDecoder ORF completeness. We find that the optimal pipeline is organ specific. For flower, CD-HIT pre-filtering followed by Cogent reconstruction produced a high-quality reference (95.3% BUSCO complete). For the leaf, the same Cogent step was detrimental, reducing BUSCO completeness from 90.1% to 78.0% by incorrectly merging distinct genes; therefore, CD-HIT at 95% identity without reconstruction was retained. We trace this divergence to organ-specific input-data characteristics: leaf transcripts show extreme full-length-read expression skew and predominantly single-isoform gene support, depriving Cogent's graph algorithm of the multi-isoform evidence it requires. We find that the concentration of full-length reads among the most highly expressed transcripts predicts pipeline suitability before reconstruction, with per-transcript read depth acting as a necessary but non-discriminating floor. Because the leaf reference lacked gene-level structure, we further recovered gene-isoform grouping using expression-aware read-clustering (Corset), which preserved completeness while restoring the paralog structure expected of a paleopolyploid genome and outperformed sequence-only clustering. Our results provide a robust, genome-free framework for constructing full-length reference transcriptomes in non-model plant species and demonstrate that pipeline choice must be evaluated per organ rather than assuming one size fits all.
Related
- SpliceCOV: Supervised Filtering of RNA-Seq-Derived Transcript Features Improves Transcriptome Assembly
- Impact of HIV-1 infection on the microglia transcriptome: induction of microRNA-23b leads to G3BP2-mediated stress granules assembly disruption
- Top2All: Benchmarking sparse-to-transcriptome reconstruction from high-expression profiles across sex- and sex-steroid-defined conditions
- Chromosome-level genome assembly and transcriptome resource for the medicinal leech Hirudo verbana
- Benchmarking Microbiome Methods for Improved Metagenomic Sequencing of Lower Respiratory Tract Samples
BibTeX
@misc{sharma2026benchmarking,
title = {{Benchmarking long-read RNA sequencing for de novo transcriptome assembly in non-model plant species: insights from Moricandia arvensis}},
author = {S. Sharma and M. Hackenberg and L. Navarro and E. Narbona and A. Gonzalez-Megias and C. Armas and J. M. Gomez and F. Perfectti},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.08.750093},
url = {https://www.biorxiv.org/content/10.64898/2026.09.08.750093v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Benchmarking long-read RNA sequencing for de novo transcriptome assembly in non-model plant species: insights from Moricandia arvensis
AU - S. Sharma
AU - M. Hackenberg
AU - L. Navarro
AU - E. Narbona
AU - A. Gonzalez-Megias
AU - C. Armas
AU - J. M. Gomez
AU - F. Perfectti
PY - 2026
DA - 2026/09/08
PB - bioRxiv
DO - 10.64898/2026.09.08.750093
UR - https://www.biorxiv.org/content/10.64898/2026.09.08.750093v1
AB - De novo transcriptome assembly is the standard approach for constructing a reference transcriptome in non-model plants that lack a high-quality genome, yet short-read assemblies struggle to resolve full-length isoforms. Long-read Iso-Seq (PacBio) captures full-length transcripts directly, but its use as a primary reference and the choice of downstream assembly pipeline remains poorly benchmarked. Here, we construct genome-free Iso-Seq reference transcriptomes for two organs, flower and leaf, of the non-model species Moricandia arvensis (L.) DC. (Brassicaceae), and systematically compare pipeline strategies combining Iso-Seq clustering, CD-HIT redundancy reduction, and Cogent graph-based reconstruction, benchmarked by BUSCO completeness, RSEM short-read mapping, and TransDecoder ORF completeness. We find that the optimal pipeline is organ specific. For flower, CD-HIT pre-filtering followed by Cogent reconstruction produced a high-quality reference (95.3% BUSCO complete). For the leaf, the same Cogent step was detrimental, reducing BUSCO completeness from 90.1% to 78.0% by incorrectly merging distinct genes; therefore, CD-HIT at 95% identity without reconstruction was retained. We trace this divergence to organ-specific input-data characteristics: leaf transcripts show extreme full-length-read expression skew and predominantly single-isoform gene support, depriving Cogent's graph algorithm of the multi-isoform evidence it requires. We find that the concentration of full-length reads among the most highly expressed transcripts predicts pipeline suitability before reconstruction, with per-transcript read depth acting as a necessary but non-discriminating floor. Because the leaf reference lacked gene-level structure, we further recovered gene-isoform grouping using expression-aware read-clustering (Corset), which preserved completeness while restoring the paralog structure expected of a paleopolyploid genome and outperformed sequence-only clustering. Our results provide a robust, genome-free framework for constructing full-length reference transcriptomes in non-model plant species and demonstrate that pipeline choice must be evaluated per organ rather than assuming one size fits all.
ER -