Accounting for Multiple Comparisons in a Genome-Wide Association Study (GWAS)

Background As we enter an era when testing millions of SNPs in a single gene association study will become the standard, consideration of multiple comparisons is an essential part of determining statistical significance. Bonferroni adjustments can be made but are conservative due to the preponderance of linkage disequilibrium (LD) between genetic markers, and permutation testing is not always a viable option. Three major classes of corrections have been proposed to correct the dependent nature of genetic data in Bonferroni adjustments: permutation testing and related alternatives, principal components analysis (PCA), and analysis of blocks of LD across the genome. We consider seven implementations of these commonly used methods using data from 1514 European American participants genotyped for 700,078 SNPs in a GWAS for AIDS. Results A Bonferroni correction using the number of LD blocks found by the three algorithms implemented by Haploview resulted in an insufficiently conservative threshold, corresponding to a genome-wide significance level of α = 0.15 - 0.20. We observed a moderate increase in power when using PRESTO, SLIDE, and simpleℳ when compared with traditional Bonferroni methods for population data genotyped on the Affymetrix 6.0 platform in European Americans (α = 0.05 thresholds between 1 × 10-7 and 7 × 10-8). Conclusions Correcting for the number of LD blocks resulted in an anti-conservative Bonferroni adjustment. SLIDE and simpleℳ are particularly useful when using a statistical test not handled in optimized permutation testing packages, and genome-wide corrected p-values using SLIDE, are much easier to interpret for consumers of GWAS studies

Pyrolysis of medium-density fiberboard: optimized search for kinetics scheme and parameters via a genetic algorithm driven by Kissinger's method

The pyrolysis kinetics of charring materials plays an important role in understanding material combustions especially for construction materials with complex degradation chemistry. Thermogravimetric analysis (TGA) is frequently used to study the heterogeneous kinetics of solid fuels; however, there is no agreed method to determine the pyrolysis scheme and kinetic parameters for charring polymers with multiple components and competing reaction pathways. This study develops a new technique to estimate the possible numbers of species and sub-reactions in pyrolysis by analyzing the second derivatives of thermogravimetry (DDTG) curves. The pyrolysis of a medium-density fiberboard (MDF) in nitrogen is studied in detail, and the DDTG curves are used to locate the temperature of the peak mass-loss rate for each sub-reaction. Then, on the basis of the TG data under multiple heating rates, Kissinger’s method is used to quickly find the possible range of values of the kinetic parameters (<i>A</i> and <i>E</i>). These ranges are used to accelerate the optimization of the inverse problem using a genetic algorithm (GA) for the kinetic and stoichiometric parameters. The proposed method and kinetic scheme found are shown to match the experimental data and are able to predict accurately results at different heating rates better than Kissinger’s method. Moreover, the search method (K–K method) is highly efficient, faster than the regular GA search alone. Modeling results show that, as the TG data available increase, the interdependence among kinetic parameters becomes weak and the accuracy of the first-order model declines. Furthermore, conducting TG experiment under multiple heating rates is found to be crucial in obtaining good kinetic parameters

The Principal Genetic Determinants for Nasopharyngeal Carcinoma in China Involve the HLA Class I Antigen Recognition Groove

Nasopharyngeal carcinoma (NPC) is an epithelial malignancy facilitated by Epstein-Barr Virus infection. Here we resolve the major genetic influences for NPC incidence using a genome-wide association study (GWAS), independent cohort replication, and high-resolution molecular HLA class I gene typing including 4,055 study participants from the Guangxi Zhuang Autonomous Region and Guangdong province of southern China. We detect and replicate strong association signals involving SNPs, HLA alleles, and amino acid (aa) variants across the major histocompatibility complex-HLA-A, HLA –B, and HLA -C class I genes (PHLA-A-aa-site-62 = 7.4×10−29; P HLA-B-aa-site-116 = 6.5×10−19; P HLA-C-aa-site-156 = 6.8×10−8 respectively). Over 250 NPC-HLA associated variants within HLA were analyzed in concert to resolve separate and largely independent HLA-A, -B, and -C gene influences. Multivariate logistical regression analysis collapsed significant associations in adjacent genes spanning 500 kb (OR2H1, GABBR1, HLA-F, and HCG9) as proxies for peptide binding motifs carried by HLA- A*11:01. A similar analysis resolved an independent association signal driven by HLA-B*13:01, B*38:02, and B*55:02 alleles together. NPC resistance alleles carrying the strongly associated amino acid variants implicate specific class I peptide recognition motifs in HLA-A and -B peptide binding groove as conferring strong genetic influence on the development of NPC in China

Genetic Factors Leading to Chronic Epstein–Barr Virus Infection and Nasopharyngeal Carcinoma in South East China: Study Design, Methods and Feasibility

Nasopharyngeal carcinoma (NPC) is a complex disease caused by a combination of Epstein-Barr virus chronic infection, the environment and host genes in a multi-step process of carcinogenesis. The identity of genetic factors involved in the development of chronic Epstein-Barr virus infection and NPC remains elusive, however. Here, we describe a two-phase, population-based, case-control study of Han Chinese from Guangxi province, where the NPC incidence rate rises to a high of 25-50 per 100,000 individuals. Phase I, powered to detect single gene associations, enrolled 984 subjects to determine feasibility, to develop infrastructure and logistics and to determine error rates in sample handling. A microsatellite screen of Phase I study participants, genotyped for 319 alleles from 34 microsatellites spanning an 18-megabase region of chromosome 4 (4p15.1-q12), previously implicated by a linkage analysis of familial NPC, found 14 alleles marginally associated with developing NPC or chronic immunoglobulin A production (p = 0.001-0.03). These associations lost significance after applying a correction for multiple tests. Although the present results await confirmation, the Phase II study population has tripled patient enrolment and has included environmental covariates, offering the potential to validate this and other genomic regions that influence the onset of NPC

Experimental review of oxygen content at mixing layer in cone calorimeter

This work aims to elucidate whether the hypothesis of zero oxygen at the mixture layer when flame takes place is assumable for every kind of material. For that purpose, we investigated the oxygen concentration there by cone calorimeter tests. A modified holder was developed in order to collect oxygen in this mixture layer. In addition, thermogravimetric tests were carried out so as to relate the possible effects of the presence of oxygen in the atmosphere where the pyrolysis process takes place, since the cone calorimeter does not allow to control the oxygen level of the atmosphere during the experiment. The reaction rates and per cent of residue in the cone calorimetric tests were measured and compared with the results from thermogravimetric tests. Six products were analysed which can be classified in three main groups: lignocellulosic, thermoplastic polymers and thermoset polymers. Cone calorimetric results showed that for some of the materials analysed (PET, Nylon and PUR foam) the oxygen level at mixture layer decreased until values close to zero. The comparison of reaction rates between cone calorimetric and thermogravimetric tests revealed the char layer created in cone calorimetric tests over the exposed face for brushed fir, Nylon and PET established an important heat barrier that modifies the thermal behaviour of these materials.Authors would like to thank to the Spanish Ministry of Economy and Competitiveness for the PYRODESIGN Project grant, Ref.: BIA2012-37890, financed jointly by FEDER funds

Genetic Variants in Nuclear-Encoded Mitochondrial Genes Influence AIDS Progression

Background: The human mitochondrial genome includes only 13 coding genes while nuclear-encoded genes account for 99% of proteins responsible for mitochondrial morphology, redox regulation, and energetics. Mitochondrial pathogenesis occurs in HIV patients and genetically, mitochondrial DNA haplogroups with presumed functional differences have been associated with differential AIDS progression. Methodology/Principal Findings: Here we explore whether single nucleotide polymorphisms (SNPs) within 904 of the estimated 1,500 genes that specify nuclear-encoded mitochondrial proteins (NEMPs) influence AIDS progression among HIV-1 infected patients. We examined NEMPs for association with the rate of AIDS progression using genotypes generated by an Affymetrix 6.0 genotyping array of 1,455 European American patients from five US AIDS cohorts. Successfully genotyped SNPs gave 50% or better haplotype coverage for 679 of known NEMP genes. With a Bonferroni adjustment for the number of genes and tests examined, multiple SNPs within two NEMP genes showed significant association with AIDS progression: acyl-CoA synthetase medium-chain family member 4 (ACSM4) on chromosome 12 and peroxisomal D3,D2-enoyl- CoA isomerase (PECI) on chromosome 6. Conclusions: Our previous studies on mitochondrial DNA showed that European haplogroups with presumed functional differences were associated with AIDS progression and HAART mediated adverse events. The modest influences of nuclearencoded mitochondrial genes found in the current study add support to the idea that mitochondrial function plays a role in AIDS pathogenesis

Analysis of genetic variation in Ashkenazi Jews by high density SNP genotyping.

BACKGROUND: Genetic isolates such as the Ashkenazi Jews (AJ) potentially offer advantages in mapping novel loci in whole genome disease association studies. To analyze patterns of genetic variation in AJ, genotypes of 101 healthy individuals were determined using the Affymetrix EAv3 500 K SNP array and compared to 60 CEPH-derived HapMap (CEU) individuals. 435,632 SNPs overlapped and met annotation criteria in the two groups. RESULTS: A small but significant global difference in allele frequencies between AJ and CEU was demonstrated by a mean FST of 0.009 (P < 0.001); large regions that differed were found on chromosomes 2 and 6. Haplotype blocks inferred from pairwise linkage disequilibrium (LD) statistics (Haploview) as well as by expectation-maximization haplotype phase inference (HAP) showed a greater number of haplotype blocks in AJ compared to CEU by Haploview (50,397 vs. 44,169) or by HAP (59,269 vs. 54,457). Average haplotype blocks were smaller in AJ compared to CEU (e.g., 36.8 kb vs. 40.5 kb HAP). Analysis of global patterns of local LD decay for closely-spaced SNPs in CEU demonstrated more LD, while for SNPs further apart, LD was slightly greater in the AJ. A likelihood ratio approach showed that runs of homozygous SNPs were approximately 20% longer in AJ. A principal components analysis was sufficient to completely resolve the CEU from the AJ. CONCLUSION: LD in the AJ versus was lower than expected by some measures and higher by others. Any putative advantage in whole genome association mapping using the AJ population will be highly dependent on regional LD structure