skip to main content

Title: Integrating Bayesian genomic cline analyses and association mapping of morphological and ecological traits to dissect reproductive isolation and introgression in a Louisiana Iris hybrid zone

Hybrid zones provide unique opportunities to examine reproductive isolation and introgression in nature. We utilized 45,384 single nucleotide polymorphism (SNP) loci to perform association mapping of 14 floral, vegetative and ecological traits that differ betweenIris hexagonaandIris fulva,and to investigate, using a Bayesian genomic cline (BGC) framework, patterns of genomic introgression in a large and phenotypically diverse hybrid zone in southern Louisiana. Many loci of small effect size were consistently found to be associated with phenotypic variation across all traits, and several individual loci were revealed to influence phenotypic variation across multiple traits. Patterns of genomic introgression were quite heterogeneous throughout the Louisiana Iris genome, withI. hexagonaalleles tending to be favoured over those ofI. fulva. Loci that were found to have exceptional patterns of introgression were also found to be significantly associated with phenotypic variation in a small number of morphological traits. However, this was the exception rather than the rule, as most loci that were associated with morphological trait variation were not significantly associated with excess ancestry. These findings provide insights into the complexity of the genomic architecture of phenotypic differences and are a first step towards identifying loci that are associated with both trait variation and reproductive isolation in nature.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Date Published:
Journal Name:
Molecular Ecology
Page Range / eLocation ID:
p. 959-978
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between thealbaandpersonatasubspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437SNPloci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome‐wide divergence. Variation in only one trait—head plumage patterning—was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome‐wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.

    more » « less
  2. Abstract

    Hybrid zones are geographic regions where isolating barriers between divergent populations are challenged by admixture. Identifying factors that facilitate or inhibit hybridization in sympatry can illuminate the processes that maintain those reproductive barriers. We analysed patterns of hybridization and phenotypic variation across two newly discovered hybrid zones between three subspecies of barn swallow (Hirundo rustica). These subspecies differ in ventral coloration and wing length, traits that are targets of sexual and natural selection, respectively, and are associated with genome‐wide differentiation in allopatry. We tested the hypothesis that the degree of divergence in these traits is associated with the extent of hybridization in secondary contact. We applied measures of population structure based on >23,000SNPs to confirm that named subspecies correspond to distinct genomic clusters, and assessed coincidence between geographic clines for ancestry and phenotype. Although gene flow was ongoing across both hybrid zones and pairwiseFSTbetween subspecies was extremely low, we found striking differences in the extent of hybridization. In the more phenotypically differentiated subspecies pair, clines for ancestry, wing length and ventral coloration were steep and coincident, suggestive of strong isolation and, potentially, selection associated with phenotype. In the less phenotypically differentiated pair, gene flow and phenotypic variation occurred over a wide geographic span, indicative of weaker isolation. Traits associated with genome‐wide differentiation in allopatry may thus also contribute to isolation in sympatry. We discuss potentially important additional roles for evolutionary history and ecology in shaping variation in the extent hybridization between closely related pairs of subspecies.

    more » « less
  3. Abstract

    A major goal of speciation research is to reveal the genomic signatures that accompany the speciation process. Genome scans are routinely used to explore genome‐wide variation and identify highly differentiated loci that may contribute to ecological divergence, but they do not incorporate spatial, phenotypic or environmental data that might enhance outlier detection. Geographic cline analysis provides a potential framework for integrating diverse forms of data in a spatially explicit framework, but has not been used to study genome‐wide patterns of divergence. Aided by a first‐draft genome assembly, we combined anFCTscan and geographic cline analysis to characterize patterns of genome‐wide divergence between divergent pollination ecotypes ofMimulus aurantiacus.FCTanalysis of 58 872SNPs generated viaRAD‐seq revealed little ecotypic differentiation (meanFCT = 0.041), although a small number of loci were moderately‐to‐highly diverged. Consistent with our previous results from the geneMaMyb2, which contributes to differences in flower colour, 130 loci have cline shapes that recapitulate the spatial pattern of trait divergence, suggesting that they may reside in or near the genomic regions that contribute to pollinator isolation. In the narrow hybrid zone between the ecotypes, extensive admixture among individuals and low linkage disequilibrium between markers indicate that most outlier loci are scattered throughout the genome, rather than being restricted to one or a few divergent regions. In addition to revealing the genomic consequences of ecological divergence in this system, we discuss how geographic cline analysis is a powerful but under‐utilized framework for studying genome‐wide patterns of divergence.

    more » « less
  4. Abstract

    Phenotypic plasticity allows organisms to change their phenotype in response to shifts in the environment. While a central topic in current discussions of evolutionary potential, a comprehensive understanding of the genetic underpinnings of plasticity is lacking in systems undergoing adaptive diversification. Here, we investigate the genetic basis of phenotypic plasticity in a textbook adaptive radiation, Lake Malawi cichlid fishes. Specifically, we crossed two divergent species to generate an F3hybrid mapping population. At early juvenile stages, hybrid families were split and reared in alternate foraging environments that mimicked benthic/scraping or limnetic/sucking modes of feeding. These alternate treatments produced a variation in morphology that was broadly similar to the major axis of divergence among Malawi cichlids, providing support for the flexible stem theory of adaptive radiation. Next, we found that the genetic architecture of several morphological traits was highly sensitive to the environment. In particular, of 22 significant quantitative trait loci (QTL), only one was shared between the environments. In addition, we identifiedQTLacting across environments with alternate alleles being differentially sensitive to the environment. Thus, our data suggest that while plasticity is largely determined by loci specific to a given environment, it may also be influenced by loci operating across environments. Finally, our mapping data provide evidence for the evolution of plasticity via genetic assimilation at an important regulatory locus,ptch1. In all, our data address long‐standing discussions about the genetic basis and evolution of plasticity. They also underscore the importance of the environment in affecting developmental outcomes, genetic architectures, morphological diversity and evolutionary potential.

    more » « less
  5. Abstract

    The amphicarpic annual legumeAmphicarpaea bracteatais unusual in producing aerial and subterranean cleistogamous flowers that always self‐fertilize and, less commonly, aerial chasmogamous flowers that outcross. Although both morphologic and genetic variants are known in this highly selfing species, debate continues over whether this variation is continuous, reflecting the segregation of standing genetic variation, or discontinuous, reflecting distinct taxa that rarely intercross. We characterizedSNPvariation in 128 individuals in southern Wisconsin to assess within‐ and among‐population variation at 3928SNPs. We also assessed genotype and leaf morphology in an additional 76 individuals to connect phenotypic variation with genetic variation. Genetic variation maps onto three strongly divergent and highly inbred genetic groups showing little relation to site location. Each group has a distinct phenotype, but the divergence of these groups differs from the varietal divisions previously identified based on morphological characters. Like previous authors, we argue that the taxonomy of this species should be revised. Despite extensive sympatry, estimates of among‐group migration rates are low, and hybrid individuals were at low frequency (<2%) in our dataset. Restricted gene flow likely results from high selfing rates and partial reproductive incompatibility as evidenced by the U‐shaped distribution of pairwiseFSTvalues reflecting “islands” of genomic divergence. These islands may be associated with hybrid incompatibility loci that arose in allopatry. The coexistence of lineages within sites may reflect density‐dependent attack by species‐specific strains of pathogenic fungi and/or root‐nodulating bacteria specializing on distinct genotypes.

    more » « less