Reproductive isolation and patterns of genetic differentiation in a cryptic butterfly species complex

Molecular studies of natural populations are often designed to detect and categorize hidden layers of cryptic diversity, and an emerging pattern suggests that cryptic species are more common and more widely distributed than previously thought. However, these studies are often decoupled from ecological and behavioural studies of species divergence. Thus, the mechanisms by which the cryptic diversity is distributed and maintained across large spatial scales are often unknown. In 1988, it was discovered that the common Eurasian Wood White butterfly consisted of two species (Leptidea sinapis and Leptidea reali), and the pair became an emerging model for the study of speciation and chromosomal evolution. In 2011, the existence of a third cryptic species (Leptidea juvernica) was proposed. This unexpected discovery raises questions about the mechanisms preventing gene flow and about the potential existence of additional species hidden in the complex. Here, we compare patterns of genetic divergence across western Eurasia in an extensive data set of mitochondrial and nuclear DNA sequences with behavioural data on inter‐ and intraspecific reproductive isolation in courtship experiments. We show that three species exist in accordance with both the phylogenetic and biological species concepts and that additional hidden diversity is unlikely to occur in Europe. The Leptidea species are now the best studied cryptic complex of butterflies in Europe and a promising model system for understanding the formation of cryptic species and the roles of local processes, colonization patterns and heterospecific interactions for ecological and evolutionary divergence.     

Population structure in two sympatric species of the Lake Tanganyika cichlid tribe Eretmodini: evidence for introgression

Patterns of genetic differentiation were analysed and compared in two sympatric species of the endemic Lake Tanganyika cichlid tribe Eretmodini by means of mitochondrial DNA (mtDNA) sequences of the control region and six microsatellite DNA loci. The sample area covers a total of 138 km of mostly uninterrupted rocky shoreline in the Democratic Republic of Congo and includes the entire distribution range of Tanganicodus cf. irsacae that stretches over a distance of 35 km. Both markers detected significant genetic differentiation within and between the two species. T. cf. irsacae contained lower overall genetic variation than Eretmoduscyanostictus, possibly due to its more restricted range of distribution and its smaller effective population sizes. Complete fixation of Tanganicodus mtDNA haplotypes was observed in Eretmodus at two localities, while at two other localities some Tanganicodus individuals possessed Eretmodus mtDNA haplotypes. Taking into account the relatively large average sequence divergence of 6.2% between the two species, as well as the geographical distribution of mtDNA haplotypes in the lake, the observed pattern is more likely to be a consequence of asymmetric introgression than of shared ancestral polymorphism. As there is significant population differentiation between sympatric Tanganicodus and Eretmodus populations, the events of introgressions may have happened after secondary contact, but our data provide no evidence for ongoing gene flow and suggest that both species are reproductively isolated at present time.     

Echoes of a distant time: effects of historical processes on contemporary genetic patterns in Galaxias platei in Patagonia

Interpreting the genetic structure of a metapopulation as the outcome of gene flow over a variety of timescales is essential for the proper understanding of how changes in landscape affect biological connectivity. Here we contrast historical and contemporary connectivity in two metapopulations of the freshwater fish Galaxias platei in northern and southernmost Patagonia where paleolakes existed during the Holocene and Pleistocene, respectively. Contemporary gene flow was mostly high and asymmetrical in the northern system while extremely reduced in the southernmost system. Historical migration patterns were high and symmetric in the northern system and high and largely asymmetric in the southern system. Both systems showed a moderate structure with a clear pattern of isolation by distance (IBD). Effective population sizes were smaller in populations with low contemporary gene flow. An approximate Bayesian computation (ABC) approach suggests a late Holocene colonization of the lakes in the northern system and recent divergence of the populations from refugial populations from east and west of the Andes. For the southern system, the ABC approach reveals that some of the extant G. platei populations most likely derive from an ancestral population inhabiting a large Pleistocene paleolake while the rest derive from a higher‐altitude lake. Our results suggest that neither historical nor contemporary processes individually fully explain the observed structure and geneflow patterns and both are necessary for a proper understanding of the factors that affect diversity and its distribution. Our study highlights the importance of a temporal perspective on connectivity to analyse the diversity of spatially complex metapopulations.     

Genetic differentiation and diversity of Acacia koa populations in the Hawaiian Islands

Acacia koa A. Gray (koa) is a leguminous tree endemic to the Hawaiian Islands and can be divided into morphologically distinguishable groups of A. koaia Hillebrand, A. koa and populations that are intermediate between these extremes. The objectives of this investigation were to distinguish among divergent groups of koa at molecular levels, and to determine genetic diversity within and among the groups. Phylogenetic analyses using the ITS/5.8S rDNA and trnK intron sequences did not separate the representative koa types into distinct clusters. An unweighted pair group method with arithmetic mean cluster analysis and principal coordinate analysis, based on allele profiles of 12 microsatellite loci for 215 individual koa samples, separated the population into three distinct clusters consistent with their morphology, A. koaia, A. koa and intermediate forms. There was an average of 8.8 alleles per polymorphic locus (AP) among all koa and koaia individuals. The intermediate populations had the highest genetic diversity (H′ = 1.599), AP (7.9) and total number of unique alleles (21), whereas A. koaia and A. koa showed similar levels of genetic diversity (H′ = 0.965 and 0.943, respectively). No correlation was observed between geographic distance and genetic distance as determined by a Mantel test (r = 0.027, P = 0.91). The data presented here support previous recommendations that morphological variation within koa should be recognized at the subspecific level rather than as distinct species.     

Climatic changes can drive the loss of genetic diversity in a Neotropical savanna tree species

The high rates of future climatic changes, compared with the rates reported for past changes, may hamper species adaptation to new climates or the tracking of suitable conditions, resulting in significant loss of genetic diversity. Trees are dominant species in many biomes and because they are long‐lived, they may not be able to cope with ongoing climatic changes. Here, we coupled ecological niche modelling (ENM) and genetic simulations to forecast the effects of climatic changes on the genetic diversity and the structure of genetic clusters. Genetic simulations were conditioned to climatic variables and restricted to plant dispersal and establishment. We used a Neotropical savanna tree as species model that shows a preference for hot and drier climates, but with low temperature seasonality. The ENM predicts a decreasing range size along the more severe future climatic scenario. Additionally, genetic diversity and allelic richness also decrease with range retraction and climatic genetic clusters are lost for both future scenarios, which will lead genetic variability to homogenize throughout the landscape. Besides, climatic genetic clusters will spatially reconfigure on the landscape following displacements of climatic conditions. Our findings indicate that climate change effects will challenge population adaptation to new environmental conditions because of the displacement of genetic ancestry clusters from their optimal conditions.     

Long-term survival despite low genetic diversity in the critically endangered Madagascar fish-eagle

The critically endangered Madagascar fish-eagle ( Haliaeetus vociferoides ) is considered to be one of the rarest birds of prey globally and at significant risk of extinction. In the most recent census, only 222 adult individuals were recorded with an estimated total breeding population of no more than 100–120 pairs. Here, levels of Madagascar fish-eagle population genetic diversity based on 47 microsatellite loci were compared with its sister species, the African fish-eagle ( Haliaeetus vocifer ), and 16 of these loci were also characterized in the white-tailed eagle ( Haliaeetus albicilla ) and the bald eagle ( Haliaeetus leucocephalus ). Overall, extremely low genetic diversity was observed in the Madagascar fish-eagle compared to other surveyed Haliaeetus species. Determining whether this low diversity is the result of a recent bottleneck or a more historic event has important implications for their conservation. Using a Bayesian coalescent-based method, we show that Madagascar fish-eagles have maintained a small effective population size for hundreds to thousands of years and that its low level of neutral genetic diversity is not the result of a recent bottleneck. Therefore, efforts made to prevent Madagascar fish-eagle extinction should place high priority on maintenance of habitat requirements and reducing direct and indirect human persecution. Given the current rate of deforestation in Madagascar, we further recommend that the population be expanded to occupy a larger geographical distribution. This will help the population persist when exposed to stochastic factors (e.g. climate and disease) that may threaten a species consisting of only 200 adult individuals while inhabiting a rapidly changing landscape.     

Parallel evolution of dwarf ecotypes in the forest tree Eucalyptus globulus

Three small populations of a dwarf ecotype of the forest tree Eucalyptus globulus are found on exposed granite headlands in south-eastern Australia. These populations are separated by at least 100 km. Here, we used 12 nuclear microsatellites and a chloroplast DNA marker to investigate the genetic affinities of the dwarf populations to one another and to their nearest populations of tall E. globulus. Cape Tourville was studied in greater detail to assess the processes enabling the maintenance of distinct ecotypes in close geographical proximity. The three dwarf populations were not related to one another and were more closely related to adjacent tall trees than to one another. At Cape Tourville the dwarf and tall ecotypes were significantly differentiated in microsatellites and in chloroplast DNA. The dwarf and tall populations differed in flowering time and no evidence of pollen dispersal from the more extensive tall to the dwarf population was found. The three dwarf populations have evolved in parallel from the local tall ecotypes. This study shows that small marginal populations of eucalypts are capable of developing reproductive isolation from nearby larger populations through differences in flowering time and/or minor spatial separation, making parapatric speciation possible.     

Genetic variation in a tropical tree species influences the associated epiphytic plant and invertebrate communities in a complex forest ecosystem

Genetic differences among tree species, their hybrids and within tree species are known to influence associated ecological communities and ecosystem processes in areas of limited species diversity. The extent to which this same phenomenon occurs based on genetic variation within a single tree species, in a diverse complex ecosystem such as a tropical forest, is unknown. The level of biodiversity and complexity of the ecosystem may reduce the impact of a single tree species on associated communities. We assessed the influence of within-species genetic variation in the tree Brosimum alicastrum (Moraceae) on associated epiphytic and invertebrate communities in a neotropical rainforest. We found a significant positive association between genetic distance of trees and community difference of the epiphytic plants growing on the tree, the invertebrates living among the leaf litter around the base of the tree, and the invertebrates found on the tree trunk. This means that the more genetically similar trees are host to more similar epiphyte and invertebrate communities. Our work has implications for whole ecosystem conservation management, since maintaining sufficient genetic diversity at the primary producer level will enhance species diversity of other plants and animals.     

Discordant patterns of genetic and phenotypic differentiation in five grasshopper species codistributed across a microreserve network

Conservation plans can be greatly improved when information on the evolutionary and demographic consequences of habitat fragmentation is available for several codistributed species. Here, we study spatial patterns of phenotypic and genetic variation among five grasshopper species that are codistributed across a network of microreserves but show remarkable differences in dispersal‐related morphology (body size and wing length), degree of habitat specialization and extent of fragmentation of their respective habitats in the study region. In particular, we tested the hypothesis that species with preferences for highly fragmented microhabitats show stronger genetic and phenotypic structure than codistributed generalist taxa inhabiting a continuous matrix of suitable habitat. We also hypothesized a higher resemblance of spatial patterns of genetic and phenotypic variability among species that have experienced a higher degree of habitat fragmentation due to their more similar responses to the parallel large‐scale destruction of their natural habitats. In partial agreement with our first hypothesis, we found that genetic structure, but not phenotypic differentiation, was higher in species linked to highly fragmented habitats. We did not find support for congruent patterns of phenotypic and genetic variability among any studied species, indicating that they show idiosyncratic evolutionary trajectories and distinctive demographic responses to habitat fragmentation across a common landscape. This suggests that conservation practices in networks of protected areas require detailed ecological and evolutionary information on target species to focus management efforts on those taxa that are more sensitive to the effects of habitat fragmentation.     

Weak population differentiation in northern European populations of the endangered anadromous clupeid Alosa fallax

Examination of genetic structure among different spawning populations of the endangered twaite shad Alosa fallax showed statistically significant but weak differentiation estimated using nine microsatellite loci. Results point to substantial gene flow among the majority of individual estuary/river spawning locations, and no evidence that habitat disruption has led to increased local genetic drift.     

Natural selection drives clinal life history patterns in the perennial sunflower species, Helianthus maximiliani

In plants, ecologically important life history traits often display clinal patterns of population divergence. Such patterns can provide strong evidence for spatially varying selection across environmental gradients but also may result from nonselective processes, such as genetic drift, population bottlenecks and spatially restricted gene flow. Comparison of population differentiation in quantitative traits (measured as Q(ST) ) with neutral molecular markers (measured as F(ST) ) provides a useful tool for understanding the relative importance of adaptive and nonadaptive processes in the formation and maintenance of clinal variation. Here, we demonstrate the existence of geographic variation in key life history traits in the diploid perennial sunflower species Helianthus maximiliani across a broad latitudinal transect in North America. Strong population differentiation was found for days to flowering, growth rate and multiple size-related traits. Differentiation in these traits greatly exceeds neutral predictions, as determined both by partial Mantel tests and by comparisons of global Q(ST) values with theoretical F(ST) distributions. These findings indicate that clinal variation in these life history traits likely results from local adaptation driven by spatially heterogeneous environments.     

Population genetic analysis reveals a homoploid hybrid origin of Stephanomeria diegensis (Asteraceae)

Homoploid hybrid speciation has generally been viewed as a rare evolutionary phenomenon, with relatively few well-documented cases in nature. Here, we investigate the origin of Stephanomeria diegensis , a diploid flowering plant species that has been proposed to have arisen as a result of hybridization between S. exigua and S. virgata . Across the range of S. diegensis , all individuals share a common chloroplast haplotype with S. virgata while showing a greater affinity for S. exigua in terms of nuclear genetic diversity. A prinicipal coordinates analysis (PCO) based on the nuclear data revealed that S. diegensis is most similar to each parent along different axes. Moreover, a Bayesian clustering analysis as well as a hybrid index-based analysis showed evidence of mixed ancestry, with approximately two thirds of the S. diegensis nuclear genome derived from S. exigua . These results provide strong support for a homoploid hybrid origin of S. diegensis . Finally, contrary to the finding that homoploid hybrid species are typically multiply-derived, our results were most consistent with a single origin of this species.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

Microsatellite analysis reveals interpopulation differentiation and gene flow in the endangered tree Changiostyrax dolichocarpa (Styracaceae) with fragmented distribution in central China

Polymorphic simple sequence repeat (SSR) markers were used to investigate the impact of habitat fragmentation on the population structure and gene flow of Changiostyrax dolichocarpa, a critically endangered tree in central China. Intrapopulation genetic diversity, population structure and gene flow in the five extant populations of this species were analysed by eight SSR markers. Intrapopulation genetic diversity results suggest that C. dolichocarpa remnants maintained a relatively high degree of genetic diversity despite severe fragmentation. Low genetic differentiation among populations was found based on Wright's F(ST) and amova analysis. Both the F(ST)-based estimate and private allele method revealed high historical gene flow among the remnant populations. Recent immigrants, detected by assignment tests, tend to decrease from the grandparent generation to the current generation. The potentially highly restricted current gene flow among fragments may render the fragmented populations of C. dolichocarpa at a higher risk of local extinction several generations after fragmentation. Both in situ and ex situ conservation management for the remnant populations of C. dolichocarpa are therefore urgently needed to rescue remaining genetic diversity.     

Divergent and reticulate processes in evolution of Ethiopian Lophuromys flavopunctatus species complex: evidence from mitochondrial and nuclear DNA differentiation patterns

Molecular study of mitochondrial and nuclear genes and cytogenetic analysis were performed to examine possible patterns of speciation in the diverse Lophuromys flavopunctatus species complex of Ethiopia. Phylogenetic analysis of mtDNA data resulted in an unresolved bush of ten deeply diverged haplotype groups corresponding to potential species either well supported by various types of character or 'cryptic'. The cytogenetic analysis showed representatives of five of these mtDNA lineages to share an identical karyotype (2 n  = 70, NFa = 84), that has not been found previously in Ethiopia. One of them, L.  cf.  sikapusi , being a member of the L. flavopunctatus species complex, demonstrates remarkable morphological similarity to representatives of another species complex, L. sikapusi s.l ., which might be considered as a result of convergent evolution in analogous environments. Analysis of RAPD data suggests that at least two mtDNA types might have been subject to interspecific transfer due to hybridization. In the case of two sympatric haplotypes of L. brunneus we may assume that the contemporary pattern of variation between them can be explained by relatively recent hybridization with another distinct species, L. flavopunctatus . The formation of two groups belonging to distinct mitochondrial lineages within northern populations could be associated with more complex processes including ancient hybridization.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 83 , 301–316.     

Development of polymorphic microsatellite markers for Dalbergia nigra (Papilionoideae), an endangered tree from the Brazilian Atlantic Forest

Dalbergia nigra is an endangered tree restricted to the Brazilian Atlantic Forest. Ten polymorphic microsatellite loci were isolated and characterized in 47 trees from two populations. The total number of alleles per locus ranged from three to 12 alleles. The levels of observed and expected heterozygosities ranged from 0.304 to 0.740 and from 0.278 to 0.872, respectively. Significant deviations from Hardy-Weinberg were detected for only three loci in each population. No pair of loci exhibited significant linkage disequilibrium. These microsatellites provide an efficient tool to investigate genetic structure in forest remnants with the purpose of conservation of this species.     

Evolutionary and demographic processes shaping geographic patterns of genetic diversity in a keystone species,the African forest elephant (Loxodonta cyclotis)

The past processes that have shaped geographic patterns of genetic diversity may be difficult to infer from current patterns. However, in species with sex differences in dispersal, differing phylogeographic patterns between mitochondrial (mt) and nuclear (nu) DNA may provide contrasting insights into past events. Forest elephants (Loxodonta cyclotis) were impacted by climate and habitat change during the Pleistocene, which likely shaped phylogeographic patterns in mitochondrial (mt) DNA that have persisted due to limited female dispersal. By contrast, the nuclear (nu) DNA phylogeography of forest elephants in Central Africa has not been determined. We therefore examined the population structure of Central African forest elephants by genotyping 94 individuals from six localities at 21 microsatellite loci. Between forest elephants in western and eastern Congolian forests, there was only modest genetic differentiation, a pattern highly discordant with that of mtDNA. Nuclear genetic patterns are consistent with isolation by distance. Alternatively, male‐mediated gene flow may have reduced the previous regional differentiation in Central Africa suggested by mtDNA patterns, which likely reflect forest fragmentation during the Pleistocene. In species like elephants, male‐mediated gene flow erases the nuclear genetic signatures of past climate and habitat changes, but these continue to persist as patterns in mtDNA because females do not disperse. Conservation implications of these results are discussed.