An objective, rapid and reproducible method for scoring AFLP peak-height data that minimizes genotyping error

Amplified fragment length polymorphism (AFLP) fingerprint data are now commonly collected using DNA sequencers. AFLP genotypes are still often scored by eye from such data - a time-consuming, error-prone and subjective process. We present a semi-automated method of genotyping sequencer-collected AFLPs at predefined fragment locations (loci) within the fingerprint. Our method uses thresholds of AFLP-polymerase chain reaction-product fluorescence intensity (peak height) in order to: (i) exclude AFLP loci that are likely to contribute high rates of error to data sets, and (ii) determine the AFLP phenotype (fragment presence or absence) at the retained loci. Error rate analysis is an integral part of this process and is used to determine optimal thresholds that minimize genotyping error, while maximizing the numbers of retained loci. We show that application of this method to a large AFLP data set allows genotype calls that are rapid, objective and repeatable, facilitating the extraction of reliable genotype data for molecular ecological studies.     

How to track and assess genotyping errors in population genetics studies

Genotyping errors occur when the genotype determined after molecular analysis does not correspond to the real genotype of the individual under consideration. Virtually every genetic data set includes some erroneous genotypes, but genotyping errors remain a taboo subject in population genetics, even though they might greatly bias the final conclusions, especially for studies based on individual identification. Here, we consider four case studies representing a large variety of population genetics investigations differing in their sampling strategies (noninvasive or traditional), in the type of organism studied (plant or animal) and the molecular markers used [microsatellites or amplified fragment length polymorphisms (AFLPs)]. In these data sets, the estimated genotyping error rate ranges from 0.8% for microsatellite loci from bear tissues to 2.6% for AFLP loci from dwarf birch leaves. Main sources of errors were allelic dropouts for microsatellites and differences in peak intensities for AFLPs, but in both cases human factors were non-negligible error generators. Therefore, tracking genotyping errors and identifying their causes are necessary to clean up the data sets and validate the final results according to the precision required. In addition, we propose the outline of a protocol designed to limit and quantify genotyping errors at each step of the genotyping process. In particular, we recommend (i) several efficient precautions to prevent contaminations and technical artefacts; (ii) systematic use of blind samples and automation; (iii) experience and rigor for laboratory work and scoring; and (iv) systematic reporting of the error rate in population genetics studies.     

Genome wide AFLP markers support cryptic species in Coniophora (Boletales)

Numerous fungal morphospecies include cryptic species that routinely are detected by sequencing a few unlinked DNA loci. However, whether the patterns observed by multi-locus sequencing are compatible with genome wide data, such as amplified fragment length polymorphisms (AFLPs), is not well known for fungi. In this study we compared the ability of three DNA loci and AFLP data to discern between cryptic fungal lineages in the three morphospecies Coniophora olivacea, Coniophora arida, and Coniophora puteana. The sequences and the AFLP data were highly congruent in delimiting the morphotaxa into multiple cryptic species. However, while the DNA sequences indicated introgression or hybridization between some of the cryptic lineages the AFLP data did not. We conclude that as few as three polymorphic DNA loci was sufficient to recognize cryptic lineages within the studied Coniophora taxa. However, based on analyses of a few (three) sequenced loci the hybridization could not easily be distinguished from incomplete lineage sorting. Hence, great caution should be taken when concluding about hybridization based on data from just a few loci.     

Genetic Analysis of Selected Strains of Eastern Oyster (Crassostrea virginica Gmelin) Using AFLP and Microsatellite Markers

Amplified fragment length polymorphisms (AFLPs) and microsatellite markers were used to examine genetic variation and divergence in 4 selected strains (DBH, NEH, FMF, and CTS) and 1 wild population (DBW) of the eastern oyster Crassostrea virginica Gmelin. Eighty-six AFLP markers (from 3 primer pairs) and 5 microsatellite loci were used for the analysis of 30 oysters from each of the 5 populations. Microsatellite loci were considerably more variable than AFLPs. The observed heterozygosity ranged from 0.560 to 0.640 across populations for microsatellites, and from 0.186 to 0.207 for AFLPs. Both F_st and _PT of microsatellite data and _PT statistics of AFLP data revealed significant divergence between all pairs of populations. There was no significant reduction in heterozygosity in all 4 selected strains; however, the number of alleles per locus was considerably lower in the selected strains than in the wild population. Two strains subjected to long-term selection for disease resistance shared frequency shifts at a few loci, which deserve further analysis to determine if they are linked to disease-resistance genes.     

Comparative analysis of genetic diversity in the mangrove species Avicennia marina (Forsk.) Vierh. (Avicenniaceae) detected by AFLPs and SSRs

Avicennia marina is an important mangrove species with a wide geographical and climatic distribution which suggests that large amounts of genetic diversity are available for conservation and breeding programs. In this study we compare the informativeness of AFLPs and SSRs for assessing genetic diversity within and among individuals, populations and subspecies of A. marina in Australia. Our comparison utilized three SSR loci and three AFLP primer sets that were known to be polymorphic, and could be run in a single analysis on a capillary electrophoresis system, using different- colored fluorescent dyes. A total of 120 individuals representing six populations and three subspecies were sampled. At the locus level, SSRs were considerably more variable than AFLPs, with a total of 52 alleles and an average heterozygosity of 0.78. Average heterozygosity for AFLPs was 0.193, but all of the 918 bands scored were polymorphic. Thus, AFLPs were considerably more efficient at revealing polymorphic loci than SSRs despite lower average heterozygosities. SSRs detected more genetic differentiation between populations (19 vs 9%) and subspecies (35 vs 11%) than AFLPs. Principal co-ordinate analysis revealed congruent patterns of genetic relationships at the individual, population and subspecific levels for both data sets. Mantel testing confirmed congruence between AFLP and SSR genetic distances among, but not within, population comparisons, indicating that the markers were segregating independently but that evolutionary groups (populations and subspecies) were similar. Three genetic criteria of importance for defining priorities for ex situ collections or in situ conservation programs (number of alleles, number of locally common alleles and number of private alleles) were correlated between the AFLP and SSR data sets. The congruence between AFLP and SSR data sets suggest that either method, or a combination, is applicable to expanded genetic studies of mangroves. The codominant nature of SSRs makes them ideal for further population-based investigations, such as mating-system analyses, for which the dominant AFLP markers are less well suited. AFLPs may be particularly useful for monitoring propagation programs and identifying duplicates within collections, since a single PCR assay can reveal many loci at once. Received: 3 October 2000 / Accepted: 19 February 2001     

Genetic diversity in cultivated carioca common beans based on molecular marker analysis

A wide array of molecular markers has been used to investigate the genetic diversity among common bean species. However, the best combination of markers for studying such diversity among common bean cultivars has yet to be determined. Few reports have examined the genetic diversity of the carioca bean, commercially one of the most important common beans in Brazil. In this study, we examined the usefulness of two molecular marker systems (simple sequence repeats - SSRs and amplified fragment length polymorphisms - AFLPs) for assessing the genetic diversity of carioca beans. The amount of information provided by Roger's modified genetic distance was used to analyze SSR data and Jaccards similarity coefficient was used for AFLP data. Seventy SSRs were polymorphic and 20 AFLP primer combinations produced 635 polymorphic bands. Molecular analysis showed that carioca genotypes were quite diverse. AFLPs revealed greater genetic differentiation and variation within the carioca genotypes (Gst = 98% and Fst = 0.83, respectively) than SSRs and provided better resolution for clustering the carioca genotypes. SSRs and AFLPs were both suitable for assessing the genetic diversity of Brazilian carioca genotypes since the number of markers used in each system provided a low coefficient of variation. However, fingerprint profiles were generated faster with AFLPs, making them a better choice for assessing genetic diversity in the carioca germplasm.     

Genetic diversity and population structure inferred from the partially duplicated genome of domesticated carp,Cyprinus carpio L.

Genetic relationships among eight populations of domesticated carp (Cyprinus carpio L.), a species with a partially duplicated genome, were studied using 12 microsatellites and 505 AFLP bands. The populations included three aquacultured carp strains and five ornamental carp (koi) variants. Grass carp (Ctenopharyngodon idella) was used as an outgroup. AFLP-based gene diversity varied from 5% (grass carp) to 32% (koi) and reflected the reasonably well understood histories and breeding practices of the populations. A large fraction of the molecular variance was due to differences between aquacultured and ornamental carps. Further analyses based on microsatellite data, including cluster analysis and neighbor-joining trees, supported the genetic distinctiveness of aquacultured and ornamental carps, despite the recent divergence of the two groups. In contrast to what was observed for AFLP-based diversity, the frequency of heterozygotes based on microsatellites was comparable among all populations. This discrepancy can potentially be explained by duplication of some loci in Cyprinus carpio L., and a model that shows how duplication can increase heterozygosity estimates for microsatellites but not for AFLP loci is discussed. Our analyses in carp can help in understanding the consequences of genotyping duplicated loci and in interpreting discrepancies between dominant and co-dominant markers in species with recent genome duplication.     

Comparative multi-locus phylogeography confirms multiple vicariance events in co-distributed rainforest frogs

Though Pleistocene refugia are frequently cited as drivers of species diversification, comparisons of molecular divergence among sister species typically indicate a continuum of divergence times from the Late Miocene, rather than a clear pulse of speciation events at the Last Glacial Maximum. Community-scale inference methods that explicitly test for multiple vicariance events, and account for differences in ancestral effective population size and gene flow, are well suited for detecting heterogeneity of species' responses to past climate fluctuations. We apply this approach to multi-locus sequence data from five co-distributed frog species endemic to the Wet Tropics rainforests of northeast Australia. Our results demonstrate at least two episodes of vicariance owing to climate-driven forest contractions: one in the Early Pleistocene and the other considerably older. Understanding how repeated cycles of rainforest contraction and expansion differentially affected lineage divergence among co-distributed species provides a framework for identifying evolutionary processes that underlie population divergence and speciation.     

Widespread amplification of amplified fragment length polymorphisms (AFLPs) in marine Antarctic animals

Although recent years have witnessed a rapid growth in the number of genetic studies of Antarctic organisms, relatively few studies have so far used nuclear markers, possibly due to the perceived cost and difficulty of isolating markers such as microsatellites. However, an often overlooked alternative is to use amplified fragment length polymorphisms (AFLPs), a versatile and low-cost method capable of generating large numbers of predominantly nuclear loci in virtually any organism. We conducted a literature review of population genetic studies of Antarctic organisms, finding that fewer than 10% used AFLPs. Moreover, a strong taxonomic bias was found, with studies employing mitochondrial DNA or microsatellites focussing predominantly on animals, while those using AFLPs were mostly of plants or lower organisms. Consequently, we explored the extent to which AFLPs amplify across a range of Antarctic marine animal taxa by genotyping eight individuals each of twelve different species, ranging from echinoderms through soft corals to pelagic fish, at four selective primer combinations. AFLPs readily amplified across all of the taxa, generating between 32 and 84 loci per species, with on average 56.5% of these being polymorphic. In general, levels of polymorphism bore little relationship with expectations based on larger populations of broadcast-spawning species being more variable, though we did find a tentative positive correlation between the number of AFLP bands amplified and a measure of effective population size. Our study lends further support for the utility and ease of use of AFLPs and their suitability for studies of Antarctic species across a wide range of taxa.     

How do SNP ascertainment schemes and population demographics affect inferences about population history?

Background

The selection of variable sites for inclusion in genomic analyses can influence results, especially when exemplar populations are used to determine polymorphic sites. We tested the impact of ascertainment bias on the inference of population genetic parameters using empirical and simulated data representing the three major continental groups of cattle: European, African, and Indian. We simulated data under three demographic models. Each simulated data set was subjected to three ascertainment schemes: (I) random selection; (II) geographically biased selection; and (III) selection biased toward loci polymorphic in multiple groups. Empirical data comprised samples of 25 individuals representing each continental group. These cattle were genotyped for 47,506 loci from the bovine 50 K SNP panel. We compared the inference of population histories for the empirical and simulated data sets across different ascertainment conditions using F_ST and principal components analysis (PCA).

Results

Bias toward shared polymorphism across continental groups is apparent in the empirical SNP data. Bias toward uneven levels of within-group polymorphism decreases estimates of F_ST between groups. Subpopulation-biased selection of SNPs changes the weighting of principal component axes and can affect inferences about proportions of admixture and population histories using PCA. PCA-based inferences of population relationships are largely congruent across types of ascertainment bias, even when ascertainment bias is strong.

Conclusions

Analyses of ascertainment bias in genomic data have largely been conducted on human data. As genomic analyses are being applied to non-model organisms, and across taxa with deeper divergences, care must be taken to consider the potential for bias in ascertainment of variation to affect inferences. Estimates of F_ST, time of separation, and population divergence as estimated by principal components analysis can be misleading if this bias is not taken into account.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1469-5) contains supplementary material, which is available to authorized users.     

A genomic scan for divergent selection in a secondary contact zone between Atlantic and Gulf of Mexico oysters, Crassostrea virginica

The degree of population structure within species often varies considerably among loci. This makes it difficult to determine whether observed variance reflects neutral-drift stochasticity or locus-specific selection at one or more loci. This uncertainty is exacerbated when evolutionary equilibrium cannot be assumed and/or admixture potentially inflates genomic variance. Thus, the value of a 'genome scan', where locus-specific summary statistics are compared with a simulated neutral distribution among loci, may be limited in secondary contact zones if the null distribution is sensitive to the timing of secondary contact. Of particular interest here is the wide variance previously observed in locus-specific patterns of population structure between Atlantic and Gulf of Mexico populations of eastern oyster, Crassostrea virginica. To test the robustness of an equilibrium null model, we compared among-locus distributions of F(ST) simulated under migration-drift equilibrium and several nonequilibrium secondary contact histories. We then tested for evidence of divergent selection between two oyster populations on either side of a secondary contact zone using 215 amplified fragment length polymorphism (AFLP) loci. Constant-migration equilibrium and nonequilibrium secondary-contact simulations produced equivalent distributions of F(ST) when anchored by the global mean F(ST) observed in oysters, 0.0917. The 99th quantile of simulated neutral F(ST) encompassed most of the variation among oyster loci. Three AFLP loci exhibited F(ST) values higher than this threshold. Although no locus was significant after correcting for multiple tests, our results show in geographically clinal organisms: AFLPs can efficiently characterize the genomic distribution of F(ST); equilibrium models can be used to evaluate outliers; these procedures help focus research on genomic regions of interest.     

Genetic structure of natural populations of the California black abalone (Haliotis cracherodii Leach, 1814), a candidate for endangered species status

The black abalone (Haliotis cracherodii) has been severely depleted across much of its historic range by a combination of overexploitation and disease. Natural recovery of extirpated populations along the southern California coast will depend on the extent to which remnant populations can serve as larval sources to geographic locations formerly supporting abalone populations. Population genetic analyses of mitochondrial cytochrome oxidase subunit one (COI) DNA sequences, four nuclear microsatellites, and 142 amplified fragment length polymorphisms (AFLPs) were used to evaluate connectivity among populations of H. cracherodii sampled from the central California coast and four islands in the Southern California Bight. Global divergence among populations was significant at COI and the AFLP loci. The Hka28 microsatellite locus and AFLP data showed significant divergence in multiple pairwise population comparisons and exhibited a signal of isolation by distance. Although estimates of gene flow based on genetic analyses must be interpreted with caution, the observed level of interpopulation genetic divergence suggests that larval dispersal is restricted, and natural recovery of decimated H. cracherodii populations along the coast of California is unlikely to occur in the near future.     

Divergent selection in trailing- versus leading-edge populations of Biscutella laevigata

Background and Aims

Knowledge on how climate-induced range shifts might affect natural selection is crucial to understand the evolution of species ranges.

Methods

Using historical demographic perspectives gathered from regional-scale phylogeography on the alpine herb Biscutella laevigata, indirect inferences on gene flow and signature of selection based on AFLP genotyping were compared between local populations persisting at the trailing edge and expanding at the leading edge.

Key Results

Spatial autocorrelation revealed that gene flow was two times more restricted at the trailing edge and genome scans indicated divergent selection in this persisting population. In contrast, no pattern of selection emerged in the expanding population at the leading edge.

Conclusions

Historical effects may determine different architecture of genetic variation and selective patterns within local populations, what is arguably important to understand evolutionary processes acting across the species ranges.     

AFLPMax: a user-friendly application for computing the optimal number of amplified fragment length polymorphism markers needed in phylogenetic reconstruction

Amplified fragment length polymorphisms (AFLPs) are widely used for phylogenetic inference especially in non-model species. Frequently, trees obtained with other nuclear or mitochondrial markers or with morphological information need additional resolution, increased branch support, or independent data sources (i.e. unlinked loci). In such cases, the use of AFLPs is a quick and cheap option. Computer simulation has shown that dominant AFLP markers lead to less accurate tree topologies than bi-allelic codominant markers such as SNPs, but this difference becomes negligible for shallow trees when using AFLP data sets that include a sufficiently large number of characters. Thus, determining how many AFLP characters are required to recover a given phylogeny is a key issue regarding the appropriateness of AFLPs for phylogenetic reconstruction. Here, we present a user-friendly, java-based graphical interface, AFLPMax, which executes an automatic pipeline of different programs providing the user with the optimal number of AFLP characters needed to recover a given phylogeny with high accuracy and support. Executables for Windows, linux and MacOS X operating systems, source code and user manual are available from: http://webs.uvigo.es/acraaj/AFLPMax.htm.     

The status of AFLP in the genomics era and a pipeline for converting AFLPs into single-locus markers

Even though next-generation sequencing (NGS) has now become the predominant state-of-the-art technique for genotyping populations, amplified fragment length polymorphism (AFLP) DNA fingerprinting is still a relevant method, thanks to its versatility, cost-effectiveness, independence of prior sequence information and broad applicability. Even though the number of AFLP studies reached its peak in 2012, it is still applied extensively for phylogenetic analysis, genotyping or identifying non-model species, which often feature complex and large genomes. For these purposes, tools continue to be developed for designing AFLP studies, scoring AFLPs or handling AFLP data. Moreover, AFLP studies embrace the NGS technology; for example, the whole-genome sequence of model species is used to design more efficient AFLP studies for non-model species. Conversely, in complexity reduction of polymorphic sequences and restriction site-associated DNA sequencing studies, polymorphisms are often found to be present in many restriction sites, which can still be studied as AFLPs. We discuss the latest advances in AFLP-based studies in the era of NGS and anticipate that AFLP will remain a relevant method in the near future, even for species with a known genome, owing to its many promising new features such as methylation-sensitive-AFLP. Here, we also present an optimized pipeline for converting AFLP markers into single-locus markers, which can be applied in both traditional AFLP and NGS studies.     

Comparative analysis of AFLPs and SSRs efficiency in resolving population genetic structure of Mediterranean Solea vulgaris

The performance of different molecular markers in the assessment of population structure was tested using samples of Solea vulgaris collected in the Mediterranean within and outside the hypothetical dispersal ability of the species. A total of 172 individuals belonging to four population samples were analysed using 15 microsatellites [simple sequence repeats (SSRs)] and 153 amplified fragment length polymorphisms (AFLPs). Considering the global qualitative patterns, we found a correlation between SSRs and AFLPs in detecting genetic differentiation among samples. However, on a small geographical scale, AFLPs were able to discriminate individuals from neighbouring populations whereas SSRs were not, and the percentage of individuals correctly assigned to their population of origin was higher with AFLPs than with SSRs. The high number of loci analysed with the AFLP technique could increase the probability to include outlier loci in the analysis; however, the neutrality test performed on our data set did not show evidence of selection acting on the S. vulgaris samples. Even if the choice of the molecular marker depends mainly on the biological question to be addressed, the higher power of discrimination and the comparative technical ease of obtaining data from AFLPs with respect to SSRs suggest the use of AFLPs for many population genetics studies.     

Challenges and pitfalls in the characterization of anonymous outlier AFLP markers in non-model species: lessons from an ocellated lizard genome scan

In the last few years, dozens of studies have documented the detection of loci influenced by selection from genome scans in a wide range of non-model species. Many of those studies used amplified fragment length polymorphism (AFLP) markers, which became popular for being easily applicable to any organism. However, because they are anonymous markers, AFLPs impose many challenges for their isolation and identification. Most recent AFLP genome scans used capillary electrophoresis (CE), which adds even more obstacles to the isolation of bands with a specific size for sequencing. These caveats might explain the extremely low number of studies that moved from the detection of outlier AFLP markers to their actual isolation and characterization. We document our efforts to characterize a set of outlier AFLP markers from a previous genome scan with CE in ocellated lizards (Lacerta lepida). Seven outliers were successfully isolated, cloned and sequenced. Their sequences are noncoding and show internal indels or polymorphic repetitive elements (microsatellites). Three outliers were converted into codominant markers by using specific internal primers to sequence and screen population variability from undigested DNA. Amplification in closely related lizard species was also achieved, revealing remarkable interspecific conservation in outlier loci sequences. We stress the importance of following up AFLP genome scans to validate selection signatures of outlier loci, but also report the main challenges and pitfalls that may be faced during the process.     

Colour pattern specification in the Mocker swallowtail Papilio dardanus: the transcription factor invected is a candidate for the mimicry locus H

The swallowtail butterfly, Papilio dardanus, is an iconic example of a polymorphic Batesian mimic. The expression of various female-limited colour forms is thought to be controlled by a single autosomal locus, termed H, whose function in determining the wing pattern remains elusive. As a step towards the physical mapping of H, we established a set of 272 polymorphic amplified fragment length polymorphism (AFLP) markers (EcoRI-MseI). Segregation patterns in a 'female-informative' brood (exploiting the absence of crossing over in female Lepidoptera) mapped these AFLPs to 30 linkage groups (putative chromosomes). The difference between the hippocoon and cenea female forms segregating in this family resides on a single one of these linkage groups, defined by 14 AFLPs. In a 'male-informative' cross (markers segregating within a linkage group), a pair of AFLPs co-segregated closely with the two female forms, except in four recombinants out of 19 female offspring. Linkage with these AFLP markers using four further female-informative families demonstrated that the genetic factor determining other morphs (poultoni, lamborni and trimeni) also maps to this same linkage group. The candidate gene invected, obtained in a screen for co-segregation of developmental genes with the colour forms, resides in a 13.9 cM interval flanked by the two AFLP markers. In the male-informative family invected co-segregated perfectly with the hippocoon/cenea factor, despite the four crossovers with the AFLPs. These findings make invected, and possibly its closely linked paralogue engrailed, strong candidates for H. This is supported by their known role in eyespot specification in nymphalid butterfly wings.     

The chromosomal location of rye AFLP bands

23 AFLP bands were assigned to different rye chromosomes by means of two different sets of wheat-rye addition lines. Only one AFLP band could be assigned to 4R, and no specific AFLPs were found on the 5R chromosome. Only one AFLP band was explicitly assigned to 4R, and no specific AFLPs were found on the 5R chromosome. At least seven co-migrating AFLPs showed the same chromosomal location in both sets of addition lines. A total of 22 AFLPs were assigned to chromosome 1R using wheat-rye substitution lines. Six of them have counterparts in one of the addition lines analyzed, but only four have the same chromosomal location. Six and four of the total AFLPs located using addition (23) and substitution (22) lines segregated in the mapping population DS2 x RXL10, but only six were simultaneously assigned to the same chromosome by both approaches. Although co-migrating AFLPs could be located on different rye chromosomes using addition and substitution lines, we believe that AFLPs can be useful as rye chromosome markers.     

Comparison of genetic diversity estimates within and among populations of maritime pine using chloroplast simple-sequence repeat and amplified fragment length polymorphism data

We compared the genetic variation of Pinus pinaster populations using amplified fragment length polymorphism (AFLP) and chloroplast simple-sequence repeat (cpSSR) loci. Populations' levels of diversity within groups were found to be similar with AFLPs, but not with cpSSRs. The high interlocus variance associated with the AFLP loci could account for the lack of differences in the former. Although AFLPs revealed much lower genetic diversity than cpSSRs, the levels of among-population differentiation found with the two types of marker were similar, provided that loci showing fewer than four null-homozygotes, in any population, were pruned from the AFLP data. Moreover, the French and Portuguese populations were clearly differentiated from each other, with both markers. The Mantel test showed that the genetic distance matrix calculated using the AFLP data was correlated with the matrix derived from the cpSSRs. Because of the concordance found between markers we conclude that gene flow was indeed the predominant force shaping nuclear and chloroplastic genetic variation of the populations within regions, at the geographical scale studied.