Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci.

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWASs) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N = 52) and nonsmokers (N = 171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and correcting for multiple testing using a two-stage procedure. We found >2 million significant meQTL variants (padj < 0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects, and five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTL variants for 958 unique eGenes (padj < 0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTN1 and ITIH4 colocalized across all data types (GWAS, meQTL, and eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci.

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWAS) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N=52) and nonsmokers (N=171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and using a two-stage multiple testing approach with eigenMT and Bonferroni corrections. We found >2 million significant meQTL variants (padj<0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects; five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTLs for 958 unique eGenes (padj<0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTIN1 and ITIH4 colocalized across all data types (GWAS + meQTL + eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

GWAS and systems biology analysis of depressive symptoms among smokers from the COPDGene cohort.

BACKGROUND: Large sample GWAS is needed to identify genetic factors associated with depression. This study used genome-wide genotypic and phenotypic data from the COPDGene study to identify genetic risk factors for depression. METHODS: Data were from 9716 COPDGene subjects with ≥10 pack-year history. Depression was defined as antidepressant use and/or a HADS depression subscale score ≥8. Non-Hispanic White (6576) and African-American (3140) subsets were analyzed. A GWAS pipeline identified SNPs associated with depression in each group. Network analysis software analyzed gene interactions through common biological pathways, genetic interactions, and tissue-specific gene expression. RESULTS: The mean age was 59.4 years (SD 9.0) with 46.5% female subjects. Depression was in 24.7% of the NHW group (1622) and 12.5% of the AA group (391). No SNPs had genome-wide significance. One of the top SNPs, rs12036147 (p = 1.28 × 10-6), is near CHRM3. Another SNP was near MDGA2 (rs17118176, p = 3.52 × 10-6). Top genes formed networks for synaptic transmission with a statistically significant level of more co-expression in brain than other tissues, particularly in the basal ganglia (p = 1.00 × 10-4). LIMITATIONS: Limitations included a depression definition based on antidepressant use and a limited HADS score subgroup, which could increase false negatives in depressed patients not on antidepressants. Antidepressants used for smoking cessation in non-depressed patients could lead to false positives. CONCLUSIONS: Systems biology analysis identified statistically significant pathways whereby multiple genes influence depression. The gene set pathway analysis and COPDGene data can help investigate depression in future studies.

GLITTER: a web-based application for gene link inspection through tissue-specific coexpression.

Accumulating evidence supports the polygenic nature of most complex diseases, suggesting the involvement of many susceptibility genes with small effect sizes. Although hundreds of genes may underlie the genetic architecture of complex diseases, those involved in a given disease are probably not randomly distributed, but likely to be functionally related. Protein-protein interaction networks have been used to evaluate the functional relatedness of susceptibility genes. However, these networks do not account for tissue specificity, are limited to protein-coding genes, and are typically biased by incomplete biological knowledge. Here, we present Gene Link Inspector Through Tissue-specific coExpRession (GLITTER), a web-based application for assessing the functional relatedness of susceptibility genes, either coding or noncoding, according to tissue-specific gene expression profiles. GLITTER can also shed light on the specific tissues in which susceptibility genes might exert their functions. We further demonstrate examples of how GLITTER can evaluate the functional relatedness of susceptibility genes underlying schizophrenia and breast cancer, and provide clues about etiology.

FLAGS: A Flexible and Adaptive Association Test for Gene Sets Using Summary Statistics.

Genome-wide association studies (GWAS) have been widely used for identifying common variants associated with complex diseases. Despite remarkable success in uncovering many risk variants and providing novel insights into disease biology, genetic variants identified to date fail to explain the vast majority of the heritability for most complex diseases. One explanation is that there are still a large number of common variants that remain to be discovered, but their effect sizes are generally too small to be detected individually. Accordingly, gene set analysis of GWAS, which examines a group of functionally related genes, has been proposed as a complementary approach to single-marker analysis. Here, we propose a FL: exible and A: daptive test for G: ene S: ets (FLAGS), using summary statistics. Extensive simulations showed that this method has an appropriate type I error rate and outperforms existing methods with increased power. As a proof of principle, through real data analyses of Crohn's disease GWAS data and bipolar disorder GWAS meta-analysis results, we demonstrated the superior performance of FLAGS over several state-of-the-art association tests for gene sets. Our method allows for the more powerful application of gene set analysis to complex diseases, which will have broad use given that GWAS summary results are increasingly publicly available.

Association of the HTR2A gene with alcohol and heroin abuse.

Positive genetic associations of rs6313 (102T/C at exon 1) and rs6311 (-1438A/G) on the 5-hydroxytryptamine (serotonin) 2A receptor gene (HTR2A or 5-HT2A) were reported for alcohol and drug abuse; however, other association studies failed to produce consistent results supporting the susceptibility of the two single nucleotide polymorphisms (SNPs). To clarify the associations of the HTR2A gene with substance use disorders, we performed a meta-analysis based on the genotypes from the available candidate gene association studies of the two SNPs with alcohol and drug abuse from multiple populations. Evidence of association was found for HTR2A rs6313 in all the combined studies (e.g., allelic P = 0.0048 and OR 0.86, 95 % CI 0.77-0.95) and also in the combined studies of alcohol dependence (abuse) (e.g., allelic P = 0.0001 and OR 0.71, 95 % CI 0.59-0.85). The same association trend was also observed in the Study of Addiction: Genetics and Environment datasets. The meta-analysis supports a contribution of the HTR2A gene to the susceptibility to substance use disorders, particularly alcohol dependence.

Eating disorder predisposition is associated with ESRRA and HDAC4 mutations.

Anorexia nervosa and bulimia nervosa are common and severe eating disorders (EDs) of unknown etiology. Although genetic factors have been implicated in the psychopathology of EDs, a clear biological pathway has not been delineated. DNA from two large families affected by EDs was collected, and mutations segregating with illness were identified by whole-genome sequencing following linkage mapping or by whole-exome sequencing. In the first family, analysis of twenty members across three generations identified a rare missense mutation in the estrogen-related receptor α (ESRRA) gene that segregated with illness. In the second family, analysis of eight members across four generations identified a missense mutation in the histone deacetylase 4 (HDAC4) gene that segregated with illness. ESRRA and HDAC4 were determined to interact both in vitro in HeLa cells and in vivo in mouse cortex. Transcriptional analysis revealed that HDAC4 potently represses the expression of known ESRRA-induced target genes. Biochemical analysis of candidate mutations revealed that the identified ESRRA mutation decreased its transcriptional activity, while the HDAC4 mutation increased transcriptional repression of ESRRA. Our findings suggest that mutations that result in decreased ESRRA activity increase the risk of developing EDs.

Admixture mapping in lupus identifies multiple functional variants within IFIH1 associated with apoptosis, inflammation, and autoantibody production.

Systemic lupus erythematosus (SLE) is an inflammatory autoimmune disease with a strong genetic component. African-Americans (AA) are at increased risk of SLE, but the genetic basis of this risk is largely unknown. To identify causal variants in SLE loci in AA, we performed admixture mapping followed by fine mapping in AA and European-Americans (EA). Through genome-wide admixture mapping in AA, we identified a strong SLE susceptibility locus at 2q22-24 (LOD=6.28), and the admixture signal is associated with the European ancestry (ancestry risk ratio ~1.5). Large-scale genotypic analysis on 19,726 individuals of African and European ancestry revealed three independently associated variants in the IFIH1 gene: an intronic variant, rs13023380 [P(meta) = 5.20×10(-14); odds ratio, 95% confidence interval = 0.82 (0.78-0.87)], and two missense variants, rs1990760 (Ala946Thr) [P(meta) = 3.08×10(-7); 0.88 (0.84-0.93)] and rs10930046 (Arg460His) [P(dom) = 1.16×10(-8); 0.70 (0.62-0.79)]. Both missense variants produced dramatic phenotypic changes in apoptosis and inflammation-related gene expression. We experimentally validated function of the intronic SNP by DNA electrophoresis, protein identification, and in vitro protein binding assays. DNA carrying the intronic risk allele rs13023380 showed reduced binding efficiency to a cellular protein complex including nucleolin and lupus autoantigen Ku70/80, and showed reduced transcriptional activity in vivo. Thus, in SLE patients, genetic susceptibility could create a biochemical imbalance that dysregulates nucleolin, Ku70/80, or other nucleic acid regulatory proteins. This could promote antibody hypermutation and auto-antibody generation, further destabilizing the cellular network. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

Association of CHRNA4 polymorphisms with smoking behavior in two populations.

CHRNA4, the gene that encodes the nicotinic acetylcholine receptor α(4) subunit, is a potential candidate gene for nicotine dependence (ND). However, studies of the association of CHNRA4 with smoking behavior have shown inconsistent results. Our meta-analysis of linkage studies of smoking behavior identified a genome-wide significant linkage of the phenotype maximum number of cigarettes smoked in a 24-hour period to a region (20q13.12-q13.32) harboring CHRNA4. This motivated us to examine the association of CHRNA4 with smoking behavior in two independent samples. In this study, we examined five single nucleotide polymorphisms (SNPs) within CHRNA4 and three smoking-related behaviors: one quantitative trait [cigarettes smoked per day (CPD)], and two binary traits [DSM-IV diagnosis of ND and dichotomized Fagerstrom test of ND (FTND)], in 1,249 unrelated European-Americans (EAs) and 1,790 unrelated African-Americans (AAs). Using the combined sample with sex, age, and race as covariates, the synonymous SNP rs1044394 was significantly associated with ND (P = 0.001) and FTND (P = 0.01). Rs2236196, which has a low correlation with rs1044394, was also significantly associated with CPD (P = 0.003). The pattern of association for these SNPs was similar in AAs and EAs. After correction for multiple testing, the association between rs1044394 and ND in the combined sample remained significant (P = 0.033). In summary, our study supports association between CHRNA4 common variation and ND in AA and EA samples. Additional studies will be necessary to evaluate the role of rare variants at CHRNA4 for ND.

Multiple independent loci at chromosome 15q25.1 affect smoking quantity: a meta-analysis and comparison with lung cancer and COPD.

Recently, genetic association findings for nicotine dependence, smoking behavior, and smoking-related diseases converged to implicate the chromosome 15q25.1 region, which includes the CHRNA5-CHRNA3-CHRNB4 cholinergic nicotinic receptor subunit genes. In particular, association with the nonsynonymous CHRNA5 SNP rs16969968 and correlates has been replicated in several independent studies. Extensive genotyping of this region has suggested additional statistically distinct signals for nicotine dependence, tagged by rs578776 and rs588765. One goal of the Consortium for the Genetic Analysis of Smoking Phenotypes (CGASP) is to elucidate the associations among these markers and dichotomous smoking quantity (heavy versus light smoking), lung cancer, and chronic obstructive pulmonary disease (COPD). We performed a meta-analysis across 34 datasets of European-ancestry subjects, including 38,617 smokers who were assessed for cigarettes-per-day, 7,700 lung cancer cases and 5,914 lung-cancer-free controls (all smokers), and 2,614 COPD cases and 3,568 COPD-free controls (all smokers). We demonstrate statistically independent associations of rs16969968 and rs588765 with smoking (mutually adjusted p-values<10(-35) and <10(-8) respectively). Because the risk alleles at these loci are negatively correlated, their association with smoking is stronger in the joint model than when each SNP is analyzed alone. Rs578776 also demonstrates association with smoking after adjustment for rs16969968 (p<10(-6)). In models adjusting for cigarettes-per-day, we confirm the association between rs16969968 and lung cancer (p<10(-20)) and observe a nominally significant association with COPD (p = 0.01); the other loci are not significantly associated with either lung cancer or COPD after adjusting for rs16969968. This study provides strong evidence that multiple statistically distinct loci in this region affect smoking behavior. This study is also the first report of association between rs588765 (and correlates) and smoking that achieves genome-wide significance; these SNPs have previously been associated with mRNA levels of CHRNA5 in brain and lung tissue.

Meta-analysis of 15 genome-wide linkage scans of smoking behavior.

BACKGROUND: A genetic contribution to smoking behavior is well-established. To identify loci that increase the risk for smoking behavior, many genome-wide linkage scans have been performed with various smoking behavior assessments. Numerous putative susceptibility loci have been identified, but only a few of these were replicated in independent studies. METHODS: We used genome search meta-analysis (GSMA) to identify risk loci by pooling all available independent genome scan results on smoking behavior. Additionally, to minimize locus heterogeneity, subgroup analyses of the smoking behavior assessed by the Fagerstrom Test for Nicotine Dependence (FTND) and maximum number of cigarettes smoked in a 24-hour period (MaxCigs24) were carried out. Samples of European ancestry were also analyzed separately. RESULTS: A total number of 15 genome scan results were available for analysis, including 3404 families with 10,253 subjects. Overall, the primary GSMA across all smoking behavior identified a genome-wide suggestive linkage in chromosome 17q24.3-q25.3 (p(SR) = .001). A secondary analysis of FTND in European-ancestry samples (625 families with 1878 subjects) detected a genome-wide suggestive linkage in 5q33.1-5q35.2 (p(SR) = .0076). Subgroup analysis of MaxCigs24 (966 families with 3273 subjects) identified a genome-wide significant linkage in 20q13.12-q13.32 (p(SR) = .00041, p(OR) = .048), where a strongly supported nicotine dependence candidate gene, CHRNA4, is located. CONCLUSIONS: The regions identified in the current study deserve close attention and will be helpful for candidate gene identification or target re-sequencing studies in the future.

Different roles of MTHFR C677T and A1298C polymorphisms in colorectal adenoma and colorectal cancer: a meta-analysis.

Association studies on the MTHFR polymorphisms (C677T and A1298C) in colorectal cancer (CRC) and colorectal adenoma have shown conflicting results. We performed a meta-analysis to better assess the purported associations. Overall, the 677T allele (10,131 patients and 15,362 controls) showed a small but significant protective effect against CRC compared to the 677C allele [P=0.0003, odds ratio (OR)=0.93; 95% confidence interval (CI) 0.89-0.98, P=0.22 (for heterogeneity)] for a worldwide population. Meta-analyses of other genetic contrasts suggested that the 677T allele is more likely to affect CRC in a recessive genetic model worldwide (P<0.0001, OR=0.86; 95% CI 0.76-0.96, P=0.06) and in Asians (P=0.0005, OR=0.75; 95% CI 0.64-0.88, P=0.71). Similarly, we found a significantly decreased risk of CRC for 1298C polymorphism (4,764 CRC patients and 6,592 controls) for a recessive genetic model worldwide (P=0.005, OR=0.81; 95% CI 0.70-0.94, P=0.40) and in Caucasians (P=0.04, OR=0.75 95% CI 0.57-0.99, P=0.35). No evidence of association of C677T (4,616 patients and 6,338 controls) and A1298C (1,272 patients and 1,684 controls) with colorectal adenoma were found. The evidence accumulated suggests that MTHFR may represent a low-penetrance susceptible gene for CRC, and that the two polymorphisms might protect against colorectal adenoma developing into cancer. A larger single study is required to further evaluate gene-gene and gene-environment interactions for MTHFR polymorphisms and the cancer risk in a specific population.

CTLA4 polymorphisms and ophthalmopathy in Graves' disease patients: association study and meta-analysis.

Studies in the past have clearly established that cytotoxic T-lymphocyte antigen-4 (CTLA4) is a susceptible gene for Graves' disease (GD). However, association studies between the CTLA4 exon-1 +49A/G polymorphism and the risk of developing Graves' ophthalmopathy (GO) in GD patients have revealed conflicting results. In this study, associations of two CTLA4 polymorphisms (+49A/G and CT60) with GD risk and GO susceptibility in GD patients were investigated in a Chinese population. In addition, a meta-analysis was performed to better assess the purported association between the +49A/G polymorphism and GO susceptibility in GD patients. Our results demonstrated that both the +49A/G and CT60 polymorphisms were associated with GD susceptibility in the Chinese population. No significant association with GO susceptibility in GD patients was confirmed regardless of which polymorphism was tested individually. Similarly, the meta-analysis results provided minimal evidence about the role of the +49A/G polymorphism and GO risk in GD patients. Interestingly, haplotypic analysis demonstrated different scenarios concerning the role of CTLA4 in GO susceptibility in the Chinese GD patients. We found that the +49A-CT60G haplotype was marginally statistically associated with the increased risk of GO in GD patients (OR = 1.63, 95%CI 1.00-2.64, p = 0.05). In conclusion, our results suggested that CTLA4 might be involved in the susceptibility to GD in the Chinese population. Although neither +49A/G nor CT60 polymorphism was associated with the risk of GO in GD patients, the haplotypic analysis provided some evidence about its role in GO susceptibility in the Chinese GD patients. We suggest that more association studies recruiting haplotypic analysis should be performed to investigate the role of CTLA4 gene in GO susceptibility in patients from different nations.

DNA repair gene XRCC3 polymorphisms and cancer risk: a meta-analysis of 48 case-control studies.

The X-ray repair cross-complementing group 3 (XRCC3) is a highly suspected candidate gene for cancer susceptibility. However, association studies on the XRCC3 polymorphisms (4541A>G, Thr(241)Met, 17893A>G) in cancer have shown conflicting results. Therefore, we performed a meta-analysis to better assess the purported associations. Forty eight eligible case-control studies including 24,975 cancer patients and 34, 209 controls were selected for our meta-analysis. Overall, individuals carrying the XRCC3 Met/Met genotype showed a small cancer risk under a recessive genetic model. The subgroup and meta-regression analysis demonstrated different scenarios concerning the XRCC3 Met/Met genotype's role in cancer susceptibility for different subgroups. Specially, there was a significantly increased risk of breast cancer (OR, 1.14; P=0.0004; 95% CI, 1.06-1.23; P=0.37 for heterogeneity), elevated but not significant risk of cancer for head and neck, bladder, surprisingly, a significantly decreased risk of non-melanoma skin cancer (OR, 0.76; P=0.007; 95% CI, 0.62-0.93; P=0.61 for heterogeneity). A significantly elevated risk of cancer was observed in population-based case-control studies but not in nested or hospital based studies. Similarly, we found a significantly increased risk of cancer for A4541G and a decreased risk for A17893G under dominant genetic models. Our meta-analysis results support that the XRCC3 might represent a low-penetrance susceptible gene especially for cancer of breast, bladder, head and neck, and non-melanoma skin cancer. A single larger study should be required to further evaluate gene-gene and gene-environment interactions on XRCC3 polymorphisms and tissue-specific cancer risk in an ethnicity specific population.

Meta-analysis of the association of CTLA-4 exon-1 +49A/G polymorphism with rheumatoid arthritis.

Rheumatoid arthritis (RA) is a common autoimmune disease. Cytotoxic T lymphocyte-associated antigen 4 (CTLA-4) is a highly suspected candidate gene for RA susceptibility. However, association studies on the polymorphism of CTLA-4 exon-1 +49A/G in RA have shown conflicting results. Therefore, we performed a meta-analysis to better assess the purported association. In order to look for ethnic effect, we performed subgroup meta-analysis in populations of European descent and Asian descent. Meta-regression analysis was also performed to explore the possible heterogeneity between the two subgroups. Ten studies (11 comparisons) with the CTLA-4 exon-1 +49A/G genotyping on 2,315 patients with RA and 2,536 controls were selected for our meta-analysis. Overall, the fixed-effects odds ratio (OR) for the G versus A allele was 1.11 (P=0.02, 95% confidence interval (CI) 1.02-1.21), with no between-study heterogeneity. Subgroup and meta-regression analysis according to the ethnicity (European or Asian) demonstrated different scenarios concerning the CTLA-4 exon-1 +49A/G polymorphism's role in RA susceptibility for the two different subgroups. No effect of G on susceptibility was seen in European descent (five comparisons; OR=1.04, P=0.30, 95% CI 0.95-1.19; no significant between-study heterogeneity). However, there is a significant association in Asian descent under both fixed [OR=1.21, 95% CI (1.06-1.39), P=0.005] and random-effect models [OR=1.19, 95% CI (1.01-1.42), P=0.04]. Meta-regression analysis also supports the heterogeneity between the two subgroups (P=0.082). We also explored the role of this polymorphism on RA risk under other various interested genetic contrasts. These results further support that this polymorphism could not be a risk factor for Europeans. Interestingly, we find that in Asians the G allele has a greater tendency to cause RA in a recessive genetic model. However, sensitivity analysis showed that the combined result of Asian populations was unstable. In conclusion, our meta-analysis results suggest that CTLA-4 exon-1 +49G allele would not be a risk factor for RA in Europeans but might play a role in RA susceptibility for Asians.