Contrasting patterns of genetic diversity at three different genetic markers in a marine mammal metapopulation

Many studies use genetic markers to explore population structure and variability within species. However, only a minority use more than one type of marker and, despite increasing evidence of a link between heterozygosity and individual fitness, few ask whether diversity correlates with population trajectory. To address these issues, we analysed data from the Steller's sea lion, Eumetiopias jubatus , where three stocks are distributed over a vast geographical range and where both genetic samples and detailed demographic data have been collected from many diverse breeding colonies. To previously published mitochondrial DNA (mtDNA) and microsatellite data sets, we have added new data for amplified fragment length polymorphism (AFLP) markers, comprising 238 loci scored in 285 sea lions sampled from 23 natal rookeries. Genotypic diversity was low relative to most vertebrates, with only 37 loci (15.5%) being polymorphic. Moreover, contrasting geographical patterns of genetic diversity were found at the three markers, with Nei's gene diversity tending to be higher for AFLPs and microsatellites in rookeries of the western and Asian stocks, while the highest mtDNA values were found in the eastern stock. Overall, and despite strongly contrasting demographic histories, after applying phylogenetic correction we found little correlation between genetic diversity and either colony size or demography. In contrast, we were able to show a highly significant positive relationship between AFLP diversity and current population size across a range of pinniped species, even though equivalent analyses did not reveal significant trends for either microsatellites or mtDNA.     

Genetic diversity and differentiation in Eryngium alpinum L. (Apiaceae): comparison of AFLP and microsatellite markers

Genetic diversity and structure of 12 populations of Eryngium alpinum L. were investigated using 63 dominant amplified fragment length polymorphism (AFLP) and seven codominant microsatellite (48 alleles) markers. Within-population diversity estimates obtained with both markers were not correlated, but the microsatellite-based fixation index Fis was correlated with both AFLP diversity indices (number of polymorphic bands and Nei's expected heterozygosity). Only AFLP diversity indices increased with the size of populations, although they did not significantly differ among them (Kruskall-Wallis test). The discrepancy between AFLPs and microsatellites may be explained by a better coverage of the genome with numerous AFLPs, the higher mutation rates of microsatellites or the absence of significant difference among within-population diversity estimates. Genetic differentiation was higher with AFLPs (theta=0.40) than with microsatellites (theta=0.23), probably due to the higher polymorphism of microsatellites. Thus, we considered global qualitative patterns rather than absolute estimates to compare the performance of both types of markers. On a large geographic scale, the Mantel test and multivariate analysis showed that genetic patterns were more congruent with the spatial arrangement of populations when inferred from microsatellites than from AFLPs, suggesting higher homoplasy of AFLP markers. On a small spatial scale, AFLPs managed to discriminate individuals from neighboring populations whereas microsatellites did not (multivariate analysis), and the percentage of individuals correctly assigned to their population of origin was higher with AFLPs than with microsatellites. However, dominant AFLPs cannot be used to study heterozygosity-related topics. Thus, distinct molecular markers should be used depending on the biological question and the geographical scale investigated.     

Assessment of the domestication state of ackee (Blighia sapida K.D. Koenig) in Benin based on AFLP and microsatellite markers

Ackee (Blighia sapida) is a native multipurpose species important for the livelihoods of the rural populations in Benin. Trees are found in natural forests or are managed by farmers in different traditional agroforestry systems. Genetic variation at amplified fragment length polymorphism (AFLP) markers, four nuclear microsatellites (nSSRs) and one chloroplast microsatellite (cpSSR) were investigated in 279 individuals from six wild and eight cultivated populations from Benin. The AFLP data revealed moderate levels of diversity of ackee in Benin (mean diversity values are proportion of polymorphic loci = 52.8% and Nei??s gene diversity = 0.157, for 375 AFLP fragments). The mean diversity values based on nSSR-markers are expected heterozygosity = 0.286, allelic richness = 2.77. Genetic variation of wild and cultivated populations did not differ markedly. AMOVA revealed that only 7.3 and 5.2% of the variation was partitioned among populations for nSSR- and AFLP-markers, respectively. A Mantel test based on these both marker-types revealed significant correlations between population pairwise geographic distance and genetic differentiation. Differentiation among cultivated populations was higher than among wild populations. The only polymorphic chloroplast microsatellite marker (ccmp7) showed three haplotypes. Cultivated populations from northeastern Benin were fixed on one haplotype which was not observed elsewhere indicating a different origin of these populations possibly from neighboring Nigeria. Farmer-led domestication had an impact on the spatial distribution of genetic variation but did not result in significant losses of diversity within populations. Measures to conserve genetic resources of ackee in each of the three main bioclimatic zones in Benin are proposed.     

The origin of cultivated <Emphasis Type="Italic">Coffea arabica</Emphasis> L. varieties revealed by AFLP and SSR markers

Molecular markers were used to assess polymorphism between and within the genetic bases of coffee (i.e. Typica and Bourbon) spread from Yemen since the early 18th century that have given rise to most arabica cultivars grown world-wide. Eleven Coffea arabica accessions derived from the disseminated bases were evaluated by amplified fragment length polymorphism (AFLP) using 37 primer combinations and simple-sequence repeats (SSRs) produced by six microsatellites. Four cultivars growing in Yemen and 11 subspontaneous accessions collected in the primary centre of diversity of the species were included in the study in order to define their relationship with the accessions derived from the genetic bases of cultivars. One hundred and seven AFLP markers were used to calculate genetic distances and construct a dendrogram. The accessions derived from the disseminated bases were grouped separately, according to their genetic origin, and were distinguished from the subspontaneous accessions. The Yemen cultivars were classified with the Typica-derived accessions. Except for one AFLP marker, all AFLP and SSR markers present in the cultivated accessions were also detected in the subspontaneous accessions. Polymorphism among the subspontaneous accessions was much higher than among the cultivated accessions. It was very low within the genetic bases, confirming the historical documentation on their dissemination. The results enabled a discussion of the genetic diversity reductions that successively occurred during the dissemination of C. arabica from its primary centre of diversity.     

Microsatellite marker diversity in common bean (Phaseolus vulgaris L.)

A diversity survey was used to estimate allelic diversity and heterozygosity of 129 microsatellite markers in a panel of 44 common bean (Phaseolus vulgaris L.) genotypes that have been used as parents of mapping populations. Two types of microsatellites were evaluated, based respectively on gene coding and genomic sequences. Genetic diversity was evaluated by estimating the polymorphism information content (PIC), as well as the distribution and range of alleles sizes. Gene-based microsatellites proved to be less polymorphic than genomic microsatellites in terms of both number of alleles (6.0 vs. 9.2) and PIC values (0.446 vs. 0.594) while greater size differences between the largest and the smallest allele were observed for the genomic microsatellites than for the gene-based microsatellites (31.4 vs. 19.1 bp). Markers that showed a high number of alleles were identified with a maximum of 28 alleles for the marker BMd1. The microsatellites were useful for distinguishing Andean and Mesoamerican genotypes, for uncovering the races within each genepool and for separating wild accessions from cultivars. Greater polymorphism and race structure was found within the Andean gene pool than within the Mesoamerican gene pool and polymorphism rate between genotypes was consistent with genepool and race identity. Comparisons between Andean genotypes had higher polymorphism (53.0%) on average than comparisons among Mesoamerican genotypes (33.4%). Within the Mesoamerican parental combinations, the intra-racial combinations between Mesoamerica and Durango or Jalisco race genotypes showed higher average rates of polymorphism (37.5%) than the within-race combinations between Mesoamerica race genotypes (31.7%). In multiple correspondance analysis we found two principal clusters of genotypes corresponding to the Mesoamerican and Andean gene pools and subgroups representing specific races especially for the Nueva Granada and Peru races of the Andean gene pool. Intra population diversity was higher within the Andean genepool than within the Mesoamerican genepool and this pattern was observed for both gene-based and genomic microsatellites. Furthermore, intra-population diversity within the Andean races (0.356 on average) was higher than within the Mesoamerican races (0.302). Within the Andean gene pool, race Peru had higher diversity compared to race Nueva Granada, while within the Mesoamerican gene pool, the races Durango, Guatemala and Jalisco had comparable levels of diversity which were below that of race Mesoamerica.     

Isolation and characterization of microsatellites from a marine foodfish species ribbonfish Trichiurus haumela

The ribbonfish, Trichiurus haumela, is an important marine foodfish species. Supply of this fish relies on capture from the wild. No information about genetic diversity and population structure of wild stocks is available. Therefore, we developed the first set of 10 polymorphic microsatellites. The average allele number was 13.2 per locus ranging from six to 20. The expected heterozygosity was 0.84, whereas the observed heterozygosity was 0.71. This set of microsatellites provides a tool for assessing genetic diversity and population structure of this species.     

Genetic diversity analysis using lowly polymorphic dominant markers: the example of AFLP in pigs

DNA markers are commonly used for large-scale evaluation of genetic diversity in farm animals, as a component of the management of animal genetic resources. AFLP markers are useful for such studies as they can be generated relatively simply; however, challenges in analysis arise from their dominant scoring and the low level of polymorphism of some markers. This paper describes the results obtained with a set of AFLP markers in a study of 59 pig breeds. AFLP fingerprints were generated using four primer combinations (PC), yielding a total of 148 marker loci, and average harmonic mean of breed sample size was 37.3. The average proportion of monomorphic populations was 63% (range across loci: 3%-98%). The moment-based method of Hill and Weir (2004, Mol Ecol 13:895-908) was applied to estimate gene frequencies, gene diversity (F(ST)), and Reynolds genetic distances. A highly significant average F(ST) of 0.11 was estimated, together with highly significant PC effects on gene diversity. The variance of F(ST) across loci also significantly exceeded the variance expected under the hypothesis of AFLP neutrality, strongly suggesting the sensitivity of AFLP to selection or other forces. Moment estimates were compared to estimates derived from the square root estimation of gene frequency, as currently applied for dominant markers, and the biases incurred in the latter method were evaluated. The paper discusses the hypotheses underlying the moment estimations and various issues relating to the biallelic, dominant, and lowly polymorphic nature of this set of AFLP markers and to their use as compared to microsatellites for measuring genetic diversity.     

Global patterns reveal strong population structure in Haemonchus contortus, a nematode parasite of domesticated ruminants

We have examined the global population genetic structure of Haemonchus contortus. The genetic variability was studied using both amplified fragment length polymorphism (AFLP) and nad4 sequences of the mitochondrial genome. To examine the performance and information content of the two different marker systems, comparative assessment of population genetic diversity was undertaken in 19 isolates of H. contortus, a parasitic nematode of small ruminants. A total of 150 individual adult worms representing 14 countries from all inhabited continents were analysed. Altogether 1,429 informative AFLP markers were generated using four different primer combinations. Also, the genetic variation was high, which agrees with results from previous AFLP studies of nematode parasites of livestock. The genetic structure was high, indicating limited gene flow between the different isolates and populations from each continent mostly formed monophyletic groups in the phylogenetic analysis. However, for isolates representing Australia, Greece and one laboratory strain that originated from South Africa (WRS), there was no clear genetic relationship between the isolates and the distance between their geographical origins. Basically the same pattern was observed for the mitochondrial marker, although the phylogenetic analysis was less resolved than for AFLP. In contrast with previous findings on the population genetic structure of H. contortus, the calculation of population structure gave high values (Nst=0.59). The strong structure was present also for the four Swedish isolates (Nst=0.16) representing a small geographical area.     

Genetic diversity and population structure of Lantana camara in India indicates multiple introductions and gene flow

Lantana camara is a highly invasive plant, which has spread over 60 countries and island groups of Asia, Africa and Australia. In India, it was introduced in the early nineteenth century, since when it has expanded and gradually established itself in almost every available ecosystem. We investigated the genetic diversity and population structure of this plant in India in order to understand its introduction, subsequent range expansion and gene flow. A total of 179 individuals were sequenced at three chloroplast loci and 218 individuals were genotyped for six nuclear microsatellites. Both chloroplasts (nine haplotypes) and microsatellites (83 alleles) showed high genetic diversity. Besides, each type of marker confirmed the presence of private polymorphism. We uncovered low to medium population structure in both markers, and found a faint signal of isolation by distance with microsatellites. Bayesian clustering analyses revealed multiple divergent genetic clusters. Taken together, these findings (i.e. high genetic diversity with private alleles and multiple genetic clusters) suggest that Lantana was introduced multiple times and gradually underwent spatial expansion with recurrent gene flow.     

Concordant genetic estimators of migration reveal anthropogenically enhanced source-sink population structure in the river sculpin, Cottus gobio

River systems are vulnerable to natural and anthropogenic habitat fragmentation and will often harbor populations deviating markedly from simplified theoretical models. We investigated fine-scale population structure in the sedentary river fish Cottus gobio using microsatellites and compared migration estimates from three FST estimators, a coalescent maximum-likelihood method and Bayesian recent migration analyses. Source-sink structure was evident via asymmetry in migration and genetic diversity with smaller upstream locations emigration biased and larger downstream subpopulations immigration biased. Patterns of isolation by distance suggested that the system was largely, but not entirely, in migration-drift equilibrium, with headwater populations harboring a signal of past colonizations and in some cases also recent population bottlenecks. Up- vs. downstream asymmetry in population structure was partly attributable to the effects of flow direction, but was enhanced by weirs prohibiting compensatory upstream migration. Estimators of migration showed strong correspondence, at least in relative terms, especially if pairwise FST was used as an indirect index of relative gene flow rather than being translated to Nm. Since true parameter values are unknown in natural systems, comparisons among estimators are important, both to determine confidence in estimates of migration and to validate the performance of different methods.     

Inheritance of amplified fragment length polymorphism markers and their utility in population genetic analysis of Plecoglossus altivelis

The reproducibility, mode of inheritance and polymorphism of amplified fragment length polymorphism (AFLP) markers were examined in ayu Plecoglossus altivelis (Salmoniformes: Plecoglossidae). The AFLP markers were highly reproducible, their inheritance following Mendelian expectations. The number of fragments amplified (34–134), polymorphic ratio (0·15–0·78) and average heterozygosity (0·02–0·25) of the AFLP markers showed significant variation among six primer pairs and among ayu populations, including a landlocked Lake‐Biwa population, two amphidromous populations ( P. a. altivelis ) and two Ryukyu‐ayu populations ( P. a. ryukyuensis ). Although AFLP analysis provided similar results in intra‐population diversity and relationships among populations to those found by analyses of allozymes, microsatellites and mitochondrial DNA sequences, AFLPs showed higher polymorphisms and hence greater distinction between genetically close populations.     

Waste not,want not: Microsatellites remain an economical and informative technology for conservation genetics

Comparisons of microsatellites and single‐nucleotide polymorphisms (SNPs) have found that SNPs outperform microsatellites in population genetic analyses, questioning the continued utility of microsatellites in population and landscape genetics. Yet, highly polymorphic markers may be of value in species that have reduced genetic variation. This study repeated previous analyses that used microsatellites with SNPs developed from ddRAD sequencing in the black‐capped vireo source‐sink system. SNPs provided greater resolution of genetic diversity, population differentiation, and migrant detection but could not reconstruct parentage relationships due to insufficient heterozygosities. The biological inferences made by both sets of markers were similar: asymmetrical gene flow from source sites to the remaining sink sites. With the landscape genetic analyses, we found different results between the two molecular markers, but associations of the top environmental features (riparian, open habitat, agriculture, and human development) with dispersal estimates were shared between marker types. Despite the higher precision of SNPs, we find that microsatellites effectively uncover population processes and patterns and are superior for parentage analyses in this species with reduced genetic diversity. This study illustrates the continued applicability and relevance of microsatellites in population genetic research.     

Comparative Assessment of SSR and AFLP Markers for Evaluation of Genetic Diversity and Conservation of Fig, <Emphasis Type="Italic">Ficus carica</Emphasis> L., Genetic Resources in Tunisia

This study characterises the genetic variability of fig, Ficus carica L., using simple sequence repeat (SSR) and amplified fragment length polymorphism (AFLP) markers. It compares the efficiency and utility of the two techniques in detecting variation and establishing genetic relationships among Tunisian fig cultivars. Our results show that using both marker systems, the Tunisian fig germ plasm is characterised by having a large genetic diversity at the deoxyribonucleic acid level, as most of AFLP bands were detected and all SSR markers were polymorphic. In fact, 351 (342 polymorphic) and 57 (57 polymorphic) bands were detected using AFLP and SSR primers, respectively. SSR markers were the most polymorphic with an average polymorphic information content value of 0.94, while AFLP markers showed the highest effective multiplex ratio (56.9) and marker index (45.2). The effective marker index was recorded highest (4.19) for AFLP markers and lowest (0.70) for the SSR ones. Our results demonstrate that (1) independent as well as combined analyses of cluster analyses of SSR and AFLP fragments showed that cultivars are clustered independently from their geographical origin, horticultural classifications and tree sex; (2) the analysis of molecular variance allowed the partitioning of genetic variation within and among fig groups and showed greater variation within groups and (3) AFLP and SSR markers datasets showed positive correlation. This study suggests the SSR and AFLP markers are suitable for diversity analysis and cultivars fingerprinting. An understanding of the genetic diversity and population structure of F. carica in Tunisia can also provide insight into the conservation and management of this species.     

Genetic evidence of local adaptation of wheat yellow rust (Puccinia striiformis f. sp. tritici) within France

Puccinia striiformis f. sp. tritici (PST), a clonal basidiomycete causing yellow rust disease on wheat, has a long record of 'overcoming' cultivar resistance introduced by breeders. Despite the long dispersal capacity of its spores, the French population of PST presents a strong geographical structure, with the presence of a specific pathotype (array of avirulence genes) at high frequencies in the south of France. The genetic diversity underlying this differentiation was analysed by microsatellite and AFLP markers. A total of 213 French isolates belonging to 10 pathotypes collected over a 15-year period were investigated. For each of the 12 microsatellites used, polymorphism resulted from a unique allelic variant associated to the south-specific pathotype. This pathotype was characterized by 40 specific markers over the total of 63 polymorphins detected using 15 AFLP primer combinations. Phylogeographical analysis indicated a strictly clonal structure of the population, and a strong genomic divergence between the northern population and a south-specific clone. Both virulence and molecular data show that the northern French population belongs to the northwestern European population, whereas the southern clone is most likely related to a Mediterranean population, the two subpopulations resulting from the ancient divergence of two clonal lineages. While the virulence complexity in the northern population may be explained by the successive introduction of corresponding resistance genes in cultivars, the maintenance of a simple virulence type in southern France, despite gene flow between the two populations, may be explained in terms of host cultivars repartition and local adaptation to specific host or climatic conditions.     

Novel polymorphic microsatellites for studying genetic diversity of red Asian arowanas

Asian arowana (Scleropages formosus) is a highly endangered fish species listed in Appendix 1 of CITES since 1980. Fourteen novel polymorphic microsatellites were isolated from a CA-enriched partial genomic DNA library, and were used for studying genetic diversity of 41 red arowanas from the wild population. Surprisingly the average allele number of the 14 polymorphic microsatellites was as high as 12.4/locus. The average observed heterozygosity was 0.78 ranging from 0.51 to 0.95, and the gene diversity was quite high (0.78). All these data suggest that high level of genetic diversity existed in the red Asian arowana population.     

Conservation genetics of endangeredBrasenia schreberi based on RAPD and AFLP markers

Brasenia schreberi J.F. Gmelin is a declared endangered species found in the lakes and ponds of South Korea. For planning its conservation strategy, we examined the genetic diversity within and among six populations, using randomly amplified polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP). Polymorphisms were more frequently detected per loci with AFLP (69.3%) than RAPD (36.8%). High genetic diversity was recognized within populations: polymorphic loci (PPL) values ranged from 36.3% in the CJM population to 74.5% in the GGT population, with a mean value of 47.8% based on AFLP markers. Great genetic differentiation (θ^B) was detected among the six populations (0.670 on RAPD and 0.196 on AFLP), and we calculated a low rate of gene flow (N_em), i.e., 0.116 on RAPD and 0.977 on AFLP. Furthermore, a Mantel test revealed that no correlation existed between genetic distances and geographical distances among the six local populations, based on RAPD or AFLP markers. These results are attributed to a number of factors, including an insufficient length of time for genetic diversity to be reduced following a natural decline in population size and isolation, adaptation of the genetic system to small population conditions, and a restricted gene flow rate. Based on both its genetic diversity and population structure, we suggest that a strategy for conserving and restoringB. schreberi must focus on maintaining historical processes, such as high levels of outbreeding, while monitoring increased gene flow among populations. This is because a reduction in genetic diversity as a result of genetic drift is undesirable.     

Population Genetic Studies Revealed Local Adaptation in a High Gene-Flow Marine Fish,the Small Yellow Croaker (Larimichthys polyactis)

The genetic differentiation of many marine fish species is low. Yet local adaptation may be common in marine fish species as the vast and changing marine environment provides more chances for natural selection. Here, we used anonymous as well as known protein gene linked microsatellites and mitochondrial DNA to detect the population structure of the small yellow croaker (Larimichthys polyactis) in the Northwest Pacific marginal seas. Among these loci, we detected at least two microsatellites, anonymous H16 and HSP27 to be clearly under diversifying selection in outlier tests. Sequence cloning and analysis revealed that H16 was located in the intron of BAHCC1 gene. Landscape genetic analysis showed that H16 mutations were significantly associated with temperature, which further supported the diversifying selection at this locus. These marker types presented different patterns of population structure: (i) mitochondrial DNA phylogeny showed no evidence of genetic divergence and demonstrated only one glacial linage; (ii) population differentiation using putatively neutral microsatellites presented a pattern of high gene flow in the L. polyactis. In addition, several genetic barriers were identified; (iii) the population differentiation pattern revealed by loci under diversifying selection was rather different from that revealed by putatively neutral loci. The results above suggest local adaptation in the small yellow croaker. In summary, population genetic studies based on different marker types disentangle the effects of demographic history, migration, genetic drift and local adaptation on population structure and also provide valuable new insights for the design of management strategies in L. polyactis.     

Isolation and characterization of dinucleotide microsatellites in greater amberjack, <Emphasis Type="Italic">Seriola dumerili</Emphasis>

Thirteen nuclear-encoded dinucleotide microsatellites were characterized from a genomic DNA library of greater amberjack, Seriola dumerili. The microsatellites include 12 perfect-repeat motifs and one imperfect-repeat motif. The number of alleles at the 13 microsatellites among a sample of 29 fish ranged from 3 to 25; gene diversity (expected heterozygosity) ranged from 0.296 to 0.948, while observed heterozygosity ranged from 0.276 to 0.897. Following Bonferroni correction, genotypes at all 13 microsatellites fit expectations of Hardy-Weinberg equilibrium. One pairwise comparison of microsatellites deviated significantly from expectations of genotypic equilibrium, suggesting that these two microsatellites may be linked. Greater amberjack support commercial and recreational fisheries along both the Atlantic and Gulf coasts of the U.S. and represent a species with potential for worldwide aquaculture. The microsatellites developed will be useful for conservation and population genetic studies of ‘wild’ and domesticated populations of greater amberjack.     

Comparing RADseq and microsatellites for estimating genetic diversity and relatedness — Implications for brown trout conservation

The conservation and management of endangered species requires information on their genetic diversity, relatedness and population structure. The main genetic markers applied for these questions are microsatellites and single nucleotide polymorphisms (SNPs), the latter of which remain the more resource demanding approach in most cases. Here, we compare the performance of two approaches, SNPs obtained by restriction‐site‐associated DNA sequencing (RADseq) and 16 DNA microsatellite loci, for estimating genetic diversity, relatedness and genetic differentiation of three, small, geographically close wild brown trout (Salmo trutta) populations and a regionally used hatchery strain. The genetic differentiation, quantified as F_ST, was similar when measured using 16 microsatellites and 4,876 SNPs. Based on both marker types, each brown trout population represented a distinct gene pool with a low level of interbreeding. Analysis of SNPs identified half‐ and full‐siblings with a higher probability than the analysis based on microsatellites, and SNPs outperformed microsatellites in estimating individual‐level multilocus heterozygosity. Overall, the results indicated that moderately polymorphic microsatellites and SNPs from RADseq agreed on estimates of population genetic structure in moderately diverged, small populations, but RADseq outperformed microsatellites for applications that required individual‐level genotype information, such as quantifying relatedness and individual‐level heterozygosity. The results can be applied to other small populations with low or moderate levels of genetic diversity.