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.

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.

Role of genetic polymorphism peroxisome proliferator-activated receptor-gamma2 Pro12Ala on ethnic susceptibility to diabetes in South-Asian and Caucasian subjects: Evidence for heterogeneity

OBJECTIVE: To determine whether the peroxisome proliferator-activated receptor (PPAR)-gamma Pro12ala polymorphism modulates susceptibility to diabetes in South Asians. RESEARCH DESIGN AND METHODS: South Asians (n = 697) and Caucasians (n = 457) living in Dallas/Forth Worth, Texas, and South Asians living in Chennai, India (n = 1,619), were enrolled for this study. PPAR-gamma Pro12Ala was determined using restriction fragment-length polymorphism. Insulin responsiveness to an oral glucose tolerance test (OGTT) was measured in nondiabetic subjects. RESULTS: The Caucasian diabetic subjects had significantly lower prevalence of PPAR-gamma 12Ala when compared with the Caucasian nondiabetic subjects (20 vs. 9%, P = 0.006). However, there were no significant differences between diabetic and nondiabetic subjects with reference to the Pro12Ala polymorphism among the South Asians living in Dallas (20 vs. 23%) and in India (19 vs. 19.3%). Although Caucasians carrying PPAR-gamma Pro12Ala had lower plasma insulin levels at 2 h of OGTT than the wild-type (Pro/Pro) carriers (76 +/- 68 and 54 +/- 33 microU/ml, respectively, P = 0.01), no differences in either fasting or 2-h plasma insulin concentrations were found between South Asians carrying the PPAR-gamma Pro12Ala polymorphism and those with the wild-type genotype at either Chennai or Dallas. CONCLUSIONS: Although further replication studies are necessary to test the validity of the described genotype-phenotype relationship, our study supports the hypothesis that the PPAR-gamma Pro12Ala polymorphism is protective against diabetes in Caucasians but not in South Asians.

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.

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.

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.