Parallel introgression, not recurrent emergence, explains apparent elevational ecotypes of polyploid Himalayan snowtrout

The recurrence of similar evolutionary patterns within different habitats often reflects parallel selective pressures acting upon either standing or independently occurring genetic variation to produce a convergence of phenotypes. This interpretation (i.e. parallel divergences within adjacent streams) has been hypothesized for drainage-specific morphological ecotypes observed in polyploid snowtrout (Cyprinidae: Schizothorax). However, parallel patterns of differential introgression during secondary contact is a viable alternative hypothesis. Here, we utilized ddRADseq (N=35,319 de-novo and N=10,884 transcriptome-aligned SNPs), as derived from Nepali/Bhutanese samples (N=48 each), to test these competing hypotheses. We first employed genome-wide allelic depths to derive appropriate ploidy models, then a Bayesian approach to yield genotypes statistically consistent under the inferred expectations. Elevational ecotypes were consistent in geometric-morphometric space, but with phylogenetic relationships at the drainage-level, sustaining an hypothesis of independent emergence. However, partitioned analyses of phylogeny and admixture identified subsets of loci under selection that retained genealogical concordance with morphology, suggesting instead that apparent patterns of morphological/phylogenetic discordance are driven by widespread genomic homogenization. Here, admixture occurring in secondary contact effectively masks previous isolation. Our results underscore two salient factors:1) Morphological adaptations are retained despite hybridization; and 2) The degree of admixture varies across tributaries, presumably concomitant with underlying environmental or anthropogenic factors.