BiForce Toolbox: powerful high-throughput computational analysis of gene-gene interactions in genome-wide association studies

Genome-wide association studies (GWAS) have discovered many loci associated with common disease and quantitative traits. However, most GWAS have not studied the gene–gene interactions (epistasis) that could be important in complex trait genetics. A major challenge in analysing epistasis in GWAS is the enormous computational demands of analysing billions of SNP combinations. Several methods have been developed recently to address this, some using computers equipped with particular graphical processing units, most restricted to binary disease traits and all poorly suited to general usage on the most widely used operating systems. We have developed the BiForce Toolbox to address the demand for high-throughput analysis of pairwise epistasis in GWAS of quantitative and disease traits across all commonly used computer systems. BiForce Toolbox is a stand-alone Java program that integrates bitwise computing with multithreaded parallelization and thus allows rapid full pairwise genome scans via a graphical user interface or the command line. Furthermore, BiForce Toolbox incorporates additional tests of interactions involving SNPs with significant marginal effects, potentially increasing the power of detection of epistasis. BiForce Toolbox is easy to use and has been applied in multiple studies of epistasis in large GWAS data sets, identifying interesting interaction signals and pathways

Association Analysis in African Americans of European-Derived Type 2 Diabetes Single Nucleotide Polymorphisms From Whole-Genome Association Studies

OBJECTIVE— Several whole-genome association studies have reported identification of type 2 diabetes susceptibility genes in various European-derived study populations. Little investigation of these loci has been reported in other ethnic groups, specifically African Americans. Striking differences exist between these populations, suggesting they may not share identical genetic risk factors. Our objective was to examine the influence of type 2 diabetes genes identified in whole-genome association studies in a large African American case-control population

Increased burden of familial-associated early-onset cancer risk among minority americans compared to non-latino whites

Background: The role of race/ethnicity in genetic predisposition of early-onset cancers can be estimated by comparing family-based cancer concordance rates among ethnic groups. Methods: We used linked California health registries to evaluate the relative cancer risks for first-degree relatives of patients diagnosed between ages 0 and 26, and the relative risks of developing distinct second primary malignancies (SPMs). From 1989 to 2015, we identified 29,631 cancer patients and 62,863 healthy family members. We calculated the standardized incident ratios (SIRs) of early-onset primary cancers diagnosed in proband siblings and mothers, as well as SPMs detected among early-onset patients. Analyses were stratified by self-identified race/ethnicity. Results: Given probands with cancer, there were increased relative risks of any cancer for siblings and mothers (SIR = 3.32; 95% confidence interval [CI]: 2.85–3.85) and of SPMs (SIR = 7.27; 95% CI: 6.56–8.03). Given a proband with solid cancer, both Latinos (SIR = 4.98; 95% CI: 3.82–6.39) and non-Latino Blacks (SIR = 7.35; 95% CI: 3.36–13.95) exhibited significantly higher relative risk of any cancer in siblings and mothers when compared to non-Latino White subjects (SIR = 3.02; 95% CI: 2.12–4.16). For hematologic cancers, higher familial risk was evident for Asian/Pacific Islanders (SIR = 7.56; 95% CI: 3.26–14.90) compared to non-Latino whites (SIR = 2.69; 95% CI: 1.62–4.20). Conclusions: The data support a need for increased attention to the genetics of early-onset cancer predisposition and environmental factors in race/ethnic minority families in the United States.</p

A complete classification of epistatic two-locus models

Background: The study of epistasis is of great importance in statistical genetics in fields such as linkage and association analysis and QTL mapping. In an effort to classify the types of epistasis in the case of two biallelic loci Li and Reich listed and described all models in the simplest case of 0/ 1 penetrance values. However, they left open the problem of finding a classification of two-locus models with continuous penetrance values. Results: We provide a complete classification of biallelic two-locus models. In addition to solving the classification problem for dichotomous trait disease models, our results apply to any instance where real numbers are assigned to genotypes, and provide a complete framework for studying epistasis in QTL data. Our approach is geometric and we show that there are 387 distinct types of two-locus models, which can be reduced to 69 when symmetry between loci and alleles is accounted for. The model types are defined by 86 circuits, which are linear combinations of genotype values, each of which measures a fundamental unit of interaction. Conclusion: The circuits provide information on epistasis beyond that contained in the additive × additive, additive × dominance, and dominance × dominance interaction terms. We discuss th

Probing the diabetes and colorectal cancer relationship using gene – environment interaction analyses

BackgroundDiabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis.MethodsWe used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test).ResultsBased on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value(3-d.f.): 5.46 x 10(-11)) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value(2-d.f.): 7.84 x 10(-09)).DiscussionThese results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship

Genetic Analysis Workshop 16: Strategies for genome-wide association study analyses

National Institutes of Health (AR44422, N01-AR-7-2232); Genome Canada and Associations AFP; Polyarctique-Groupe Taitbout; Rhumatisme et Travail; Arthritis Research Campaign; National Institute of General Medical Sciences (R01 GM31575

Increased Oxidative Burden Associated with Traffic Component of Ambient Particulate Matter at Roadside and Urban Background Schools Sites in London

As the incidence of respiratory and allergic symptoms has been reported to be increased in children attending schools in close proximity to busy roads, it was hypothesised that PM from roadside schools would display enhanced oxidative potential (OP). Two consecutive one-week air quality monitoring campaigns were conducted at seven school sampling sites, reflecting roadside and urban background in London. Chemical characteristics of size fractionated particulate matter (PM) samples were related to the capacity to drive biological oxidation reactions in a synthetic respiratory tract lining fluid. Contrary to hypothesised contrasts in particulate OP between school site types, no robust size-fractionated differences in OP were identified due high temporal variability in concentrations of PM components over the one-week sampling campaigns. For OP assessed both by ascorbate (OPAA m−3) and glutathione (OPGSH m−3) depletion, the highest OP per cubic metre of air was in the largest size fraction, PM1.9–10.2. However, when expressed per unit mass of particles OPAA µg−1 showed no significant dependence upon particle size, while OPGSH µg−1 had a tendency to increase with increasing particle size, paralleling increased concentrations of Fe, Ba and Cu. The two OP metrics were not significantly correlated with one another, suggesting that the glutathione and ascorbate depletion assays respond to different components of the particles. Ascorbate depletion per unit mass did not show the same dependence as for GSH and it is possible that other trace metals (Zn, Ni, V) or organic components which are enriched in the finer particle fractions, or the greater surface area of smaller particles, counter-balance the redox activity of Fe, Ba and Cu in the coarse particles. Further work with longer-term sampling and a larger suite of analytes is advised in order to better elucidate the determinants of oxidative potential, and to fuller explore the contrasts between site types.\ud \u

Comparison of Oxidative Properties, Light Absorbance, and Total and Elemental Mass Concentration of Ambient PM(2.5) Collected at 20 European Sites

OBJECTIVE: It has been proposed that the redox activity of particles may represent a major determinant of their toxicity. We measured the in vitro ability of ambient fine particles [particulate matter with aerodynamic diameters ≤2.5 μm (PM(2.5))] to form hydroxyl radicals ((•)OH) in an oxidant environment, as well as to deplete physiologic antioxidants (ascorbic acid, glutathione) in the naturally reducing environment of the respiratory tract lining fluid (RTLF). The objective was to examine how these toxicologically relevant measures were related to other PM characteristics, such as total and elemental mass concentration and light absorbance. DESIGN: Gravimetric PM(2.5) samples (n = 716) collected over 1 year from 20 centers participating in the European Community Respiratory Health Survey were available. Light absorbance of these filters was measured with reflectometry. PM suspensions were recovered from filters by vortexing and sonication before dilution to a standard concentration. The oxidative activity of these particle suspensions was then assessed by measuring their ability to generate (•)OH in the presence of hydrogen peroxide, using electron spin resonance and 5,5-dimethyl-1-pyrroline-N-oxide as spin trap, or by establishing their capacity to deplete antioxidants from a synthetic model of the RTLF. RESULTS AND CONCLUSION: PM oxidative activity varied significantly among European sampling sites. Correlations between oxidative activity and all other characteristics of PM were low, both within centers (temporal correlation) and across communities (annual mean). Thus, no single surrogate measure of PM redox activity could be identified. Because these novel measures are suggested to reflect crucial biologic mechanisms of PM, their use may be pertinent in epidemiologic studies. Therefore, it is important to define the appropriate methods to determine oxidative activity of PM