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.

Characterization of para-Nitrophenol-degrading bacterial communities in river water by using functional markers and stable isotope probing

Microbial degradation is a major determinant of the fate of pollutants in the environment. para-Nitrophenol (PNP) is an EPA listed priority pollutant with a wide environmental distribution, but little is known about the microorganisms that degrade it in the environment. We studied the diversity of active PNP-degrading bacterial populations in river water using a novel functional marker approach coupled with [13C6]PNP stable isotope probing (SIP). Culturing together with culture-independent terminal restriction fragment length polymorphism analysis of 16S rRNA gene amplicons identified Pseudomonas syringae to be the major driver of PNP degradation in river water microcosms. This was confirmed by SIP-pyrosequencing of amplified 16S rRNA. Similarly, functional gene analysis showed that degradation followed the Gram-negative bacterial pathway and involved pnpA from Pseudomonas spp. However, analysis of maleylacetate reductase (encoded by mar), an enzyme common to late stages of both Gram-negative and Gram-positive bacterial PNP degradation pathways, identified a diverse assemblage of bacteria associated with PNP degradation, suggesting that mar has limited use as a specific marker of PNP biodegradation. Both the pnpA and mar genes were detected in a PNP-degrading isolate, P. syringae AKHD2, which was isolated from river water. Our results suggest that PNP-degrading cultures of Pseudomonas spp. are representative of environmental PNP-degrading populations.

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.

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.

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.

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.

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.

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.

Recent developments in genetic data analysis: what can they tell us about human demographic history?

Over the last decade, a number of new methods of population genetic analysis based on likelihood have been introduced. This review describes and explains the general statistical techniques that have recently been used, and discusses the underlying population genetic models. Experimental papers that use these methods to infer human demographic and phylogeographic history are reviewed. It appears that the use of likelihood has hitherto had little impact in the field of human population genetics, which is still primarily driven by more traditional approaches. However, with the current uncertainty about the effects of natural selection, population structure and ascertainment of single-nucleotide polymorphism markers, it is suggested that likelihood-based methods may have a greater impact in the future.