Colocalization of GWAS and eQTL signals at loci with multiple signals identifies additional candidate genes for body fat distribution.

Integration of genome-wide association study (GWAS) signals with expression quantitative trait loci (eQTL) studies enables identification of candidate genes. However, evaluating whether nearby signals may share causal variants, termed colocalization, is affected by the presence of allelic heterogeneity, different variants at the same locus impacting the same phenotype. We previously identified eQTL in subcutaneous adipose tissue from 770 participants in the Metabolic Syndrome in Men (METSIM) study and detected 15 eQTL signals that colocalized with GWAS signals for waist-hip ratio adjusted for body mass index (WHRadjBMI) from the Genetic Investigation of Anthropometric Traits consortium. Here, we reevaluated evidence of colocalization using two approaches, conditional analysis and the Bayesian test COLOC, and show that providing COLOC with approximate conditional summary statistics at multi-signal GWAS loci can reconcile disagreements in colocalization classification between the two tests. Next, we performed conditional analysis on the METSIM subcutaneous adipose tissue data to identify conditionally distinct or secondary eQTL signals. We used the two approaches to test for colocalization with WHRadjBMI GWAS signals and evaluated the differences in colocalization classification between the two tests. Through these analyses, we identified four GWAS signals colocalized with secondary eQTL signals for FAM13A, SSR3, GRB14 and FMO1. Thus, at loci with multiple eQTL and/or GWAS signals, analyzing each signal independently enabled additional candidate genes to be identified.

Genetic Regulation of Adipose Gene Expression and Cardio-Metabolic Traits.

Subcutaneous adipose tissue stores excess lipids and maintains energy balance. We performed expression quantitative trait locus (eQTL) analyses by using abdominal subcutaneous adipose tissue of 770 extensively phenotyped participants of the METSIM study. We identified cis-eQTLs for 12,400 genes at a 1% false-discovery rate. Among an approximately 680 known genome-wide association study (GWAS) loci for cardio-metabolic traits, we identified 140 coincident cis-eQTLs at 109 GWAS loci, including 93 eQTLs not previously described. At 49 of these 140 eQTLs, gene expression was nominally associated (p < 0.05) with levels of the GWAS trait. The size of our dataset enabled identification of five loci associated (p < 5 × 10-8) with at least five genes located >5 Mb away. These trans-eQTL signals confirmed and extended the previously reported KLF14-mediated network to 55 target genes, validated the CIITA regulation of class II MHC genes, and identified ZNF800 as a candidate master regulator. Finally, we observed similar expression-clinical trait correlations of genes associated with GWAS loci in both humans and a panel of genetically diverse mice. These results provide candidate genes for further investigation of their potential roles in adipose biology and in regulating cardio-metabolic traits.

Genome-wide transcriptome analysis identifies novel dysregulated genes implicated in Alzheimer's pathology.

INTRODUCTION: Abnormal gene expression patterns may contribute to the onset and progression of late-onset Alzheimer's disease (LOAD). METHODS: We performed transcriptome-wide meta-analysis (N = 1440) of blood-based microarray gene expression profiles as well as neuroimaging and cerebrospinal fluid (CSF) endophenotype analysis. RESULTS: We identified and replicated five genes (CREB5, CD46, TMBIM6, IRAK3, and RPAIN) as significantly dysregulated in LOAD. The most significantly altered gene, CREB5, was also associated with brain atrophy and increased amyloid beta (Aß) accumulation, especially in the entorhinal cortex region. cis-expression quantitative trait loci mapping analysis of CREB5 detected five significant associations (P < 5 × 10-8 ), where rs56388170 (most significant) was also significantly associated with global cortical Aß deposition measured by [18 F]Florbetapir positron emission tomography and CSF Aß1-42 . DISCUSSION: RNA from peripheral blood indicated a differential gene expression pattern in LOAD. Genes identified have been implicated in biological processes relevant to Alzheimer's disease. CREB, in particular, plays a key role in nervous system development, cell survival, plasticity, and learning and memory.

Hypoxia and transforming growth factor-beta1 pathway activation promote Chinese Hamster Ovary cell aggregation.

Suspension cultivation is the preferred mode of operation for the large-scale production of many biologics. Chinese Hamster Ovary (CHO) cells are anchorage-dependent in origin, but they have been widely adapted to suspension culture. In suspension culture, formation of CHO cell aggregates is a common phenomenon and compromises cell culture performance in multiple ways. To better understand the underlying mechanisms that regulate cell aggregation, we utilized CHO-specific transcriptome profiling as a screening tool and demonstrated that many genes encoding extracellular matrix (ECM) proteins were upregulated in the cultures with increased cell aggregation. Significantly, hypoxia was identified to be a cause for promoting CHO cell aggregation, and transforming growth factor beta1 (TGFß1) pathway activation served as an intermediate step mediating this biological cascade. These transcriptomics findings were confirmed by cell culture experiments, and it was further demonstrated that adding recombinant TGFß1 to the culture significantly increased ECM protein fibronectin expression and cell aggregation. The results of this study emphasize the importance of adequate mixing and oxygen supply for suspension cultures from a new angle, and regulating the TGFß1 pathway is proposed as a new strategy for mitigating cell aggregation to improve cell culture performance.

Genetic regulation of mouse liver metabolite levels.

We profiled and analyzed 283 metabolites representing eight major classes of molecules including Lipids, Carbohydrates, Amino Acids, Peptides, Xenobiotics, Vitamins and Cofactors, Energy Metabolism, and Nucleotides in mouse liver of 104 inbred and recombinant inbred strains. We find that metabolites exhibit a wide range of variation, as has been previously observed with metabolites in blood serum. Using genome-wide association analysis, we mapped 40% of the quantified metabolites to at least one locus in the genome and for 75% of the loci mapped we identified at least one candidate gene by local expression QTL analysis of the transcripts. Moreover, we validated 2 of 3 of the significant loci examined by adenoviral overexpression of the genes in mice. In our GWAS results, we find that at significant loci the peak markers explained on average between 20 and 40% of variation in the metabolites. Moreover, 39% of loci found to be regulating liver metabolites in mice were also found in human GWAS results for serum metabolites, providing support for similarity in genetic regulation of metabolites between mice and human. We also integrated the metabolomic data with transcriptomic and clinical phenotypic data to evaluate the extent of co-variation across various biological scales.

HuR is required for IL-17-induced Act1-mediated CXCL1 and CXCL5 mRNA stabilization.

IL-17, a major inflammatory cytokine plays a critical role in the pathogenesis of many autoimmune inflammatory diseases. In this study, we report a new function of RNA-binding protein HuR in IL-17-induced Act1-mediated chemokine mRNA stabilization. HuR deficiency markedly reduced IL-17-induced chemokine expression due to increased mRNA decay. Act1-mediated HuR polyubiquitination was required for the binding of HuR to CXCL1 mRNA, leading to mRNA stabilization. Although IL-17 induced the coshift of Act1 and HuR to the polysomal fractions in a sucrose gradient, HuR deficiency reduced the ratio of translation-active/translation-inactive IL-17-induced chemokine mRNAs. Furthermore, HuR deletion in distal lung epithelium attenuated IL-17-induced neutrophilia. In summary, HuR functions to couple receptor-proximal signaling to posttranscriptional machinery, contributing to IL-17-induced inflammation.

Systems genetics analysis of gene-by-environment interactions in human cells.

Gene by environment (GxE) interactions are clearly important in many human diseases, but they have proven to be difficult to study on a molecular level. We report genetic analysis of thousands of transcript abundance traits in human primary endothelial cell (EC) lines in response to proinflammatory oxidized phospholipids implicated in cardiovascular disease. Of the 59 most regulated transcripts, approximately one-third showed evidence of GxE interactions. The interactions resulted primarily from effects of distal-, trans-acting loci, but a striking example of a local-GxE interaction was also observed for FGD6. Some of the distal interactions were validated by siRNA knockdown experiments, including a locus involved in the regulation of multiple transcripts involved in the ER stress pathway. Our findings add to the understanding of the overall architecture of complex human traits and are consistent with the possibility that GxE interactions are responsible, in part, for the failure of association studies to more fully explain common disease variation.

A high-resolution association mapping panel for the dissection of complex traits in mice.

Systems genetics relies on common genetic variants to elucidate biologic networks contributing to complex disease-related phenotypes. Mice are ideal model organisms for such approaches, but linkage analysis has been only modestly successful due to low mapping resolution. Association analysis in mice has the potential of much better resolution, but it is confounded by population structure and inadequate power to map traits that explain less than 10% of the variance, typical of mouse quantitative trait loci (QTL). We report a novel strategy for association mapping that combines classic inbred strains for mapping resolution and recombinant inbred strains for mapping power. Using a mixed model algorithm to correct for population structure, we validate the approach by mapping over 2500 cis-expression QTL with a resolution an order of magnitude narrower than traditional QTL analysis. We also report the fine mapping of metabolic traits such as plasma lipids. This resource, termed the Hybrid Mouse Diversity Panel, makes possible the integration of multiple data sets and should prove useful for systems-based approaches to complex traits and studies of gene-by-environment interactions.

Identification and validation of G protein-coupled receptors modulating flow-dependent signaling pathways in vascular endothelial cells.

Vascular endothelial cells are exposed to mechanical forces due to their presence at the interface between the vessel wall and flowing blood. The patterns of these mechanical forces (laminar vs. turbulent) regulate endothelial cell function and play an important role in determining endothelial phenotype and ultimately cardiovascular health. One of the key transcriptional mediators of the positive effects of laminar flow patterns on endothelial cell phenotype is the zinc-finger transcription factor, krüppel-like factor 2 (KLF2). Given its importance in maintaining a healthy endothelium, we sought to identify endothelial regulators of the KLF2 transcriptional program as potential new therapeutic approaches to treating cardiovascular disease. Using an approach that utilized both bioinformatics and targeted gene knockdown, we identified endothelial GPCRs capable of modulating KLF2 expression. Genetic screening using siRNAs directed to these GPCRs identified 12 potential GPCR targets that could modulate the KLF2 program, including a subset capable of regulating flow-induced KLF2 expression in primary endothelial cells. Among these targets, we describe the ability of several GPCRs (GPR116, SSTR3, GPR101, LGR4) to affect KLF2 transcriptional activation. We also identify these targets as potential validated targets for the development of novel treatments targeting the endothelium. Finally, we highlight the initiation of drug discovery efforts for LGR4 and report the identification of the first known synthetic ligands to this receptor as a proof-of-concept for pathway-directed phenotypic screening to identify novel drug targets.

Biochemical and transcriptional profiling to triage additional activities in a series of IGF-1R/IR inhibitors.

Therapeutic development of a targeted agent involves a series of decisions over additional activities that may be ignored, eliminated or pursued. This paper details the concurrent application of two methods that provide a spectrum of information about the biological activity of a compound: biochemical profiling on a large panel of kinase assays and transcriptional profiling of mRNA responses. Our mRNA profiling studies used a full dose range, identifying subsets of transcriptional responses with differing EC(50)s which may reflect distinct targets. Profiling data allowed prioritization for validation in xenograft models, generated testable hypotheses for active compounds, and informed decisions on the general utility of the series.

Differential Changes in ACPA Fine Specificity and Gene Expression in a Randomized Trial of Abatacept and Adalimumab in Rheumatoid Arthritis.

INTRODUCTION: The biologics abatacept and adalimumab have different mechanisms of action (MoAs). We analyzed data from patients with rheumatoid arthritis treated in AMPLE (NCT00929864) to explore the pharmacodynamic effects of abatacept or adalimumab on anti-citrullinated protein antibodies (ACPAs) and gene expression. METHODS: AMPLE was a phase IIIb, 2-year, randomized, head-to-head trial of abatacept versus adalimumab. Post hoc analyses of baseline anti-cyclic citrullinated peptide-2 (anti-CCP2, an ACPA surrogate) positive (+) status and ACPA fine-specificity profiles over time, as well as transcriptional profiling (peripheral whole blood), were performed. RESULTS: Of 646 patients treated (abatacept, n = 318; adalimumab, n = 328), ACPA and gene expression data were available from 508 and 566 patients, respectively. In anti-CCP2+ patients (n = 388), baseline fine specificities for most ACPAs were highly correlated; over 2 years, levels decreased with abatacept but not adalimumab. By year 2, expression of genes associated with T cell co-stimulation and antibody production was lower for abatacept versus adalimumab; expression of genes associated with proinflammatory signaling was lower for adalimumab versus abatacept. Treatment modulated the expression of T- and B-cell gene signatures, with differences in CD8+ T cells, activated T cells, plasma cells, B cells, natural killer cells (all lower with abatacept versus adalimumab), and polymorphonuclear leukocytes (higher with abatacept versus adalimumab). CONCLUSIONS: In AMPLE, despite similar clinical outcomes, data showed that pharmacodynamic/genetic changes after 2 years of abatacept or adalimumab were consistent with drug MoAs. Further assessment of the relationship between such changes and clinical outcomes, including prediction of response, is warranted. TRIAL REGISTRATION: ClinicalTrials.gov identifier, NCT00929864.

Transcriptional profiling of the dose response: a more powerful approach for characterizing drug activities.

The dose response curve is the gold standard for measuring the effect of a drug treatment, but is rarely used in genomic scale transcriptional profiling due to perceived obstacles of cost and analysis. One barrier to examining transcriptional dose responses is that existing methods for microarray data analysis can identify patterns, but provide no quantitative pharmacological information. We developed analytical methods that identify transcripts responsive to dose, calculate classical pharmacological parameters such as the EC50, and enable an in-depth analysis of coordinated dose-dependent treatment effects. The approach was applied to a transcriptional profiling study that evaluated four kinase inhibitors (imatinib, nilotinib, dasatinib and PD0325901) across a six-logarithm dose range, using 12 arrays per compound. The transcript responses proved a powerful means to characterize and compare the compounds: the distribution of EC50 values for the transcriptome was linked to specific targets, dose-dependent effects on cellular processes were identified using automated pathway analysis, and a connection was seen between EC50s in standard cellular assays and transcriptional EC50s. Our approach greatly enriches the information that can be obtained from standard transcriptional profiling technology. Moreover, these methods are automated, robust to non-optimized assays, and could be applied to other sources of quantitative data.

Abatacept does not induce direct gene expression changes in antigen-presenting cells.

BACKGROUND: It has been proposed that ligation of CD80 and CD86 induces reverse signaling into antigen-presenting cells. In this study, we tested the ability of abatacept, a soluble human fusion protein comprising the extracellular domain of cytotoxic T lymphocyte antigen 4 and a fragment of the Fc domain of IgG(1), to activate antigen-presenting cells by measuring changes in global transcriptional responses. METHODS: Affymetrix chips were used to measure gene expression levels using mRNA isolated from immature and mature human dendritic cells and a B cell line following 6 h of treatment with abatacept. RESULTS: In contrast to robust transcriptional responses induced by the control treatment phorbol-12-myristate-13-acetate, abatacept induced minimal gene changes in three different populations of antigen-presenting cells. Furthermore, no gene changes were observed in response to belatacept, a modified version of abatacept that binds with higher avidity to CD80 and CD86. CONCLUSIONS: We conclude that reverse signaling in antigen-presenting cells is unlikely to occur in response to either abatacept or belatacept, thereby supporting the modulation of CD28 signaling on T cells as the main mechanism of action for these therapeutics.

Salivary Gland Stem Cells Age Prematurely in Primary Sjögren's Syndrome.

OBJECTIVE: A major characteristic of the autoimmune disease primary Sjögren's syndrome (SS) is salivary gland (SG) hypofunction. The inability of resident SG stem cells (SGSCs) to maintain homeostasis and saliva production has never been explained and limits our comprehension of mechanisms underlying primary SS. The present study was undertaken to investigate the role of salivary gland stem cells in hyposalivation in primary SS. METHODS: SGSCs were isolated from parotid biopsy samples from controls and patients classified as having primary SS or incomplete primary SS, according to the American College of Rheumatology/European League Against Rheumatism criteria. Self-renewal and differentiation assays were used to determine SGSC regenerative potential, RNA was extracted for sequencing analysis, single telomere length analysis was conducted to determine telomere length, and frozen tissue samples were used for immunohistochemical analysis. RESULTS: SGSCs isolated from primary SS parotid gland biopsy samples were regeneratively inferior to healthy control specimens. We demonstrated that SGSCs from samples from patients with primary SS are not only lower in number and less able to differentiate, but are likely to be senescent, as revealed by telomere length analysis, RNA sequencing, and immunostaining. We further found that SGSCs exposed to primary SS-associated proinflammatory cytokines we induced to proliferate, express senescence-associated genes, and subsequently differentiate into intercalated duct cells. We also localized p16+ senescent cells to the intercalated ducts in primary SS SG tissue, suggesting a block in SGSC differentiation into acinar cells. CONCLUSION: This study represents the first characterization of SGSCs in primary SS, and also the first demonstration of a linkage between an autoimmune disease and a parenchymal premature-aging phenotype. The knowledge garnered in this study indicates that disease-modifying antirheumatic drugs used to treat primary SS are not likely to restore saliva production, and should be supplemented with fresh SGSCs to recover saliva production.

The unfolded protein response is an important regulator of inflammatory genes in endothelial cells.

OBJECTIVE: Oxidized 1-palmitoyl-2-arachidonyl-sn-3-glycero-phosphorylcholine (oxPAPC) accumulates in atherosclerotic lesions and in vitro studies suggest that it mediates chronic inflammatory response in endothelial cells (ECs). The goal of our studies was to identify pathways mediating the induction of inflammatory genes by oxPAPC. METHODS AND RESULTS: Using expression arrays, quantitative polymerase chain reaction (PCR), and immunoblotting we demonstrate that oxPAPC leads to endoplasmic reticulum stress and activation of the unfolded protein response (UPR) in human aortic ECs. Immunohistochemistry analysis of human atherosclerotic lesions indicated that UPR is induced in areas containing oxidized phospholipids. Using the UPR inducing agent tunicamycin and selective siRNA targeting of the ATF4 and XBP1 branches of the UPR, we demonstrate that these transcription factors are essential mediators of IL8, IL6, and MCP1 expression in human aortic ECs required for maximal inflammatory gene expression in the basal state and after oxPAPC treatment. We also identify a novel oxPAPC-induced chemokine, the CXC motif ligand 3 (CXCL3), and show that its expression requires XBP1. CONCLUSIONS: These data suggest that the UPR pathway is a general mediator of vascular inflammation and EC dysfunction in atherosclerosis, and, likely, other inflammatory disorders.

Evaluation and selection of a non-PCR based technology for improved gene expression profiling from clinical formalin-fixed, paraffin-embedded samples.

AIM: Formalin-fixed, paraffin-embedded (FFPE) clinical tissue samples have the potential to provide valuable gene-expression data for the development of cancer biomarkers. However, FFPE RNA is extensively fragmented, presenting a significant challenge for reliably detecting gene expression using traditional qPCR methods. RESULTS: We evaluated three novel methodologies along with the traditional qPCR method for their ability to detect Notch pathway gene expression in colorectal cancer FFPE tissue RNAs. We found that quantitative nuclease protection assay-detected gene expression in high-quality RNAs as sensitively as qPCR, and consistently detected mRNAs in highly fragmented FFPE tissue RNAs. CONCLUSION: Quantitative nuclease protection assay represents an improved methodology for detecting gene expression in FFPE tissue and has the potential to advance the development of cancer biomarkers.

Beneficial and Adverse Effects of an LXR Agonist on Human Lipid and Lipoprotein Metabolism and Circulating Neutrophils.

The development of LXR agonists for the treatment of coronary artery disease has been challenged by undesirable properties in animal models. Here we show the effects of an LXR agonist on lipid and lipoprotein metabolism and neutrophils in human subjects. BMS-852927, a novel LXRß-selective compound, had favorable profiles in animal models with a wide therapeutic index in cynomolgus monkeys and mice. In healthy subjects and hypercholesterolemic patients, reverse cholesterol transport pathways were induced similarly to that in animal models. However, increased plasma and hepatic TG, plasma LDL-C, apoB, apoE, and CETP and decreased circulating neutrophils were also evident. Furthermore, similar increases in LDL-C were observed in normocholesterolemic subjects and statin-treated patients. The primate model markedly underestimated human lipogenic responses and did not predict human neutrophil effects. These studies demonstrate both beneficial and adverse LXR agonist clinical responses and emphasize the importance of further translational research in this area.

Hypoxia influences protein transport and epigenetic repression of CHO cell cultures in shake flasks.

Shake flasks and bench-top bioreactors are widely used for cell culture process development, however, culture performances significantly differ between them. In order to apply the results received from small-scale cultures to production scale, it is important to understand the mechanisms underlying the differences between various culture systems. This study analyzes the expression patterns of Chinese hamster ovary (CHO) cells producing IgG-fusion protein B0 cultured in shake flasks and 5-L bench-top bioreactors by CHO-specific DNA microarrays. The data show that hypoxia was present in shake flask cultures but not in controlled, bench-top bioreactors. Hypoxic conditions appeared to be associated with epigenetic repression resulting in decreased cell culture performance and protein productivity, which is also present during large-scale bioreactor operations due to oxygen gradients. High protein productivity was associated with increased cellular machinery for protein transport and secretion in conjunction with decreased epigenetic repression in bench-top bioreactor cultivation. Metal ions could improve cell growth and protein production under hypoxia and this condition could be mimicked in small-scale bioreactors to facilitate cell culture process scale-up.

Tau overexpression impacts a neuroinflammation gene expression network perturbed in Alzheimer's disease.

Filamentous inclusions of the microtubule-associated protein, tau, define a variety of neurodegenerative diseases known as tauopathies, including Alzheimer's disease (AD). To better understand the role of tau-mediated effects on pathophysiology and global central nervous system function, we extensively characterized gene expression, pathology and behavior of the rTg4510 mouse model, which overexpresses a mutant form of human tau that causes Frontotemporal dementia and parkinsonism linked to chromosome 17 (FTDP-17). We found that the most predominantly altered gene expression pathways in rTg4510 mice were in inflammatory processes. These results closely matched the causal immune function and microglial gene-regulatory network recently identified in AD. We identified additional gene expression changes by laser microdissecting specific regions of the hippocampus, which highlighted alterations in neuronal network activity. Expression of inflammatory genes and markers of neuronal activity changed as a function of age in rTg4510 mice and coincided with behavioral deficits. Inflammatory changes were tau-dependent, as they were reversed by suppression of the tau transgene. Our results suggest that the alterations in microglial phenotypes that appear to contribute to the pathogenesis of Alzheimer's disease may be driven by tau dysfunction, in addition to the direct effects of beta-amyloid.

MicroRNA changes associated with atypical CYP1A1 inducer BMS-764459.

The corticotrophin releasing factor (CRF) receptor I antagonist, BMS-764459 (evaluated as a potential treatment of affective disorders), was orally dosed to female Sprague-Dawley rats once daily for 2 weeks (vehicle control or 175mg/kg/day). To investigate the mechanism of BMS-764459-related liver weight increases, total liver RNA was isolated and evaluated for mRNA gene expression by microarray analysis (assessing the expression of approximately 24,000 genes) from snap-frozen tissue. Subsequently, mRNA and miRNA (microRNA) were also analyzed 5 years later from FFPE (Formalin Fixed Paraffin Embedded) samples via RT-PCR (about 800 miRNA evaluated). Genomic analyses showed that BMS-764459 induces AhR target genes with additional inductions of CYP2B, CYP3A, and Abcc3 consistent with the gene expression pattern of atypical CYP1A1 inducers. Analysis of miRNA expression identified a number of significantly affected miRNAs. To further evaluate their role in atypical CYP1A1 induction, an in silico evaluation of differentially expressed miRNA was performed and their putative mRNA 3'-UTR (untranslated region) binding sequences were evaluated. MiR-680 and miR-29a were identified as potential regulators and biomarkers of atypical CYP1A1 induction by regulating Abcc3, CYP3A and CYP2B as well as a number of AhR targeted genes.