Hidden Biodiversity in an Ecologically Important Freshwater Amphipod: Differences in Genetic Structure between Two Cryptic Species

Cryptic species, i.e. species that are morphologically hard to distinguish, have been detected repeatedly in various taxa and ecosystems. In order to evaluate the importance of this finding, we have to know in how far cryptic species differ in various aspects of their biology. The amphipod Gammarus fossarum is a key invertebrate in freshwater streams and contains several cryptic species. We examined the population genetic structure, genetic diversity and demographic history of two of them (type A and type B) using microsatellite markers and asked whether they show significant differences. We present results of population genetic analyses based on a total of 37 populations from the headwaters of two major European drainages, Rhine and Rhone. We found that, in both species, genetic diversity was geographically structured among and within drainages. For type A in the Rhine and type B in the Rhone, we detected significant patterns of isolation by distance. The increase of genetic differentiation with geographical distance, however, was much higher in type A than in type B. This result indicates substantial interspecific differences in population history and/or the extent of current gene flow between populations. In the Rhine, type B does not show evidence of isolation by distance, and population differentiation is relatively low across hundreds of kilometres. The majority of these populations also show signatures of recent bottlenecks. These patterns are consistent with a recent expansion of type B into the Rhine drainage. In summary, our results suggest considerable and previously unrecognized interspecific differences in the genetic structure of these cryptic keystone species.

Complexities in the genetic structure of Anopheles gambiae populations in west Africa as revealed by microsatellite DNA analysis

Chromosomal forms of Anopheles gambiae, given the informal designations Bamako, Mopti, and Savannah, have been recognized by the presence or absence of four paracentric inversions on chromosome 2. Studies of karyotype frequencies at sites where the forms occur in sympatry have led to the suggestion that these forms represent species. We conducted a study of the genetic structure of populations of An. gambiae from two villages in Mali, west Africa. Populations at each site were composed of the Bamako and Mopti forms and the sibling species, Anopheles arabiensis. Karyotypes were determined for each individual mosquito and genotypes at 21 microsatellite loci determined. A number of the microsatellites have been physically mapped to polytene chromosomes, making it possible to select loci based on their position relative to the inversions used to define forms. We found that the chromosomal forms differ at all loci on chromosome 2, but there were few differences for loci on other chromosomes. Geographic variation was small. Gene flow appears to vary among different regions within the genome, being lowest on chromosome 2, probably due to hitchhiking with the inversions. We conclude that the majority of observed genetic divergence between chromosomal forms can be explained by forces that need not involve reproductive isolation, although reproductive isolation is not ruled out. We found low levels of gene flow between the sibling species Anopheles gambiae and Anopheles arabiensis, similar to estimates based on observed frequencies of hybrid karyotypes in natural populations.

Effects of Interspecific Gene Flow on the Phenotypic Variance-Covariance Matrix in Lake Victoria Cichlids

Quantitative genetics theory predicts adaptive evolution to be constrained along evolutionary lines of least resistance. In theory, hybridization and subsequent interspecific gene flow may however rapidly change the evolutionary constraints of a population and eventually change its evolutionary potential, but empirical evidence is still scarce. Using closely related species pairs of Lake Victoria cichlids sampled from four different islands with different levels of interspecific gene flow, we tested for potential effects of introgressive hybridization on phenotypic evolution in wild populations. We found that these effects differed among our study species. Constraints measured as the eccentricity of phenotypic variance-covariance matrices declined significantly with increasing gene flow in the less abundant species for matrices that have a diverged line of least resistance. In contrast we find no such decline for the more abundant species. Overall our results suggest that hybridization can change the underlying phenotypic variance-covariance matrix, potentially increasing the adaptive potential of such populations.

Mitochondrial control region diversity of the houbara bustard Chlamydotis undulata complex and genetic structure along the Atlantic seaboard of North Africa

The houbara bustard, Chlamydotis undulata, is a declining cryptic desert bird whose range extends from North Africa to Central Asia. Three subspecies are currently recognized by geographical distribution and morphology: C.u.fuertaventurae, C.u.undulata and C.u.macqueenii. We have sequenced 854 bp of mitochondrial control region from 73 birds to describe their population genetic structure with a particular sampling focus on the connectivity between C.u.fuertaventurae and C.u.undulata along the Atlantic seaboard of North Africa. Nucleotide and haplotypic diversity varied among the subspecies being highest in C.u.undulata, lowest in C.u.fuertaventurae and intermediate in C.u.macqueenii. C.u.fuertaventurae and C.u.undulata are paraphyletic and an average nucleotide divergence of 2.08% splits the later from C.u.macqueenii. We estimate that C.u.fuertaventurae and C.u.undulata split from C.u.macqueenii approximately 430 000 years ago. C.u.fuertaventurae and C.u.undulata are weakly differentiated (F-ST = 0.27, N-m = 1.3), indicative of a recent shared history. Archaeological evidence indicates that houbara bustards have been present on the Canary Islands for 130-170 000 years. However, our genetic data point to a more recent separation of C.u.fuertaventurae and C.u.undulata at around 20-25 000 years. Concordant archaeological, climatic opportunities for colonization and genetic data point to a scenario of: (i) initial colonization of the Canary Islands about 130 000 years ago; (ii) a period of secondary contact 19-30 000 years ago homogenizing any pre-existing genetic structure followed by; (iii) a period of relative isolation that persists today.

Optimal sampling strategy for estimation of spatial genetic structure in tree populations

Fine-scale spatial genetic structure (SGS) in natural tree populations is largely a result of restricted pollen and seed dispersal. Understanding the link between limitations to dispersal in gene vectors and SGS is of key interest to biologists and the availability of highly variable molecular markers has facilitated fine-scale analysis of populations. However, estimation of SGS may depend strongly on the type of genetic marker and sampling strategy (of both loci and individuals). To explore sampling limits, we created a model population with simulated distributions of dominant and codominant alleles, resulting from natural regeneration with restricted gene flow. SGS estimates from subsamples (simulating collection and analysis with amplified fragment length polymorphism (AFLP) and microsatellite markers) were correlated with the 'real' estimate (from the full model population). For both marker types, sampling ranges were evident, with lower limits below which estimation was poorly correlated and upper limits above which sampling became inefficient. Lower limits (correlation of 0.9) were 100 individuals, 10 loci for microsatellites and 150 individuals, 100 loci for AFLPs. Upper limits were 200 individuals, five loci for microsatellites and 200 individuals, 100 loci for AFLPs. The limits indicated by simulation were compared with data sets from real species. Instances where sampling effort had been either insufficient or inefficient were identified. The model results should form practical boundaries for studies aiming to detect SGS. However, greater sample sizes will be required in cases where SGS is weaker than for our simulated population, for example, in species with effective pollen/seed dispersal mechanisms.

Genetic structure of a reef-building coral from thermally distinct environments on the Great Barrier Reef

Adaptation to localised thermal regimes is facilitated by restricted gene flow, ultimately leading to genetic divergence among populations and differences in their physiological tolerances. Allozyme analysis of six polymorphic loci was used to assess genetic differentiation between nine populations of the reef-building coral Acropora millepora over a latitudinal temperature gradient on the inshore regions of the Great Barrier Reef (GBR). Small but significant genetic differentiation indicative of moderate levels of gene flow (pairwise F-ST 0.023 to 0.077) was found between southern populations of A. millepora in cooler regions of the GBR and the warmer, central or northern GBR populations. Patterns of genetic differentiation at these putatively neutral allozyme loci broadly matched experimental variation in thermal tolerance and were consistent with local thermal regimes (warmest monthly-averages) for the A. millepora populations examined. It is therefore hypothesized that natural selection has influenced the thermal tolerance of the A. millepora populations examined and greater genetic divergence is likely to be revealed by examination of genetic markers under the direct effects of natural selection.

Environmental variables, habitat discontinuity and life history shaping the genetic structure of Pomatoschistus marmoratus

Coastal lagoons are semi-isolated ecosystems exposed to wide fluctuations of environmental conditions and showing habitat fragmentation. These features may play an important role in separating species into different populations, even at small spatial scales. In this study, we evaluate the concordance between mitochondrial (previous published data) and nuclear data analyzing the genetic variability of Pomatoschistus marmoratus in five localities, inside and outside the Mar Menor coastal lagoon (SE Spain) using eight microsatellites. High genetic diversity and similar levels of allele richness were observed across all loci and localities, although significant genic and genotypic differentiation was found between populations inside and outside the lagoon. In contrast to the FST values obtained from previous mitochondrial DNA analyses (control region), the microsatellite data exhibited significant differentiation among samples inside the Mar Menor and between lagoonal and marine samples. This pattern was corroborated using Cavalli-Sforza genetic distances. The habitat fragmentation inside the coastal lagoon and among lagoon and marine localities could be acting as a barrier to gene flow and contributing to the observed genetic structure. Our results from generalized additive models point a significant link between extreme lagoonal environmental conditions (mainly maximum salinity) and P. marmoratus genetic composition. Thereby, these environmental features could be also acting on genetic structure of coastal lagoon populations of P. marmoratus favoring their genetic divergence. The mating strategy of P. marmoratus could be also influencing our results obtained from mitochondrial and nuclear DNA. Therefore, a special consideration must be done in the selection of the DNA markers depending on the reproductive strategy of the species.

Small-scale genetic structure of Cerastoderma glaucum in a lagoonal environment: potential significance of habitat discontinuity and unstable population dynamics

Environmental heterogeneity in coastal lagoons is expected to facilitate local adaptation in response to different ecological conditions, causing significant genetic structuring within lagoon populations at a small scale and also differentiation between lagoons. However, these patterns and processes of genetic structuring are still poorly understood. The aims of our study were (1) to seek genetic structure at a small scale in Cerastoderma glaucum inside the Mar Menor coastal lagoon using a mitochondrial DNA marker (COI) that has previously detected genetic differentiation inside the lagoon in other species and (2) to evaluate the influence of extreme environmental conditions and habitat discontinuity on its genetic composition. The results indicate high levels of haplotype diversity and low values of nucleotide diversity. COI data provide evidence of significant population differentiation among some localities within the lagoon. Limited gene flow and unstable population dynamics (i.e. fluctuations in population size caused by local extinction and recolonization), probably due to the high environmental heterogeneity, could generate the small-scale genetic divergence detected between populations within the lagoon.

Connectivity patterns inferred from the genetic structure of white seabream (Diplodus sargus L.)

Themarine environment seems, at first sight, to be a homogeneousmediumlacking barriers to species dispersal. Nevertheless, populations of marine species show varying levels of gene flow and population differentiation, so barriers to gene flow can often be detected. Weaimto elucidate the role of oceanographical factors ingenerating connectivity among populations shaping the phylogeographical patterns in the marine realm, which is not only a topic of considerable interest for understanding the evolution ofmarine biodiversity but also formanagement and conservation of marine life. For this proposal,we investigate the genetic structure and connectivity between continental and insular populations ofwhite seabreamin North East Atlantic (NEA) and Mediterranean Sea (MS) aswell as the influence of historical and contemporary factors in this scenario using mitochondrial (cytochrome b) and nuclear (a set of 9 microsatellite) molecular markers. Azores population appeared genetically differentiated in a single cluster using Structure analysis. This result was corroborated by Principal Component Analysis (PCA) and Monmonier algorithm which suggested a boundary to gene flow, isolating this locality. Azorean population also shows the highest significant values of FST and genetic distances for both molecular markers (microsatellites and mtDNA). We suggest that the breakdown of effective genetic exchange between Azores and the others' samples could be explained simultaneously by hydrographic (deep water) and hydrodynamic (isolating current regimes) factors acting as barriers to the free dispersal of white seabream(adults and larvae) and by historical factors which could be favoured for the survival of Azorean white seabream population at the last glaciation. Mediterranean islands show similar genetic diversity to the neighbouring continental samples and nonsignificant genetic differences. Proximity to continental coasts and the current system could promote an optimal larval dispersion among Mediterranean islands (Mallorca and Castellamare) and coasts with high gene flow.

Association of kinase insert domain-containing receptor (KDR) gene polymorphism/ haplotypes with recurrent spontaneous abortion and genetic structure

Background: Recurrent spontaneous abortion is one of the diseases that can lead to physical, psychological, and, economical problems for both individuals and society. Recently a few numbers of genetic polymorphisms in kinase insert domain-containing receptor (KDR) gene are examined that can endanger the life of the fetus in pregnant women. Objective: The risk of KDR gene polymorphisms was investigated in Iranian women with idiopathic recurrent spontaneous abortion (RSA). Materials and Methods: A case controlled study was performed. One hundred idiopathic recurrent spontaneous abortion patients with at least two consecutive pregnancy losses before 20 weeks of gestational age with normal karyotypes were included in the study. Also, 100 healthy women with at least one natural pregnancy were studied as control group. Two functional SNPs located in KDR gene; rs1870377 (Q472H), and rs2305948 (V297I) as well as one tag SNP in the intron region (rs6838752) were genotyped by using PCR based restriction fragment length polymorphism (PCR-RFLP) technique. Haplotype frequency was determined for these three SNPs’ genotypes. Analysis of genetic STRUCTURE and K means clustering were performed to study genetic variation. Results: Functional SNP (rs1870377) was highly linked to tag SNP (rs6838752) (D´ value=0. 214; χ2 = 16.44, p<0. 001). K means clustering showed that k = 8 as the best fit for the optimal number of genetic subgroups in our studied materials. This result was in agreement with Neighbor Joining cluster analysis. Conclusion: In our study, the allele and genotype frequencies were not associated with RSA between patient and control individuals. Inconsistent results in different populations with different allele frequencies among RSA patients and controls may be due to ethnic variation and used sample size.

Population genetic structuring in Acanthopagrus butcheri (Pisces: Sparidae): does low gene flow among estuaries apply to both sexes?

Acanthopagrus butcheri completes its entire life history within estuaries and coastal lakes of southern Australia, although adults occasionally move between estuaries via the sea. Consequently, it is expected that populations of A. butcheri in different estuaries will be genetically distinct, with the magnitude of genetic divergence increasing with geographic isolation. However, previous genetic studies of A. butcheri from southeast Australia yielded conflicting results; allozyme variation exhibited minimal spatial structuring (θ = 0.012), whereas mitochondrial DNA distinguished the majority of populations analyzed (θ = 0.263) and genetic divergence was positively correlated with geographic isolation. This discrepancy could reflect high male gene flow, which impacts nuclear but not mitochondrial markers. Here we estimated allele frequencies at five nuclear microsatellite loci across 11 southeast Australian populations (595 individuals). Overall structuring of microsatellite variation was weaker (θ = 0.088) than that observed for mitochondrial DNA, but was able to distinguish a greater number of populations and was positively correlated with geographic distance. Therefore, we reject high male gene flow and invoke a stepping-stone model of infrequent gene flow among estuaries for both sexes. Likewise, management of A. butcheri within the study range should be conducted at the scale of individual or geographically proximate estuaries for both sexes. The lack of allozyme structuring in southeast Australia reflects either the large variance in structuring expected among loci under neutral conditions and the low number of allozymes surveyed or a recent colonization of estuaries such that some but not all nuclear loci have approached migration-drift equilibrium.

Population genetic structure of serotine bats (Eptesicus serotinus) across Europe and implications for the potential spread of bat rabies (European bat lyssavirus EBLV-1).

Understanding of the movements of species at multiple scales is essential to appreciate patterns of population connectivity and in some cases, the potential for pathogen transmission. The serotine bat (Eptesicus serotinus) is a common and widely distributed species in Europe where it frequently harbours European bat lyssavirus type 1 (EBLV-1), a virus causing rabies and transmissible to humans. In the United Kingdom, it is rare, with a distribution restricted to south of the country and so far the virus has never been found there. We investigated the genetic structure and gene flow of E. serotinus across the England and continental Europe. Greater genetic structuring was found in England compared with continental Europe. Nuclear data suggest a single population on the continent, although further work with more intensive sampling is required to confirm this, while mitochondrial sequences indicate an east-west substructure. In contrast, three distinct populations were found in England using microsatellite markers, and mitochondrial diversity was very low. Evidence of nuclear admixture indicated strong male-mediated gene flow among populations. Differences in connectivity could contribute to the high viral prevalence on the continent in contrast with the United Kingdom. Although the English Channel was previously thought to restrict gene flow, our data indicate relatively frequent movement from the continent to England highlighting the potential for movement of EBLV-1 into the United Kingdom.

Population genetic structure of the red fox (Vulpes vulpes) in the UK

The red fox (Vulpes vulpes) is common and widely distributed within the UK. It is a carrier or potential carrier of numerous zoonotic diseases. Despite this, there are no published reports on the population genetics of foxes in Britain. In this study, we aim to provide an insight into recent historical movement of foxes within Britain, as well as a current assessment of the genetic diversity and gene flow within British populations. We used 14 microsatellite markers to analyse 501 red fox samples originating from England, southern Scotland and northern France. High genetic diversity was evident within the sample set as a whole and limited population genetic structure was present in British samples analysed. Notably, STRUCTURE analysis found support of four population clusters, one of which grouped two southern England sampling areas with the nearby French samples from Calais, indicating recent (post-formation of the Channel) mixing of British and French populations. This may coincide with reports of large-scale translocations of foxes into Britain during the nineteenth century for sport hunting. Other STRUCTURE populations may be related to geographic features or to cultural practices such as fox hunting. In addition, the two British urban populations analysed showed some degree of differentiation from their local rural counterparts.

Genetic structure and sex-biased dispersal of a declining cooperative-breeder, the Grey-crowned Babbler, Pomatostomus temporalis, at the southern edge of its range

Loss and fragmentation of habitat can disrupt genetic exchange between populations, which is reflected in changes to the genetic structure of populations. The Grey-crowned Babbler (Pomatostomus temporalis) is a cooperatively breeding woodland bird, once common and widespread in south-eastern Australia. The species has suffered population declines of >90% across its southern distribution as a result of loss and fragmentation of habitat. We investigated patterns of genetic diversity and population structure of Grey-crowned Babblers in fragmented habitats at the southernmost extent of its range. We sampled blood from 135 individual Babblers from 39 groups stratified into six subpopulations in three regions. Genotypic data were used to estimate genetic diversity, population substructure, local relatedness and dispersal patterns. Individuals showed high heterozygosity within regions, and varying numbers of private alleles among regions suggested differences in levels of connectivity between regions. Four genetic clusters revealed population substructure consistent with treeless landscapes acting as strong barriers to gene flow. In contrast to previous studies,we identified a male-biased dispersal pattern and significant isolation-by-distance patterns for females at fine spatial scales. We recommend that conservation plans for this species incorporate opportunities to increase and enhance corridor areas to facilitate genetic exchange among subpopulations.

European lampreys: new insights on postglacial colonization, gene flow and speciation.

Ice ages are known to be the most dominant palaeoclimatic feature occurring on Earth, producing severe climatic oscillations and consequently shaping the distribution and the population structure of several species. Lampreys constitute excellent models to study the colonization of freshwater systems, as they commonly appear in pairs of closely related species of anadromous versus freshwater resident adults, thus having the ability to colonize new habitats, through the anadromous species, and establish freshwater resident derivates. We used 10 microsatellite loci to investigate the spatial structure, patterns of gene flow and migration routes of Lampetra populations in Europe. We sampled 11 populations including the migratory L. fluviatilis and four resident species, L. planeri, L. alavariensis, L. auremensis and L. lusitanica, the last three endemic to the Iberian Peninsula. In this southern glacial refugium almost all sampled populations represent a distinct genetic cluster, showing high levels of allopatric differentiation, reflecting long periods of isolation. As result of their more recent common ancestor, populations from northern Europe are less divergent among them, they are represented by fewer genetic clusters, and there is evidence of strong recent gene flow among populations. These previously glaciated areas from northern Europe may have been colonized from lampreys expanding out of the Iberian refugia. The pair L. fluviatilis/L. planeri is apparently at different stages of speciation in different locations, showing evidences of high reproductive isolation in the southern refugium, and low differentiation in the north.