Variations in DNA elucidate molecular networks that cause disease.

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.

Pancreatic Cancer Action Network's SPARK: A Cloud-Based Patient Health Data and Analytics Platform for Pancreatic Cancer.

PURPOSE: Pancreatic cancer currently holds the position of third deadliest cancer in the United States and the 5-year survival rate is among the lowest for major cancers at just 12%. Thus, continued research efforts to better understand the clinical and molecular underpinnings of pancreatic cancer are critical to developing both early detection methodologies as well as improved therapeutic options. This study introduces Pancreatic Cancer Action Network's (PanCAN's) SPARK, a cloud-based data and analytics platform that integrates patient health data from the PanCAN's research initiatives and aims to accelerate pancreatic cancer research by making real-world patient health data and analysis tools easier to access and use. MATERIALS AND METHODS: The SPARK platform integrates clinical, molecular, multiomic, imaging, and patient-reported data generated from PanCAN's research initiatives. The platform is built on a cloud-based infrastructure powered by Velsera. Cohort exploration and browser capabilities are built using Velsera ARIA, a specialized product for leveraging clinicogenomic data to build cohorts, query variant information, and drive downstream association analyses. Data science and analytic capabilities are also built into the platform allowing researchers to perform simple to complex analysis. RESULTS: Version 1 of the SPARK platform was released to pilot users, who represented diverse end users, including molecular biologists, clinicians, and bioinformaticians. Included in the pilot release of SPARK are deidentified clinical (including treatment and outcomes data), molecular, multiomic, and whole-slide pathology images for over 600 patients enrolled in PanCAN's Know Your Tumor molecular profiling service. CONCLUSION: The pilot release of the SPARK platform introduces qualified researchers to PanCAN real-world patient health data and analytical resources in a centralized location.

Impact of KRAS mutations and co-mutations on clinical outcomes in pancreatic ductal adenocarcinoma.

The relevance of KRAS mutation alleles to clinical outcome remains inconclusive in pancreatic adenocarcinoma (PDAC). We conducted a retrospective study of 803 patients with PDAC (42% with metastatic disease) at MD Anderson Cancer Center. Overall survival (OS) analysis demonstrated that KRAS mutation status and subtypes were prognostic (p < 0.001). Relative to patients with KRAS wildtype tumors (median OS 38 months), patients with KRASG12R had a similar OS (median 34 months), while patients with KRASQ61 and KRASG12D mutated tumors had shorter OS (median 20 months [HR: 1.9, 95% CI 1.2-3.0, p = 0.006] and 22 months [HR: 1.7, 95% CI 1.3-2.3, p < 0.001], respectively). There was enrichment of KRASG12D mutation in metastatic tumors (34% vs 24%, OR: 1.7, 95% CI 1.2-2.4, p = 0.001) and enrichment of KRASG12R in well and moderately differentiated tumors (14% vs 9%, OR: 1.7, 95% CI 1.05-2.99, p = 0.04). Similar findings were observed in the external validation cohort (PanCAN's Know Your Tumor® dataset, n = 408).

Impact of KRAS Mutations and Co-mutations on Clinical Outcomes in Pancreatic Ductal Adenocarcinoma.

The relevance of KRAS mutation alleles to clinical outcome remains inconclusive in pancreatic adenocarcinoma (PDAC). We conducted a retrospective study of 803 PDAC patients (42% with metastatic disease) at MD Anderson Cancer Center. Overall survival (OS) analysis demonstrated that KRAS mutation status and subtypes were prognostic (p<0.001). Relative to patients with KRAS wildtype tumors (median OS 38 months), patients with KRASG12R had a similar OS (median 34 months), while patients with KRASQ61 and KRASG12D mutated tumors had shorter OS (median 20 months [HR: 1.9, 95% CI 1.2-3.0, p=0.006] and 22 months [HR: 1.7, 95% CI 1.3-2.3, p<0.001], respectively). There was enrichment of KRASG12D mutation in metastatic tumors (34% vs 24%, OR: 1.7, 95% CI 1.2-2.4, p=0.001) and enrichment of KRASG12R in well and moderately differentiated tumors (14% vs 9%, OR: 1.7, 95% CI 1.05-2.99, p=0.04). Similar findings were observed in the external validation cohort (PanCAN's Know Your Tumor® dataset, n=408).

High-resolution mapping of gene expression using association in an outbred mouse stock.

Quantitative trait locus (QTL) analysis is a powerful tool for mapping genes for complex traits in mice, but its utility is limited by poor resolution. A promising mapping approach is association analysis in outbred stocks or different inbred strains. As a proof of concept for the association approach, we applied whole-genome association analysis to hepatic gene expression traits in an outbred mouse population, the MF1 stock, and replicated expression QTL (eQTL) identified in previous studies of F2 intercross mice. We found that the mapping resolution of these eQTL was significantly greater in the outbred population. Through an example, we also showed how this precise mapping can be used to resolve previously identified loci (in intercross studies), which affect many different transcript levels (known as eQTL "hotspots"), into distinct regions. Our results also highlight the importance of correcting for population structure in whole-genome association studies in the outbred stock.

Leveraging patient-reported outcomes (PROs) in patients with pancreatic cancer: The Pancreatic Cancer Action Network (PanCAN) online patient registry experience.

BACKGROUND: The Pancreatic Cancer Action Network (PanCAN) Patient Registry is an online, pancreatic cancer-specific, global registry enabling patients to self-report sociodemographics, disease/management characteristics, and patient-reported outcomes (PROs). We sought to describe the creation, user experience, and research potential of the PanCAN Registry. METHODS: We obtained data to describe (1) the creation of the Registry (questionnaire development, marketing efforts, and regulatory considerations); (2) the user experience (user characteristics and interactions with the registry following inception); and (3) the research potential of the registry (comparing PROs and treatment patterns by age [±65 years] and treatment site [community or academic] for users with de novo metastatic disease). RESULTS: The Registry was conceived as part of PanCAN's strategic plan for a personalized therapy initiative. PanCAN staff and disease expert consultants developed questionnaires hosted on an electronic PRO platform. Users had the option to include their data in research efforts, and the Registry platform received institutional review board approval. From 7/2015 to 12/2020, 2187 patients visited the registry and 1697 (77.6%) completed at least one survey (median age = 64 years [range: 24-90], 47.9% women, 88.7% White, 34.0% metastatic disease). Among patients with metastatic disease (N = 567), 46.0% were ≥65 years old and 67.5% received treatment at community sites. Patients ≥65 years reported feeling less hopeful about the treatment plan (12.4% vs. 24.3%, p = 0.003), and patients treated at community sites reported more frequent treatment breaks of >2 weeks (58.2% vs. 28.1%, p < 0.001). CONCLUSIONS: Our findings demonstrate the feasibility, usability, and research potential of an online PRO registry for patients with cancer. This description of the PanCAN Registry should inform future registry-building efforts to facilitate standardized PRO reporting and provide a valuable research database. CLINICAL TRIAL REGISTRATION NUMBER: Not applicable.

Integrating genetic and network analysis to characterize genes related to mouse weight.

Systems biology approaches that are based on the genetics of gene expression have been fruitful in identifying genetic regulatory loci related to complex traits. We use microarray and genetic marker data from an F2 mouse intercross to examine the large-scale organization of the gene co-expression network in liver, and annotate several gene modules in terms of 22 physiological traits. We identify chromosomal loci (referred to as module quantitative trait loci, mQTL) that perturb the modules and describe a novel approach that integrates network properties with genetic marker information to model gene/trait relationships. Specifically, using the mQTL and the intramodular connectivity of a body weight-related module, we describe which factors determine the relationship between gene expression profiles and weight. Our approach results in the identification of genetic targets that influence gene modules (pathways) that are related to the clinical phenotypes of interest.

An integrative genetics approach to identify candidate genes regulating BMD: combining linkage, gene expression, and association.

Numerous quantitative trait loci (QTLs) affecting bone traits have been identified in the mouse; however, few of the underlying genes have been discovered. To improve the process of transitioning from QTL to gene, we describe an integrative genetics approach, which combines linkage analysis, expression QTL (eQTL) mapping, causality modeling, and genetic association in outbred mice. In C57BL/6J x C3H/HeJ (BXH) F(2) mice, nine QTLs regulating femoral BMD were identified. To select candidate genes from within each QTL region, microarray gene expression profiles from individual F(2) mice were used to identify 148 genes whose expression was correlated with BMD and regulated by local eQTLs. Many of the genes that were the most highly correlated with BMD have been previously shown to modulate bone mass or skeletal development. Candidates were further prioritized by determining whether their expression was predicted to underlie variation in BMD. Using network edge orienting (NEO), a causality modeling algorithm, 18 of the 148 candidates were predicted to be causally related to differences in BMD. To fine-map QTLs, markers in outbred MF1 mice were tested for association with BMD. Three chromosome 11 SNPs were identified that were associated with BMD within the Bmd11 QTL. Finally, our approach provides strong support for Wnt9a, Rasd1, or both underlying Bmd11. Integration of multiple genetic and genomic data sets can substantially improve the efficiency of QTL fine-mapping and candidate gene identification.

Cis-acting expression quantitative trait loci in mice.

We previously reported the analysis of genome-wide expression profiles and various diabetes-related traits in a segregating cross between inbred mouse strains C57BL/6J (B6) and DBA/2J (DBA). By considering transcript levels as quantitative traits, we identified several thousand expression quantitative trait loci (eQTL) with LOD score >4.3. We now experimentally address the problem of multiple comparisons by estimating the fraction of false-positive eQTL that are under cis-acting regulation. For this, we have utilized a classic cis-trans test with (B6 x DBA)F(1) mice to determine the relative levels of transcripts from the B6 and DBA alleles. The results suggest that at least 64% of cis-acting eQTL with LOD >4.3 are true positives, while the remaining 36% could not be confirmed as truly cis-acting. Moreover, we find that >96% of apparent cis-acting eQTL occur in regions that do not share SNP haplotypes. Cis-acting eQTL serve as an important new resource for the identification of positional candidates in QTL studies in mice. Also, we use the analysis of the correlation structures between genotypes, gene expression traits, and phenotypic traits to further characterize genes expressed in liver that are under cis-acting control, and highlight the advantages and disadvantages of integrating genetics and gene expression data in segregating populations.

Genomic analysis of metabolic pathway gene expression in mice.

BACKGROUND: A segregating population of (C57BL/6J x DBA/2J)F2 intercross mice was studied for obesity-related traits and for global gene expression in liver. Quantitative trait locus analyses were applied to the subcutaneous fat-mass trait and all gene-expression data. These data were then used to identify gene sets that are differentially perturbed in lean and obese mice. RESULTS: We integrated global gene-expression data with phenotypic and genetic segregation analyses to evaluate metabolic pathways associated with obesity. Using two approaches we identified 13 metabolic pathways whose genes are coordinately regulated in association with obesity. Four genomic regions on chromosomes 3, 6, 16, and 19 were found to control the coordinated expression of these pathways. Using criteria that included trait correlation, differential gene expression, and linkage to genomic regions, we identified novel genes potentially associated with the identified pathways. CONCLUSION: This study demonstrates that genetic and gene-expression data can be integrated to identify pathways associated with clinical traits and their underlying genetic determinants.

A genome-wide set of congenic mouse strains derived from DBA/2J on a C57BL/6J background.

In the analysis of complex traits, congenic strains are powerful tools because they allow characterization of a single locus in the absence of genetic variation throughout the remainder of the genome. Here, we report the construction and initial characterization of a genome-wide panel of congenic strains derived from the donor strain DBA/2J on the background strain C57BL/6J. For many strains, we have carried out high-density SNP genotyping to precisely map the congenic interval and to identify any contaminating regions. Certain strains exhibit striking variation in litter size and in the ratio of females to males. We illustrate the utility of the set by "Mendelizing" the complex trait of myocardial calcification. These 65 strains cover more than 95% of the autosomal genome and should facilitate the analysis of the many genetic trait differences that have been reported between these parental strains.

Thematic review series: The pathogenesis of atherosclerosis. Toward a biological network for atherosclerosis.

The goal of systems biology is to define all of the elements present in a given system and to create an interaction network between these components so that the behavior of the system, as a whole and in parts, can be explained under specified conditions. The elements constituting the network that influences the development of atherosclerosis could be genes, pathways, transcript levels, proteins, or physiologic traits. In this review, we discuss how the integration of genetics and technologies such as transcriptomics and proteomics, combined with mathematical modeling, may lead to an understanding of such networks.