Genetic analysis of biological pathway data through genomic randomization

Statistical driver genes as a means to uncover missing heritability for age-related macular degeneration

Background

Age-related macular degeneration (AMD) is a progressive retinal disease contributing to blindness worldwide. Multiple estimates for AMD heritability (h²) exist; however, a substantial proportion of h² is not attributable to known genomic loci. The International AMD Genomics Consortium (IAMDGC) gathered the largest dataset of advanced AMD (ADV) cases and controls available and identified 34 loci containing 52 independent risk variants defining known AMD h². To better define AMD heterogeneity, we used Pathway Analysis by Randomization Incorporating Structure (PARIS) on the IAMDGC data and identified 8 statistical driver genes (SDGs), including 2 novel SDGs not discovered by the IAMDGC. We chose to further investigate these pathway-based risk genes and determine their contribution to ADV h², as well as the differential ADV subtype h².

Methods

We performed genomic-relatedness-based restricted maximum-likelihood (GREML) analyses on ADV, geographic atrophy (GA), and choroidal neovascularization (CNV) subtypes to investigate the h² of genotyped variants on the full DNA array chip, 34 risk loci (n = 2758 common variants), 52 variants from the IAMDGC 2016 GWAS, and the 8 SDGs, specifically the novel 2 SDGs, PPARA and PLCG2.

Results

Via GREML, full chip h² was 44.05% for ADV, 46.37% for GA, and 62.03% for CNV. The lead 52 variants’ h² (ADV: 14.52%, GA: 8.02%, CNV: 13.62%) and 34 loci h² (ADV: 13.73%, GA: 8.81%, CNV: 12.89%) indicate that known variants contribute ~ 14% to ADV h². SDG variants account for a small percentage of ADV, GA, and CNV heritability, but estimates based on the combination of SDGs and the 34 known loci are similar to those calculated for known loci alone. We identified modest epistatic interactions among variants in the 2 SDGs and the 52 IAMDGC variants, including modest interactions between variants in PPARA and PLCG2.

Conclusions

Pathway analyses, which leverage biological relationships among genes in a pathway, may be useful in identifying additional loci that contribute to the heritability of complex disorders in a non-additive manner. Heritability analyses of these loci, especially amongst disease subtypes, may provide clues to the importance of specific genes to the genetic architecture of AMD.

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.

Pathway Analysis Integrating Genome-Wide and Functional Data Identifies PLCG2 as a Candidate Gene for Age-Related Macular Degeneration

Purpose

Age-related macular degeneration (AMD) is the worldwide leading cause of blindness among the elderly. Although genome-wide association studies (GWAS) have identified AMD risk variants, their roles in disease etiology are not well-characterized, and they only explain a portion of AMD heritability.

Methods

We performed pathway analyses using summary statistics from the International AMD Genomics Consortium's 2016 GWAS and multiple pathway databases to identify biological pathways wherein genetic association signals for AMD may be aggregating. We determined which genes contributed most to significant pathway signals across the databases. We characterized these genes by constructing protein-protein interaction networks and performing motif analysis.

Results

We determined that eight genes (C2, C3, LIPC, MICA, NOTCH4, PLCG2, PPARA, and RAD51B) "drive" the statistical signals observed across pathways curated in the Kyoto Encyclopedia of Genes and Genomes (KEGG), Reactome, and Gene Ontology (GO) databases. We further refined our definition of statistical driver gene to identify PLCG2 as a candidate gene for AMD due to its significant gene-level signals (P < 0.0001) across KEGG, Reactome, GO, and NetPath pathways.

Conclusions

We performed pathway analyses on the largest available collection of advanced AMD cases and controls in the world. Eight genes strongly contributed to significant pathways from the three larger databases, and one gene (PLCG2) was central to significant pathways from all four databases. This is, to our knowledge, the first study to identify PLCG2 as a candidate gene for AMD based solely on genetic burden. Our findings reinforce the utility of integrating in silico genetic and biological pathway data to investigate the genetic architecture of AMD.

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 Review on Methods for Detecting SNP Interactions in High-Dimensional Genomic Data

In this era of genome-wide association studies (GWAS), the quest for understanding the genetic architecture of complex diseases is rapidly increasing more than ever before. The development of high throughput genotyping and next generation sequencing technologies enables genetic epidemiological analysis of large scale data. These advances have led to the identification of a number of single nucleotide polymorphisms (SNPs) responsible for disease susceptibility. The interactions between SNPs associated with complex diseases are increasingly being explored in the current literature. These interaction studies are mathematically challenging and computationally complex. These challenges have been addressed by a number of data mining and machine learning approaches. This paper reviews the current methods and the related software packages to detect the SNP interactions that contribute to diseases. The issues that need to be considered when developing these models are addressed in this review. The paper also reviews the achievements in data simulation to evaluate the performance of these models. Further, it discusses the future of SNP interaction analysis.

Testosterone Pathway Genetic Polymorphisms in Relation to Primary Open-Angle Glaucoma: An Analysis in Two Large Datasets

Purpose

Sex hormones may be associated with primary open-angle glaucoma (POAG), although the mechanisms are unclear. We previously observed that gene variants involved with estrogen metabolism were collectively associated with POAG in women but not men; here we assessed gene variants related to testosterone metabolism collectively and POAG risk.

Methods

We used two datasets: one from the United States (3853 cases and 33,480 controls) and another from Australia (1155 cases and 1992 controls). Both datasets contained densely called genotypes imputed to the 1000 Genomes reference panel. We used pathway- and gene-based approaches with Pathway Analysis by Randomization Incorporating Structure (PARIS) software to assess the overall association between a panel of single nucleotide polymorphisms (SNPs) in testosterone metabolism genes and POAG. In sex-stratified analyses, we evaluated POAG overall and POAG subtypes defined by maximum IOP (high-tension [HTG] or normal tension glaucoma [NTG]).

Results

In the US dataset, the SNP panel was not associated with POAG (permuted P = 0.77), although there was an association in the Australian sample (permuted P = 0.018). In both datasets, the SNP panel was associated with POAG in men (permuted P ≤ 0.033) and not women (permuted P ≥ 0.42), but in gene-based analyses, there was no consistency on the main genes responsible for these findings. In both datasets, the testosterone pathway association with HTG was significant (permuted P ≤ 0.011), but again, gene-based analyses showed no consistent driver gene associations.

Conclusions

Collectively, testosterone metabolism pathway SNPs were consistently associated with the high-tension subtype of POAG in two datasets.

Genome-wide association and pathway analysis of left ventricular function after anthracycline exposure in adults

BACKGROUND

Anthracyclines are important chemotherapeutic agents, but their use is limited by cardiotoxicity. Candidate gene and genome-wide studies have identified putative risk loci for overt cardiotoxicity and heart failure, but there has been no comprehensive assessment of genomic variation influencing the intermediate phenotype of anthracycline-related changes in left ventricular (LV) function. The purpose of this study was to identify genetic factors influencing changes in LV function after anthracycline chemotherapy.

METHODS

We conducted a genome-wide association study (GWAS) of change in LV function after anthracycline exposure in 385 patients identified from BioVU, a resource linking DNA samples to de-identified electronic medical record data. Variants with P values less than 1×10 were independently tested for replication in a cohort of 181 anthracycline-exposed patients from a prospective clinical trial. Pathway analysis was performed to assess combined effects of multiple genetic variants.

RESULTS

Both cohorts were middle-aged adults of predominantly European descent. Among 11 candidate loci identified in discovery GWAS, one single nucleotide polymorphism near PR domain containing 2, with ZNF domain (PRDM2), rs7542939, had a combined P value of 6.5×10 in meta-analysis. Eighteen Kyoto Encyclopedia of Gene and Genomes pathways showed strong enrichment for variants associated with the primary outcome. Identified pathways related to DNA repair, cellular metabolism, and cardiac remodeling.

CONCLUSION

Using genome-wide association we identified a novel candidate susceptibility locus near PRDM2. Variation in genes belonging to pathways related to DNA repair, metabolism, and cardiac remodeling may influence changes in LV function after anthracycline exposure.

Multiple analyses of large-scale genome-wide association study highlight new risk pathways in lumbar spine bone mineral density

Osteoporosis is a common human complex disease. It is mainly characterized by low bone mineral density (BMD) and low-trauma osteoporotic fractures (OF). Until now, a large proportion of heritability has yet to be explained. The existing large-scale genome-wide association studies (GWAS) provide strong support for the investigation of osteoporosis mechanisms using pathway analysis. Recent findings showed that different risk pathways may be involved in BMD in different tissues. Here, we conducted multiple pathway analyses of a large-scale lumbar spine BMD GWAS dataset (2,468,080 SNPs and 31,800 samples) using two published gene-based analysis software including ProxyGeneLD and the PLINK. Using BMD genes from ProxyGeneLD, we identified 51 significant KEGG pathways with adjusted P<0.01. Using BMD genes from PLINK, we identified 38 significant KEGG pathways with adjusted P<0.01. Interestingly, 33 pathways are shared in both methods. In summary, we not only identified the known risk pathway such as Wnt signaling, in which the top GWAS variants are significantly enriched, but also highlight some new risk pathways. Interestingly, evidence from further supports the involvement of these pathways in MBD.

Bio-collections in autism research

Autism spectrum disorder (ASD) is a group of complex neurodevelopmental disorders with diverse clinical manifestations and symptoms. In the last 10 years, there have been significant advances in understanding the genetic basis for ASD, critically supported through the establishment of ASD bio-collections and application in research. Here, we summarise a selection of major ASD bio-collections and their associated findings. Collectively, these include mapping ASD candidate genes, assessing the nature and frequency of gene mutations and their association with ASD clinical subgroups, insights into related molecular pathways such as the synapses, chromatin remodelling, transcription and ASD-related brain regions. We also briefly review emerging studies on the use of induced pluripotent stem cells (iPSCs) to potentially model ASD in culture. These provide deeper insight into ASD progression during development and could generate human cell models for drug screening. Finally, we provide perspectives concerning the utilities of ASD bio-collections and limitations, and highlight considerations in setting up a new bio-collection for ASD research.

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.

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.

Common variants at PVT1, ATG13-AMBRA1, AHI1 and CLEC16A are associated with selective IgA deficiency

Selective immunoglobulin A deficiency (IgAD) is the most common primary immunodeficiency in Europeans. Our genome-wide association study (GWAS) meta-analysis of 1,635 patients with IgAD and 4,852 controls identified four new significant (P < 5 × 10-8) loci and association with a rare IFIH1 variant (p.Ile923Val). Peak new variants (PVT1, P = 4.3 × 10-11; ATG13-AMBRA1, P = 6.7 × 10-10; AHI1, P = 8.4 × 10-10; CLEC16A, P = 1.4 × 10-9) overlapped with autoimmune markers (3/4) and correlated with 21 putative regulatory variants, including expression quantitative trait loci (eQTLs) for AHI1 and DEXI and DNase hypersensitivity sites in FOXP3+ regulatory T cells. Pathway analysis of the meta-analysis results showed striking association with the KEGG pathway for IgA production (pathway P < 0.0001), with 22 of the 30 annotated pathway genes containing at least one variant with P ≤ 0.05 in the IgAD meta-analysis. These data suggest that a complex network of genetic effects, including genes known to influence the biology of IgA production, contributes to IgAD.

Shared Genetic Etiology of Autoimmune Diseases in Patients from a Biorepository Linked to De-identified Electronic Health Records

Autoimmune diseases represent a significant medical burden affecting up to 5–8% of the U.S. population. While genetics is known to play a role, studies of common autoimmune diseases are complicated by phenotype heterogeneity, limited sample sizes, and a single disease approach. Here we performed a targeted genetic association study for cases of multiple sclerosis (MS), rheumatoid arthritis (RA), and Crohn's disease (CD) to assess which common genetic variants contribute individually and pleiotropically to disease risk. Joint modeling and pathway analysis combining the three phenotypes were performed to identify common underlying mechanisms of risk of autoimmune conditions. European American cases of MS, RA, and CD, (n = 119, 53, and 129, respectively) and 1924 controls were identified using de-identified electronic health records (EHRs) through a combination of International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) billing codes, Current Procedural Terminology (CPT) codes, medication lists, and text matching. As expected, hallmark SNPs in MS, such as DQA1 rs9271366 (OR = 1.91; p = 0.008), replicated in the present study. Both MS and CD were associated with TIMMDC1 rs2293370 (OR = 0.27, p = 0.01; OR = 0.25, p = 0.02; respectively). Additionally, PDE2A rs3781913 was significantly associated with both CD and RA (OR = 0.46, p = 0.02; OR = 0.32, p = 0.02; respectively). Joint modeling and pathway analysis identified variants within the KEGG NOD-like receptor signaling pathway and Shigellosis pathway as being correlated with the combined autoimmune phenotype. Our study replicated previously-reported genetic associations for MS and CD in a population derived from de-identified EHRs. We found evidence to support a shared genetic etiology between CD/MS and CD/RA outside of the major histocompatibility complex region and identified KEGG pathways indicative of a bacterial pathogenesis risk for autoimmunity in a joint model. Future work to elucidate this shared etiology will be key in the development of risk models as envisioned in the era of precision medicine.

Assessing the Association of Mitochondrial Genetic Variation With Primary Open-Angle Glaucoma Using Gene-Set Analyses

PURPOSE

Recent studies indicate that mitochondrial proteins may contribute to the pathogenesis of primary open-angle glaucoma (POAG). In this study, we examined the association between POAG and common variations in gene-encoding mitochondrial proteins.

METHODS

We examined genetic data from 3430 POAG cases and 3108 controls derived from the combination of the GLAUGEN and NEIGHBOR studies. We constructed biological-system coherent mitochondrial nuclear-encoded protein gene-sets by intersecting the MitoCarta database with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We examined the mitochondrial gene-sets for association with POAG and with normal-tension glaucoma (NTG) and high-tension glaucoma (HTG) subsets using Pathway Analysis by Randomization Incorporating Structure.

RESULTS

We identified 22 KEGG pathways with significant mitochondrial protein-encoding gene enrichment, belonging to six general biological classes. Among the pathway classes, mitochondrial lipid metabolism was associated with POAG overall (P = 0.013) and with NTG (P = 0.0006), and mitochondrial carbohydrate metabolism was associated with NTG (P = 0.030). Examining the individual KEGG pathway mitochondrial gene-sets, fatty acid elongation and synthesis and degradation of ketone bodies, both lipid metabolism pathways, were significantly associated with POAG (P = 0.005 and P = 0.002, respectively) and NTG (P = 0.0004 and P < 0.0001, respectively). Butanoate metabolism, a carbohydrate metabolism pathway, was significantly associated with POAG (P = 0.004), NTG (P = 0.001), and HTG (P = 0.010).

CONCLUSIONS

We present an effective approach for assessing the contributions of mitochondrial genetic variation to open-angle glaucoma. Our findings support a role for mitochondria in POAG pathogenesis and specifically point to lipid and carbohydrate metabolism pathways as being important.

Genetic Association Analysis of Drusen Progression

PURPOSE

Age-related macular degeneration is a common form of vision loss affecting older adults. The etiology of AMD is multifactorial and is influenced by environmental and genetic risk factors. In this study, we examine how 19 common risk variants contribute to drusen progression, a hallmark of AMD pathogenesis.

METHODS

Exome chip data was made available through the International AMD Genomics Consortium (IAMDGC). Drusen quantification was carried out with color fundus photographs using an automated drusen detection and quantification algorithm. A genetic risk score (GRS) was calculated per subject by summing risk allele counts at 19 common genetic risk variants weighted by their respective effect sizes. Pathway analysis of drusen progression was carried out with the software package Pathway Analysis by Randomization Incorporating Structure.

RESULTS

We observed significant correlation with drusen baseline area and the GRS in the age-related eye disease study (AREDS) dataset (ρ = 0.175, P = 0.006). Measures of association were not statistically significant between drusen progression and the GRS (P = 0.54). Pathway analysis revealed the cell adhesion molecules pathway as the most highly significant pathway associated with drusen progression (corrected P = 0.02).

CONCLUSIONS

In this study, we explored the potential influence of known common AMD genetic risk factors on drusen progression. Our results from the GRS analysis showed association of increasing genetic burden (from 19 AMD associated loci) to baseline drusen load but not drusen progression in the AREDS dataset while pathway analysis suggests additional genetic contributors to AMD risk.

Pathway analysis by randomization incorporating structure-PARIS: an update

MOTIVATION

We present an update to the pathway enrichment analysis tool 'Pathway Analysis by Randomization Incorporating Structure (PARIS)' that determines aggregated association signals generated from genome-wide association study results. Pathway-based analyses highlight biological pathways associated with phenotypes. PARIS uses a unique permutation strategy to evaluate the genomic structure of interrogated pathways, through permutation testing of genomic features, thus eliminating many of the over-testing concerns arising with other pathway analysis approaches.

RESULTS

We have updated PARIS to incorporate expanded pathway definitions through the incorporation of new expert knowledge from multiple database sources, through customized user provided pathways, and other improvements in user flexibility and functionality.

AVAILABILITY AND IMPLEMENTATION

PARIS is freely available to all users at https://ritchielab.psu.edu/software/paris-download

CONTACT

jnc43@case.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The genetics of human autoimmune disease: A perspective on progress in the field and future directions

Progress in defining the genetics of autoimmune disease has been dramatically enhanced by large scale genetic studies. Genome-wide approaches, examining hundreds or for some diseases thousands of cases and controls, have been implemented using high throughput genotyping and appropriate algorithms to provide a wealth of data over the last decade. These studies have identified hundreds of non-HLA loci as well as further defining HLA variations that predispose to different autoimmune diseases. These studies to identify genetic risk loci are also complemented by progress in gene expression studies including definition of expression quantitative trait loci (eQTL), various alterations in chromatin structure including histone marks, DNase I sensitivity, repressed chromatin regions as well as transcript factor binding sites. Integration of this information can partially explain why particular variations can alter proclivity to autoimmune phenotypes. Despite our incomplete knowledge base with only partial definition of hereditary factors and possible functional connections, this progress has and will continue to facilitate a better understanding of critical pathways and critical changes in immunoregulation. Advances in defining and understanding functional variants potentially can lead to both novel therapeutics and personalized medicine in which therapeutic approaches are chosen based on particular molecular phenotypes and genomic alterations.

Gene set analysis: A step-by-step guide

To maximize the potential of genome-wide association studies, many researchers are performing secondary analyses to identify sets of genes jointly associated with the trait of interest. Although methods for gene-set analyses (GSA), also called pathway analyses, have been around for more than a decade, the field is still evolving. There are numerous algorithms available for testing the cumulative effect of multiple SNPs, yet no real consensus in the field about the best way to perform a GSA. This paper provides an overview of the factors that can affect the results of a GSA, the lessons learned from past studies, and suggestions for how to make analysis choices that are most appropriate for different types of data. © 2015 Wiley Periodicals, Inc.

Phenome-Wide Association Studies: Embracing Complexity for Discovery

The inherent complexity of biological systems can be leveraged for a greater understanding of the impact of genetic architecture on outcomes, traits, and pharmacological response. The genome-wide association study (GWAS) approach has well-developed methods and relatively straight-forward methodologies; however, the bigger picture of the impact of genetic architecture on phenotypic outcome still remains to be elucidated even with an ever-growing number of GWAS performed. Greater consideration of the complexity of biological processes, using more data from the phenome, exposome, and diverse -omic resources, including considering the interplay of pleiotropy and genetic interactions, may provide additional leverage for making the most of the incredible wealth of information available for study. Here, we describe how incorporating greater complexity into analyses through the use of additional phenotypic data and widespread deployment of phenome-wide association studies may provide new insights into genetic factors influencing diseases, traits, and pharmacological response.

Heterogeneous Network Edge Prediction: A Data Integration Approach to Prioritize Disease-Associated Genes

The first decade of Genome Wide Association Studies (GWAS) has uncovered a wealth of disease-associated variants. Two important derivations will be the translation of this information into a multiscale understanding of pathogenic variants and leveraging existing data to increase the power of existing and future studies through prioritization. We explore edge prediction on heterogeneous networks—graphs with multiple node and edge types—for accomplishing both tasks. First we constructed a network with 18 node types—genes, diseases, tissues, pathophysiologies, and 14 MSigDB (molecular signatures database) collections—and 19 edge types from high-throughput publicly-available resources. From this network composed of 40,343 nodes and 1,608,168 edges, we extracted features that describe the topology between specific genes and diseases. Next, we trained a model from GWAS associations and predicted the probability of association between each protein-coding gene and each of 29 well-studied complex diseases. The model, which achieved 132-fold enrichment in precision at 10% recall, outperformed any individual domain, highlighting the benefit of integrative approaches. We identified pleiotropy, transcriptional signatures of perturbations, pathways, and protein interactions as influential mechanisms explaining pathogenesis. Our method successfully predicted the results (with AUROC = 0.79) from a withheld multiple sclerosis (MS) GWAS despite starting with only 13 previously associated genes. Finally, we combined our network predictions with statistical evidence of association to propose four novel MS genes, three of which (JAK2, REL, RUNX3) validated on the masked GWAS. Furthermore, our predictions provide biological support highlighting REL as the causal gene within its gene-rich locus. Users can browse all predictions online (http://het.io). Heterogeneous network edge prediction effectively prioritized genetic associations and provides a powerful new approach for data integration across multiple domains.

For complex human diseases, identifying the genes harboring susceptibility variants has taken on medical importance. Disease-associated genes provide clues for elucidating disease etiology, predicting disease risk, and highlighting therapeutic targets. Here, we develop a method to predict whether a given gene and disease are associated. To capture the multitude of biological entities underlying pathogenesis, we constructed a heterogeneous network, containing multiple node and edge types. We built on a technique developed for social network analysis, which embraces disparate sources of data to make predictions from heterogeneous networks. Using the compendium of associations from genome-wide studies, we learned the influential mechanisms underlying pathogenesis. Our findings provide a novel perspective about the existence of pervasive pleiotropy across complex diseases. Furthermore, we suggest transcriptional signatures of perturbations are an underutilized resource amongst prioritization approaches. For multiple sclerosis, we demonstrated our ability to prioritize future studies and discover novel susceptibility genes. Researchers can use these predictions to increase the statistical power of their studies, to suggest the causal genes from a set of candidates, or to generate evidence-based experimental hypothesis.

Bioinformatics tools for discovery and functional analysis of single nucleotide polymorphisms

With the high speed DNA sequencing of genome, databases of genome data continue to grow, and the understanding of genetic variation between individuals grows as well. Single nucleotide polymorphisms (SNPs), a main type of genetic variation, are increasingly important resource for understanding the structure and function of the human genome and become a valuable resource for investigating the genetic basis of disease. During the past years, in addition to experimental approaches to characterize specific variants, intense bioinformatics techniques were applied to understand effects of these genetic changes. In the genetics studies, one intends to understand the molecular basis of disease, and computational methods are becoming increasingly important for SNPs selection, prediction and understanding the downstream effects of genetic variation. The review provides systematic information on the available resources and methods for SNPs discovery and analysis. We also report some new results on DNA sequence-based prediction of SNPs in human cytochrome P450, which serves as an example of computational methods to predict and discovery SNPs. Additionally, annotation and prediction of functional SNPs, as well as a comprehensive list of existing tools and online recourses, are reviewed and described.

Integrated systems analysis reveals a molecular network underlying autism spectrum disorders

Autism is a complex disease whose etiology remains elusive. We integrated previously and newly generated data and developed a systems framework involving the interactome, gene expression and genome sequencing to identify a protein interaction module with members strongly enriched for autism candidate genes. Sequencing of 25 patients confirmed the involvement of this module in autism, which was subsequently validated using an independent cohort of over 500 patients. Expression of this module was dichotomized with a ubiquitously expressed subcomponent and another subcomponent preferentially expressed in the corpus callosum, which was significantly affected by our identified mutations in the network center. RNA-sequencing of the corpus callosum from patients with autism exhibited extensive gene mis-expression in this module, and our immunochemical analysis showed that the human corpus callosum is predominantly populated by oligodendrocyte cells. Analysis of functional genomic data further revealed a significant involvement of this module in the development of oligodendrocyte cells in mouse brain. Our analysis delineates a natural network involved in autism, helps uncover novel candidate genes for this disease and improves our understanding of its molecular pathology.

Genomic architecture of pharmacological efficacy and adverse events

The pharmacokinetic and pharmacodynamic disciplines address pharmacological traits, including efficacy and adverse events. Pharmacogenomics studies have identified pervasive genetic effects on treatment outcomes, resulting in the development of genetic biomarkers for optimization of drug therapy. Pharmacogenomics-based tests are already being applied in clinical decision making. However, despite substantial progress in identifying the genetic etiology of pharmacological response, current biomarker panels still largely rely on single gene tests with a large portion of the genetic effects remaining to be discovered. Future research must account for the combined effects of multiple genetic variants, incorporate pathway-based approaches, explore gene-gene interactions and nonprotein coding functional genetic variants, extend studies across ancestral populations, and prioritize laboratory characterization of molecular mechanisms. Because genetic factors can play a key role in drug response, accurate biomarker tests capturing the main genetic factors determining treatment outcomes have substantial potential for improving individual clinical care.

Protein interaction networks reveal novel autism risk genes within GWAS statistical noise

Genome-wide association studies (GWAS) for Autism Spectrum Disorder (ASD) thus far met limited success in the identification of common risk variants, consistent with the notion that variants with small individual effects cannot be detected individually in single SNP analysis. To further capture disease risk gene information from ASD association studies, we applied a network-based strategy to the Autism Genome Project (AGP) and the Autism Genetics Resource Exchange GWAS datasets, combining family-based association data with Human Protein-Protein interaction (PPI) data. Our analysis showed that autism-associated proteins at higher than conventional levels of significance (P<0.1) directly interact more than random expectation and are involved in a limited number of interconnected biological processes, indicating that they are functionally related. The functionally coherent networks generated by this approach contain ASD-relevant disease biology, as demonstrated by an improved positive predictive value and sensitivity in retrieving known ASD candidate genes relative to the top associated genes from either GWAS, as well as a higher gene overlap between the two ASD datasets. Analysis of the intersection between the networks obtained from the two ASD GWAS and six unrelated disease datasets identified fourteen genes exclusively present in the ASD networks. These are mostly novel genes involved in abnormal nervous system phenotypes in animal models, and in fundamental biological processes previously implicated in ASD, such as axon guidance, cell adhesion or cytoskeleton organization. Overall, our results highlighted novel susceptibility genes previously hidden within GWAS statistical "noise" that warrant further analysis for causal variants.

Hypothesis-independent pathway analysis implicates GABA and acetyl-CoA metabolism in primary open-angle glaucoma and normal-pressure glaucoma

Primary open-angle glaucoma (POAG) is a leading cause of blindness worldwide. Using genome-wide association single-nucleotide polymorphism data from the Glaucoma Genes and Environment study and National Eye Institute Glaucoma Human Genetics Collaboration comprising 3,108 cases and 3,430 controls, we assessed biologic pathways as annotated in the KEGG database for association with risk of POAG. After correction for genic overlap among pathways, we found 4 pathways, butanoate metabolism (hsa00650), hematopoietic cell lineage (hsa04640), lysine degradation (hsa00310) and basal transcription factors (hsa03022) related to POAG with permuted p < 0.001. In addition, the human leukocyte antigen (HLA) gene family was significantly associated with POAG (p < 0.001). In the POAG subset with normal-pressure glaucoma (NPG), the butanoate metabolism pathway was also significantly associated (p < 0.001) as well as the MAPK and Hedgehog signaling pathways (hsa04010 and hsa04340), glycosaminoglycan biosynthesis-heparan sulfate pathway (hsa00534) and the phenylalanine, tyrosine and tryptophan biosynthesis pathway (hsa0400). The butanoate metabolism pathway overall, and specifically the aspects of the pathway that contribute to GABA and acetyl-CoA metabolism, was the only pathway significantly associated with both POAG and NPG. Collectively these results implicate GABA and acetyl-CoA metabolism in glaucoma pathogenesis, and suggest new potential therapeutic targets.

PINBPA: cytoscape app for network analysis of GWAS data

Protein interaction network-based pathway analysis (PINBPA) for genome-wide association studies (GWAS) has been developed as a Cytoscape app, to enable analysis of GWAS data in a network fashion. Users can easily import GWAS summary-level data, draw Manhattan plots, define blocks, prioritize genes with random walk with restart, detect enriched subnetworks and test the significance of subnetworks via a user-friendly interface.

AVAILABILITY AND IMPLEMENTATION

PINBPA app is freely available in Cytoscape app store.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide association studies: getting to pathogenesis, the role of inflammation/complement in age-related macular degeneration

Age-related macular degeneration (AMD) is a chronic, degenerative, and significant cause of visual impairment and blindness in the elderly. Genetic and epidemiological studies have confirmed that AMD has a strong genetic component, which has encouraged the application of increasingly sophisticated genetic techniques to uncover the important underlying genetic variants. Although various genes and pathways have been implicated in the risk for AMD, complement activation has been emphasized repeatedly throughout the literature as having a major role both physiologically and genetically in susceptibility to and pathogenesis of this disease. This article explores the research efforts that brought about the discovery and characterization of the role of inflammatory and immune processes (specifically complement) in AMD. The focus herein is on the genetic evidence for the role of complement in AMD as supported specifically by genome-wide association (GWA) studies, which interrogate hundreds of thousands of variants across the genome in a hypothesis-free approach, and other genetic interrogation methods.

Functional and genomic context in pathway analysis of GWAS data

Gene set analysis (GSA) is a promising tool for uncovering the polygenic effects associated with complex diseases. However, the available techniques reflect a wide variety of hypotheses about how genetic effects interact to contribute to disease susceptibility. The lack of consensus about the best way to perform GSA has led to confusion in the field and has made it difficult to compare results across methods. A clear understanding of the various choices made during GSA - such as how gene sets are defined, how single-nucleotide polymorphisms (SNPs) are assigned to genes, and how individual SNP-level effects are aggregated to produce gene- or pathway-level effects - will improve the interpretability and comparability of results across methods and studies. In this review we provide an overview of the various data sources used to construct gene sets and the statistical methods used to test for gene set association, as well as provide guidelines for ensuring the comparability of results.

The impact of the human genome project on complex disease

In the decade that has passed since the initial release of the Human Genome, numerous advancements in science and technology within and beyond genetics and genomics have been encouraged and enhanced by the availability of this vast and remarkable data resource. Progress in understanding three common, complex diseases: age-related macular degeneration (AMD), Alzheimer's disease (AD), and multiple sclerosis (MS), are three exemplars of the incredible impact on the elucidation of the genetic architecture of disease. The approaches used in these diseases have been successfully applied to numerous other complex diseases. For example, the heritability of AMD was confirmed upon the release of the first genome-wide association study (GWAS) along with confirmatory reports that supported the findings of that state-of-the art method, thus setting the foundation for future GWAS in other heritable diseases. Following this seminal discovery and applying it to other diseases including AD and MS, the genetic knowledge of AD expanded far beyond the well-known APOE locus and now includes more than 20 loci. MS genetics saw a similar increase beyond the HLA loci and now has more than 100 known risk loci. Ongoing and future efforts will seek to define the remaining heritability of these diseases; the next decade could very well hold the key to attaining this goal.

Vascular tone pathway polymorphisms in relation to primary open-angle glaucoma

AIMS

Vascular perfusion may be impaired in primary open-angle glaucoma (POAG); thus, we evaluated a panel of markers in vascular tone-regulating genes in relation to POAG.

METHODS

We used Illumina 660W-Quad array genotype data and pooled P-values from 3108 POAG cases and 3430 controls from the combined National Eye Institute Glaucoma Human Genetics Collaboration consortium and Glaucoma Genes and Environment studies. Using information from previous literature and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, we compiled single-nucleotide polymorphisms (SNPs) in 186 vascular tone-regulating genes. We used the 'Pathway Analysis by Randomization Incorporating Structure' analysis software, which performed 1000 permutations to compare the overall pathway and selected genes with comparable randomly generated pathways and genes in their association with POAG.

RESULTS

The vascular tone pathway was not associated with POAG overall or POAG subtypes, defined by the type of visual field loss (early paracentral loss (n=224 cases) or only peripheral loss (n=993 cases)) (permuted P≥0.20). In gene-based analyses, eight were associated with POAG overall at permuted P<0.001: PRKAA1, CAV1, ITPR3, EDNRB, GNB2, DNM2, HFE, and MYL9. Notably, six of these eight (the first six listed) code for factors involved in the endothelial nitric oxide synthase activity, and three of these six (CAV1, ITPR3, and EDNRB) were also associated with early paracentral loss at P<0.001, whereas none of the six genes reached P<0.001 for peripheral loss only.

DISCUSSION

Although the assembled vascular tone SNP set was not associated with POAG, genes that code for local factors involved in setting vascular tone were associated with POAG.

Nutrigenetics, nutrigenomics and inflammatory bowel diseases

Inflammatory bowel disease includes ulcerative colitis and Crohn's disease, which are both inflammatory disorders of the gastrointestinal tract. Both types of inflammatory bowel disease have a complex etiology, resulting from a genetically determined susceptibility interacting with environmental factors, including the diet and gut microbiota. Genome Wide Association Studies have implicated more than 160 single-nucleotide polymorphisms in disease susceptibility. Consideration of the different pathways suggested to be involved implies that specific dietary interventions are likely to be appropriate, dependent upon the nature of the genes involved. Epigenetics and the gut microbiota are also responsive to dietary interventions. Nutrigenetics may lead to personalized nutrition for disease prevention and treatment, while nutrigenomics may help to understand the nature of the disease and individual response to nutrients.

Incorporating prior knowledge to increase the power of genome-wide association studies

Typical methods of analyzing genome-wide single nucleotide variant (SNV) data in cases and controls involve testing each variant's genotypes separately for phenotype association, and then using a substantial multiple-testing penalty to minimize the rate of false positives. This approach, however, can result in low power for modestly associated SNVs. Furthermore, simply looking at the most associated SNVs may not directly yield biological insights about disease etiology. SNVset methods attempt to address both limitations of the traditional approach by testing biologically meaningful sets of SNVs (e.g., genes or pathways). The number of tests run in a SNVset analysis is typically much lower (hundreds or thousands instead of millions) than in a traditional analysis, so the false-positive rate is lower. Additionally, by testing SNVsets that are biologically meaningful finding a significant set may more quickly yield insights into disease etiology.In this chapter we summarize the short history of SNVset testing and provide an overview of the many recently proposed methods. Furthermore, we provide detailed step-by-step instructions on how to perform a SNVset analysis, including a substantial number of practical tips and questions that researchers should consider before undertaking a SNVset analysis. Lastly, we describe a companion R package (snvset) that implements recently proposed SNVset methods. While SNVset testing is a new approach, with many new methods still being developed and many open questions, the promise of the approach is worth serious consideration when considering analytic methods for GWAS.

Pathway-PDT: a flexible pathway analysis tool for nuclear families

Background

Pathway analysis based on Genome-Wide Association Studies (GWAS) data has become popular as a secondary analysis strategy. Although many pathway analysis tools have been developed for case–control studies, there is no tool that can use all information from raw genotypes in general nuclear families. We developed Pathway-PDT, which uses the framework of Pedigree Disequilibrium Test (PDT) for general family data, to perform pathway analysis based on raw genotypes in family-based GWAS.

Results

Simulation results showed that Pathway-PDT is more powerful than the p-value based method, ALIGATOR. Pathway-PDT also can be more powerful than the PLINK set-based test when analyzing general nuclear families with multiple siblings or missing parents. Additionally, Pathway-PDT has a flexible and convenient user interface, which allows users to modify their analysis parameters as well as to apply various types of gene and pathway definitions.

Conclusions

The Pathway-PDT method is implemented in C++ with POSIX threads and is computationally feasible for pathway analysis with large scale family GWAS datasets. The Windows binary along with Makefile and source codes for the Linux are available at https://sourceforge.net/projects/pathway-pdt/.

The genetic landscapes of autism spectrum disorders

The autism spectrum disorders (ASD) are characterized by impairments in social interaction and stereotyped behaviors. For the majority of individuals with ASD, the causes of the disorder remain unknown; however, in up to 25% of cases, a genetic cause can be identified. Chromosomal rearrangements as well as rare and de novo copy-number variants are present in ∼10-20% of individuals with ASD, compared with 1-2% in the general population and/or unaffected siblings. Rare and de novo coding-sequence mutations affecting neuronal genes have also been identified in ∼5-10% of individuals with ASD. Common variants such as single-nucleotide polymorphisms seem to contribute to ASD susceptibility, but, taken individually, their effects appear to be small. Despite a heterogeneous genetic landscape, the genes implicated thus far-which are involved in chromatin remodeling, metabolism, mRNA translation, and synaptic function-seem to converge in common pathways affecting neuronal and synaptic homeostasis. Animal models developed to study these genes should lead to a better understanding of the diversity of the genetic landscapes of ASD.

Estrogen pathway polymorphisms in relation to primary open angle glaucoma: an analysis accounting for gender from the United States

PURPOSE

Circulating estrogen levels are relevant in glaucoma phenotypic traits. We assessed the association between an estrogen metabolism single nucleotide polymorphism (SNP) panel in relation to primary open angle glaucoma (POAG), accounting for gender.

METHODS

We included 3,108 POAG cases and 3,430 controls of both genders from the Glaucoma Genes and Environment (GLAUGEN) study and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium genotyped on the Illumina 660W-Quad platform. We assessed the relation between the SNP panels representative of estrogen metabolism and POAG using pathway- and gene-based approaches with the Pathway Analysis by Randomization Incorporating Structure (PARIS) software. PARIS executes a permutation algorithm to assess statistical significance relative to the pathways and genes of comparable genetic architecture. These analyses were performed using the meta-analyzed results from the GLAUGEN and NEIGHBOR data sets. We evaluated POAG overall as well as two subtypes of POAG defined as intraocular pressure (IOP) ≥22 mmHg (high-pressure glaucoma [HPG]) or IOP <22 mmHg (normal pressure glaucoma [NPG]) at diagnosis. We conducted these analyses for each gender separately and then jointly in men and women.

RESULTS

Among women, the estrogen SNP pathway was associated with POAG overall (permuted p=0.006) and HPG (permuted p<0.001) but not NPG (permuted p=0.09). Interestingly, there was no relation between the estrogen SNP pathway and POAG when men were considered alone (permuted p>0.99). Among women, gene-based analyses revealed that the catechol-O-methyltransferase gene showed strong associations with HTG (permuted gene p≤0.001) and NPG (permuted gene p=0.01).

CONCLUSIONS

The estrogen SNP pathway was associated with POAG among women.

Neurobiology meets genomic science: the promise of human-induced pluripotent stem cells

The recent introduction of the induced pluripotent stem cell technology has made possible the derivation of neuronal cells from somatic cells obtained from human individuals. This in turn has opened new areas of investigation that can potentially bridge the gap between neuroscience and psychopathology. For the first time we can study the cell biology and genetics of neurons derived from any individual. Furthermore, by recapitulating in vitro the developmental steps whereby stem cells give rise to neuronal cells, we can now hope to understand factors that control typical and atypical development. We can begin to explore how human genes and their variants are transcribed into messenger RNAs within developing neurons and how these gene transcripts control the biology of developing cells. Thus, human-induced pluripotent stem cells have the potential to uncover not only what aspects of development are uniquely human but also variations in the series of events necessary for normal human brain development that predispose to psychopathology.

Individual common variants exert weak effects on the risk for autism spectrum disorders

While it is apparent that rare variation can play an important role in the genetic architecture of autism spectrum disorders (ASDs), the contribution of common variation to the risk of developing ASD is less clear. To produce a more comprehensive picture, we report Stage 2 of the Autism Genome Project genome-wide association study, adding 1301 ASD families and bringing the total to 2705 families analysed (Stages 1 and 2). In addition to evaluating the association of individual single nucleotide polymorphisms (SNPs), we also sought evidence that common variants, en masse, might affect the risk. Despite genotyping over a million SNPs covering the genome, no single SNP shows significant association with ASD or selected phenotypes at a genome-wide level. The SNP that achieves the smallest P-value from secondary analyses is rs1718101. It falls in CNTNAP2, a gene previously implicated in susceptibility for ASD. This SNP also shows modest association with age of word/phrase acquisition in ASD subjects, of interest because features of language development are also associated with other variation in CNTNAP2. In contrast, allele scores derived from the transmission of common alleles to Stage 1 cases significantly predict case status in the independent Stage 2 sample. Despite being significant, the variance explained by these allele scores was small (Vm< 1%). Based on results from individual SNPs and their en masse effect on risk, as inferred from the allele score results, it is reasonable to conclude that common variants affect the risk for ASD but their individual effects are modest.

Loci nominally associated with autism from genome-wide analysis show enrichment of brain expression quantitative trait loci but not lymphoblastoid cell line expression quantitative trait loci

Background

Autism spectrum disorder is a severe early onset neurodevelopmental disorder with high heritability but significant heterogeneity. Traditional genome-wide approaches to test for an association of common variants with autism susceptibility risk have met with limited success. However, novel methods to identify moderate risk alleles in attainable sample sizes are now gaining momentum.

Methods

In this study, we utilized publically available genome-wide association study data from the Autism Genome Project and annotated the results (P <0.001) for expression quantitative trait loci present in the parietal lobe (GSE35977), cerebellum (GSE35974) and lymphoblastoid cell lines (GSE7761). We then performed a test of enrichment by comparing these results to simulated data conditioned on minor allele frequency to generate an empirical P-value indicating statistically significant enrichment of expression quantitative trait loci in top results from the autism genome-wide association study.

Results

Our findings show a global enrichment of brain expression quantitative trait loci, but not lymphoblastoid cell line expression quantitative trait loci, among top single nucleotide polymorphisms from an autism genome-wide association study. Additionally, the data implicates individual genes SLC25A12, PANX1 and PANX2 as well as pathways previously implicated in autism.

Conclusions

These findings provide supportive rationale for the use of annotation-based approaches to genome-wide association studies.

Common variants at 9p21 and 8q22 are associated with increased susceptibility to optic nerve degeneration in glaucoma

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63-0.75], p = 1.86×10⁻¹⁸), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21-1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50-0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53-0.72], p = 8.88×10⁻¹⁰). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41-0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54-1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

Infant siblings and the investigation of autism risk factors

Infant sibling studies have been at the vanguard of autism spectrum disorders (ASD) research over the past decade, providing important new knowledge about the earliest emerging signs of ASD and expanding our understanding of the developmental course of this complex disorder. Studies focused on siblings of children with ASD also have unrealized potential for contributing to ASD etiologic research. Moving targeted time of enrollment back from infancy toward conception creates tremendous opportunities for optimally studying risk factors and risk biomarkers during the pre-, peri- and neonatal periods. By doing so, a traditional sibling study, which already incorporates close developmental follow-up of at-risk infants through the third year of life, is essentially reconfigured as an enriched-risk pregnancy cohort study. This review considers the enriched-risk pregnancy cohort approach of studying infant siblings in the context of current thinking on ASD etiologic mechanisms. It then discusses the key features of this approach and provides a description of the design and implementation strategy of one major ASD enriched-risk pregnancy cohort study: the Early Autism Risk Longitudinal Investigation (EARLI).

Pathway analysis of genomic data: concepts, methods, and prospects for future development

Genome-wide data sets are increasingly being used to identify biological pathways and networks underlying complex diseases. In particular, analyzing genomic data through sets defined by functional pathways offers the potential of greater power for discovery and natural connections to biological mechanisms. With the burgeoning availability of next-generation sequencing, this is an opportune moment to revisit strategies for pathway-based analysis of genomic data. Here, we synthesize relevant concepts and extant methodologies to guide investigators in study design and execution. We also highlight ongoing challenges and proposed solutions. As relevant analytical strategies mature, pathways and networks will be ideally placed to integrate data from diverse -omics sources to harness the extensive, rich information related to disease and treatment mechanisms.

Strategies for pathway analysis from GWAS data

Genome-wide association studies (GWAS) 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 error concerns. 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, as they are difficult to separate from a sea of false-positive associations. Organizing these individual SNPs into biologically meaningful groups to look at 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.

Links between genetics and pathophysiology in the autism spectrum disorders

Autism spectrum disorders (ASD) are important neuropsychiatric disorders, currently estimated to affect approximately 1% of children, with considerable emotional and financial costs. Significant collaborative effort has been made over the last 15 years in an attempt to unravel the genetic mechanisms underlying these conditions. This has led to important discoveries, both of the roles of specific genes, as well as larger scale chromosomal copy number changes. Here, we summarize some of the latest genetic findings in the field of ASD and attempt to link them with the results of pathophysiological studies to provide an overall picture of at least one of the mechanisms by which ASD may develop.