A versatile gene-based test for genome-wide association studies

Beyond genome-wide significance: integrative approaches to the interpretation and extension of GWAS findings for alcohol use disorder

Alcohol use disorder (AUD) is a heritable complex behavior. Due to the highly polygenic nature of AUD, identifying genetic variants that comprise this heritable variation has proved to be challenging. With the exception of functional variants in alcohol metabolizing genes (e.g. ADH1B and ALDH2), few other candidate loci have been confidently linked to AUD. Genome-wide association studies (GWAS) of AUD and other alcohol-related phenotypes have either produced few hits with genome-wide significance or have failed to replicate on further study. These issues reinforce the complex nature of the genetic underpinnings for AUD and suggest that both GWAS studies with larger samples and additional analysis approaches that better harness the nominally significant loci in existing GWAS are needed. Here, we review approaches of interest in the post-GWAS era, including in silico functional analyses; functional partitioning of single nucleotide polymorphism heritability; aggregation of signal into genes and gene networks; and validation of identified loci, genes and gene networks in postmortem brain tissue and across species. These integrative approaches hold promise to illuminate our understanding of the biological basis of AUD; however, we recognize that the main challenge continues to be the extremely polygenic nature of AUD, which necessitates large samples to identify multiple loci associated with AUD liability.

Candidate gene analyses for acute pain and morphine analgesia after pediatric day surgery: African American versus European Caucasian ancestry and dose prediction limits

Acute pain and opioid analgesia demonstrate inter-individual variability and polygenic influence. In 241 children of African American and 277 of European Caucasian ancestry, we sought to replicate select candidate gene associations with morphine dose and postoperative pain and then to estimate dose prediction limits. Twenty-seven single-nucleotide polymorphisms (SNPs) from nine genes (ABCB1, ARRB2, COMT, DRD2, KCNJ6, MC1R, OPRD1, OPRM1, and UGT2B7) met selection criteria and were analyzed along with TAOK3. Few associations replicated: morphine dose (mcg/kg) in African American children and ABCB1 rs1045642 (A allele, β = -9.30, 95% CI: -17.25 to -1.35, p = 0.02) and OPRM1 rs1799971 (G allele, β = 23.19, 95% CI: 3.27-43.11, p = 0.02); KCNJ6 rs2211843 and high pain in African American subjects (T allele, OR 2.08, 95% CI: 1.17-3.71, p = 0.01) and in congruent European Caucasian pain phenotypes; and COMT rs740603 for high pain in European Caucasian subjects (A allele, OR: 0.69, 95% CI: 0.48-0.99, p = 0.046). With age, body mass index, and physical status as covariates, simple top SNP candidate gene models could explain theoretical maximums of 24.2% (European Caucasian) and 14.6% (African American) of morphine dose variances.

Eleven loci with new reproducible genetic associations with allergic disease risk

BACKGROUND

A recent genome-wide association study (GWAS) identified 99 loci that contain genetic risk variants shared between asthma, hay fever, and eczema. Many more risk loci shared between these common allergic diseases remain to be discovered, which could point to new therapeutic opportunities.

OBJECTIVE

We sought to identify novel risk loci shared between asthma, hay fever, and eczema by applying a gene-based test of association to results from a published GWAS that included data from 360,838 subjects.

METHODS

We used approximate conditional analysis to adjust the results from the published GWAS for the effects of the top risk variants identified in that study. We then analyzed the adjusted GWAS results with the EUGENE gene-based approach, which combines evidence for association with disease risk across regulatory variants identified in different tissues. Novel gene-based associations were followed up in an independent sample of 233,898 subjects from the UK Biobank study.

RESULTS

Of the 19,432 genes tested, 30 had a significant gene-based association at a Bonferroni-corrected P value of 2.5 × 10-6. Of these, 20 were also significantly associated (P < .05/30 = .0016) with disease risk in the replication sample, including 19 that were located in 11 loci not reported to contain allergy risk variants in previous GWASs. Among these were 9 genes with a known function that is directly relevant to allergic disease: FOSL2, VPRBP, IPCEF1, PRR5L, NCF4, APOBR, IL27, ATXN2L, and LAT. For 4 genes (eg, ATXN2L), a genetically determined decrease in gene expression was associated with decreased allergy risk, and therefore drugs that inhibit gene expression or function are predicted to ameliorate disease symptoms. The opposite directional effect was observed for 14 genes, including IL27, a cytokine known to suppress TH2 responses.

CONCLUSION

Using a gene-based approach, we identified 11 risk loci for allergic disease that were not reported in previous GWASs. Functional studies that investigate the contribution of the 19 associated genes to the pathophysiology of allergic disease and assess their therapeutic potential are warranted.

Prioritizing Crohn's disease genes by integrating association signals with gene expression implicates monocyte subsets

Genome-wide association studies have identified ~170 loci associated with Crohn's disease (CD) and defining which genes drive these association signals is a major challenge. The primary aim of this study was to define which CD locus genes are most likely to be disease related. We developed a gene prioritization regression model (GPRM) by integrating complementary mRNA expression datasets, including bulk RNA-Seq from the terminal ileum of 302 newly diagnosed, untreated CD patients and controls, and in stimulated monocytes. Transcriptome-wide association and co-expression network analyses were performed on the ileal RNA-Seq datasets, identifying 40 genome-wide significant genes. Co-expression network analysis identified a single gene module, which was substantially enriched for CD locus genes and most highly expressed in monocytes. By including expression-based and epigenetic information, we refined likely CD genes to 2.5 prioritized genes per locus from an average of 7.8 total genes. We validated our model structure using cross-validation and our prioritization results by protein-association network analyses, which demonstrated significantly higher CD gene interactions for prioritized compared with non-prioritized genes. Although individual datasets cannot convey all of the information relevant to a disease, combining data from multiple relevant expression-based datasets improves prediction of disease genes and helps to further understanding of disease pathogenesis.

Comparison of methods for multivariate gene-based association tests for complex diseases using common variants

Complex diseases are usually associated with multiple correlated phenotypes, and the analysis of composite scores or disease status may not fully capture the complexity (or multidimensionality). Joint analysis of multiple disease-related phenotypes in genetic tests could potentially increase power to detect association of a disease with common SNPs (or genes). Gene-based tests are designed to identify genes containing multiple risk variants that individually are weakly associated with a univariate trait. We combined three multivariate association tests (O'Brien method, TATES, and MultiPhen) with two gene-based association tests (GATES and VEGAS) and compared performance (type I error and power) of six multivariate gene-based methods using simulated data. Data (n = 2000) for genetic sequence and correlated phenotypes were simulated by varying causal variant proportions and phenotype correlations for various scenarios. These simulations showed that two multivariate association tests (TATES and MultiPhen, but not O'Brien) paired with VEGAS have inflated type I error in all scenarios, while the three multivariate association tests paired with GATES have correct type I error. MultiPhen paired with GATES has higher power than competing methods if the correlations among phenotypes are low (r < 0.57). We applied these gene-based association methods to a GWAS dataset from the Alzheimer's Disease Genetics Consortium containing three neuropathological traits related to Alzheimer disease (neuritic plaque, neurofibrillary tangles, and cerebral amyloid angiopathy) measured in 3500 autopsied brains. Gene-level significant evidence (P < 2.7 × 10^-6) was identified in a region containing three contiguous genes (TRAPPC12, TRAPPC12-AS1, ADI1) using O'Brien and VEGAS. Gene-wide significant associations were not observed in univariate gene-based tests.

Genome-Wide Gene-Environment Interaction Analysis Using Set-Based Association Tests

The identification of gene-environment interactions (G × E) may eventually guide health-related choices and medical interventions for complex diseases. More powerful methods must be developed to identify G × E. The “adaptive combination of Bayes factors method” (ADABF) has been proposed as a powerful genome-wide polygenic approach to detect G × E. In this work, we evaluate its performance when serving as a gene-based G × E test. We compare ADABF with six tests including the “Set-Based gene-EnviRonment InterAction test” (SBERIA), “gene-environment set association test” (GESAT), etc. With extensive simulations, SBERIA and ADABF are found to be more powerful than other G × E tests. However, SBERIA suffers from a power loss when 50% SNP main effects are in the same direction with the SNP × E interaction effects while 50% are in the opposite direction. We further applied these seven G × E methods to the Taiwan Biobank data to explore gene× alcohol interactions on blood pressure levels. The ADAMTS7P1 gene at chromosome 15q25.2 was detected to interact with alcohol consumption on diastolic blood pressure (p = 9.5 × 10⁻⁷, according to the GESAT test). At this gene, the P-values provided by other six tests all reached the suggestive significance level (p < 5 × 10⁻⁵). Regarding the computation time required for a genome-wide G × E analysis, SBERIA is the fastest method, followed by ADABF. Considering the validity, power performance, robustness, and computation time, ADABF is recommended for genome-wide G × E analyses.

GWAS and network analysis of co-occurring nicotine and alcohol dependence identifies significantly associated alleles and network

Alcohol dependence (AD) and nicotine dependence (ND) co-occur frequently (AD+ND). We integrated SNP-based, gene-based, and protein-protein interaction network analyses to identify shared risk genes or gene subnetworks for AD+ND in African Americans (AAs, N = 2,094) and European Americans (EAs, N = 1,207). The DSM-IV criterion counts for AD and ND were modeled as two dependent variables in a multivariate linear mixed model, and analyzed separately for the two populations. The most significant SNP was rs6579845 in EAs (p < 1.29 × 10^-8 ) in GM2A, which encodes GM2 ganglioside activator, and is a cis-expression quantitative locus that affects GM2A expression in blood and brain tissues. However, this SNP was not replicated in our another small sample (N = 678). We identified a subnetwork of 24 genes that contributed to the AD+ND criterion counts. In the gene-set analysis for the subnetwork in an independent sample, the Study of Addiction: Genetics and Environment project (predominately EAs), these 24 genes as a set differed in AD+ND versus control subjects in EAs (p = .041). Functional enrichment analysis for this subnetwork revealed that the gene enrichment involved primarily nerve growth factor pathways, and cocaine and amphetamine addiction. In conclusion, we identified a genome-wide significant variant at GM2A and a gene subnetwork underlying the genetic trait of shared AD+ND. These results increase our understanding of the shared (pleiotropic) genetic risk that underlies AD+ND.

Sex-specific gene and pathway modeling of inherited glioma risk

Background

To date, genome-wide association studies (GWAS) have identified 25 risk variants for glioma, explaining 30% of heritable risk. Most histologies occur with significantly higher incidence in males, and this difference is not explained by currently known risk factors. A previous GWAS identified sex-specific glioma risk variants, and this analysis aims to further elucidate risk variation by sex using gene- and pathway-based approaches.

Methods

Results from the Glioma International Case-Control Study were used as a testing set, and results from 3 GWAS were combined via meta-analysis and used as a validation set. Using summary statistics for nominally significant autosomal SNPs (P < 0.01 in a previous meta-analysis) and nominally significant X-chromosome SNPs (P < 0.01), 3 algorithms (Pascal, BimBam, and GATES) were used to generate gene scores, and Pascal was used to generate pathway scores. Results were considered statistically significant in the discovery set when P < 3.3 × 10-6 and in the validation set when P < 0.001 in 2 of 3 algorithms.

Results

Twenty-five genes within 5 regions and 19 genes within 6 regions reached statistical significance in at least 2 of 3 algorithms in males and females, respectively. EGFR was significantly associated with all glioma and glioblastoma in males only and a female-specific association in TERT, all of which remained nominally significant after conditioning on known risk loci. There were nominal associations with the BioCarta telomeres pathway in both males and females.

Conclusions

These results provide additional evidence that there may be differences by sex in genetic risk for glioma. Additional analyses may further elucidate the biological processes through which this risk is conferred.

Picomolar concentrations of oligomeric alpha-synuclein sensitizes TLR4 to play an initiating role in Parkinson's disease pathogenesis

Despite the wealth of genomic and transcriptomic data in Parkinson’s disease (PD), the initial molecular events are unknown. Using LD score regression analysis, we show significant enrichment in PD heritability within regulatory sites for LPS-activated monocytes and that TLR4 expression is highest within human substantia nigra, the most affected brain region, suggesting a role for TLR4 inflammatory responses. We then performed extended incubation of cells with physiological concentrations of small alpha-synuclein oligomers observing the development of a TLR4-dependent sensitized inflammatory response with time, including TNF-α production. ROS and cell death in primary neuronal cultures were significantly reduced by TLR4 antagonists revealing that an indirect inflammatory mechanism involving cytokines produced by glial cells makes a major contribution to neuronal death. Prolonged exposure to low levels of alpha-synuclein oligomers sensitizes TLR4 responsiveness in astrocytes and microglial, explaining how they become pro-inflammatory, and may be an early causative event in PD.

Strategies for Pathway Analysis Using GWAS and WGS Data

Single-allele study designs, commonly used in genome-wide association studies (GWAS) as well as the more recently developed whole genome sequencing (WGS) studies, are a standard approach for investigating the relationship of common variation within the human genome to a given phenotype of interest. However, single-allele association results published for many GWAS studies represent only the tip of the iceberg for the information that can be extracted from these datasets. The primary analysis strategy for GWAS entails association analysis in which only the single nucleotide polymorphisms (SNPs) with the strongest p-values are declared statistically significant due to issues arising from multiple testing and type I errors. Factors such as locus heterogeneity, epistasis, and multiple genes conferring small effects contribute to the complexity of the genetic models underlying phenotype expression. Thus, many biologically meaningful associations having lower effect sizes at individual genes are overlooked, making it difficult to separate true associations from a sea of false-positive associations. Organizing these individual SNPs into biologically meaningful groups to look at the overall effects of minor perturbations to genes and pathways is desirable. This pathway-based approach provides researchers with insight into the functional foundations of the phenotype being studied and allows testing of various genetic scenarios. © 2018 by John Wiley & Sons, Inc.

A combined analysis of genetically correlated traits identifies 187 loci and a role for neurogenesis and myelination in intelligence

Intelligence, or general cognitive function, is phenotypically and genetically correlated with many traits, including a wide range of physical, and mental health variables. Education is strongly genetically correlated with intelligence (r_g = 0.70). We used these findings as foundations for our use of a novel approach-multi-trait analysis of genome-wide association studies (MTAG; Turley et al. 2017)-to combine two large genome-wide association studies (GWASs) of education and intelligence, increasing statistical power and resulting in the largest GWAS of intelligence yet reported. Our study had four goals: first, to facilitate the discovery of new genetic loci associated with intelligence; second, to add to our understanding of the biology of intelligence differences; third, to examine whether combining genetically correlated traits in this way produces results consistent with the primary phenotype of intelligence; and, finally, to test how well this new meta-analytic data sample on intelligence predicts phenotypic intelligence in an independent sample. By combining datasets using MTAG, our functional sample size increased from 199,242 participants to 248,482. We found 187 independent loci associated with intelligence, implicating 538 genes, using both SNP-based and gene-based GWAS. We found evidence that neurogenesis and myelination-as well as genes expressed in the synapse, and those involved in the regulation of the nervous system-may explain some of the biological differences in intelligence. The results of our combined analysis demonstrated the same pattern of genetic correlations as those from previous GWASs of intelligence, providing support for the meta-analysis of these genetically-related phenotypes.

Gene-Based Nonparametric Testing of Interactions Using Distance Correlation Coefficient in Case-Control Association Studies

Among the various statistical methods for identifying gene⁻gene interactions in qualitative genome-wide association studies (GWAS), gene-based methods have recently grown in popularity because they confer advantages in both statistical power and biological interpretability. However, most of these methods make strong assumptions about the form of the relationship between traits and single-nucleotide polymorphisms, which result in limited statistical power. In this paper, we propose a gene-based method based on the distance correlation coefficient called gene-based gene-gene interaction via distance correlation coefficient (GBDcor). The distance correlation (dCor) is a measurement of the dependency between two random vectors with arbitrary, and not necessarily equal, dimensions. We used the difference in dCor in case and control datasets as an indicator of gene⁻gene interaction, which was based on the assumption that the joint distribution of two genes in case subjects and in control subjects should not be significantly different if the two genes do not interact. We designed a permutation-based statistical test to evaluate the difference between dCor in cases and controls for a pair of genes, and we provided the p-value for the statistic to represent the significance of the interaction between the two genes. In experiments with both simulated and real-world data, our method outperformed previous approaches in detecting interactions accurately.

Proper joint analysis of summary association statistics requires the adjustment of heterogeneity in SNP coverage pattern

As meta-analysis results published by consortia of genome-wide association studies (GWASs) become increasingly available, many association summary statistics-based multi-locus tests have been developed to jointly evaluate multiple single-nucleotide polymorphisms (SNPs) to reveal novel genetic architectures of various complex traits. The validity of these approaches relies on the accurate estimate of z-score correlations at considered SNPs, which in turn requires knowledge on the set of SNPs assessed by each study participating in the meta-analysis. However, this exact SNP coverage information is usually unavailable from the meta-analysis results published by GWAS consortia. In the absence of the coverage information, researchers typically estimate the z-score correlations by making oversimplified coverage assumptions. We show through real studies that such a practice can generate highly inflated type I errors, and we demonstrate the proper way to incorporate correct coverage information into multi-locus analyses. We advocate that consortia should make SNP coverage information available when posting their meta-analysis results, and that investigators who develop analytic tools for joint analyses based on summary data should pay attention to the variation in SNP coverage and adjust for it appropriately.

Heritability and Genome-Wide Association Study of Plasma Cholesterol in Chinese Adult Twins

Dyslipidemia represents a strong and independent risk factor for cardiovascular disease. Plasma cholesterol, such as total cholesterol (TC), low density lipoprotein cholesterol (LDL-C), and high density lipoprotein cholesterol (HDL-C), is the common indicator of diagnosing dyslipidemia. Here based on 382 Chinese twin pairs, we explored the magnitude of genetic impact on TC, HDL-C, LDL-C variation and further searched for genetic susceptibility loci for them using genome-wide association study (GWAS). The ACE model was the best fit model with additive genetic parameter (A) accounting for 26.6%, common or shared environmental parameter (C) accounting for 47.8%, unique/non-shared environmental parameter (E) accounting for 25.6% for the variance in HDL-C. The AE model was the best fit model for TC (A: 61.4%; E: 38.6%) and LDL-C (A: 65.5%; E: 34.5%). While no SNPs reached the genome-wide significance level (P < 5 × 10-8), 8, 14, 9 SNPs exceeded the suggestive significance level (P < 1 × 10-5) for TC, HDL-C, LDL-C, respectively. The promising genetic regions for TC, HDL-C, LDL-C were on chromosome 11 around rs7107698, chromosome 5 around rs12518218, chromosome 2 around rs10490120, respectively. Gene-based analysis found 1038, 1033 and 1090 genes nominally associated with TC, HDL-C, LDL-C (P < 0.05), especially FAF1, KLKB1 for TC, KLKB1 for HDL-C, and NTRK1, FAF1, SNTB2 for LDL-C, respectively. The number of common related genes among TC, HDL-C and LDL-C was 71, including FAF1, KLKB1, etc. Pathway enrichment analysis discovered known related pathways-zinc transporters, metal ion SLC transporters for TC, cell adhesion molecules CAMs, IL-6 signaling for HDL, FC epsilon RI signaling pathway, NFAT pathway for LDL, respectively. In conclusion, the TC and LDL-C level is moderately heritable and the HDL-C level is lowly heritable in Chinese population. The genomic loci, functional genes and pathways are identified to account for the heritability of plasma cholesterol level. Our findings provide important insights into plasma cholesterol molecular physiology and expect future research to replicate and validate our results.

Ancestry and genetic associations with bronchopulmonary dysplasia in preterm infants

Bronchopulmonary dysplasia in premature infants is a common and often severe lung disease with long-term sequelae. A genetic component is suspected but not fully defined. We performed an ancestry and genome-wide association study to identify variants, genes, and pathways associated with survival without bronchopulmonary dysplasia in 387 high-risk infants treated with inhaled nitric oxide in the Trial of Late Surfactant study. Global African genetic ancestry was associated with increased survival without bronchopulmonary dysplasia among infants of maternal self-reported Hispanic white race/ethnicity [odds ratio (OR) = 4.5, P = 0.01]. Admixture mapping found suggestive outcome associations with local African ancestry at chromosome bands 18q21 and 10q22 among infants of maternal self-reported African-American race/ethnicity. For all infants, the top individual variant identified was within the intron of NBL1, which is expressed in midtrimester lung and is an antagonist of bone morphogenetic proteins ( rs372271081 , OR = 0.17, P = 7.4 × 10-7). The protective allele of this variant was significantly associated with lower nitric oxide metabolites in the urine of non-Hispanic white infants ( P = 0.006), supporting a role in the racial differential response to nitric oxide. Interrogating genes upregulated in bronchopulmonary dysplasia lungs indicated association with variants in CCL18, a cytokine associated with fibrosis and interstitial lung disease, and pathway analyses implicated variation in genes involved in immune/inflammatory processes in response to infection and mechanical ventilation. Our results suggest that genetic variation related to lung development, drug metabolism, and immune response contribute to individual and racial/ethnic differences in respiratory outcomes following inhaled nitric oxide treatment of high-risk premature infants.

HisCoM-GGI: Hierarchical structural component analysis of gene-gene interactions

Although genome-wide association studies (GWAS) have successfully identified thousands of single nucleotide polymorphisms (SNPs) associated with common diseases, these observations are limited for fully explaining "missing heritability". Determining gene-gene interactions (GGI) are one possible avenue for addressing the missing heritability problem. While many statistical approaches have been proposed to detect GGI, most of these focus primarily on SNP-to-SNP interactions. While there are many advantages of gene-based GGI analyses, such as reducing the burden of multiple-testing correction, and increasing power by aggregating multiple causal signals across SNPs in specific genes, only a few methods are available. In this study, we proposed a new statistical approach for gene-based GGI analysis, "Hierarchical structural CoMponent analysis of Gene-Gene Interactions" (HisCoM-GGI). HisCoM-GGI is based on generalized structured component analysis, and can consider hierarchical structural relationships between genes and SNPs. For a pair of genes, HisCoM-GGI first effectively summarizes all possible pairwise SNP-SNP interactions into a latent variable, from which it then performs GGI analysis. HisCoM-GGI can evaluate both gene-level and SNP-level interactions. Through simulation studies, HisCoM-GGI demonstrated higher statistical power than existing gene-based GGI methods, in analyzing a GWAS of a Korean population for identifying GGI associated with body mass index. Resultantly, HisCoM-GGI successfully identified 14 potential GGI, two of which, (NCOR2 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mo>×</mml:mo></mml:math> SPOCK1) and (LINGO2 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mo>×</mml:mo></mml:math> ZNF385D) were successfully replicated in independent datasets. We conclude that HisCoM-GGI method may be a valuable tool for genome to identify GGI in missing heritability, allowing us to better understand the biological genetic mechanisms of complex traits. We conclude that HisCoM-GGI method may be a valuable tool for genome to identify GGI in missing heritability, allowing us to better understand biological genetic mechanisms of complex traits. An implementation of HisCoM-GGI can be downloaded from the website ( http://statgen.snu.ac.kr/software/hiscom-ggi ).

V-ATPases and osteoclasts: ambiguous future of V-ATPases inhibitors in osteoporosis

Vacuolar ATPases (V-ATPases) play a critical role in regulating extracellular acidification of osteoclasts and bone resorption. The deficiencies of subunit a3 and d2 of V-ATPases result in increased bone density in humans and mice. One of the traditional drug design strategies in treating osteoporosis is the use of subunit a3 inhibitor. Recent findings connect subunits H and G1 with decreased bone density. Given the controversial effects of ATPase subunits on bone density, there is a critical need to review the subunits of V-ATPase in osteoclasts and their functions in regulating osteoclasts and bone remodeling. In this review, we comprehensively address the following areas: information about all V-ATPase subunits and their isoforms; summary of V-ATPase subunits associated with human genetic diseases; V-ATPase subunits and osteopetrosis/osteoporosis; screening of all V-ATPase subunits variants in GEFOS data and in-house data; spectrum of V-ATPase subunits during osteoclastogenesis; direct and indirect roles of subunits of V-ATPases in osteoclasts; V-ATPase-associated signaling pathways in osteoclasts; interactions among V-ATPase subunits in osteoclasts; osteoclast-specific V-ATPase inhibitors; perspective of future inhibitors or activators targeting V-ATPase subunits in the treatment of osteoporosis.

Large-scale genome-wide enrichment analyses identify new trait-associated genes and pathways across 31 human phenotypes

Genome-wide association studies (GWAS) aim to identify genetic factors associated with phenotypes. Standard analyses test variants for associations individually. However, variant-level associations are hard to identify and can be difficult to interpret biologically. Enrichment analyses help address both problems by targeting sets of biologically related variants. Here we introduce a new model-based enrichment method that requires only GWAS summary statistics. Applying this method to interrogate 4,026 gene sets in 31 human phenotypes identifies many previously-unreported enrichments, including enrichments of endochondral ossification pathway for height, NFAT-dependent transcription pathway for rheumatoid arthritis, brain-related genes for coronary artery disease, and liver-related genes for Alzheimer’s disease. A key feature of our method is that inferred enrichments automatically help identify new trait-associated genes. For example, accounting for enrichment in lipid transport genes highlights association between MTTP and low-density lipoprotein levels, whereas conventional analyses of the same data found no significant variants near this gene.

In genome-wide association studies, variant-level associations are hard to identify and can be difficult to interpret biologically. Here, the authors develop a new model-based enrichment analysis method, and apply it to identify new associated genes, pathways and tissues across 31 human phenotypes.

Interactive effect between ATPase-related genes and early-life tobacco smoke exposure on bronchial hyper-responsiveness detected in asthma-ascertained families

BACKGROUND

A positional cloning study of bronchial hyper-responsiveness (BHR) at the 17p11 locus in the French Epidemiological study on the Genetics and Environment of Asthma (EGEA) families showed significant interaction between early-life environmental tobacco smoke (ETS) exposure and genetic variants located in DNAH9. This gene encodes the heavy chain subunit of axonemal dynein, which is involved with ATP in the motile cilia function.Our goal was to identify genetic variants at other genes interacting with ETS in BHR by investigating all genes belonging to the 'ATP-binding' and 'ATPase activity' pathways which include DNAH9, are targets of cigarette smoke and play a crucial role in the airway inflammation.

METHODS

Family-based interaction tests between ETS-exposed and unexposed BHR siblings were conducted in 388 EGEA families. Twenty single-nucleotide polymorphisms (SNP) showing interaction signals (p≤5.10^-3) were tested in the 253 Saguenay-Lac-Saint-Jean (SLSJ) families.

RESULTS

One of these SNPs was significantly replicated for interaction with ETS in SLSJ families (p=0.003). Another SNP reached the significance threshold after correction for multiple testing in the combined analysis of the two samples (p=10^-5). Results were confirmed using both a robust log-linear test and a gene-based interaction test.

CONCLUSION

The SNPs showing interaction with ETS belong to the ATP8A1 and ABCA1 genes, which play a role in the maintenance of asymmetry and homeostasis of lung membrane lipids.

Protein network analysis reveals selectively vulnerable regions and biological processes in FTD

OBJECTIVE

The neuroanatomical profile of behavioral variant frontotemporal dementia (bvFTD) suggests a common biological etiology of disease despite disparate pathologic causes; we investigated the genetic underpinnings of this selective regional vulnerability to identify new risk factors for bvFTD.

METHODS

We used recently developed analytical techniques designed to address the limitations of genome-wide association studies to generate a protein interaction network of 63 bvFTD risk genes. We characterized this network using gene expression data from healthy and diseased human brain tissue, evaluating regional network expression patterns across the lifespan as well as the cell types and biological processes most affected in bvFTD.

RESULTS

We found that bvFTD network genes show enriched expression across the human lifespan in vulnerable neuronal populations, are implicated in cell signaling, cell cycle, immune function, and development, and are differentially expressed in pathologically confirmed frontotemporal lobar degeneration cases. Five of the genes highlighted by our differential expression analyses, BAIAP2, ERBB3, POU2F2, SMARCA2, and CDC37, appear to be novel bvFTD risk loci.

CONCLUSIONS

Our findings suggest that the cumulative burden of common genetic variation in an interacting protein network expressed in specific brain regions across the lifespan may influence susceptibility to bvFTD.

Recent Advances in the Systems Biology of Aging

SIGNIFICANCE

Reductionist studies have contributed greatly to our understanding of the basic biology of aging in recent years but we still do not understand fundamental mechanisms for many identified drugs and pathways. Use of systems approaches will help us move forward in our understanding of aging. Recent Advances: Recent work described here has illustrated the power of systems biology to inform our understanding of aging through the study of (i) diet restriction, (ii) neurodegenerative disease, and (iii) biomarkers of aging.

CRITICAL ISSUES

Although we do not understand all of the individual genes and pathways that affect aging, as we continue to uncover more of them, we have now also begun to synthesize existing data using systems-level approaches, often to great effect. The three examples noted here all benefit from computational approaches that were unknown a few years ago, and from biological insights gleaned from multiple model systems, from aging laboratories as well as many other areas of biology.

FUTURE DIRECTIONS

Many new technologies, such as single-cell sequencing, advances in epigenetics beyond the methylome (specifically, assay for transposase-accessible chromatin with high throughput sequencing ), and multiomic network studies, will increase the reach of systems biologists. This suggests that approaches similar to those described here will continue to lead to striking findings, and to interventions that may allow us to delay some of the many age-associated diseases in humans; perhaps sooner that we expect. Antioxid. Redox Signal. 29, 973-984.

Common and Rare Genetic Risk Factors Converge in Protein Interaction Networks Underlying Schizophrenia

Hundreds of genomic loci have been identified with the recent advances of schizophrenia in genome-wide association studies (GWAS) and sequencing studies. However, the functional interactions among those genes remain largely unknown. We developed a network-based approach to integrate multiple genetic risk factors, which lead to the discovery of new susceptibility genes and causal sub-networks, or pathways in schizophrenia. We identified significantly and consistently over-represented pathways in the largest schizophrenia GWA studies, which are highly relevant to synaptic plasticity, neural development and signaling transduction, such as long-term potentiation, neurotrophin signaling pathway, and the ERBB signaling pathway. We also demonstrated that genes targeted by common SNPs are more likely to interact with genes harboring de novo mutations (DNMs) in the protein-protein interaction (PPI) network, suggesting a mutual interplay of both common and rare variants in schizophrenia. We further developed an edge-based search algorithm to identify the top-ranked gene modules associated with schizophrenia risk. Our results suggest that the N-methyl-D-aspartate receptor (NMDAR) interactome may play a leading role in the pathology of schizophrenia, as it is highly targeted by multiple types of genetic risk factors.

Integration of Molecular Interactome and Targeted Interaction Analysis to Identify a COPD Disease Network Module

The polygenic nature of complex diseases offers potential opportunities to utilize network-based approaches that leverage the comprehensive set of protein-protein interactions (the human interactome) to identify new genes of interest and relevant biological pathways. However, the incompleteness of the current human interactome prevents it from reaching its full potential to extract network-based knowledge from gene discovery efforts, such as genome-wide association studies, for complex diseases like chronic obstructive pulmonary disease (COPD). Here, we provide a framework that integrates the existing human interactome information with experimental protein-protein interaction data for FAM13A, one of the most highly associated genetic loci to COPD, to find a more comprehensive disease network module. We identified an initial disease network neighborhood by applying a random-walk method. Next, we developed a network-based closeness approach (C_AB) that revealed 9 out of 96 FAM13A interacting partners identified by affinity purification assays were significantly close to the initial network neighborhood. Moreover, compared to a similar method (local radiality), the C_AB approach predicts low-degree genes as potential candidates. The candidates identified by the network-based closeness approach were combined with the initial network neighborhood to build a comprehensive disease network module (163 genes) that was enriched with genes differentially expressed between controls and COPD subjects in alveolar macrophages, lung tissue, sputum, blood, and bronchial brushing datasets. Overall, we demonstrate an approach to find disease-related network components using new laboratory data to overcome incompleteness of the current interactome.

HSP90 as a global genetic modifier for male genital morphology in Drosophila melanogaster

The molecular chaperone protein HSP90 has been proposed to modulate genotype-phenotype relationship in a broad range of organisms. We explore the proposed genetic modifier effect of HSP90 through a genomewide analysis. Here, we show that HSP90 functions as a genetic modifier of genital morphology in Drosophila melanogaster. We identified a large number of single-nucleotide polymorphisms (SNPs) with an HSP90-dependent effect by using genome wide association analysis. We classified the SNPs into the ones under capacitance effect (smaller allelic effect under HSP90 inhibition) or the ones under potentiation effect (larger allelic effect under HSP90 inhibition). Although the majority of SNPs are under capacitance, there are a large number of SNPs under potentiation. This observation provides support for a model in which Hsp90 is not described exclusively as a "genetic capacitor," but is described more broadly as a "genetic modifier." Because the majority of the candidate genes estimated from SNPs with an HSP90-dependent effect in the current study has never been reported to interact with HSP90 directly, the global genetic modifier effect of HSP90 may be exhibited through epistatic interactions in gene regulatory networks.

Genomic Analyses of Visual Cognition: Perceptual Rivalry and Top-Down Control

Visual cognition in humans has traditionally been studied with cognitive behavioral methods and brain imaging, but much less with genetic methods. Perceptual rivalry, an important phenomenon in visual cognition, is the spontaneous perceptual alternation that occurs between two distinct interpretations of a physically constant visual stimulus (e.g., binocular rivalry stimuli) or a perceptually ambiguous stimulus (e.g., the Necker cube). The switching rate varies dramatically across individuals and can be voluntarily modulated by observers. Here, we adopted a genomic approach to systematically investigate the genetics underlying binocular rivalry, Necker cube rivalry and voluntary modulation of Necker cube rivalry in young Chinese adults (Homo sapiens, 81% female, 20 ± 1 years old) at multiple levels, including common single nucleotide polymorphism (SNP)-based heritability estimation, SNP-based genome-wide association study (GWAS), gene-based analysis, and pathway analysis. We performed a pilot GWAS in 2441 individuals and replicated it in an independent cohort of 943 individuals. Common SNP-based heritability was estimated to be 25% for spontaneous perceptual rivalry. SNPs rs184765639 and rs75595941 were associated with voluntary modulation, and imaging data suggested genotypic difference of rs184765639 in the surface area of the left caudal-middle frontal cortex. Additionally, converging evidence from multilevel analyses associated genes such as PRMT1 with perceptual switching rate, and MIR1178 with voluntary modulation strength. In summary, this study discovered specific genetic contributions to perceptual rivalry and its voluntary modulation in human beings. These findings may promote our understanding of psychiatric disorders, as perceptual rivalry is a potential psychiatric biomarker.SIGNIFICANCE STATEMENT Perceptual rivalry is an important visual phenomenon in which our perception of a physically constant visual input spontaneously switches between two different states. There are individual variations in perceptual switching rate and voluntary modulation strength. Our genomic analyses reveal several loci associated with these two kinds of variation. Because perceptual rivalry is thought to be relevant to and potentially an endophenotype for psychiatric disorders, these results may help understand not only visual cognition, but also psychiatric disorders.

Shared Biological Pathways Between Alzheimer's Disease and Ischemic Stroke

Alzheimer's disease (AD) and ischemic stroke (IS) are an immense socioeconomic burden worldwide. There is a possibility that shared genetic factors lead to their links at epidemiological and pathophysiological levels. Although recent genome-wide association studies (GWAS) have provided profound insights into the genetics of AD and IS, no shared genetic variants have been identified to date. This prompted us to initiate this study, which sought to identify shared pathways linking AD and IS. We took advantage of large-scale GWAS summary data of AD (17,008 AD cases and 37,154 controls) and IS (10,307 cases and 19,326 controls) to conduct pathway analyses using genetic pathways from multiple well-studied databases, including GO, KEGG, PANTHER, Reactome, and Wikipathways. Collectively, we discovered that AD and IS shared 179 GO categories (56 biological processes, 95 cellular components, and 28 molecular functions); and the following pathways: six KEGG pathways; two PANTHER pathways; four Reactome pathways; and one in Wikipathways pathway. The more fine-grained GO terms were mainly summarized into different functional categories: transcriptional and post-transcriptional regulation, synapse, endocytic membrane traffic through the endosomal system, signaling transduction, immune process, multi-organism process, protein catabolic metabolism, and cell adhesion. The shared pathways were roughly classified into three categories: immune system; cancer (NSCLC and glioma); and signal transduction pathways involving the cadherin signaling pathway, Wnt signaling pathway, G-protein signaling and downstream signaling mediated by phosphoinositides (PIPs). The majority of these common pathways linked to both AD and IS were supported by convincing evidence from the literature. In conclusion, our findings contribute to a better understanding of common biological mechanisms underlying AD and IS and serve as a guide to direct future research.

Prioritizing candidate genes post-GWAS using multiple sources of data for mastitis resistance in dairy cattle

Background

Improving resistance to mastitis, one of the costliest diseases in dairy production, has become an important objective in dairy cattle breeding. However, mastitis resistance is influenced by many genes involved in multiple processes, including the response to infection, inflammation, and post-infection healing. Low genetic heritability, environmental variations, and farm management differences further complicate the identification of links between genetic variants and mastitis resistance. Consequently, studies of the genetics of variation in mastitis resistance in dairy cattle lack agreement about the responsible genes.

Results

We associated 15,552,968 imputed whole-genome sequencing markers for 5147 Nordic Holstein cattle with mastitis resistance in a genome-wide association study (GWAS). Next, we augmented P-values for markers in genes in the associated regions using Gene Ontology terms, Kyoto Encyclopedia of Genes and Genomes pathway analysis, and mammalian phenotype database. To confirm results of gene-based analyses, we used gene expression data from E. coli-challenged cow udders. We identified 22 independent quantitative trait loci (QTL) that collectively explained 14% of the variance in breeding values for resistance to clinical mastitis (CM). Using association test statistics with multiple pieces of independent information on gene function and differential expression during bacterial infection, we suggested putative causal genes with biological relevance for 12 QTL affecting resistance to CM in dairy cattle.

Conclusion

Combining information on the nearest positional genes, gene-based analyses, and differential gene expression data from RNA-seq, we identified putative causal genes (candidate genes with biological evidence) in QTL for mastitis resistance in Nordic Holstein cattle. The same strategy can be applied for other traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5050-x) contains supplementary material, which is available to authorized users.

Multivariate analysis of genomics data to identify potential pleiotropic genes for type 2 diabetes, obesity and dyslipidemia using Meta-CCA and gene-based approach

Previous studies have demonstrated the genetic correlations between type 2 diabetes, obesity and dyslipidemia, and indicated that many genes have pleiotropic effects on them. However, these pleiotropic genes have not been well-defined. It is essential to identify pleiotropic genes using systematic approaches because systematically analyzing correlated traits is an effective way to enhance their statistical power. To identify potential pleiotropic genes for these three disorders, we performed a systematic analysis by incorporating GWAS (genome-wide associated study) datasets of six correlated traits related to type 2 diabetes, obesity and dyslipidemia using Meta-CCA (meta-analysis using canonical correlation analysis). Meta-CCA is an emerging method to systematically identify potential pleiotropic genes using GWAS summary statistics of multiple correlated traits. 2,720 genes were identified as significant genes after multiple testing (Bonferroni corrected p value < 0.05). Further, to refine the identified genes, we tested their relationship to the six correlated traits using VEGAS-2 (versatile gene-based association study-2). Only the genes significantly associated (Bonferroni corrected p value < 0.05) with more than one trait were kept. Finally, 25 genes (including two confirmed pleiotropic genes and eleven novel pleiotropic genes) were identified as potential pleiotropic genes. They were enriched in 5 pathways including the statin pathway and the PPAR (peroxisome proliferator-activated receptor) Alpha pathway. In summary, our study identified potential pleiotropic genes and pathways of type 2 diabetes, obesity and dyslipidemia, which may shed light on the common biological etiology and pathogenesis of these three disorders and provide promising insights for new therapies.

FDR-Corrected Sparse Canonical Correlation Analysis With Applications to Imaging Genomics

Reducing the number of false discoveries is presently one of the most pressing issues in the life sciences. It is of especially great importance for many applications in neuroimaging and genomics, where data sets are typically high-dimensional, which means that the number of explanatory variables exceeds the sample size. The false discovery rate (FDR) is a criterion that can be employed to address that issue. Thus it has gained great popularity as a tool for testing multiple hypotheses. Canonical correlation analysis (CCA) is a statistical technique that is used to make sense of the cross-correlation of two sets of measurements collected on the same set of samples (e.g., brain imaging and genomic data for the same mental illness patients), and sparse CCA extends the classical method to high-dimensional settings. Here, we propose a way of applying the FDR concept to sparse CCA, and a method to control the FDR. The proposed FDR correction directly influences the sparsity of the solution, adapting it to the unknown true sparsity level. Theoretical derivation as well as simulation studies show that our procedure indeed keeps the FDR of the canonical vectors below a user-specified target level. We apply the proposed method to an imaging genomics data set from the Philadelphia Neurodevelopmental Cohort. Our results link the brain connectivity profiles derived from brain activity during an emotion identification task, as measured by functional magnetic resonance imaging, to the corresponding subjects' genomic data.

Gene-Centric Analysis of Preeclampsia Identifies Maternal Association at PLEKHG1

The genetic susceptibility to preeclampsia, a pregnancy-specific complication with significant maternal and fetal morbidity, has been poorly characterized. To identify maternal genes associated with preeclampsia risk, we assembled 498 cases and 1864 controls of European ancestry from preeclampsia case-control collections in 5 different US sites (with additional matched population controls), genotyped samples on a cardiovascular gene-centric array composed of variants from ≈2000 genes selected based on prior genetic studies of cardiovascular and metabolic diseases and performed case-control genetic association analysis on 27 429 variants passing quality control. In silico replication testing of 9 lead signals with P<10-4 was performed in independent European samples from the SOPHIA (Study of Pregnancy Hypertension in Iowa) and Inova cohorts (212 cases, 456 controls). Multiethnic assessment of lead signals was then performed in samples of black (26 cases, 136 controls), Hispanic (132 cases, 468 controls), and East Asian (9 cases, 80 controls) ancestry. Multiethnic meta-analysis (877 cases, 3004 controls) revealed a study-wide statistically significant association of the rs9478812 variant in the pleiotropic PLEKHG1 gene (odds ratio, 1.40 [1.23-1.60]; Pmeta=5.90×10-7). The rs9478812 effect was even stronger in the subset of European cases with known early-onset preeclampsia (236 cases diagnosed <37 weeks, 1864 controls; odds ratio, 1.59 [1.27-1.98]; P=4.01×10-5). PLEKHG1 variants have previously been implicated in genome-wide association studies of blood pressure, body weight, and neurological disorders. Although larger studies are required to further define maternal preeclampsia heritability, this study identifies a novel maternal risk locus for further investigation.

Novel Myopia Genes and Pathways Identified From Syndromic Forms of Myopia

Purpose

To test the hypothesis that genes known to cause clinical syndromes featuring myopia also harbor polymorphisms contributing to nonsyndromic refractive errors.

Methods

Clinical phenotypes and syndromes that have refractive errors as a recognized feature were identified using the Online Mendelian Inheritance in Man (OMIM) database. One hundred fifty-four unique causative genes were identified, of which 119 were specifically linked with myopia and 114 represented syndromic myopia (i.e., myopia and at least one other clinical feature). Myopia was the only refractive error listed for 98 genes and hyperopia and the only refractive error noted for 28 genes, with the remaining 28 genes linked to phenotypes with multiple forms of refractive error. Pathway analysis was carried out to find biological processes overrepresented within these sets of genes. Genetic variants located within 50 kb of the 119 myopia-related genes were evaluated for involvement in refractive error by analysis of summary statistics from genome-wide association studies (GWAS) conducted by the CREAM Consortium and 23andMe, using both single-marker and gene-based tests.

Results

Pathway analysis identified several biological processes already implicated in refractive error development through prior GWAS analyses and animal studies, including extracellular matrix remodeling, focal adhesion, and axon guidance, supporting the research hypothesis. Novel pathways also implicated in myopia development included mannosylation, glycosylation, lens development, gliogenesis, and Schwann cell differentiation. Hyperopia was found to be linked to a different pattern of biological processes, mostly related to organogenesis. Comparison with GWAS findings further confirmed that syndromic myopia genes were enriched for genetic variants that influence refractive errors in the general population. Gene-based analyses implicated 21 novel candidate myopia genes (ADAMTS18, ADAMTS2, ADAMTSL4, AGK, ALDH18A1, ASXL1, COL4A1, COL9A2, ERBB3, FBN1, GJA1, GNPTG, IFIH1, KIF11, LTBP2, OCA2, POLR3B, POMT1, PTPN11, TFAP2A, ZNF469).

Conclusions

Common genetic variants within or nearby genes that cause syndromic myopia are enriched for variants that cause nonsyndromic, common myopia. Analysis of syndromic forms of refractive errors can provide new insights into the etiology of myopia and additional potential targets for therapeutic interventions.

Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals

Here we conducted a large-scale genetic association analysis of educational attainment in a sample of approximately 1.1 million individuals and identify 1,271 independent genome-wide-significant SNPs. For the SNPs taken together, we found evidence of heterogeneous effects across environments. The SNPs implicate genes involved in brain-development processes and neuron-to-neuron communication. In a separate analysis of the X chromosome, we identify 10 independent genome-wide-significant SNPs and estimate a SNP heritability of around 0.3% in both men and women, consistent with partial dosage compensation. A joint (multi-phenotype) analysis of educational attainment and three related cognitive phenotypes generates polygenic scores that explain 11-13% of the variance in educational attainment and 7-10% of the variance in cognitive performance. This prediction accuracy substantially increases the utility of polygenic scores as tools in research.

A genome scan for genes underlying adult body size differences between Central African hunter-gatherers and farmers

The evolutionary and biological bases of the Central African "pygmy" phenotype, a characteristic of rainforest hunter-gatherers defined by reduced body size compared with neighboring farmers, remain largely unknown. Here, we perform a joint investigation in Central African hunter-gatherers and farmers of adult standing height, sitting height, leg length, and body mass index (BMI), considering 358 hunter-gatherers and 169 farmers with genotypes for 153,798 SNPs. In addition to reduced standing heights, hunter-gatherers have shorter sitting heights and leg lengths and higher sitting/standing height ratios than farmers and lower BMI for males. Standing height, sitting height, and leg length are strongly correlated with inferred levels of farmer genetic ancestry, whereas BMI is only weakly correlated, perhaps reflecting greater contributions of non-genetic factors to body weight than to height. Single- and multi-marker association tests identify one region and eight genes associated with hunter-gatherer/farmer status, and 24 genes associated with the height-related traits. Many of these genes have putative functions consistent with roles in determining their associated traits and the pygmy phenotype, and they include three associated with standing height in non-Africans (PRKG1, DSCAM, MAGI2). We find evidence that European height-associated SNPs or variants in linkage disequilibrium with them contribute to standing- and sitting-height determination in Central Africans, but not to the differential status of hunter-gatherers and farmers. These findings provide new insights into the biological basis of the pygmy phenotype, and they highlight the potential of cross-population studies for exploring the genetic basis of phenotypes that vary naturally across populations.

Socio-demographic, Clinical, and Genetic Determinants of Quality of Life in Lung Cancer Patients

Patient reported health-related quality of life (QOL) is a major component of the overall well-being of cancer patients, with links to prognosis. In 6,420 lung cancer patients, we identified patient characteristics and genetic determinants of QOL. Patient responses from the SF-12 questionnaire was used to calculate normalized Physical Component Summary (PCS) and Mental Component Summary (MCS) scores. Further, we analyzed 218 single nucleotide polymorphisms (SNPs) in the p38 MAPK signaling pathway, a key mediator of response to cellular and environmental stress, as genetic determinants of QOL in a subset of the study population (N = 641). Trends among demographic factors for mean PCS and MCS included smoking status (PCS P_trend < 0.001, MCS P_trend < 0.001) and education (PCS P_trend < 0.001, MCS P_trend < 0.001). Similar relationships were seen for MCS. The homozygous rare genotype of MEF2B: rs2040562 showed an increased risk of a poor MCS (OR: 3.06, 95% CI: 1.05–8.92, P = 0.041). Finally, survival analysis showed that a low PCS or a MCS was associated with increased risks of five-year mortality (HR = 1.63, 95% CI: 1.51–1.77, HR = 1.23, 95% CI: 1.16–1.32, respectively) and there was a significant reduction in median survival time (P_log-rank < 0.001). These findings suggest that multiple factors contribute to QOL in lung cancer patients, and baseline QOL can impact survival.

Multiscale Analysis of Independent Alzheimer's Cohorts Finds Disruption of Molecular, Genetic, and Clinical Networks by Human Herpesvirus

Investigators have long suspected that pathogenic microbes might contribute to the onset and progression of Alzheimer's disease (AD) although definitive evidence has not been presented. Whether such findings represent a causal contribution, or reflect opportunistic passengers of neurodegeneration, is also difficult to resolve. We constructed multiscale networks of the late-onset AD-associated virome, integrating genomic, transcriptomic, proteomic, and histopathological data across four brain regions from human post-mortem tissue. We observed increased human herpesvirus 6A (HHV-6A) and human herpesvirus 7 (HHV-7) from subjects with AD compared with controls. These results were replicated in two additional, independent and geographically dispersed cohorts. We observed regulatory relationships linking viral abundance and modulators of APP metabolism, including induction of APBB2, APPBP2, BIN1, BACE1, CLU, PICALM, and PSEN1 by HHV-6A. This study elucidates networks linking molecular, clinical, and neuropathological features with viral activity and is consistent with viral activity constituting a general feature of AD.

Multi-level genomic analyses suggest new genetic variants involved in human memory

Development of high-throughput genotyping platforms provides an opportunity to identify new genetic elements related to complex cognitive functions. Taking advantage of multi-level genomic analysis, here we studied the genetic basis of human short-term (STM, n = 1623) and long-term (LTM, n = 1522) memory functions. Heritability estimation based on single nucleotide polymorphism showed moderate (61%, standard error 35%) heritability of short-term memory but almost zero heritability of long-term memory. We further performed a two-step genome-wide association study, but failed to find any SNPs that could pass genome-wide significance and survive replication at the same time. However, suggestive significance for rs7011450 was found in the shared component of the two STM tasks. Further inspections on its nearby gene zinc finger and at-hook domain containing and SNPs around this gene showed suggestive association with STM. In LTM, a polymorphism within branched chain amino acid transaminase 2 showed suggestive significance in the discovery cohort and has been replicated in another independent population of 1862. Furthermore, we performed a pathway analysis based on the current genomic data and found pathways including mTOR signaling and axon guidance significantly associated with STM capacity. These findings warrant further replication in other larger populations.

Controllability in an islet specific regulatory network identifies the transcriptional factor NFATC4, which regulates Type 2 Diabetes associated genes

Probing the dynamic control features of biological networks represents a new frontier in capturing the dysregulated pathways in complex diseases. Here, using patient samples obtained from a pancreatic islet transplantation program, we constructed a tissue-specific gene regulatory network and used the control centrality (Cc) concept to identify the high control centrality (HiCc) pathways, which might serve as key pathobiological pathways for Type 2 Diabetes (T2D). We found that HiCc pathway genes were significantly enriched with modest GWAS p-values in the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) study. We identified variants regulating gene expression (expression quantitative loci, eQTL) of HiCc pathway genes in islet samples. These eQTL genes showed higher levels of differential expression compared to non-eQTL genes in low, medium, and high glucose concentrations in rat islets. Among genes with highly significant eQTL evidence, NFATC4 belonged to four HiCc pathways. We asked if the expressions of T2D-associated candidate genes from GWAS and literature are regulated by Nfatc4 in rat islets. Extensive in vitro silencing of Nfatc4 in rat islet cells displayed reduced expression of 16, and increased expression of four putative downstream T2D genes. Overall, our approach uncovers the mechanistic connection of NFATC4 with downstream targets including a previously unknown one, TCF7L2, and establishes the HiCc pathways' relationship to T2D.

Integrated analysis of human genetic association study and mouse transcriptome suggests LBH and SHF genes as novel susceptible genes for amyloid-β accumulation in Alzheimer's disease

Alzheimer's disease (AD) is a common neurological disease that causes dementia in humans. Although the reports of associated pathological genes have been increasing, the molecular mechanism leading to the accumulation of amyloid-β (Aβ) in human brain is still not well understood. To identify novel genes that cause accumulation of Aβ in AD patients, we conducted an integrative analysis by combining a human genetic association study and transcriptome analysis in mouse brain. First, we examined genome-wide gene expression levels in the hippocampus, comparing them to amyloid Aβ level in mice with mixed genetic backgrounds. Next, based on a GWAS statistics obtained by a previous study with human AD subjects, we obtained gene-based statistics from the SNP-based statistics. We combined p values from the two types of analysis across orthologous gene pairs in human and mouse into one p value for each gene to evaluate AD susceptibility. As a result, we found five genes with significant p values in this integrated analysis among the 373 genes analyzed. We also examined the gene expression level of these five genes in the hippocampus of independent human AD cases and control subjects. Two genes, LBH and SHF, showed lower expression levels in AD cases than control subjects. This is consistent with the gene expression levels of both the genes in mouse which were negatively correlated with Aβ accumulation. These results, obtained from the integrative approach, suggest that LBH and SHF are associated with the AD pathogenesis.

Arc Requires PSD95 for Assembly into Postsynaptic Complexes Involved with Neural Dysfunction and Intelligence

Arc is an activity-regulated neuronal protein, but little is known about its interactions, assembly into multiprotein complexes, and role in human disease and cognition. We applied an integrated proteomic and genetic strategy by targeting a tandem affinity purification (TAP) tag and Venus fluorescent protein into the endogenous Arc gene in mice. This allowed biochemical and proteomic characterization of native complexes in wild-type and knockout mice. We identified many Arc-interacting proteins, of which PSD95 was the most abundant. PSD95 was essential for Arc assembly into 1.5-MDa complexes and activity-dependent recruitment to excitatory synapses. Integrating human genetic data with proteomic data showed that Arc-PSD95 complexes are enriched in schizophrenia, intellectual disability, autism, and epilepsy mutations and normal variants in intelligence. We propose that Arc-PSD95 postsynaptic complexes potentially affect human cognitive function.

•

TAP tag and purification of endogenous Arc protein complexes from the mouse brain

•

PSD95 is the major Arc binding protein, and both assemble into 1.5-MDa supercomplexes

•

PSD95 is essential for recruitment of Arc to synapses

•

Mutations and genetic variants in Arc-PSD95 are linked to cognition

Novel Methods for Family-Based Genetic Studies

The recent development of microarray and sequencing technology allows identification of disease susceptibility genes. Although the genome-wide association studies (GWAS) have successfully identified many genetic markers related to human diseases, the traditional statistical methods are not powerful to detect rare genetic markers. The rare genetic markers are usually grouped together and tested at the set level. One of such methods is the sequence kernel association test (SKAT), which has been commonly used in the rare genetic marker analysis. In recent publications, SKAT has been extended to be applicable for family-based rare variant analysis. Here, I present three published statistical approaches for family-based rare variant analysis for: 1. continuous traits, 2. binary traits, and 3. multiple correlated traits.

IGESS: a statistical approach to integrating individual-level genotype data and summary statistics in genome-wide association studies

Motivation

Results from genome-wide association studies (GWAS) suggest that a complex phenotype is often affected by many variants with small effects, known as 'polygenicity'. Tens of thousands of samples are often required to ensure statistical power of identifying these variants with small effects. However, it is often the case that a research group can only get approval for the access to individual-level genotype data with a limited sample size (e.g. a few hundreds or thousands). Meanwhile, summary statistics generated using single-variant-based analysis are becoming publicly available. The sample sizes associated with the summary statistics datasets are usually quite large. How to make the most efficient use of existing abundant data resources largely remains an open question.

Results

In this study, we propose a statistical approach, IGESS, to increasing statistical power of identifying risk variants and improving accuracy of risk prediction by i ntegrating individual level ge notype data and s ummary s tatistics. An efficient algorithm based on variational inference is developed to handle the genome-wide analysis. Through comprehensive simulation studies, we demonstrated the advantages of IGESS over the methods which take either individual-level data or summary statistics data as input. We applied IGESS to perform integrative analysis of Crohns Disease from WTCCC and summary statistics from other studies. IGESS was able to significantly increase the statistical power of identifying risk variants and improve the risk prediction accuracy from 63.2% ( ±0.4% ) to 69.4% ( ±0.1% ) using about 240 000 variants.

Availability and implementation

The IGESS software is available at https://github.com/daviddaigithub/IGESS .

Contact

zbxu@xjtu.edu.cn or xwan@comp.hkbu.edu.hk or eeyang@hkbu.edu.hk.

Supplementary information

Supplementary data are available at Bioinformatics online.

POLARIS: Polygenic LD-adjusted risk score approach for set-based analysis of GWAS data

Polygenic risk scores (PRSs) are a method to summarize the additive trait variance captured by a set of SNPs, and can increase the power of set-based analyses by leveraging public genome-wide association study (GWAS) datasets. PRS aims to assess the genetic liability to some phenotype on the basis of polygenic risk for the same or different phenotype estimated from independent data. We propose the application of PRSs as a set-based method with an additional component of adjustment for linkage disequilibrium (LD), with potential extension of the PRS approach to analyze biologically meaningful SNP sets. We call this method POLARIS: POlygenic Ld-Adjusted RIsk Score. POLARIS identifies the LD structure of SNPs using spectral decomposition of the SNP correlation matrix and replaces the individuals' SNP allele counts with LD-adjusted dosages. Using a raw genotype dataset together with SNP effect sizes from a second independent dataset, POLARIS can be used for set-based analysis. MAGMA is an alternative set-based approach employing principal component analysis to account for LD between markers in a raw genotype dataset. We used simulations, both with simple constructed and real LD-structure, to compare the power of these methods. POLARIS shows more power than MAGMA applied to the raw genotype dataset only, but less or comparable power to combined analysis of both datasets. POLARIS has the advantages that it produces a risk score per person per set using all available SNPs, and aims to increase power by leveraging the effect sizes from the discovery set in a self-contained test of association in the test dataset.

A molecule-based genetic association approach implicates a range of voltage-gated calcium channels associated with schizophrenia

Traditional genome-wide association studies (GWAS) have successfully detected genetic variants associated with schizophrenia. However, only a small fraction of heritability can be explained. Gene-set/pathway-based methods can overcome limitations arising from single nucleotide polymorphism (SNP)-based analysis, but most of them place constraints on size which may exclude highly specific and functional sets, like macromolecules. Voltage-gated calcium (Ca_v ) channels, belonging to macromolecules, are composed of several subunits whose encoding genes are located far away or even on different chromosomes. We combined information about such molecules with GWAS data to investigate how functional channels associated with schizophrenia. We defined a biologically meaningful SNP-set based on channel structure and performed an association study by using a validated method: SNP-set (sequence) kernel association test. We identified eight subtypes of Ca_v channels significantly associated with schizophrenia from a subsample of published data (N = 56,605), including the L-type channels (Ca_v 1.1, Ca_v 1.2, Ca_v 1.3), P-/Q-type Ca_v 2.1, N-type Ca_v 2.2, R-type Ca_v 2.3, T-type Ca_v 3.1, and Ca_v 3.3. Only genes from Ca_v 1.2 and Ca_v 3.3 have been implicated by the largest GWAS (N = 82,315). Each subtype of Ca_v channels showed relatively high chip heritability, proportional to the size of its constituent gene regions. The results suggest that abnormalities of Ca_v channels may play an important role in the pathophysiology of schizophrenia and these channels may represent appropriate drug targets for therapeutics. Analyzing subunit-encoding genes of a macromolecule in aggregate is a complementary way to identify more genetic variants of polygenic diseases. This study offers the potential of power for discovery the biological mechanisms of schizophrenia.

Analysis of the Influence of microRNAs in Lithium Response in Bipolar Disorder

Bipolar disorder (BD) is a common, highly heritable neuropsychiatric disease characterized by recurrent episodes of mania and depression. Lithium is the best-established long-term treatment for BD, even though individual response is highly variable. Evidence suggests that some of this variability has a genetic basis. This is supported by the largest genome-wide association study (GWAS) of lithium response to date conducted by the International Consortium on Lithium Genetics (ConLiGen). Recently, we performed the first genome-wide analysis of the involvement of miRNAs in BD and identified nine BD-associated miRNAs. However, it is unknown whether these miRNAs are also associated with lithium response in BD. In the present study, we therefore tested whether common variants at these nine candidate miRNAs contribute to the variance in lithium response in BD. Furthermore, we systematically analyzed whether any other miRNA in the genome is implicated in the response to lithium. For this purpose, we performed gene-based tests for all known miRNA coding genes in the ConLiGen GWAS dataset (n = 2,563 patients) using a set-based testing approach adapted from the versatile gene-based test for GWAS (VEGAS2). In the candidate approach, miR-499a showed a nominally significant association with lithium response, providing some evidence for involvement in both development and treatment of BD. In the genome-wide miRNA analysis, 71 miRNAs showed nominally significant associations with the dichotomous phenotype and 106 with the continuous trait for treatment response. A total of 15 miRNAs revealed nominal significance in both phenotypes with miR-633 showing the strongest association with the continuous trait (p = 9.80E-04) and miR-607 with the dichotomous phenotype (p = 5.79E-04). No association between miRNAs and treatment response to lithium in BD in either of the tested conditions withstood multiple testing correction. Given the limited power of our study, the investigation of miRNAs in larger GWAS samples of BD and lithium response is warranted.

Association Study and Fine-Mapping Major Histocompatibility Complex Analysis of Pemphigus Vulgaris in a Han Chinese Population

To identify possible additional genetic susceptibility loci for pemphigus vulgaris (PV), we performed a genome-wide association study of 240 PV patients and 1,031 control individuals, and we selected the top single nucleotide polymorphisms for replication in independent samples, with 252 patient samples and 1,852 control samples. We identified rs11218708 (P = 3.1 × 10^-8, odds ratio = 1.54) at chromosome locus 11q24.1 as significantly associated with PV. A fine-mapping analysis of PV risk in the major histocompatibility complex region showed three independent variants predisposed to PV using stepwise analysis: HLA-DRB1*14:04 (P = 2.47 × 10^-38, odds ratio = 6.28), rs7454108 at the TAP2 gene (P = 2.78 × 10^-12, odds ratio = 3.25), and rs1051336 at the HLA-DRA gene (P = 3.06 × 10^-6, odds ratio = 0.33). A systematic evaluation using gene- and pathway-based analyses showed a high tendency for PV susceptibility genes to be associated with autoimmunity. Our study highlights the involvement of immune-mediated processes in the pathophysiology of PV and illustrates the value of imputation to identify variants in the major histocompatibility complex region.

Study of 300,486 individuals identifies 148 independent genetic loci influencing general cognitive function

General cognitive function is a prominent and relatively stable human trait that is associated with many important life outcomes. We combine cognitive and genetic data from the CHARGE and COGENT consortia, and UK Biobank (total N = 300,486; age 16-102) and find 148 genome-wide significant independent loci (P < 5 × 10-8) associated with general cognitive function. Within the novel genetic loci are variants associated with neurodegenerative and neurodevelopmental disorders, physical and psychiatric illnesses, and brain structure. Gene-based analyses find 709 genes associated with general cognitive function. Expression levels across the cortex are associated with general cognitive function. Using polygenic scores, up to 4.3% of variance in general cognitive function is predicted in independent samples. We detect significant genetic overlap between general cognitive function, reaction time, and many health variables including eyesight, hypertension, and longevity. In conclusion we identify novel genetic loci and pathways contributing to the heritability of general cognitive function.

Cross-ancestry genome-wide association analysis of corneal thickness strengthens link between complex and Mendelian eye diseases

Central corneal thickness (CCT) is a highly heritable trait associated with complex eye diseases such as keratoconus and glaucoma. We perform a genome-wide association meta-analysis of CCT and identify 19 novel regions. In addition to adding support for known connective tissue-related pathways, pathway analyses uncover previously unreported gene sets. Remarkably, >20% of the CCT-loci are near or within Mendelian disorder genes. These included FBN1, ADAMTS2 and TGFB2 which associate with connective tissue disorders (Marfan, Ehlers-Danlos and Loeys-Dietz syndromes), and the LUM-DCN-KERA gene complex involved in myopia, corneal dystrophies and cornea plana. Using index CCT-increasing variants, we find a significant inverse correlation in effect sizes between CCT and keratoconus (r = -0.62, P = 5.30 × 10-5) but not between CCT and primary open-angle glaucoma (r = -0.17, P = 0.2). Our findings provide evidence for shared genetic influences between CCT and keratoconus, and implicate candidate genes acting in collagen and extracellular matrix regulation.

Integrative omics for health and disease

Advances in omics technologies - such as genomics, transcriptomics, proteomics and metabolomics - have begun to enable personalized medicine at an extraordinarily detailed molecular level. Individually, these technologies have contributed medical advances that have begun to enter clinical practice. However, each technology individually cannot capture the entire biological complexity of most human diseases. Integration of multiple technologies has emerged as an approach to provide a more comprehensive view of biology and disease. In this Review, we discuss the potential for combining diverse types of data and the utility of this approach in human health and disease. We provide examples of data integration to understand, diagnose and inform treatment of diseases, including rare and common diseases as well as cancer and transplant biology. Finally, we discuss technical and other challenges to clinical implementation of integrative omics.

Integrated analysis supports ATXN1 as a schizophrenia risk gene

Protein-protein interaction (PPI) is informative in identifying hidden disease risk genes that tend to interact with known risk genes usually working together in the same disease module. With the use of an integrated approach combining PPI information with pathway and expression analysis as well as genome-wide association study (GWAS), we intended to find new risk genes for schizophrenia (SCZ). We showed that ATXN1 was the only direct PPI partner of the know SCZ risk gene ZNF804A, and it also had direct PPIs with other 18 known SCZ risk genes. ATXN1 serves as one of the hub genes in the PPI network containing many known SCZ risk genes, and this network is significantly enriched for the MAPK signaling pathway. Further gene expression analysis indicated that ATXN1 is highly expressed in prefrontal cortex, and SCZ patients had significantly decreased expression compared with healthy controls. Finally, the published GWAS data supports an association of ATXN1 with SCZ as well as other psychiatric disorders though not reaching genome-wide significance. These convergent evidences support ATXN1 as a promising risk gene for SCZ, and the integrated approach serves as a useful tool for dissecting the genetic basis of psychiatric disorders.