Highly interconnected genes in disease-specific networks are enriched for disease-associated polymorphisms

Background

Complex diseases are associated with altered interactions between thousands of genes. We developed a novel method to identify and prioritize disease genes, which was generally applicable to complex diseases.

Results

We identified modules of highly interconnected genes in disease-specific networks derived from integrating gene-expression and protein interaction data. We examined if those modules were enriched for disease-associated SNPs, and could be used to find novel genes for functional studies. First, we analyzed publicly available gene expression microarray and genome-wide association study (GWAS) data from 13, highly diverse, complex diseases. In each disease, highly interconnected genes formed modules, which were significantly enriched for genes harboring disease-associated SNPs. To test if such modules could be used to find novel genes for functional studies, we repeated the analyses using our own gene expression microarray and GWAS data from seasonal allergic rhinitis. We identified a novel gene, FGF2, whose relevance was supported by functional studies using combined small interfering RNA-mediated knock-down and gene expression microarrays. The modules in the 13 complex diseases analyzed here tended to overlap and were enriched for pathways related to oncological, metabolic and inflammatory diseases. This suggested that this union of the modules would be associated with a general increase in susceptibility for complex diseases. Indeed, we found that this union was enriched with GWAS genes for 145 other complex diseases.

Conclusions

Modules of highly interconnected complex disease genes were enriched for disease-associated SNPs, and could be used to find novel genes for functional studies.     

Genetics of sputum gene expression in chronic obstructive pulmonary disease

Previous expression quantitative trait loci (eQTL) studies have performed genetic association studies for gene expression, but most of these studies examined lymphoblastoid cell lines from non-diseased individuals. We examined the genetics of gene expression in a relevant disease tissue from chronic obstructive pulmonary disease (COPD) patients to identify functional effects of known susceptibility genes and to find novel disease genes. By combining gene expression profiling on induced sputum samples from 131 COPD cases from the ECLIPSE Study with genomewide single nucleotide polymorphism (SNP) data, we found 4315 significant cis-eQTL SNP-probe set associations (3309 unique SNPs). The 3309 SNPs were tested for association with COPD in a genomewide association study (GWAS) dataset, which included 2940 COPD cases and 1380 controls. Adjusting for 3309 tests (p<1.5e-5), the two SNPs which were significantly associated with COPD were located in two separate genes in a known COPD locus on chromosome 15: CHRNA5 and IREB2. Detailed analysis of chromosome 15 demonstrated additional eQTLs for IREB2 mapping to that gene. eQTL SNPs for CHRNA5 mapped to multiple linkage disequilibrium (LD) bins. The eQTLs for IREB2 and CHRNA5 were not in LD. Seventy-four additional eQTL SNPs were associated with COPD at p<0.01. These were genotyped in two COPD populations, finding replicated associations with a SNP in PSORS1C1, in the HLA-C region on chromosome 6. Integrative analysis of GWAS and gene expression data from relevant tissue from diseased subjects has located potential functional variants in two known COPD genes and has identified a novel COPD susceptibility locus.     

Targeted sequencing of the 9p21.3 region reveals association with reduced disease risks in Ashkenazi Jewish centenarians

Genome-wide association studies (GWAS) have pinpointed the chromosomal locus 9p21.3 as a genetic hotspot for various age-related disorders. Common genetic variants in this locus are linked to multiple traits, including coronary artery diseases, cancers, and diabetes. Centenarians are known for their reduced risk and delayed onset of these conditions. To investigate whether this evasion of disease risks involves diminished genetic risks in the 9p21.3 locus, we sequenced this region in an Ashkenazi Jewish centenarian cohort (centenarians: n = 450, healthy controls: n = 500). Risk alleles associated with cancers, glaucoma, CAD, and T2D showed a significant depletion in centenarians. Furthermore, the risk and non-risk genotypes are linked to two distinct low-frequency variant profiles, enriched in controls and centenarians, respectively. Our findings provide evidence that the extreme longevity cohort is associated with collectively lower risks of multiple age-related diseases in the 9p21.3 locus.     

Sex differences in disease risk from reported genome-wide association study findings

Men and women differ in susceptibility to many diseases and in responses to treatment. Recent advances in genome-wide association studies (GWAS) provide a wealth of data for associating genetic profiles with disease risk; however, in general, these data have not been systematically probed for sex differences in gene-disease associations. Incorporating sex into the analysis of GWAS results can elucidate new relationships between single nucleotide polymorphisms (SNPs) and human disease. In this study, we performed a sex-differentiated analysis on significant SNPs from GWAS data of the seven common diseases studied by the Wellcome Trust Case Control Consortium. We employed and compared three methods: logistic regression, Woolf’s test of heterogeneity, and a novel statistical metric that we developed called permutation method to assess sex effects (PMASE). After correction for false discovery, PMASE finds SNPs that are significantly associated with disease in only one sex. These sexually dimorphic SNP-disease associations occur in Coronary Artery Disease and Crohn’s Disease. GWAS analyses that fail to consider sex-specific effects may miss discovering sexual dimorphism in SNP-disease associations that give new insights into differences in disease mechanism between men and women.     

Genome-wide association analysis of canine atopic dermatitis and identification of disease related SNPs

In humans, genome-wide association studies (GWAS) have been shown to be an effective and thorough approach for identifying polymorphisms associated with disease phenotypes. Here, we describe the first study to perform a genome-wide association study in canine atopic dermatitis (cAD) using the Illumina Canine SNP20 array, containing 22,362 single-nucleotide polymorphisms (SNPs). The aim of the study was to identify SNPs associated with cAD using affected and unaffected Golden Retrievers. Further validation studies were performed for potentially associated SNPs using Sequenom genotyping of larger numbers of cases and controls across eight breeds (Boxer, German Shepherd Dog, Labrador, Golden Retriever, Shiba Inu, Shih Tzu, Pit Bull, and West Highland White Terriers). Using meta-analysis, two SNPs were associated with cAD in all breeds tested. RS22114085 was identified as a susceptibility locus (p?=?0.00014, odds ratio?=?2) and RS23472497 as a protective locus (p?=?0.0015, odds ratio?=?0.6). Both of these SNPs were located in intergenic regions, and their effects have been demonstrated to be independent of each other, highlighting that further fine mapping and resequencing is required of these areas. Further, 12 SNPs were validated by Sequenom genotyping as associated with cAD, but these were not associated with all breeds. This study suggests that GWAS will be a useful approach for identifying genetic risk factors for cAD. Given the clinical heterogeneity within this condition and the likelihood that the relative genetic effect sizes are small, greater sample sizes and further studies will be required.     

Genome‐wide association study discovered candidate genes of Verticillium wilt resistance in upland cotton (Gossypium hirsutum L.)

Verticillium wilt (VW), caused by infection by Verticillium dahliae, is considered one of the most yield‐limiting diseases in cotton. To examine the genetic architecture of cotton VW resistance, we performed a genome‐wide association study (GWAS) using a panel of 299 accessions and 85 630 single nucleotide polymorphisms (SNPs) detected using the specific‐locus amplified fragment sequencing (SLAF‐seq) approach. Trait–SNP association analysis detected a total of 17 significant SNPs at P < 1.17 × 10^–5 (P = 1/85 630, –log₁₀P = 4.93); the peaks of SNPs associated with VW resistance on A10 were continuous and common in three environments (RDIG2015, RDIF2015 and RDIF2016). Haplotype block structure analysis predicted 22 candidate genes for VW resistance based on A10_99672586 with a minimum P‐value (–log₁₀P = 6.21). One of these genes (CG02) was near the significant SNP A10_99672586 (0.26 Mb), located in a 372‐kb haplotype block, and its Arabidopsis AT3G25510 homologues contain TIR‐NBS‐LRR domains that may be involved in disease resistance response. Real‐time quantitative PCR and virus‐induced gene silencing (VIGS) analysis showed that CG02 was specific to up‐regulation in the resistant (R) genotype Zhongzhimian2 (ZZM2) and that silenced plants were more susceptible to V. dahliae. These results indicate that CG02 is likely the candidate gene for resistance against V. dahliae in cotton. The identified locus or gene may serve as a promising target for genetic engineering and selection for improving resistance to VW in cotton.     

Elimination of p19ARF‐expressing cells protects against pulmonary emphysema in mice

Senescent cells accumulate in tissues during aging and are considered to underlie several aging‐associated phenotypes and diseases. We recently reported that the elimination of p19^ARF‐expressing senescent cells from lung tissue restored tissue function and gene expression in middle‐aged (12‐month‐old) mice. The aging of lung tissue increases the risk of pulmonary diseases such as emphysema, and cellular senescence is accelerated in emphysema patients. However, there is currently no direct evidence to show that cellular senescence promotes the pathology of emphysema, and the involvement of senescence in the development of this disease has yet to be clarified. We herein demonstrated that p19^ARF facilitated the development of pulmonary emphysema in mice. The elimination of p19^ARF‐expressing cells prevented lung tissue from elastase‐induced lung dysfunction. These effects appeared to depend on reduced pulmonary inflammation, which is enhanced after elastase stimulation. Furthermore, the administration of a senolytic drug that selectively kills senescent cells attenuated emphysema‐associated pathologies. These results strongly suggest the potential of senescent cells as therapeutic/preventive targets for pulmonary emphysema.     

Genome-Wide Scan Informed by Age-Related Disease Identifies Loci for Exceptional Human Longevity

We developed a new statistical framework to find genetic variants associated with extreme longevity. The method, informed GWAS (iGWAS), takes advantage of knowledge from large studies of age-related disease in order to narrow the search for SNPs associated with longevity. To gain support for our approach, we first show there is an overlap between loci involved in disease and loci associated with extreme longevity. These results indicate that several disease variants may be depleted in centenarians versus the general population. Next, we used iGWAS to harness information from 14 meta-analyses of disease and trait GWAS to identify longevity loci in two studies of long-lived humans. In a standard GWAS analysis, only one locus in these studies is significant (APOE/TOMM40) when controlling the false discovery rate (FDR) at 10%. With iGWAS, we identify eight genetic loci to associate significantly with exceptional human longevity at FDR < 10%. We followed up the eight lead SNPs in independent cohorts, and found replication evidence of four loci and suggestive evidence for one more with exceptional longevity. The loci that replicated (FDR < 5%) included APOE/TOMM40 (associated with Alzheimer’s disease), CDKN2B/ANRIL (implicated in the regulation of cellular senescence), ABO (tags the O blood group), and SH2B3/ATXN2 (a signaling gene that extends lifespan in Drosophila and a gene involved in neurological disease). Our results implicate new loci in longevity and reveal a genetic overlap between longevity and age-related diseases and traits, including coronary artery disease and Alzheimer’s disease. iGWAS provides a new analytical strategy for uncovering SNPs that influence extreme longevity, and can be applied more broadly to boost power in other studies of complex phenotypes.     

The emerging role of alternative splicing in senescence and aging

Deregulation of precursor mRNA splicing is associated with many illnesses and has been linked to age‐related chronic diseases. Here we review recent progress documenting how defects in the machinery that performs intron removal and controls splice site selection contribute to cellular senescence and organismal aging. We discuss the functional association linking p53, IGF‐1, SIRT1, and ING‐1 splice variants with senescence and aging, and review a selection of splicing defects occurring in accelerated aging (progeria), vascular aging, and Alzheimer's disease. Overall, it is becoming increasingly clear that changes in the activity of splicing factors and in the production of key splice variants can impact cellular senescence and the aging phenotype.     

Genetic evidence for common pathways in human age‐related diseases

Aging is the single largest risk factor for chronic disease. Studies in model organisms have identified conserved pathways that modulate aging rate and the onset and progression of multiple age‐related diseases, suggesting that common pathways of aging may influence age‐related diseases in humans as well. To determine whether there is genetic evidence supporting the notion of common pathways underlying age‐related diseases, we analyzed the genes and pathways found to be associated with five major categories of age‐related disease using a total of 410 genomewide association studies (GWAS). While only a small number of genes are shared among all five disease categories, those found in at least three of the five major age‐related disease categories are highly enriched for apoliprotein metabolism genes. We found that a more substantial number of gene ontology (GO) terms are shared among the 5 age‐related disease categories and shared GO terms include canonical aging pathways identified in model organisms, such as nutrient‐sensing signaling, translation, proteostasis, stress responses, and genome maintenance. Taking advantage of the vast amount of genetic data from the GWAS, our findings provide the first direct evidence that conserved pathways of aging simultaneously influence multiple age‐related diseases in humans as has been demonstrated in model organisms.     

Genome-wide association studies reveal the role of polymorphisms affecting factor H binding protein expression in host invasion by Neisseria meningitidis

Many invasive bacterial diseases are caused by organisms that are ordinarily harmless components of the human microbiome. Effective interventions against these microbes require an understanding of the processes whereby symbiotic or commensal relationships transition into pathology. Here, we describe bacterial genome-wide association studies (GWAS) of Neisseria meningitidis, a common commensal of the human respiratory tract that is nevertheless a leading cause of meningitis and sepsis. An initial GWAS discovered bacterial genetic variants, including single nucleotide polymorphisms (SNPs), associated with invasive meningococcal disease (IMD) versus carriage in several loci across the meningococcal genome, encoding antigens and other extracellular components, confirming the polygenic nature of the invasive phenotype. In particular, there was a significant peak of association around the fHbp locus, encoding factor H binding protein (fHbp), which promotes bacterial immune evasion of human complement by recruiting complement factor H (CFH) to the meningococcal surface. The association around fHbp with IMD was confirmed by a validation GWAS, and we found that the SNPs identified in the validation affected the 5’ region of fHbp mRNA, altering secondary RNA structures, thereby increasing fHbp expression and enhancing bacterial escape from complement-mediated killing. This finding is consistent with the known link between complement deficiencies and CFH variation with human susceptibility to IMD. These observations demonstrate the importance of human and bacterial genetic variation across the fHbp:CFH interface in determining IMD susceptibility, the transition from carriage to disease.     

Blood cis-eQTL Analysis Fails to Identify Novel Association Signals among Sub-Threshold Candidates from Genome-Wide Association Studies in Restless Legs Syndrome

Restless legs syndrome (RLS) is a common neurologic disorder characterized by nightly dysesthesias affecting the legs primarily during periods of rest and relieved by movement. RLS is a complex genetic disease and susceptibility factors in six genomic regions have been identified by means of genome-wide association studies (GWAS). For some complex genetic traits, expression quantitative trait loci (eQTLs) are enriched among trait-associated single nucleotide polymorphisms (SNPs). With the aim of identifying new genetic susceptibility factors for RLS, we assessed the 332 best-associated SNPs from the genome-wide phase of the to date largest RLS GWAS for cis-eQTL effects in peripheral blood from individuals of European descent. In 740 individuals belonging to the KORA general population cohort, 52 cis-eQTLs with p_nominal<10⁻³ were identified, while in 976 individuals belonging to the SHIP-TREND general population study 53 cis-eQTLs with p_nominal<10⁻³ were present. 23 of these cis-eQTLs overlapped between the two cohorts. Subsequently, the twelve of the 23 cis-eQTL SNPs, which were not located at an already published RLS-associated locus, were tested for association in 2449 RLS cases and 1462 controls. The top SNP, located in the DET1 gene, was nominally significant (p<0.05) but did not withstand correction for multiple testing (p = 0.42). Although a similar approach has been used successfully with regard to other complex diseases, we were unable to identify new genetic susceptibility factor for RLS by adding this novel level of functional assessment to RLS GWAS data.     

Alzheimer’s Disease Risk Polymorphisms Regulate Gene Expression in the ZCWPW1 and the CELF1 Loci

Late onset Alzheimer’s disease (LOAD) is a genetically complex and clinically heterogeneous disease. Recent large-scale genome wide association studies (GWAS) have identified more than twenty loci that modify risk for AD. Despite the identification of these loci, little progress has been made in identifying the functional variants that explain the association with AD risk. Thus, we sought to determine whether the novel LOAD GWAS single nucleotide polymorphisms (SNPs) alter expression of LOAD GWAS genes and whether expression of these genes is altered in AD brains. The majority of LOAD GWAS SNPs occur in gene dense regions under large linkage disequilibrium (LD) blocks, making it unclear which gene(s) are modified by the SNP. Thus, we tested for brain expression quantitative trait loci (eQTLs) between LOAD GWAS SNPs and SNPs in high LD with the LOAD GWAS SNPs in all of the genes within the GWAS loci. We found a significant eQTL between rs1476679 and PILRB and GATS, which occurs within the ZCWPW1 locus. PILRB and GATS expression levels, within the ZCWPW1 locus, were also associated with AD status. Rs7120548 was associated with MTCH2 expression, which occurs within the CELF1 locus. Additionally, expression of several genes within the CELF1 locus, including MTCH2, were highly correlated with one another and were associated with AD status. We further demonstrate that PILRB, as well as other genes within the GWAS loci, are most highly expressed in microglia. These findings together with the function of PILRB as a DAP12 receptor supports the critical role of microglia and neuroinflammation in AD risk.     

Association between putative functional variants in the PSMB9 gene and risk of melanoma – re‐analysis of published melanoma genome‐wide association studies

To mine possibly hidden causal single‐nucleotide polymorphisms (SNPs) of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome‐wide association study (GWAS) data set with 1804 melanoma cases and 1026 cancer‐free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in two other GWAS data sets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL data set (P = 0.013 and 0.004, respectively), but failed in validation using the Australian data set. Genotype–phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression level of PSMB9. Further experiments revealed that SNP rs2071480, which is in high LD with rs1351383 and rs2127675, may have a weak effect on the promoter activity of PSMB9. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility.