An approximate Bayesian computation approach to overcome biases that arise when using amplified fragment length polymorphism markers to study population structure

There is great interest in using amplified fragment length polymorphism (AFLP) markers because they are inexpensive and easy to produce. It is, therefore, possible to generate a large number of markers that have a wide coverage of species genotnes. Several statistical methods have been proposed to study the genetic structure using AFLP's but they assume Hardy-Weinberg equilibrium and do not estimate the inbreeding coefficient, F-IS. A Bayesian method has been proposed by Holsinger and colleagues that relaxes these simplifying assumptions but we have identified two sources of bias that can influence estimates based on these markers: (i) the use of a uniform prior on ancestral allele frequencies and (ii) the ascertainment bias of AFLP markers. We present a new Bayesian method that avoids these biases by using an implementation based on the approximate Bayesian computation (ABC) algorithm. This new method estimates population-specific F-IS and F-ST values and offers users the possibility of taking into account the criteria for selecting the markers that are used in the analyses. The software is available at our web site (http://www-leca.uif-grenoble.fi-/logiciels.htm). Finally, we provide advice on how to avoid the effects of ascertainment bias.

Amplified fragment length polymorphism genetic fingerprinting challenges the taxonomic status of the near-endemic species Chara curta Nolte ex Kutz. (Characeae)

Amplified fragment length polymorphism (AFLP) genetic fingerprinting of 14 accessions of Chara curta and Chara aspera Willd., sampled across a range of habitats and morphologies in Britain, suggests that these taxa are part of the variation within a single species complex. Two primer combinations generating 397 fragments (97% of which were polymorphic), analysed by Jaccard's similarity coefficient and principal co-ordinate analysis, did not recover groups which reflect the current taxonomy. By contrast with the genetic study, a Gower general similarity coefficient and principal co-ordinate analysis of 52 morphological characters recovered the currently recognized species groups. A Mantel test showed no significant correlation between the genetic data and the morphological data, supporting the hypothesis that phenotypic variability in Chara L. is either to some extent environmentally induced or represents developmental stages. Implications for the conservation status of C. curta in Britain are discussed. (c) 2007 The Linnean Society of London, Botanical Journal of the Linnean Society, 2007, 155, 467-476.

Genotyping of human and porcine Yersinia enterocolitica, Yersinia intermedia, and Yersinia bercovieri strains from Switzerland by amplified fragment length polymorphism analysis

In this study, 231 strains of Yersinia enterocolitica, 25 strains of Y. intermedia, and 10 strains of Y. bercovieri from human and porcine sources (including reference strains) were analyzed using amplified fragment length polymorphism (AFLP), a whole-genome fingerprinting method for subtyping bacterial isolates. AFLP typing distinguished the different Yersinia species examined. Representatives of Y. enterocolitica biotypes 1A, 1B, 2, 3, and 4 belonged to biotype-related AFLP clusters and were clearly distinguished from each other. Y. enterocolitica biotypes 2, 3, and 4 appeared to be more closely related to each other (83% similarity) than to biotypes 1A (11%) and 1B (47%). Biotype 1A strains exhibited the greatest genetic heterogeneity of the biotypes studied. The biotype 1A genotypes were distributed among four major clusters, each containing strains from both human and porcine sources, confirming the zoonotic potential of this organism. The AFLP technique is a valuable genotypic method for identification and typing of Y. enterocolitica and other Yersinia spp.

Genotyping of Francisella tularensis strains by pulsed-field gel electrophoresis, amplified fragment length polymorphism fingerprinting, and 16S rRNA gene sequencing

We evaluated three molecular methods for identification of Francisella strains: pulsed-field gel electrophoresis (PFGE), amplified fragment length polymorphism (AFLP) analysis, and 16S rRNA gene sequencing. The analysis was performed with 54 Francisella tularensis subsp. holarctica, 5 F. tularensis subsp. tularensis, 2 F. tularensis subsp. novicida, and 1 F. philomiragia strains. On the basis of the combination of results obtained by PFGE with the restriction enzymes XhoI and BamHI, PFGE revealed seven pulsotypes, which allowed us to discriminate the strains to the subspecies level and which even allowed us to discriminate among some isolates of F. tularensis subsp. holarctica. The AFLP analysis technique produced some degree of discrimination among F. tularensis subsp. holarctica strains (one primary cluster with three major subclusters and minor variations within subclusters) when EcoRI-C and MseI-A, EcoRI-T and MseI-T, EcoRI-A and MseI-C, and EcoRI-0 and MseI-CA were used as primers. The degree of similarity among the strains was about 94%. The percent similarities of the AFLP profiles of this subspecies compared to those of F. tularensis subsp. tularensis, F. tularensis subsp. novicida, and F. philomiragia were less than 90%, about 72%, and less than 24%, respectively, thus permitting easy differentiation of this subspecies. 16S rRNA gene sequencing revealed 100% similarity for all F. tularensis subsp. holarctica isolates compared in this study. These results suggest that although limited genetic heterogeneity among F. tularensis subsp. holarctica isolates was observed, PFGE and AFLP analysis appear to be promising tools for the diagnosis of infections caused by different subspecies of F. tularensis and suitable techniques for the differentiation of individual strains.

Molecular characterization by amplified fragment length polymorphism and aflatoxin production of Aspergillus flavus isolated from freshly harvested peanut in Brazil

Brazil contributes substantially to the global peanut production, and the state of Sao Paulo is the largest producer in the country. Peanut crops can be contaminated by Aspergillus flavus strains producing aflatoxins, which are highly toxic and carcinogenic. Thus, the production of high-quality peanuts is crucial both for the commercial peanut industry and as a matter of public health. In this study, we used amplified fragment length polymorphism analysis (AFLP) to investigate the genetic variability among A. flavus strains isolated from fresh peanuts harvested in four different regions in the state of Sao Paulo, and to determine whether the molecular genetic profiles correlated with aflatoxin production or sclerotia formation. AFLP analysis generated 78 fragments ranging from 27 to 365 base pairs in length. Thirteen percent were not polymorphic. Genotyping identified twelve groups of A. flavus. On the basis of the polymorphisms identified, similarity between the isolates ranged from 37% to 100%. Of all isolates collected, 91.7% produced aflatoxins and 83.9% produced small sclerotia. Statistical analysis failed to suggest any relationship between the presence of sclerotia and mean levels of aflatoxins B-1 and B-2. Furthermore, a dendrogram based on AFLP data revealed substantial genetic variability among the A. flavus strains, but showed no correlation between dendrogram groups separated by molecular genetic features and production of aflatoxins B-1 or B-2 or the formation of sclerotia.

Genotyping of Brazilian Erysipelothrix spp. strains by amplified fragment length polymorphism

One hundred and fifty-one Erysipelothrix spp. isolates from diseased and carrier swine from Brazil were identified by PCR, submitted to serotyping and analyzed by amplified fragment length polymorphism with a single enzyme (AFLP). Reference strains from Australia and the United Kingdom were also examined. The 151 strains were classified into 18 different serotypes (1a, 1b, 2a, 2b, 4, 5, 6, 7, 8, 10, 11, 12, 15, 17, 19, 21, 24 and 25), being serotype 2b the most frequent (39.7%). By associating serotyping and PCR results, it was possible to identify 146 strains as E. rhusiopathiae and five strains as E. tonsillarum. Despite the fact that for this genus AFLP did not cluster all isolates according to serotype, origin, disease or isolation data, the execution of the technique was easy and fast, demonstrating high discriminatory power. The results produced by the AFLP analysis of Erysipelothrix spp. could also support its use as a discriminatory tool for E. rhusiopathiae and E. tonsillarum species. (C) 2010 Elsevier B.V. All rights reserved.

Genetic structure of populations of Pissodes castaneus (De Geer) (Coleoptera, Curculionidae) using amplified fragment length polymorphism

Genetic structure of populations of Pissodes castaneus (De Geer) (Coleoptera, Curculionidae) using amplified fragment length polymorphism. The objective of this study was to determine the genetic structure of populations of Pissodes castaneus from different areas and on different species of Pinus using the PCR-AFLP technique. Twenty samples were analyzed, representing 19 populations from Brazil and one from Florence, Italy, which is the region of origin of P. castaneus. The four combinations of primers generated a total of 367 fragments of DNA, and 100% of polymorphic loci, indicating high degree of molecular polymorphism. The dendrogram did not reveal trends for grouping the populations in relation to origin. The low genetic similarity (0.11 between the most distant groups) and genetic distances of 0.13 and 0.44 for 10 out of the 20 samples may indicate several founding events or multiple introductions of heterogeneous strains into Brazil. The allelic fixation index (Fst) was 0.3851, considered high, and the number of migrants (Nm) was 0.3991, indicating low gene flow among populations. The highest genetic distances were between the population from Irani, SC and Cambará do Sul, RS and Bituruna, PR, indicating an independent founding event or a particular allelic fixation in the former location. The high genetic diversity among populations points out that the populations are genetically heterogeneous with a diverse gene pool in the surveyed areas, what makes them to respond differently to control measures.

Fluorescent amplified fragment length polymorphism genotyping of human and animal Staphylococcus aureus isolates from dairy farms with manual milking

Fluorescence amplified fragment length polymorphism (fAFLP) was used to assess the genetic relatedness of 40 Staphylococcus aureus strains isolated from human and animal skin samples in seven dairy farms with manual milking. S. aureus was isolated from 11 out of 30 (36%) human skin samples and from 29 out of 100 (29%) teat skin samples from apparently healthy cows. Genomic DNA from each isolate was double-digested with EcoRI and MseI and complementary oligonucleotide adaptors were ligated to the restriction fragments. Pre-selective and selective, amplification reactions were performed, the amplified fragments were separated by electrophoresis in an ABI377 sequencer and analysed using GeneScan 3.1 and Genotyper 2.5. Three single isolates (a-c), a predominant cluster with 35 isolates (d) and another cluster with two isolates (e) were identified. Both clusters d and e included human and animal isolates genetically related, because the profiles had 90-100% homology. Since no cluster was comprised uniquely of human or animal isolates and given the close genetic relatedness among human and animal samples in the farms, the present findings support the. hypothesis that dairy workers can spread S. aureus through manual milking. (C) 2005 Elsevier B.V. All rights reserved.

Fluorescent amplified fragment length polymorphism (fAFLP) analyses and genetic diversity in Litopenaeus vannamei (Penaeidae)

The Pacific white shrimp, Litopenaeus vannamei (Penaeidae), represents about 95% of all Brazilian shrimp production. The Brazilian L. vannamei foundation broodstock was made up of specimens collected from different American Pacific sites, but little information was collected on the genetic structure of the broodstock. We used the fluorescence amplified fragment length polymorphism (fAFLP) method to study the genetic diversity of L. vannamei broodstock lines 03CMF1 and 03CBF1 originally produced by breeder-shrimps imported mainly from Panama and Ecuador, although wild individuals from other localities may also have been used in producing these two lines. Our results showed a total of 93 polymorphic bands ranging from 50 to 500 bp, the mean Nei's genetic diversity calculated for the total sample was 13.4% and identity and genetic distance analyses indicated high genetic homogeneity within and between both the broodstock lineages studied which suggests that they had similar genetic structure. These results may represent an important tool for the appropriate management of L. vannamei broodstocks. Copyright by the Brazilian Society of Genetics.

Combining niche modelling and landscape genetics to study local adaptation: A novel approach illustrated using alpine plants

Understanding the factors that shape adaptive genetic variation across species niches has become of paramount importance in evolutionary ecology, especially to understand how adaptation to changing climate affects the geographic range of species. The distribution of adaptive alleles in the ecological niche is determined by the emergence of novel mutations, their fitness consequences and gene flow that connects populations across species niches. Striking demographical differences and source sink dynamics of populations between the centre and the margin of the niche can play a major role in the emergence and spread of adaptive alleles. Although some theoretical predictions have long been proposed, the origin and distribution of adaptive alleles within species niches remain untested. In this paper, we propose and discuss a novel empirical approach that combines landscape genetics with species niche modelling, to test whether alleles that confer local adaptation are more likely to occur in either marginal or central populations of species niches. We illustrate this new approach by using a published data set of 21 alpine plant species genotyped with a total of 2483 amplified fragment length polymorphisms (AFLP), distributed over more than 1733 sampling sites across the Alps. Based on the assumption that alleles that were statistically associated with environmental variables were adaptive, we found that adaptive alleles in the margin of a species niche were also present in the niche centre, which suggests that adaptation originates in the niche centre. These findings corroborate models of species range evolution, in which the centre of the niche contributes to the emergence of novel adaptive alleles, which diffuse towards niche margins and facilitate niche and range expansion through subsequent local adaptation. Although these results need to be confirmed via fitness measurements in natural populations and functionally characterised genetic sequences, this study provides a first step towards understanding how adaptive genetic variation emerges and shapes species niches and geographic ranges along environmental gradients.

A second generation radiation hybrid map to aid the assembly of the bovine genome sequence

BACKGROUND: Several approaches can be used to determine the order of loci on chromosomes and hence develop maps of the genome. However, all mapping approaches are prone to errors either arising from technical deficiencies or lack of statistical support to distinguish between alternative orders of loci. The accuracy of the genome maps could be improved, in principle, if information from different sources was combined to produce integrated maps. The publicly available bovine genomic sequence assembly with 6x coverage (Btau_2.0) is based on whole genome shotgun sequence data and limited mapping data however, it is recognised that this assembly is a draft that contains errors. Correcting the sequence assembly requires extensive additional mapping information to improve the reliability of the ordering of sequence scaffolds on chromosomes. The radiation hybrid (RH) map described here has been contributed to the international sequencing project to aid this process. RESULTS: An RH map for the 30 bovine chromosomes is presented. The map was built using the Roslin 3000-rad RH panel (BovGen RH map) and contains 3966 markers including 2473 new loci in addition to 262 amplified fragment-length polymorphisms (AFLP) and 1231 markers previously published with the first generation RH map. Sequences of the mapped loci were aligned with published bovine genome maps to identify inconsistencies. In addition to differences in the order of loci, several cases were observed where the chromosomal assignment of loci differed between maps. All the chromosome maps were aligned with the current 6x bovine assembly (Btau_2.0) and 2898 loci were unambiguously located in the bovine sequence. The order of loci on the RH map for BTA 5, 7, 16, 22, 25 and 29 differed substantially from the assembled bovine sequence. From the 2898 loci unambiguously identified in the bovine sequence assembly, 131 mapped to different chromosomes in the BovGen RH map. CONCLUSION: Alignment of the BovGen RH map with other published RH and genetic maps showed higher consistency in marker order and chromosome assignment than with the current 6x sequence assembly. This suggests that the bovine sequence assembly could be significantly improved by incorporating additional independent mapping information.

Occurrence of Mycobacterium avium subspecies paratuberculosis across host species and European countries with evidence for transmission between wildlife and domestic ruminants

BACKGROUND Mycobacterium avium subspecies paratuberculosis (Map) causes an infectious chronic enteritis (paratuberculosis or Johne's disease) principally of ruminants. The epidemiology of Map is poorly understood, particularly with respect to the role of wildlife reservoirs and the controversial issue of zoonotic potential (Crohn's disease). Genotypic discrimination of Map isolates is pivotal to descriptive epidemiology and resolving these issues. This study was undertaken to determine the genetic diversity of Map, enhance our understanding of the host range and distribution and assess the potential for interspecies transmission. RESULTS 164 Map isolates from seven European countries representing 19 different host species were genotyped by standardized IS900--restriction fragment length polymorphism (IS900-RFLP), pulsed-field gel electrophoresis (PFGE), amplified fragment length polymorphisms (AFLP) and mycobacterial interspersed repeat unit-variable number tandem repeat (MIRU-VNTR) analyses. Six PstI and 17 BstEII IS900-RFLP, 31 multiplex [SnaBI-SpeI] PFGE profiles and 23 MIRU-VNTR profiles were detected. AFLP gave insufficient discrimination of isolates for meaningful genetic analysis. Point estimates for Simpson's index of diversity calculated for the individual typing techniques were in the range of 0.636 to 0.664 but a combination of all three methods increased the discriminating power to 0.879, sufficient for investigating transmission dynamics. Two predominant strain types were detected across Europe with all three typing techniques. Evidence for interspecies transmission between wildlife and domestic ruminants on the same property was demonstrated in four cases, between wildlife species on the same property in two cases and between different species of domestic livestock on one property. CONCLUSION The results of this study showed that it is necessary to use multiple genotyping techniques targeting different sources of genetic variation to obtain the level of discrimination necessary to investigate transmission dynamics and trace the source of Map infections. Furthermore, the combination of genotyping techniques may depend on the geographical location of the population to be tested. Identical genotypes were obtained from Map isolated from different host species co-habiting on the same property strongly suggesting that interspecies transmission occurs. Interspecies transmission of Map between wildlife species and domestic livestock on the same property provides further evidence to support a role for wildlife reservoirs of infection.

A population genetics study of Anopheles darlingi (Diptera: Culicidae) from Colombia based on random amplified polymorphic DNA-polymerase chain reaction and amplified fragment lenght polymorphism markers

The genetic variation and population structure of three populations of Anopheles darlingi from Colombia were studied using random amplified polymorphic markers (RAPDs) and amplified fragment length polymorphism markers (AFLPs). Six RAPD primers produced 46 polymorphic fragments, while two AFLP primer combinations produced 197 polymorphic fragments from 71 DNA samples. Both of the evaluated genetic markers showed the presence of gene flow, suggesting that Colombian An. darlingi populations are in panmixia. Average genetic diversity, estimated from observed heterozygosity, was 0.374 (RAPD) and 0.309 (AFLP). RAPD and AFLP markers showed little evidence of geographic separation between eastern and western populations; however, the F ST values showed high gene flow between the two western populations (RAPD: F ST = 0.029; Nm: 8.5; AFLP: F ST = 0.051; Nm: 4.7). According to molecular variance analysis (AMOVA), the genetic distance between populations was significant (RAPD:phiST = 0.084; AFLP:phiST = 0.229, P < 0.001). The F ST distances and AMOVAs using AFLP loci support the differentiation of the Guyana biogeographic province population from those of the Chocó-Magdalena. In this last region, Chocó and Córdoba populations showed the highest genetic flow.

An integrated molecular map of yellow passion fruit based on simultaneous maximum-likelihood estimation of linkage and linkage phases

The development of genetic maps for auto-incompatible species, such as the yellow passion fruit (Passiflora edulis Sims f.flavicarpa Deg.) is restricted due to the unfeasibility of obtaining traditional mapping populations based on inbred lines. For this reason, yellow passion fruit linkage maps were generally constructed using a strategy known as two-way pseudo-testeross, based on monoparental dominant markers segregating in a 1:1 fashion. Due to the lack of information from these markers in one of the parents, two individual (parental) maps were obtained. However, integration of these maps is essential, and biparental markers can be used for such an operation. The objective of our study was to construct an integrated molecular map for a full-sib population of yellow passion fruit combining different loci configuration generated from amplified fragment length polymorphisms (AFLPs) and microsatellite markers and using a novel approach based on simultaneous maximum-likelihood estimation of linkage and linkage phases, specially designed for outcrossing species. Of the total number of loci, approximate to 76%, 21%, 0.7%, and 2.3% did segregate in 1:1, 3:1, 1:2:1, and 1:1:1:1 ratios, respectively. Ten linkage groups (LGs) were established with a logarithm of the odds (LOD) score >= 5.0 assuming a recombination fraction : <= 0.35. On average, 24 markers were assigned per LG, representing a total map length of 1687 cM, with a marker density of 6.9 cM. No markers were placed as accessories on the map as was done with previously constructed individual maps.

Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation.

Identifying adaptive genetic variation is a challenging task, in particular in non-model species for which genomic information is still limited or absent. Here, we studied distribution patterns of amplified fragment length polymorphisms (AFLPs) in response to environmental variation, in 13 alpine plant species consistently sampled across the entire European Alps. Multiple linear regressions were performed between AFLP allele frequencies per site as dependent variables and two categories of independent variables, namely Moran's eigenvector map MEM variables (to account for spatial and unaccounted environmental variation, and historical demographic processes) and environmental variables. These associations allowed the identification of 153 loci of ecological relevance. Univariate regressions between allele frequency and each environmental factor further showed that loci of ecological relevance were mainly correlated with MEM variables. We found that precipitation and temperature were the best environmental predictors, whereas topographic factors were rarely involved in environmental associations. Climatic factors, subject to rapid variation as a result of the current global warming, are known to strongly influence the fate of alpine plants. Our study shows, for the first time for a large number of species, that the same environmental variables are drivers of plant adaptation at the scale of a whole biome, here the European Alps.