Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y‐chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least‐divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi‐absence of Y‐chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non‐neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.     

Phylogeography and Population Structure of Glossina fuscipes fuscipes in Uganda: Implications for Control of Tsetse

Background

Glossina fuscipes fuscipes, a riverine species of tsetse, is the main vector of both human and animal trypanosomiasis in Uganda. Successful implementation of vector control will require establishing an appropriate geographical scale for these activities. Population genetics can help to resolve this issue by characterizing the extent of linkage among apparently isolated groups of tsetse.

Methodology/Principal Findings

We conducted genetic analyses on mitochondrial and microsatellite data accumulated from approximately 1000 individual tsetse captured in Uganda and neighboring regions of Kenya and Sudan. Phylogeographic analyses suggested that the largest scale genetic structure in G. f. fuscipes arose from an historical event that divided two divergent mitochondrial lineages. These lineages are currently partitioned to northern and southern Uganda and co-occur only in a narrow zone of contact extending across central Uganda. Bayesian assignment tests, which provided evidence for admixture between northern and southern flies at the zone of contact and evidence for northerly gene flow across the zone of contact, indicated that this structure may be impermanent. On the other hand, microsatellite structure within the southern lineage indicated that gene flow is currently limited between populations in western and southeastern Uganda. Within regions, the average F_ST between populations separated by less than 100 km was less than ∼0.1. Significant tests of isolation by distance suggested that gene flow is ongoing between neighboring populations and that island populations are not uniformly more isolated than mainland populations.

Conclusions/Significance

Despite the presence of population structure arising from historical colonization events, our results have revealed strong signals of current gene flow within regions that should be accounted for when planning tsetse control in Uganda. Populations in southeastern Uganda appeared to receive little gene flow from populations in western or northern Uganda, supporting the feasibility of area wide control in the Lake Victoria region by the Pan African Tsetse and Trypanosomiasis Eradication Campaign.     

The tsetse fly Glossina fuscipes fuscipes (Diptera: Glossina) harbours a surprising diversity of bacteria other than symbionts

Three different bacterial species are regularly described from tsetse flies. However, no broad screens have been performed to investigate the existence of other bacteria in this medically and agriculturally important vector insect. Utilising both culture dependent and independent methods we show that Kenyan populations of Glossina fuscipes fuscipes harbour a surprising diversity of bacteria. Bacteria were isolated from 72% of flies with 23 different bacterial species identified. The Firmicutes phylum dominated with 16 species of which seven belong to the genus Bacillus. The tsetse fly primary symbiont, Wigglesworthia glossinidia, was identified by the culture independent pathway. However, neither the secondary symbiont Sodalis nor Wolbachia was detected with either of the methods used. Two other bacterial species were identified with the DNA based method, Bacillus subtilis and Serratia marcescens. Further studies are needed to determine how tsetse flies, which only ever feed on vertebrate blood, pick up bacteria and to investigate the possible impact of these bacteria on Glossina longevity and vector competence.     

Hybrid swarm between divergent lineages of mule deer (Odocoileus hemionus)

Studies of hybrid zones have revealed an array of evolutionary outcomes, yet the underlying structure is typically characterized as one of three types: a hybrid zone, a hybrid swarm or a hybrid taxon. Our primary objective was to determine which of these three structures best characterizes a zone of hybridization between two divergent lineages of mule deer (Odocoileus hemionus), mule deer and black-tailed deer. These lineages are morphologically, ecologically and genetically distinct, yet hybridize readily along a zone of secondary contact between the east and west slopes of the Cascade Mountains (Washington and Oregon, USA). Using microsatellite and mitochondrial DNA, we found clear evidence for extensive hybridization and introgression between lineages, with varying degrees of admixture across the zone of contact. The pattern of hybridization in this region closely resembles a hybrid swarm; based on data from 10 microsatellite loci, we detected hybrids that extend well beyond the F1 generation, did not detect linkage disequilibrium at the centre of the zone and found that genotypes were associated randomly within the zone of contact. Introgression was characterized as bidirectional and symmetric, which is surprising given that the zone of contact occurs along a sharp ecotone and that lineages are characterized by large differences in body size (a key component of mating success). Regardless of the underlying mechanisms promoting hybrid swarm maintenance, it is clear that the persistence of a hybrid swarm presents unique challenges for management in this region.     

Geographic variation in hybridization across a reinforcement contact zone of chorus frogs (Pseudacris)

Reinforcement contact zones, which are secondary contact zones where species are diverging in reproductive behaviors due to selection against hybridization, represent natural laboratories for studying speciation‐in‐action. Here, we examined replicate localities across the entire reinforcement contact zone between North American chorus frogs Pseudacris feriarum and P. nigrita to investigate geographic variation in hybridization frequencies and to assess whether reinforcement may have contributed to increased genetic divergence within species. Previous work indicated these species have undergone reproductive character displacement (RCD) in male acoustic signals and female preferences due to reinforcement. We also examined acoustic signal variation across the contact zone to assess whether signal characteristics reliably predict hybrid index and to elucidate whether the degree of RCD predicts hybridization rate. Using microsatellites, mitochondrial sequences, and acoustic signal information from >1,000 individuals across >50 localities and ten sympatric focal regions, we demonstrate: (1) hybridization occurs and (2) varies substantially across the geographic range of the contact zone, (3) hybridization is asymmetric and in the direction predicted from observed patterns of asymmetric RCD, (4) in one species, genetic distance is higher between conspecific localities where one or both have been reinforced than between nonreinforced localities, after controlling for geographic distance, (5) acoustic signal characters strongly predict hybrid index, and (6) the degree of RCD does not strongly predict admixture levels. By showing that hybridization occurs in all sympatric localities, this study provides the fifth and final line of evidence that reproductive character displacement is due to reinforcement in the chorus frog contact zone. Furthermore, this work suggests that the dual action of cascade reinforcement and partial geographic isolation is promoting genetic diversification within one of the reinforced species.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

Inferring the degree of incipient speciation in secondary contact zones of closely related lineages of Palearctic green toads (Bufo viridis subgroup)

Reproductive isolation between lineages is expected to accumulate with divergence time, but the time taken to speciate may strongly vary between different groups of organisms. In anuran amphibians, laboratory crosses can still produce viable hybrid offspring >20 My after separation, but the speed of speciation in closely related anuran lineages under natural conditions is poorly studied. Palearctic green toads (Bufo viridis subgroup) offer an excellent system to address this question, comprising several lineages that arose at different times and form secondary contact zones. Using mitochondrial and nuclear markers, we previously demonstrated that in Sicily, B. siculus and B. balearicus developed advanced reproductive isolation after Plio-Pleistocene divergence (2.6 My, 3.3–1.9), with limited historic mtDNA introgression, scarce nuclear admixture, but low, if any, current gene flow. Here, we study genetic interactions between younger lineages of early Pleistocene divergence (1.9 My, 2.5–1.3) in northeastern Italy (B. balearicus, B. viridis). We find significantly more, asymmetric nuclear and wider, differential mtDNA introgression. The population structure seems to be molded by geographic distance and barriers (rivers), much more than by intrinsic genomic incompatibilities. These differences of hybridization between zones may be partly explained by differences in the duration of previous isolation. Scattered research on other anurans suggests that wide hybrid zones with strong introgression may develop when secondary contacts occur <2 My after divergence, whereas narrower zones with restricted gene flow form when divergence exceeds 3 My. Our study strengthens support for this rule of thumb by comparing lineages with different divergence times within the same radiation.     

Postmating reproductive barriers contribute to the incipient sexual isolation of the United States and Caribbean Drosophila melanogaster

The nascent stages of speciation start with the emergence of sexual isolation. Understanding the influence of reproductive barriers in this evolutionary process is an ongoing effort. We present a study of Drosophila melanogaster admixed populations from the southeast United States and the Caribbean islands known to be a secondary contact zone of European‐ and African‐derived populations undergoing incipient sexual isolation. The existence of premating reproductive barriers has been previously established, but these types of barriers are not the only source shaping sexual isolation. To assess the influence of postmating barriers, we investigated putative postmating barriers of female remating and egg‐laying behavior, as well as hatchability of eggs laid and female longevity after mating. In the central region of our putative hybrid zone of American and Caribbean populations, we observed lower hatchability of eggs laid accompanied by increased resistance to harm after mating to less‐related males. These results illustrate that postmating reproductive barriers act alongside premating barriers and genetic admixture such as hybrid incompatibilities and influence early phases of sexual isolation.     

Characterizing a hybrid zone between a cryptic species pair of freshwater snails

Characterizing hybrid zones and their dynamics is a central goal in evolutionary biology, but this is particularly challenging for morphologically cryptic species. The lack of conspicuous divergence between parental types means intermediate hybrid forms often go undetected. We aimed to detect and characterize a suspected hybrid zone between a pair of morphologically cryptic lineages of the freshwater snail, Radix. We sampled Radix from across a contact zone between two mitochondrial lineages (Radix balthica and an undescribed lineage termed ‘MOTU3’) and detected admixture between two nuclear genotype clusters, which were significantly but not categorically associated with the mitochondrial lineages. Using a model selection approach, we show that the admixture cline is best explained by an interaction between precipitation and temperature gradients over the area, rather than geographic distance. We thus hypothesize that the correlation with climatic gradients suggests environmental selection has played a role in maintaining the hybrid zone. In a 2050 climate change scenario, we furthermore predict an expansion of one of the nuclear clusters and a widening of the hybrid zone as the climate warms and dries.     

The diurnal activity, movement and trypanosome infection rates of Glossina fuscipes fuscipes （Diptera： Glossinidae） in Buvuma Island, Lake Victoria, Uganda

The diurnal activity patterns, trypanosome infection rates and movement of Glossinafuscipesfuscipes （Diptera： Glossinidae） were investigated in Buvuma Island, Lake Victoria, Uganda. Hourly trapping of tsetse flies was undertaken to determine their activity rhythm while a capture-mark-release-recapture method was conducted to assess the movement and dispersal of tsetse flies between lakeshore, hinterland and further inland sites along a transected area. Dissection of tsetse flies was also undertaken to determine the trypanosome infection rates in salivary glands, proboscis and mid-gut. Results indicated a bimodal diurnal activity profile for G. f fuscipes on the Island, both on the lakeshore and in the hinterland. Movement and dispersal of G. f fuscipes tsetse flies occurred between lakeshore, hinterland and further inland sites with a greater tendency of flies to move to the lakeshore. Trypanosome infection rates of 4.32% for Trypasoma vivax and 1.15% for 7. congolense were found in G. f. fuscipes.     

Genetic panmixia within a narrow contact zone between chromosomally and ecologically distinct black fly sibling species (Diptera: Simuliidae)

Hybrid zones are windows into the speciation process, and their study can give clues into the maintenance and breakdown of species boundaries. Using both genetic and ecological tools, we investigate lineage diversification across a contact zone characterized by chromosome rearrangements. We show that black fly sibling species, Simulium arcticum sensu stricto (s.s.) and Simulium saxosum, lack genetic differentiation at both microsatellite and mtDNA loci in allopatry and sympatry, as well as exhibit high levels of gene flow and continuous chromosome variation in sympatry. Furthermore, hybrid frequencies at the contact zone are similar to those seen between races, rather than species. In contrast, S. arcticum s.s. and S. saxosum maintain ecological differences and distinct habitat associations ‐ the contact zone situated at the margin of suitable habitat for each sibling species. Moreover, gene flow occurs only in a narrow band along an ecological transition. Except for the contact zone, S. arcticum s.s. and S. saxosum hybrids do not occur elsewhere within the sibling species' ranges. Although S. arcticum s.s. and S. saxosum maintain the potential to interbreed freely, we conclude that habitat associations and, perhaps, chromosome systems prevent expansion of ranges and assimilation of lineages.     

Vegetation and the importance of insecticide-treated target siting for control of Glossina fuscipes fuscipes

Control of tsetse flies using insecticide-treated targets is often hampered by vegetation re-growth and encroachment which obscures a target and renders it less effective. Potentially this is of particular concern for the newly developed small targets (0.25 high × 0.5 m wide) which show promise for cost-efficient control of Palpalis group tsetse flies. Consequently the performance of a small target was investigated for Glossina fuscipes fuscipes in Kenya, when the target was obscured following the placement of vegetation to simulate various degrees of natural bush encroachment. Catches decreased significantly only when the target was obscured by more than 80%. Even if a small target is underneath a very low overhanging bush (0.5 m above ground), the numbers of G. f. fuscipes decreased by only about 30% compared to a target in the open. We show that the efficiency of the small targets, even in small (1 m diameter) clearings, is largely uncompromised by vegetation re-growth because G. f. fuscipes readily enter between and under vegetation. The essential characteristic is that there should be some openings between vegetation. This implies that for this important vector of HAT, and possibly other Palpalis group flies, a smaller initial clearance zone around targets can be made and longer interval between site maintenance visits is possible both of which will result in cost savings for large scale operations. We also investigated and discuss other site features e.g. large solid objects and position in relation to the water's edge in terms of the efficacy of the small targets.     

Prospects of a fungus-contamination device for the control of tsetse fly Glossina fuscipes fuscipes

The prospects of the fungus Metarhizium anisopliae (Metsch.) Sorok. applied in contamination devices (Cds) to control tsetse fly Glossina fuscipes fuscipes Newstead was tested in a field experiment in Lake Victoria from 2 March 1999 to 31 August 2000. One hundred and sixty pyramidal traps mounted with Cds were deployed along the lakeshore and rivers on Mfangano Island. Contamination devices were loaded with 1.5-2.0 g of dry conidia/Cd. On the second island, Nzenze Island, four pyramidal traps fitted with plastic bags were deployed and served as the conventional 'trap and kill' population suppression method. A third island, Ngodhe Island, remained untreated and served as a control. Cds were recharged monthly with fresh conidia; plastic bags were also changed monthly. The apparent changes in population density were monitored weekly using biconical traps set at random on the three islands. To assess the incidence of M. anisopliae in tsetse flies on Mfangano Island, flies captured during monitoring were maintained in the laboratory and their mortality recorded. Fly population was reduced to 82.4 and 95.8% relative to untreated control on Mfangano and Nzenze islands, respectively, during the experimental period. Compared to the fungus-treated island, the number of flies caught in monitoring traps increased considerably in 'trap kill' treatment at 5 months after the treatments were removed. The incidence of M. anisopliae in fly populations was low during the 12 weeks following the initiation of the experiment but increased afterward until termination of the treatment. M. anisopliae could still be recovered from fly populations at 3 months after termination of the treatment, although the incidence was low. The results of this study have shown that application of M. anisopliae in a contamination device can suppress the population of G. fuscipes fuscipes comparable to the 'trap and kill' technology.     

Patterns of genetic,phenotypic, and acoustic variation across a chiffchaff (Phylloscopus collybita abietinus/tristis) hybrid zone

Characterizing patterns of evolution of genetic and phenotypic divergence between incipient species is essential to understand how evolution of reproductive isolation proceeds. Hybrid zones are excellent for studying such processes, as they provide opportunities to assess trait variation in individuals with mixed genetic background and to quantify gene flow across different genomic regions. Here, we combine plumage, song, mtDNA and whole‐genome sequence data and analyze variation across a sympatric zone between the European and the Siberian chiffchaff (Phylloscopus collybita abietinus/tristis) to study how gene exchange between the lineages affects trait variation. Our results show that chiffchaff within the sympatric region show more extensive trait variation than allopatric birds, with a large proportion of individuals exhibiting intermediate phenotypic characters. The genomic differentiation between the subspecies is lower in sympatry than in allopatry and sympatric birds have a mix of genetic ancestry indicating extensive ongoing and past gene flow. Patterns of phenotypic and genetic variation also vary between regions within the hybrid zone, potentially reflecting differences in population densities, age of secondary contact, or differences in mate recognition or mate preference. The genomic data support the presence of two distinct genetic clades corresponding to allopatric abietinus and tristis and that genetic admixture is the force underlying trait variation in the sympatric region—the previously described subspecies (“fulvescens”) from the region is therefore not likely a distinct taxon. In addition, we conclude that subspecies identification based on appearance is uncertain as an individual with an apparently distinct phenotype can have a considerable proportion of the genome composed of mixed alleles, or even a major part of the genome introgressed from the other subspecies. Our results provide insights into the dynamics of admixture across subspecies boundaries and have implications for understanding speciation processes and for the identification of specific chiffchaff individuals based on phenotypic characters.     

Contrasting patterns of mitochondrial DNA and microsatellite introgressive hybridization between lineages of lake whitefish (Coregonus clupeaformis); relevance for speciation

We performed a combined analysis of mitochondrial DNA (mtDNA) and microsatellite loci among lake whitefish (Coregonus clupeaformis) populations in order to assess the levels of congruence between both types of markers in defining patterns of genetic structuring, introgressive hybridization and inferring population origins in the hybrid zone of the St. John River basin. A second objective was to test the hypothesis that secondary contact between glacial lineages always resulted in the occurrence of sympatric dwarf and normal whitefish ecotypes. Fish were sampled from 35 populations and polymorphism was screened at mtDNA and six microsatellite loci for a total of 688 and 763 whitefish, respectively. Four lakes harbouring a single whitefish population of normal ecotype admixed with mtDNA haplotypes of different lineages were found. This confirmed that secondary contact between whitefish evolutionary lineages did not always result in the persistence of reproductively isolated ecotypes. Microsatellites further supported the definition of distinct glacial lineages by identifying lineage-specific allelic size groups. They also further supported the hypothesis that ecotypes originated from either a single founding lineage (sympatric divergence) or following secondary contacts between lineages (allopatric divergence), depending on the lake. In general, however, the pattern of population differentiation and introgressive hybridization observed at microsatellites was in sharp contrast with that depicted by mtDNA variation. Both factorial correspondence analysis and analysis of admixture proportion revealed a much more pronounced pattern of introgressive hybridization than depicted by mtDNA analyses. Variable levels of introgression indicated that environmental differences may be as important as the historical contingency of secondary contact in explaining the persistence of sympatric ecotypes and the differential pattern of introgressive hybridization among lakes. Whitefish populations from the St. John River basin hybrid zone represent a rare illustration of a continuum of both morphological and genetic differentiation within a given taxon, spanning from complete introgression to possibly complete reproductive isolation, depending on lakes. Thus, each lake may be viewed as a different temporal snapshot taken throughout the gradual process of speciation.     

Simple identification of mitochondrial lineages in contact zones based on lineage-selective primers

A variety of research projects focus on genetic variation among and within maternal lineages as encompassed by mitochondrial DNA (mtDNA). While mtDNA often differs substantially between species, large differences may also be found within species. The evaluation of such divergent lineages, for example in intraspecific contact zones (hybrid zones), commonly involves sequencing numerous individuals. Large‐scale sequencing is both expensive and labour‐intensive. Based on sequences from 15 individuals, we devised a simple and quick polymerase chain reaction assay for identification of divergent mtDNA lineages in a secondary contact zone of the side‐blotched lizard (Uta stansburiana). The application uses lineage‐selective primers to amplify a lineage‐diagnostic product, and is based on each group of mtDNA haplotypes being a monophyletic assemblage of haplotypes sharing the same maternal ancestry, deeply divergent from the other group. The assay was tested on a larger sample (n = 147) of specimens from the contact zone, confirming its usefulness in quick and reliable identification of mtDNA lineages. This approach can be modified for other species, provided diagnostic lineage variation is available, and may also be performed in simple laboratory settings while conducting fieldwork.     

Polymorphic microsatellite markers for the tsetse fly Glossina fuscipes fuscipes (Diptera: Glossinidae), a vector of human African trypanosomiasis

Our understanding of Glossina fuscipes fuscipes, a major vector of sleeping sickness, has been severely constrained by a lack of genetic markers for mapping and population genetic studies. Here we present 10 newly developed microsatellite loci for this tsetse species. Heterozygosity levels in Moyo, an Ugandan population, averaged 0.57, with only two loci showing very low heterozygosity. Five loci carried more than six alleles. Together with five recently published microsatellite loci, this brings the number of available microsatellite loci for this species to 15. Their availability will greatly facilitate future studies on the genetics of this important human disease vector.     

Can alternative mating tactics facilitate introgression across a hybrid zone by circumventing female choice?

Reproductive barriers and divergence in species’ mate recognition systems underlie major models of speciation. However, hybridization between divergent species is common, and classic mechanisms to explain permeable reproductive barriers rarely consider how an individual may attain reproductive success. Alternative mating tactics (AMTs) exist in various forms across animal taxa. Such tactics may allow poorer quality individuals to gain mating opportunities and facilitate introgression either through asymmetrical positive selection or by circumventing female choice altogether in areas of secondary contact. One such tactic is satellite behaviour in frogs, where silent males perch near advertisers in an attempt to intercept females. To test whether such satellite male tactics are context‐dependent and favoured by hybrids, we genotyped and quantified the morphology of 80 male spring peeper (Pseudacris crucifer) individuals involved in caller–satellite associations from a secondary contact zone between two intraspecific mitochondrial lineages. Irrespective of population, satellite behaviour was best predicted by size but not body condition. Within the contact zone, pure individuals showed a significantly greater probability of being active callers, whereas hybrids of one lineage were more likely to adopt the satellite tactic. We suggest that satellite behaviour in P. crucifer promotes introgression, breaks down reproductive isolating barriers and contributes to asymmetrical introgression in this secondary contact zone. AMTs may thus be an underexplored but important alternative to oft‐discussed causes of genetic discordance found in hybrid zones.     

Evolution of a hybrid zone of two willow species (Salix L.) in the European Alps analyzed by RAD-seq and morphometrics

Natural hybridization of plants can result in many outcomes with several evolutionary consequences, such as hybrid speciation and introgression. Natural hybrid zones can arise in mountain systems as a result of fluctuating climate during the exchange of glacial and interglacial periods, where species retract and expand their territories, resulting in secondary contacts. Willows are a large genus of woody plants with an immense capability of interspecific crossing. In this study, the sympatric area of two diploid sister species, S. foetida and S. waldsteiniana in the eastern European Alps, was investigated to study the genomic structure of populations within and outside their contact zone and to analyze congruence of morphological phenotypes with genetic data. Eleven populations of the two species were sampled across the Alps and examined using phylogenetic network and population genetic structure analyses of RAD Seq data and morphometric analyses of leaves. The results showed that a homoploid hybrid zone between the two species was established within their sympatric area. Patterns of genetic admixture in homoploid hybrids indicated introgression with asymmetric backcrossing to not only one of the parental species but also one hybrid population forming a separate lineage. The lack of F1 hybrids indicated a long-term persistence of the hybrid populations. Insignificant isolation by distance suggests that gene flow can act over large geographical scales. Morphometric characteristics of hybrids supported the molecular data and clearly separated populations of the parental species, but showed intermediacy in the hybrid zone populations with a bias toward S. waldsteiniana. The homoploid hybrid zone might have been established via secondary contact hybridization, and its establishment was fostered by the low genetic divergence of parental species and a lack of strong intrinsic crossing barriers. Incomplete ecological separation and the ability of long-distance dispersal of willows could have contributed to the spatial expansion of the hybrid zone.     

Secondary contact seeds phenotypic novelty in cichlid fishes

Theory proposes that genomic admixture between formerly reproductively isolated populations can generate phenotypic novelty for selection to act upon. Secondary contact may therefore be a significant promoter of phenotypic novelty that allows species to overcome environmental challenges and adapt to novel environments, including during adaptive radiation. To date, this has largely been considered from the perspective of interspecific hybridization at contact zones. However, it is also possible that this process occurs more commonly between natural populations of a single species, and thus its importance in adaptive evolution may have been underestimated. In this study, we tested the consequences of genomic introgression during apparent secondary contact between phenotypically similar lineages of the riverine cichlid fish Astatotilapia calliptera. We provide population genetic evidence of a secondary contact zone in the wild, and then demonstrate using mate-choice experiments that both lineages can reproduce together successfully in laboratory conditions. Finally, we show that genomically admixed individuals display extreme phenotypes not observed in the parental lineages. Collectively, the evidence shows that secondary contact can drive the evolution of phenotypic novelty, suggesting that pulses of secondary contact may repeatedly seed genetic novelty, which when coupled with ecological opportunity could promote rapid adaptive evolution in natural circumstances.