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.

Restriction fragment length polymorphism of polymerase chain reaction products applied to the differentiation of poultry campylobacters for epidemiological investigations

A technique for subtyping Camplobacter jejuni isolates has been developed by using the restriction fragment length polymorphism (Rnp) of polymerase chain reaction (PCR) products of the fluA and flaB genes. The technique was validated by using strains representing 28 serotypes of C jejuni and it may also be applied to C coli. From these strains 12 distinct RFLP profiles were observed but there was no direct relationship between the RFLP profile and the serotype. One hundred and thirty-five campylobacter isolates from 15 geographically distinct broiler flocks were investigated. All the isolates could be subtyped by using the RFLP method. Isolates from most of the flocks had a single RFLP profile despite data indicating that several serotypes were involved. Although it is possible that further restriction analysis may have demonstrated profile variations in these strains, it is more likely that antigenic variation can occur within genotypically related campylobacters. As a result, serotyping may give conflicting information for veterinary epidemiological purposes. This RFLP typing scheme appears to provide a suitable tool for the investigation of the sources and routes of transmission of campylobacters in chickens.

Suitability of cryopreservation for the long-term storage of rare and endangered plant species: a case history for Cosmos atrosanguineus

The suitability of cryopreservation for the secure, long-term storage of the rare and endangered species Cosmos atrosanguineus was investigated. Using encapsulation/dehydration of shoot tips in alginate strips, survival rates of up to 100 % and shoot regeneration of up to 35 % were achieved. Light and electron microscopy studies indicated that cellular damage to some regions of the shoot tip during the freeze/thaw procedure was high, although cell survival in and around the meristematic region allowed shoot tip regeneration. The genetic fingerprinting technique, amplified fragment length polymorphisms (AFLPs), showed that no detectable genetic variation was present between material of C. atrosanguineus at the time of initiation into tissue culture and that which had been cryopreserved, stored in liquid nitrogen for 12 months and regenerated. Wearied plantlets that were grown under glasshouse conditions exhibited no morphological variation from non-frozen controls. (C) 2003 Annals of Botany Company.

Friends or relatives? Phylogenetics and species delimitation in the controversial European orchid Genus ophrys

Background and Aims Highly variable, yet possibly convergent, morphology and lack of sequence variation have severely hindered production of a robust phylogenetic framework for the genus Ophrys. The aim of this study is to produce this framework as a basis for more rigorous species delimitation and conservation recommendations. Methods Nuclear and plastid DNA sequencing and amplified fragment length polymorphism (AFLP) were performed on 85 accessions of Ophrys, spanning the full range of species aggregates currently recognized. Data were analysed using a combination of parsimony and Bayesian tree-building techniques and by principal coordinates analysis. Key Results Complementary phylogenetic analyses and ordinations using nuclear, plastid and AFLP datasets identify ten genetically distinct groups (six robust) within the genus that may in turn be grouped into three sections (treated as subgenera by some authors). Additionally, genetic evidence is provided for a close relationship between the O. tenthredinifera, O. bombyliflora and O. speculum groups. The combination of these analytical techniques provides new insights into Ophrys systematics, notably recognition of the novel O. umbilicata group. Conclusions Heterogeneous copies of the nuclear ITS region show that some putative Ophrys species arose through hybridization rather than divergent speciation. The supposedly highly specific pseudocopulatory pollination syndrome of Ophrys is demonstrably 'leaky', suggesting that the genus has been substantially over-divided at the species level.

Intercontinental dispersal prior to human translocation revealed in a cryptogenic invasive tree

In this study, complementary species-level and intraspecific phylogenies were used to better circumscribe the original native range and history of translocation of the invasive tree Parkinsonia aculeata. Species-level phylogenies were reconstructed using three chloroplast gene regions, and amplified fragment length polymorphism (AFLP) markers were used to reconstruct the intraspecific phylogeny. Together, these phylogenies revealed the timescale of transcontinental lineage divergence and the likely source of recent introductions of the invasive. The sequence data showed that divergence between North American and Argentinean P. aculeata occurred at least 5.7 million years ago, refuting previous hypotheses of recent dispersal between North and South America. AFLP phylogenies revealed the most likely sources of naturalized populations. The AFLP data also identified putatively introgressed plants, underlining the importance of wide sampling of AFLPs and of comparison with uniparentally inherited marker data when investigating hybridizing groups. Although P. aculeata has generally been considered North American, these data show that the original native range of P. aculeata included South America; recent introductions to Africa and Australia are most likely to have occurred from South American populations.

Genomic analysis of recombinant Sabin clinical isolates

Recombination in Poliovirus vaccine strains is a very frequent phenomenon. In this report 23 polio/Sabin strains isolated from healthy vaccinees or from VAPP patients after OPV administration, were investigated in order to identify recombination sites from 2C to 3D regions of the poliovirus genome. RT-PCR, followed by Restriction Fragment Length Polymorphism (RFLP) screening analysis were applied in four distant genomic regions (5' UTR, VP1, 2C and 3C-3D) in order to detect any putative recombinant. The detected recombinants were sequenced from 2C to the end of the genome (3' UTR) and the exact recombination sites were determined with computational analysis. Five of the 23 isolated strains were recombinant in one genomic region, two of them in 2C, isolates EP16:S3/S2, EP23:S3/S1, two in 3D isolates EP6:S2/S1, EP12:S2/S1 and one in 3A isolate EP9:S2/Sl. Point mutations were found in strains EP3, EP6, EP9 and EP12. Recombination specific types and sites re-occurrence along with point mutations are discussed concerning the polioviruses evolution.

Ty-1 copia retrotransposon-based SSAP marker development in Cashew

The most popular retrotransposon-based molecular marker system in use at the present time is the sequence-specific amplification polymorphism (SSAP) system . This system exploits the insertional polymorphism of long terminal repeat (LTR) retrotransposons around the genome. Because the LTR sequence is used to design primers for this method, its successful application requires sequence information from the terminal region of the mobile elements . In this study, two LTR sequences were isolated from the cashew genome and used successfully to develop SSAP marker systems. These were shown to have higher levels of polymorphism than amplified fragment length polymorphic markers for this species.

Comparison of the euryarchaeal microbial community in guts and food-soil of the soil-feeding termite Cubitermes fungifaber across different soil types

Termites are an important component of tropical soil communities and have a significant affect on the structure and nutrient content of soil. Digestion in termites is related to gut structure, gut physico-chemical conditions and gut symbiotic microbiota. Here we describe the use of 16S rRNA gene sequencing and Terminal-restriction Fragment Length Polymorphism (T-RFLP) analysis to examine methanogenic Archaea (MA) in the guts and food-soil of the soil-feeder Cubitermes fungifaber Sjostedt across a range of soil types. If they are strictly vertically inherited, then MA in guts should be the same in all individuals even if the soils differ across sites. In contrast, gut MA should reflect what is present in soil if populations are merely a reflection of what is ingested as the insects forage. We show clear differences between the euryarchaeal communities in termite guts and in food-soils from five different sites. Analysis of 16S rRNA gene clones indicated little overlap between the gut and soil communities. Gut clones were related to a termite-derived Methanomicrobiales cluster, to Methanobrevibacter and, surprisingly, to the haloalkaliphile Natronococcus. Soil clones clustered with Methanosarcina, Methanomicrococcus or Rice Cluster I. T-RFLP analysis indicated that the archaeal communities in the soil samples differed from site to site, whereas those in termite guts were similar between sites. There was some overlap between the gut and soil communities but these may represent transient populations in either guts or soil. Our data does not support the hypothesis that termite gut MA are derived from their food soil but also does not support a purely vertical transmission of gut microflora.

Comparison of Euryarchaea strains in the guts and food-soil of the soil-feeding termite Cubitermes fungifaber across different soil types

Termites are an important component of tropical soil communities and have a significant effect on the structure and nutrient content of soil. Digestion in termites is related to gut structure, gut physicochemical conditions, and gut symbiotic microbiota. Here we describe the use of 16S rRNA gene sequencing and terminal-restriction fragment length polymorphism (T-RFLP) analysis to examine methanogenic archaea (MA) in the guts and food-soil of the soil-feeder Cubitermes fungifaber Sjostedt across a range of soil types. If these MA are strictly vertically inherited, then the MA in guts should be the same in all individuals even if the soils differ across sites. In contrast, gut MA should reflect what is present in soil if populations are merely a reflection of what is ingested as the insects forage. We show clear differences between the euryarchaeal communities in termite guts and in food-soils from five different sites. Analysis of 16S rRNA gene clones indicated little overlap between the gut and soil communities. Gut clones were related to a termite-derived Methanomicrobiales cluster, to Methanobrevibacter and, surprisingly, to the haloalkaliphile Natronococcus. Soil clones clustered with Methanosarcina, Methanomicrococcus, or rice cluster I. T-RFLP analysis indicated that the archaeal communities in the soil samples differed from site to site, whereas those in termite guts were similar between sites. There was some overlap between the gut and soil communities, but these may represent transient populations in either guts or soil. Our data do not support the hypothesis that termite gut MA are derived from their food-soil but also do not support a purely vertical transmission of gut microflora.

Reversible phase variation in the phnE gene, which is required for phosphonate metabolism in Escherichia coli K-12

It is known that Escherichia coli K-12 is cryptic (Phn(-)) for utilization of methyl phosphonate (MePn) and that Phn(+) variants can be selected for growth on MePn as the sole P source. Variants arise from deletion via a possible slip strand mechanism of one of three direct 8-bp repeat sequences in phnE, which restores function to a component of a putative ABC type transporter. Here we show that Phn(+) variants are present at the surprisingly high frequency of >10(-2) in K-12 strains. Amplified-fragment length polymorphism analysis was used to monitor instability in phnE in various strains growing under different conditions. This revealed that, once selection for growth on MePn is removed, Phn(+) revertants reappear and accumulate at high levels through reinsertion of the 8-bp repeat element sequence. It appears that, in K-12, phnE contains a high-frequency reversible gene switch, producing phase variation which either allows ("on" form) or blocks ("off" form) MePn utilization. The switch can also block usage of other metabolizable alkyl phosphonates, including the naturally occurring 2-aminoethylphosphonate. All K-12 strains, obtained from collections, appear in the "off" form even when bearing mutations in mutS, mutD, or dnaQ which are known to enhance slip strand events between repetitive sequences. The ability to inactivate the phnE gene appears to be unique to K-12 strains since the B strain is naturally Phn(+) and lacks the inactivating 8-bp insertion in phnE, as do important pathogenic strains for which genome sequences are known and also strains isolated recently from environmental sources.

Mapping quantitative trait loci for flag leave senescence as a yield determinant in winter wheat under optimal and drought stressed environments

The timing of flag leaf senescence (FLS) is an important determinant of yield under stress and optimal environments. A doubled haploid population derived from crossing the photo period-sensitive variety Beaver,with the photo period-insensitive variety Soissons, varied significantly for this trait, measured as the percent green flag leaf area remaining at 14 days and 35 days after anthesis. This trait also showed a significantly positive correlation with yield under variable environmental regimes. QTL analysis based on a genetic map derived from 48 doubled haploid lines using amplified fragment length polymorphism (AFLP) and simple sequence repeat (SSR) markers, revealed the genetic control of this trait. The coincidence of QTL for senescence on chromosomes 2B and 2D under drought-stressed and optimal environments, respectively, indicate a complex genetic mechanism of this trait involving the re-mobilisation of resources from the source to the sink during senescence.

DNA variation in a 13-Mb region including the F9 gene: inferring the genealogical history and causal role of a hemophilia B mutation (IVS 5+13 A -> G)

About 5.5% of all UK hemophilia B patients have the base substitution IVS 5+13 A-->G as the only change in their factor (F)IX gene (F9). This generates a novel donor splice site which fits the consensus better than the normal intron 5 donor splice. Use of the novel splice site should result in a missense mutation followed by the abnormal addition of four amino acids to the patients' FIX. In order to explain the prevalence of this mutation, its genealogical history is examined. Analysis of restriction fragment length polymorphism in the 21 reference UK individuals (from different families) with the above mutation showed identical haplotypes in 19 while two differed from the rest and from each other. In order to investigate the history of the mutation and to verify that it had occurred independently more than once, the sequence variation in 1.5-kb segments scattered over a 13-Mb region including F9 was examined in 18 patients and 15 controls. This variation was then analyzed with a recently developed Bayesian approach that reconstructs the genealogy of the gene investigated while providing evidence of independent mutations that contribute disconnected branches to the genealogical tree. The method also provides minimum estimates of the age of the mutation inherited by the members of coherent trees. This revealed that 17 or 18 mutant genes descend from a founder who probably lived 450 years ago, while one patient carries an independent mutation. The independent recurrence of the IVS5+13 A-->G mutation strongly supports the conclusion that it is the cause of these patients' mild hemophilia.

Suitability of cryopreservation for the long-term storage of rare and endangered plant species: a case history for cosmos atrosanguineus

The suitability of cryopreservation for the secure, long-term storage of the rare and endangered species Cosmos atrosanguineus was investigated. Using encapsulation/dehydration of shoot tips in alginate strips, survival rates of up to 100 % and shoot regeneration of up to 35 % were achieved. Light and electron microscopy studies indicated that cellular damage to some regions of the shoot tip during the freeze/thaw procedure was high, although cell survival in and around the meristematic region allowed shoot tip regeneration. The genetic fingerprinting technique, amplified fragment length polymorphisms (AFLPs), showed that no detectable genetic variation was present between material of C. atrosanguineus at the time of initiation into tissue culture and that which had been cryopreserved, stored in liquid nitrogen for 12 months and regenerated. Weaned plantlets that were grown under glasshouse conditions exhibited no morphological variation from non-frozen controls.

Adaptive divergence and speciation among sexual and pseudoviviparous populations of Festuca

Pseudovivipary is an environmentally induced flowering abnormality in which vegetative shoots replace seminiferous (sexual) inflorescences. Pseudovivipary is usually retained in transplantation experiments, indicating that the trait is not solely induced by the growing environment. Pseudovivipary is the defining characteristic of Festuca vivipara, and arguably the only feature separating this species from its closest seminiferous relative, Festuca ovina. We performed phylogenetic and population genetic analysis on sympatric F. ovina and F. vivipara samples to establish whether pseudovivipary is an adaptive trait that accurately defines the separation of genetically distinct Festuca species. Chloroplast and nuclear marker-based analyses revealed that variation at a geographical level can exceed that between F. vivipara and F. ovina. We deduced that F. vivipara is a recent species that frequently arises independently within F. ovina populations and has not accumulated significant genetic differentiation from its progenitor. We inferred local gene flow between the species. We identified one amplified fragment length polymorphism marker that may be linked to a pseudovivipary-related region of the genome, and several other markers provide evidence of regional local adaptation in Festuca populations. We conclude that F. vivipara can only be appropriately recognized as a morphologically and ecologically distinct species; it lacks genetic differentiation from its relatives. This is the first report of a ‘failure in normal flowering development’ that repeatedly appears to be adaptive, such that the trait responsible for species recognition constantly reappears on a local basis.