A note on the use of the generalized odds ratio in meta-analysis of association studies involving bi- and tri-allelic polymorphisms

Abstract Background The generalized odds ratio (GOR) was recently suggested as a genetic model-free measure for association studies. However, its properties were not extensively investigated. We used Monte Carlo simulations to investigate type-I error rates, power and bias in both effect size and between-study variance estimates of meta-analyses using the GOR as a summary effect, and compared these results to those obtained by usual approaches of model specification. We further applied the GOR in a real meta-analysis of three genome-wide association studies in Alzheimer's disease. Findings For bi-allelic polymorphisms, the GOR performs virtually identical to a standard multiplicative model of analysis (e.g. per-allele odds ratio) for variants acting multiplicatively, but augments slightly the power to detect variants with a dominant mode of action, while reducing the probability to detect recessive variants. Although there were differences among the GOR and usual approaches in terms of bias and type-I error rates, both simulation- and real data-based results provided little indication that these differences will be substantial in practice for meta-analyses involving bi-allelic polymorphisms. However, the use of the GOR may be slightly more powerful for the synthesis of data from tri-allelic variants, particularly when susceptibility alleles are less common in the populations (≤10%). This gain in power may depend on knowledge of the direction of the effects. Conclusions For the synthesis of data from bi-allelic variants, the GOR may be regarded as a multiplicative-like model of analysis. The use of the GOR may be slightly more powerful in the tri-allelic case, particularly when susceptibility alleles are less common in the populations.

Testing allele homogeneity: the problem of nested hypotheses

Background: The evaluation of associations between genotypes and diseases in a case-control framework plays an important role in genetic epidemiology. This paper focuses on the evaluation of the homogeneity of both genotypic and allelic frequencies. The traditional test that is used to check allelic homogeneity is known to be valid only under Hardy-Weinberg equilibrium, a property that may not hold in practice. Results: We first describe the flaws of the traditional (chi-squared) tests for both allelic and genotypic homogeneity. Besides the known problem of the allelic procedure, we show that whenever these tests are used, an incoherence may arise: sometimes the genotypic homogeneity hypothesis is not rejected, but the allelic hypothesis is. As we argue, this is logically impossible. Some methods that were recently proposed implicitly rely on the idea that this does not happen. In an attempt to correct this incoherence, we describe an alternative frequentist approach that is appropriate even when Hardy-Weinberg equilibrium does not hold. It is then shown that the problem remains and is intrinsic of frequentist procedures. Finally, we introduce the Full Bayesian Significance Test to test both hypotheses and prove that the incoherence cannot happen with these new tests. To illustrate this, all five tests are applied to real and simulated datasets. Using the celebrated power analysis, we show that the Bayesian method is comparable to the frequentist one and has the advantage of being coherent. Conclusions: Contrary to more traditional approaches, the Full Bayesian Significance Test for association studies provides a simple, coherent and powerful tool for detecting associations.

Differentiation of African Components of Ancestry to Stratify Groups in a Case-Control Study of a Brazilian Urban Population

Background: Balancing the subject composition of case and control groups to create homogenous ancestries between each group is essential for medical association studies. Methods: We explored the applicability of single-tube 34-plex ancestry informative markers (AIM) single nucleotide polymorphisms (SNPs) to estimate the African Component of Ancestry (ACA) to design a future case-control association study of a Brazilian urban sample. Results: One hundred eighty individuals (107 case group; 73 control group) self-described as white, brown-intermediate or black were selected. The proportions of the relative contribution of a variable number of ancestral population components were similar between case and control groups. Moreover, the case and control groups demonstrated similar distributions for ACA <0.25 and >0.50 categories. Notably a high number of outlier values (23 samples) were observed among individuals with ACA <0.25. These individuals presented a high probability of Native American and East Asian ancestral components; however, no individuals originally giving these self-described ancestries were observed in this study. Conclusions: The strategy proposed for the assessment of ancestry and adjustment of case and control groups for an association study is an important step for the proper construction of the study, particularly when subjects are taken from a complex urban population. This can be achieved using a straight forward multiplexed AIM-SNPs assay of highly discriminatory ancestry markers.

MYLIP p.N342S polymorphism is not associated with lipid profile in the Brazilian population

Background: A recent study investigated the MYLIP region in the Mexican population in order to fine-map the actual susceptibility variants of this locus. The p.N342S polymorphism was identified as the underlying functional variant accounting for one of the previous signals of genome-wide association studies and the N342 allele was associated with higher cholesterol concentrations in Mexican dyslipidemic individuals. To date, there is no further evaluation on this genotype-phenotype association in the literature. In this scenario, and because of a possible pharmacotherapeutic target of dyslipidemia, the main aim of this study was to assess the influence of the MYLIP p.N342S polymorphism on lipid profile in Brazilian individuals. Methods: 1295 subjects of the general population and 1425 consecutive patients submitted to coronary angiography were selected. General characteristics, biochemical tests, blood pressures, pulse wave velocity, and coronary artery disease scores were analyzed. Genotypes for the MYLIP rs9370867 (p.N342S, c.G1025A) polymorphism were detected by high resolution melting analysis. Results: No association of the MYLIP rs9370867 genotypes with lipid profile, hemodynamic data, and coronary angiographic data was found. Analysis stratified by hyperlipidemia, gender, and ethnicity was also performed and the sub-groups presented similar results. In both general population and patient samples, the MYLIP rs9370867 polymorphism was differently distributed according to ethnicity. In the general population, subjects carrying GG genotypes had higher systolic blood pressure (BP), diastolic BP, and mean BP values (129.0 +/- 23.3; 84.9 +/- 14.6; 99.5 +/- 16.8 mmHg) compared with subjects carrying AA genotypes (123.7 +/- 19.5; 81.6 +/- 11.8; 95.6 +/- 13.6 mmHg) (p = 0.01; p = 0.02; p = 0.01, respectively), even after adjustment for covariates. However, in analysis stratified by ethnicity, this finding was not found and there is no evidence that the polymorphism influences BP. Conclusion: Our findings indicate that association studies involving this MYLIP variant can present distinct results according to the studied population. In this moment, further studies are needed to reaffirm if the MYLIP p.N342S polymorphism is functional or not, and to identify other functional markers within this gene.

Using Prior Information from the Medical Literature in GWAS of Oral Cancer Identifies Novel Susceptibility Variant on Chromosome 4-the AdAPT Method

Background: Genome-wide association studies (GWAS) require large sample sizes to obtain adequate statistical power, but it may be possible to increase the power by incorporating complementary data. In this study we investigated the feasibility of automatically retrieving information from the medical literature and leveraging this information in GWAS. Methods: We developed a method that searches through PubMed abstracts for pre-assigned keywords and key concepts, and uses this information to assign prior probabilities of association for each single nucleotide polymorphism (SNP) with the phenotype of interest - the Adjusting Association Priors with Text (AdAPT) method. Association results from a GWAS can subsequently be ranked in the context of these priors using the Bayes False Discovery Probability (BFDP) framework. We initially tested AdAPT by comparing rankings of known susceptibility alleles in a previous lung cancer GWAS, and subsequently applied it in a two-phase GWAS of oral cancer. Results: Known lung cancer susceptibility SNPs were consistently ranked higher by AdAPT BFDPs than by p-values. In the oral cancer GWAS, we sought to replicate the top five SNPs as ranked by AdAPT BFDPs, of which rs991316, located in the ADH gene region of 4q23, displayed a statistically significant association with oral cancer risk in the replication phase (per-rare-allele log additive p-value [p(trend)] = 2.5 x 10(-3)). The combined OR for having one additional rare allele was 0.83 (95% CI: 0.76-0.90), and this association was independent of previously identified susceptibility SNPs that are associated with overall UADT cancer in this gene region. We also investigated if rs991316 was associated with other cancers of the upper aerodigestive tract (UADT), but no additional association signal was found. Conclusion: This study highlights the potential utility of systematically incorporating prior knowledge from the medical literature in genome-wide analyses using the AdAPT methodology. AdAPT is available online (url: http://services.gate.ac.uk/lld/gwas/service/config).

Brazilian urban population genetic structure reveals a high degree of admixture

Advances in genotyping technologies have contributed to a better understanding of human population genetic structure and improved the analysis of association studies. To analyze patterns of human genetic variation in Brazil, we used SNP data from 1129 individuals - 138 from the urban population of Sao Paulo, Brazil, and 991 from 11 populations of the HapMap Project. Principal components analysis was performed on the SNPs common to these populations, to identify the composition and the number of SNPs needed to capture the genetic variation of them. Both admixture and local ancestry inference were performed in individuals of the Brazilian sample. Individuals from the Brazilian sample fell between Europeans, Mexicans, and Africans. Brazilians are suggested to have the highest internal genetic variation of sampled populations. Our results indicate, as expected, that the Brazilian sample analyzed descend from Amerindians, African, and/or European ancestors, but intermarriage between individuals of different ethnic origin had an important role in generating the broad genetic variation observed in the present-day population. The data support the notion that the Brazilian population, due to its high degree of admixture, can provide a valuable resource for strategies aiming at using admixture as a tool for mapping complex traits in humans. European Journal of Human Genetics (2012) 20, 111-116; doi:10.1038/ejhg.2011.144; published online 24 August 2011

ALDH1A2 (RALDH2) genetic variation in human congenital heart disease

Abstract Background Signaling by the vitamin A-derived morphogen retinoic acid (RA) is required at multiple steps of cardiac development. Since conversion of retinaldehyde to RA by retinaldehyde dehydrogenase type II (ALDH1A2, a.k.a RALDH2) is critical for cardiac development, we screened patients with congenital heart disease (CHDs) for genetic variation at the ALDH1A2 locus. Methods One-hundred and thirty-three CHD patients were screened for genetic variation at the ALDH1A2 locus through bi-directional sequencing. In addition, six SNPs (rs2704188, rs1441815, rs3784259, rs1530293, rs1899430) at the same locus were studied using a TDT-based association approach in 101 CHD trios. Observed mutations were modeled through molecular mechanics (MM) simulations using the AMBER 9 package, Sander and Pmemd programs. Sequence conservation of observed mutations was evaluated through phylogenetic tree construction from ungapped alignments containing ALDH8 s, ALDH1Ls, ALDH1 s and ALDH2 s. Trees were generated by the Neighbor Joining method. Variations potentially affecting splicing mechanisms were cloned and functional assays were designed to test splicing alterations using the pSPL3 splicing assay. Results We describe in Tetralogy of Fallot (TOF) the mutations Ala151Ser and Ile157Thr that change non-polar to polar residues at exon 4. Exon 4 encodes part of the highly-conserved tetramerization domain, a structural motif required for ALDH oligomerization. Molecular mechanics simulation studies of the two mutations indicate that they hinder tetramerization. We determined that the SNP rs16939660, previously associated with spina bifida and observed in patients with TOF, does not affect splicing. Moreover, association studies performed with classical models and with the transmission disequilibrium test (TDT) design using single marker genotype, or haplotype information do not show differences between cases and controls. Conclusion In summary, our screen indicates that ALDH1A2 genetic variation is present in TOF patients, suggesting a possible causal role for this gene in rare cases of human CHD, but does not support the hypothesis that variation at the ALDH1A2 locus is a significant modifier of the risk for CHD in humans.

Relationship between Vitamin D Receptor gene polymorphisms and the components of metabolic syndrome

Abstract Background The Vitamin D Receptor gene (VDR) is expressed in many tissues and modulates the expression of several other genes. The purpose of this study was to investigate the association between metabolic syndrome (MetSyn) with the presence of VDR 2228570 C > T and VDR 1544410 A > G polymorphisms in Brazilian adults. Methods Two hundred forty three (243) individuals were included in a cross-sectional study. MetSyn was classified using the criteria proposed by National Cholesterol Educational Program - Adult Treatment Panel III. Insulin resistance and β cell secretion were estimated by the mathematical models of HOMA IR and β, respectively. The VDR 2228570 C > T and VDR 1544410 A > G polymorphisms were detected by enzymatic digestion and confirmed by allele specific PCR or amplification of refractory mutation. Results Individuals with MetSyn and heterozygosis for VDR 2228570 C > T have higher concentrations of iPTH and HOMA β than those without this polymorphism, and subjects with recessive homozygosis for the same polymorphisms presented higher insulin resistance than those with the heterozygous genotype. There is no association among VDR 1544410 A > G and components of MetSyn, HOMA IR and β, serum vitamin D (25(OH)D3) and intact parathormone (iPTH) levels in patients with MetSyn. A significant lower concentration of 25(OH)D3 was observed only in individuals without MetSyn in the VDR 1544410 A > G genotype. Additionally, individuals without MetSyn and heterozygosis for VDR 2228570 C > T presented higher concentration of triglycerides and lower HDL than those without this polymorphism. Conclusions Using two common VDR polymorphism data suggests they may influence insulin secretion, insulin resistance an serum HDL-cholesterol in our highly heterogeneous population. Whether VDR polymorphism may influence the severity of MetSyn component disorder, warrants examination in larger cohorts used for genome-wide association studies.

The humankind genome: from genetic diversity to the origin of human diseases

Genome-wide association studies have failed to establish common variant risk for the majority of common human diseases. The underlying reasons for this failure are explained by recent studies of resequencing and comparison of over 1200 human genomes and 10 000 exomes, together with the delineation of DNA methylation patterns (epigenome) and full characterization of coding and noncoding RNAs (transcriptome) being transcribed. These studies have provided the most comprehensive catalogues of functional elements and genetic variants that are now available for global integrative analysis and experimental validation in prospective cohort studies. With these datasets, researchers will have unparalleled opportunities for the alignment, mining, and testing of hypotheses for the roles of specific genetic variants, including copy number variations, single nucleotide polymorphisms, and indels as the cause of specific phenotypes and diseases. Through the use of next-generation sequencing technologies for genotyping and standardized ontological annotation to systematically analyze the effects of genomic variation on humans and model organism phenotypes, we will be able to find candidate genes and new clues for disease’s etiology and treatment. This article describes essential concepts in genetics and genomic technologies as well as the emerging computational framework to comprehensively search websites and platforms available for the analysis and interpretation of genomic data.

Revealing archaeal diversity patterns and methane fluxes in Admiralty Bay, King George Island, and their association to Brazilian Antarctic Station activities

The study of Antarctic archaeal communities adds information on the biogeography of this group and helps understanding the dynamics of biogenic methane production in such extreme habitats. Molecular methods were combined to methane flux determinations in Martel Inlet, Admiralty Bay, to assess archaeal diversity, to obtain information about contribution of the area to atmospheric methane budget and to detect possible interferences of the Antarctic Brazilian Station Comandante Ferraz (EACF) wastewater discharge on local archaeal communities and methane emissions. Methane fluxes in Martel Inlet ranged from 3.2 to 117.9 mu mol CH(4) m(-2) d(-1), with an average of 51.3 +/- 8.5 mu mol CH(4) m(-2) d(-1) and a median of 57.6 mu mol CH(4) m(-2)d(-1). However, three negative fluxes averaging -11.3 mu mol CH(4) m(-2) d(-1) were detected in MacKellar Inlet, indicating that Admiralty Bay can be either a source or sink of atmospheric methane. Denaturing gradient gel electrophoresis (DGGE) showed that archaeal communities at EACF varied with depth and formed a group separated from the reference sites. Granulometric analysis indicated that differences observed may be mostly related to sediment type. However, an influence of wastewater input could not be discarded, since higher methane fluxes were found at CF site. suggesting stimulation of local methanogenesis. DGGE profile of the wastewater sample grouped separated from all other samples, suggesting that methanogenesis stimulation may be due to changes in environmental conditions rather than to the input of allochtonous species from the wastewater. 16S ribosomal DNA clone libraries analysis showed that all wastewater sequences were related to known methanogenic groups belonging to the hydrogenotrophic genera Methanobacterium and Methanobrevibacter and the aceticlastic genus Methanosaeta. EACF and Botany Point sediment clone libraries retrieved only groups of uncultivated Archaea, with predominance of Crenarchaeota representatives (MCG, MG1, MBG-B, MBG-C and MHVG groups). Euryarchaeota sequences found were mostly related to the LDS and RC-V groups, but MBG-D and DHVE-5 were also present. No representatives of cultivated methanogenic groups were found, but coverage estimates suggest that a higher number of clones would have to be analyzed in order to cover the greater archaeal diversity of Martel Inlet sediment. Nevertheless, the analysis of the libraries revealed groups not commonly found by other authors in Antarctic habitats and also indicated the presence of groups of uncultivated archaea previously associated to methane rich environments or to the methane cycle. (C) 2010 Elsevier Ltd. All rights reserved.

Quantitative proteomic analysis and functional studies reveal that nucleophosmin is involved in cell death in glioblastoma cell line transfected with siRNA

Previously, we reported that nucleophosmin (NPM) was increased in glioblastoma multiforme (GBM). NPM is a phosphoprotein related to apoptosis, ribosome biogenesis, mitosis, and DNA repair, but details about its function remain unclear. We treated U87MG and A172 cells with small interference RNA (siRNA) and obtained a reduction of 80% in NPM1 expression. Knockdown at the protein level was evident after the 4th day and was maintained until the 7th day of transfection that was investigated by quantitative proteomic analysis using isobaric tags. The comparison of proteomic analysis of NPM1-siRNA against controls allowed the identification of 14 proteins, two proteins showed increase and 12 presented a reduction of expression levels. Gene ontology assigned most of the hypoexpressed proteins to apoptosis regulation, including GRP78. NPM1 silencing did not impair cell proliferation until the 7th day after transfection, but sensitized U87MG cells to temozolomide (TMZ), culminating with an increase in cell death and provoking at a later period a reduction of colony formation. In a large data set of GBM patients, both GRP78 and NPM1 genes were upregulated and presented a tendency to shorter overall survival time. In conclusion, NPM proved to participate in the apoptotic process, sensitizing TMZ-treated U87MG and A172 cells to cell death, and in association with upregulation of GRP78 may be helpful as a predictive factor of poor prognosis in GBM patients.

The association between low alcohol use and traffic risk behaviors among Brazilian college students

Although there are a large number of studies focused on binge drinking and traffic risk behaviors (TRB), little is known regarding low levels of alcohol consumption and its association to TRB. The aim of this cross-sectional study is to examine the association of low to moderate alcohol intake pattern and TRB in college students in Brazil. 7037 students from a National representative sample were selected under rigorous inclusion criteria. All study participants voluntarily fulfilled a structured, anonymous, and self-questionnaire regarding alcohol and drug use, social-demographic data, and TRB. Alcohol was assessed according to the average number of alcoholic units consumed on standard occasions over the past 12 months. The associations between alcohol intake and TRB were summarized with odds ratio and their confidence interval obtained from logistic regression. Compared with abstainers students who consumed only one alcohol unit had the risk of being a passenger in a car driven by a drunk driver increased by almost four times, students who reported using five or more units were increased by almost five times the risk of being involved in a car crash. Compared with students who consumed one alcohol unit, the risk of driving under the influence of alcohol increased four times in students using three alcohol units. Age group, use of illicit drugs, employment status, gender, and marital status significantly influenced occurrence of TRB among college students. Our study highlights the potential detrimental effects of low and moderate pattern of alcohol consumption and its relation to riding with an intoxicated driver and other TRB. These data suggest that targeted interventions should be implemented in order to prevent negative consequences due to alcohol use in this population. (C) 2012 Elsevier Inc. All rights reserved,

Pure red cell aplasia and primary antiphospholipid syndrome: a unique association

Pure red cell aplasia (PRCA) is a disease with important relationships to autoimmune mechanisms. Although some autoimmune diseases, such as rheumatoid arthritis and systemic lupus erythematosus, have been associated with PRCA, until this point no studies have described the association between PRCA and primary antiphospholipid syndrome (APS). This is the first case report of PRCA associated with primary APS in a 39-year-old man with acute heart failure secondary to an anaemic condition that was diagnosed as pure red cell aplasia. The patient was later diagnosed with retinal artery and vein thromboses and bilateral deep venous thromboses of the femoral and popliteal veins. The most common causes of PRCA and other thrombophilias were ruled out in this investigation through complementary tests. This association with APS adds a new possibility to the study of PRCA pathophysiology.

Motor Neuron Disease and Acquired Axonal Neuropathy Association in HIV Infection: Case Report and Update

Background: A possible viral etiology has been documented in the genesis of motor neuron disorders and acquired peripheral neuropathies, mainly due to the vulnerability of peripheral nerves and the anterior horn to certain viruses. In recent years, several reports show association of HIV infection with Amyotrophic Lateral Sclerosis Syndrome, Motor Neuron Diseases and peripheral neuropathies. Objective: To report a case of an association between Motor Neuron Disease and Acquired Axonal neuropathy in HIV infection, and describe the findings of neurological examination, cerebrospinal fluid, neuroimaging and electrophysiology. Methods: The patient underwent neurological examination. General medical examinations were performed, including, specific neuromuscular tests, analysis of cerebrospinal fluid, muscle biopsy and imaging studies. Results and Discussion: The initial clinical presentation of our case was marked by cramps and fasciculations with posterior distal paresis and atrophy in the left arm. We found electromyography tracings with deficits in the anterior horn of the spinal cord and peripheral nerves. Dysphagia and release of primitive reflexes were also identified. At the same time, the patient was informed to be HIV positive with high viral load. He received antiretroviral therapy, with load control but with no clinical remission. Conclusion: Motor Neuron disorders and peripheral neuropathy may occur in association with HIV infection. However, a causal relationship remains uncertain. It is noteworthy that the antiretroviral regimen may be implicated in some cases.