FastEpistasis: a high performance computing solution for quantitative trait epistasis.

Motivation: Genome-wide association studies have become widely used tools to study effects of genetic variants on complex diseases. While it is of great interest to extend existing analysis methods by considering interaction effects between pairs of loci, the large number of possible tests presents a significant computational challenge. The number of computations is further multiplied in the study of gene expression quantitative trait mapping, in which tests are performed for thousands of gene phenotypes simultaneously. Results: We present FastEpistasis, an efficient parallel solution extending the PLINK epistasis module, designed to test for epistasis effects when analyzing continuous phenotypes. Our results show that the algorithm scales with the number of processors and offers a reduction in computation time when several phenotypes are analyzed simultaneously. FastEpistasis is capable of testing the association of a continuous trait with all single nucleotide polymorphism ( SNP) pairs from 500 000 SNPs, totaling 125 billion tests, in a population of 5000 individuals in 29, 4 or 0.5 days using 8, 64 or 512 processors.

Crohn's disease-associated polymorphism within the PTPN2 gene affects muramyl-dipeptide-induced cytokine secretion and autophagy.

BACKGROUND: The single nucleotide polymorphism (SNP) rs2542151 within the gene locus region encoding protein tyrosine phosphatase non-receptor type 2 (PTPN2) has been associated with Crohn's disease (CD), ulcerative colitis (UC), type-I diabetes, and rheumatoid arthritis. We have previously shown that PTPN2 regulates mitogen-activated protein kinase (MAPK) signaling and cytokine secretion in human THP-1 monocytes and intestinal epithelial cells (IEC). Here, we studied whether intronic PTPN2 SNP rs1893217 regulates immune responses to the nucleotide-oligomerization domain 2 (NOD2) ligand, muramyl-dipeptide (MDP). MATERIALS AND METHODS: Genomic DNA samples from 343 CD and 663 non-IBD control patients (male and female) from a combined German, Swiss, and Polish cohort were genotyped for the presence of the PTPN2 SNPs, rs2542151, and rs1893217. PTPN2-variant rs1893217 was introduced into T(84) IEC or THP-1 cells using a lentiviral vector. RESULTS: We identified a novel association between the genetic variant, rs1893217, located in intron 7 of the PTPN2 gene and CD. Human THP-1 monocytes carrying this variant revealed increased MAPK activation as well as elevated mRNA expression of T-bet transcription factor and secretion of interferon-γ in response to the bacterial wall component, MDP. In contrast, secretion of interleukin-8 and tumor necrosis factor were reduced. In both, T(84) IEC and THP-1 monocytes, autophagosome formation was impaired. CONCLUSIONS: We identified a novel CD-associated PTPN2 variant that modulates innate immune responses to bacterial antigens. These findings not only provide key insights into the effects of a functional mutation on a clinically relevant gene, but also reveal how such a mutation could contribute to the onset of disease.

Improving the prediction of the functional impact of cancer mutations by baseline tolerance transformation

High-throughput prioritization of cancer-causing mutations (drivers) is a key challenge of cancer genome projects, due to the number of somatic variants detected in tumors. One important step in this task is to assess the functional impact of tumor somatic mutations. A number of computational methods have been employed for that purpose, although most were originally developed to distinguish disease-related nonsynonymous single nucleotide variants (nsSNVs) from polymorphisms. Our new method, transformed Functional Impact score for Cancer (transFIC), improves the assessment of the functional impact of tumor nsSNVs by taking into account the baseline tolerance of genes to functional variants.

Normal variation in fronto-occipital circuitry and cerebellar structure with an autism-associated polymorphism of CNTNAP2.

Recent genetic studies have implicated a number of candidate genes in the pathogenesis of Autism Spectrum Disorder (ASD). Polymorphisms of CNTNAP2 (contactin-associated like protein-2), a member of the neurexin family, have already been implicated as a susceptibility gene for autism by at least 3 separate studies. We investigated variation in white and grey matter morphology using structural MRI and diffusion tensor imaging. We compared volumetric differences in white and grey matter and fractional anisotropy values in control subjects characterised by genotype at rs7794745, a single nucleotide polymorphism in CNTNAP2. Homozygotes for the risk allele showed significant reductions in grey and white matter volume and fractional anisotropy in several regions that have already been implicated in ASD, including the cerebellum, fusiform gyrus, occipital and frontal cortices. Male homozygotes for the risk alleles showed greater reductions in grey matter in the right frontal pole and in FA in the right rostral fronto-occipital fasciculus compared to their female counterparts who showed greater reductions in FA of the anterior thalamic radiation. Thus a risk allele for autism results in significant cerebral morphological variation, despite the absence of overt symptoms or behavioural abnormalities. The results are consistent with accumulating evidence of CNTNAP2's function in neuronal development. The finding suggests the possibility that the heterogeneous manifestations of ASD can be aetiologically characterised into distinct subtypes through genetic-morphological analysis.

Impact of donor and recipient IL28B rs12979860 genotypes on hepatitis C virus liver graft reinfection.

BACKGROUND & AIMS: Recent studies have described a major impact of genetic variations near the IL28B gene on the natural course and outcome of antiviral therapy in chronic hepatitis C. We therefore, aimed to explore the impact of donor and recipient genotypes of these polymorphisms on hepatitis C virus (HCV) liver graft reinfection. METHODS: Donor and recipient genotypes of IL28B rs12979860C>T single nucleotide polymorphism were determined in 91 patients with HCV liver graft reinfection, 47 of whom were treated with pegylated interferon-α (PEG-IFN-α) and ribavirin. IL28B genetic polymorphisms were correlated with the natural course and treatment outcome of recurrent hepatitis C. RESULTS: Patients requiring liver transplantation due to end-stage chronic hepatitis C appeared to be selected toward the adverse genotypes rs12979860 CT/TT compared to non-transplanted HCV-infected patients (p=0.046). Patients with the donor genotype rs12979860 CC had higher peak ALT and HCV RNA serum concentrations than those with CT/TT (p=0.04 and 0.06, respectively). No association was observed between ALT/HCV RNA serum concentrations and recipient genotypes (p>0.3). More important, donor IL28B rs12979860 CC vs. CT/TT genotypes were associated with rapid, complete early, and sustained virologic response (RVR, cEVR, SVR) to treatment with PEG-IFN-α and ribavirin (p=0.003, 0.0012, 0.008, respectively), but weaker associations of recipient genotypes with RVR, cEVR, and SVR were observed as well (p=0.0046, 0.115, 0.118, respectively). CONCLUSIONS: We provide evidence for a dominant, but not exclusive impact of the donor rather than the recipient IL28B genetic background on the natural course and treatment outcome of HCV liver graft reinfection.

GuavaH: a compendium of host genomic data in HIV biology and disease.

BACKGROUND: There is an ever-increasing volume of data on host genes that are modulated during HIV infection, influence disease susceptibility or carry genetic variants that impact HIV infection. We created GuavaH (Genomic Utility for Association and Viral Analyses in HIV, http://www.GuavaH.org), a public resource that supports multipurpose analysis of genome-wide genetic variation and gene expression profile across multiple phenotypes relevant to HIV biology. FINDINGS: We included original data from 8 genome and transcriptome studies addressing viral and host responses in and ex vivo. These studies cover phenotypes such as HIV acquisition, plasma viral load, disease progression, viral replication cycle, latency and viral-host genome interaction. This represents genome-wide association data from more than 4,000 individuals, exome sequencing data from 392 individuals, in vivo transcriptome microarray data from 127 patients/conditions, and 60 sets of RNA-seq data. Additionally, GuavaH allows visualization of protein variation in ~8,000 individuals from the general population. The publicly available GuavaH framework supports queries on (i) unique single nucleotide polymorphism across different HIV related phenotypes, (ii) gene structure and variation, (iii) in vivo gene expression in the setting of human infection (CD4+ T cells), and (iv) in vitro gene expression data in models of permissive infection, latency and reactivation. CONCLUSIONS: The complexity of the analysis of host genetic influences on HIV biology and pathogenesis calls for comprehensive motors of research on curated data. The tool developed here allows queries and supports validation of the rapidly growing body of host genomic information pertinent to HIV research.

No interactions between previously associated 2-hour glucose gene variants and physical activity or BMI on 2-hour glucose levels.

Gene-lifestyle interactions have been suggested to contribute to the development of type 2 diabetes. Glucose levels 2 h after a standard 75-g glucose challenge are used to diagnose diabetes and are associated with both genetic and lifestyle factors. However, whether these factors interact to determine 2-h glucose levels is unknown. We meta-analyzed single nucleotide polymorphism (SNP) × BMI and SNP × physical activity (PA) interaction regression models for five SNPs previously associated with 2-h glucose levels from up to 22 studies comprising 54,884 individuals without diabetes. PA levels were dichotomized, with individuals below the first quintile classified as inactive (20%) and the remainder as active (80%). BMI was considered a continuous trait. Inactive individuals had higher 2-h glucose levels than active individuals (β = 0.22 mmol/L [95% CI 0.13-0.31], P = 1.63 × 10(-6)). All SNPs were associated with 2-h glucose (β = 0.06-0.12 mmol/allele, P ≤ 1.53 × 10(-7)), but no significant interactions were found with PA (P > 0.18) or BMI (P ≥ 0.04). In this large study of gene-lifestyle interaction, we observed no interactions between genetic and lifestyle factors, both of which were associated with 2-h glucose. It is perhaps unlikely that top loci from genome-wide association studies will exhibit strong subgroup-specific effects, and may not, therefore, make the best candidates for the study of interactions.

Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index.

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation.

Signature polymorphisms in the internal transcribed spacer region relevant for the differentiation of zoophilic and anthropophilic strains of Trichophyton interdigitale and other species of T. mentagrophytes sensu lato.

Summary Background Dermatophytes are the main cause of superficial mycoses in humans and animals. Molecular research has given useful insights into the phylogeny and taxonomy of the dermatophytes to overcome the difficulties with conventional diagnostics. Objectives The Trichophyton mentagrophytes complex consists of anthropophilic as well as zoophilic species. Although several molecular markers have been developed for the differentiation of strains belonging to T. mentagrophytes sensu lato, correct identification still remains problematic, especially concerning the delineation of anthropophilic and zoophilic strains of T. interdigitale. This differentiation is not academic but is essential for selection of the correct antimycotic therapy to treat infected patients. Methods One hundred and thirty isolates identified by morphological characteristics as T. mentagrophytes sensu lato were investigated using restriction fragment length polymorphism (RFLP) and sequence analysis of the polymerase chain reaction-amplified internal transcribed spacer (ITS) region of the rDNA. Results Species of this complex produced individual RFLP patterns obtained by the restriction enzyme MvaI. Subsequent sequence analysis of the ITS1, 5.8S and ITS2 region of all strains, but of T. interdigitale in particular, revealed single unique polymorphisms in anthropophilic and zoophilic strains. Conclusions Signature polymorphisms were observed to be useful for the differentiation of these strains and epidemiological data showed a host specificity among zoophilic strains of T. interdigitale/Arthroderma vanbreuseghemii compared with A. benhamiae as well as characteristic clinical pictures in humans when caused by zoophilic or anthropophilic strains. The delineation is relevant because it helps in determining the correct treatment and provides clues regarding the source of the infection.

Is the adiposity-associated FTO gene variant related to all-cause mortality independent of adiposity? Meta-analysis of data from 169,551 Caucasian adults.

Previously, a single nucleotide polymorphism (SNP), rs9939609, in the FTO gene showed a much stronger association with all-cause mortality than expected from its association with body mass index (BMI), body fat mass index (FMI) and waist circumference (WC). This finding implies that the SNP has strong pleiotropic effects on adiposity and adiposity-independent pathological pathways that leads to increased mortality. To investigate this further, we conducted a meta-analysis of similar data from 34 longitudinal studies including 169,551 adult Caucasians among whom 27,100 died during follow-up. Linear regression showed that the minor allele of the FTO SNP was associated with greater BMI (n = 169,551; 0.32 kg m(-2) ; 95% CI 0.28-0.32, P < 1 × 10(-32) ), WC (n = 152,631; 0.76 cm; 0.68-0.84, P < 1 × 10(-32) ) and FMI (n = 48,192; 0.17 kg m(-2) ; 0.13-0.22, P = 1.0 × 10(-13) ). Cox proportional hazard regression analyses for mortality showed that the hazards ratio (HR) for the minor allele of the FTO SNPs was 1.02 (1.00-1.04, P = 0.097), but the apparent excess risk was eliminated after adjustment for BMI and WC (HR: 1.00; 0.98-1.03, P = 0.662) and for FMI (HR: 1.00; 0.96-1.04, P = 0.932). In conclusion, this study does not support that the FTO SNP is associated with all-cause mortality independently of the adiposity phenotypes.

The Impact of Catechol-O-Methyltransferase and Dopamine D4 Receptor Genotypes on Neurophysiological Markers of Performance Monitoring

Dynamic adaptations of one"s behavior by means of performance monitoring are a central function of the human executive system, that underlies considerable interindividual variation. Converging evidence from electrophysiological and neuroimaging studies in both animals and humans hints atthe importance ofthe dopaminergic system forthe regulation of performance monitoring. Here, we studied the impact of two polymorphisms affecting dopaminergic functioning in the prefrontal cortex [catechol-O-methyltransferase (COMT) Val108/158Met and dopamine D4 receptor (DRD4) single-nucleotide polymorphism (SNP)-521] on neurophysiological correlates of performance monitoring. We applied a modified version of a standard flanker task with an embedded stop-signal task to tap into the different functions involved, particularly error monitoring, conflict detection and inhibitory processes. Participants homozygous for the DRD4 T allele produced an increased error-related negativity after both choice errors and failed inhibitions compared with C-homozygotes. This was associated with pronounced compensatory behavior reflected in higher post-error slowing. No group differences were seen in the incompatibility N2, suggesting distinct effects of the DRD4 polymorphism on error monitoring processes. Additionally, participants homozygous for the COMTVal allele, with a thereby diminished prefrontal dopaminergic level, revealed increased prefrontal processing related to inhibitory functions, reflected in the enhanced stop-signal-related components N2 and P3a. The results extend previous findings from mainly behavioral and neuroimaging data on the relationship between dopaminergic genes and executive functions and present possible underlying mechanisms for the previously suggested association between these dopaminergic polymorphisms and psychiatric disorders as schizophrenia or attention deficit hyperactivity disorder.

Association of PCK1 with Body Mass Index and Other Metabolic Features in Patients With Psychotropic Treatments.

Weight gain is a major health problem among psychiatric populations. It implicates several receptors and hormones involved in energy balance and metabolism. Phosphoenolpyruvate carboxykinase 1 is a rate-controlling enzyme involved in gluconeogenesis, glyceroneogenesis and cataplerosis and has been related to obesity and diabetes phenotypes in animals and humans. The aim of this study was to investigate the association of phosphoenolpyruvate carboxykinase 1 polymorphisms with metabolic traits in psychiatric patients treated with psychotropic drugs inducing weight gain and in general population samples. One polymorphism (rs11552145G > A) significantly associated with body mass index in the psychiatric discovery sample (n = 478) was replicated in 2 other psychiatric samples (n1 = 168, n2 = 188), with AA-genotype carriers having lower body mass index as compared to G-allele carriers. Stronger associations were found among women younger than 45 years carrying AA-genotype as compared to G-allele carriers (-2.25 kg/m, n = 151, P = 0.009) and in the discovery sample (-2.20 kg/m, n = 423, P = 0.0004). In the discovery sample for which metabolic parameters were available, AA-genotype showed lower waist circumference (-6.86 cm, P = 0.008) and triglycerides levels (-5.58 mg/100 mL, P < 0.002) when compared to G-allele carriers. Finally, waist-to-hip ratio was associated with rs6070157 (proxy of rs11552145, r = 0.99) in a population-based sample (N = 123,865, P = 0.022). Our results suggest an association of rs11552145G > A polymorphism with metabolic-related traits, especially in psychiatric populations and in women younger than 45 years.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study.

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHRadjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHRadjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHRadjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHRadjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHRadjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHRadjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Relationship between obesity and the risk of clinically significant depression: Mendelian randomisation study.

BACKGROUND: Obesity has been shown to be associated with depression and it has been suggested that higher body mass index (BMI) increases the risk of depression and other common mental disorders. However, the causal relationship remains unclear and Mendelian randomisation, a form of instrumental variable analysis, has recently been employed to attempt to resolve this issue. AIMS: To investigate whether higher BMI increases the risk of major depression. METHOD: Two instrumental variable analyses were conducted to test the causal relationship between obesity and major depression in RADIANT, a large case-control study of major depression. We used a single nucleotide polymorphism (SNP) in FTO and a genetic risk score (GRS) based on 32 SNPs with well-established associations with BMI. RESULTS: Linear regression analysis, as expected, showed that individuals carrying more risk alleles of FTO or having higher score of GRS had a higher BMI. Probit regression suggested that higher BMI is associated with increased risk of major depression. However, our two instrumental variable analyses did not support a causal relationship between higher BMI and major depression (FTO genotype: coefficient -0.03, 95% CI -0.18 to 0.13, P = 0.73; GRS: coefficient -0.02, 95% CI -0.11 to 0.07, P = 0.62). CONCLUSIONS: Our instrumental variable analyses did not support a causal relationship between higher BMI and major depression. The positive associations of higher BMI with major depression in probit regression analyses might be explained by reverse causality and/or residual confounding.

Multiplicative interaction of functional inflammasome genetic variants in determining the risk of gout.

INTRODUCTION: The acute gout flare results from a localised self-limiting innate immune response to monosodium urate (MSU) crystals deposited in joints in hyperuricaemic individuals. Activation of the caspase recruitment domain-containing protein 8 (CARD8) NOD-like receptor pyrin-containing 3 (NLRP3) inflammasome by MSU crystals and production of mature interleukin-1β (IL-1β) is central to acute gouty arthritis. However very little is known about genetic control of the innate immune response involved in acute gouty arthritis. Therefore our aim was to test functional single nucleotide polymorphism (SNP) variants in the toll-like receptor (TLR)-inflammasome-IL-1β axis for association with gout. METHODS: 1,494 gout cases of European and 863 gout cases of New Zealand (NZ) Polynesian (Māori and Pacific Island) ancestry were included. Gout was diagnosed by the 1977 ARA gout classification criteria. There were 1,030 Polynesian controls and 10,942 European controls including from the publicly-available Atherosclerosis Risk in Communities (ARIC) and Framingham Heart (FHS) studies. The ten SNPs were either genotyped by Sequenom MassArray or by Affymetrix SNP array or imputed in the ARIC and FHS datasets. Allelic association was done by logistic regression adjusting by age and sex with European and Polynesian data combined by meta-analysis. Sample sets were pooled for multiplicative interaction analysis, which was also adjusted by sample set. RESULTS: Eleven SNPs were tested in the TLR2, CD14, IL1B, CARD8, NLRP3, MYD88, P2RX7, DAPK1 and TNXIP genes. Nominally significant (P < 0.05) associations with gout were detected at CARD8 rs2043211 (OR = 1.12, P = 0.007), IL1B rs1143623 (OR = 1.10, P = 0.020) and CD14 rs2569190 (OR = 1.08; P = 0.036). There was significant multiplicative interaction between CARD8 and IL1B (P = 0.005), with the IL1B risk genotype amplifying the risk effect of CARD8. CONCLUSION: There is evidence for association of gout with functional variants in CARD8, IL1B and CD14. The gout-associated allele of IL1B increases expression of IL-1β - the multiplicative interaction with CARD8 would be consistent with a synergy of greater inflammasome activity (resulting from reduced CARD8) combined with higher levels of pre-IL-1β expression leading to increased production of mature IL-1β in gout.