A Proof-of-Concept Study: Reference Transcriptome Choice as a Hidden Analytical Variable in Breast Cancer RNA-Seq Analysis
- 1Pathomorphology and Clinical Immunology, Poznan University of Medical Sciences, 49 Przybyszewskiego St., 60-356 Poznań, Poland
- 2Department of Genetics, Wroclaw University of Environmental and Life Sciences, 7 Kożuchowska St., 51-631 Wroclaw, Poland
Abstract
Reference genome and transcriptome are fundamental components of RNA-seq analysis and can influence transcript quantification and downstream differential expression results. Although GRCh38 is the current human reference build, GRCh37 remains in use in research and clinical bioinformatics workflows. The extent to which reference transcriptome choice affects transcriptomic inference in breast cancer remains insufficiently characterized. We performed a comparative analysis of publicly available RNA-seq data comprising two MCF7 breast cancer cell-line samples and two healthy breast tissue samples from the European Nucleotide Archive (PRJEB4829). Identical preprocessing, transcript quantification, and differential expression workflows were applied using GRCh37- and GRCh38-based reference transcriptomes. Reference transcriptome-dependent variability was further examined by comparing the same samples quantified under the two reference transcriptomes. Concordance of differentially expressed genes and pathway enrichment was assessed, and sensitivity to reference transcriptome choice was evaluated for genes included in the MammaPrint and Naderi expression signatures and selected hereditary cancer gene panels. The MCF7-versus-normal-tissue comparison identified 15,950 differentially expressed genes using GRCh37 and 16,348 using GRCh38, with 13,899 genes shared between the two analyses. The most significantly enriched Hallmark and KEGG pathways were also concordant between annotations. Principal component analysis showed that the dominant pattern of variation in the MCF7-versus-normal-tissue comparison was associated with biological status. In contrast, when the same MCF7 samples were quantified using different reference transcriptomes, reference transcriptome choice represented a major source of variation. Expression of genes comprising the evaluated prognostic signatures and hereditary cancer panels showed no significant differences in median expression between reference transcriptomes. However, individual cancer-relevant genes, including PTEN, RAD51D, NANOG, and SOX11, demonstrated reference transcriptome-dependent expression differences. In this proof-of-concept analysis, reference transcriptome choice produced measurable differences in RNA-seq quantification and differential expression, while broad biological patterns and the aggregate expression of established gene signatures were relatively stable. Reference transcriptome-dependent variability became particularly apparent when biological heterogeneity was limited, suggesting that reference transcriptome choice can act as a hidden analytical variable in transcriptomic studies. These findings support explicit reporting and version control of reference annotations and highlight the need for validation in larger, independently replicated breast cancer cohorts before drawing clinical conclusions.
Related
- Reference Genome Choice Shapes RNA-Seq Quantification and Downstream Interpretation in Chinese Breast Cancer Patients
- De novo transcriptome of the red alga Gracilariopsis tenuifrons (Gracilariales, Rhodophyta) in response to prolonged temperature stress
- A conserved thermodynamic signature links transcriptome organization to growth state
- De novo transcriptome analysis reveals the effects of abiotic and biotic elicitation on naphthoquinone and flavonoid production in Lawsonia inermis L
- Integrated transcriptome and single-cell RNA sequencing identifies transient receptor potential channel-related biomarkers in ulcerative colitis
BibTeX
@article{zielinski2026proofofconcept,
title = {{A Proof-of-Concept Study: Reference Transcriptome Choice as a Hidden Analytical Variable in Breast Cancer RNA-Seq Analysis}},
author = {Paweł Zieliński and Magdalena Frąszczak and Bartosz Czech and Paula Dobosz and Kinga Adamiak and Patryk Mierzejewski and Rafał Stępień and Magda Mielczarek},
journal = {International Journal of Molecular Sciences},
year = {2026},
month = oct,
doi = {10.3390/ijms27219398},
url = {https://www.mdpi.com/1422-0067/27/21/9398},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - A Proof-of-Concept Study: Reference Transcriptome Choice as a Hidden Analytical Variable in Breast Cancer RNA-Seq Analysis
AU - Paweł Zieliński
AU - Magdalena Frąszczak
AU - Bartosz Czech
AU - Paula Dobosz
AU - Kinga Adamiak
AU - Patryk Mierzejewski
AU - Rafał Stępień
AU - Magda Mielczarek
PY - 2026
DA - 2026/10/22
JO - International Journal of Molecular Sciences
DO - 10.3390/ijms27219398
UR - https://www.mdpi.com/1422-0067/27/21/9398
AB - Reference genome and transcriptome are fundamental components of RNA-seq analysis and can influence transcript quantification and downstream differential expression results. Although GRCh38 is the current human reference build, GRCh37 remains in use in research and clinical bioinformatics workflows. The extent to which reference transcriptome choice affects transcriptomic inference in breast cancer remains insufficiently characterized. We performed a comparative analysis of publicly available RNA-seq data comprising two MCF7 breast cancer cell-line samples and two healthy breast tissue samples from the European Nucleotide Archive (PRJEB4829). Identical preprocessing, transcript quantification, and differential expression workflows were applied using GRCh37- and GRCh38-based reference transcriptomes. Reference transcriptome-dependent variability was further examined by comparing the same samples quantified under the two reference transcriptomes. Concordance of differentially expressed genes and pathway enrichment was assessed, and sensitivity to reference transcriptome choice was evaluated for genes included in the MammaPrint and Naderi expression signatures and selected hereditary cancer gene panels. The MCF7-versus-normal-tissue comparison identified 15,950 differentially expressed genes using GRCh37 and 16,348 using GRCh38, with 13,899 genes shared between the two analyses. The most significantly enriched Hallmark and KEGG pathways were also concordant between annotations. Principal component analysis showed that the dominant pattern of variation in the MCF7-versus-normal-tissue comparison was associated with biological status. In contrast, when the same MCF7 samples were quantified using different reference transcriptomes, reference transcriptome choice represented a major source of variation. Expression of genes comprising the evaluated prognostic signatures and hereditary cancer panels showed no significant differences in median expression between reference transcriptomes. However, individual cancer-relevant genes, including PTEN, RAD51D, NANOG, and SOX11, demonstrated reference transcriptome-dependent expression differences. In this proof-of-concept analysis, reference transcriptome choice produced measurable differences in RNA-seq quantification and differential expression, while broad biological patterns and the aggregate expression of established gene signatures were relatively stable. Reference transcriptome-dependent variability became particularly apparent when biological heterogeneity was limited, suggesting that reference transcriptome choice can act as a hidden analytical variable in transcriptomic studies. These findings support explicit reporting and version control of reference annotations and highlight the need for validation in larger, independently replicated breast cancer cohorts before drawing clinical conclusions.
SN - 1422-0067
ER -