Ecology predicts parapatric distributions in two closely related Antirrhinum majus subspecies

Using a species distribution model, we reconstructed the environmental niches of Antirrhinum majus pseudomajus and Antirrhinum majus striatum, two closely related species with parapatric distributions. We tested whether retention of ancestral environmental niche (i.e. niche conservatism) or adaptation to different ecological conditions (i.e. niche divergence) could explain the maintenance of their non-overlapping geographic ranges. We found that the environmental niche of A. m. pseudomajus is almost twice as large as that of A. m. striatum, with substantial overlap indicating that A. m. pseudomajus and A. m. striatum should co-occur frequently within the geographic range of A. m. striatum. By analysing contact zones where both subspecies are geographically close, we found that the presence of one subspecies instead of the other was significantly influenced by particular combinations of climatic factors. Since independent genetic evidence indicates that the two subspecies have experienced phases of range overlap at or near contact zones over the course of their evolutionary history, we propose that ecological niche displacement might be an important factor in explaining the absence of current range overlap between A. majus subspecies.     

Introgression at differentially aged hybrid zones in red-tailed chipmunks

Hybrid zones allow us to investigate the maintenance and the break down of reproductive isolation; they are a window into the speciation process. Tamias ruficaudus (red-tailed chipmunk) has a roughly ring-like distribution in the Inland Northwest and includes two morphologically well-differentiated subspecies, T. r. ruficaudus (in the eastern portion of its range) and T. r. simulans (in the western portion). These taxa meet at two contact zones: the Lochsa River in Idaho and 200 km to the north, near Whitefish, Montana. The Lochsa Zone is encompassed within the Clearwater River Drainage, which has been proposed as a glacial refugium for many taxa throughout the Pleistocene, whereas the Whitefish Zone was under the Cordilleran ice sheet during the most recent glacial maxima approximately 10,000 years ago. Mitochondrial DNA introgression has been documented at both contact zones, yet the subspecies remain significantly distinct with respect to bacular morphology and no intermediate morphologies have ever been reported. Here, we elucidate differentiation and introgression of the nuclear genome using ten microsatellite loci and compare findings to previously described mitochondrial DNA haplotype distribution and introgression. We found significant substructure in the nuclear data; each subspecies is divided into at least two genetically distinct demes. At the Lochsa contact zone, individuals restricted to the mtDNA zone of introgression form a distinct deme at microsatellite loci whereas in the younger, Whitefish contact zone, there is no hybrid-zone specific group. The genetic distances of the demes within these two subspecies indicate recent northward expansion.     

Cytonuclear discordance across a leopard frog contact zone

We conducted a genetic analysis of the extent of hybridization within and outside a contact zone between two North American leopard frogs, Rana blairi and Rana pipiens by comparing distribution patterns of mitochondrial and nuclear haplotypes. The contact zone, located between South Dakota, Nebraska, and Iowa, USA, was previously examined using morphological data in the 1970s, leading to the conclusion that hybridization was rare between R. pipiens and R. blairi. Our genetic analysis of 51 populations (611 samples) shows strong cytonuclear discordance. Mitochondrial-haplotype distribution matches the same pattern as the documented species spatial distributions based on morphology. However, the geographic distribution of the nuclear haplotypes reveals asymmetrical swamping of the R. pipiens nuclear haplotypes by R. blairi haplotypes. Phylogenetic analyses of both mitochondrial and nuclear markers provide strong evidence for the presence of R. blairi-R. pipiens introgression for the nuclear marker. A pattern of mitochondrial isolation with nuclear introgression is extremely unusual, and predicted to occur much less often than the reverse.     

Phylogeography of white-spotted Charr (Salvelinus leucomaenis) inferred from mitochondrial DNA sequences

The white-spotted charr (Salvelinus leucomaenis) is a coldwater-adapted fish distributed in far-eastern Asia. To assess phylogeographic patterns of this species over most of its range in the Japanese archipelago and Sakhalin Island, Russia, we examined nucleotide sequences of the mitochondrial DNA (mtDNA) cytochrome b region (557 bp) in 141 individuals from 50 populations. A total of 33 (5.5%) nucleotide positions were polymorphic and defined 29 haplotypes. Phylogenetic analysis assigned the observed haplotypes to four main clades, which were characterized by the idiosyncrasies and discontinuity of geographic distributions. The nested clade analyses revealed that the geographical distribution patterns of some haplotypes and clades were explained by historical event such as past fragmentation. Although substantial genetic differentiation was found among the four main clades, their geographic distributions overlapped extensively in several regions. Since white-spotted charr can potentially use both freshwater and marine environments, coexistence among different lineages can be attributed to secondary contact through range expansion by migratory individuals during multiple glacial periods after interglacial isolation. Finally, our data demonstrate that the current subspecies designation does not reflect the phylogeography of this species based on mtDNA analysis. Hierarchical analysis (AMOVA) also showed that genetic variation was far more pronounced within subspecies than among subspecies (i.e., among discrete regions). These results suggest that each population, rather than each subspecies, must be treated as an evolutionarily significant unit.     

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.     

Population genetic structure and hybridization patterns in the Mediterranean endemics Phlomis lychnitis and P. crinita (Lamiaceae)

BACKGROUND AND AIMS: The historical influence of gene flow and genetic drift after the last glacial phase of the Quaternary Period is reflected in current levels of genetic diversity and population structure of plant species. Moreover, hybridization after secondary contact might also affect population genetic diversity and structure. An assessment was made of the genetic variation and hybrid zone structure in Iberian populations of the Mediterranean Phlomis lychnitis and P. crinita, for which phylogenetic relationships are controversial, and hybridization and introgression are common. METHODS: Allozyme variation at 13 loci was analysed in 1723 individual plants sampled from 35 natural locations of P. lychnitis, P. crinita subsp. malacitana and P. crinita subsp. crinita in southern and eastern Spain. Standard genetic diversity parameters were calculated and patterns of genetic structure in each taxon were tested to fit the equilibrium between gene flow and genetic drift. Individual multilocus genotypes were subjected to Bayesian clustering analysis to estimate hybridization and introgression rates for both geographic regions. KEY RESULTS: Contrasting patterns in the distribution of genetic variation among the three taxa were found. Phlomis lychnitis showed no significant inbreeding, low genetic differentiation among populations and no evidence of isolation by distance. Phlomis crinita subsp. malacitana and P. crinita subsp. crinita showed high levels of genetic structure consistent with a pattern of gene flow-drift equilibrium. Higher instances of hybridization and introgression were detected in locations from southern Spain compared with locations from eastern Spain, matching unimodal and bimodal hybrid zones, respectively. CONCLUSIONS: High instances of historical gene flow, range expansion and altitudinal movement during the Quaternary Period, and lineage sorting can explain the diversity of patterns observed. The results suggest that P. lychnitis is the most differentiated lineage in the group; however, the relationship between the three taxa remains unclear.     

Assessing multilocus introgression patterns: a case study on the mouse X chromosome in central Europe

Multilocus hybrid zone (HZ) studies predate genomics by decades. The power of early methods is becoming apparent and now large datasets are commonplace. Relating introgression along a chromosome to evolutionary process is challenging: although reduced introgression regions may indicate speciation genes, this pattern may be obscured by asymmetric introgression of linked invasive genes. Further, HZ movement may form salients and leave islands in its wake. Barton's concordance was proposed 24 years ago for assessing introgression where geographic patterns are complex. The geographic axis of introgression is replaced with the hybrid index. We compare this, a recently proposed genomic clines approach, and two-dimensional (2D) geographic analyses, for 24 X chromosome loci of 2873 mice from the central-European house mouse HZ. In 2D, 14 loci show linear contact, seven precisely matching previous studies. Four show introgression islands to the east of the zone, suggesting past westward zone movement, two show westward salients. Barton's concordance both recovers and refines this information. A region of reduced introgression on the central X is supported, despite X centromere-proximal male-biased westward introgression matching a westward 2D geographic salient. Genomic clines results are consistent regarding introgression asymmetries, but otherwise more difficult to interpret. Evidence for genetic conflict is discussed.     

Geographic structure of genetic and phenotypic variation in the hybrid zone between Quercus affinis and Q. laurina in Mexico

Analyzing the structure of hybrid zones is important for inferring their origin, dynamics and evolutionary significance. We examined the geographic structure of phenotypic and genetic variation in the contact zone between two Mexican red oaks, Quercus affinis and Q. laurina. A total of 105 individuals from seven populations were sampled along a 600‐km latitudinal gradient representing the distribution area of the two species and their contact zone. Individuals were genotyped for nine nuclear and four chloroplast DNA microsatellite loci (ncSSR and cpSSR, respectively), and characterized for several leaf and acorn traits. The cpSSR data revealed extensive haplotype sharing among populations of the two species, while a Bayesian assignment analysis based on ncSSRs identified two main genetic groups, each corresponding to one of the species, and two populations in the contact zone showing evidence of admixture. The proportion of genetic ancestry in the populations was strongly associated with latitude and showed a pattern of variation with the shape of a narrow sigmoidal cline. The variation in three of the seven phenotypic traits was partially congruent with molecular variation, while the other traits did not conform to a geographic cline but instead were correlated with environmental variables. In conclusion, the hybrid zone between the two oak species has some of the characteristics of a tension zone, but heterogeneous variation across traits suggests differential introgression and the action of extrinsic selection.     

Sympatric serpentine endemic Monardella (Lamiaceae) species maintain habitat differences despite hybridization

Ecological differentiation and genetic isolation are thought to be critical in facilitating coexistence between related species, but the relative importance of these phenomena and the interactions between them are not well understood. Here, we examine divergence in abiotic habitat affinity and the extent of hybridization and introgression between two rare species of Monardella (Lamiaceae) that are both restricted to the same serpentine soil exposure in California. Although broadly sympatric, they are found in microhabitats that differ consistently in soil chemistry, slope, rockiness and vegetation. We identify one active hybrid zone at a site with intermediate soil and above‐ground characteristics, and we document admixture patterns indicative of extensive and asymmetric introgression from one species into the other. We find that genetic distance among heterospecific populations is related to geographic distance, such that the extent of apparent introgression is partly explained by the spatial proximity to the hybrid zone. Our work shows that plant species can maintain morphological and ecological integrity in the face of weak genetic isolation, intermediate habitats can facilitate the establishment of hybrids, and that the degree of apparent introgression a population experiences is related to its geographic location rather than its local habitat characteristics.     

THE VARIABLE GENOMIC ARCHITECTURE OF ISOLATION BETWEEN HYBRIDIZING SPECIES OF HOUSE MICE

Studies of the genetics of hybrid zones can provide insight into the genomic architecture of species boundaries. By examining patterns of introgression of multiple loci across a hybrid zone, it may be possible to identify regions of the genome that have experienced selection. Here, we present a comparison of introgression in two replicate transects through the house mouse hybrid zone through central Europe, using data from 41 single nucleotide markers. Using both genomic and geographic clines, we found many differences in patterns of introgression between the two transects, as well as some similarities. We found that many loci may have experienced the effects of selection at linked sites, including selection against hybrid genotypes, as well as positive selection in the form of genotypes introgressed into a foreign genetic background. We also found many positive associations of conspecific alleles among unlinked markers, which could be caused by epistatic interactions. Different patterns of introgression in the two transects highlight the challenge of using hybrid zones to identify genes underlying isolation and raise the possibility that the genetic basis of isolation between these species may be dependent on the local population genetic make-up or the local ecological setting.     

Sharp acoustic boundaries across an altitudinal avian hybrid zone despite asymmetric introgression

Birdsong is a sexually selected trait that could play an important evolutionary role when related taxa come into secondary contact. Many songbird species, however, learn their songs through copying one or more tutors, which complicates the evolutionary outcome of such contact. Two subspecies of a presumed vocal learner, the grey‐breasted wood‐wren (Henicorhina leucophrys), replace each other altitudinally across the western slope of the Ecuadorian Andes. These subspecies are morphologically very similar, but show striking differences in their song. We examined variation in acoustic traits and genetic composition across the altitudinal range covered by both subspecies and between two allopatric populations. The acoustic boundary between the subspecies was found to be highly abrupt across a narrow elevational range with virtually no evidence of song convergence. Mixed singing and use of hetero‐subspecific song occurred in the contact zone and was biased towards the use of leucophrys song types. Hetero‐subspecific song copying by hilaris and not by leucophrys reflected a previously found asymmetric pattern of response to song playback. Using amplified fragment length polymorphisms (AFLP) markers, we detected hybridization in the contact zone and asymmetric introgression in parapatric populations, with more leucophrys alleles present in hilaris populations than vice versa. This pattern may be a trail of introgression due to upslope displacement of leucophrys by hilaris. Our data suggest that song learning may impact speciation and hybridization in contrasting ways at different spatial scales: although learning may speed up population divergence in songs, thereby enhancing assortative mating and reducing gene flow, it may at a local level also lead to the copying of heterospecific songs, therefore allowing some level of hybridization and introgression.     

Genetic structure, introgression, and a narrow hybrid zone between northern and California spotted owls (Strix occidentalis)

The northern spotted owl (Strix occidentalis caurina) is a threatened subspecies and the California spotted owl (Strix occidentalis occidentalis) is a subspecies of special concern in the western United States. Concern for their continued viability has arisen because of habitat loss caused by timber harvesting. The taxonomic status of the northern subspecies has been the subject of continuing controversy. We investigated the phylogeographical and population genetic structure of northern and California spotted owls with special reference to their region of contact. Mitochondrial DNA (mtDNA) control region sequences confirmed the existence of two well-differentiated lineages connected by a narrow hybrid zone in a region of low population density in north central California. Maximum-likelihood estimates indicated bidirectional gene flow between the lineages but limited introgression outside the region of contact. The lengths of both the mtDNA hybrid zone and the reduced density patch were similar and slightly exceeded estimates of natal dispersal distances. This suggests that the two subspecies were in secondary contact in a hybrid zone trapped by a population density trough. Consequently, the zone of interaction is expected to be geographically stable. We discovered a third, rare clade of haplotypes, which we interpreted to be a result of incomplete lineage sorting; those haplotypes result in a paraphyletic northern spotted owl with respect to the California spotted owl. A congeneric species, the barred owl (Strix varia), occasionally hybridizes with spotted owls; our results indicated an upper bound for the frequency of barred owl mtDNA haplotypes in northern spotted owl populations of 3%.     

Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula

The role of the Y chromosome in speciation is unclear. Hybrid zones provide natural arenas for studying speciation, as differential introgression of markers may reveal selection acting against incompatibilities. Two subspecies of the European rabbit (Oryctolagus cuniculus) form a hybrid zone in the Iberian Peninsula. Previous work on mitochondrial DNA (mtDNA), Y- and X-linked loci revealed the existence of two divergent lineages in the rabbit genome and that these lineages are largely subspecies-specific for mtDNA and two X-linked loci. Here we investigated the geographic distribution of the two Y chromosome lineages by genotyping two diagnostic single nucleotide polymorphisms in a sample of 353 male rabbits representing both subspecies, and found that Y chromosome lineages are also largely subspecies-specific. We then sequenced three autosomal loci and discovered considerable variation in levels of differentiation at these loci. Finally, we compared estimates of population differentiation between rabbit subspecies at 26 markers and found a surprising bimodal distribution of F(ST)values. The vast majority of loci showed little or no differentiation between rabbit subspecies while a few loci, including the SRY gene, showed little or no introgression across the hybrid zone. Estimates of population differentiation for the Y chromosome were surprisingly high given that there is male-biased dispersal in rabbits. Taken together, these data indicate that there is a clear dichotomy in the rabbit genome and that some loci remain highly differentiated despite extensive gene flow following secondary contact.     

Low intraspecific variation for genomic isolation between hybridizing sunflower species

Barriers to gene flow between species result from selection against foreign linkage blocks in hybrids. When the geographic ranges of taxa meet at multiple locations, the opportunity exists for variation in the genetic architecture of isolating barriers. Hybrid zones between two sunflower species (Helianthus annuus and H. petiolaris) in Nebraska and California exhibited remarkably similar patterns of introgression of mapped molecular markers. Congruence among hybrid zones may result from limited intraspecific variation at loci contributing to isolation and from similar selective effects of alleles in the heterospecific genetic background. The observed consistency of introgression patterns across distantly separated hybrid zones suggests that intrinsic forces predominate in determining hybrid zone dynamics and boundaries between these sunflower species.     

Phylogeography of speciation: allopatric divergence and secondary contact between outcrossing and selfing Clarkia

The origins of hybrid zones between parapatric taxa have been of particular interest for understanding the evolution of reproductive isolation and the geographic context of species divergence. One challenge has been to distinguish between allopatric divergence (followed by secondary contact) versus primary intergradation (parapatric speciation) as alternative divergence histories. Here, we use complementary phylogeographic and population genetic analyses to investigate the recent divergence of two subspecies of Clarkia xantiana and the formation of a hybrid zone within the narrow region of sympatry. We tested alternative phylogeographic models of divergence using approximate Bayesian computation (ABC) and found strong support for a secondary contact model and little support for a model allowing for gene flow throughout the divergence process (i.e. primary intergradation). Two independent methods for inferring the ancestral geography of each subspecies, one based on probabilistic character state reconstructions and the other on palaeo-distribution modelling, also support a model of divergence in allopatry and range expansion leading to secondary contact. The membership of individuals to genetic clusters suggests geographic substructure within each taxon where allopatric and sympatric samples are primarily found in separate clusters. We also observed coincidence and concordance of genetic clines across three types of molecular markers, which suggests that there is a strong barrier to gene flow. Taken together, our results provide evidence for allopatric divergence followed by range expansion leading to secondary contact. The location of refugial populations and the directionality of range expansion are consistent with expectations based on climate change since the last glacial maximum. Our approach also illustrates the utility of combining phylogeographic hypothesis testing with species distribution modelling and fine-scale population genetic analyses for inferring the geography of the divergence process.     

SIGNIFICANT ROLE FOR HISTORICAL EFFECTS IN THE EVOLUTION OF REPRODUCTIVE ISOLATION: EVIDENCE FROM PATTERNS OF INTROGRESSION BETWEEN THE CYPRINID FISHES,LUXILUS CORNUTUS AND LUXILUS CHRYSOCEPHALUS

Samples of Luxilus cornutus, Luxilus chrysocephalus, and their hybrids were collected along hypothesized routes of dispersal from Pleistocene refugia to examine the significance of geographic variation in patterns of introgression between these species. Patterns of allozyme and mitochondrial DNA (mtDNA) variation were generally consistent with those from previous studies. Tests of Hardy-Weinberg equilibrium revealed significant deficiencies of heterozygotes in all samples, indicating some form of reproductive isolation. Mitochondrial DNAs of each species were not equally represented in F₁ hybrids; however, this bias was eliminated when the two largest samples were excluded from the analysis. Backcross hybrids exhibited biased mtDNA introgression, as samples from Lake Erie (eastern) and Lake Michigan (western) drainages showed significant excesses of mtDNAs from L. chrysocephalus and L. cornutus, respectively, relative to frequencies of diagnostic allozyme markers. The extent and direction of allozyme and mtDNA introgression was quantified by calculating isolation index values from morphologically “pure” individuals of each species from each locality. Analysis of variance of these measures identified limited introgression of allozyme variants with no geographic pattern, but significant differences in direction of mtDNA introgression between drainages (i.e., postglacial dispersal route). Association between patterns of mtDNA introgression and dispersal route across the latitudinal width of the contact zone is best explained by genetic divergence during past isolation of ancestral populations from these drainages. These results identify a significant role for historical effects in the evolution of reproductive isolation and the process of speciation.     

Micro-spatial genetic structure in song sparrows (<Emphasis Type="Italic">Melospiza melodia</Emphasis>)

The spatial genetic structure of populations is strongly influenced by current and historical patterns of gene flow and drift, which in the simplest case, is limited by geographic distance. We examined the microspatial genetic structure within 33 populations of song sparrows (Melospiza melodia) which included eight subspecies located across coastal areas in southern British Columbia (BC) and California. We also examined the effect of water barriers and local density estimates on genetic structuring. Across both regions, positive genetic structure was detectable at distances of less than 10 km. Genetic divergence was highest in Californian subspecies, perhaps due to reduced gene flow across sub-specific contact zones. In BC, populations distributed across islands displayed greater genetic structuring over similar spatial scales than those across mainland sites, supporting the prediction that water barriers reduce gene flow in this species. Our results confirm both the expectation for fine-scale genetic structure in these generally sedentary subspecies, and the role of landscape features in generating geographic variation in genetic structure.     

Distortion of symmetrical introgression in a hybrid zone: evidence for locus-specific selection and uni-directional range expansion

The fate of species integrity upon natural hybridization depends on the interaction between selection and dispersal. The relative significance of these processes may be studied in the initial phase of contact before selection and gene flow reach equilibrium. Here we study a hybrid zone of two salamander species, Lyciasalamandra antalyana and Lyciasalamandra billae, at the initial phase of hybridization. We quantify the degree and mode of introgression using nuclear and mtDNA markers. The hybrid zone can be characterized as an abrupt transition zone, the central hybrid zone being only c. 400 m, but introgressed genes were traced up to 3 km. Introgression was traced in both sexes but gene flow may be slightly male-biased. Indirect evidence suggests that hybrid males are less viable than females. Introgression occurred at two levels: (1) locus-specific selection led to different allelic introgression patterns independent of species, while (2) asymmetrical species-level introgression occurred predominately from L. antalyana to L. billae due to range expansion of the former. This indicates that foreign genes can be incorporated into novel genomic environments, which in turn may contribute to the great diversity of morphological variants in Lyciasalamandra.     

Asymmetric admixture and morphological variability at a suture zone: parapatric burbot subspecies (Pisces) in the Mackenzie River basin,Canada

The burbot (Lota lota L.) is a northern freshwater fish with a circumpolar distribution. Two subspecies diverged due to isolation during glacial maxima: Lota lota lota is the Eurasian-Beringian subspecies and Lota lota maculosa the North American subspecies. We sampled burbot from the Great Slave Lake and Mackenzie River area, Canada, the only known contact zone of these two lineages. Using molecular methods (microsatellite loci and mtDNA sequence) we found that the subspecies’ distributions abut in the Mackenzie River delta, with L. l. lota in the lower delta and L. l. maculosa in all upstream rivers and lakes. Admixture between subspecies was minimal, decreased with increasing geographic distance, and was asymmetrical: mitochondrial and nuclear genetic introgression was from L. l. lota into L. l. maculosa but not the reverse. Within subspecies, there was low inter-population genetic differentiation, no isolation-by-distance, and no evidence for sex-biased dispersal. We did not identify a difference in body length between subspecies per se, though mean lengths differed among localities. Thus, genetic data demonstrate that burbot subspecies are reproductively isolated though the extent to which morphologically variability relates to local versus subspecific variation remains unclear.