A temporal classifier predicts histopathology state and parses acute-chronic phasing in inflammatory bowel disease patients.

Previous studies have conducted time course characterization of murine colitis models through transcriptional profiling of differential expression. We characterize the transcriptional landscape of acute and chronic models of dextran sodium sulfate (DSS) and adoptive transfer (AT) colitis to derive temporal gene expression and splicing signatures in blood and colonic tissue in order to capture dynamics of colitis remission and relapse. We identify sub networks of patient-derived causal networks that are enriched in these temporal signatures to distinguish acute and chronic disease components within the broader molecular landscape of IBD. The interaction between the DSS phenotype and chronological time-point naturally defines parsimonious temporal gene expression and splicing signatures associated with acute and chronic phases disease (as opposed to ordinary time-specific differential expression/splicing). We show these expression and splicing signatures are largely orthogonal, i.e. affect different genetic bodies, and that using machine learning, signatures are predictive of histopathological measures from both blood and intestinal data in murine colitis models as well as an independent cohort of IBD patients. Through access to longitudinal multi-scale profiling from disease tissue in IBD patient cohorts, we can apply this machine learning pipeline to generation of direct patient temporal multimodal regulatory signatures for prediction of histopathological outcomes.

Evolving use of real-world evidence in the regulatory process: a focus on immuno-oncology treatment and outcomes.

In recent years, regulatory bodies have increasingly recognized the utility of real-world evidence (RWE) for supplementing and supporting clinical trial data in new drug applications. Nevertheless, the integration of RWE into established regulatory processes is complex and the generation of 'regulatory-grade' real-world data faces operational, methodological, data-related and policy-related challenges. In parallel with this evolving role for RWE, immuno-oncology therapies have emerged as leading cancer treatments and are expected to continue to play a central role in the future. In this article, we review the current literature on the use of RWE for regulatory submissions, with a focus on novel anticancer immunotherapies, and discuss the utility and current limitations of RWE in the context of drug development and regulatory approvals.

STK11 and KEAP1 mutations as prognostic biomarkers in an observational real-world lung adenocarcinoma cohort.

INTRODUCTION: Somatic mutations in STK11 and KEAP1, frequently comutated in non-squamous non-small cell lung cancer (NSQ NSCLC), have been associated with poor response to immune checkpoint blockade (ICB). However, previous reports lack non-ICB controls needed to properly ascertain the predictive nature of those biomarkers. The objective of this study was to evaluate the predictive versus prognostic effect of STK11 or KEAP1 mutations in NSQ NSCLC. METHODS: Patients diagnosed with stage IIIB, IIIC, IVA or IVB NSQ NSCLC from a real-world data cohort from the Flatiron Health Network linked with genetic testing from Foundation Medicine were retrospectively assessed. Real-world, progression-free survival (rwPFS) and overall survival (OS) were calculated from time of initiation of first-line treatment. RESULTS: We analysed clinical and mutational data for 2276 patients including patients treated with anti-programmed death-1 (PD-1)/anti-programmed death ligand 1 (PD-L1) inhibitors at first line (n=574). Mutations in STK11 or KEAP1 were associated with poor outcomes across multiple therapeutic classes and were not specifically associated with poor outcomes in ICB cohorts. There was no observable interaction between STK11 mutations and anti-PD-1/anti-PD-L1 treatment on rwPFS (HR, 1.05; 95% CI 0.76 to 1.44; p=0.785) or OS (HR, 1.13; 95% CI 0.76 to 1.67; p=0.540). Similarly, there was no observable interaction between KEAP1 mutations and treatment on rwPFS (HR, 0.93; 95% CI 0.67 to 1.28; p=0.653) or OS (HR, 0.98; 95% CI 0.66 to 1.45; p=0.913). CONCLUSION: Our results show that STK11-KEAP1 mutations are prognostic, not predictive, biomarkers for anti-PD-1/anti-PD-L1 therapy.

High-Throughput Identification of the Plasma Proteomic Signature of Inflammatory Bowel Disease.

BACKGROUND: The molecular aetiology of inflammatory bowel disease [IBD] and its two subtypes, ulcerative colitis [UC] and Crohn's disease [CD], have been carefully investigated at genome and transcriptome levels. Recent advances in high-throughput proteome quantification has enabled comprehensive large-scale plasma proteomics studies of IBD. METHODS: The study used two cohorts: [1] The CERTIFI-cohort: 42 samples from the CERTIFI trial of anti-TNFα-refractory CD patients; [2] the PROgECT-UNITI-HCs cohort: 46 UC samples of the PROgECT study, 84 CD samples of the UNITI I and UNITI II studies, and 72 healthy controls recruited in Mount Sinai Hospital, New York, USA. The plasma proteome for these two cohorts was quantified using high-throughput platforms. RESULTS: For the PROgECT-UNITI-HCs cohort, we measured a total of 1310 proteins. Of these, 493 proteins showed different plasma levels in IBD patients to the plasma levels in controls at 10% false discovery rate [FDR], among which 11 proteins had a fold change greater than 2. The proteins upregulated in IBD were associated with immunity functionality, whereas the proteins downregulated in IBD were associated with nutrition and metabolism. The proteomic profiles were very similar between UC and CD. In the CERTIFI cohort, 1014 proteins were measured, and it was found that the plasma protein level had little correlation with the blood or intestine transcriptomes. CONCLUSIONS: We report the largest proteomics study to date on IBD and controls. A large proportion of plasma proteins are altered in IBD, which provides insights into the disease aetiology and indicates a potential for biomarker discovery.

A functional genomics predictive network model identifies regulators of inflammatory bowel disease.

A major challenge in inflammatory bowel disease (IBD) is the integration of diverse IBD data sets to construct predictive models of IBD. We present a predictive model of the immune component of IBD that informs causal relationships among loci previously linked to IBD through genome-wide association studies (GWAS) using functional and regulatory annotations that relate to the cells, tissues, and pathophysiology of IBD. Our model consists of individual networks constructed using molecular data generated from intestinal samples isolated from three populations of patients with IBD at different stages of disease. We performed key driver analysis to identify genes predicted to modulate network regulatory states associated with IBD, prioritizing and prospectively validating 12 of the top key drivers experimentally. This validated key driver set not only introduces new regulators of processes central to IBD but also provides the integrated circuits of genetic, molecular, and clinical traits that can be directly queried to interrogate and refine the regulatory framework defining IBD.

Integrating personalized gene expression profiles into predictive disease-associated gene pools.

Gene expression data are routinely used to identify genes that on average exhibit different expression levels between a case and a control group. Yet, very few of such differentially expressed genes are detectably perturbed in individual patients. Here, we develop a framework to construct personalized perturbation profiles for individual subjects, identifying the set of genes that are significantly perturbed in each individual. This allows us to characterize the heterogeneity of the molecular manifestations of complex diseases by quantifying the expression-level similarities and differences among patients with the same phenotype. We show that despite the high heterogeneity of the individual perturbation profiles, patients with asthma, Parkinson and Huntington's disease share a broadpool of sporadically disease-associated genes, and that individuals with statistically significant overlap with this pool have a 80-100% chance of being diagnosed with the disease. The developed framework opens up the possibility to apply gene expression data in the context of precision medicine, with important implications for biomarker identification, drug development, diagnosis and treatment.

High-Throughput Characterization of Blood Serum Proteomics of IBD Patients with Respect to Aging and Genetic Factors.

To date, no large scale, systematic description of the blood serum proteome has been performed in inflammatory bowel disease (IBD) patients. By using microarray technology, a more complete description of the blood proteome of IBD patients is feasible. It may help to achieve a better understanding of the disease. We analyzed blood serum profiles of 1128 proteins in IBD patients of European descent (84 Crohn's Disease (CD) subjects and 88 Ulcerative Colitis (UC) subjects) as well as 15 healthy control subjects, and linked protein variability to patient age (all cohorts) and genetic components (genotype data generated from CD patients). We discovered new, previously unreported aging-associated proteomic traits (such as serum Albumin level), confirmed previously reported results from different tissues (i.e., upregulation of APOE with aging), and found loss of regulation of MMP7 in CD patients. In carrying out a genome wide genotype-protein association study (proteomic Quantitative Trait Loci, pQTL) within the CD patients, we identified 41 distinct proteomic traits influenced by cis pQTLs (underlying SNPs are referred to as pSNPs). Significant overlaps between pQTLs and cis eQTLs corresponding to the same gene were observed and in some cases the QTL were related to inflammatory disease susceptibility. Importantly, we discovered that serum protein levels of MST1 (Macrophage Stimulating 1) were regulated by SNP rs3197999 (p = 5.96E-10, FDR<5%), an accepted GWAS locus for IBD. Filling the knowledge gap of molecular mechanisms between GWAS hits and disease susceptibility requires systematically dissecting the impact of the locus at the cell, mRNA expression, and protein levels. The technology and analysis tools that are now available for large-scale molecular studies can elucidate how alterations in the proteome driven by genetic polymorphisms cause or provide protection against disease. Herein, we demonstrated this directly by integrating proteomic and pQTLs with existing GWAS, mRNA expression, and eQTL datasets to provide insights into the biological processes underlying IBD and pinpoint causal genetic variants along with their downstream molecular consequences.

Blood and Intestine eQTLs from an Anti-TNF-Resistant Crohn's Disease Cohort Inform IBD Genetic Association Loci.

OBJECTIVES: Genome-wide association studies (GWAS) have identified loci reproducibly associated with inflammatory bowel disease (IBD) and other immune-mediated diseases; however, the molecular mechanisms underlying most of genetic susceptibility remain undefined. Expressional quantitative trait loci (eQTL) of disease-relevant tissue can be employed in order to elucidate the genes and pathways affected by disease-specific genetic variance. METHODS: In this study, we derived eQTLs for human whole blood and intestine tissues of anti-tumor necrosis factor-resistant Crohn's disease (CD) patients. We interpreted these eQTLs in the context of published IBD GWAS hits to inform on the disease process. RESULTS: At 10% false discovery rate, we discovered that 5,174 genes in blood and 2,063 genes in the intestine were controlled by a nearby single-nucleotide polymorphism (SNP) (i.e., cis-eQTL), among which 1,360 were shared between the two tissues. A large fraction of the identified eQTLs were supported by the regulomeDB database, showing that the eQTLs reside in regulatory elements (odds ratio; OR=3.44 and 3.24 for blood and intestine eQTLs, respectively) as opposed to protein-coding regions. Published IBD GWAS hits as a whole were enriched for blood and intestine eQTLs (OR=2.88 and 2.05; and P value=2.51E-9 and 0.013, respectively), thereby linking genetic susceptibility to control of gene expression in these tissues. Through a systematic search, we used eQTL data to inform 109 out of 372 IBD GWAS SNPs documented in National Human Genome Research Institute catalog, and we categorized the genes influenced by eQTLs according to their functions. Many of these genes have experimentally validated roles in specific cell types contributing to intestinal inflammation. CONCLUSIONS: The blood and intestine eQTLs described in this study represent a powerful tool to link GWAS loci to a regulatory function and thus elucidate the mechanisms underlying the genetic loci associated with IBD and related conditions. Overall, our eQTL discovery approach empirically identifies the disease-associated variants including their impact on the direction and extent of expression changes in the context of disease-relevant cellular pathways in order to infer the functional outcome of this aspect of genetic susceptibility.

Integrative genomic deconvolution of rheumatoid arthritis GWAS loci into gene and cell type associations.

BACKGROUND: Although genome-wide association studies (GWAS) have identified over 100 genetic loci associated with rheumatoid arthritis (RA), our ability to translate these results into disease understanding and novel therapeutics is limited. Most RA GWAS loci reside outside of protein-coding regions and likely affect distal transcriptional enhancers. Furthermore, GWAS do not identify the cell types where the associated causal gene functions. Thus, mapping the transcriptional regulatory roles of GWAS hits and the relevant cell types will lead to better understanding of RA pathogenesis. RESULTS: We combine the whole-genome sequences and blood transcription profiles of 377 RA patients and identify over 6000 unique genes with expression quantitative trait loci (eQTLs). We demonstrate the quality of the identified eQTLs through comparison to non-RA individuals. We integrate the eQTLs with immune cell epigenome maps, RA GWAS risk loci, and adjustment for linkage disequilibrium to propose target genes of immune cell enhancers that overlap RA risk loci. We examine 20 immune cell epigenomes and perform a focused analysis on primary monocytes, B cells, and T cells. CONCLUSIONS: We highlight cell-specific gene associations with relevance to RA pathogenesis including the identification of FCGR2B in B cells as possessing both intragenic and enhancer regulatory GWAS hits. We show that our RA patient cohort derived eQTL network is more informative for studying RA than that from a healthy cohort. While not experimentally validated here, the reported eQTLs and cell type-specific RA risk associations can prioritize future experiments with the goal of elucidating the regulatory mechanisms behind genetic risk associations.

Variants in TRIM22 That Affect NOD2 Signaling Are Associated With Very-Early-Onset Inflammatory Bowel Disease.

BACKGROUND & AIMS: Severe forms of inflammatory bowel disease (IBD) that develop in very young children can be caused by variants in a single gene. We performed whole-exome sequence (WES) analysis to identify genetic factors that might cause granulomatous colitis and severe perianal disease, with recurrent bacterial and viral infections, in an infant of consanguineous parents. METHODS: We performed targeted WES analysis of DNA collected from the patient and her parents. We validated our findings by a similar analysis of DNA from 150 patients with very-early-onset IBD not associated with known genetic factors analyzed in Toronto, Oxford, and Munich. We compared gene expression signatures in inflamed vs noninflamed intestinal and rectal tissues collected from patients with treatment-resistant Crohn's disease who participated in a trial of ustekinumab. We performed functional studies of identified variants in primary cells from patients and cell culture. RESULTS: We identified a homozygous variant in the tripartite motif containing 22 gene (TRIM22) of the patient, as well as in 2 patients with a disease similar phenotype. Functional studies showed that the variant disrupted the ability of TRIM22 to regulate nucleotide binding oligomerization domain containing 2 (NOD2)-dependent activation of interferon-beta signaling and nuclear factor-κB. Computational studies demonstrated a correlation between the TRIM22-NOD2 network and signaling pathways and genetic factors associated very early onset and adult-onset IBD. TRIM22 is also associated with antiviral and mycobacterial effectors and markers of inflammation, such as fecal calprotectin, C-reactive protein, and Crohn's disease activity index scores. CONCLUSIONS: In WES and targeted exome sequence analyses of an infant with severe IBD characterized by granulomatous colitis and severe perianal disease, we identified a homozygous variant of TRIM22 that affects the ability of its product to regulate NOD2. Combined computational and functional studies showed that the TRIM22-NOD2 network regulates antiviral and antibacterial signaling pathways that contribute to inflammation. Further study of this network could lead to new disease markers and therapeutic targets for patients with very early and adult-onset IBD.

NDEx, the Network Data Exchange.

Networks are a powerful and flexible methodology for expressing biological knowledge for computation and communication. Network-encoded information can include systematic screens for molecular interactions, biological relationships curated from literature, and outputs from analysis of Big Data. NDEx, the Network Data Exchange (www.ndexbio.org), is an online commons where scientists can upload, share, and publicly distribute networks. Networks in NDEx receive globally unique accession IDs and can be stored for private use, shared in pre-publication collaboration, or released for public access. Standard and novel data formats are accommodated in a flexible storage model. Organizations can use NDEx as a distribution channel for networks they generate or curate. Developers of bioinformatic applications can store and query NDEx networks via a common programmatic interface. NDEx helps expand the role of networks in scientific discourse and facilitates the integration of networks as data in publications. It is a step towards an ecosystem in which networks bearing data, hypotheses, and findings flow easily between scientists.

A disease module in the interactome explains disease heterogeneity, drug response and captures novel pathways and genes in asthma.

Recent advances in genetics have spurred rapid progress towards the systematic identification of genes involved in complex diseases. Still, the detailed understanding of the molecular and physiological mechanisms through which these genes affect disease phenotypes remains a major challenge. Here, we identify the asthma disease module, i.e. the local neighborhood of the interactome whose perturbation is associated with asthma, and validate it for functional and pathophysiological relevance, using both computational and experimental approaches. We find that the asthma disease module is enriched with modest GWAS P-values against the background of random variation, and with differentially expressed genes from normal and asthmatic fibroblast cells treated with an asthma-specific drug. The asthma module also contains immune response mechanisms that are shared with other immune-related disease modules. Further, using diverse omics (genomics, gene-expression, drug response) data, we identify the GAB1 signaling pathway as an important novel modulator in asthma. The wiring diagram of the uncovered asthma module suggests a relatively close link between GAB1 and glucocorticoids (GCs), which we experimentally validate, observing an increase in the level of GAB1 after GC treatment in BEAS-2B bronchial epithelial cells. The siRNA knockdown of GAB1 in the BEAS-2B cell line resulted in a decrease in the NFkB level, suggesting a novel regulatory path of the pro-inflammatory factor NFkB by GAB1 in asthma.

Expression quantitative trait Loci acting across multiple tissues are enriched in inherited risk for coronary artery disease.

BACKGROUND: Despite recent discoveries of new genetic risk factors, the majority of risk for coronary artery disease (CAD) remains elusive. As the most proximal sensor of DNA variation, RNA abundance can help identify subpopulations of genetic variants active in and across tissues mediating CAD risk through gene expression. METHODS AND RESULTS: By generating new genomic data on DNA and RNA samples from the Stockholm Atherosclerosis Gene Expression (STAGE) study, 8156 cis-acting expression quantitative trait loci (eQTLs) for 6450 genes across 7 CAD-relevant tissues were detected. The inherited risk enrichments of tissue-defined sets of these eQTLs were assessed using 2 independent genome-wide association data sets. eQTLs acting across increasing numbers of tissues were found increasingly enriched for CAD risk and resided at regulatory hot spots. The risk enrichment of 42 eQTLs acting across 5 to 6 tissues was particularly high (≤7.3-fold) and confirmed in the combined genome-wide association data from Coronary Artery Disease Genome Wide Replication And Meta-Analysis Consortium. Sixteen of the 42 eQTLs associated with 19 master regulatory genes and 29 downstream gene sets (n>30) were further risk enriched comparable to that of the 153 genome-wide association risk single-nucleotide polymorphisms established for CAD (8.4-fold versus 10-fold). Three gene sets, governed by the master regulators FLYWCH1, PSORSIC3, and G3BP1, segregated the STAGE patients according to extent of CAD, and small interfering RNA targeting of these master regulators affected cholesterol-ester accumulation in foam cells of the THP1 monocytic cell line. CONCLUSIONS: eQTLs acting across multiple tissues are significant carriers of inherited risk for CAD. FLYWCH1, PSORSIC3, and G3BP1 are novel master regulatory genes in CAD that may be suitable targets.

Synthesis of 53 tissue and cell line expression QTL datasets reveals master eQTLs.

BACKGROUND: Gene expression genetic studies in human tissues and cells identify cis- and trans-acting expression quantitative trait loci (eQTLs). These eQTLs provide insights into regulatory mechanisms underlying disease risk. However, few studies systematically characterized eQTL results across cell and tissues types. We synthesized eQTL results from >50 datasets, including new primary data from human brain, peripheral plaque and kidney samples, in order to discover features of human eQTLs. RESULTS: We find a substantial number of robust cis-eQTLs and far fewer trans-eQTLs consistent across tissues. Analysis of 45 full human GWAS scans indicates eQTLs are enriched overall, and above nSNPs, among positive statistical signals in genetic mapping studies, and account for a significant fraction of the strongest human trait effects. Expression QTLs are enriched for gene centricity, higher population allele frequencies, in housekeeping genes, and for coincidence with regulatory features, though there is little evidence of 5' or 3' positional bias. Several regulatory categories are not enriched including microRNAs and their predicted binding sites and long, intergenic non-coding RNAs. Among the most tissue-ubiquitous cis-eQTLs, there is enrichment for genes involved in xenobiotic metabolism and mitochondrial function, suggesting these eQTLs may have adaptive origins. Several strong eQTLs (CDK5RAP2, NBPFs) coincide with regions of reported human lineage selection. The intersection of new kidney and plaque eQTLs with related GWAS suggest possible gene prioritization. For example, butyrophilins are now linked to arterial pathogenesis via multiple genetic and expression studies. Expression QTL and GWAS results are made available as a community resource through the NHLBI GRASP database [http://apps.nhlbi.nih.gov/grasp/]. CONCLUSIONS: Expression QTLs inform the interpretation of human trait variability, and may account for a greater fraction of phenotypic variability than protein-coding variants. The synthesis of available tissue eQTL data highlights many strong cis-eQTLs that may have important biologic roles and could serve as positive controls in future studies. Our results indicate some strong tissue-ubiquitous eQTLs may have adaptive origins in humans. Efforts to expand the genetic, splicing and tissue coverage of known eQTLs will provide further insights into human gene regulation.

Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer's disease.

The genetics of complex disease produce alterations in the molecular interactions of cellular pathways whose collective effect may become clear through the organized structure of molecular networks. To characterize molecular systems associated with late-onset Alzheimer's disease (LOAD), we constructed gene-regulatory networks in 1,647 postmortem brain tissues from LOAD patients and nondemented subjects, and we demonstrate that LOAD reconfigures specific portions of the molecular interaction structure. Through an integrative network-based approach, we rank-ordered these network structures for relevance to LOAD pathology, highlighting an immune- and microglia-specific module that is dominated by genes involved in pathogen phagocytosis, contains TYROBP as a key regulator, and is upregulated in LOAD. Mouse microglia cells overexpressing intact or truncated TYROBP revealed expression changes that significantly overlapped the human brain TYROBP network. Thus the causal network structure is a useful predictor of response to gene perturbations and presents a framework to test models of disease mechanisms underlying LOAD.

Inverse regulation of inflammation and mitochondrial function in adipose tissue defines extreme insulin sensitivity in morbidly obese patients.

Obesity is associated with insulin resistance, a major risk factor for type 2 diabetes and cardiovascular disease. However, not all obese individuals are insulin resistant, which confounds our understanding of the mechanistic link between these conditions. We conducted transcriptome analyses on 835 obese subjects with mean BMI of 48.8, on which we have previously reported genetic associations of gene expression. Here, we selected ~320 nondiabetic (HbA(1c) <7.0) subjects and further stratified the cohort into insulin-resistant versus insulin-sensitive subgroups based on homeostasis model assessment-insulin resistance. An unsupervised informatics analysis revealed that immune response and inflammation-related genes were significantly downregulated in the omental adipose tissue of obese individuals with extreme insulin sensitivity and, to a much lesser extent, in subcutaneous adipose tissue. In contrast, genes related to ß-oxidation and the citric acid cycle were relatively overexpressed in adipose of insulin-sensitive patients. These observations were verified by querying an independent cohort of our published dataset of 37 subjects whose subcutaneous adipose tissue was sampled before and after treatment with thiazolidinediones. Whereas the immune response and inflammation pathway genes were downregulated by thiazolidinedione treatment, ß-oxidation and citric acid cycle genes were upregulated. This work highlights the critical role that omental adipose inflammatory pathways might play in the pathophysiology of insulin resistance, independent of body weight.

The arrestin domain-containing 3 protein regulates body mass and energy expenditure.

A human genome-wide linkage scan for obesity identified a linkage peak on chromosome 5q13-15. Positional cloning revealed an association of a rare haplotype to high body-mass index (BMI) in males but not females. The risk locus contains a single gene, "arrestin domain-containing 3" (ARRDC3), an uncharacterized α-arrestin. Inactivating Arrdc3 in mice led to a striking resistance to obesity, with greater impact on male mice. Mice with decreased ARRDC3 levels were protected from obesity due to increased energy expenditure through increased activity levels and increased thermogenesis of both brown and white adipose tissues. ARRDC3 interacted directly with ß-adrenergic receptors, and loss of ARRDC3 increased the response to ß-adrenergic stimulation in isolated adipose tissue. These results demonstrate that ARRDC3 is a gender-sensitive regulator of obesity and energy expenditure and reveal a surprising diversity for arrestin family protein functions.

Dissecting cis regulation of gene expression in human metabolic tissues.

Complex diseases such as obesity and type II diabetes can result from a failure in multiple organ systems including the central nervous system and tissues involved in partitioning and disposal of nutrients. Studying the genetics of gene expression in tissues that are involved in the development of these diseases can provide insights into how these tissues interact within the context of disease. Expression quantitative trait locus (eQTL) studies identify mRNA expression changes linked to proximal genetic signals (cis eQTLs) that have been shown to affect disease. Given the high impact of recent eQTL studies, it is important to understand what role sample size and environment plays in identification of cis eQTLs. Here we show in a genotyped obese human population that the number of cis eQTLs obey precise scaling laws as a function of sample size in three profiled tissues, i.e. omental adipose, subcutaneous adipose and liver. Also, we show that genes (or transcripts) with cis eQTL associations detected in a small population are detected at approximately 90% rate in the largest population available for our study, indicating that genes with strong cis acting regulatory elements can be identified with relatively high confidence in smaller populations. However, by increasing the sample size we allow for better detection of weaker and more distantly located cis-regulatory elements. Yet, we determined that the number of tissue specific cis eQTLs saturates in a modestly sized cohort while the number of cis eQTLs common to all tissues fails to reach a maximum value. Understanding the power laws that govern the number and specificity of eQTLs detected in different tissues, will allow a better utilization of genetics of gene expression to inform the molecular mechanism underlying complex disease traits.

A survey of the genetics of stomach, liver, and adipose gene expression from a morbidly obese cohort.

To map the genetics of gene expression in metabolically relevant tissues and investigate the diversity of expression SNPs (eSNPs) in multiple tissues from the same individual, we collected four tissues from approximately 1000 patients undergoing Roux-en-Y gastric bypass (RYGB) and clinical traits associated with their weight loss and co-morbidities. We then performed high-throughput genotyping and gene expression profiling and carried out a genome-wide association analyses for more than 100,000 gene expression traits representing four metabolically relevant tissues: liver, omental adipose, subcutaneous adipose, and stomach. We successfully identified 24,531 eSNPs corresponding to about 10,000 distinct genes. This represents the greatest number of eSNPs identified to our knowledge by any study to date and the first study to identify eSNPs from stomach tissue. We then demonstrate how these eSNPs provide a high-quality disease map for each tissue in morbidly obese patients to not only inform genetic associations identified in this cohort, but in previously published genome-wide association studies as well. These data can aid in elucidating the key networks associated with morbid obesity, response to RYGB, and disease as a whole.

Liver and adipose expression associated SNPs are enriched for association to type 2 diabetes.

Genome-wide association studies (GWAS) have demonstrated the ability to identify the strongest causal common variants in complex human diseases. However, to date, the massive data generated from GWAS have not been maximally explored to identify true associations that fail to meet the stringent level of association required to achieve genome-wide significance. Genetics of gene expression (GGE) studies have shown promise towards identifying DNA variations associated with disease and providing a path to functionally characterize findings from GWAS. Here, we present the first empiric study to systematically characterize the set of single nucleotide polymorphisms associated with expression (eSNPs) in liver, subcutaneous fat, and omental fat tissues, demonstrating these eSNPs are significantly more enriched for SNPs that associate with type 2 diabetes (T2D) in three large-scale GWAS than a matched set of randomly selected SNPs. This enrichment for T2D association increases as we restrict to eSNPs that correspond to genes comprising gene networks constructed from adipose gene expression data isolated from a mouse population segregating a T2D phenotype. Finally, by restricting to eSNPs corresponding to genes comprising an adipose subnetwork strongly predicted as causal for T2D, we dramatically increased the enrichment for SNPs associated with T2D and were able to identify a functionally related set of diabetes susceptibility genes. We identified and validated malic enzyme 1 (Me1) as a key regulator of this T2D subnetwork in mouse and provided support for the association of this gene to T2D in humans. This integration of eSNPs and networks provides a novel approach to identify disease susceptibility networks rather than the single SNPs or genes traditionally identified through GWAS, thereby extracting additional value from the wealth of data currently being generated by GWAS.