Migration through a major Andean ecogeographic disruption as a driver of genotypic and phenotypic diversity in a wild tomato species

The large number of species on our planet arises from the phenotypic variation and reproductive isolation occurring at the population level. In this study, we sought to understand the origins of such population-level variation in defensive acylsugar chemistry and mating systems in Solanum habrochaites – a wild tomato species found in diverse Andean habitats in Ecuador and Peru. Using Restriction-Associated-Digestion Sequencing (RAD-seq) of 50 S. habrochaites accessions, we identified eight population clusters generated via isolation and hybridization dynamics of 4-6 ancestral populations. Estimation of heterozygosity, fixation index, isolation by distance, and migration probabilities, allowed identification of multiple barriers to gene flow leading to the establishment of extant populations. One major barrier is the Amotape-Huancabamba Zone (AHZ) – a geographical feature in the Andes with high endemism, where the mountainous range breaks up into isolated microhabitats. The AHZ was associated with emergence of alleles for novel reproductive and acylsugar phenotypes. These alleles led to the evolution of self-compatibility in the northern populations, where alleles for novel defense-related enzyme variants were also found to be fixed. We identified geographical distance as a major force causing population differentiation in the central/southern part of the range, where S. habrochaites was also inferred to have originated. Findings presented here highlight the role of the diverse ecogeography of Peru and Ecuador in generating new, reproductively isolated populations, and enhance our understanding of the microevolutionary processes that lay a path to speciation.