Mapping Bias Overestimates Reference Allele Frequencies at the HLA Genes in the 1000 Genomes Project Phase I Data.

Next-generation sequencing (NGS) technologies have become the standard for data generation in studies of population genomics, as the 1000 Genomes Project (1000G). However, these techniques are known to be problematic when applied to highly polymorphic genomic regions, such as the human leukocyte antigen (HLA) genes. Because accurate genotype calls and allele frequency estimations are crucial to population genomics analyses, it is important to assess the reliability of NGS data. Here, we evaluate the reliability of genotype calls and allele frequency estimates of the single-nucleotide polymorphisms (SNPs) reported by 1000G (phase I) at five HLA genes (HLA-A, -B, -C, -DRB1, and -DQB1). We take advantage of the availability of HLA Sanger sequencing of 930 of the 1092 1000G samples and use this as a gold standard to benchmark the 1000G data. We document that 18.6% of SNP genotype calls in HLA genes are incorrect and that allele frequencies are estimated with an error greater than ±0.1 at approximately 25% of the SNPs in HLA genes. We found a bias toward overestimation of reference allele frequency for the 1000G data, indicating mapping bias is an important cause of error in frequency estimation in this dataset. We provide a list of sites that have poor allele frequency estimates and discuss the outcomes of including those sites in different kinds of analyses. Because the HLA region is the most polymorphic in the human genome, our results provide insights into the challenges of using of NGS data at other genomic regions of high diversity.

Allele and genotype frequencies of polymorphic cytochromes P4502D6, 2C19 and 2E1 in Aborigines from Western Australia

The polymorphisms of the important xenobiotic metabolizing enzymes CYP2D6, CYP2C19 and CYP2E1 have been studied extensively in a large number of populations and show significant heterogeneity in the frequency of different alleles/genotypes and in the prevalence of the extensive and poor metabolizer phenotypes, Understanding of inter-ethnic differences in genotypes is important in prediction of either beneficial or adverse effects from therapeutic agents and other xenobiotics. Since no data were available for Australian Aborigines, we investigated the frequencies of alleles and genotypes for CYP2D6, CYP2C19 and CYP2E1 in a population living in the far north of Western Australia. Because of its geographical isolation, this population can serve as a model to study the impact of evolutionary forces on the distribution of different alleles for xenobiotic metabolizing enzymes. Twelve CYP2D6 alleles were analysed, The wild-type allele *1 was the most frequent (85.8%) and the non-functional alleles (*4, *5, *16) had an overall frequency of less than 10%. Only one subject (0.4%) was a poor metabolizer for CYP2D6 because of the genotype *5/*5, For CYP2C19, the frequencies of the *1 (wild-type) and the non-functional (*2 and *3) alleles were 50.2%, 35.5% and 14.3%, respectively. The combined CYP2C19 genotypes (*2/*2, *2/*3 or *3/*3) correspond to a predicted frequency of 25.6% for the CYP2C19 poor metabolizer phenotype, For CYP2E1, only one subject had the rare c2 allele giving an overall allele frequency of 0.2%. For CYP2D6 and CYP2C19, allele frequencies and predicted phenotypes differed significantly from those for Caucasians but were similar to those for Orientals indicating a close relationship to East Asian populations. Differences between Aborigines and Orientals in allele frequencies for CYP2D6*10 and CYP2E1 c2 may have arisen through natural selection, or genetic drift, respectively, Pharmacogenetics 11:69-76 (C) 2001 Lippincott Williams & Wilkins.

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.

Absence of the E2 allele of apolipoprotein in Amerindians

Determination of the ApoE allele distribution in five South American Amerindian tribes revealed absence of the ApoE2 allele, accompanied by high ApoE3 and low ApoE4 allele frequencies for most tribes, a distribution only previously reported for the Inuit Eskimo from Greenland.

Allele frequency distributions of six hypervariable loci (D1S80, APOB, D4S43, vW1, F13A and DYS19) in two African-Brazilian communities from the Amazon region

ABSTRACT: The allele frequency distributions of three VNTR (D1S80, APOB and D4S43) and three STR (vW1, F13A1 and DYS19) loci were investigated in two Afro-Brazilian populations from the Amazon: Curiau and Pacoval. Exact tests for population differentiation revealed significant differences in allele frequency between populations only for the D1S80 and APOB loci. A statistically significant deviation from the Hardy-Weinberg equilibrium was observed only in the D1S80 locus of the Pacoval sample. A neighbor-joining tree was constructed based on DA genetic distances of allele frequencies in four Afro-Brazilian populations from the Amazon (Pacoval, Curiau, Trombetas, and Cametá), along with those from Congo, Cameroon, Brazilian Amerindians, and Europeans. This analysis revealed the usefulness of these Amp-FLPs for population studies - African and African-derived populations were closely grouped, and clearly separated from Amerindians and Europeans. Estimates of admixture components based on the gene identity method revealed the prevalence of the African component in both populations studied, amounting to 51% in Pacoval, and to 43% in Curiau. The Amerindian component was also important in both populations (37% in Pacoval, and 24% in Curiau). The European component reached 33% in Curiau.

Allelic frequencies and statistical data obtained from 12 codis STR loci in an admixed population of the Brazilian Amazon

The allelic frequencies of 12 short tandem repeat loci were obtained from a sample of 307 unrelated individuals living in Macapá, a city in the northern Amazon region, Brazil. These loci are the most commonly used in forensics and paternity testing. Based on the allele frequency obtained for the population of Macapá, we estimated an interethnic admixture for the three parental groups (European, Native American and African) of, respectively, 46%, 35% and 19%. Comparing these allele frequencies with those of other Brazilian populations and of the Iberian Peninsula population, no significant distances were observed. The interpopulation genetic distances (F_ST coefficients) to the present database ranged from F_ST = 0.0016 between Macapá and Belém to F_ST = 0.0036 between Macapá and the Iberian Peninsula.

Distribution of the CCR5delta32 allele (gene variant CCR5) in Rondonia, Western Amazonian region, Brazil

Since around 1723, on the occasion of its initial colonization by Europeans, Rondonia has received successive waves of immigrants. This has been further swelled by individuals from northeastern Brazil, who began entering at the beginning of the twentieth century. The ethnic composition varies across the state according to the various sites of settlement of each wave of immigrants. We analyzed the frequency of the CCR5 Delta 32 allele of the CCR5 chemokine receptor, which is considered a Caucasian marker, in five sample sets from the population. Four were collected in Porto Velho, the state capital and the site of several waves of migration. Of these, two, from the Hospital de Base were comprised of HB Mothers and HB Newborns presenting allele frequencies of 3.5% and 3.1%, respectively, a third from the peri-urban neighborhoods of Candelaria/Bate-Estaca (1.8%), whereas a fourth, from the Research Center on Tropical Medicine/CEPEM (0.6%), was composed of malaria patients under treament. The fifth sample (3.4%) came from the inland Quilombola village of Pedras Negras. Two homozygous individuals (CCR5 Delta 32/CCR5 Delta 32) were detected among the HB Mother samples. The frequency of this allele was heterogeneous and higher where the European inflow was more pronounced. The presence of the allele in Pedras Negras revealed European miscegenation in a community largely comprising Quilombolas.

Associations between prefrontal γ-aminobutyric acid concentration and the tryptophan hydroxylase isoform 2 gene, a panic disorder risk allele in women

Associations between the central serotonergic and γ-aminobutyric acid (GABA) systems play key roles in the prefrontal cortical regulation of emotion and cognition and in the pathophysiology and pharmacotherapy of highly prevalent psychiatric disorders. The goal of this study was to test the effects of common variants of the tryptophan hydroxylase isoform 2 (TPH2) gene on GABA concentration in the prefrontal cortex (PFC) using magnetic resonance spectroscopy. In this study involving 64 individuals, we examined the associations between prefrontal cortical GABA concentration and 12 single nucleotide polymorphisms (SNPs) spanning the TPH2 gene, including rs4570625 (−703 G/T SNP), a potentially functional TPH2 polymorphism that has been associated with decreased TPH2 mRNA expression and panic disorder. Our results revealed a significant association between increased GABA concentration in the PFC and the T-allele frequencies of two TPH2 SNPs, namely rs4570625 (−703 G/T) and rs2129575 (p≤0.0004) and the C-allele frequency of one TPH2 SNP, namely rs1386491 (p = 0.0003) in female subjects. We concluded that rs4570625 (−703 G/T), rs2129575 and rs1386491 play a significant role in GABAergic neurotransmission and may contribute to the sex-specific dysfunction of the GABAergic system in the PFC.

Bovine gene polymorphisms related to fat deposition and meat tenderness

Leptin, thyroglobulin and diacylglycerol O-acyltransferase play important roles in fat metabolism. Fat deposition has an influence on meat quality and consumers' choice. The aim of this study was to determine allele and genotype frequencies of polymorphisms of the bovine genes, which encode leptin (LEP), thyroglobulin (TG) and diacylglycerol O-acyltransferase (DGAT1). A further objective was to establish the effects of these polymorphisms on meat characteristics. We genotyped 147 animals belonging to the Nelore (Bos indicus), Canchim (5/8 Bos taurus + 3/8 Bos indicus), Rubia Gallega X Nelore (1/2 Bos taurus + 1/2 Bos indicus), Brangus Three-way cross (9/16 Bos taurus + 7/16 Bos indicus) and Braunvieh Three-way cross (3/4 Bos taurus + 1/4 Bos indicus) breeds. Backfat thickness, total lipids, marbling score, ribeye area and shear force were fitted, using the General Linear Model (GLM) procedure of the SAS software. The least square means of genotypes and genetic groups were compared using Tukey's test. Allele frequencies vary among the genetic groups, depending on Bos indicus versus Bos taurus influence. The LEP polymorphism segregates in pure Bos indicus Nelore animals, which is a new finding. The T allele of TG is fixed in Nelore, and DGAT1 segregates in all groups, but the frequency of allele A is lower in Nelore animals. The results showed no association between the genotypes and traits studied, but a genetic group effect on these traits was found. So, the genetic background remains relevant for fat deposition and meat tenderness, but the gene markers developed for Bos taurus may be insufficient for Bos indicus.

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.

CYP2C19 and ABCB1 gene polymorphisms are differently distributed according to ethnicity in the Brazilian general population

Background: Recent studies have reported the clinical importance of CYP2C19 and ABCB1 polymorphisms in an individualized approach to clopidogrel treatment. The aims of this study were to evaluate the frequencies of CYP2C19 and ABCB1 polymorphisms and to identify the clopidogrel-predicted metabolic phenotypes according to ethnic groups in a sample of individuals representative of a highly admixtured population. Methods: One hundred and eighty-three Amerindians and 1,029 subjects of the general population of 4 regions of the country were included. Genotypes for the ABCB1c.C3435T (rs1045642), CYP2C19*2 (rs4244285), CYP2C19*3 (rs4986893), CYP2C19*4 (rs28399504), CYP2C19*5 (rs56337013), and CYP2C19*17 (rs12248560) polymorphisms were detected by polymerase chain reaction followed by high resolution melting analysis. The CYP2C19*3, CYP2C19*4 and CYP2C19*5 variants were genotyped in a subsample of subjects (300 samples randomly selected). Results: The CYP2C19*3 and CYP2C19*5 variant alleles were not detected and the CYP2C19*4 variant allele presented a frequency of 0.3%. The allelic frequencies for the ABCB1c.C3435T, CYP2C19*2 and CYP2C19*17 polymorphisms were differently distributed according to ethnicity: Amerindian (51.4%, 10.4%, 15.8%); Caucasian descent (43.2%, 16.9%, 18.0%); Mulatto (35.9%, 16.5%, 21.3%); and African descent (32.8%, 20.2%, 26.3%) individuals, respectively. As a result, self-referred ethnicity was able to predict significantly different clopidogrel-predicted metabolic phenotypes prevalence even for a highly admixtured population. Conclusion: Our findings indicate the existence of inter-ethnic differences in the ABCB1 and CYP2C19 variant allele frequencies in the Brazilian general population plus Amerindians. This information could help in stratifying individuals from this population regarding clopidogrel-predicted metabolic phenotypes and design more cost-effective programs towards individualization of clopidogrel therapy.

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.