The Baetic mountains, located in the southern Iberian Peninsula, is a major hotspot of biodiversity in the Mediterranean Basin, constituting one of the most important glacial refugia for vascular plants in Europe. Despite their relatively limited extension, the Baetic Mountains contain almost 50% of the total endemic Erysimum species in the Iberian Peninsula. The broadly distributed Erysimum genus has diversified profusely in the Mediterranean region, with more than a hundred species described in the area, out of a total of c. 200 species included in the genus. We used two plastid DNA regions (ndhF and trnT-L) and one nuclear DNA region (ITS1-5.8S rDNA-ITS2), with 3,556 bp total length, to carry out phylogenetic analysis by Bayesian inference, maximum likelihood and maximum parsimony, in order to explore the evolutionary relationships between the Erysimum species inhabiting these ranges. Analyses of concatenated sequences from the two genomes identified two main clades with no overlap in species composition so that samples from the same species fell within the same major clade. The phylogenetic relationships depicted by those two clades do not give support to the E. nevadense group, previously proposed on taxonomic grounds. In addition, our results indicated recurrent changes in flower colour in the Baetic Erysimum species although, alternatively, reticulate evolution, which is suggested by incongruent position of taxa in the different trees, may have also affected this trait. Las cordilleras Béticas, localizadas en el sudeste de la Península Ibérica, representan una importante zona para la biodiversidad de la cuenca mediterránea, constituyendo uno de los refugios glaciares más destacados de plantas vasculares en Europa. A pesar de su extensión relativamente limitada, las cordilleras Béticas albergan casi el 50% del total de las especies endémicas de Erysimum de la Península Ibérica. Erysimum es un género ampliamente distribuido, que se ha diversificado profusamente en la región mediterránea, con más de un centenar de especies descritas en dicha área, del total de las alrededor de 200 especies incluidas en el género. Usamos dos regiones de ADN plastidial (ndhF y trnT-L) y una región de ADN nuclear (ITS1-5.8S rDNA-ITS2), con una longitud total de 3.556 pb, para llevar a cabo análisis filogenéticos mediante inferencia bayesiana, máxima verosimilitud y máxima parsimonia, con el fin de explorar las relaciones evolutivas entre las especies de Erysimum que habitan en éstas cordilleras. El análisis de secuencias concatenadas de los dos genomas identifica dos clados principales y las muestras de la misma especie aparecen siempre en un clado o en el otro pero no en ambos. Las relaciones filogenéticas que indican esos dos clados no apoyan la circunscripción del grupo taxonómico E. nevadense, previamente definido en dicho género. Además, los análisis indican cambios recurrentes en el color floral de las especies de Erysimum que habitan los Sistemas Béticos aunque la evolución reticulada, que sugiere la posición de algunos ...