GSA-SNP: a general approach for gene set analysis of polymorphisms

Clinical effects of novel susceptibility genes for beta-amyloid: a gene-based association study in the Korean population

Amyloid-beta (Aβ) is a pathological hallmark of Alzheimer's disease (AD). We aimed to identify genes related to Aβ uptake in the Korean population and investigate the effects of these novel genes on clinical outcomes, including neurodegeneration and cognitive impairments. We recruited a total of 759 Korean participants who underwent neuropsychological tests, brain magnetic resonance imaging, 18F-flutemetamol positron emission tomography, and microarray genotyping data. We performed gene-based association analysis, and also performed expression quantitative trait loci and network analysis. In genome-wide association studies, no single nucleotide polymorphism (SNP) passed the genome-wide significance threshold. In gene-based association analysis, six genes (LCMT1, SCRN2, LRRC46, MRPL10, SP6, and OSBPL7) were significantly associated with Aβ standardised uptake value ratio in the brain. The three most significant SNPs (rs4787307, rs9903904, and rs11079797) on these genes are associated with the regulation of the LCMT1, OSBPL7, and SCRN2 genes, respectively. These SNPs are involved in decreasing hippocampal volume and cognitive scores by mediating Aβ uptake. The 19 enriched gene sets identified by pathway analysis included axon and chemokine activity. Our findings suggest novel susceptibility genes associated with the uptake of Aβ, which in turn leads to worse clinical outcomes. Our findings might lead to the discovery of new AD treatment targets.

Genome-wide analysis identifies a novel LINC-PINT splice variant associated with vascular amyloid pathology in Alzheimer's disease

Cerebral amyloid angiopathy (CAA) contributes to accelerated cognitive decline in Alzheimer’s disease (AD) dementia and is a common finding at autopsy. The APOEε4 allele and male sex have previously been reported to associate with increased CAA in AD. To inform biomarker and therapeutic target discovery, we aimed to identify additional genetic risk factors and biological pathways involved in this vascular component of AD etiology. We present a genome-wide association study of CAA pathology in AD cases and report sex- and APOE-stratified assessment of this phenotype. Genome-wide genotypes were collected from 853 neuropathology-confirmed AD cases scored for CAA across five brain regions, and imputed to the Haplotype Reference Consortium panel. Key variables and genome-wide genotypes were tested for association with CAA in all individuals and in sex and APOEε4 stratified subsets. Pathway enrichment was run for each of the genetic analyses. Implicated loci were further investigated for functional consequences using brain transcriptome data from 1,186 samples representing seven brain regions profiled as part of the AMP-AD consortium. We confirmed association of male sex, AD neuropathology and APOEε4 with increased CAA, and identified a novel locus, LINC-PINT, associated with lower CAA amongst APOEε4-negative individuals (rs10234094-C, beta = −3.70 [95% CI −0.49—−0.24]; p = 1.63E-08). Transcriptome profiling revealed higher LINC-PINT expression levels in AD cases, and association of rs10234094-C with altered LINC-PINT splicing. Pathway analysis indicates variation in genes involved in neuronal health and function are linked to CAA in AD patients. Further studies in additional and diverse cohorts are needed to assess broader translation of our findings.

Genome-wide association study identifies susceptibility loci of brain atrophy to NFIA and ST18 in Alzheimer's disease

To identify genetic variants influencing cortical atrophy in Alzheimer's disease (AD), we performed genome-wide association studies (GWAS) of mean cortical thicknesses in 17 AD-related brain. In this study, we used neuroimaging and genetic data of 919 participants from the Alzheimer's Disease Neuroimaging Initiative cohort, which include 268 cognitively normal controls, 488 mild cognitive impairment, 163 AD individuals. We performed GWAS with 3,041,429 single nucleotide polymorphisms (SNPs) for cortical thickness. The results of GWAS indicated that rs10109716 in ST18 (ST18 C2H2C-type zinc finger transcription factor) and rs661526 in NFIA (nuclear factor I A) genes are significantly associated with mean cortical thicknesses of the left inferior frontal gyrus and left parahippocampal gyrus, respectively. The rs661526 regulates the expression levels of NFIA in the substantia nigra and frontal cortex and rs10109716 regulates the expression levels of ST18 in the thalamus. These results suggest a crucial role of identified genes for cortical atrophy and could provide further insights into the genetic basis of AD.

Regional imaging genetic enrichment analysis

MOTIVATION

Brain imaging genetics aims to reveal genetic effects on brain phenotypes, where most studies examine phenotypes defined on anatomical or functional regions of interest (ROIs) given their biologically meaningful interpretation and modest dimensionality compared with voxelwise approaches. Typical ROI-level measures used in these studies are summary statistics from voxelwise measures in the region, without making full use of individual voxel signals.

RESULTS

In this article, we propose a flexible and powerful framework for mining regional imaging genetic associations via voxelwise enrichment analysis, which embraces the collective effect of weak voxel-level signals and integrates brain anatomical annotation information. Our proposed method achieves three goals at the same time: (i) increase the statistical power by substantially reducing the burden of multiple comparison correction; (ii) employ brain annotation information to enable biologically meaningful interpretation and (iii) make full use of fine-grained voxelwise signals. We demonstrate our method on an imaging genetic analysis using data from the Alzheimer's Disease Neuroimaging Initiative, where we assess the collective regional genetic effects of voxelwise FDG-positron emission tomography measures between 116 ROIs and 565 373 single-nucleotide polymorphisms. Compared with traditional ROI-wise and voxelwise approaches, our method identified 2946 novel imaging genetic associations in addition to 33 ones overlapping with the two benchmark methods. In particular, two newly reported variants were further supported by transcriptome evidences from region-specific expression analysis. This demonstrates the promise of the proposed method as a flexible and powerful framework for exploring imaging genetic effects on the brain.

AVAILABILITY AND IMPLEMENTATION

The R code and sample data are freely available at https://github.com/lshen/RIGEA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Models of the Gene Must Inform Data-Mining Strategies in Genomics

The gene is a fundamental concept of genetics, which emerged with the Mendelian paradigm of heredity at the beginning of the 20th century. However, the concept has since diversified. Somewhat different narratives and models of the gene developed in several sub-disciplines of genetics, that is in classical genetics, population genetics, molecular genetics, genomics, and, recently, also, in systems genetics. Here, I ask how the diversity of the concept impacts data-integration and data-mining strategies for bioinformatics, genomics, statistical genetics, and data science. I also consider theoretical background of the concept of the gene in the ideas of empiricism and experimentalism, as well as reductionist and anti-reductionist narratives on the concept. Finally, a few strategies of analysis from published examples of data-mining projects are discussed. Moreover, the examples are re-interpreted in the light of the theoretical material. I argue that the choice of an optimal level of abstraction for the gene is vital for a successful genome analysis.

A Bioinformatics Crash Course for Interpreting Genomics Data

Reductions in genotyping costs and improvements in computational power have made conducting genome-wide association studies (GWAS) standard practice for many complex diseases. GWAS is the assessment of genetic variants across the genome of many individuals to determine which, if any, genetic variants are associated with a specific trait. As with any analysis, there are evolving best practices that should be followed to ensure scientific rigor and reliability in the conclusions. This article presents a brief summary for many of the key bioinformatics considerations when either planning or evaluating GWAS. This review is meant to serve as a guide to those without deep expertise in bioinformatics and GWAS and give them tools to critically evaluate this popular approach to investigating complex diseases. In addition, a checklist is provided that can be used by investigators to evaluate whether a GWAS has appropriately accounted for the many potential sources of bias and generally followed current best practices.

Genome-wide association study identifies zonisamide responsive gene in Parkinson's disease patients

Long-term treatment of Parkinson's disease (PD) by levodopa leads to motor complication "wearing-off". Zonisamide is a nondopaminergic antiparkinsonian drug that can improve "wearing-off" although response to the treatment varies between individuals. To clarify the genetic basis of zonisamide responsiveness, we conducted a genome-wide association study (GWAS) on 200 PD patients from a placebo-controlled clinical trial, including 67 responders whose "off" time decreased ≥1.5 h after 12 weeks of zonisamide treatment and 133 poor responders. We genotyped and evaluated the association between 611,492 single nucleotide polymorphisms (SNPs) and "off" time reduction. We also performed whole-genome imputation, gene- and pathway-based analyses of GWAS data. For promising SNPs, we examined single-tissue expression quantitative trait loci (eQTL) data in the GTEx database. SNP rs16854023 (Mouse double minute 4, MDM4) showed genome-wide significant association with reduced "off" time (P_Adjusted = 4.85 × 10^-9). Carriers of responsive genotype showed >7-fold decrease in mean "off" time compared to noncarriers (1.42 h vs 0.19 h; P = 2.71 × 10^-7). In silico eQTL data indicated that zonisamide sensitivity is associated with higher MDM4 expression. Among the 37 pathways significantly influencing "off" time, calcium and glutamate signaling have also been associated with anti-epileptic effect of zonisamide. MDM4 encodes a negative regulator of p53. The association between improved motor fluctuation and MDM4 upregulation implies that p53 inhibition may prevent dopaminergic neuron loss and consequent motor symptoms. This is the first genome-wide pharmacogenetics study on antiparkinsonian drug. The findings provide a basis for improved management of "wearing-off" in PD by genotype-guided zonisamide treatment.

Fifteen Years of Gene Set Analysis for High-Throughput Genomic Data: A Review of Statistical Approaches and Future Challenges

Over the last decade, gene set analysis has become the first choice for gaining insights into underlying complex biology of diseases through gene expression and gene association studies. It also reduces the complexity of statistical analysis and enhances the explanatory power of the obtained results. Although gene set analysis approaches are extensively used in gene expression and genome wide association data analysis, the statistical structure and steps common to these approaches have not yet been comprehensively discussed, which limits their utility. In this article, we provide a comprehensive overview, statistical structure and steps of gene set analysis approaches used for microarrays, RNA-sequencing and genome wide association data analysis. Further, we also classify the gene set analysis approaches and tools by the type of genomic study, null hypothesis, sampling model and nature of the test statistic, etc. Rather than reviewing the gene set analysis approaches individually, we provide the generation-wise evolution of such approaches for microarrays, RNA-sequencing and genome wide association studies and discuss their relative merits and limitations. Here, we identify the key biological and statistical challenges in current gene set analysis, which will be addressed by statisticians and biologists collectively in order to develop the next generation of gene set analysis approaches. Further, this study will serve as a catalog and provide guidelines to genome researchers and experimental biologists for choosing the proper gene set analysis approach based on several factors.

Bioinformatics services for analyzing massive genomic datasets

The explosive growth of next-generation sequencing data has resulted in ultra-large-scale datasets and ensuing computational problems. In Korea, the amount of genomic data has been increasing rapidly in the recent years. Leveraging these big data requires researchers to use large-scale computational resources and analysis pipelines. A promising solution for addressing this computational challenge is cloud computing, where CPUs, memory, storage, and programs are accessible in the form of virtual machines. Here, we present a cloud computing-based system, Bio-Express, that provides user-friendly, cost-effective analysis of massive genomic datasets. Bio-Express is loaded with predefined multi-omics data analysis pipelines, which are divided into genome, transcriptome, epigenome, and metagenome pipelines. Users can employ predefined pipelines or create a new pipeline for analyzing their own omics data. We also developed several web-based services for facilitating downstream analysis of genome data. Bio-Express web service is freely available at https://www.bioexpress.re.kr/.

Identification of 34 genes conferring genetic and pharmacological risk for the comorbidity of schizophrenia and smoking behaviors

The prevalence of smoking is significantly higher in persons with schizophrenia (SCZ) than in the general population. However, the biological mechanisms of the comorbidity of smoking and SCZ are largely unknown. This study aimed to reveal shared biological pathways for the two diseases by analyzing data from two genome-wide association studies with a total sample size of 153,898. With pathway-based analysis, we first discovered 18 significantly enriched pathways shared by SCZ and smoking, which were classified into five groups: postsynaptic density, cadherin binding, dendritic spine, long-term depression, and axon guidance. Then, by using an integrative analysis of genetic, epigenetic, and expression data, we found not only 34 critical genes (e.g., PRKCZ, ARHGEF3, and CDKN1A) but also various risk-associated SNPs in these genes, which convey susceptibility to the comorbidity of the two disorders. Finally, using both in vivo and in vitro data, we demonstrated that the expression profiles of the 34 genes were significantly altered by multiple psychotropic drugs. Together, this multi-omics study not only reveals target genes for new drugs to treat SCZ but also reveals new insights into the shared genetic vulnerabilities of SCZ and smoking behaviors.

Brain Imaging Genomics: Integrated Analysis and Machine Learning

Brain imaging genomics is an emerging data science field, where integrated analysis of brain imaging and genomics data, often combined with other biomarker, clinical and environmental data, is performed to gain new insights into the phenotypic, genetic and molecular characteristics of the brain as well as their impact on normal and disordered brain function and behavior. It has enormous potential to contribute significantly to biomedical discoveries in brain science. Given the increasingly important role of statistical and machine learning in biomedicine and rapidly growing literature in brain imaging genomics, we provide an up-to-date and comprehensive review of statistical and machine learning methods for brain imaging genomics, as well as a practical discussion on method selection for various biomedical applications.

Gene-set association and epistatic analyses reveal complex gene interaction networks affecting flowering time in a worldwide barley collection

Single-marker genome-wide association studies (GWAS) have successfully detected associations between single nucleotide polymorphisms (SNPs) and agronomic traits such as flowering time and grain yield in barley. However, the analysis of individual SNPs can only account for a small proportion of genetic variation, and can only provide limited knowledge on gene network interactions. Gene-based GWAS approaches provide enormous opportunity both to combine genetic information and to examine interactions among genetic variants. Here, we revisited a previously published phenotypic and genotypic data set of 895 barley varieties grown in two years at four different field locations in Australia. We employed statistical models to examine gene-phenotype associations, as well as two-way epistasis analyses to increase the capability to find novel genes that have significant roles in controlling flowering time in barley. Genetic associations were tested between flowering time and corresponding genotypes of 174 putative flowering time-related genes. Gene-phenotype association analysis detected 113 genes associated with flowering time in barley, demonstrating the unprecedented power of gene-based analysis. Subsequent two-way epistasis analysis revealed 19 pairs of gene×gene interactions involved in controlling flowering time. Our study demonstrates that gene-based association approaches can provide higher capacity for future crop improvement to increase crop performance and adaptation to different environments.

Gene set analysis methods for the functional interpretation of non-mRNA data—Genomic range and ncRNA data

Gene set analysis (GSA) is one of the methods of choice for analyzing the results of current omics studies; however, it has been mainly developed to analyze mRNA (microarray, RNA-Seq) data. The following review includes an update regarding general methods and resources for GSA and then emphasizes GSA methods and tools for non-mRNA omics datasets, specifically genomic range data (ChIP-Seq, SNP and methylation) and ncRNA data (miRNAs, lncRNAs and others). In the end, the state of the GSA field for non-mRNA datasets is discussed, and some current challenges and trends are highlighted, especially the use of network approaches to face complexity issues.

Genome-wide association analysis of hippocampal volume identifies enrichment of neurogenesis-related pathways

Adult neurogenesis occurs in the dentate gyrus of the hippocampus during adulthood and contributes to sustaining the hippocampal formation. To investigate whether neurogenesis-related pathways are associated with hippocampal volume, we performed gene-set enrichment analysis using summary statistics from a large-scale genome-wide association study (N = 13,163) of hippocampal volume from the Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Consortium and two year hippocampal volume changes from baseline in cognitively normal individuals from Alzheimer's Disease Neuroimaging Initiative Cohort (ADNI). Gene-set enrichment analysis of hippocampal volume identified 44 significantly enriched biological pathways (FDR corrected p-value < 0.05), of which 38 pathways were related to neurogenesis-related processes including neurogenesis, generation of new neurons, neuronal development, and neuronal migration and differentiation. For genes highly represented in the significantly enriched neurogenesis-related pathways, gene-based association analysis identified TESC, ACVR1, MSRB3, and DPP4 as significantly associated with hippocampal volume. Furthermore, co-expression network-based functional analysis of gene expression data in the hippocampal subfields, CA1 and CA3, from 32 normal controls showed that distinct co-expression modules were mostly enriched in neurogenesis related pathways. Our results suggest that neurogenesis-related pathways may be enriched for hippocampal volume and that hippocampal volume may serve as a potential phenotype for the investigation of human adult neurogenesis.

Mining Regional Imaging Genetic Associations via Voxel-wise Enrichment Analysis

Brain imaging genetics aims to reveal genetic effects on brain phenotypes, where most studies examine phenotypes defined on anatomical or functional regions of interest (ROIs) given their biologically meaningful annotation and modest dimensionality compared with voxel-wise approaches. Typical ROI-level measures used in these studies are summary statistics from voxel-wise measures in the region, without making full use of individual voxel signals. In this paper, we propose a flexible and powerful framework for mining regional imaging genetic associations via voxel-wise enrichment analysis, which embraces the collective effect of weak voxel-level signals within an ROI. We demonstrate our method on an imaging genetic analysis using data from the Alzheimers Disease Neuroimaging Initiative, where we assess the collective regional genetic effects of voxel-wise FDGPET measures between 116 ROIs and 19 AD candidate SNPs. Compared with traditional ROI-wise and voxel-wise approaches, our method identified 102 additional significant associations, some of which were further supported by evidences in brain tissue-specific expression analysis. This demonstrates the promise of the proposed method as a flexible and powerful framework for exploring imaging genetic effects on the brain.

Efficient pathway enrichment and network analysis of GWAS summary data using GSA-SNP2

Pathway-based analysis in genome-wide association study (GWAS) is being widely used to uncover novel multi-genic functional associations. Many of these pathway-based methods have been used to test the enrichment of the associated genes in the pathways, but exhibited low powers and were highly affected by free parameters. We present the novel method and software GSA-SNP2 for pathway enrichment analysis of GWAS P-value data. GSA-SNP2 provides high power, decent type I error control and fast computation by incorporating the random set model and SNP-count adjusted gene score. In a comparative study using simulated and real GWAS data, GSA-SNP2 exhibited high power and best prioritized gold standard positive pathways compared with six existing enrichment-based methods and two self-contained methods (alternative pathway analysis approach). Based on these results, the difference between pathway analysis approaches was investigated and the effects of the gene correlation structures on the pathway enrichment analysis were also discussed. In addition, GSA-SNP2 is able to visualize protein interaction networks within and across the significant pathways so that the user can prioritize the core subnetworks for further studies. GSA-SNP2 is freely available at https://sourceforge.net/projects/gsasnp2.

Risk of spontaneous preterm birth and fetal growth associates with fetal SLIT2

Spontaneous preterm birth (SPTB) is the leading cause of neonatal death and morbidity worldwide. Both maternal and fetal genetic factors likely contribute to SPTB. We performed a genome-wide association study (GWAS) on a population of Finnish origin that included 247 infants with SPTB (gestational age [GA] < 36 weeks) and 419 term controls (GA 38-41 weeks). The strongest signal came within the gene encoding slit guidance ligand 2 (SLIT2; rs116461311, minor allele frequency 0.05, p = 1.6×10-6). Pathway analysis revealed the top-ranking pathway was axon guidance, which includes SLIT2. In 172 very preterm-born infants (GA <32 weeks), rs116461311 was clearly overrepresented (odds ratio 4.06, p = 1.55×10-7). SLIT2 variants were associated with SPTB in another European population that comprised 260 very preterm infants and 9,630 controls. To gain functional insight, we used immunohistochemistry to visualize SLIT2 and its receptor ROBO1 in placentas from spontaneous preterm and term births. Both SLIT2 and ROBO1 were located in villous and decidual trophoblasts of embryonic origin. Based on qRT-PCR, the mRNA levels of SLIT2 and ROBO1 were higher in the basal plate of SPTB placentas compared to those from term or elective preterm deliveries. In addition, in spontaneous term and preterm births, placental SLIT2 expression was correlated with variations in fetal growth. Knockdown of ROBO1 in trophoblast-derived HTR8/SVneo cells by siRNA indicated that it regulate expression of several pregnancy-specific beta-1-glycoprotein (PSG) genes and genes involved in inflammation. Our results show that the fetal SLIT2 variant and both SLIT2 and ROBO1 expression in placenta and trophoblast cells may be correlated with susceptibility to SPTB. SLIT2-ROBO1 signaling was linked with regulation of genes involved in inflammation, PSG genes, decidualization and fetal growth. We propose that this receptor-ligand couple is a component of the signaling network that promotes SPTB.

Genome-Wide Interaction and Pathway Association Studies for Body Mass Index

Objective: We investigated gene interactions (epistasis) for body mass index (BMI) in a European-American adult female cohort via genome-wide interaction analyses (GWIA) and pathway association analyses. Methods: Genome-wide pairwise interaction analyses were carried out for BMI in 493 extremely obese cases (BMI > 35 kg/m2) and 537 never-overweight controls (BMI < 25 kg/m2). To further validate the results, specific SNPs were selected based on the GWIA results for haplotype-based association studies. Pathway-based association analyses were performed using a modified Gene Set Enrichment Algorithm (GSEA) (GenGen program) to further explore BMI-related pathways using our genome wide association study (GWAS) data set, GIANT, ENGAGE, and DIAGRAM Consortia. Results: The EXOC4-1q23.1 interaction was associated with BMI, with the most significant epistasis between rs7800006 and rs10797020 (P = 2.63 × 10-11). In the pathway-based association analysis, Tob1 pathway showed the most significant association with BMI (empirical P < 0.001, FDR = 0.044, FWER = 0.040). These findings were further validated in different populations. Conclusion: Genome-wide pairwise SNP-SNP interaction and pathway analyses suggest that EXOC4 and TOB1-related pathways may contribute to the development of obesity.

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.

NKG2D gene variation and susceptibility to viral bronchiolitis in childhood

BACKGROUND

Genetic factors associated with bronchiolitis are inadequately characterized. We therefore inspected a selected subpopulation of our previous genome-wide association study (GWAS) of bronchiolitis for overlap with known quantitative trait loci (QTLs) to identify susceptibility loci that potentially affect mRNA and protein levels.

METHODS

GWAS included a Finnish-Swedish case-control population (n = 187), matched for age and site. We integrated GWAS variants (p < 10^-4) with QTL data. We subsequently verified allele-specific expression of identified QTLs by flow cytometry. Association of the resulting candidate loci with bronchiolitis was tested in three additional cohorts from Finland and Denmark (n = 1201).

RESULTS

Bronchiolitis-susceptibility variant rs10772271 resided within QTLs previously associated with NKG2D (NK group 2, member D) mRNA and protein levels. Flow cytometric analysis confirmed the association with protein level in NK cells. The GWAS susceptibility allele (A) of rs10772271 (odds ratio [OR] = 2.34) corresponded with decreased NKG2D expression. The allele was nominally associated with bronchiolitis in one Finnish replicate (OR = 1.50), and the other showed directional consistency (OR = 1.43). No association was detected in Danish population CONCLUSIONS: The bronchiolitis GWAS susceptibility allele was linked to decreased NKG2D expression in the QTL data and in our expression analysis. We propose that reduced NKG2D expression predisposes infants to severe bronchiolitis.

Identification of susceptibility pathways for the role of chromosome 15q25.1 in modifying lung cancer risk

Genome-wide association studies (GWAS) identified the chromosome 15q25.1 locus as a leading susceptibility region for lung cancer. However, the pathogenic pathways, through which susceptibility SNPs within chromosome 15q25.1 affects lung cancer risk, have not been explored. We analyzed three cohorts with GWAS data consisting 42,901 individuals and lung expression quantitative trait loci (eQTL) data on 409 individuals to identify and validate the underlying pathways and to investigate the combined effect of genes from the identified susceptibility pathways. The KEGG neuroactive ligand receptor interaction pathway, two Reactome pathways, and 22 Gene Ontology terms were identified and replicated to be significantly associated with lung cancer risk, with P values less than 0.05 and FDR less than 0.1. Functional annotation of eQTL analysis results showed that the neuroactive ligand receptor interaction pathway and gated channel activity were involved in lung cancer risk. These pathways provide important insights for the etiology of lung cancer.

Gene expression profiles for predicting antibody‑mediated kidney allograft rejection: Analysis of GEO datasets

Antibody‑mediated rejections (AMRs) are one of the most challenging complications that result in the deterioration of renal allograft function and graft loss in a large majority of cases. The purpose of the present study was to characterize a meta‑signature of differentially expressed RNAs associated with AMR in cases of kidney transplantation. Gene Expression Omnibus (GEO) dataset searches up to September 11, 2017, using Medical Subject Heading terms and keywords associated with kidney transplantation, AMR and mRNA arrays were downloaded from the GEO dataset. Using a computational analysis, a meta‑signature was determined that characterized the significant intersection of differentially expressed genes (DEGs). Gene‑set and network analyses were also performed to identify gene sets and sub‑networks associated with the AMR‑related traits. A statistically significant mRNA meta‑signature of upregulated and downregulated gene expression levels that were significantly associated with AMR was identified. C‑X‑C motif chemokine ligand 10 (CXCL10), CXCL9 and guanylate binding protein 1 were the most significantly associated with AMR. DEGs were efficiently identified and were found to be able to predict the occurrence of AMR according to a meta‑analysis approach from publicly available datasets. These methods and results can be applied for a more accurate diagnosis of AMR in transplant cases.

Tissue-specific network-based genome wide study of amygdala imaging phenotypes to identify functional interaction modules

Motivation

Network-based genome-wide association studies (GWAS) aim to identify functional modules from biological networks that are enriched by top GWAS findings. Although gene functions are relevant to tissue context, most existing methods analyze tissue-free networks without reflecting phenotypic specificity.

Results

We propose a novel module identification framework for imaging genetic studies using the tissue-specific functional interaction network. Our method includes three steps: (i) re-prioritize imaging GWAS findings by applying machine learning methods to incorporate network topological information and enhance the connectivity among top genes; (ii) detect densely connected modules based on interactions among top re-prioritized genes; and (iii) identify phenotype-relevant modules enriched by top GWAS findings. We demonstrate our method on the GWAS of [18F]FDG-PET measures in the amygdala region using the imaging genetic data from the Alzheimer's Disease Neuroimaging Initiative, and map the GWAS results onto the amygdala-specific functional interaction network. The proposed network-based GWAS method can effectively detect densely connected modules enriched by top GWAS findings. Tissue-specific functional network can provide precise context to help explore the collective effects of genes with biologically meaningful interactions specific to the studied phenotype.

Availability and implementation

The R code and sample data are freely available at http://www.iu.edu/shenlab/tools/gwasmodule/.

Contact

shenli@iu.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Gene set enrichment analysis and expression pattern exploration implicate an involvement of neurodevelopmental processes in bipolar disorder

BACKGROUND

Bipolar disorder (BD) is a common and highly heritable disorder of mood. Genome-wide association studies (GWAS) have identified several independent susceptibility loci. In order to extract more biological information from GWAS data, multi-locus approaches represent powerful tools since they utilize knowledge about biological processes to integrate functional sets of genes at strongly to moderately associated loci.

METHODS

We conducted gene set enrichment analyses (GSEA) using 2.3 million single-nucleotide polymorphisms, 397 Reactome pathways and 24,025 patients with BD and controls. RNA expression of implicated individual genes and gene sets were examined in post-mortem brains across lifespan.

RESULTS

Two pathways showed a significant enrichment after correction for multiple comparisons in the GSEA: GRB2 events in ERBB2 signaling, for which 6 of 21 genes were BD associated (PFDR = 0.0377), and NCAM signaling for neurite out-growth, for which 11 out of 62 genes were BD associated (PFDR = 0.0451). Most pathway genes showed peaks of RNA co-expression during fetal development and infancy and mapped to neocortical areas and parts of the limbic system.

LIMITATIONS

Pathway associations were technically reproduced by two methods, although they were not formally replicated in independent samples. Gene expression was explored in controls but not in patients.

CONCLUSIONS

Pathway analysis in large GWAS data of BD and follow-up of gene expression patterns in healthy brains provide support for an involvement of neurodevelopmental processes in the etiology of this neuropsychiatric disease. Future studies are required to further evaluate the relevance of the implicated genes on pathway functioning and clinical aspects of BD.

Shared activity patterns arising at genetic susceptibility loci reveal underlying genomic and cellular architecture of human disease

Genetic variants underlying complex traits, including disease susceptibility, are enriched within the transcriptional regulatory elements, promoters and enhancers. There is emerging evidence that regulatory elements associated with particular traits or diseases share similar patterns of transcriptional activity. Accordingly, shared transcriptional activity (coexpression) may help prioritise loci associated with a given trait, and help to identify underlying biological processes. Using cap analysis of gene expression (CAGE) profiles of promoter- and enhancer-derived RNAs across 1824 human samples, we have analysed coexpression of RNAs originating from trait-associated regulatory regions using a novel quantitative method (network density analysis; NDA). For most traits studied, phenotype-associated variants in regulatory regions were linked to tightly-coexpressed networks that are likely to share important functional characteristics. Coexpression provides a new signal, independent of phenotype association, to enable fine mapping of causative variants. The NDA coexpression approach identifies new genetic variants associated with specific traits, including an association between the regulation of the OCT1 cation transporter and genetic variants underlying circulating cholesterol levels. NDA strongly implicates particular cell types and tissues in disease pathogenesis. For example, distinct groupings of disease-associated regulatory regions implicate two distinct biological processes in the pathogenesis of ulcerative colitis; a further two separate processes are implicated in Crohn's disease. Thus, our functional analysis of genetic predisposition to disease defines new distinct disease endotypes. We predict that patients with a preponderance of susceptibility variants in each group are likely to respond differently to pharmacological therapy. Together, these findings enable a deeper biological understanding of the causal basis of complex traits.

Genome-wide analysis of differential RNA editing in epilepsy

The recoding of genetic information through RNA editing contributes to proteomic diversity, but the extent and significance of RNA editing in disease is poorly understood. In particular, few studies have investigated the relationship between RNA editing and disease at a genome-wide level. Here, we developed a framework for the genome-wide detection of RNA sites that are differentially edited in disease. Using RNA-sequencing data from 100 hippocampi from mice with epilepsy (pilocarpine–temporal lobe epilepsy model) and 100 healthy control hippocampi, we identified 256 RNA sites (overlapping with 87 genes) that were significantly differentially edited between epileptic cases and controls. The degree of differential RNA editing in epileptic mice correlated with frequency of seizures, and the set of genes differentially RNA-edited between case and control mice were enriched for functional terms highly relevant to epilepsy, including “neuron projection” and “seizures.” Genes with differential RNA editing were preferentially enriched for genes with a genetic association to epilepsy. Indeed, we found that they are significantly enriched for genes that harbor nonsynonymous de novo mutations in patients with epileptic encephalopathy and for common susceptibility variants associated with generalized epilepsy. These analyses reveal a functional convergence between genes that are differentially RNA-edited in acquired symptomatic epilepsy and those that contribute risk for genetic epilepsy. Taken together, our results suggest a potential role for RNA editing in the epileptic hippocampus in the occurrence and severity of epileptic seizures.

Real-life helping behaviours in North America: A genome-wide association approach

In humans, prosocial behaviour is essential for social functioning. Twin studies suggest this distinct human trait to be partly hardwired. In the last decade research on the genetics of prosocial behaviour focused on neurotransmitters and neuropeptides, such as oxytocin, dopamine, and their respective pathways. Recent trends towards large scale medical studies targeting the genetic basis of complex diseases such as Alzheimer’s disease and schizophrenia pave the way for new directions also in behavioural genetics.

Based on data from 10,713 participants of the American Health and Retirement Study we estimated heritability of helping behaviour–its total variance explained by 1.2 million single nucleotide polymorphisms–to be 11%. Both, fixed models and mixed linear models identified rs11697300, an intergene variant on chromosome 20, as a candidate variant moderating this particular helping behaviour. We assume that this so far undescribed area is worth further investigation in association with human prosocial behaviour.

Multivariate association between single-nucleotide polymorphisms in Alzgene linkage regions and structural changes in the brain: discovery, refinement and validation

Using publicly-available data from the Alzheimer’s Disease Neuroimaging Initiative, we investigate the joint association between single-nucleotide polymorphisms (SNPs) in previously established linkage regions for Alzheimer’s disease (AD) and rates of decline in brain structure. In an initial, discovery stage of analysis, we applied a weighted

Genome-wide association study of language performance in Alzheimer's disease

Language impairment is common in prodromal stages of Alzheimer's disease (AD) and progresses over time. However, the genetic architecture underlying language performance is poorly understood. To identify novel genetic variants associated with language performance, we analyzed brain MRI and performed a genome-wide association study (GWAS) using a composite measure of language performance from the Alzheimer's Disease Neuroimaging Initiative (ADNI; n=1560). The language composite score was associated with brain atrophy on MRI in language and semantic areas. GWAS identified GLI3 (GLI family zinc finger 3) as significantly associated with language performance (p<5×10^-8). Enrichment of GWAS association was identified in pathways related to nervous system development and glutamate receptor function and trafficking. Our results, which warrant further investigation in independent and larger cohorts, implicate GLI3, a developmental transcription factor involved in patterning brain structures, as a putative gene associated with language dysfunction in AD.

Influence maximization in time bounded network identifies transcription factors regulating perturbed pathways

Motivation: To understand the dynamic nature of the biological process, it is crucial to identify perturbed pathways in an altered environment and also to infer regulators that trigger the response. Current time-series analysis methods, however, are not powerful enough to identify perturbed pathways and regulators simultaneously. Widely used methods include methods to determine gene sets such as differentially expressed genes or gene clusters and these genes sets need to be further interpreted in terms of biological pathways using other tools. Most pathway analysis methods are not designed for time series data and they do not consider gene-gene influence on the time dimension.

Results: In this article, we propose a novel time-series analysis method TimeTP for determining transcription factors (TFs) regulating pathway perturbation, which narrows the focus to perturbed sub-pathways and utilizes the gene regulatory network and protein–protein interaction network to locate TFs triggering the perturbation. TimeTP first identifies perturbed sub-pathways that propagate the expression changes along the time. Starting points of the perturbed sub-pathways are mapped into the network and the most influential TFs are determined by influence maximization technique. The analysis result is visually summarized in TF-Pathway map in time clock. TimeTP was applied to PIK3CA knock-in dataset and found significant sub-pathways and their regulators relevant to the PIP3 signaling pathway.

Availability and Implementation: TimeTP is implemented in Python and available at http://biohealth.snu.ac.kr/software/TimeTP/.

Supplementary information:Supplementary data are available at Bioinformatics online.

Contact:sunkim.bioinfo@snu.ac.kr

Genome-wide association study of borderline personality disorder reveals genetic overlap with bipolar disorder, major depression and schizophrenia

Borderline personality disorder (BOR) is determined by environmental and genetic factors, and characterized by affective instability and impulsivity, diagnostic symptoms also observed in manic phases of bipolar disorder (BIP). Up to 20% of BIP patients show comorbidity with BOR. This report describes the first case-control genome-wide association study (GWAS) of BOR, performed in one of the largest BOR patient samples worldwide. The focus of our analysis was (i) to detect genes and gene sets involved in BOR and (ii) to investigate the genetic overlap with BIP. As there is considerable genetic overlap between BIP, major depression (MDD) and schizophrenia (SCZ) and a high comorbidity of BOR and MDD, we also analyzed the genetic overlap of BOR with SCZ and MDD. GWAS, gene-based tests and gene-set analyses were performed in 998 BOR patients and 1545 controls. Linkage disequilibrium score regression was used to detect the genetic overlap between BOR and these disorders. Single marker analysis revealed no significant association after correction for multiple testing. Gene-based analysis yielded two significant genes: DPYD (P=4.42 × 10-7) and PKP4 (P=8.67 × 10-7); and gene-set analysis yielded a significant finding for exocytosis (GO:0006887, PFDR=0.019; FDR, false discovery rate). Prior studies have implicated DPYD, PKP4 and exocytosis in BIP and SCZ. The most notable finding of the present study was the genetic overlap of BOR with BIP (rg=0.28 [P=2.99 × 10-3]), SCZ (rg=0.34 [P=4.37 × 10-5]) and MDD (rg=0.57 [P=1.04 × 10-3]). We believe our study is the first to demonstrate that BOR overlaps with BIP, MDD and SCZ on the genetic level. Whether this is confined to transdiagnostic clinical symptoms should be examined in future studies.

NETWORK-BASED GENOME WIDE STUDY OF HIPPOCAMPAL IMAGING PHENOTYPE IN ALZHEIMER'S DISEASE TO IDENTIFY FUNCTIONAL INTERACTION MODULES

Identification of functional modules from biological network is a promising approach to enhance the statistical power of genome-wide association study (GWAS) and improve biological interpretation for complex diseases. The precise functions of genes are highly relevant to tissue context, while a majority of module identification studies are based on tissue-free biological networks that lacks phenotypic specificity. In this study, we propose a module identification method that maps the GWAS results of an imaging phenotype onto the corresponding tissue-specific functional interaction network by applying a machine learning framework. Ridge regression and support vector machine (SVM) models are constructed to re-prioritize GWAS results, followed by exploring hippocampus-relevant modules based on top predictions using GWAS top findings. We also propose a GWAS top-neighbor-based module identification approach and compare it with Ridge and SVM based approaches. Modules conserving both tissue specificity and GWAS discoveries are identified, showing the promise of the proposal method for providing insight into the mechanism of complex diseases.

Improving the detection of pathways in genome-wide association studies by combined effects of SNPs from Linkage Disequilibrium blocks

Genome-wide association studies (GWAS) have successfully identified single variants associated with diseases. To increase the power of GWAS, gene-based and pathway-based tests are commonly employed to detect more risk factors. However, the gene- and pathway-based association tests may be biased towards genes or pathways containing a large number of single-nucleotide polymorphisms (SNPs) with small P-values caused by high linkage disequilibrium (LD) correlations. To address such bias, numerous pathway-based methods have been developed. Here we propose a novel method, DGAT-path, to divide all SNPs assigned to genes in each pathway into LD blocks, and to sum the chi-square statistics of LD blocks for assessing the significance of the pathway by permutation tests. The method was proven robust with the type I error rate >1.6 times lower than other methods. Meanwhile, the method displays a higher power and is not biased by the pathway size. The applications to the GWAS summary statistics for schizophrenia and breast cancer indicate that the detected top pathways contain more genes close to associated SNPs than other methods. As a result, the method identified 17 and 12 significant pathways containing 20 and 21 novel associated genes, respectively for two diseases. The method is available online by http://sparks-lab.org/server/DGAT-path.

Genomic and transcriptomic predictors of response levels to endurance exercise training

Predicting the responsiveness to regular exercise is a topic of great relevance due to its potential role in personalized exercise medicine applications. The present review focuses on cardiorespiratory fitness (commonly measured by maximal oxygen uptake, V̇O2 max ), a trait with wide-ranging impact on health and performance indicators. Gains in V̇O2 max demonstrate large inter-individual variation even in response to standardized exercise training programmes. The estimated ΔVO2 max heritability of 47% suggests that genomic-based predictors alone are insufficient to account for the total trainability variance. Candidate gene and genome-wide linkage studies have not significantly contributed to our understanding of the molecular basis of trainability. A genome-wide association study suggested that V̇O2 max trainability is influenced by multiple genes of small effects, but these findings still await rigorous replication. Valuable evidence, however, has been obtained by combining skeletal muscle transcript abundance profiles with common DNA variants for the prediction of the V̇O2 max response to exercise training. Although the physiological determinants of V̇O2 max measured at a given time are largely enunciated, what is poorly understood are the details of tissue-specific molecular mechanisms that limit V̇O2 max and related signalling pathways in response to exercise training. Bioinformatics explorations based on thousands of variants have been used to interrogate pathways and systems instead of single variants and genes, and the main findings, along with those from exercise experimental studies, have been summarized here in a working model of V̇O2 max trainability.

PINTnet: construction of condition-specific pathway interaction network by computing shortest paths on weighted PPI

Background

Identifying perturbed pathways in a given condition is crucial in understanding biological phenomena. In addition to identifying perturbed pathways individually, pathway analysis should consider interactions among pathways. Currently available pathway interaction prediction methods are based on the existence of overlapping genes between pathways, protein-protein interaction (PPI) or functional similarities. However, these approaches just consider the pathways as a set of genes, thus they do not take account of topological features. In addition, most of the existing approaches do not handle the explicit gene expression quantity information that is routinely measured by RNA-sequecing.

Results

To overcome these technical issues, we developed a new pathway interaction network construction method using PPI, closeness centrality and shortest paths. We tested our approach on three different high-throughput RNA-seq data sets: pregnant mice data to reveal the role of serotonin on beta cell mass, bone-metastatic breast cancer data and autoimmune thyroiditis data to study the role of IFN- α. Our approach successfully identified the pathways reported in the original papers. For the pathways that are not directly mentioned in the original papers, we were able to find evidences of pathway interactions by the literature search. Our method outperformed two existing approaches, overlapping gene-based approach (OGB) and protein-protein interaction-based approach (PB), in experiments with the three data sets.

Conclusion

Our results show that PINTnet successfully identified condition-specific perturbed pathways and the interactions between the pathways. We believe that our method will be very useful in characterizing biological mechanisms at the pathway level. PINTnet is available at http://biohealth.snu.ac.kr/software/PINTnet/.

Rare and common epilepsies converge on a shared gene regulatory network providing opportunities for novel antiepileptic drug discovery

BACKGROUND

The relationship between monogenic and polygenic forms of epilepsy is poorly understood and the extent to which the genetic and acquired epilepsies share common pathways is unclear. Here, we use an integrated systems-level analysis of brain gene expression data to identify molecular networks disrupted in epilepsy.

RESULTS

We identified a co-expression network of 320 genes (M30), which is significantly enriched for non-synonymous de novo mutations ascertained from patients with monogenic epilepsy and for common variants associated with polygenic epilepsy. The genes in the M30 network are expressed widely in the human brain under tight developmental control and encode physically interacting proteins involved in synaptic processes. The most highly connected proteins within the M30 network were preferentially disrupted by deleterious de novo mutations for monogenic epilepsy, in line with the centrality-lethality hypothesis. Analysis of M30 expression revealed consistent downregulation in the epileptic brain in heterogeneous forms of epilepsy including human temporal lobe epilepsy, a mouse model of acquired temporal lobe epilepsy, and a mouse model of monogenic Dravet (SCN1A) disease. These results suggest functional disruption of M30 via gene mutation or altered expression as a convergent mechanism regulating susceptibility to epilepsy broadly. Using the large collection of drug-induced gene expression data from Connectivity Map, several drugs were predicted to preferentially restore the downregulation of M30 in epilepsy toward health, most notably valproic acid, whose effect on M30 expression was replicated in neurons.

CONCLUSIONS

Taken together, our results suggest targeting the expression of M30 as a potential new therapeutic strategy in epilepsy.

snpGeneSets: An R Package for Genome-Wide Study Annotation

Genome-wide studies (GWS) of SNP associations and differential gene expressions have generated abundant results; next-generation sequencing technology has further boosted the number of variants and genes identified. Effective interpretation requires massive annotation and downstream analysis of these genome-wide results, a computationally challenging task. We developed the snpGeneSets package to simplify annotation and analysis of GWS results. Our package integrates local copies of knowledge bases for SNPs, genes, and gene sets, and implements wrapper functions in the R language to enable transparent access to low-level databases for efficient annotation of large genomic data. The package contains functions that execute three types of annotations: (1) genomic mapping annotation for SNPs and genes and functional annotation for gene sets; (2) bidirectional mapping between SNPs and genes, and genes and gene sets; and (3) calculation of gene effect measures from SNP associations and performance of gene set enrichment analyses to identify functional pathways. We applied snpGeneSets to type 2 diabetes (T2D) results from the NHGRI genome-wide association study (GWAS) catalog, a Finnish GWAS, and a genome-wide expression study (GWES). These studies demonstrate the usefulness of snpGeneSets for annotating and performing enrichment analysis of GWS results. The package is open-source, free, and can be downloaded at: https://www.umc.edu/biostats_software/.

Pathway analysis of complex diseases for GWAS, extending to consider rare variants, multi-omics and interactions

BACKGROUND

Genome-wide association studies (GWAS) is a major method for studying the genetics of complex diseases. Finding all sequence variants to explain fully the aetiology of a disease is difficult because of their small effect sizes. To better explain disease mechanisms, pathway analysis is used to consolidate the effects of multiple variants, and hence increase the power of the study. While pathway analysis has previously been performed within GWAS only, it can now be extended to examining rare variants, other "-omics" and interaction data.

SCOPE OF REVIEW

1. Factors to consider in the choice of software for GWAS pathway analysis. 2. Examples of how pathway analysis is used to analyse rare variants, other "-omics" and interaction data.

MAJOR CONCLUSIONS

To choose appropriate software tools, factors for consideration include covariate compatibility, null hypothesis, one- or two-step analysis required, curation method of gene sets, size of pathways, and size of flanking regions to define gene boundaries. For rare variants, analysis performance depends on consistency between assumed and actual effect distribution of variants. Integration of other "-omics" data and interaction can better explain gene functions.

GENERAL SIGNIFICANCE

Pathway analysis methods will be more readily used for integration of multiple sources of data, and enable more accurate prediction of phenotypes.

Improving Gene-Set Enrichment Analysis of RNA-Seq Data with Small Replicates

Deregulated pathways identified from transcriptome data of two sample groups have played a key role in many genomic studies. Gene-set enrichment analysis (GSEA) has been commonly used for pathway or functional analysis of microarray data, and it is also being applied to RNA-seq data. However, most RNA-seq data so far have only small replicates. This enforces to apply the gene-permuting GSEA method (or preranked GSEA) which results in a great number of false positives due to the inter-gene correlation in each gene-set. We demonstrate that incorporating the absolute gene statistic in one-tailed GSEA considerably improves the false-positive control and the overall discriminatory ability of the gene-permuting GSEA methods for RNA-seq data. To test the performance, a simulation method to generate correlated read counts within a gene-set was newly developed, and a dozen of currently available RNA-seq enrichment analysis methods were compared, where the proposed methods outperformed others that do not account for the inter-gene correlation. Analysis of real RNA-seq data also supported the proposed methods in terms of false positive control, ranks of true positives and biological relevance. An efficient R package (AbsFilterGSEA) coded with C++ (Rcpp) is available from CRAN.

Genetic Influences on Plasma Homocysteine Levels in African Americans and Yoruba Nigerians

Plasma homocysteine, a metabolite involved in key cellular methylation processes seems to be implicated in cognitive functions and cardiovascular health with its high levels representing a potential modifiable risk factor for Alzheimer's disease (AD) and other dementias. A better understanding of the genetic factors regulating homocysteine levels, particularly in non-white populations, may help in risk stratification analyses of existing clinical trials and may point to novel targets for homocysteine-lowering therapy. To identify genetic influences on plasma homocysteine levels in individuals with African ancestry, we performed a targeted gene and pathway-based analysis using a priori biological information and then to identify new association performed a genome-wide association study. All analyses used combined data from the African American and Yoruba cohorts from the Indianapolis-Ibadan Dementia Project. Targeted analyses demonstrated significant associations of homocysteine and variants within the CBS (Cystathionine beta-Synthase) gene. We identified a novel genome-wide significant association of the AD risk gene CD2AP (CD2-associated protein) with plasma homocysteine levels in both cohorts. Minor allele (T) carriers of identified CD2AP variant (rs6940729) exhibited decreased homocysteine level. Pathway enrichment analysis identified several interesting pathways including the GABA receptor activation pathway. This is noteworthy given the known antagonistic effect of homocysteine on GABA receptors. These findings identify several new targets warranting further investigation in relation to the role of homocysteine in neurodegeneration.

Genome-Wide Association of CKD Progression: The Chronic Renal Insufficiency Cohort Study

The rate of decline of renal function varies significantly among individuals with CKD. To understand better the contribution of genetics to CKD progression, we performed a genome-wide association study among participants in the Chronic Renal Insufficiency Cohort Study. Our outcome of interest was CKD progression measured as change in eGFR over time among 1331 blacks and 1476 whites with CKD. We stratified all analyses by race and subsequently, diabetes status. Single-nucleotide polymorphisms (SNPs) that surpassed a significance threshold of P<1×10^-6 for association with eGFR slope were selected as candidates for follow-up and secondarily tested for association with proteinuria and time to ESRD. We identified 12 such SNPs among black patients and six such SNPs among white patients. We were able to conduct follow-up analyses of three candidate SNPs in similar (replication) cohorts and eight candidate SNPs in phenotype-related (validation) cohorts. Among blacks without diabetes, rs653747 in LINC00923 replicated in the African American Study of Kidney Disease and Hypertension cohort (discovery P=5.42×10^-7; replication P=0.039; combined P=7.42×10^-9). This SNP also associated with ESRD (hazard ratio, 2.0 (95% confidence interval, 1.5 to 2.7); P=4.90×10^-6). Similarly, rs931891 in LINC00923 associated with eGFR decline (P=1.44×10^-4) in white patients without diabetes. In summary, SNPs in LINC00923, an RNA gene expressed in the kidney, significantly associated with CKD progression in individuals with nondiabetic CKD. However, the lack of equivalent cohorts hampered replication for most discovery loci. Further replication of our findings in comparable study populations is warranted.

Pathway analysis of bladder cancer genome-wide association study identifies novel pathways involved in bladder cancer development

Genome-wide association studies (GWAS) are designed to identify individual regions associated with cancer risk, but only explain a small fraction of the inherited variability. Alternative approach analyzing genetic variants within biological pathways has been proposed to discover networks of susceptibility genes with additional effects. The gene set enrichment analysis (GSEA) may complement and expand traditional GWAS analysis to identify novel genes and pathways associated with bladder cancer risk. We selected three GSEA methods: Gen-Gen, Aligator, and the SNP Ratio Test to evaluate cellular signaling pathways involved in bladder cancer susceptibility in a Texas GWAS population. The candidate genetic polymorphisms from the significant pathway selected by GSEA were validated in an independent NCI GWAS. We identified 18 novel pathways (P < 0.05) significantly associated with bladder cancer risk. Five of the most promising pathways (P ≤ 0.001 in any of the three GSEA methods) among the 18 pathways included two cell cycle pathways and neural cell adhesion molecule (NCAM), platelet-derived growth factor (PDGF), and unfolded protein response pathways. We validated the candidate polymorphisms in the NCI GWAS and found variants of RAPGEF1, SKP1, HERPUD1, CACNB2, CACNA1C, CACNA1S, COL4A2, SRC, and CACNA1C were associated with bladder cancer risk. Two CCNE1 variants, rs8102137 and rs997669, from cell cycle pathways showed the strongest associations; the CCNE1 signal at 19q12 has already been reported in previous GWAS. These findings offer additional etiologic insights highlighting the specific genes and pathways associated with bladder cancer development. GSEA may be a complementary tool to GWAS to identify additional loci of cancer susceptibility.

Genetics of Coronary Artery Disease

Genetic factors contribute importantly to the risk of coronary artery disease (CAD), and in the past decade, there has been major progress in this area. The tools applied include genome-wide association studies encompassing >200,000 individuals complemented by bioinformatic approaches, including 1000 Genomes imputation, expression quantitative trait locus analyses, and interrogation of Encyclopedia of DNA Elements, Roadmap, and other data sets. close to 60 common SNPs (minor allele frequency>0.05) associated with CAD risk and reaching genome-wide significance (P<5 × 10(-8)) have been identified. Furthermore, a total of 202 independent signals in 109 loci have achieved a false discovery rate (q<0.05) and together explain 28% of the estimated heritability of CAD. These data have been used successfully to create genetic risk scores that can improve risk prediction beyond conventional risk factors and identify those individuals who will benefit most from statin therapy. Such information also has important applications in clinical medicine and drug discovery by using a Mendelian randomization approach to interrogate the causal nature of many factors found to associate with CAD risk in epidemiological studies. In contrast to genome-wide association studies, whole-exome sequencing has provided valuable information directly relevant to genes with known roles in plasma lipoprotein metabolism but has, thus far, failed to identify other rare coding variants linked to CAD. Overall, recent studies have led to a broader understanding of the genetic architecture of CAD and demonstrate that it largely derives from the cumulative effect of multiple common risk alleles individually of small effect size rather than rare variants with large effects on CAD risk. Despite this success, there has been limited progress in understanding the function of the novel loci; the majority of which are in noncoding regions of the genome.

Integrated pathway analysis of nasopharyngeal carcinoma implicates the axonemal dynein complex in the Malaysian cohort

Nasopharyngeal carcinoma (NPC) is an epithelial squamous cell carcinoma on the mucosal lining of the nasopharynx. The etiology of NPC remains elusive despite many reported studies. Most studies employ a single platform approach, neglecting the cumulative influence of both the genome and transcriptome toward NPC development. We aim to employ an integrated pathway approach to identify dysregulated pathways linked to NPC. Our approach combines imputation NPC GWAS data from a Malaysian cohort as well as published expression data GSE12452 from both NPC and non-NPC nasopharynx tissues. Pathway association for GWAS data was performed using MAGENTA while for expression data, GSA-SNP was used with gene p values derived from differential expression values from GEO2R. Our study identified NPC association in the gene ontology (GO) axonemal dynein complex pathway (pGWAS-GSEA  = 1.98 × 10(-2) ; pExpr-GSEA  = 1.27 × 10(-24) ; pBonf-Combined  = 4.15 × 10(-21) ). This association was replicated in a separate cohort using gene expression data from NPC and non-NPC nasopharynx tissues (pAmpliSeq-GSEA  = 6.56 × 10(-4) ). Loss of function in the axonemal dynein complex causes impaired cilia function, leading to poor mucociliary clearance and subsequently upper or lower respiratory tract infection, the former of which includes the nasopharynx. Our approach illustrates the potential use of integrated pathway analysis in detecting gene sets involved in the development of NPC in the Malaysian cohort.

Two-dimensional enrichment analysis for mining high-level imaging genetic associations

Enrichment analysis has been widely applied in the genome-wide association studies, where gene sets corresponding to biological pathways are examined for significant associations with a phenotype to help increase statistical power and improve biological interpretation. In this work, we expand the scope of enrichment analysis into brain imaging genetics, an emerging field that studies how genetic variation influences brain structure and function measured by neuroimaging quantitative traits (QT). Given the high dimensionality of both imaging and genetic data, we propose to study Imaging Genetic Enrichment Analysis (IGEA), a new enrichment analysis paradigm that jointly considers meaningful gene sets (GS) and brain circuits (BC) and examines whether any given GS–BC pair is enriched in a list of gene–QT findings. Using gene expression data from Allen Human Brain Atlas and imaging genetics data from Alzheimer’s Disease Neuroimaging Initiative as test beds, we present an IGEA framework and conduct a proof-of-concept study. This empirical study identifies 25 significant high-level two-dimensional imaging genetics modules. Many of these modules are relevant to a variety of neurobiological pathways or neurodegenerative diseases, showing the promise of the proposal framework for providing insight into the mechanism of complex diseases.

Genome-wide association study of antidepressant response: involvement of the inorganic cation transmembrane transporter activity pathway

BACKGROUND

Genome-wide association studies (GWAS) represent the current frontier in pharmacogenomics. Thousands of subjects of Caucasian ancestry have been included in previous GWAS investigating antidepressant response. GWAS focused on this phenotype are lacking in Asian populations.

METHODS

A sample of 109 major depressive disorder (MDD) patients of Korean origin in antidepressant treatment was collected. Phenotypes were response and remission according to the Hamilton Rating Scale for Depression (HRSD). Genome-wide genotyping was performed using the Illumina Human Omni2.5-8 platform. The same phenotypes were used in the STAR*D level 1 (n = 1677) for independent replication. In order to corroborate findings and increase the comparability between the two datasets, three levels of analysis (SNPs, genes and pathways) were carried out. Bonferroni correction, permutations, and replication across samples were used to reduce the risk of false positives.

RESULTS

Among the genes replicated across the two samples (permutated p < 0.05 in both of them), CTNNA3 appeared promising. The inorganic cation transmembrane transporter activity pathway (GO:0022890) was associated with antidepressant response in both samples (p = 2.9e-5 and p = 0.001 in the Korean and STAR*D samples, respectively) and this pathway included CACNA1A, CACNA1C, and CACNB2 genes.

CONCLUSIONS

The present study supported the involvement of genes coding for subunits of L-type voltage-gated calcium channel in antidepressant efficacy across different ethnicities but replication of findings is required before any definitive statement.