Convergent evolution, where independent lineages evolve similar traits when adapting to similar habitats, is a common phenomenon and testament to the repeatability of evolutionary processes. Still, non-convergence is also common, and a major question is whether apparently idiosyncratic, lineage-specific evolutionary changes are reflective of chance events inherent to evolutionary processes, or whether they are also influenced by deterministic genetic or ecological factors. To address this question, we quantified the degree of convergence in genome-wide patterns of gene expression across lineages of livebearing fishes (family Poeciliidae) that span 40 million years of evolution and have colonized extreme environments in the form of toxic, hydrogen-sulfide-rich springs. We specifically asked whether the degree of convergence across lineage pairs was related to their phylogenetic relatedness or the ecological similarity of the habitats they inhabit. Using phylogenetic comparative analyses, we showed that the degree of convergence was highly variable across lineage pairs residing in sulfide springs. While closely related lineages did not exhibit higher degrees of convergence than distantly related ones, we uncovered a strong relationship between degree of convergence and ecological similarity. Our results indicate that variation in the degree of convergence is not merely noise associated with evolutionary contingency. Rather, cryptic environmental variation that is frequently ignored when we employ reductionist approaches can significantly contribute to adaptive evolution. This study highlights the importance of multivariate approaches that capture the complexities of both selective regimes and organismal design when assessing the roles of determinism and contingency in evolution.