SUMMARY 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 1647 post-mortem brain tissues from LOAD patients and non-demented subjects, and 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 dominated by genes involved in pathogen phagocytosis, containing TYROBP as a key regulator and up-regulated in LOAD. Mouse microglia cells over-expressing 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 novel framework to test models of disease mechanisms underlying LOAD.

INTRODUCTION: Immunotherapy targeting the programmed cell death protein–1 (PD-1) axis elicits durable antitumor responses in multiple cancer types. However, clinical responses vary, and biomarkers predictive of response may help to identify patients who will derive the greatest therapeutic benefit. Clinically validated biomarkers predictive of response to the anti–PD-1 monoclonal antibody pembrolizumab include PD-1 ligand 1 (PD-L1) expression in specific cancers and high microsatellite instability (MSI-H) regardless of tumor type. Tumor mutational burden (TMB) and T cell–inflamed gene expression profile (GEP) are emerging predictive biomarkers for pembrolizumab. Both PD-L1 and GEP are inflammatory biomarkers indicative of a T cell–inflamed tumor microenvironment (TME), whereas TMB and MSI-H are indirect measures of tumor antigenicity generated by somatic tumor mutations. However, the relationship between these two categories of biomarkers is not well characterized. RATIONALE: This study assessed the potential for TMB and a T cell–inflamed GEP to jointly predict clinical response to pembrolizumab in >300 patient samples with advanced solid tumors and melanoma across 22 tumor types from four KEYNOTE clinical trials. To assess the individual and joint clinical utility of TMB and GEP, patients were stratified in four biomarker–defined clinical response groups [GEP low and TMB low (GEPlo TMBlo), GEP low and TMB high (GEPlo TMBhi), GEPhi TMBlo, and GEPhi TMBhi] based on predefined cutoffs for TMB and GEP. These patient–defined biomarker groups were further used to guide transcriptome and exome analyses of tumors in a large molecular database [The Cancer Genome Atlas (TCGA)] (n = 6384 tumors) to identify targetable patterns of biology that may modulate response and resistance. RESULTS: TMB and GEP exhibited only modest correlation and were independently predictive of response across the KEYNOTE clinical datasets. We found that objective response rates were strongest in patients with GEPhi TMBhi (37 to 57%), moderate in those with GEPhi TMBlo (12 to 35%) and GEPlo TMBhi (11 to 42%), and reduced or absent in those with GEPlo TMBlo (0 to 9%) (see the figure). Additionally, longer progression–free survival times were seen in patients with higher levels of both TMB and GEP. Findings were comparable when TMB and PD-L1 expression were jointly assessed. Within TCGA database, GEP and TMB again had a low correlation, demonstrating the potential to jointly stratify transcriptomic and genomic features across cancer types. Specific gene expression patterns reflective of TME biology showed significant associations with TMB, GEP, or both. In particular, gene set enrichment analysis identified proliferative and stromal, myeloid, and vascular biology corresponding to specific TMB-defined subgroups within GEPhi tumors. In TMBhi tumors, indication-dependent somatic DNA alterations in key cancer driver genes showed a strong negative association with GEP. CONCLUSION: This analysis shows that TMB and inflammatory biomarkers (T cell–in...

A key goal of biomedical research is to elucidate the complex network of gene interactions underlying complex traits such as common human diseases. Here we detail a multistep procedure for identifying potential key drivers of complex traits that integrates DNA-variation and gene-expression data with other complex trait data in segregating mouse populations. Ordering gene expression traits relative to one another and relative to other complex traits is achieved by systematically testing whether variations in DNA that lead to variations in relative transcript abundances statistically support an independent, causative or reactive function relative to the complex traits under consideration. We show that this approach can predict transcriptional responses to single gene-perturbation experiments using gene-expression data in the context of a segregating mouse population. We also demonstrate the utility of this approach by identifying and experimentally validating the involvement of three new genes in susceptibility to obesity.In the past few years, gene-expression microarrays and other general molecular profiling technologies have been applied to a wide range of biological problems and have contributed to discoveries about the complex network of biochemical processes underlying living Correspondence should be addressed to E.E.S. (eric_schadt@merck.com). Note: Supplementary information is available on the Nature Genetics website. COMPETING INTERESTS STATEMENTThe authors declare that they have no competing financial interests. NIH Public Access Author ManuscriptNat Genet. Author manuscript; available in PMC 2010 March 18. Published in final edited form as:Nat Genet. 2005 July ; 37(7): 710-717. doi:10.1038/ng1589. NIH-PA Author ManuscriptNIH-PA Author Manuscript NIH-PA Author Manuscript systems 1 , common human diseases 2,3 and gene discovery and structure determination [4][5][6] . Microarrays have also helped to identify biomarkers 7 , disease subtypes 3,8,9 and mechanisms of toxicity 10 and, more recently, to elucidate the genetics of gene expression in human populations 11,12 and to reconstruct gene networks by integrating gene-expression and genetic data 13 . The use of molecular profiling technologies as tools to identify genes underlying common, polygenic diseases has been less successful. Hundreds or even thousands of genes whose expression changes are associated with disease traits have been identified, but determining which of the genes cause disease rather than respond to the disease state has proven difficult.Microarray data have recently been combined with other experimental approaches to facilitate identification of key mechanistic drivers of complex traits 3,[13][14][15][16][17] . One such technique involves treating relative transcript abundances as quantitative traits in segregating populations. In this method, chromosomal regions that control the level of expression of a particular gene are mapped as expression quantitative trait loci (eQTLs). Gene-expression QTLs that contain the gene encoding t...

To survey hepatitis B virus (HBV) integration in liver cancer genomes, we conducted massively parallel sequencing of 81 HBV-positive and 7 HBV-negative hepatocellular carcinomas (HCCs) and adjacent normal tissues. We found that HBV integration is observed more frequently in the tumors (86.4%) than in adjacent liver tissues (30.7%). Copy-number variations (CNVs) were significantly increased at HBV breakpoint locations where chromosomal instability was likely induced. Approximately 40% of HBV breakpoints within the HBV genome were located within a 1,800-bp region where the viral enhancer, X gene and core gene are located. We also identified recurrent HBV integration events (in ≥ 4 HCCs) that were validated by RNA sequencing (RNA-seq) and Sanger sequencing at the known and putative cancer-related TERT, MLL4 and CCNE1 genes, which showed upregulated gene expression in tumor versus normal tissue. We also report evidence that suggests that the number of HBV integrations is associated with patient survival.

Identifying variations in DNA that increase susceptibility to disease is one of the primary aims of genetic studies using a forward genetics approach. However, identification of disease-susceptibility genes by means of such studies provides limited functional information on how genes lead to disease. In fact, in most cases there is an absence of functional information altogether, preventing a definitive identification of the susceptibility gene or genes. Here we develop an alternative to the classic forward genetics approach for dissecting complex disease traits where, instead of identifying susceptibility genes directly affected by variations in DNA, we identify gene networks that are perturbed by susceptibility loci and that in turn lead to disease. Application of this method to liver and adipose gene expression data generated from a segregating mouse population results in the identification of a macrophage-enriched network supported as having a causal relationship with disease traits associated with metabolic syndrome. Three genes in this network, lipoprotein lipase (Lpl), lactamase beta (Lactb) and protein phosphatase 1-like (Ppm1l), are validated as previously unknown obesity genes, strengthening the association between this network and metabolic disease traits. Our analysis provides direct experimental support that complex traits such as obesity are emergent properties of molecular networks that are modulated by complex genetic loci and environmental factors.

Hepatocellular carcinoma (HCC) is one of the most deadly cancers worldwide and has no effective treatment, yet the molecular basis of hepatocarcinogenesis remains largely unknown. Here we report findings from a whole-genome sequencing (WGS) study of 88 matched HCC tumor/normal pairs, 81 of which are Hepatitis B virus (HBV) positive, seeking to identify genetically altered genes and pathways implicated in HBV-associated HCC. We find beta-catenin to be the most frequently mutated oncogene (15.9%) and TP53 the most frequently mutated tumor suppressor (35.2%). The Wnt/beta-catenin and JAK/STAT pathways, altered in 62.5% and 45.5% of cases, respectively, are likely to act as two major oncogenic drivers in HCC. This study also identifies several prevalent and potentially actionable mutations, including activating mutations of Janus kinase 1 (JAK1), in 9.1% of patients and provides a path toward therapeutic intervention of the disease.

Summary Tumor associated macrophages (TAM) contribute to all aspects of tumor progression. Use of CSF1R inhibitors to target TAM is therapeutically appealing, but has had very limited antitumor effects. Here, we have identified the mechanism that limited the effect of CSF1R targeted therapy. We demonstrated that carcinoma associated fibroblasts (CAF) are major sources of chemokines that recruit granulocytes to tumors. CSF1 produced by tumor cells caused HDAC2-mediated down-regulation of granulocyte-specific chemokine expression in CAF, which limited migration of these cells to tumors. Treatment with CSF1R inhibitors disrupted this cross talk and triggered a profound increase in granulocyte recruitment to tumors. Combining CSF1R inhibitor with a CXCR2 antagonist blocked granulocyte infiltration of tumors and showed strong anti-tumor effects.

Genome-wide association studies (GWAS) have identified loci reproducibly associated with pulmonary diseases; however, the molecular mechanism underlying these associations are largely unknown. The objectives of this study were to discover genetic variants affecting gene expression in human lung tissue, to refine susceptibility loci for asthma identified in GWAS studies, and to use the genetics of gene expression and network analyses to find key molecular drivers of asthma. We performed a genome-wide search for expression quantitative trait loci (eQTL) in 1,111 human lung samples. The lung eQTL dataset was then used to inform asthma genetic studies reported in the literature. The top ranked lung eQTLs were integrated with the GWAS on asthma reported by the GABRIEL consortium to generate a Bayesian gene expression network for discovery of novel molecular pathways underpinning asthma. We detected 17,178 cis- and 593 trans- lung eQTLs, which can be used to explore the functional consequences of loci associated with lung diseases and traits. Some strong eQTLs are also asthma susceptibility loci. For example, rs3859192 on chr17q21 is robustly associated with the mRNA levels of GSDMA (P = 3.55×10−151). The genetic-gene expression network identified the SOCS3 pathway as one of the key drivers of asthma. The eQTLs and gene networks identified in this study are powerful tools for elucidating the causal mechanisms underlying pulmonary disease. This data resource offers much-needed support to pinpoint the causal genes and characterize the molecular function of gene variants associated with lung diseases.