Many a slip: dissecting the causes of reproductive isolation in two species of Tegenaria spiders (Agelenidae)

Key to our understanding of the mechanisms underlying the process of speciation is the determination of the nature of the barriers to gene flow between related taxa. Species that show zero gene flow in sympatry or parapatry are of little use in this respect. In the present study, we used two closely‐related species of large house spider, Tegenaria saeva and Tegenaria gigantea, which hybridize to a limited extent along a natural contact zone in southern Britain. The species are apparently indistinguishable with respect to habitat utilization and phenology. Laboratory crosses using individuals from both allopatric and parapatric populations suggest that, although male and female courtship, as well as web and cuticular‐borne pheromones, are conserved between the species, mechanical difficulties are experienced during interspecific copulation. Copulation bouts are, on average, significantly shorter during interspecific matings because of these difficulties, and are probably not sufficiently long for effective sperm transfer to take place. In the two cases of successful interspecific crossing, and in subsequent F₁ and backcross generations, there are few indications of differential fertility, fecundity or viability, suggesting little post‐zygotic incompatibility. The high success rate of crosses between F₁ hybrids and both parental species underlines the principally mechanical barrier to gene flow between these taxa. Once this is breached, there appears to be little impediment to continuing introgression, which could, in some geographical areas at least, ultimately lead to the fusion of the two species. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 355–367.     

The role of hybridization and introgression in maintaining species integrity and cohesion in naturally isolated inselberg bromeliad populations

• Hybridization is a widespread phenomenon present in numerous lineages across the tree of life. Its evolutionary consequences range from effects on the origin and maintenance, to the loss of biodiversity.
• We studied genetic diversity and intra‐ and interspecific gene flow between two sympatric populations of closely‐related species, Pitcairnia flammea and P. corcovadensis (Bromeliaceae), which are adapted to naturally fragmented Neotropical inselbergs, based on nuclear and plastidial DNA.
• Our main results indicate a strong reproductive isolation barrier, although low levels of interspecific gene flow were observed in both sympatric populations. The low rates of intraspecific gene flow observed for both P. corcovadensis and P. flammea populations corroborate the increasing body of evidence that inselberg bromeliad species are maintained as discrete evolutionary units despite the presence of low genetic connectivity. Nuclear patterns of genetic diversity and gene flow revealed that hybridization and introgression might not cause species extinction via genetic assimilation of the rare P. corcovadensis.
• In the face of reduced intraspecific gene exchange, hybridization and introgression may be important aspects of the Pitcairnia diversification process, with a positive evolutionary impact at the bromeliad community level, and thus contribute to increasing and maintaining genetic diversity in local isolated inselberg populations.

     

Isolation and characterization of microsatellite markers in the tropical tree species Dicorynia guianensis (Caesalpinaceae)

Dicorynia guianensis is a canopy tree endemic to French Guyana and the first to be exploited as timber wood. To investigate levels of diversity and dynamics of gene flow in a 40‐ha study plot, highly polymorphic microsatellite markers were developed. An enriched microsatellite library was constructed and seven microsatellite primer pairs were developed. Total genomic DNA was extracted and amplified from cambium tissue from 172 adult individuals. The level of genetic diversity was found to be low. A large heterozygote deficiency was found at one locus. This was resolved by redesigning the primer set.     

Hybridization between ecotypes in a phenotypically and ecologically heterogeneous population of Iris savannarum (Iridaceae) in Florida

Iris series Hexagonae is a small, monophyletic complex of five species and associated hybrid populations, popularly known as the ‘Louisiana irises’. Series Hexagonae has been recognized as a textbook case of introgressive hybridization based on numerous studies in Louisiana. We previously explored patterns of genetic structure and diversity in populations of the complex in Florida. Populations that occupy high, dry habitats have plants with a distinctive floral ‘highlands’ phenotype. Aquatic populations of I. savannarum also have a consistent floral phenotype (‘coastal’). Jacks Branch slough in Glades County harbors the largest population of series Hexagonae sampled in our previous studies. It is an ecologically heterogeneous site, punctuated by sandy uplands that remain above high water and harbor Iris with the highlands floral phenotype as well as coastal types in the wetter areas of the slough. We hypothesized that mesic, but non‐inundated, sites in the slough would host hybrids of the two phenotypic groups, and tested these hypotheses with 19 microsatellite loci. All data analyses support our hypothesis of hybridization between the upland and hydric subpopulations. Although two methods of introgression analysis identify some of the individuals from the admixture zone as first‐generation hybrids, subsequent admixture was asymmetric with introgression from the ‘highland’ parent. This is the first report of such assortative processes occurring in the confines of a single population of Iris. We suggest that the evolutionary processes described interspecifically for series Hexagonae commence early at the population level.     

Phylogeographic patterns and high levels of chloroplast DNA diversity in four Packera (asteraceae) species in southwestern Alberta

Chloroplast DNA (cpDNA) haplotype variation is compared among alpine and prairie/montane species of Packera from a region in southwestern Alberta that straddles the boundary of Pleistocene glaciation. The phylogeny of the 15 haplotypes identified reveals the presence of two groups: one generally found in coastal and northern species and the other from species in drier habitats. The presence of both groups in all four species and most populations from southwestern Alberta is evidence of past hybridization involving species or lineages that may no longer be present in the region. With the exception of the alpine P. subnuda (phiST = 1.0), interpopulational subdivision of haplotype variation is low (phiST < 0.350), suggesting that interpopulational gene flow is high. However, based on haplotype distribution patterns, we propose that Pleistocene hybridization and incomplete lineage sorting have resulted in reduced subdivision of interpopulational variation so that gene flow may not be as high as indicated. Drift has been more important in the alpine species populations, especially P. subnuda.     

Genetic isolation between coastal and fishery‐impacted,offshore bottlenose dolphin (Tursiops spp.) populations

The identification of species and population boundaries is important in both evolutionary and conservation biology. In recent years, new population genetic and computational methods for estimating population parameters and testing hypotheses in a quantitative manner have emerged. Using a Bayesian framework and a quantitative model‐testing approach, we evaluated the species status and genetic connectedness of bottlenose dolphin (Tursiops spp.) populations off remote northwestern Australia, with a focus on pelagic ‘offshore’ dolphins subject to incidental capture in a trawl fishery. We analysed 71 dolphin samples from three sites beyond the 50 m depth contour (the inshore boundary of the fishery) and up to 170 km offshore, including incidentally caught and free‐ranging individuals associating with trawl vessels, and 273 dolphins sampled at 12 coastal sites inshore of the 50 m depth contour and within 10 km of the coast. Results from 19 nuclear microsatellite markers showed significant population structure between dolphins from within the fishery and coastal sites, but also among dolphins from coastal sites, identifying three coastal populations. Moreover, we found no current or historic gene flow into the offshore population in the region of the fishery, indicating a complete lack of recruitment from coastal sites. Mitochondrial DNA corroborated our findings of genetic isolation between dolphins from the offshore population and coastal sites. Most offshore individuals formed a monophyletic clade with common bottlenose dolphins (T. truncatus), while all 273 individuals sampled coastally formed a well‐supported clade of Indo‐Pacific bottlenose dolphins (T. aduncus). By including a quantitative modelling approach, our study explicitly took evolutionary processes into account for informing the conservation and management of protected species. As such, it may serve as a template for other, similarly inaccessible study populations.     

Microsatellite analysis of genetic structure in the mangrove species Avicennia marina (Forsk.) Vierh. (Avicenniaceae)

The level of genetic variation throughout the entire worldwide range of the mangrove species Avicennia marina (Forsk.) Vierh. was examined using microsatellite markers. Three microsatellite loci detected high levels of allelic diversity (70 alleles in total), essential for an accurate estimation of population genetic parameters. The informativeness of the microsatellite loci tended to increase with increasing average number of repeats. The levels of heterozygosity detected for each population, over all loci, ranged from 0.0 to 0.8, with an average of 0.407, indicating that some populations had little or no genetic variation, whereas others had a large amount. Populations at the extremes of the distribution range showed reduced levels of heterozygosity, and significant levels of inbreeding. This is not unexpected as these populations may be subject to founder effects and environmental constraints. The presence of genetic structure was tested in A. marina populations using three models: (i) a single panmictic model; (ii) the discrete subpopulation model; and (iii) the isolation by distance model. The discrete subpopulations model was supported by the overall measures of population differentiation based on the infinite alleles model (F-statistics), and the stepwise mutation model (R statistics). In addition, an analysis of molecular variance (AMOVA), using both theoretical models, found that most of the variation was between populations (41-71%), and within individuals in the total population (31-49%). There was little variation among individuals within populations (0-10%). There was no significant isolation by distance. The high levels of genetic differentiation observed among populations of A. marina may be due to environmental and ecological factors, particularly past sea level and climatic changes.     

Population differentiation and species cohesion in two closely related plants adapted to neotropical high-altitude 'inselbergs', Alcantarea imperialis and Alcantarea geniculata (Bromeliaceae)

Isolated granitic rock outcrops or 'inselbergs' may provide a window into the molecular ecology and genetics of continental radiations under simplified conditions, in analogy to the use of oceanic islands in studies of species radiations. Patterns of variability and gene flow in inselberg species have never been thoroughly evaluated in comparison to related taxa with more continuous distribution ranges, or to other species in the same kingdom in general. We use nuclear microsatellites to study population differentiation and gene flow in two diploid, perennial plants adapted to high-altitude neotropical inselbergs, Alcantarea imperialis and Alcantarea geniculata (Bromeliaceae). Population differentiation is pronounced in both taxa, especially in A. imperialis. Gene flow in this species is considerably lower than expected from the literature on plants in general and Bromeliaceae in particular, and too low to prevent differentiation due to drift (N(e)m < 1), unless selection coefficients/effect sizes of favourable alleles are great enough to maintain species cohesion. Low gene flow in A. imperialis indicates that the ability of pollinating bats to promote gene exchange between inselbergs is smaller than previously assumed. Population subdivision in one inselberg population of A. imperialis appears to be associated with the presence of two colour morphs that differ in the coloration of rosettes and bracts. Our results indicate a high potential for inselbergs as venues for studies of the molecular ecology and genetics of continental radiations, such as the one that gave rise to the extraordinary diversity of adaptive strategies and phenotypes seen in Bromeliaceae.     

Development and characterization of chloroplast microsatellite markers for Cryptomeria japonica D. Don

We developed 17 chloroplast microsatellite markers, which consisted of seven mononucleotide microsatellites with a minimum repeat number of 10 and 10 dinucleotide microsatellites with a minimum repeat number of six, from the complete chloroplast genomic sequence of Cryptomeria japonica. A survey of 45 C. japonica individuals showed the number of alleles ranging from two to 11 alleles and a diversity index ranging from 0.085 to 0.895. Consequently, the 45 C. japonica individuals were divided into 39 haplotypes. These markers will be useful genetic markers in the gene flow analysis and population genetics of C. japonica.     

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.     

Development and characterization of 11 microsatellite markers in a widespread Neotropical seasonally dry forest tree species,Geoffroea spinosa Jacq. (Leguminosae)

Eleven dinucleotide microsatellites were developed in Geoffroea spinosa (Leguminosae), a widespread tree of the seasonally dry Neotropical forests, and characterized on six populations from Peru, Argentina and Paraguay. Four of them amplified on the Peruvian populations only, probably because of mutations in the microsatellite flanking regions in the other populations. Ten microsatellites were found polymorphic, with within population gene diversities ranging from 0.17 to 0.95, and a number of alleles varying from seven to 19. A significant overall genetic differentiation was also found (θ = 0.212; P < 0.01).     

Differential gene flow of mitochondrial and nuclear DNA markers among chromosomal races of Australian morabine grasshoppers (Vandiemenella, viatica species group)

Recent theoretical developments have led to a renewed interest in the potential role of chromosomal rearrangements in speciation. Australian morabine grasshoppers (genus Vandiemenella, viatica species group) provide an excellent study system to test this potential role of chromosomal rearrangements because they show extensive chromosomal variation and formed the basis of a classic chromosomal speciation model. There are three chromosomal races, viatica19, viatica17, and P24(XY), on Kangaroo Island, South Australia, forming five parapatric populations with four putative contact zones among them. We investigate the extent to which chromosomal variation among these populations may be associated with barriers to gene flow. Population genetic and phylogeographical analyses using 15 variable allozyme loci and the elongation factor-1alpha (EF-1alpha) gene indicate that the three races represent genetically distinct taxa. In contrast, analyses of the mitochondrial cytochrome c oxidase subunit I (COI) gene show the presence of three distinctive and geographically localized groups that do not correspond with the distribution of the chromosomal races. These discordant population genetic patterns are likely to result from introgressive hybridization between the chromosomal races and range expansions/contractions. Overall, these results suggest that reduction of nuclear gene flow may be associated with chromosomal variation, or underlying genetic variation linked with chromosomal variation, whereas mitochondrial gene flow appears to be independent of this variation in these morabine grasshoppers. The identification of an intact contact zone between P24(XY) and viatica17 offers considerable potential for further investigation of molecular mechanisms that maintain distinct nuclear genomes among the chromosomal races.     

Even population differentiation for maternal and biparental gene markers in Eugenia uniflora, a widely distributed species from the Brazilian coastal Atlantic rain forest

Brazilian cherry (Eugenia uniflora L.) occurs in the Brazilian Atlantic Coastal Tropical Forest from the north‐eastern state of Ceara, to the southernmost state of Rio Grande do Sul. E. uniflora plays an important role in the maintenance of ‘Restinga’ ecosystems, at the interface between forest and strand vegetation. Here, we characterize the genetic diversity, its allocation within and between populations and the possible components of nuclear and cytoplasmic gene flow that determine the spatial distribution of the genetic variability. Five E. uniflora populations encompassing distinct biogeographical components were sampled: populations at extremes of latitudinal distribution, at an oceanic island and close to the major urban centre of Rio de Janeiro. AFLP markers showed that genetic variance among studied populations was moderate (F_st = 0.211) with 78.9% of variability residing within populations. The estimated seed and pollen components of gene flow among populations were approximately equal, with a ratio between 1.003 and 0.713 in function of the cpDNA marker used. The five populations present a considerable genetic structure, as assigned by both nuclear and chloroplastic DNA markers. Our data suggest the existence of different glacial refugia and a limited pollen and seed gene flow, mainly between the southern and the other regions enclosed in the Atlantic rain forest. Therefore, any strategy of conservation and management has to assure the preservation of several populations along the Atlantic coast, to maintain the majority of the intra specific level of diversity.     

Hybridization and sequential components of reproductive isolation between parapatric walnut‐infesting sister species Rhagoletis completa and Rhagoletis zoqui

Hybridization can provide a window into how populations diverge to form new species. Here, we confirm hybridization between Rhagoletis completa Cresson, 1929 and Rhagoletis zoqui Bush, 1966, two species of walnut husk‐infesting flies that geographically overlap in a narrow area of parapatry in Northeastern Mexico. Rhagoletis completa and R. zoqui are members of a species group (Rhagoletis suavis) that has been hypothesized to speciate in allopatry. Sexual selection has been argued to be a potentially important factor for generating pre‐mating isolation among walnut husk flies, because of the differences in wing morphology and body coloration, and the existence of sexual dimorphism within species. However, there was no evidence for pre‐mating isolation between R. completa and R. zoqui, based on choice and no‐choice mating experiments conducted on adults of fly populations outside the contact zone. There was also no support for reduced fertility of hybrid matings or for F₁ inviability; however, F₁ hybrids appeared to have lower fertility and F₂ offspring have reduced survivorship. Postzygotic isolation in later generation hybrids of mixed ancestry therefore appears to be the first intrinsic barrier to gene flow evolving between R. completa and R. zoqui. We discuss the implications of our results for allopatric speciation in walnut flies and the potential evolutionary fate of R. zoqui and R. completa if they were to come into broad geographic contact in the future. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Chloroplast DNA variation and postglacial recolonization of common ash (Fraxinus excelsior L.) in Europe

We used chloroplast polymerase chain reaction-restriction-fragment length polymorphism (PCR-RFLP) and chloroplast microsatellites to assess the structure of genetic variation and postglacial history across the entire natural range of the common ash (Fraxinus excelsior L.), a broad-leaved wind-pollinated and wind-dispersed European forest tree. A low level of polymorphism was observed, with only 12 haplotypes at four polymorphic microsatellites in 201 populations, and two PCR-RFLP haplotypes in a subset of 62 populations. The clear geographical pattern displayed by the five most common haplotypes was in agreement with glacial refugia for ash being located in Iberia, Italy, the eastern Alps and the Balkan Peninsula, as had been suggested from fossil pollen data. A low chloroplast DNA mutation rate, a low effective population size in glacial refugia related to ash's life history traits, as well as features of postglacial expansion were put forward to explain the low level of polymorphism. Differentiation among populations was high (GST= 0.89), reflecting poor mixing among recolonizing lineages. Therefore, the responsible factor for the highly homogeneous genetic pattern previously identified at nuclear microsatellites throughout western and central Europe (Heuertz et al. 2004) must have been efficient postglacial pollen flow. Further comparison of variation patterns at both marker systems revealed that nuclear microsatellites identified complex differentiation patterns in south-eastern Europe which remained undetected with chloroplast microsatellites. The results suggest that data from different markers should be combined in order to capture the most important genetic patterns in a species.     

Fine-scale population structure in a desert amphibian: landscape genetics of the black toad (Bufo exsul)

Environmental variables can strongly influence a variety of intra- and inter-population processes, including demography, population structure and gene flow. When environmental conditions are particularly harsh for a certain species, investigating these effects is important to understanding how populations persist under difficult conditions. Furthermore, species inhabiting challenging environments present excellent opportunities to examine the effects of complex landscapes on population processes because these effects will often be more pronounced. In this study, I use 16 microsatellite loci to examine population structure, gene flow and demographic history in the black toad, Bufo exsul , which has one of the most restricted natural ranges of any amphibian. Bufo exsul inhabits four springs in the Deep Springs Valley high desert basin and has never been observed more than several meters from any source of water. My results reveal limited gene flow and moderately high levels of population structure ( F _ST = 0.051–0.063) between all but the two closest springs. I found that the geographic distance across the arid scrub habitat between springs is significantly correlated with genetic structure when distance accounts for topography and barriers to dispersal. I also found very low effective population sizes ( N _e = 7–30) and substantial evidence for historical population bottlenecks in all four populations. Together, these results suggest that the desert landscape and B.   exsul 's high habitat specificity contribute significantly to population structure and demography in this species and emphasize the importance of considering behavioural and landscape data in conservation genetic studies of natural systems.     

Properties and unbiased estimation of F- and D-statistics in samples containing related and inbred individuals

The Patterson F- and D-statistics are commonly used measures for quantifying population relationships and for testing hypotheses about demographic history. These statistics make use of allele frequency information across populations to infer different aspects of population history, such as population structure and introgression events. Inclusion of related or inbred individuals can bias such statistics, which may often lead to the filtering of such individuals. Here, we derive statistical properties of the F- and D-statistics, including their biases due to the inclusion of related or inbred individuals, their variances, and their corresponding mean squared errors. Moreover, for those statistics that are biased, we develop unbiased estimators and evaluate the variances of these new quantities. Comparisons of the new unbiased statistics to the originals demonstrates that our newly derived statistics often have lower error across a wide population parameter space. Furthermore, we apply these unbiased estimators using several global human populations with the inclusion of related individuals to highlight their application on an empirical dataset. Finally, we implement these unbiased estimators in open-source software package funbiased for easy application by the scientific community.     

Genetic analysis of a chromosomal hybrid zone in the Australian morabine grasshoppers (Vandiemenella, viatica species group)

Whether chromosomal rearrangements promote speciation by providing barriers to gene exchange between populations is one of the long-standing debates in evolutionary biology. This question can be addressed by studying patterns of gene flow and selection in hybrid zones between chromosomally diverse taxa. Here we present results of the first study of the genetic structure of a hybrid zone between chromosomal races of morabine grasshoppers Vandiemenella viatica , P24(XY) and viatica 17, on Kangaroo Island, Australia. Chromosomal and 11 nuclear markers revealed a narrow hybrid zone with strong linkage disequilibrium and heterozygote deficits, most likely maintained by a balance between dispersal and selection. Widths and positions of clines for these markers are concordant and coincident, suggesting that selection is unlikely to be concentrated on a few chromosomes. In contrast, a mitochondrial marker showed a significantly wider cline with centre offset toward the P24(XY) side. We argue that the discordance between the mitochondrial and nuclear/chromosomal clines and overall asymmetry of the clines suggest a secondary origin of the contact zone and potential movement of the zone after contact. Genome-wide scans using many genetic markers and chromosomal mapping of these markers are needed to investigate whether chromosomal differences directly reduce gene flow after secondary contact.