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.

Genome-wide association mapping to candidate polymorphism resolution in the unsequenced barley genome

Although commonplace in human disease genetics, genome-wide association (GWA) studies have only relatively recently been applied to plants. Using 32 phenotypes in the inbreeding crop barley, we report GWA mapping of 15 morphological traits across ∼500 cultivars genotyped with 1,536 SNPs. In contrast to the majority of human GWA studies, we observe high levels of linkage disequilibrium within and between chromosomes. Despite this, GWA analysis readily detected common alleles of high penetrance. To investigate the potential of combining GWA mapping with comparative analysis to resolve traits to candidate polymorphism level in unsequenced genomes, we fine-mapped a selected phenotype (anthocyanin pigmentation) within a 140-kb interval containing three genes. Of these, resequencing the putative anthocyanin pathway gene HvbHLH1 identified a deletion resulting in a premature stop codon upstream of the basic helix-loop-helix domain, which was diagnostic for lack of anthocyanin in our association and biparental mapping populations. The methodology described here is transferable to species with limited genomic resources, providing a paradigm for reducing the threshold of map-based cloning in unsequenced crops.

A novel association of a polymorphism in the first intron of adiponectin gene with type 2 diabetes, obesity and hypoadiponectinemia in Asian Indians

Adiponectin is an adipose tissue specific protein that is decreased in subjects with obesity and type 2 diabetes. The objective of the present study was to examine whether variants in the regulatory regions of the adiponectin gene contribute to type 2 diabetes in Asian Indians. The study comprised of 2,000 normal glucose tolerant (NGT) and 2,000 type 2 diabetic, unrelated subjects randomly selected from the Chennai Urban Rural Epidemiology Study (CURES), in southern India. Fasting serum adiponectin levels were measured by radioimmunoassay. We identified two proximal promoter SNPs (-11377C-->G and -11282T-->C), one intronic SNP (+10211T-->G) and one exonic SNP (+45T-->G) by SSCP and direct sequencing in a pilot study (n = 500). The +10211T-->G SNP alone was genotyped using PCR-RFLP in 4,000 study subjects. Logistic regression analysis revealed that subjects with TG genotype of +10211T-->G had significantly higher risk for diabetes compared to TT genotype [Odds ratio 1.28; 95% Confidence Interval (CI) 1.07-1.54; P = 0.008]. However, no association with diabetes was observed with GG genotype (P = 0.22). Stratification of the study subjects based on BMI showed that the odds ratio for obesity for the TG genotype was 1.53 (95%CI 1.3-1.8; P < 10(-7)) and that for GG genotype, 2.10 (95% CI 1.3-3.3; P = 0.002). Among NGT subjects, the mean serum adiponectin levels were significantly lower among the GG (P = 0.007) and TG (P = 0.001) genotypes compared to TT genotype. Among Asian Indians there is an association of +10211T-->G polymorphism in the first intron of the adiponectin gene with type 2 diabetes, obesity and hypoadiponectinemia.

A novel parallel approach to the likelihood-based estimation of admixture in population genetics

Inferring population admixture from genetic data and quantifying it is a difficult but crucial task in evolutionary and conservation biology. Unfortunately state-of-the-art probabilistic approaches are computationally demanding. Effectively exploiting the computational power of modern multiprocessor systems can thus have a positive impact to Monte Carlo-based simulation of admixture modeling. A novel parallel approach is briefly described and promising results on its message passing interface (MPI)-based C implementation are reported.