Antonio Di Narzo

Cardiometabolic risk loci share downstream cis- and trans-gene regulation across tissues and diseases

Genetic variation and coronary artery disease Most genetic variants lie outside protein-coding genes, but their effects, especially in human health, are not well understood. Franzén et al. examined gene expression in tissues affected by coronary artery disease (CAD). They found that individuals with loci that have been associated with CAD in genome-wide analyses had different patterns of tissue-specific gene expression than individuals without these genetic variants. Similarly, tissues not associated with CAD did not have CAD-like expression patterns. Thus, tissue-specific data can be used to dissect the genetic effects that predispose individuals to CAD. Science, this issue p. 827 A gene expression survey in patients with coronary artery disease reveals how genetic variation affects the risk of heart failure. Genome-wide association studies (GWAS) have identified hundreds of cardiometabolic disease (CMD) risk loci. However, they contribute little to genetic variance, and most downstream gene-regulatory mechanisms are unknown. We genotyped and RNA-sequenced vascular and metabolic tissues from 600 coronary artery disease patients in the Stockholm-Tartu Atherosclerosis Reverse Networks Engineering Task study (STARNET). Gene expression traits associated with CMD risk single-nucleotide polymorphism (SNPs) identified by GWAS were more extensively found in STARNET than in tissue- and disease-unspecific gene-tissue expression studies, indicating sharing of downstream cis-/trans-gene regulation across tissues and CMDs. In contrast, the regulatory effects of other GWAS risk SNPs were tissue-specific; abdominal fat emerged as an important gene-regulatory site for blood lipids, such as for the low-density lipoprotein cholesterol and coronary artery disease risk gene PCSK9. STARNET provides insights into gene-regulatory mechanisms for CMD risk loci, facilitating their translation into opportunities for diagnosis, therapy, and prevention.

A common haplotype lowers PU.1 expression in myeloid cells and delays onset of Alzheimer's disease

A genome-wide survival analysis of 14,406 Alzheimers disease (AD) cases and 25,849 controls identified eight previously reported AD risk loci and 14 novel loci associated with age at onset. Linkage disequilibrium score regression of 220 cell types implicated the regulation of myeloid gene expression in AD risk. The minor allele of rs1057233 (G), within the previously reported CELF1 AD risk locus, showed association with delayed AD onset and lower expression of SPI1 in monocytes and macrophages. SPI1 encodes PU.1, a transcription factor critical for myeloid cell development and function. AD heritability was enriched within the PU.1 cistrome, implicating a myeloid PU.1 target gene network in AD. Finally, experimentally altered PU.1 levels affected the expression of mouse orthologs of many AD risk genes and the phagocytic activity of mouse microglial cells. Our results suggest that lower SPI1 expression reduces AD risk by regulating myeloid gene expression and cell function.

Prediction of Causal Candidate Genes in Coronary Artery Disease Loci

Objective—Genome-wide association studies have to date identified 159 significant and suggestive loci for coronary artery disease (CAD). We now report comprehensive bioinformatics analyses of sequence variation in these loci to predict candidate causal genes. Approach and Results—All annotated genes in the loci were evaluated with respect to protein-coding single-nucleotide polymorphism and gene expression parameters. The latter included expression quantitative trait loci, tissue specificity, and miRNA binding. High priority candidate genes were further identified based on literature searches and our experimental data. We conclude that the great majority of causal variations affecting CAD risk occur in noncoding regions, with 41% affecting gene expression robustly versus 6% leading to amino acid changes. Many of these genes differed from the traditionally annotated genes, which was usually based on proximity to the lead single-nucleotide polymorphism. Indeed, we obtained evidence that genetic variants at CAD loci affect 98 genes which had not been linked to CAD previously. Conclusions—Our results substantially revise the list of likely candidates for CAD and suggest that genome-wide association studies efforts in other diseases may benefit from similar bioinformatics analyses.

Pericytes impair capillary blood flow and motor function after chronic spinal cord injury

Blood vessels in the central nervous system (CNS) are controlled by neuronal activity. For example, widespread vessel constriction (vessel tone) is induced by brainstem neurons that release the monoamines serotonin and noradrenaline, and local vessel dilation is induced by glutamatergic neuron activity. Here we examined how vessel tone adapts to the loss of neuron-derived monoamines after spinal cord injury (SCI) in rats. We find that, months after the imposition of SCI, the spinal cord below the site of injury is in a chronic state of hypoxia owing to paradoxical excess activity of monoamine receptors (5-HT1) on pericytes, despite the absence of monoamines. This monoamine-receptor activity causes pericytes to locally constrict capillaries, which reduces blood flow to ischemic levels. Receptor activation in the absence of monoamines results from the production of trace amines (such as tryptamine) by pericytes that ectopically express the enzyme aromatic L-amino acid decarboxylase (AADC), which synthesizes trace amines directly from dietary amino acids (such as tryptophan). Inhibition of monoamine receptors or of AADC, or even an increase in inhaled oxygen, produces substantial relief from hypoxia and improves motoneuron and locomotor function after SCI.

Molecularly defined unfolded protein response subclasses have distinct correlations with fatty liver disease in zebrafish

The unfolded protein response (UPR) is a complex network of sensors and target genes that ensure efficient folding of secretory proteins in the endoplasmic reticulum (ER). UPR activation is mediated by three main sensors, which regulate the expression of hundreds of targets. UPR activation can result in outcomes ranging from enhanced cellular function to cell dysfunction and cell death. How this pathway causes such different outcomes is unknown. Fatty liver disease (steatosis) is associated with markers of UPR activation and robust UPR induction can cause steatosis; however, in other cases, UPR activation can protect against this disease. By assessing the magnitude of activation of UPR sensors and target genes in the liver of zebrafish larvae exposed to three commonly used ER stressors (tunicamycin, thapsigargin and Brefeldin A), we have identified distinct combinations of UPR sensors and targets (i.e. subclasses) activated by each stressor. We found that only the UPR subclass characterized by maximal induction of UPR target genes, which we term a stressed-UPR, induced steatosis. Principal component analysis demonstrated a significant positive association between UPR target gene induction and steatosis. The same principal component analysis showed significant correlation with steatosis in samples from patients with fatty liver disease. We demonstrate that an adaptive UPR induced by a short exposure to thapsigargin prior to challenging with tunicamycin reduced both the induction of a stressed UPR and steatosis incidence. We conclude that a stressed UPR causes steatosis and an adaptive UPR prevents it, demonstrating that this pathway plays dichotomous roles in fatty liver disease.

Shared genetic etiology underlying Alzheimer's disease and type 2 diabetes

Epidemiological evidence supports the observation that subjects with type 2 diabetes (T2D) are at higher risk to develop Alzheimers disease (AD). However, whether and how these two conditions are causally linked is unknown. Possible mechanisms include shared genetic risk factors, which were investigated in this study based on recent genome wide association study (GWAS) findings. In order to achieve our goal, we retrieved single nucleotide polymorphisms (SNPs) associated with T2D and AD from large-scale GWAS meta-analysis consortia and tested for overlap among the T2D- and AD-associated SNPs at various p-value thresholds. We then explored the function of the shared T2D/AD GWAS SNPs by leveraging expressional quantitative trait loci, pathways, gene ontology data, and co-expression networks. We found 927 SNPs associated with both AD and T2D with p-value ≤0.01, an overlap significantly larger than random chance (overlapping p-value of 6.93E-28). Among these, 395 of the shared GWAS SNPs have the same risk allele for AD and T2D, suggesting common pathogenic mechanisms underlying the development of both AD and T2D. Genes influenced by shared T2D/AD SNPs with the same risk allele were first identified using a SNP annotation variation (ANNOVAR) software, followed by using Association Protein-Protein Link Evaluator (DAPPLE) software to identify additional proteins that are known to physically interact with the ANNOVAR gene annotations. We found that gene annotations from ANNOVAR and DAPPLE significantly enriched specific KEGG pathways pertaining to immune responses, cell signaling and neuronal plasticity, cellular processes in which abnormalities are known to contribute to both T2D and AD pathogenesis. Thus, our observation suggests that among T2D subjects with common genetic predispositions (e.g., SNPs with consistent risk alleles for T2D and AD), dysregulation of these pathogenic pathways could contribute to the elevated risks for AD in subjects. Interestingly, we found that 532 of the shared T2D/AD GWAS SNPs had divergent risk alleles in the two diseases. For individual shared T2D/AD SNPs with divergent alleles, one of the allelic forms may contribute to one of the diseases (e.g., T2D), but not necessarily to the other (e.g., AD), or vice versa. Collectively, our GWAS studies tentatively support the epidemiological observation of disease concordance between T2D and AD. Moreover, the studies provide the much needed information for the design of future novel therapeutic approaches, for a subpopulation of T2D subjects with genetic disposition to AD, that could benefit T2D and reduce the risk for subsequent development of AD.

A unique gene expression signature associated with serotonin 2C receptor RNA editing in the prefrontal cortex and altered in suicide

Editing of the pre-mRNA for the serotonin receptor 2C (5-HT2CR) by site-specific adenosine deamination (A-to-I pre-mRNA editing) substantially increases the functional plasticity of this key neurotransmitter receptor and is thought to contribute to homeostatic mechanisms in neurons. 5-HT2CR mRNA editing generates up to 24 different receptor isoforms. The extent of editing correlates with 5-HT2CR functional activity: more highly edited isoforms exhibit the least function. Altered 5-HT2CR editing has been reported in postmortem brains of suicide victims. We report a comparative analysis of the connections among 5-HT2CR editing, genome-wide gene expression and DNA methylation in suicide victims, individuals with major depressive disorder and non-psychiatric controls. The results confirm previous findings of an overrepresentation of highly edited mRNA variants (which encode hypoactive 5-HT2CR receptors) in the brains of suicide victims. A large set of genes for which the expression level is associated with editing was detected. This signature set of editing-associated genes is significantly enriched for genes that are involved in synaptic transmission, genes that are preferentially expressed in neurons, and genes whose expression is correlated with the level of DNA methylation. Notably, we report that the link between 5-HT2CR editing and gene expression is disrupted in suicide victims. The results suggest that the postulated homeostatic function of 5-HT2CR editing is dysregulated in individuals who committed suicide.

Decrease of mRNA Editing after Spinal Cord Injury is Caused by Down-regulation of ADAR2 that is Triggered by Inflammatory Response

We recently showed that spinal cord injury (SCI) leads to a decrease in mRNA editing of serotonin receptor 2C (5-HT2CR) contributing to post-SCI spasticity. Here we study post-SCI mRNA editing and global gene expression using massively parallel sequencing. Evidence is presented that the decrease in 5-HT2CR editing is caused by down-regulation of adenosine deaminase ADAR2 and that editing of at least one other ADAR2 target, potassium channel Kv1.1, is decreased after SCI. Bayesian network analysis of genome-wide transcriptome data indicates that down-regulation of ADAR2 (1) is triggered by persistent inflammatory response to SCI that is associated with activation of microglia and (2) results in changes in neuronal gene expression that are likely to contribute both to post-SCI restoration of neuronal excitability and muscle spasms. These findings have broad implications for other diseases of the Central Nervous System and could open new avenues for developing efficacious antispastic treatments.

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.

Airway Epithelial Expression Quantitative Trait Loci Reveal Genes Underlying Asthma and Other Airway Diseases.

Genome-wide association studies (GWASs) have identified loci that are robustly associated with asthma and related phenotypes; however, the molecular mechanisms underlying these associations need to be explored. The most relevant tissues to study the functional consequences of asthma are the airways. We used publically available data to derive expression quantitative trait loci (eQTLs) for human epithelial cells from small and large airways and applied the eQTLs in the interpretation of GWAS results of asthma and related phenotypes. For the small airways (n = 105), we discovered 660 eQTLs at a 10% false discovery rate (FDR), among which 315 eQTLs were not previously reported in a large-scale eQTL study of whole lung tissue. A large fraction of the identified eQTLs is supported by data from Encyclopedia of DNA Elements (ENCODE) showing that the eQTLs reside in regulatory elements (57.5 and 67.6% of cis- and trans-eQTLs, respectively). Published pulmonary GWAS hits were enriched as airway epithelial eQTLs (9.2-fold). Further, genes regulated by asthma GWAS loci in epithelium are significantly enriched in immune response pathways, such as IL-4 signaling (FDR, 5.2 × 10(-4)). The airway epithelial eQTLs described in this study are complementary to previously reported lung eQTLs and represent a powerful resource to link GWAS-associated variants to their regulatory function and thus elucidate the molecular mechanisms underlying asthma and airway-related conditions.