Single-nucleotide polymorphism, linkage disequilibrium and geographic structure in the malaria parasite Plasmodium vivax: prospects for genome-wide association studies

Background: The ideal malaria parasite populations for initial mapping of genomic regions contributing to phenotypes such as drug resistance and virulence, through genome-wide association studies, are those with high genetic diversity, allowing for numerous informative markers, and rare meiotic recombination, allowing for strong linkage disequilibrium (LD) between markers and phenotype-determining loci. However, levels of genetic diversity and LD in field populations of the major human malaria parasite P. vivax remain little characterized. Results: We examined single-nucleotide polymorphisms (SNPs) and LD patterns across a 100-kb chromosome segment of P. vivax in 238 field isolates from areas of low to moderate malaria endemicity in South America and Asia, where LD tends to be more extensive than in holoendemic populations, and in two monkey-adapted strains (Salvador-I, from El Salvador, and Belem, from Brazil). We found varying levels of SNP diversity and LD across populations, with the highest diversity and strongest LD in the area of lowest malaria transmission. We found several clusters of contiguous markers with rare meiotic recombination and characterized a relatively conserved haplotype structure among populations, suggesting the existence of recombination hotspots in the genome region analyzed. Both silent and nonsynonymous SNPs revealed substantial between-population differentiation, which accounted for similar to 40% of the overall genetic diversity observed. Although parasites clustered according to their continental origin, we found evidence for substructure within the Brazilian population of P. vivax. We also explored between-population differentiation patterns revealed by loci putatively affected by natural selection and found marked geographic variation in frequencies of nucleotide substitutions at the pvmdr-1 locus, putatively associated with drug resistance. Conclusion: These findings support the feasibility of genome-wide association studies in carefully selected populations of P. vivax, using relatively low densities of markers, but underscore the risk of false positives caused by population structure at both local and regional levels.

Mendelian inheritance, linkage and linkage disequilibrium in microsatellite loci of Copaifera langsdorffii Desf.

Copaifera langsdorffii is a Neotropical tree with wide distribution in the Brazilian Atlantic rain forest and savanna. Although eight microsatellite loci (SSR) were developed in 2000 and have been widely used since then, there is yet no information about their inheritance, linkage and linkage disequilibrium (LD). Through the analysis of 28 open-pollinated (OP) progenies, the SSR loci revealed Mendelian inheritance and independent assortment. Using these progenies, young and adult trees LD was mainly detected in OP progenies. Our results show clear evidence that the eight SSR loci can be used without restriction in genetic diversity, mating system and parentage analysis.

Pharmacogenetics of Warfarin: Development of a Dosing Algorithm for Brazilian Patients

A dosing algorithm including genetic (VKORC1 and CYP2C9 genotypes) and nongenetic factors (age, weight, therapeutic indication, and cotreatment with amiodarone or simvastatin) explained 51% of the variance in stable weekly warfarin doses in 390 patients attending an anticoagulant clinic in a Brazilian public hospital. The VKORC1 3673G>A genotype was the most important predictor of warfarin dose, with a partial R(2) value of 23.9%. Replacing the VKORC1 3673G>A genotype with VKORC1 diplotype did not increase the algorithm`s predictive power. We suggest that three other single-nucleotide polymorphisms (SNPs) (5808T>G, 6853G>C, and 9041G>A) that are in strong linkage disequilibrium (LD) with 3673G>A would be equally good predictors of the warfarin dose requirement. The algorithm`s predictive power was similar across the self-identified ""race/color"" subsets. ""Race/color"" was not associated with stable warfarin dose in the multiple regression model, although the required warfarin dose was significantly lower (P = 0.006) in white (29 +/- 13 mg/week, n = 196) than in black patients (35 +/- 15 mg/week, n = 76).

An empirical evaluation of imputation accuracy for association statistics reveals increased type-I error rates in genome-wide associations

Background: Genome wide association studies (GWAS) are becoming the approach of choice to identify genetic determinants of complex phenotypes and common diseases. The astonishing amount of generated data and the use of distinct genotyping platforms with variable genomic coverage are still analytical challenges. Imputation algorithms combine directly genotyped markers information with haplotypic structure for the population of interest for the inference of a badly genotyped or missing marker and are considered a near zero cost approach to allow the comparison and combination of data generated in different studies. Several reports stated that imputed markers have an overall acceptable accuracy but no published report has performed a pair wise comparison of imputed and empiric association statistics of a complete set of GWAS markers. Results: In this report we identified a total of 73 imputed markers that yielded a nominally statistically significant association at P < 10(-5) for type 2 Diabetes Mellitus and compared them with results obtained based on empirical allelic frequencies. Interestingly, despite their overall high correlation, association statistics based on imputed frequencies were discordant in 35 of the 73 (47%) associated markers, considerably inflating the type I error rate of imputed markers. We comprehensively tested several quality thresholds, the haplotypic structure underlying imputed markers and the use of flanking markers as predictors of inaccurate association statistics derived from imputed markers. Conclusions: Our results suggest that association statistics from imputed markers showing specific MAF (Minor Allele Frequencies) range, located in weak linkage disequilibrium blocks or strongly deviating from local patterns of association are prone to have inflated false positive association signals. The present study highlights the potential of imputation procedures and proposes simple procedures for selecting the best imputed markers for follow-up genotyping studies.

A Genome-Wide Association Study of Upper Aerodigestive Tract Cancers Conducted within the INHANCE Consortium

Genome-wide association studies (GWAS) have been successful in identifying common genetic variation involved in susceptibility to etiologically complex disease. We conducted a GWAS to identify common genetic variation involved in susceptibility to upper aero-digestive tract (UADT) cancers. Genome-wide genotyping was carried out using the Illumina HumanHap300 beadchips in 2,091 UADT cancer cases and 3,513 controls from two large European multi-centre UADT cancer studies, as well as 4,821 generic controls. The 19 top-ranked variants were investigated further in an additional 6,514 UADT cancer cases and 7,892 controls of European descent from an additional 13 UADT cancer studies participating in the INHANCE consortium. Five common variants presented evidence for significant association in the combined analysis (p <= 5 x 10(-7)). Two novel variants were identified, a 4q21 variant (rs1494961, p = 1 x 10(-8)) located near DNA repair related genes HEL308 and FAM175A (or Abraxas) and a 12q24 variant (rs4767364, p = 2 x 10(-8)) located in an extended linkage disequilibrium region that contains multiple genes including the aldehyde dehydrogenase 2 (ALDH2) gene. Three remaining variants are located in the ADH gene cluster and were identified previously in a candidate gene study involving some of these samples. The association between these three variants and UADT cancers was independently replicated in 5,092 UADT cancer cases and 6,794 controls non-overlapping samples presented here (rs1573496-ADH7, p = 5 x 10(-8); rs1229984-ADH1B, p = 7 x 10(-9); and rs698-ADH1C, p = 0.02). These results implicate two variants at 4q21 and 12q24 and further highlight three ADH variants in UADT cancer susceptibility.

Association of SNPs on CAPN1 and CAST genes with tenderness in Nellore cattle

We examined whether single-nucleotide polymorphisms (SNPs) in the calpain (CAPN) and calpastatin (CAST) genes, described from Bos primigenius taurus, are polymorphic in Nellore cattle. We also looked for a possible association of linkage disequilibrium of this polymorphism with tenderness of the longissimus dorsi muscle after 7, 14 and 21 days of postmortem aging in 638 purebred Nellore bulls. Meat tenderness was measured as Warner-Bratzler shear force. Additive and dominance effects were tested for SNPs of the three genotypic classes; the substitution effect was tested for SNPs with missing genotypic classes. Genotypic and gene frequencies were also calculated for the different SNPs. An increase in tenderness was observed from 7 to 21 days; the average values for shear force at 7, 14 and 21 days of aging were 5.92 +/- 0.06, 4.92 +/- 0.05, and 4.38 +/- 0.04 kg, respectively. All markers showed polymorphism, but there was no CC genotype for CAPN316, and few animals showed the AA genotype for CAPN530. The alleles CAPN4751, UOGCAST1, and WSUCAST were found to have additive and dominance effects for shear force at 7, 14 and 21 days, while CAPN316 showed a substitution effect for shear force at 7 and 21 days. An additive-by-additive epistatic interaction was observed between CAPN4751 and markers on the CAST gene. In conclusion, these markers should be considered for use in breeding programs.

Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence

Background: The malaria parasite Plasmodium falciparum exhibits abundant genetic diversity, and this diversity is key to its success as a pathogen. Previous efforts to study genetic diversity in P. falciparum have begun to elucidate the demographic history of the species, as well as patterns of population structure and patterns of linkage disequilibrium within its genome. Such studies will be greatly enhanced by new genomic tools and recent large-scale efforts to map genomic variation. To that end, we have developed a high throughput single nucleotide polymorphism (SNP) genotyping platform for P. falciparum. Results: Using an Affymetrix 3,000 SNP assay array, we found roughly half the assays (1,638) yielded high quality, 100% accurate genotyping calls for both major and minor SNP alleles. Genotype data from 76 global isolates confirm significant genetic differentiation among continental populations and varying levels of SNP diversity and linkage disequilibrium according to geographic location and local epidemiological factors. We further discovered that nonsynonymous and silent (synonymous or noncoding) SNPs differ with respect to within-population diversity, interpopulation differentiation, and the degree to which allele frequencies are correlated between populations. Conclusions: The distinct population profile of nonsynonymous variants indicates that natural selection has a significant influence on genomic diversity in P. falciparum, and that many of these changes may reflect functional variants deserving of follow-up study. Our analysis demonstrates the potential for new high-throughput genotyping technologies to enhance studies of population structure, natural selection, and ultimately enable genome-wide association studies in P. falciparum to find genes underlying key phenotypic traits.

Genetic variability and natural selection at the ligand domain of the Duffy binding protein in brazilian Plasmodium vivax populations

Background: Plasmodium vivax malaria is a major public health challenge in Latin America, Asia and Oceania, with 130-435 million clinical cases per year worldwide. Invasion of host blood cells by P. vivax mainly depends on a type I membrane protein called Duffy binding protein (PvDBP). The erythrocyte-binding motif of PvDBP is a 170 amino-acid stretch located in its cysteine-rich region II (PvDBP(II)), which is the most variable segment of the protein. Methods: To test whether diversifying natural selection has shaped the nucleotide diversity of PvDBP(II) in Brazilian populations, this region was sequenced in 122 isolates from six different geographic areas. A Bayesian method was applied to test for the action of natural selection under a population genetic model that incorporates recombination. The analysis was integrated with a structural model of PvDBP(II), and T-and B-cell epitopes were localized on the 3-D structure. Results: The results suggest that: (i) recombination plays an important role in determining the haplotype structure of PvDBP(II), and (ii) PvDBP(II) appears to contain neutrally evolving codons as well as codons evolving under natural selection. Diversifying selection preferentially acts on sites identified as epitopes, particularly on amino acid residues 417, 419, and 424, which show strong linkage disequilibrium. Conclusions: This study shows that some polymorphisms of PvDBP(II) are present near the erythrocyte-binding domain and might serve to elude antibodies that inhibit cell invasion. Therefore, these polymorphisms should be taken into account when designing vaccines aimed at eliciting antibodies to inhibit erythrocyte invasion.

Modeling associations between genetic markers using Bayesian networks

Motivation: Understanding the patterns of association between polymorphisms at different loci in a population ( linkage disequilibrium, LD) is of fundamental importance in various genetic studies. Many coefficients were proposed for measuring the degree of LD, but they provide only a static view of the current LD structure. Generative models (GMs) were proposed to go beyond these measures, giving not only a description of the actual LD structure but also a tool to help understanding the process that generated such structure. GMs based in coalescent theory have been the most appealing because they link LD to evolutionary factors. Nevertheless, the inference and parameter estimation of such models is still computationally challenging. Results: We present a more practical method to build GM that describe LD. The method is based on learning weighted Bayesian network structures from haplotype data, extracting equivalence structure classes and using them to model LD. The results obtained in public data from the HapMap database showed that the method is a promising tool for modeling LD. The associations represented by the learned models are correlated with the traditional measure of LD D`. The method was able to represent LD blocks found by standard tools. The granularity of the association blocks and the readability of the models can be controlled in the method. The results suggest that the causality information gained by our method can be useful to tell about the conservability of the genetic markers and to guide the selection of subset of representative markers.

Effects of Perilipin (PLIN) Gene Variation on Metabolic Syndrome Risk and Weight Loss in Obese Children and Adolescents

Context: Genetic polymorphisms at the perilipin (PLIN) locus have been investigated for their potential utility as markers for obesity and metabolic syndrome (MS). We examined in obese children and adolescents (OCA) aged 7-14 yr the association of single-nucleotide polymorphisms (SNP) at the PLIN locus with anthropometric, metabolic traits, and weight loss after 20-wk multi-disciplinary behavioral and nutritional treatment without medication. Design: A total of 234 OCA [body mass index (BMI = 30.4 +/- 4.4 kg/m(2); BMI Z-score = 2.31 +/- 0.4) were evaluated at baseline and after intervention. We genotyped four SNPs (PLIN1 6209T -> C, PLIN4 11482G -> A, PLIN5 13041A -> G, and PLIN6 14995A -> T). Results: Allele frequencies were similar to other populations, PLIN1 and PLIN4 were in linkage disequilibrium (D` = 0.999; P < 0.001). At baseline, no anthropometric differences were observed, but minor allele A at PLIN4 was associated with higher triglycerides (111 +/- 49 vs. 94 +/- 42 mg/dl; P = 0.003), lower high-density lipoprotein cholesterol (40 +/- 9 vs. 44 +/- 10 mg/dl; P = 0.003) and higher homeostasis model assessment for insulin resistance (4.0 +/- 2.3 vs. 3.5 +/- 2.1; P +/- 0.015). Minor allele A at PLIN4 was associated with MS risk (age and sex adjusted) hazard ratio 2.4 (95% confidence interval = 1.1-4.9) for genotype GA and 3.5 (95% confidence interval = 1.2-9.9) for AA. After intervention, subjects carrying minor allele T at PLIN6 had increased weight loss (3.3 +/- 3.7 vs. 1.9 +/- 3.4 kg; P = 0.002) and increased loss of the BMI Z-score (0.23 +/- 0.18 vs. 0.18 +/- 0.15; P +/- 0.003). Due to group size, risk of by-chance findings cannot be excluded. Conclusion: The minor A allele at PLIN4 was associated with higher risk of MS at baseline, whereas the PLIN6 SNP was associated with better weight loss, suggesting that these polymorphisms may predict outcome strategies based on multidisciplinary treatment for OCA. (J Clin Endocrinol Metab 93: 4933-4940, 2008)

Polymorphisms of APOE and LRP genes in Brazilian individuals with Alzheimer disease

Alzheimer disease (AD) is the most frequent cause of dementia in Western countries. Putative genetic risk factors for AD are polymorphisms in the apolipoprotein E (APOE) gene and in the low-density lipoprotein receptor-related protein (LRP) gene. Our objective was to investigate the role of the APOE coding region polymorphisms epsilon 2, epsilon 3, and epsilon 4 and APOE promoter variants A/T at position -491 and G/T at -219, as well as LRP polymorphism C/T, as risk factors for AD in Brazilian individuals. One hundred and twenty patients with probable AD, along with 120 controls were analyzed. A significant difference between patients and controls for 64 alleles was observed: frequency of this allele in AD was 0.31, and 0.10 in controls. Individuals with 2 FA alleles had a higher risk for AD than subjects with only 1 such allele; presence of 1 epsilon 2 allele proved protective. The presence of the T allele of the -219 polymorphism was also associated with an increased risk of AD, but this polymorphism is in linkage disequilibrium with APOE F polymorphisms. No significant differences between patients and controls were observed for -491 APOE or LRP polymorphisms. In this Brazilian population, both the epsilon 4 allele and T -219 polymorphism were associated with an increased risk for AD.

Further delineation of the 15q13 microdeletion and duplication syndromes: a clinical spectrum varying from non-pathogenic to a severe outcome

Background: Recurrent 15q13.3 microdeletions were recently identified with identical proximal (BP4) and distal (BP5) breakpoints and associated with mild to moderate mental retardation and epilepsy. Methods: To assess further the clinical implications of this novel 15q13.3 microdeletion syndrome, 18 new probands with a deletion were molecularly and clinically characterised. In addition, we evaluated the characteristics of a family with a more proximal deletion between BP3 and BP4. Finally, four patients with a duplication in the BP3-BP4-BP5 region were included in this study to ascertain the clinical significance of duplications in this region. Results: The 15q13.3 microdeletion in our series was associated with a highly variable intra-and inter-familial phenotype. At least 11 of the 18 deletions identified were inherited. Moreover, 7 of 10 siblings from four different families also had this deletion: one had a mild developmental delay, four had only learning problems during childhood, but functioned well in daily life as adults, whereas the other two had no learning problems at all. In contrast to previous findings, seizures were not a common feature in our series (only 2 of 17 living probands). Three patients with deletions had cardiac defects and deletion of the KLF13 gene, located in the critical region, may contribute to these abnormalities. The limited data from the single family with the more proximal BP3-BP4 deletion suggest this deletion may have little clinical significance. Patients with duplications of the BP3-BP4-BP5 region did not share a recognisable phenotype, but psychiatric disease was noted in 2 of 4 patients. Conclusions: Overall, our findings broaden the phenotypic spectrum associated with 15q13.3 deletions and suggest that, in some individuals, deletion of 15q13.3 is not sufficient to cause disease. The existence of microdeletion syndromes, associated with an unpredictable and variable phenotypic outcome, will pose the clinician with diagnostic difficulties and challenge the commonly used paradigm in the diagnostic setting that aberrations inherited from a phenotypically normal parent are usually without clinical consequences.

HLA-G polymorphism and transitional cell carcinoma of the bladder in a Brazilian population

The morphologic appearance and clinical behavior of the human urinary bladder papillary transitional cell carcinoma (TCC) probably result from a complex interaction between carcinogenic insults and host resistance during the patient`s life. While the main recognized risk factors are of environmental origin (e.g. smoking), relatively little information exists about the susceptibility to TCC development. The human leukocyte antigen G (HLA-G) molecule plays an important role in immune response regulation and has been implicated in the inhibition of the cytolytic function of natural killer and cytotoxic T cells. Several lines of evidence indicate that HLA-G polymorphisms influence the expression level and production of different HLA-G isoforms. The aim of this study was to explore a possible influence of the HLA-G polymorphism on the susceptibility to urinary bladder TCC development and progression in smokers and nonsmokers Brazilian subjects. The HLA-G locus was found to be associated with susceptibility to TCC development and progression. The G*0104 allelic group (specially the G*010404 allele) and the G*0103 allele were associated with a tobacco-dependent influence on TCC development. The G*0104 group was associated with progression to high-grade tumors, irrespective of smoking habit, while the G*0103 allele was associated to high-grade tumor only in smoking patients. Our results are an evidence that the HLA-G locus itself, or as part of an extended haplotype encompassing this chromosome region (particularly the HLA-A given the high linkage disequilibrium observed between them in this data series), may be associated with TCC susceptibility and tumor progression, suggesting a tobacco-dependent influence of these polymorphisms.

Evaluation of computational methods for the reconstruction of HLA haplotypes

Human leukocyte antigen (HLA) haplotypes are frequently evaluated for population history inferences and association studies. However, the available typing techniques for the main HLA loci usually do not allow the determination of the allele phase and the constitution of a haplotype, which may be obtained by a very time-consuming and expensive family-based segregation study. Without the family-based study, computational inference by probabilistic models is necessary to obtain haplotypes. Several authors have used the expectation-maximization (EM) algorithm to determine HLA haplotypes, but high levels of erroneous inferences are expected because of the genetic distance among the main HLA loci and the presence of several recombination hotspots. In order to evaluate the efficiency of computational inference methods, 763 unrelated individuals stratified into three different datasets had their haplotypes manually defined in a family-based study of HLA-A, -B, -DRB1 and -DQB1 segregation, and these haplotypes were compared with the data obtained by the following three methods: the Expectation-Maximization (EM) and Excoffier-Laval-Balding (ELB) algorithms using the arlequin 3.11 software, and the PHASE method. When comparing the methods, we observed that all algorithms showed a poor performance for haplotype reconstruction with distant loci, estimating incorrect haplotypes for 38%-57% of the samples considering all algorithms and datasets. We suggest that computational haplotype inferences involving low-resolution HLA-A, HLA-B, HLA-DRB1 and HLA-DQB1 haplotypes should be considered with caution.

The genetic structure of 3 ` untranslated region of the HLA-G gene: polymorphisms and haplotypes

The HLA-G gene is predominantly expressed at the maternal-fetal interface. It has been associated with maternal-fetal tolerance and in the inhibition of cytotoxic T lymphocyte and natural killer cytolytic functions. At least two variations in the 3` untranslated region (UTR) of HLA-G locus are associated with HLA-G expression levels, the 14-bp deletion/insertion polymorphism and the +3142 single-nucleotide polymorphism (SNP). However, this region has not been completely characterized yet. The variability of the 3`UTR of HLA-G gene and its haplotype structure were characterized in 155 individuals from Brazil, as well as HLA-G alleles associated with each of the 3`UTR haplotype. The following eight variation sites were detected: the 14-bp polymorphism and SNPs at the positions +3003T/C, +3010C/G, +3027A/C, +3035C/T, +3142G/C, +3187A/G and +3196C/G. Similarly, 11 different 3`UTR haplotypes were identified and several HLA-G alleles presented only one 3`UTR haplotype. In addition, a high linkage disequilibrium among the variation sites was detected, especially among the 14-bp insertion and the alleles +3142G and +3187A, all previously associated with low mRNA availability, demonstrating that their effects are not independent. The detailed analyses of 3`UTR of the HLA-G locus may shed some light into mechanisms underlying the regulation of HLA-G expression. Genes and Immunity (2010) 11, 134-141; doi: 10.1038/gene.2009.74; published online 1 October 2009