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.

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.

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.

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.

Piezotolerant small-colony variants with increased thermotolerance, antibiotic susceptibility, and low invasiveness in a clonal Staphylococcus aureus population

Following a pressure treatment of a clonal Staphylococcus aureus culture with 400 MPa for 30 min, piezotolerant variants were isolated. Among 21 randomly selected survivors, 9 were piezotolerant and all formed small colonies on several agar media. The majority of the isolates showed increased thermotolerance, impaired growth, and reduced antibiotic resistance compared to the wild type. However, several nonpiezotolerant isolates also demonstrated impaired growth and the small-colony phenotype. In agglutination tests for the detection of protein A and fibrinogen, the piezotolerant variants showed weaker agglutination reactions than the wild type and the other isolates. All variants also showed defective production of the typical S. aureus golden color, a characteristic which has previously been linked with virulence. They were also less able to invade intestinal epithelial cells than the wild type. These S. aureus variants showed phenotypic similarities to previously isolated Listeria monocytogenes piezotolerant mutants that contained mutations in ctsR. Because of these similarities, possible alterations in the ctsR hypermutable regions of the S. aureus variants were investigated through amplified fragment length polymorphism analysis. No mutations were identified, and subsequently we sequenced the ctsR and hrcA genes of three representative variants, finding no mutations. This work demonstrates that S. aureus probably possesses a strategy resulting in an abundance of multiple-stressresistant variants within clonal populations. This strategy, however, seems to involve genes and regulatory mechanisms different from those previously reported for L. monocytogenes. We are in the process of identifying these mechanisms.

Construction of a Festuca pratensis BAC library for map-based cloning in Festulolium substitution lines

Introgression in Festulolium is a potentially powerful tool to isolate genes for a large number of traits which differ between Festuca pratensis Huds. and Lolium perenne L. Not only are hybrids between the two species fertile, but the two genomes can be distinguished by genomic in situ hybridisation and a high frequency of recombination occurs between homoeologous chromosomes and chromosome segments. By a programme of introgression and a series of backcrosses, L. perenne lines have been produced which contain small F. pratensis substitutions. This material is a rich source of polymorphic markers targeted towards any trait carried on the F. pratensis substitution not observed in the L. perenne background. We describe here the construction of an F. pratensis BAC library, which establishes the basis of a map-based cloning strategy in L. perenne. The library contains 49,152 clones, with an average insert size of 112 kbp, providing coverage of 2.5 haploid genome equivalents. We have screened the library for eight amplified fragment length polymorphism (AFLP) derived markers known to be linked to an F. pratensis gene introgressed into L. perenne and conferring a staygreen phenotype as a consequence of a mutation in primary chlorophyll catabolism. While for four of the markers it was possible to identify bacterial artificial chromosome (BAC) clones, the other four AFLPs were too repetitive to enable reliable identification of locus-specific BACs. Moreover, when the four BACs were partially sequenced, no obvious coding regions could be identified. This contrasted to BACs identified using cDNA sequences, when multiple genes were identified on the same BAC.

Genetic diversity and specialisation of Eudarluca caricis on some graminaceous Puccinia species

Eudarluca caricis is a common hyperparasite of rusts. A total of 100 cultures were isolated from six Puccinia species or forms growing on 10 species of British grasses at two sites approximately 3 km apart. 82 isolates collected in 2005 were partially sequenced at the ITS locus, and amplified fragment length polymorphism profiles generated for 86 isolates from 2005 and 12 from 2007. Partial ITS sequences of most isolates grouped closely, in a clade with previously reported graminaceous Puccinia isolates and a number of Melampsora isolates. A second clade was very distinct and contained mostly isolates from P. poarum on Poa trivialis. All isolates had distinct AFLP haplotypes. The P. poarum isolates were very distinct from isolates collected from other rusts at the same site. Isolates from P. brachypodii f. sp. arrehenatheri growing on Arrhenatherum elatius in 2005 and 2007 at the same location were distinct (P < 0.001). Isolates from each rust or grass in one year and site were more similar than expected from overall variation between isolates (P<0.001). Isolates from P. coronata on different grasses clustered together (with isolates from P. brachypodii f. sp. poae-nemoralis), suggesting partial host rust specialisation in E. caricis.

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.

Detection and identification of species-specific bacteria associated with synanthropic mites

Internal bacterial communities of synanthropic mites Acarus siro, Dermatophagoides farinae, Lepidoglyphus destructor, and Tyrophagus putrescentiae (Acari: Astigmata) were analyzed by culturing and culture-independent approaches from specimens obtained from laboratory colonies. Homogenates of surface-sterilized mites were used for cultivation on non-selective agar and DNA extraction. Isolated bacteria were identified by sequencing of the 16S rRNA gene. PCR amplified 16S rRNA genes were analyzed by terminal restriction fragment length polymorphism analysis (T-RFLP) and cloning sequencing. Fluorescence in situ hybridization using universal bacterial probes was used for direct bacterial localization. T-RFLP analysis of 16S rRNA gene revealed distinct species-specific bacterial communities. The results were further confirmed by cloning and sequencing (284 clones). L. destructor and D. farinae showed more diverse communities then A. siro and T. putrescentiae. In the cultivated part of the community, the mean CFUs from four mite species ranged from 5.2 × 102 to 1.4 × 103 per mite. D. farinae had significantly higher CFUs than the other species. Bacteria were located in the digestive and reproductive tract, parenchymatical tissue, and in bacteriocytes. Among the clones, Bartonella-like bacteria occurring in A. siro and T. putresecentiae represented a distinct group related to Bartonellaceae and to Bartonella-like symbionts of ants. The clones of high similarity to Xenorhabdus cabanillasii were found in L. destructor and D. farinae, and one clone related to Photorhabdus temperata in A. siro. Members of Sphingobacteriales cloned from D. farinae and A. siro clustered with the sequences of “Candidatus Cardinium hertigii” and as a separate novel cluster.

Association of the PPARGC1A gene polymorphism with diabetic nephropathy in an Asian Indian population (CURES-41)

BACKGROUND: The aim of this study was to evaluate the association of polymorphisms of the peroxisome proliferator-activated receptor gamma (PPARG) gene and peroxisome proliferators-activated receptor gamma co-activator 1 alpha (PPARGC1A) gene with diabetic nephropathy (DN) in Asian Indians. METHODS: Six common polymorphisms, 3 of the PPARG gene [-1279G/A, Pro12Ala, and His478His (C/T)] and 3 of the PPARGC1A gene (Thr394Thr, Gly482Ser, and +A2962G) were studied in 571 normal glucose-tolerant (NGT) subjects, 255 type 2 diabetic (T2D) subjects without nephropathy, and 141 DN subjects. Genotypes were determined by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and direct sequencing. Logistic regression analysis was performed to assess the covariables associated with DN. RESULTS: Among the 6 polymorphisms examined, only the Gly482Ser of the PPARGC1A gene was significantly associated with DN. The genotype frequency of Ser/Ser genotype of the PPARGC1A gene was 8.8% (50/571) in NGT subjects, 7.8% (20/255) in T2D subjects, and 29.8% (42/141) in DN subjects. The odds ratios (ORs) for DN for the susceptible Gly/Ser and Ser/Ser genotype after adjusting for age, sex, body mass index, and duration of diabetes were 2.14 [95% confidence interval (CI), 1.23-3.72; P = 0.007] and 8.01 (95% CI, 3.89-16.47; P < 0.001), respectively. The unadjusted OR for DN for the XA genotype of the Thr394Thr polymorphism was 1.87 (95% CI, 1.20-2.92; P = 0.006) compared to T2D subjects. However, the significance was lost (P = 0.061) when adjusted for age, sex, BMI, and duration of diabetes. The +A2962G of PPARGC1A and the 3 polymorphisms of PPARG were not associated with DN. CONCLUSION: The Gly482Ser polymorphism of the PPARGC1A gene is associated with DN in Asian Indians.

Lack of association between serum adiponectin levels and the Pro12Ala polymorphism in Asian Indians

AIMS: The aim of the study was to investigate the association of serum adiponectin levels with the Pro12Ala polymorphism of the peroxisome proliferator activated receptor-gamma (PPARG) gene in Asian Indians. METHODS: We selected 400 diabetic subjects, 200 with the Pro12Pro genotype (100 male and 100 female) and 200 with the Pro12Ala genotype (100 male and 100 female) and 400 age- and sex-matched normal glucose tolerance subjects with similar genotype profiles from the Chennai Urban Rural Epidemiology Study. Fasting serum adiponection levels were measured using radioimmunoassay. The Pro12Ala polymorphism was genotyped by PCR-restriction fragment length polymorphism using BstUI. RESULTS: All clinical and biochemical parameters were similar in the subjects with the Pro12Pro and Pro12Ala genotypes. There was no significant difference in serum adiponectin values between subjects with the Pro12Pro and Pro12Ala genotypes (males 5.4 vs. 5.8 microg/ml, P = 0.546; females 6.9 vs. 7.2 microg/ml, P = 0.748). Adiponectin values did not differ among these two genotypes even when categorized based on their diabetes status (normal glucose tolerance Pro12Pro 7.9 vs. Pro12Ala 7.7 microg/ml, P = 0.994; diabetes Pro12Pro 4.7 vs. Pro12Ala 5.4 microg/ml, P = 0.622). CONCLUSION: The Pro12Ala polymorphism of the PPARG gene is not associated with serum adiponectin levels in Asian Indians.