Genetic structure and signature of population decrease in the critically endangered freshwater cyprinid <Emphasis Type="Italic">Chondrostoma lusitanicum</Emphasis>

The endemic and critically endangered cyprinid Chondrostoma lusitanicum has a very restricted distribution range. In order to estimate genetic diversity, characterize population structure and infer the demographic history, we examined six microsatellite loci and cytochrome b (mtDNA) sequences from samples taken throughout C. lusitanicum’s geographical range. Estimates of genetic diversity were low in all samples (average He < 0.35). The microsatellite data pointed to a major difference between northern (Samarra and Tejo drainages) and southern (Sado and Sines drainages) samples. This separation was not so clear with mtDNA, since one sample from the Tejo drainage grouped with the southern samples. This could be related with ancestral polymorphism or with admixture events between northern and southern sites during the late Pleistocene. Nevertheless, both markers indicate high levels of population differentiation in the north (for microsatellites F _ST >  0.23; and for mtDNA Φ_ST > 0.74) and lower levels in the south (F _ST < 0.05; Φ_ST < 0.40). With microsatellites we detected strong signals of a recent population decrease in effective size, by more than one order of magnitude, starting in the last centuries. This is consistent with field observations reporting a severe anthropogenic-driven population decline in the last decades. On the contrary mtDNA suggested a much older expansion. Overall, these results suggest that the distribution of genetic diversity in C. lusitanicum is the result of both ancient events related with drainage system formation, and recent human activities. The potential effect of population substructure generating genetic patterns similar to a population decrease is discussed, as well as the implications of these results for the conservation of C. lusitanicum. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Isolation and characterization of polymorphic microsatellite loci in Castanopsis chinensis Hance (Fagaceae)

Twelve novel microsatellite loci were isolated and characterized from enriched genomic libraries of Castanopsis chinensis. Four previously reported microsatellites from Castanopsis cuspidata were cross-amplified in C. chinensis. Forty-two sample trees from a wild population were tested for polymorphism using a set of the 16 polymorphic microsatellites. The average allele number of these microsatellites was 4.6 per locus, ranging from 2 to 7. The ranges of observed and expected heterozygosity were 0.262–1.000 and 0.238–0.818, respectively. Significant deviations from Hardy–Weinberg equilibrium were detected at five loci and no linkage disequilibrium was observed.     

Population genetic inferences using immune gene SNPs mirror patterns inferred by microsatellites

Single nucleotide polymorphisms (SNPs) are replacing microsatellites for population genetic analyses, but it is not apparent how many SNPs are needed or how well SNPs correlate with microsatellites. We used data from the gopher tortoise, Gopherus polyphemus—a species with small populations, to compare SNPs and microsatellites to estimate population genetic parameters. Specifically, we compared one SNP data set (16 tortoises from four populations sequenced at 17 901 SNPs) to two microsatellite data sets, a full data set of 101 tortoises and a partial data set of 16 tortoises previously genotyped at 10 microsatellites. For the full microsatellite data set, observed heterozygosity, expected heterozygosity and F_ST were correlated between SNPs and microsatellites; however, allelic richness was not. The same was true for the partial microsatellite data set, except that allelic richness, but not observed heterozygosity, was correlated. The number of clusters estimated by structure differed for each data set (SNPs = 2; partial microsatellite = 3; full microsatellite = 4). Principle component analyses (PCA) showed four clusters for all data sets. More than 800 SNPs were needed to correlate with allelic richness, observed heterozygosity and expected heterozygosity, but only 100 were needed for F_ST. The number of SNPs typically obtained from next‐generation sequencing (NGS) far exceeds the number needed to correlate with microsatellite parameter estimates. Our study illustrates that diversity, F_ST and PCA results from microsatellites can mirror those obtained with SNPs. These results may be generally applicable to small populations, a defining feature of endangered and threatened species, because theory predicts that genetic drift will tend to outweigh selection in small populations.     

Characterization of 14 novel microsatellite loci in the endangered <Emphasis Type="Italic">Liriodendron chinense</Emphasis> (Magnoliaceae) and cross-species amplification in closely related taxa

Fourteen polymorphic microsatellite loci were isolated and characterized from an AC-enriched genomic library of Liriodendron chinense (Magnoliaceae). The average allele number of these microsatellites was 5.2 per locus, ranging from three to eight. The observed and expected heterozygosity at population level were 0.07–1.00 and 0.10–0.83, respectively. These sets of microsatellites will be useful for studies of population genetic structure of L. chinense and L. tulipifera as well as to estimate fine-scale gene flow rates.     

Polymorphic microsatellites identified by cross-species amplifications in the European Coot Fulica atra

We tested the cross-amplification of 26 microsatellites developed for passerines and an additional three developed for Gallinula species in eight European Coots from two populations. Sixteen microsatellite markers successfully amplified, of which nine were polymorphic with 2–6 alleles (mean 3.7 alleles) and an expected heterozygosity (H _e) ranging from 0.375 to 0.805 (mean H _e = 0.589). On average, we found 2.22 alleles/locus and a mean H _e of 0.440 in one nest, and 2.56 alleles/locus and a mean H _e of 0.494 in the other one. These nine polymorphic markers could be of potential use in studies of genetic variability, population structure and reproductive strategy of European Coots.     

Development and characterization of microsatellite loci in the endangered Tricyrtis ishiiana (Convallariaceae), a local endemic plant in Japan

Nine microsatellite loci were isolated and characterized for the endangered and tertiary relict perennial herb, Tricyrtis ishiiana. The number of alleles ranged from 2 to 33. The expected (H _E) and observed (H _O) heterozygosities were 0.207–0.944 and 0.215–0.813, respectively, from 96 individuals on one population. Five loci exhibited significantly fewer heterozygotes than expected under Hardy–Weinberg equilibrium (P < 0.05). These primers amplifying microsatellites in this species may provide a useful tool for population genetics to establish conservation strategy.     

Isolation and characterisation of microsatellite loci in Calystegia soldanella (Convolvulaceae), an endangered coastal plant isolated in Lake Biwa,Japan

Eight microsatellite loci of Calystegia soldanella useful for comparisons of the genetic structure of isolated populations in the ancient Lake Biwa and coastal populations in Japan were isolated and characterised. The number of alleles ranged from 2 to 5. The expected (H _E) and observed (H _O) heterozygosities were 0.097–0.583 and 0.000–0.380, respectively, from 100 individuals from Lake Biwa and coastal populations. Seven of the eight loci exhibited significantly fewer heterozygotes than expected based on the Hardy–Weinberg equilibrium (P < 0.05). These primers amplifying microsatellites in C. soldanella may provide a population genetics tool useful in the establishment of a conservation strategy.     

Polymorphic microsatellites in largemouth bronze gudgeon (Coreius guichenoti) developed from repeat‐enriched libraries and cross‐species amplifications

Largemouth bronze gudgeon (Coreius guichenoti) is a medium‐sized fish endemic from the upper Yangtze River of China and its survival is threatened by the construction of the Three Gorges Dam. This study reports 20 new polymorphic microsatellites from a repeat‐enriched genomic library with a mean number allele of 5.2, and observed and expected heterozygosities ranging from 0.035 to 1, and from 0.13 to 0.917, respectively. In a cross‐species amplification test, nine of the 37 tested loci were found to be also polymorphic in a congeneric species, brass gudgeon (C. heterodon). In addition, other four loci from common carp (Cyprinus carpio) were also polymorphic in C. guichenoti. Out of these 24 polymorphic microsatellites, only three loci significantly deviated from Hardy–Weinberg equilibrium in the sampled population (P < 0.0025), and all pairwise tests for linkage disequilibrium among loci were nonsignificant after applying sequential Bonferroni correction (P > 0.0026). These novel microsatellites provide sufficient levels of polymorphism for studies on population genetics and conservation in C. guichenoti and its related species.     

Comparing RADseq and microsatellites to infer complex phylogeographic patterns,an empirical perspective in the Crucian carp,Carassius carassius,L.

The conservation of threatened species must be underpinned by phylogeographic knowledge. This need is epitomized by the freshwater fish Carassius carassius, which is in decline across much of its European range. Restriction site‐associated DNA sequencing (RADseq) is increasingly used for such applications; however, RADseq is expensive, and limitations on sample number must be weighed against the benefit of large numbers of markers. This trade‐off has previously been examined using simulation studies; however, empirical comparisons between these markers, especially in a phylogeographic context, are lacking. Here, we compare the results from microsatellites and RADseq for the phylogeography of C. carassius to test whether it is more advantageous to genotype fewer markers (microsatellites) in many samples, or many markers (SNPs) in fewer samples. These data sets, along with data from the mitochondrial cytochrome b gene, agree on broad phylogeographic patterns, showing the existence of two previously unidentified C. carassius lineages in Europe: one found throughout northern and central‐eastern European drainages and a second almost exclusively confined to the Danubian catchment. These lineages have been isolated for approximately 2.15 m years and should be considered separate conservation units. RADseq recovered finer population structure and stronger patterns of IBD than microsatellites, despite including only 17.6% of samples (38% of populations and 52% of samples per population). RADseq was also used along with approximate Bayesian computation to show that the postglacial colonization routes of C. carassius differ from the general patterns of freshwater fish in Europe, likely as a result of their distinctive ecology.     

Genetic diversity and gene flow in a Caribbean tree <Emphasis Type="Italic">Pterocarpus officinalis</Emphasis> Jacq.: a study based on chloroplast and nuclear microsatellites

We analysed the molecular diversity of Pterocarpus officinalis, a tree species distributed in Caribbean islands, South and Central America to quantify the genetic variation within island, to assess the pattern of differentiation and infer levels of gene flow; with the overall goal of defining a strategy of conservation. Two hundred two individuals of 9 populations were analysed using three chloroplast and six nuclear microsatellite markers. The observed heterozygosity varied markedly among the populations for nuclear (H _Onuc = 0.20–0.50) and chloroplast microsatellites (H _cp = 0.22–0.68). The continental population from French Guyana showed a higher value of H _Onuc than island populations, and the differences were significant in some cases. The fixation index F _IS ranged from −0.043 to 0.368; a significant heterozygote deficit was detected in 7 populations. The heterozygosity excess method suggested that two populations in Guadeloupe have undergone a recent bottleneck. Global and pairwise F _ST were high for both nuclear (F _STnuc = 0.29) and chloroplast microsatellites (F _STcp = 0.58). The neighbour-joining tree based on both markers, presented a differentiation pattern that can be explained by the seed dispersal by flotation and marine stream. The comparison of Bayesian approach and the method based on allelic frequency demonstrate a very limited number of migrants between populations.     

Isolation and characterization of polymorphic microsatellite loci in <Emphasis Type="Italic">Ardisia crenata</Emphasis> (Myrsinaceae)

Ardisia crenata, an evergreen shrub native to East Asia, has been a serious invasive plant to the southeastern USA. Here 13 polymorphic microsatellite loci were isolated and characterized from an enrichment genomic library of A. crenata. The average allele number of these microsatellites was four per locus, ranging from two to seven. The ranges of observed and expected heterozygosity were 0.000–1.000 and 0.239–0.789, respectively. These microsatellite markers will be useful for investigating population genetics and reproductive ecology of A. crenata.     

Morphological and molecular characterization of the Landes honey bee (<Emphasis Type="Italic">Apis mellifera</Emphasis> L.) ecotype for genetic conservation

A population of honey bees (Apis mellifera mellifera L.) with an annual colony brood cycle adapted to a locally abundant floral source in the Landes region of Southwest France is the subject of genetic conservation efforts. This population is maintained by local beekeepers in an area that experiences both an annual seasonal influx of non-local colonies and the permanent culture of imported stock. However, some colonies native to the Landes do not have the adapted brood cycle and their status as ecotypic are in question. The present study used morphology, mitochondrial DNA and microsatellites to characterize the endemic population and suggests further genetic conservation strategies. These methods yielded different degrees of discrimination of native and imported colonies and provided a powerful suite of tools for local resource managers. Colonies from the Landes could be differentiated from non-local French A. m. mellifera populations using morphometric analysis, and from non-native and reference populations using mtDNA and microsatellites. Seven morphological characters were identified by discriminant analysis as informative for delineating the Landes ecotype from other A. m. mellifera populations. Mitochondrial haplotypes for the population were characterized and five microsatellite loci were found to be informative in characterizing the Landes population. Asymmetric gene flow detected with microsatellite alleles was observed to be 5.5–5.9% from imported to native stocks of honey bees while introgression of native microsatellite alleles into imported colonies was 21.6%.     

Genetic signs of multiple colonization events in Baltic ciscoes with radiation into sympatric spring‐ and autumn‐spawners confined to early postglacial arrival

Presence of sympatric populations may reflect local diversification or secondary contact of already distinct forms. The Baltic cisco (Coregonus albula) normally spawns in late autumn, but in a few lakes in Northern Europe sympatric autumn and spring‐ or winter‐spawners have been described. So far, the evolutionary relationships and taxonomic status of these main life history forms have remained largely unclear. With microsatellites and mtDNA sequences, we analyzed extant and extinct spring‐ and autumn‐spawners from a total of 23 Swedish localities, including sympatric populations. Published sequences from Baltic ciscoes in Germany and Finland, and Coregonus sardinella from North America were also included together with novel mtDNA sequences from Siberian C. sardinella. A clear genetic structure within Sweden was found that included two population assemblages markedly differentiated at microsatellites and apparently fixed for mtDNA haplotypes from two distinct clades. All sympatric Swedish populations belonged to the same assemblage, suggesting parallel evolution of spring‐spawning rather than secondary contact. The pattern observed further suggests that postglacial immigration to Northern Europe occurred from at least two different refugia. Previous results showing that mtDNA in Baltic cisco is paraphyletic with respect to North American C. sardinella were confirmed. However, the inclusion of Siberian C. sardinella revealed a more complicated pattern, as these novel haplotypes were found within one of the two main C. albula clades and were clearly distinct from those in North American C. sardinella. The evolutionary history of Northern Hemisphere ciscoes thus seems to be more complex than previously recognized.     

Isolation and characterization of 14 polymorphic microsatellite loci in the ark shell Scapharca subcrenata (Bivalvia: Arcidae)

The ark shell Scapharca subcrenata is an important fishery resource, but has been suffering from severe population decline due to the deterioration of the coastal environment and over-exploitation. To investigate the genetic variation and structure among populations of S. subcrenata, 14 polymorphic microsatellite loci were developed. The number of alleles per polymorphic locus ranged from 11 to 29 and the observed and expected heterozygosity varied from 0.400 to 0.923 and from 0.705 to 0.965, respectively. These microsatellites will help advance the investigation of the genetic population structure and genetic diversity in this species.     

Characterization of microsatellites from cacao–<Emphasis Type="Italic">Moniliophthora perniciosa</Emphasis> interaction expressed sequence tags

Theobroma cacao L.–Moniliophthora perniciosa expressed sequence tags (ESTs) were converted into useful satellite markers for population analysis and genetic mapping. Forty-nine flanking primer pairs from TSH1188 (a resistant genotype) and Catongo (a susceptible genotype) ESTs were designed and screened for polymorphism analysis. Eleven were polymorphic, with an average of 3.81 alleles per locus and a total of 42 alleles. The satellite markers were tested on 21 cacao accessions and two bulked DNAs generated from 6 resistant and 6 susceptible plants from a segregating F₂ (SCA6 × ICS1) population for witches’ broom resistance. These results show that EST-derived microsatellites (short sequence repeats, SSRs) in Theobroma cacao have many potential applications in linkage mapping and the planning of crosses.     

Isolation and characterization of microsatellite loci in Shortia rotundifolia (Diapensiaceae), an endangered relict plant on the Ryukyu Islands and Taiwan

Ten microsatellite loci were isolated and characterized for Shortia rotundifolia, a relict temperate plant on the subtropical Ryukyu Islands and Taiwan. The number of alleles ranged from 2 to 21. The expected (H _E) and observed (H _O) heterozygosities were 0.027–0.877 and 0.027–0.622, respectively, from 74 individuals on Iriomotejima Island of Ryukyus. One locus exhibited significantly fewer heterozygotes than expected under Hardy–Weinberg equilibrium (P < 0.05). These primers amplifying microsatellites in this species may provide a useful tool for population genetics to establish conservation strategy.     

Depleted genetic variation of the European ground squirrel in Central Europe in both microsatellites and the major histocompatibility complex gene: implications for conservation

Habitat fragmentation may influence the genetic make-up and adaptability of endangered populations. To facilitate genetic monitoring of the endangered European ground squirrel (EGS), we analyzed 382 individuals from 16 populations in Central Europe, covering almost half of its natural range. We tested how fragmentation affects the genetic architecture of presumably selectively neutral (12 microsatellites) and non-neutral (the major histocompatibility class II DRB gene) loci. Spatial genetic analyses defined two groups of populations, “western” and “eastern”, with a significantly higher level of habitat fragmentation in the former group. The highly fragmented western populations had significantly lower genetic diversity in both types of markers. Only one allele of the DRB gene predominated in populations of the western group, while four alleles were evenly distributed across the eastern populations. Coefficient of inbreeding values (F _IS) calculated from microsatellites were significantly higher in the western (0.27–0.79) than in eastern populations (−0.060–0.119). Inter-population differentiation was very high, but similar in both groups (western F _ST = 0.23, eastern F _ST = 0.25). The test of isolation by distance was significant for the whole dataset, as well as for the two groups analyzed separately. Comparison of genetic variability and structure on microsatellites and the DRB gene does not provide any evidence for contemporary selection on MHC genes. We suggest that genetic drift in small bottlenecked and fragmented populations may overact the role of balancing selection. Based on the resulting risk of inbreeding depression in the western populations, we support population management by crossbreeding between the western and eastern populations.     

Genetic variability and the effect of habitat fragmentation in spotted suslik Spermophilus suslicus populations from two different regions

Endangered species worldwide exist in remnant populations, often within fragmented landscapes. Although assessment of genetic diversity in fragmented habitats is very important for conservation purposes, it is usually impossible to evaluate the amount of diversity that has actually been lost. Here, we compared population structure and levels of genetic diversity within populations of spotted suslik Spermophilus suslicus, inhabiting two different parts of the species range characterized by different levels of habitat connectivity. We used microsatellites to analyze 10 critically endangered populations located at the western part of the range, where suslik habitat have been severely devastated due to agriculture industrialization. Their genetic composition was compared with four populations from the eastern part of the range where the species still occupies habitat with reasonable levels of connectivity. In the western region, we detected extreme population structure (F _ST = 0.20) and levels of genetic diversity (Allelic richness ranged from 1.45 to 3.07) characteristic for highly endangered populations. Alternatively, in the eastern region we found significantly higher allelic richness (from 5.09 to 5.81) and insignificant population structure (F _ST = 0.03). As we identified a strong correlation between genetic and geographic distance and a lack of private alleles in the western region, we conclude that extreme population structure and lower genetic diversity is due to recent habitat loss. Results from this study provide guidelines for conservation and management of this highly endangered species.     

Microsatellite loci isolation in the endangered Gran Canarian blue chaffinch (Fringilla teydea polatzeki) and their utility in closely related taxa

The blue chaffinch, Fringilla teydea, is an endemic species of the Canary Islands. This species is formed by two subspecies: The Teneriffean blue chaffinch (F. t. teydea), and the endangered Gran Canarian blue chaffinch, (F. t. polatzeki). Here we report the isolation and characterization of nine tetranucleotide microsatellites (AAAG and AAAT) from the Gran Canarian subspecies, using an enrichment protocol. An average of 7.8 alleles per locus and an average observed heterozygosity of 0.773 were found (n = 28). The loci were tested for their ability to cross amplify in the Teneriffean subspecies and in the common chaffinch (Fringilla coelebs). These microsatellites will be used to manage a captive breeding programme for the endangered Gran Canarian subspecies.     

Genomewide gene‐associated microsatellite markers for the model invasive ascidian,Ciona intestinalis species complex

The vase tunicate, Ciona intestinalis species complex, has become a good model for ecological and evolutionary studies, especially those focusing on microevolution associated with rapidly changing environments. However, genomewide genetic markers are still lacking. Here, we characterized a large set of genomewide gene‐associated microsatellite markers for C. intestinalis spA (=C. robusta). Bioinformatic analysis identified 4654 microsatellites from expressed sequence tags (ESTs), 2126 of which successfully assigned to chromosomes were selected for further analysis. Based on the distribution evenness on chromosomes, function annotation and suitability for primer design, we chose 545 candidate microsatellites for further characterization. After amplification validation and variation assessment, 218 loci were polymorphic in at least one of the two populations collected from the coast of Arenys de Mar, Spain (N = 24–48), and Cape Town, South Africa (N = 24–33). The number of alleles, observed heterozygosity and expected heterozygosity ranged from 2 to 11, 0 to 0.833 and 0.021 to 0.818, and from 2 to 10, 0 to 0.879 and 0.031 to 0.845 for the Spanish and African populations, respectively. When all microsatellites were tested for cross‐species utility, only 60 loci (25.8%) could be successfully amplified and all loci were polymorphic in C. intestinalis spB. A high level of genomewide polymorphism is likely responsible for the low transferability. The large set of microsatellite markers characterized here is expected to provide a useful genomewide resource for ecological and evolutionary studies using C. intestinalis as a model.