Müllerian mimicry of a quantitative trait despite contrasting levels of genomic divergence and selection

Hybrid zones, where distinct populations meet and interbreed, give insight into how differences between populations are maintained despite gene flow. Studying clines in genetic loci and adaptive traits across hybrid zones is a powerful method for understanding how selection drives differentiation within a single species, but can also be used to compare parallel divergence in different species responding to a common selective pressure. Here, we study parallel divergence of wing colouration in the butterflies Heliconius erato and H. melpomene, which are distantly related Müllerian mimics which show parallel geographic variation in both discrete variation in pigmentation, and quantitative variation in structural colour. Using geographic cline analysis, we show that clines in these traits are positioned in roughly the same geographic region for both species, which is consistent with direct selection for mimicry. However, the width of the clines varies markedly between species. This difference is explained in part by variation in the strength of selection acting on colour traits within each species, but may also be influenced by differences in the dispersal rate and total strength of selection against hybrids between the species. Genotyping‐by‐sequencing also revealed weaker population structure in H. melpomene, suggesting the hybrid zones may have evolved differently in each species, which may also contribute to the patterns of phenotypic divergence in this system. Overall, we conclude that multiple factors are needed to explain patterns of clinal variation within and between these species, although mimicry has probably played a central role.     

Variation in complex mating signals in an “island” hybrid zone between Stenobothrus grasshopper species

Two grasshopper species Stenobothrus rubicundus and S. clavatus were previously shown to meet in a narrow hybrid zone on Mount Tomaros in northern Greece. The species are remarkable for their complex courtship songs accompanied by conspicuous movements of antennae and wings. We analyzed variations in forewing morphology, antenna shape, and courtship song across the hybrid zone using a geographic information system, and we documented three contact zones on Mount Tomaros. All male traits and female wings show abrupt transitions across the contact zones, suggesting that these traits are driven by selection rather than by drift. Male clines in antennae are displaced toward S. clavatus, whereas all clines in wings are displaced toward S. rubicundus. We explain cline discordance as depending on sexual selection via female choice. The high covariance between wings and antennae found in the centers of all contact zones results from high levels of linkage disequilibria among the underlying loci, which in turn more likely results from assortative mating than from selection against hybrids. The covariance is found to be higher in clavatus‐like than rubicundus‐like populations, which implies asymmetric assortative mating in parental‐like sites of the hybrid zone and a movement of the hybrid zone in favor of S. clavatus.     

A fine-scale spatial analysis of the mosaic hybrid zone between Gryllus firmus and Gryllus pennsylvanicus

Abstract.— The pattern of character variation within a hybrid zone, the hybrid zone structure, has been used to infer the processes that maintain hybrid zones. Unfortunately it is difficult to infer process from structure alone because many different processes can produce the same pattern of character variation. Mosaic hybrid zones may be maintained by exogenous selection in a heterogeneous environment and/or endogenous selection against hybrid individuals; habitat preference, premating isolating barriers and/or fertility selection can also contribute. The spatial scale at which a hybrid zone is sampled affects its apparent structure; a hybrid zone may appear clinal at one scale and mosaic at another. Here, we sample the mosaic hybrid zone between two field crickets, Gryllus firmus and G. pennsylvanicus , at a scale that spans the boundaries between individual soil-habitat patches. From our analysis, we find that at fine scales, the mosaic hybrid zone resolves into a set of steep clines across patch boundaries. Both morphological and molecular traits exhibit sharp and generally concordant clines. However, clines for mitochondrial DNA and one anonymous nuclear marker are clearly displaced as a result of current hybridization or past introgression (the "ghost of hybridization past"). Thus, scale is important for the structure of this and probably other hybrid zones. The extremely sharp, concordant clines across patch boundaries indicate that the cricket hybrid zone is undoubtedly structured by selection. However, the detailed mechanisms responsible for the maintenance of the hybrid zone–whether endogenous selection against hybrids, exogenous selection by the environment, and/or behavioral preferences for mates or habitats– remain to be elucidated. Determining these mechanisms will depend on closer inspection of the organisms themselves and their interactions, as is the case for all hybrid zones.     

Distinct sources of gene flow produce contrasting population genetic dynamics at different range boundaries of a Choristoneura budworm

Populations are often exposed to multiple sources of gene flow, but accounts are lacking of the population genetic dynamics that result from these interactions or their effects on local evolution. Using a genomic clines framework applied to 1,195 single nucleotide polymorphisms, we documented genomewide, locus‐specific patterns of introgression between Choristoneura occidentalis biennis spruce budworms and two ecologically divergent relatives, C. o. occidentalis and Choristoneura fumiferana, that it interacts with at alternate boundaries of its range. We observe contrasting hybrid indexes between the two hybrid zones, no overlap in “gene‐flow outliers” (clines showing relatively extreme extents or rates of locus‐specific introgression) and variable linkage disequilibrium among those outliers. At the same time, correlated genomewide rates of introgression between zones suggest the presence of processes common to both boundaries. These findings highlight the contrasting population genetic dynamics that can occur at separate frontiers of a single population, while also suggesting that shared patterns may frequently accompany cases of divergence‐with‐gene‐flow that involve a lineage in common. Our results point to potentially complex evolutionary outcomes for populations experiencing multiple sources of gene flow.     

Comparing clines on molecular and phenotypic traits in hybrid zones: a window on tension zone models

The study of zones of secondary contact provides insight into the maintenance of reproductive isolation. Tension zone theory supplies powerful tools for assessing how dispersal and selection shape hybrid zones. We present a multimodal analysis of phenotypic clines in conjunction with clines at molecular markers in a hybrid zone between Larus glaucescens and Larus occidentalis. We developed a new method to analyze simultaneously clines of quantitative traits and molecular data. Low linkage disequilibrium and the lack of coincidence between clines at six microsatellites, a mitochondrial DNA region, and two phenotypic traits indicated introgression. However, the hypothesis of neutral diffusion was rejected based on evidence that all of the clines were concordant and narrower than expected for neutral clines, indicating some indirect selection. The analysis of phenotypic variance gave evidence of restricted phenotypic introgression and together with the bimodal distribution of phenotypes suggested that disruptive selection is acting across the hybrid zone, especially on the coloration of bare parts. Multimodal analysis of phenotypic clines also highlighted a shift between the peak of intermediates and the cline center, left behind by hybrid zone motion. High-resolution analysis of phenotypes distribution thus proved useful for detecting hybrid zone movement even without temporal data.     

PLANT RESPONSES TO CLIMATE IN THE CAPE FLORISTIC REGION OF SOUTH AFRICA: EVIDENCE FOR ADAPTIVE DIFFERENTIATION IN THE PROTEACEAE

Local adaptation along environmental gradients may drive plant species radiation within the Cape Floristic Region (CFR), yet few studies examine the role of ecologically based divergent selection within CFR clades. In this study, we ask whether populations within the monophyletic white protea clade (Protea section Exsertae, Proteaceae) differ in key functional traits along environmental gradients and whether differences are consistent with local adaptation. Using seven taxa, we measured trait–environment associations and selection gradients across 35 populations of wild adults and their offspring grown in two common gardens. Focal traits were leaf size and shape, specific leaf area (SLA), stomatal density, growth, and photosynthetic rate. Analyses on wild and common garden plants revealed heritable trait differences that were associated with gradients in rainfall seasonality, drought stress, cold stress, and less frequently, soil fertility. Divergent selection between gardens generally matched trait–environment correlations and literature‐based predictions, yet variation in selection regimes among wild populations generally did not. Thus, selection via seedling survival may promote gradient‐wide differences in SLA and leaf area more than does selection via adult fecundity. By focusing on the traits, life stages, and environmental clines that drive divergent selection, our study uniquely demonstrates adaptive differentiation among plant populations in the CFR.     

Natural selection drives leaf divergence in experimental populations of Senecio lautus under natural conditions

相似文献   

Cline coupling and uncoupling in a stickleback hybrid zone

Strong ecological selection on a genetic locus can maintain allele frequency differences between populations in different environments, even in the face of hybridization. When alleles at divergent loci come into tight linkage disequilibrium, selection acts on them as a unit and can significantly reduce gene flow. For populations interbreeding across a hybrid zone, linkage disequilibria between loci can force clines to share the same slopes and centers. However, strong ecological selection on a locus can also pull its cline away from the others, reducing linkage disequilibrium and weakening the barrier to gene flow. We looked for this “cline uncoupling” effect in a hybrid zone between stream resident and anadromous sticklebacks at two genes known to be under divergent natural selection (Eda and ATP1a1) and five morphological traits that repeatedly evolve in freshwater stickleback. These clines were all steep and located together at the top of the estuary, such that we found no evidence for cline uncoupling. However, we did not observe the stepped shape normally associated with steep concordant clines. It thus remains possible that these clines cluster together because their individual selection regimes are identical, but this would be very surprising given their diverse roles in osmoregulation, body armor, and swimming performance.     

Hybridization between Townsend's Dendroica townsendi and black‐throated green warblers D. virens in an avian suture zone

Hybrid zones between species provide natural systems for the study of processes involved in divergence, reproductive isolation and speciation. Townsend's Dendroica townsendi and black‐throated green D. virens warblers are phenotypically and genetically divergent groups that occur in western and eastern North America respectively, with potential for range contact in the Rocky Mountains of British Columbia, where other west–east avian pairs come into contact. Although one potential hybrid (a phenotypic Townsend's warbler with the black‐throated green mitochondrial DNA) has been previously reported, there have been no studies of interactions between the taxa in potential areas of sympatry. To determine whether interbreeding between these species is a regular occurrence we examined variation in individuals across the area of putative range overlap. Analysis of plumage, morphology, and mitochondrial (COI) and nuclear molecular markers (CHD1Z and numt‐Dco1) shows surprisingly extensive hybridization between these species, with at least 38% of individuals in the hybrid zone being either hybrids or backcrosses. Each of the traits displays a sigmoidal cline centred along the eastern slope of the Rocky Mountains (molecular cline centres averaging 50 km east of the crest of the Rockies, ranging from 41 to 56 km). The clines are narrow (average molecular cline width is 60 km, ranging from 40 to 87 km) relative to the dispersal distance of related warbler species, suggesting that selection is maintaining the hybrid zone; we discuss possible sources of selection. Given the narrowness of the zone we recommend the two forms should continue to be treated as separate taxonomic species. Townsend's warblers also form an extensively studied hybrid zone with their more closely related southern relative, the hermit warbler D. occidentalis. The combined system of three discrete forms separated by narrow hybrid zones provides an excellent system for the study of hybridization, reproductive isolation and speciation.     

Differential adaptation to a harsh granite outcrop habitat between sympatric Mimulus species

Understanding which environmental variables and traits underlie adaptation to harsh environments is difficult because many traits evolve simultaneously as populations or species diverge. Here, we investigate the ecological variables and traits that underlie Mimulus laciniatus’ adaptation to granite outcrops compared to its sympatric, mesic‐adapted progenitor, Mimulus guttatus. We use fine‐scale measurements of soil moisture and herbivory to examine differences in selective forces between the species’ habitats, and measure selection on flowering time, flower size, plant height, and leaf shape in a reciprocal transplant using M. laciniatus × M. guttatus F₄ hybrids. We find that differences in drought and herbivory drive survival differences between habitats, that M. laciniatus and M. guttatus are each better adapted to their native habitat, and differential habitat selection on flowering time, plant stature, and leaf shape. Although early flowering time, small stature, and lobed leaf shape underlie plant fitness in M. laciniatus’ seasonally dry environment, increased plant size is advantageous in a competitive mesic environment replete with herbivores like M. guttatus’. Given that we observed divergent selection between habitats in the direction of species differences, we conclude that adaptation to different microhabitats is an important component of reproductive isolation in this sympatric species pair.     

Variation in complex courtship traits across a hybrid zone between grasshopper species of the Chorthippus albomarginatus group

Closely related grasshopper species of the Chorthippus albomarginatus group are notable for their extremely complex courtship songs, accompanied by a visual display. Two species of this group, Ch. albomarginatus and Ch. oschei, were previously shown to hybridize in a wide mosaic hybrid zone in Ukraine and Moldova. In this paper, variation in five courtship song characters, one character of visual display and the number of stridulatory pegs were analysed across the hybrid zone to estimate selection against hybrids and strength of assortative mating. Comparison of cline width and position across the hybrid zone showed concordant and coincident clines in four traits, such as three song characters and one morphological character, and discordant and non‐coincident clines in two other song characters and the character of visual display. Concordance of clines in different characters suggests an equal strength of selection acting on underlying loci. Increase of variance and covariance between phenotypic traits at the cline centre could more likely result from assortative mating than from selection against hybrids. Most pairwise cases showed the highest covariance for the oschei‐like, than for the albomarginatus‐like hybrid populations. This indicates that introgression of the oschei genes into the albomarginatus genome is stronger than vice versa, and may be evidence of the movement of the hybrid zone in favour of Ch. albomarginatus. Analysis of associations between phenotype and local vegetation showed that mosaic structure of the hybrid zone is explained to a great extent by habitat–phenotype associations. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 275–291.     

Clinal variation in a brown lemur (Eulemur spp.) hybrid zone: Combining morphological,genetic and climatic data to examine stability

Studies of hybrid zones can inform our understanding of reproductive isolation and speciation. Two species of brown lemur (Eulemur rufifrons and E. cinereiceps) form an apparently stable hybrid zone in the Andringitra region of south‐eastern Madagascar. The aim of this study was to identify factors that contribute to this stability. We sampled animals at 11 sites along a 90‐km transect through the hybrid zone and examined variation in 26 microsatellites, the D‐loop region of mitochondrial DNA, six pelage and nine morphological traits; we also included samples collected in more distant allopatric sites. Clines in these traits were noncoincident, and there was no increase in either inbreeding coefficients or linkage disequilibrium at the centre of the zone. These results could suggest that the hybrid zone is maintained by weak selection against hybrids, conforming to either the tension zone or geographical selection‐gradient model. However, a closer examination of clines in pelage and microsatellites indicates that these clines are not sigmoid or stepped in shape but instead plateau at their centre. Sites within the hybrid zone also occur in a distinct habitat, characterized by greater seasonality in precipitation and lower seasonality in temperature. Together, these findings suggest that the hybrid zone may follow the bounded superiority model, with exogenous selection favouring hybrids within the transitional zone. These findings are noteworthy, as examples supporting the bounded superiority model are rare and may indicate a process of ecologically driven speciation without geographical isolation.     

Genomic variation across two barn swallow hybrid zones reveals traits associated with divergence in sympatry and allopatry

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,000 SNPs 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 pairwise F_ST between 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.     

EVOLUTIONARY IMPLICATIONS OF DIVERGENT CLINES IN AN AVIAN (MANACUS: AVES) HYBRID ZONE

Abstract A previous study of the hybrid zone in western Panama between white‐collared (Manacus candei) and golden‐collared manakins (M. vitellinus) documented the unidirectional introgression of vitellinus male secondary sexual traits across the zone. Here, we examine the hybrid zone in greater genetic and morphological detail. Statistical comparisons of clines are performed using maximum‐likelihood and nonparametric bootstrap methods. Our results demonstrate that an array of six molecular and two morphometric markers agree in cline position and width. Clines for male collar and belly color are similar in width to the first eight clines, but are shifted in position by at least five cline widths. The result is that birds in intervening populations are genetically and morphometrically very like parental candei, but males have the plumage color of parental vitellinus. Neither neutral diffusion nor nonlinearity of color scales appear to be viable explanations for the large cline shifts. Genetic dominance of vitellinus plumage traits is another potential explanation that will require breeding experiments to test. Sexual selection remains a plausible explanation for the observed introgression of vitellinus color traits in these highly dimorphic, polygynous, lek‐mating birds. Two other clines, including a nondiagnostic isozyme locus, are similar in position to the main cluster of clines, but are broader in width. Thus, introgression at some loci is greater than that detected with diagnostic markers. Assuming that narrow clines are maintained by selection, variation in cline width indicates that selection is not uniform throughout the genome and that diagnostic markers are under more intense selective pressure. The traditional focus on diagnostic markers in studies of hybrid zones may therefore lead to underestimates of average introgression. This effect may be more pronounced in organisms with low levels of genetic divergence between hybridizing taxa.     

Genetic and phenotypic variation across a hybrid zone between ecologically divergent tree squirrels (Tamiasciurus)

A hybrid zone along an environmental gradient should contain a clinal pattern of genetic and phenotypic variation. This occurs because divergent selection in the two parental habitats is typically strong enough to overcome the homogenizing effects of gene flow across the environmental transition. We studied hybridization between two parapatric tree squirrels (Tamiasciurus spp.) across a forest gradient over which the two species vary in coloration, cranial morphology and body size. We sampled 397 individuals at 29 locations across a 600‐km transect to seek genetic evidence for hybridization; upon confirming hybridization, we examined levels of genetic admixture in relation to maintenance of phenotypic divergence despite potentially homogenizing gene flow. Applying population assignment analyses to microsatellite data, we found that Tamiasciurus douglasii and T. hudsonicus form two distinct genetic clusters but also hybridize, mostly within transitional forest habitat. Overall, based on this nuclear analysis, 48% of the specimens were characterized as T. douglasii, 9% as hybrids and 43% as T. hudsonicus. Hybrids appeared to be reproductively viable, as evidenced by the presence of later‐generation hybrid genotypes. Observed clines in ecologically important phenotypic traits—fur coloration and cranial morphology—were sharper than the cline of putatively neutral mtDNA, which suggests that divergent selection may maintain phenotypic distinctiveness. The relatively recent divergence of these two species (probably late Pleistocene), apparent lack of prezygotic isolating mechanisms and geographic coincidence of cline centres for both genetic and phenotypic variation suggest that environmental factors play a large role in maintaining the distinctiveness of these two species across the hybrid zone.     

The effect of marker choice on estimated levels of introgression across an avian (Pipridae: Manacus) hybrid zone

Hybrid zones are often characterized by narrow, coincident clines for diverse traits, suggesting that little introgression occurs across them. However, this pattern may result from a bias in focussing on traits that are diagnostic of parental populations. Such choice of highly differentiated traits may cause us to overlook differential introgression in nondiagnostic traits and to distort our perception of hybrid zones. We tested this hypothesis in an avian hybrid zone by comparing cline structure in two sets of molecular markers: isozyme and restriction fragment length polymorphism markers chosen for differentiation between parental forms, and microsatellite markers chosen for polymorphism. Two cline‐fitting methods showed that cline centre positions of microsatellite alleles were more variable than those of isozyme and restriction fragment length polymorphism markers, and several were significantly shifted from those of the diagnostic markers. Cline widths of microsatellite alleles were also variable and two‐ to eightfold wider than those of the diagnostic markers. These patterns are consistent with the idea that markers chosen for differentiation are more likely to be under purifying selection, and studies focussed on these markers will underestimate overall introgression across hybrid zones. Our results suggest that neutral and positively selected alleles may introgress freely across many hybrid zones without altering perceived boundaries between hybridizing forms.     

Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome‐wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly‐differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs’ allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16–0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species – pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.     

Stable Heliconius butterfly hybrid zones are correlated with a local rainfall peak at the edge of the Amazon basin

Multilocus clines between Müllerian mimetic races of Heliconius butterflies provide a classic example of the maintenance of hybrid zones and their importance in speciation. Concordant hybrid zones in the mimics Heliconius erato and H. melpomene in northern Peru were carefully documented in the 1980s, and this prior work now permits a historical analysis of the movement or stasis of the zones. Previous work predicted that these zones might be moving toward the Andes due to selective asymmetry. Extensive deforestation and climate change might also be expected to affect the positions and widths of the hybrid zones. We show that the positions and shapes of these hybrid zones have instead remained remarkably stable between 1985 and 2012. The stability of this interaction strongly implicates continued selection, rather than neutral mixing following secondary contact. The stability of cline widths and strong linkage disequilibria (gametic correlation coefficients R_max = 0.35–0.56 among unlinked loci) over 25 years suggest that mimetic selection pressures on each color pattern locus have remained approximately constant (s ≈ 0.13–0.40 per locus in both species). Exceptionally high levels of precipitation at the edge of the easternmost Andes may act as a population density trough for butterflies, trapping the hybrid zones at the foot of the mountains, and preventing movement. As such, our results falsify one prediction of the Pleistocene Refugium theory: That the ranges of divergent species or subspecies should be centered on regions characterized by maxima of rainfall, with hybrid zones falling in more arid regions between them.     

Preliminary analysis of a hybrid zone between two subspecies of Zootermopsis nevadensis

For species largely allopatric in distribution, zones of contact provide an opportunity for hybridization, testing grounds for species boundaries, and may result in the formation of a new species. Thus, hybrid zones have the potential to provide important insights into speciation. In this study, we performed a preliminary analysis of a hybrid zone between two subspecies of the dampwood termite, Zootermopsis nevadensis (Z. n. nuttingi and Z. n. nevadensis) near Bartle, CA, using 12 microsatellite loci and a mitochondrial gene. Fifty-seven colonies collected in 36 locations were analyzed. The analysis of genetic markers revealed a large hybrid zone approximately 104 km in width. Although stepped clines best explained the data, we are unable to rule out the existence of a mosaic hybrid zone. We inferred weak selection (s* < 0.05%) against hybrids, but the data also suggested the existence of a barrier to gene flow from Z. n. nevadensis to Z. n. nuttingi, but not in the other direction. Given the large zone of contact, extensive sampling is needed to obtain a more complete characterization of the hybrid zone. However, the results of this study suggest that despite the accumulation of phenotypic and genetic differences and time since divergence (~2 million years), Z. n. nuttingi and Z. n. nevadensis are capable of extensive hybridization.