Patterns of population genetic variation in sympatric chiltoniid amphipods within a calcrete aquifer reveal a dynamic subterranean environment

Calcrete aquifers from the Yilgarn region of arid central Western Australia contain an assemblage of obligate groundwater invertebrate species that are each endemic to single aquifers. Fine-scale phylogeographic and population genetic analyses of three sympatric and independently derived species of amphipod (Chiltoniidae) were carried out to determine whether there were common patterns of population genetic structure or evidence for past geographic isolation of populations within a single calcrete aquifer. Genetic diversity in amphipod mitochondrial DNA (cytochrome c oxidase subunit I gene) and allozymes were examined across a 3.5 km² region of the Sturt Meadows calcrete, which contains a grid of 115 bore holes (=wells). Stygobiont amphipods were found to have high levels of mitochondrial haplotype diversity coupled with low nucleotide diversity. Mitochondrial phylogeographic structuring was found between haplogroups for one of the chiltoniid species, which also showed population structuring for nuclear markers. Signatures of population expansion in two of the three species, match previous findings for diving beetles at the same site, indicating that the system is dynamic. We propose isolation of populations in refugia within the calcrete, followed by expansion events, as the most likely source of intraspecific genetic diversity, due to changes in water level influencing gene flow across the calcrete.     

Evolution of blind beetles in isolated aquifers: a test of alternative modes of speciation

Evidence is growing that not only allopatric but also sympatric speciation can be important in the evolution of species. Sympatric speciation has most convincingly been demonstrated in laboratory experiments with bacteria, but field-based evidence is limited to a few cases. The recently discovered plethora of subterranean diving beetle species in isolated aquifers in the arid interior of Australia offers a unique opportunity to evaluate alternative modes of speciation. This naturally replicated evolutionary experiment started 10-5 million years ago, when climate change forced the surface species to occupy geographically isolated subterranean aquifers. Using phylogenetic analysis, we determine the frequency of aquifers containing closely related sister species. By comparing observed frequencies with predictions from different statistical models, we show that it is very unlikely that the high number of sympatrically occurring sister species can be explained by a combination of allopatric evolution and repeated colonisations alone. Thus, diversification has occurred within the aquifers and likely involved sympatric, parapatric and/or microallopatric speciation.     

Low reproductive isolation and highly variable levels of gene flow reveal limited progress towards speciation between European river and brook lampreys

Ecologically based divergent selection is a factor that could drive reproductive isolation even in the presence of gene flow. Population pairs arrayed along a continuum of divergence provide a good opportunity to address this issue. Here, we used a combination of mating trials, experimental crosses and population genetic analyses to investigate the evolution of reproductive isolation between two closely related species of lampreys with distinct life histories. We used microsatellite markers to genotype over 1000 individuals of the migratory parasitic river lamprey (Lampetra fluviatilis) and freshwater‐resident nonparasitic brook lamprey (Lampetra planeri) distributed in 10 sympatric and parapatric population pairs in France. Mating trials, parentage analyses and artificial fertilizations demonstrated a low level of reproductive isolation between species even though size‐assortative mating may contribute to isolation. Most parapatric population pairs were strongly differentiated due to the joint effects of geographic distance and barriers to migration. In contrast, we found variable levels of gene flow between sympatric populations ranging from panmixia to moderate differentiation, which indicates a gradient of divergence with some population pairs that may correspond to alternative morphs or ecotypes of a single species and others that remain partially isolated. Ecologically based divergent selection may explain these variable levels of divergence among sympatric population pairs, but incomplete genome swamping following secondary contact could have also played a role. Overall, this study illustrates how highly differentiated phenotypes can be maintained despite high levels of gene flow that limit the progress towards speciation.     

Evolution of subterranean diving beetles (Coleoptera: Dytiscidae: Hydroporini, Bidessini) in the arid zone of Australia

Calcrete aquifers in arid inland Australia have recently been found to contain the world's most diverse assemblage of subterranean diving beetles (Coleoptera: Dytiscidae). In this study we test whether the adaptive shift hypothesis (ASH) or the climatic relict hypothesis (CRH) is the most likely mode of evolution for the Australian subterranean diving beetles by using a phylogeny based on two sequenced fragments of mitochondrial genes (CO1 and 16S-tRNA-ND1) and linearized using a relaxed molecular clock method. Most individual calcrete aquifers contain an assemblage of diving beetle species of distantly related lineages and/or a single pair of sister species that significantly differ in size and morphology. Evolutionary transitions from surface to subterranean life took place in a relatively small time frame between nine and four million years ago. Most of the variation in divergence times of the sympatric sister species is explained by the variation in latitude of the localities, which correlates with the onset of aridity from the north to the south and with an aridity maximum in the Early Pliocene (five mya). We conclude that individual calcrete aquifers were colonized by several distantly related diving beetle lineages. Several lines of evidence from molecular clock analyses support the CRH, indicating that all evolutionary transitions took place during the Late Miocene and Early Pliocene as a result of aridification.     

Evidence for speciation underground in diving beetles (Dytiscidae) from a subterranean archipelago

Most subterranean animals are assumed to have evolved from surface ancestors following colonization of a cave system; however, very few studies have raised the possibility of “subterranean speciation” in underground habitats (i.e., obligate cave‐dwelling organisms [troglobionts] descended from troglobiotic ancestors). Numerous endemic subterranean diving beetle species from spatially discrete calcrete aquifers in Western Australia (stygobionts) have evolved independently from surface ancestors; however, several cases of sympatric sister species raise the possibility of subterranean speciation. We tested this hypothesis using vision (phototransduction) genes that are evolving under neutral processes in subterranean species and purifying selection in surface species. Using sequence data from 32 subterranean and five surface species in the genus Paroster (Dytiscidae), we identified deleterious mutations in long wavelength opsin (lwop), arrestin 1 (arr1), and arrestin 2 (arr2) shared by a sympatric sister‐species triplet, arr1 shared by a sympatric sister‐species pair, and lwop and arr2 shared among closely related species in adjacent calcrete aquifers. In all cases, a common ancestor possessed the function‐altering mutations, implying they were already adapted to aphotic environments. Our study represents one of the first confirmed cases of subterranean speciation in cave insects. The assessment of genes undergoing pseudogenization provides a novel way of testing modes of speciation and the history of diversification in blind cave animals.     

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.     

Contrasting diversity and demographic signals in sympatric narrow‐range endemic shrubs of the south‐west Western Australian semi‐arid zone

Comparative studies of sympatric species that integrate both phylogeographical and population genetic approaches provide insight into how demographic events and life history traits shape adaptive potential and drive species persistence. Such studies are rare for species‐rich and strongly structured environments, especially those of the southern hemisphere. For two sympatric, perennial shrubs of the south‐west Western Australian semi‐arid zone, Grevillea globosa and Mirbelia sp. Bursarioides, we assessed historical and contemporary genetic diversity and structure, demographic processes and ratios of pollen to seed dispersal. Phylogeographical structure was not detected and haplotype networks were star‐like. Number of haplotypes, nucleotide diversity, haplotype diversity, and allelic diversity were statistically significantly lower for G. globosa than for M. sp. Bursarioides. Levels of haplotype divergence and more contemporary genetic divergence and expected heterozygosity were lower for G. globosa than for M. sp. Bursarioides, but differences were not statistically significant. Both species exhibited signals of isolation by distance and low pollen to seed dispersal ratios (5.26:1 and 6.88:1). Grevillea globosa displayed signals of historical and contemporary demographic expansion. Results imply an important role for aspects of seed ecology that impact population demography, as well as direct dispersal and a significant contribution of seed dispersal to genetic connectivity in a semi‐arid landscape.     

Contrasting genetic diversity and population structure among three sympatric Madagascan shorebirds: parallels with rarity,endemism, and dispersal

Understanding the relative contributions of intrinsic and extrinsic factors to population structure and genetic diversity is a central goal of conservation and evolutionary genetics. One way to achieve this is through comparative population genetic analysis of sympatric sister taxa, which allows evaluation of intrinsic factors such as population demography and life history while controlling for phylogenetic relatedness and geography. We used ten conserved microsatellites to explore the population structure and genetic diversity of three sympatric and closely related plover species in southwestern Madagascar: Kittlitz's plover (Charadrius pecuarius), white‐fronted plover (C. marginatus), and Madagascar plover (C. thoracicus). Bayesian clustering revealed strong population structure in the rare and endemic Madagascar plover, intermediate population structure in the white‐fronted plover, and no detectable population structure in the geographically widespread Kittlitz's plover. In contrast, allelic richness and heterozygosity were highest for the Kittlitz's plover, intermediate for the white‐fronted plover and lowest for the Madagascar plover. No evidence was found in support of the “watershed mechanism” proposed to facilitate vicariant divergence of Madagascan lemurs and reptiles, which we attribute to the vagility of birds. However, we found a significant pattern of genetic isolation by distance among populations of the Madagascar plover, but not for the other two species. These findings suggest that interspecific variation in rarity, endemism, and dispersal propensity may influence genetic structure and diversity, even in highly vagile species.     

Molecular evidence of a peripatric origin for two sympatric species of field crickets (Gryllus rubens and G. texensis) revealed from coalescent simulations and population genetic tests

Species pairs that differ primarily in characters involved in mating interactions and are largely sympatric raise intriguing questions about the mode of speciation. When species divergence is relatively recent, the footprint of the demographic history during speciation might be preserved and used to reconstruct the biogeography of species divergence. In this study, patterns of genetic variation were examined throughout the geographical range of two cryptic sister taxa of field crickets, Gryllus texensis and G. rubens; mitochondrial cytochrome oxidase I (COI) was sequenced in 365 individuals sampled from 48 localities. Despite significant molecular divergence between the species, they were not reciprocally monophyletic. We devised several analyses to statistically explore what historical processes might have given rise to this genealogical structure. The analyses indicated that the biogeographical pattern of genetic variation does not support a model of recent gene flow between species. Instead, coalescent simulations suggested that the genealogical structure within G. texensis, namely a deep split between two geographically overlapping clades, reflects historical substructure within G. texensis. Additional tests that consider the concentration of G. rubens haplotypes in one of the two G. texensis genetic clusters suggest a model of speciation in which G. rubens was derived from one lineage of a geographically subdivided ancestor. These results indicate that, despite the contemporary sympatry of G. texensis and G. rubens, the data are indicative of an peripatric origin in which G. rubens was derived from one of the two historical partitions in the species currently recognized as G. texensis. This proposed model of species divergence suggests how the interplay of geography and selection may give rise to new species, although this requires testing with multilocus data. Specifically, the model highlights how that geographical partitioning of ancestral variation in the past may augment the selectively driven divergence of characters involved in the reproductive isolation of the species today.     

Isolation by distance and gene flow in the Eurasian badger (Meles meles) at both a local and broad scale

Eurasian badgers, Meles meles, have been shown to possess limited genetic population structure within Europe; however, field studies have detected high levels of philopatry, which are expected to increase population structure. Population structure will be a consequence of both contemporary dispersal and historical processes, each of which is expected to be evident at a different scale. Therefore, to gain a greater understanding of gene flow in the badger, we examined microsatellite diversity both among and within badger populations, focusing on populations from the British Isles and western Europe. We found that while populations differed in their allelic diversity, the British Isles displayed a similar degree of diversity to the rest of western Europe. The lower genetic diversity occurring in Ireland, Norway and Scotland was more likely to have resulted from founder effects rather than contemporary population density. While there was significant population structure (F ST = 0.19), divergence among populations was generally well explained by geographic distance (P < 0.0001) across the entire range studied of more than 3000 km. Transient effects from the Pleistocene appear to have been replaced by a strong pattern of genetic isolation by distance across western Europe, suggestive of colonization from a single refugium. Analysis of individuals within British populations through Mantel tests and spatial autocorrelation demonstrated that there was significant local population structure across 3-30 km, confirming that dispersal is indeed restricted. The isolation by distance observed among badger populations across western Europe is likely to be a consequence of this restricted local dispersal.     

Apparent 'sympatric' speciation in ecologically similar herbivorous beetles facilitated by multiple colonizations of an island

Coexistence of recently diverged and ecologically similar sister species in complete sympatry represents a particularly compelling case for sympatric speciation. This study investigates the possible sympatric origin of two coexisting bark beetle taxa that utilize the same host plant on the island of La Palma in the Canary Islands. Aphanarthrum subglabrum and Aphanarthrum glabrum ssp. nudum breed inside dead twigs of Euphorbia lamarckii plants and are closely related to the allopatric A. glabrum ssp. glabrum in Tenerife, El Hierro and La Gomera. We tested the various speciation hypotheses in a genealogical context, using mitochondrial gene fragments from Cytochrome Oxidase I and 16S, and nuclear gene fragments from Enolase, Elongation Factor 1alpha and Histone H3. Phylogenetic analyses of the combined nuclear DNA data strongly supported a sister relationship between two sympatric and reproductively isolated taxa in La Palma. However, network analyses of subdivided nonrecombinant segments of the Enolase locus indicated a closer relationship between the two allopatric A. glabrum subspecies, suggesting multiple colonizations of this island. A bimodal distribution of mtDNA haplotypes in La Palma further documented the independent colonization of this island, with asymmetric introgression of mtDNA between two lineages. Consequently, the sympatric origin of the La Palma species is concluded to have involved allopatric phases before the parallel colonization of this island and subsequent introgression at some loci. The clear genetic and morphological evidence for reproductive isolation between these species suggests that the sympatric completion of divergence was either due to initial genetic incompatibility, morphological character displacement in male genitalia, or a combination of these factors.     

Comparative molecular phylogeography of two Xenopus species, X. gilli and X. laevis, in the south-western Cape Province, South Africa

Xenopus gilli is a vulnerable anuran with a patchy distribution along the south-western coast of the Cape Province, South Africa. This species is sympatric with Xenopus laevis laevis , a widespread relative found over much of southern Africa. We examined the molecular phylogeography and population structure of the contact zone between these species to obtain information about historical biogeography and conservation management of this region. Analyses of the distribution, frequency, and cladistic and phenetic relationships among mitochondrial DNA haplotypes indicate that population subdivision is present in both taxa but that long-term isolation of sets of populations has occurred in X. gilli only. Haplotype and nucleotide diversity are also considerably higher within and among X. gilli ponds than X. l. laevis ponds in this region. We attribute the genetic segregation of X. gilli populations to ancient habitat fragmentation by ocean transgression into X. gilli habitat and to continued habitat alteration by human activity. The lower level of genetic diversity in X. l. laevis in this region is likely a result of a recent arrival of this taxon to the south-western Cape region relative to X. gilli . Population structure in X. l. laevis may be a result of isolation by distance. Clear evidence exists for at least two management units within X. gilli and strongly supports the establishment of protective measures east of False Bay in order to conserve a substantial portion of this species' extant genetic diversity.     

Subterranean archipelago in the Australian arid zone: mitochondrial DNA phylogeography of amphipods from central Western Australia

In 1998, a unique subterranean ecosystem was discovered in numerous isolated calcrete (carbonate) aquifers in the arid Yilgarn region of Western Australia. Previous morphological and genetic analyses of a subterranean water beetle fauna suggest that calcrete aquifers are equivalent to closed island habitats that have been isolated for millions of years. We tested this hypothesis further by phylogeographic analyses of subterranean amphipods (Crangonyctoidea: Paramelitidae and Hyalidae) using mitochondrial DNA sequence data derived from the cytochrome oxidase I gene. Phylogenetic analyses and population genetic analyses (samova) provided strong evidence for the existence of at least 16 crangonyctoid and six hyalid divergent mitochondrial lineages, each restricted in their distribution to a single calcrete aquifer, in support of the 'subterranean island (archipelago) hypothesis' and extending its scope to include entirely water respiring invertebrates. Sequence divergence estimates between proximate calcrete populations suggest that calcretes have been isolated at least since the Pliocene, coinciding with a major aridity phase that led to the intermittent drying of surface water. The distribution of calcretes along palaeodrainage channels and on either side of drainage divides, have had less influence on the overall phylogeographic structure of populations, with evidence that ancestral crangonyctoid and hyalid species moved between catchments multiple times prior to their isolation within calcretes. At least two potential modes of evolution may account for the diversity of subterranean amphipod populations: dispersal/vicariance of stygobitic species or colonization of calcretes by surface species and independent evolution of stygobitic characteristics.     

Genital divergence in sympatric sister snails

A pattern of greater divergence in mating traits between sister‐species pairs with overlapping ranges than between allopatric species pairs is expected if reinforcement commonly contributes to speciation. Few large‐scale comparative analyses have addressed this prediction, especially for genital form. Here, we show that penial morphology follows the predicted pattern in 40 robustly identified sister‐species pairs in the marine gastropod subfamily Littorininae. Further work is needed to exclude other processes that may contribute to genital divergence between sympatric species, but the clear pattern we observe strongly suggests a role for genital form in reproductive isolation in this large clade.     

Population genomic evidence for multiple Pliocene refugia in a montane‐restricted harvestman (Arachnida,Opiliones, Sclerobunus robustus) from the southwestern United States

The integration of ecological niche modelling into phylogeographic analyses has allowed for the identification and testing of potential refugia under a hypothesis‐based framework, where the expected patterns of higher genetic diversity in refugial populations and evidence of range expansion of nonrefugial populations are corroborated with empirical data. In this study, we focus on a montane‐restricted cryophilic harvestman, Sclerobunus robustus, distributed throughout the heterogeneous Southern Rocky Mountains and Intermontane Plateau of southwestern North America. We identified hypothetical refugia using ecological niche models (ENMs) across three time periods, corroborated these refugia with population genetic methods using double‐digest RAD‐seq data and conducted population‐level phylogenetic and divergence dating analyses. ENMs identify two large temporally persistent regions in the mid‐latitude highlands. Genetic patterns support these two hypothesized refugia with higher genetic diversity within refugial populations and evidence for range expansion in populations found outside hypothesized refugia. Phylogenetic analyses identify five to six genetically divergent, geographically cohesive clades of S. robustus. Divergence dating analyses suggest that these separate refugia date to the Pliocene and that divergence between clades pre‐dates the late Pleistocene glacial cycles, while diversification within clades was likely driven by these cycles. Population genetic analyses reveal effects of both isolation by distance (IBD) and isolation by environment (IBE), with IBD more important in the continuous mountainous portion of the distribution, while IBE was stronger in the populations inhabiting the isolated sky islands of the south. Using model‐based coalescent approaches, we find support for postdivergence migration between clades from separate refugia.     

Population dynamics of the endangered plant, Phaedranassa tunguraguae, from the Tropical Andean hotspot

The Tropical Andes is a diversity hotspot for plants, but there is a scant knowledge about patterns of genetic variation within its constituent species. Phaedranassa tunguraguae is an IUCN endangered plant species endemic to a single valley in the Ecuadorian Andes. We estimate the levels of genetic differentiation across the geographic distribution of P. tunguraguae using 12 microsatellite loci. We discuss factors that might influence the genetic structure of this species. Genetic distance was used to evaluate relationship among populations and geographic patterns. Bayesian methods were used to investigate population structure, migration, evidence of recent bottlenecks, and time of divergence. The 7 populations form 2 genetic clusters. These clusters show highly significant differentiation between them, along with isolation by distance. Allele richness decreases from the most diverse westernmost population to the least diverse easternmost population. The species overall shows an excess of homozygotes, with highest levels of inbreeding in the easternmost population. We found evidence of recent bottleneck events. Migration rates were in general low but were higher between populations within each of the clusters. Time of divergence between populations was related to historical volcanic activity in the area. Based on our results, we propose 2 management units for P. tunguraguae.     

Comprehensive analysis of population genetics of Phoxinus phoxinus ujmonensis in the Irtysh River: Abiotic and biotic factors

As a widely distributed species along the Irtysh River, Phoxinus phoxinus ujmonensis (Kaschtschenko, 1899) was used as a model to investigate genetic diversity and population structure as well as the influence of environmental factors on population genetics. In this study, we specifically developed 12 polymorphic microsatellite loci. The analysis of microsatellite and mtDNA markers revealed a high and a moderate genetic diversity across seven populations, respectively. Moderate differentiation was also detected among several populations, indicating the impact of habitat fragmentation and divergence. The absence of isolation by distance implied an extensive gene flow, while the presence of isolation by adaptation implied that these populations might be in the process of adapting to divergent habitats. Correlation analysis showed that abiotic factors like dissolved oxygen, pH, total dissolved solids, and conductivity in water as well as biotic factors like plankton diversity and fish species diversity had impact on genetic diversity and divergence in P. phoxinus ujmonensis populations. The results of this study will provide an insight into the effect of environmental factors on genetic diversity and contribute to the study of population genetics of sympatric species.     

ANALYSIS OF MECHANISMS OF MICROEVOLUTIONARY CHANGE IN CEPPHUS GUILLEMOTS USING PATTERNS OF CONTROL REGION VARIATION

We surveyed population-level sequence variation in part of the mitochondrial control region for three species including eight subspecies of Cepphus guillemots (Charadriiformes: Alcidae) to test specific predictions about mechanisms of population differentiation. We found that sequences of spectacled guillemots (C. carbo) were more closely related to those of pigeon guillemots (C. columba; both found in the Pacific Ocean) than to those of black guillemots (C. grylle; Arctic and Atlantic Oceans), despite dissimilarities in plumage between spectacled guillemots and the other species. Distributions of species and timing of divergence events suggest that speciation involved allopatric and microallopatric populations isolated by Pleistocene glaciers. Control region sequences were significantly differentiated among populations within species and suggest that gene flow is low; however, populations are probably not in genetic equilibrium, so these results probably reflect historical isolation of colonies. In contrast, phylogenetic relationships among sequences within species were poorly resolved, probably because of a combination of incomplete lineage sorting and contemporary gene flow. Indices of genetic diversity provided no suggestion of recent bottlenecks in most populations, although two populations apparently underwent recent severe bottlenecks. Genetic divergence among populations was not correlated with geographic distance, which argues against isolation by distance. Results of these analyses, combined with breeding distributions and timing of divergence events, suggest that populations diverged during isolation in glacial refugia. Our results are consistent with earlier hypotheses posed by Storer and Udvardy.     

MULTIPLE MODES OF SPECIATION INVOLVED IN THE PARALLEL EVOLUTION OF SYMPATRIC MORPHOTYPES OF LAKE WHITEFISH (COREGONUS CLUPEAFORMIS,SALMONIDAE)

We performed a phylogenetic analysis of mtDNA variation among seven sympatric pairs of dwarf and normal morphotypes of whitefish from northern Québec and the St. John River drainage to address three questions relevant to understanding their radiation. Are all sympatric pairs reproductively isolated? Do phylogenetic analyses confirm that sympatric whitefish morphotypes found in eastern North America represent the outcome of polyphyletic evolutionary events? If so, did all sympatric pairs from the St. John River drainage originate from the same scenario of allopatric divergence and secondary contact? The hypothesis of genetic differentiation was supported for all sympatric pairs from the St. John River drainage, whereas lack of mtDNA diversity precluded any test of reproductive isolation for northern Québec populations. Patterns of mtDNA variation confirmed that dwarf and normal morphotypes evolved in parallel among independent, yet closely related, lineages, thus providing indirect evidence for the role of natural selection in promoting phenotypic radiation in whitefish. Patterns of mtDNA diversity among sympatric pairs of the St. John River indicated a complex picture of whitefish evolution that implied sympatric divergence and multiple allopatric divergence/secondary contact events on a small geographic scale. These results suggests that ecological opportunities, namely trophic niche availability, may promote population divergence in whitefish.     

Phenotypic divergence but not genetic distance predicts assortative mating among species of a cichlid fish radiation

The hypothesis of ecological divergence giving rise to premating isolation in the face of gene flow is controversial. However, this may be an important mechanism to explain the rapid multiplication of species during adaptive radiation following the colonization of a new environment when geographical barriers to gene flow are largely absent but underutilized niche space is abundant. Using cichlid fish, we tested the prediction of ecological speciation that the strength of premating isolation among species is predicted by phenotypic rather than genetic distance. We conducted mate choice experiments between three closely related, sympatric species of a recent radiation in Lake Mweru (Zambia/DRC) that differ in habitat use and phenotype, and a distantly related population from Lake Bangweulu that resembles one of the species in Lake Mweru. We found significant assortative mating among all closely related, sympatric species that differed phenotypically, but none between the distantly related allopatric populations of more similar phenotype. Phenotypic distance between species was a good predictor of the strength of premating isolation, suggesting that assortative mating can evolve rapidly in association with ecological divergence during adaptive radiation. Our data also reveals that distantly related allopatric populations that have not diverged phenotypically, may hybridize when coming into secondary contact, e.g. upon river capture because of diversion of drainage systems.