RNA-seq meta-analysis and machine learning identify stress-responsive genes and improve genomic prediction in common bean (Phaseolus vulgaris L.) with cross-species application in cowpea (Vigna unguiculata L.)

Common bean (Phaseolus vulgaris L.) is exposed to a broad spectrum of abiotic and biotic stresses that impose severe constraints on productivity, yet the molecular basis of stress tolerance remains poorly resolved, with independent studies yielding inconsistent and incomplete conclusions. To establish a comprehensive picture of the common bean stress transcriptome, we conducted a systematic meta-analysis of publicly available RNA-sequencing datasets spanning abiotic and biotic stress conditions across leaf and root tissues. Integrating statistical meta-analysis with machine-learning approaches, we identified stress-responsive gene sets whose robustness was verified through rigorous statistical approaches including independent dataset validation. Beyond confirming established stress-responsive genes, the machine-learning framework uncovered candidates overlooked by standard significance thresholds in individual studies yet carrying consistent transcriptional signals across studies. Co-expression and protein-protein network analyses further resolved these candidates into functionally coherent modules linked to specific stress-response programs. Notably, ethylene-responsive transcription factors were identified as hub genes in three of four stress-tissue groups, with NAC domain transcription factors emerging as additional hub genes in biotic stress contexts. Importantly, the biological significance of the identified gene sets was validated genomically: marker panels targeting consensus meta-analysis-derived and machine-learning-discovered gene regions improved genomic prediction accuracy for disease resistance traits in common bean and abiotic stress tolerance traits in cowpea relative to a baseline model with equivalent-sized random marker sets. Overall, these findings revealed conserved stress transcriptome signatures in common bean and provided a cross-species, evidence-based framework for prioritizing candidate genes and constructing biologically informed genomic selection tools to advance stress-resilient legume breeding.