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.

Genetic Contribution for Non-Syndromic Cleft Lip With or Without Cleft Palate (NS CL/P) in Different Regions of Brazil and Implications for Association Studies

Non-syndromic cleft lip with or without cleft palate (NS CL/P) is a complex disease in which heritability estimates vary widely depending on the population studied. To evaluate the importance of genetic contribution to NS CL/P in the Brazilian population, we conducted a study with 1,042 families from five different locations (Santarem, Fortaleza, Barbalha, Maceio, and Rio de Janeiro). We also evaluated the role of consanguinity and ethnic background. The proportion of familial cases varied significantly across locations, with the highest values found in Santarem (44%) and the lowest in Maceio (23%). Heritability estimates showed a higher genetic contribution to NS CL/P in Barbalha (85%), followed by Santarem (71%), Rio de Janeiro (70%), Fortaleza (64%), and Maceio (45%). Ancestry was not correlated with the occurrence of NS CL/P or with the variability in heritability. Only in Rio de Janeiro was the coefficient of inbreeding significantly larger in NS CL/P families than in the local population. Recurrence risk for the total sample was approximately 1.5-1.6%, varying according to the location studied (0.6-0.7% in Maceio to 2.2-2.8% in Barbalha). Our findings show that the degree of genetic contribution to NS CL/P varies according to the geographic region studied, and this difference cannot be attributed to consanguinity or ancestry. These findings suggest that Barbalha is a promising region for genetic studies. The data presented here will be useful in interpreting results from molecular analyses and show that care must be taken when pooling samples from different populations for association studies. (C) 2011 Wiley-Liss, Inc.

Association of polymorphisms at the ADIPOR1 regulatory region with type 2 diabetes and body mass index in a Brazilian population with European or African ancestry

Association studies between ADIPOR1 genetic variants and predisposition to type 2 diabetes (DM2) have provided contradictory results. We determined if two single nucleotide polymorphisms (SNP c.-8503G>A and SNP c.10225C>G) in regulatory regions of ADIPOR1 in 567 Brazilian individuals of European (EA; N = 443) or African (AfA; N = 124) ancestry from rural (quilombo remnants; N = 439) and urban (N = 567) areas. We detected a significant effect of ethnicity on the distribution of the allelic frequencies of both SNPs in these populations (EA: -8503A = 0.27; AfA: -8503A = 0.16; P = 0.001 and EA: 10225G = 0.35; AfA: 10225G = 0.51; P < 0.001). Neither of the polymorphisms were associated with DM2 in the case-control study in EA (SNP c.-8503G>A: DM2 group -8503A = 0.26; control group -8503A = 0.30; P = 0.14/SNP 10225C>G: DM2 group 10225G = 0.37; control group 10225G = 0.32; P = 0.40) and AfA populations (SNP c.-8503G>A: DM2 group -8503A = 0.16; control group -8503A = 0.15; P = 0.34/SNP 10225C>G: DM2 group 10225G = 0.51; control group 10225G = 0.52; P = 0.50). Similarly, none of the polymorphisms were associated with metabolic/anthropometric risk factors for DM2 in any of the three populations, except for HDL cholesterol, which was significantly higher in AfA heterozygotes (GC = 53.75 ± 17.26 mg/dL) than in homozygotes. We conclude that ADIPOR1 polymorphisms are unlikely to be major risk factors for DM2 or for metabolic/anthropometric measurements that represent risk factors for DM2 in populations of European and African ancestries.

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.

Genetic Background of Patients from a University Medical Center in Manhattan: Implications for Personalized Medicine

Background: The rapid progress currently being made in genomic science has created interest in potential clinical applications; however, formal translational research has been limited thus far. Studies of population genetics have demonstrated substantial variation in allele frequencies and haplotype structure at loci of medical relevance and the genetic background of patient cohorts may often be complex. Methods and Findings: To describe the heterogeneity in an unselected clinical sample we used the Affymetrix 6.0 gene array chip to genotype self-identified European Americans (N = 326), African Americans (N = 324) and Hispanics (N = 327) from the medical practice of Mount Sinai Medical Center in Manhattan, NY. Additional data from US minority groups and Brazil were used for external comparison. Substantial variation in ancestral origin was observed for both African Americans and Hispanics; data from the latter group overlapped with both Mexican Americans and Brazilians in the external data sets. A pooled analysis of the African Americans and Hispanics from NY demonstrated a broad continuum of ancestral origin making classification by race/ethnicity uninformative. Selected loci harboring variants associated with medical traits and drug response confirmed substantial within-and between-group heterogeneity. Conclusion: As a consequence of these complementary levels of heterogeneity group labels offered no guidance at the individual level. These findings demonstrate the complexity involved in clinical translation of the results from genome-wide association studies and suggest that in the genomic era conventional racial/ethnic labels are of little value.

Prognostic value of TP53 Pro47Ser and Arg72Pro single nucleotide polymorphisms and the susceptibility to gliomas in individuals from Southeast Brazil

The TP53 tumor suppressor gene codifies a protein responsible for preventing cells with genetic damage from growing and dividing by blocking cell growth or apoptosis pathways. A common single nucleotide polymorphism (SNP) in TP53 codon 72 (Arg72Pro) induces a 15-fold decrease of apoptosis-inducing ability and has been associated with susceptibility to human cancers. Recently, another TP53 SNP at codon 47 (Pro47Ser) was reported to have a low apoptosis-inducing ability; however, there are no association studies between this SNP and cancer. Aiming to study the role of TP53 Pro47Ser and Arg72Pro on glioma susceptibility and oncologic prognosis of patients, we investigated the genotype distribution of these SNPs in 94 gliomas (81 astrocytomas, 8 ependymomas and 5 oligodendrogliomas) and in 100 healthy subjects by the polymerase chain reaction-restriction fragment length polymorphism approach. Chi-square and Fisher exact test comparisons for genotype distributions and allele frequencies did not reveal any significant difference between patients and control groups. Overall and disease-free survivals were calculated by the Kaplan-Meier method, and the log-rank test was used for comparisons, but no significant statistical difference was observed between the two groups. Our data suggest that TP53 Pro47Ser and Arg72Pro SNPs are not involved either in susceptibility to developing gliomas or in patient survival, at least in the Brazilian population.

Endothelial Nitric Oxide Synthase Polymorphisms and Haplotypes in Amerindians

Interethnic disparities in the distribution of endothelial nitric oxide synthase (eNOS) polymorphisms may affect nitric oxide (NO)-mediated effects of and responses to drugs. While there are differences between black and white subjects there is no information regarding the distribution of eNOS gene alleles and haplotypes in Amerindians. We studied three clinically relevant eNOS polymorphisms (T(-786) C in the promoter, a variable number of tandem repeats in intron 4, and the Glu298Asp in exon 7) and eNOS haplotypes in 170 Amerindians from three tribes of the Brazilian Amazon. The results were compared with previous findings for black and white Brazilians. The Asp298, C(-786), and 4a alleles were much less common in Amerindians (5.0%, 3.2%, and 4.1%, respectively) than in blacks (15.1%, 19.5%, and 32.0%, respectively) or whites (32.8%, 41.9%, and 17.9%, respectively) (p<0.001). The haplotype including the most common alleles for each polymorphism was much more common in Amerindians (89%) than in blacks (45%) or whites (41%). Our findings are consistent with a lower genetic diversity in Amerindians compared with blacks and whites. These striking differences may be of major relevance for case-control association studies focusing on eNOS gene polymorphisms and may explain, at least in part, differences in the responses to cardiovascular drugs.

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.

Detection of the TLR4 1196C > T Polymorphism by Mismatched-Polymerase Chain Reaction Using Plasmid DNA as Internal Control in Restriction Fragment Length Polymorphism Assays

Background: Restriction fragment length polymorphism (RFLP) is a common molecular assay used for genotyping, and it requires validated quality control procedures to prevent mistyping caused by impaired endonuclease activity. We have evaluated the usefulness of a plasmid-based internal control in RFLP assays. Results: Blood samples were collected from 102 individuals with acute myocardial infarction (AMI) and 108 non-AMI individuals (controls) for DNA extraction and laboratory analyses. The 1196C> T polymorphism in the toll-like receptor 4 (TLR4) gene was amplified by mismatched-polymerase chain reaction (PCR). Amplicons and pBluescript II SK-plasmid were simultaneously digested with endonuclease HincII. Fragments were separated on 2% agarose gels. Plasmid was completely digested using up to 55.2 nmL/L DNA solutions and 1 mu L PCR product. Nevertheless, plasmid DNA with 41.4 nM or higher concentrations was incompletely digested in the presence of 7 mL PCR product. In standardized conditions, TLR4 1196C> T variant was accurately genotyped. TLR4 1196T allele frequency was similar between AMI (3.1%) and controls (2.0%, p = 0.948). TLR4 SNP was not associated with AMI in this sample population. In conclusion, the plasmid-based control is a useful approach to prevent mistyping in RFLP assays, and it is validate for genetic association studies such as TLR4 1196C> T.

Strategies for genetic model specification in the screening of genome-wide meta-analysis signals for further replication

Background Meta-analysis is increasingly being employed as a screening procedure in large-scale association studies to select promising variants for follow-up studies. However, standard methods for meta-analysis require the assumption of an underlying genetic model, which is typically unknown a priori. This drawback can introduce model misspecifications, causing power to be suboptimal, or the evaluation of multiple genetic models, which augments the number of false-positive associations, ultimately leading to waste of resources with fruitless replication studies. We used simulated meta-analyses of large genetic association studies to investigate naive strategies of genetic model specification to optimize screenings of genome-wide meta-analysis signals for further replication. Methods Different methods, meta-analytical models and strategies were compared in terms of power and type-I error. Simulations were carried out for a binary trait in a wide range of true genetic models, genome-wide thresholds, minor allele frequencies (MAFs), odds ratios and between-study heterogeneity (tau(2)). Results Among the investigated strategies, a simple Bonferroni-corrected approach that fits both multiplicative and recessive models was found to be optimal in most examined scenarios, reducing the likelihood of false discoveries and enhancing power in scenarios with small MAFs either in the presence or in absence of heterogeneity. Nonetheless, this strategy is sensitive to tau(2) whenever the susceptibility allele is common (MAF epsilon 30%), resulting in an increased number of false-positive associations compared with an analysis that considers only the multiplicative model. Conclusion Invoking a simple Bonferroni adjustment and testing for both multiplicative and recessive models is fast and an optimal strategy in large meta-analysis-based screenings. However, care must be taken when examined variants are common, where specification of a multiplicative model alone may be preferable.

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.

Assessing Individual Interethnic Admixture and Population Substructure Using a 48-Insertion-Deletion (INSEL) Ancestry-Informative Marker (AIM) Panel

Estimating the proportions of different ancestries in admixed populations is very important in population genetics studies, and it is particularly important for detecting population substructure effects in case-control association studies. In this work, a set of 48 ancestry, informative insertion, deletion polymorphisms (INDELs) were selected with the goal of efficiently measuring the proportions of three different ancestries (sub-Saharan African, European, and Native American) in mixed populations. All selected markers can be easily analyzed via multiplex PCR and detected with standard capillary electrophoresis. A total of 593 unrelated individuals representative of European, African, and Native American parental populations were typed, as were 380 individuals from three Brazilian populations with known admixture patterns. As expected, the interethnic admixture estimates show that individuals from southern Brazil present an almost exclusively European ancestry; Afro-descendant communities in the Amazon region, apart from the major African contribution, present some degree of admixture with Europeans and Native Americans; and a sample from Belem, in the northeastern Amazon, shows a significant contribution of the three ethnic groups, although with a greater European proportion. In summary, a panel of ancestry-informative INDELs was optimized and proven to be a variable tool for estimating individual and global ancestry proportions in admixed populations. The ability to accurately infer interethnic admixtures highlights the usefulness of this marker set for assessing population substructure in association studies, particularly those conducted in Brazilian and other Latin American populations sharing trihybrid ancestry patterns. Hum Mutat 31:184-190, 2010. (C) 2009 Wiley-Liss, Inc.

T allele of-344C/T polymorphism in aldosterone synthase gene is not associated with resistant hypertension

Resistant hypertension (RH) is the maintenance of elevated blood pressure concurrent with the use of three different antihypertensive drugs, one of which is a diuretic. The Renin-Angiotensin-Aldosterone System plays a major role in volume-dependent hypertension. Therefore, its components are interesting targets for genetic association studies. This work focused on the -344 C/T polymorphism in the CYP11b2 gene, which encodes aldosterone synthase. This work evaluates the association between T allele and resistance to anti-hypertensive treatment. Genotyping analysis included 88 subjects with RH, 142 who were responsive to anti-hypertensive treatment and 110 subjects as a control group. Plasmatic concentrations of aldosterone, renin and cortisol, carotid intima-media thickness and carotid-femoral pulse wave velocity were assessed in a smaller subset of hypertensive patients. An association was found between T allele and hypertension (P < 0.005), but there was no difference in allele frequencies between both hypertensive groups. There was no difference in plasmatic parameters either, in remodeling indicators between the genotypic groups.

Protector mechanisms of the association between gastroesophageal reflux disease and asthma: experimental study in rats

BACKGROUND: It is well known the association between gastroesophageal reflux disease and asthma. The hyperreactivity of the airways is a characteristic of an asthmatic. Many studies associate the increase of the airways reactivity with gastroesophageal reflux disease. AIM: In this study we have evaluated the effect of the intraluminal exposition to gastric juice of trachea on the reactivity to methacholine from rats submitted to a pulmonary allergic inflammation. METHODS: Group of rats were sensitized and challenged with ovalbumin. After 24 hours the animals were sacrificed, and their tracheae were removed to be cultured with gastric juice. The gastric juice was obtained from a donor rat. Subsequently the segments were placed into plastic plates with RPMI-1640 for incubation, under suitable atmosphere and time. After the period of incubation the segments were put into chambers for the analysis of the contractile response to methacholine. RESULTS: We observed reduction in the contractile response of trachea cultured with gastric juice from allergic rats. This result was confirmed by the pharmacological treatments with compound 48/80 and dissodium cromoglicate (mast cells blockade), L-NAME (nitric oxide inhibitor, NO), capsaicin (neuropeptides depletion) and indomethacin (ciclooxigenase inhibitor). CONCLUSIONS: Our results highlight to the existence of a complex interaction between pulmonary allergy and gastric juice in the airways. The involvement of the non-adrenergic non-cholinergic system, NO, prostanoids and mast cells are directly related to this interaction. We suggest that the reduced contractile response observed in vitro may represent a protector mechanism of the airways. Despite its presence in the human body it can not be observed due to the predominant effects of excitatory the non-adrenergic non-cholinergic system.