Cis Expression Research Articles

Inside the USCAP Journals

Inside the USCAP Journals

Long Non-Coding RNAs: A Novel Paradigm for Toxicology.

Long non-coding RNAs (lncRNAs) are over 200 nucleotides in length and are transcribed from the mammalian genome in a tissue-specific and developmentally regulated pattern. There is growing recognition that lncRNAs are novel biomarkers and/or key regulators of toxicological responses in humans and animal models. Lacking protein-coding capacity, the numerous types of lncRNAs possess a myriad of transcriptional regulatory functions that include cis and trans gene expression, transcription factor activity, chromatin remodeling, imprinting, and enhancer up-regulation. LncRNAs also influence mRNA processing, post-transcriptional regulation, and protein trafficking. Dysregulation of lncRNAs has been implicated in various human health outcomes such as various cancers, Alzheimer's disease, cardiovascular disease, autoimmune diseases, as well as intermediary metabolism such as glucose, lipid, and bile acid homeostasis. Interestingly, emerging evidence in the literature over the past five years has shown that lncRNA regulation is impacted by exposures to various chemicals such as polycyclic aromatic hydrocarbons, benzene, cadmium, chlorpyrifos-methyl, bisphenol A, phthalates, phenols, and bile acids. Recent technological advancements, including next-generation sequencing technologies and novel computational algorithms, have enabled the profiling and functional characterizations of lncRNAs on a genomic scale. In this review, we summarize the biogenesis and general biological functions of lncRNAs, highlight the important roles of lncRNAs in human diseases and especially during the toxicological responses to various xenobiotics, evaluate current methods for identifying aberrant lncRNA expression and molecular target interactions, and discuss the potential to implement these tools to address fundamental questions in toxicology.

Autoimmune risk variants in ERAP2 are associated with gene-expression levels in thymus.

Genetic polymorphisms in the endoplasmic reticulum aminopeptidase (ERAP)1 and ERAP2 genes have been associated with several autoimmune diseases (AIDs) at a genome-wide significance level. In this study, we performed a cis expression quantitative trait locus (eQTL) screen to investigate whether seven fine-mapped AID single-nucleotide polymorphisms (SNPs) in the ERAP-region influence the gene-expression levels of ERAP1 and ERAP2 in thymus. After quality control, we identified six significant eQTLs. We further assessed the peak eQTL signals, and both genes showed highly significant and independent thymic eQTL signals (P=2.16 × 10-15 and P=8.22 × 10-23, respectively). Interestingly, the peak eQTL signal overlapped with the AID risk loci in ERAP2 (r2>0.94), but were distinct in ERAP1 (r2<0.4). Finally, among the SNPs showing the most significant eQTL associations with ERAP2 (P<3.4 × 10-20), six were located within transcription factor motifs in an enhancer region in thymus. Our study therefore reveals the fine-mapped AID risk variants that act as eQTLs with ERAP2 in thymus, and highlights the potential causal regulatory variants.

Variation in the glucose transporter gene SLC2A2 is associated with glycemic response to metformin.

Metformin is the first-line antidiabetic drug with over 100 million users worldwide, yet its mechanism of action remains unclear1. Here the Metformin Genetics (MetGen) Consortium reports a three-stage genome wide association study (GWAS), consisting of 13,123 participants of different ancestries. The C-allele of rs8192675 in the intron of SLC2A2, which encodes the facilitated glucose transporter GLUT2, was associated with a 0.17% (p=6.6x10-14) greater metformin induced HbA1c reduction in 10,577 participants of European ancestry. rs8192675 is the top cis-eQTL for SLC2A2 in 1,226 human liver samples, suggesting a key role for hepatic GLUT2 in regulation of metformin action. In obese individuals C-allele homozygotes at rs8192675 had a 0.33% (3.6mmol/mol) greater absolute HbA1c reduction than T-allele homozygotes.This is about half the effect seen with the addition of a DPP-4 inhibitor, and equates to a dose difference of 550mg of metformin, suggesting rs8192675 as a potential biomarker for stratified medicine.

Preliminary evidence for association of genetic variants in pri-miR-34b/c and abnormal miR-34c expression with attention deficit and hyperactivity disorder.

Attention deficit and hyperactivity disorder (ADHD) is a prevalent neurodevelopmental disorder characterized by impairment to sustain attention and inability to control impulses and activity level. The etiology of ADHD is complex, with an estimated heritability of 70–80%. Under the hypothesis that alterations in the processing or target binding of microRNAs (miRNAs) may result in functional alterations predisposing to ADHD, we explored whether common polymorphisms potentially affecting miRNA-mediated regulation are involved in this psychiatric disorder. We performed a comprehensive association study focused on 134 miRNAs in 754 ADHD subjects and 766 controls and found association between the miR-34b/c locus and ADHD. Subsequently, we provided preliminary evidence for overexpression of the miR-34c-3p mature form in peripheral blood mononuclear cells of ADHD subjects. Next, we tested the effect on gene expression of single-nucleotide polymorphisms within the ADHD-associated region and found that rs4938923 in the promoter of the pri-miR-34b/c tags cis expression quantitative trait loci for both miR-34b and miR-34c and has an impact on the expression levels of 681 transcripts in trans, including genes previously associated with ADHD. This gene set was enriched for miR-34b/c binding sites, functional categories related to the central nervous system, such as axon guidance or neuron differentiation, and serotonin biosynthesis and signaling canonical pathways. Our results provide preliminary evidence for the contribution to ADHD of a functional variant in the pri-miR-34b/c promoter, possibly through dysregulation of the expression of mature forms of miR-34b and miR-34c and some target genes. These data highlight the importance of abnormal miRNA function as a potential epigenetic mechanism contributing to ADHD.

Genetic risk variants for autoimmune diseases that influence gene expression in thymus.

Genome-wide association studies (GWAS) have boosted our knowledge of genetic risk variants in autoimmune diseases (AIDs). Most risk variants are located within or near genes with immunological functions, and the majority is found to be non-coding, pointing towards a regulatory role. In this study, we performed a cis expression quantitative trait locus (eQTL) screen restricted to 353 AID associated risk variants selected from the GWAS catalog to investigate whether these single nucleotide polymorphisms (SNPs) influence gene expression in thymus. Genotypes were obtained by Immunochip (Ichip) and tested against expression of surrounding genes (±1 Mb) in human thymic tissue (n = 42). We identified eight significant eQTLs located within seven genetic regions (FCRL3, RNASET2, C2orf74, NPIPB8, SIRPG, SYS1 and AJ006998.2) where the expression was associated with AID risk SNPs at a study-wide level of significance (P < 2.7 × 10-5). In NPIPB8 and AJ006998.2, the eQTL signals appeared to be thymus-specific. Furthermore, many AID risk SNPs from GWAS have been subsequently fine-mapped in recent Ichip projects, and fine-mapped AID SNPs overlapped with the thymic eQTLs within RNASET2 and SIRPG Finally, in all the eQTL regions, except C2orf74, SNPs underlying the thymic eQTLs were predicted to interfere with transcription factors important in T cell development. Our study therefore reveals autoimmune risk variants that act as eQTLs in thymus, and suggest that thymic gene regulation may play a functional role at some AID risk loci.

Identification of a functional variant for colorectal cancer risk mapping to chromosome 5q31.1.

Genome-wide association studies (GWASs) have established chromosome 5q31.1 as a risk locus for colorectal cancer (CRC). We previously identified a potentially regulatory single nucleotide polymorphism (SNP) rs17716310 within 5q31.1. Now, we extended our study with another independent Chinese population, functional assays and analyses of TCGA (The Cancer Genome Atlas) data. Significant associations between rs17716310 and CRC risk were found in Present Study including 1075 CRC cases and 1999 controls (additive model: OR = 1.149, 95% CI = 1.027–1.286, P = 0.016), and in Combined Study including 1766 cases and 2708 controls (additive model: OR = 1.145, 95% CI = 1.045–1.254, P = 0.004). Dual luciferase reporter gene assays indicated that the variant C allele obviously increased transcriptional activity. Using TCGA datasets, we indicated rs17716310 as a cis expression quantitative trait locus (eQTL) for the gene SMAD5, whose expression was significantly higher in CRC tissues. These findings suggested that the functional polymorphism rs17716310 A > C might be a genetic modifier for CRC, promoting the expression of SMAD5 that belonged to the transforming growth factor beta (TGF-β) signaling pathway.

Integrative Multi-omic Analysis of Human Platelet eQTLs Reveals Alternative Start Site in Mitofusin 2

Platelets play a central role in ischemic cardiovascular events. Cardiovascular disease (CVD) is a major cause of death worldwide. Numerous genome-wide association studies (GWASs) have identified loci associated with CVD risk. However, our understanding of how these variants contribute to disease is limited. Using data from the platelet RNA and expression 1 (PRAX1) study, we analyzed cis expression quantitative trait loci (eQTLs) in platelets from 154 normal human subjects. We confirmed these results insilico by performing allele-specific expression (ASE) analysis, which demonstrated that the allelic directionality of eQTLs and ASE patterns correlate significantly. Comparison of platelet eQTLs with data from the Genotype-Tissue Expression (GTEx) project revealed that a number of platelet eQTLs are platelet specific and that platelet eQTL peaks localize to the gene body at a higher rate than eQTLs from other tissues. Upon integration with data from previously published GWASs, we found that the trait-associated variant rs1474868 coincides with the eQTL peak for mitofusin 2 (MFN2). Additional experimental and computational analyses revealed that this eQTL is linked to an unannotated alternate MFN2 start site preferentially expressed in platelets. Integration of phenotype data from the PRAX1 study showed that MFN2 expression levels were significantly associated with platelet count. This study links the variant rs1474868 to a platelet-specific regulatory role for MFN2 and demonstrates the utility of integrating multi-omic data with eQTL analysis in disease-relevant tissues for interpreting GWAS results.

Quantitative Trait Locus and Brain Expression of HLA-DPA1 Offers Evidence of Shared Immune Alterations in Psychiatric Disorders.

Genome-wide association studies of schizophrenia encompassing the major histocompatibility locus (MHC) were highly significant following genome-wide correction. This broad region implicates many genes including the MHC complex class II. Within this interval we examined the expression of two MHC II genes (HLA-DPA1 and HLA-DRB1) in brain from individual subjects with schizophrenia (SZ), bipolar disorder (BD), major depressive disorder (MDD), and controls by differential gene expression methods. A third MHC II mRNA, CD74, was studied outside of the MHC II locus, as it interacts within the same immune complex. Exon microarrays were performed in anterior cingulate cortex (ACC) in BD compared to controls, and both HLA-DPA1 and CD74 were decreased in expression in BD. The expression of HLA-DPA1 and CD74 were both reduced in hippocampus, amygdala, and dorsolateral prefrontal cortex regions in SZ and BD compared to controls by specific qPCR assay. We found several novel HLA-DPA1 mRNA variants spanning HLA-DPA1 exons 2-3-4 as suggested by exon microarrays. The intronic rs9277341 SNP was a significant cis expression quantitative trait locus (eQTL) that was associated with the total expression of HLA-DPA1 in five brain regions. A biomarker study of MHC II mRNAs was conducted in SZ, BD, MDD, and control lymphoblastic cell lines (LCL) by qPCR assay of 87 subjects. There was significantly decreased expression of HLA-DPA1 and CD74 in BD, and trends for reductions in SZ in LCLs. The discovery of multiple splicing variants in brain for HLA-DPA1 is important as the HLA-DPA1 gene is highly conserved, there are no reported splicing variants, and the functions in brain are unknown. Future work on the function and localization of MHC Class II proteins in brain will help to understand the role of alterations in neuropsychiatric disorders. The HLA-DPA1 eQTL is located within a large linkage disequilibrium block that has an irrefutable association with schizophrenia. Future tests in a larger cohort are needed to determine the significance of this eQTL association with schizophrenia. Our findings support the long-held hypothesis that alterations in immune function are associated with the pathophysiology of psychiatric disorders.

DRD2 Schizophrenia-Risk Allele Is Associated With Impaired Striatal Functioning in Unaffected Siblings of Schizophrenia Patients.

A recent Genome-Wide Association Study showed that the rs2514218 single nucleotide polymorphism (SNP) in close proximity to dopamine receptor D2 is strongly associated with schizophrenia. Further, an in silico experiment showed that rs2514218 has a cis expression quantitative trait locus effect in the basal ganglia. To date, however, the functional consequence of this SNP is unknown. Here, we used functional Magnetic resonance imaging to investigate the impact of this risk allele on striatal activation during proactive and reactive response inhibition in 45 unaffected siblings of schizophrenia patients. We included siblings to circumvent the illness specific confounds affecting striatal functioning independent from gene effects. Behavioral analyses revealed no differences between the carriers (n= 21) and noncarriers (n= 24). Risk allele carriers showed a diminished striatal response to increasing proactive inhibitory control demands, whereas overall level of striatal activation in carriers was elevated compared to noncarriers. Finally, risk allele carriers showed a blunted striatal response during successful reactive inhibition compared to the noncarriers. These data are consistent with earlier reports showing similar deficits in schizophrenia patients, and point to a failure to flexibly engage the striatum in response to contextual cues. This is the first study to demonstrate an association between impaired striatal functioning and the rs2514218 polymorphism. We take our findings to indicate that striatal functioning is impaired in carriers of the DRD2 risk allele, likely due to dopamine dysregulation at the DRD2 location.

Abstract 12905: Gene Expression and Expression Quantifying Trait Loci Analysis of Left Atrium of Patients With Post-operative Atrial Fibrillation

Introduction: The majority of information on the genetic background of atrial fibrillation (AF) results from genomic DNA variant analysis without consideration of tissue expression. Hypothesis: Analysis of tissue-specific gene expression in left atrium (LA) can further understanding of the molecular mechanism of identified AF risk variants, and identify novel genes and gene variants associated with AF. Methods: We isolated mRNA from samples of the LA free wall taken during mitral valve surgery in 62 Caucasian individuals. Gene expression in the LA was compared between patients who did and did not have post-operative AF (poAF) using high-throughput RNA expression. Using genotypes of 1.4 million single nucleotide polymorphisms (SNP) we performed cis expression quantifying trait loci (eQTL) analysis, correlating gene expression of each gene with the genotypes of adjacent (&lt;1Mbp) SNPs. Results: We identified 23 differentially expressed genes in the LA of patients with poAF, including three potassium channel genes (KCNA7, KCNH8 and KCNK17). The largest expression difference was in LOC645323, a long non-coding RNA. The expression of PITX2, ZFHX3 and KCNN3, previously shown to be associated with AF, did not differ between patients with and without poAF. We identified 12,476 cis eQTL relationships in the LA, several of those included genetic regions and genes previously associated with AF. We confirmed an eQTL relationship between rs3744029 genotype and the expression of MYOZ1. Furthermore we describe a novel eQTL relationship between rs6795970 genotype and the expression of the SCN10A gene. Conclusions: We have analysed the human LA expression via high-throughput RNA sequencing, and identified novel genes and gene variants likely involved in the molecular pathophysiology of AF.

Genome-wide significant association with seven novel multiple sclerosis risk loci

ObjectiveA recent large-scale study in multiple sclerosis (MS) using the ImmunoChip platform reported on 11 loci that showed suggestive genetic association with MS. Additional data in sufficiently sized and independent...

Identification of quantitative trait loci influencing inflammation-mediated alveolar bone loss: insights into polygenic inheritance of host-biofilm disequilibria in periodontitis.

The relative contribution of genetic and environmental factors to the onset and progression of periodontitis is inconclusive. Despite the high prevalence, phenotypic heterogeneity and significant local and systemic implications of this disease, early detection and individualized therapy are problematic. Using a murine model of periodontitis in a panel of 17 recombinant inbred mice, the current study addressed the heritability of, and oral dysbiosis associated with, inflammation-mediated alveolar bone loss (iABL), the hallmark of periodontitis. Quantitative trait locus (QTL) genomics and quantitative PCR for over 99% of known murine oral microbiota were used. It was found that iABL is a polygenic trait with 32.7% heritability. One suggestive QTL, nicknamed inflammation-mediated alveolar bone loss locus (iABLL), was identified on chromosome 2. Eleven genes involved in innate immune responses and bone metabolism, particularly related to macrophage and osteoblast function, namely Etl4, Pdss1, Cobll1, 9330158F14Rik, Xirp2, Stk39, Mettl5, Metapl1, Itga6, Pdk1 and Sp3, were found in the iABLL using cis expression QTL and nonsynonymous single nucleotide polymorphism analyses. Specific oral microbiome shifts in saliva and tongue mucosa are associated with disease in this model. Our results indicate that complex host-biofilm interactions generate pathogenic states that extend beyond subgingival biofilms and periodontal tissues. Although no temporal relationship between the onset of iABL and microbiome changes were established, our findings suggest that host factors may be responsible for pathogenic shifts in subgingival biofilms when persistent and undisturbed.

A common variant in the CLDN7/ELP5 locus predicts adiponectin change with lifestyle intervention and improved fitness in obese individuals with diabetes.

Overweight/obese individuals with Type 2 diabetes have low adiponectin levels, which may improve with lifestyle changes. We investigated whether genetic variants associated with adiponectin levels in genome-wide association studies (GWAS) would also be related with adiponectin changes in response to an intensive lifestyle intervention (ILI), potentially through mechanisms altering the adipose microenvironment via weight loss and/or improved cardiorespiratory fitness. Look AHEAD was a randomized trial comparing the cardiovascular benefits of ILI-induced weight loss and physical activity compared with diabetes support and education among overweight/obese individuals with Type 2 diabetes. In a subsample of Look AHEAD with adiponectin data and genetic consent (n=1,351), we evaluated the effects of 24 genetic variants, demonstrated by GWAS to be cross-sectionally associated with adiponectin, on adiponectin change 1-yr postintervention. We explored via mediational analyses whether any differential effects by treatment arm were occurring through weight loss and/or improved fitness. A variant, rs222857, in the CLDN7 locus, potentially associated with epithelial barrier integrity and tight junction physiology, and a putative cis expression quantitative trail locus for elongator acetyltransferase complex subunit 5 (ELP5), predicted adiponectin increases within ILI (log-adiponectin in overall sample per copy: β±SE=0.05±0.02, P=0.008; in non-Hispanic whites: 0.06±0.02, P=0.009). The favorable effects of rs222857 (minor allele frequency 45.5%) appeared to be mediated by mechanisms associated with improved fitness, and not weight loss. This is the first study to identify a genetic variant that modifies adiponectin response to lifestyle intervention in overweight/obese diabetic individuals.

Loss of fibroblast growth factor 21 action induces insulin resistance, pancreatic islet hyperplasia and dysfunction in mice.

Fibroblast growth factor (FGF) 21 is an endocrine factor that normalizes glucose homeostasis and reduces insulin resistance in diabetes. Although the pancreas is an FGF21 target organ, its role in pancreatic islets remains obscure. This study aimed to elucidate the physiological role of FGF21 in pancreatic islets using FGF21-knockout (FGF21-KO) mice. Twenty-four-week-old male global FGF21-KO mice were used in this study. Glucose and insulin tolerance were assessed. Expression of genes and proteins related to islet function and underlying mechanisms were also examined. Islet morphology and insulin-secreting capacity were further evaluated. FGF21-KO mice exhibited insulin resistance while being normoglycemic, associated with increases in beta-cell proliferation and insulin synthesis, acting as compensatory responses. This phenotype probably results from enhanced growth hormone (GH) sensitivity in FGF21-KO mouse islets. In addition, ex vivo FGF21 treatment in normal C57BL/6J mouse islets reduced GH signaling, probably via upregulation of peroxisome proliferator-activated receptor gamma (PPARγ) and cytokine-inducible SH-2 containing (CIS) protein, whereas KO mouse islets displayed reduced PPARγ and CIS expression. FGF21 treatment also reversed GH-induced insulin expression, beta-cell proliferation and GH-impaired glucose-stimulated insulin secretion (GSIS) in islets. Furthermore, distorted islet morphology and impaired GSIS were observed in KO mice, suggestive of islet dysfunction, whereas the enhanced insulin expression and impaired GSIS in FGF21-KO mouse islets could be reversed by blockade of GH signaling. Our data indicate that FGF21 is important in the regulation of beta-cell proliferation and insulin synthesis, probably via modulation of GH signaling. These findings provide evidence that FGF21 is an obligatory metabolic regulator in pancreatic islets and shed new light onto the role of endogenous FGF21 in the pathogenesis of insulin resistance and islet dysfunction.

Relationship of DNA methylation and gene expression in idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is an untreatable and often fatal lung disease that is increasing in prevalence and is caused by complex interactions between genetic and environmental factors. Epigenetic mechanisms control gene expression and are likely to regulate the IPF transcriptome. To identify methylation marks that modify gene expression in IPF lung. We assessed DNA methylation (comprehensive high-throughput arrays for relative methylation arrays [CHARM]) and gene expression (Agilent gene expression arrays) in 94 patients with IPF and 67 control subjects, and performed integrative genomic analyses to define methylation-gene expression relationships in IPF lung. We validated methylation changes by a targeted analysis (Epityper), and performed functional validation of one of the genes identified by our analysis. We identified 2,130 differentially methylated regions (DMRs; <5% false discovery rate), of which 738 are associated with significant changes in gene expression and enriched for expected inverse relationship between methylation and expression (P < 2.2 × 10(-16)). We validated 13/15 DMRs by targeted analysis of methylation. Methylation-expression quantitative trait loci (methyl-eQTL) identified methylation marks that control cis and trans gene expression, with an enrichment for cis relationships (P < 2.2 × 10(-16)). We found five trans methyl-eQTLs where a methylation change at a single DMR is associated with transcriptional changes in a substantial number of genes; four of these DMRs are near transcription factors (castor zinc finger 1 [CASZ1], FOXC1, MXD4, and ZDHHC4). We studied the in vitro effects of change in CASZ1 expression and validated its role in regulation of target genes in the methyl-eQTL. These results suggest that DNA methylation may be involved in the pathogenesis of IPF.

Gene Age Predicts the Strength of Purifying Selection Acting on Gene Expression Variation in Humans

Gene expression levels can be subject to selection. We hypothesized that the age of gene origin is associated with expression constraints, given that it affects the level of gene integration into the functional cellular environment. By studying the genetic variation affecting gene expression levels (cis expression quantitative trait loci [cis-eQTLs]) and protein levels (cis protein QTLs [cis-pQTLs]), we determined that young, primate-specific genes are enriched in cis-eQTLs and cis-pQTLs. Compared to cis-eQTLs of old genes originating before the zebrafish divergence, cis-eQTLs of young genes have a higher effect size, are located closer to the transcription start site, are more significant, and tend to influence genes in multiple tissues and populations. These results suggest that the expression constraint of each gene increases throughout its lifespan. We also detected a positive correlation between expression constraints (approximated by cis-eQTL properties) and coding constraints (approximated by Ka/Ks) and observed that this correlation might be driven by gene age. To uncover factors associated with the increase in gene-age-related expression constraints, we demonstrated that gene connectivity, gene involvement in complex regulatory networks, gene haploinsufficiency, and the strength of posttranscriptional regulation increase with gene age. We also observed an increase in heritability of gene expression levels with age, implying a reduction of the environmental component. In summary, we show that gene age shapes key gene properties during evolution and is therefore an important component of genome function.

Linking polymorphic p53 response elements with gene expression in airway epithelial cells of smokers and cancer risk.

Chronic cigarette smoking exposes airway epithelial cells to thousands of carcinogens, oxidants and DNA-damaging agents, creating a field of molecular injury in the airway and altering gene expression. Studies of cytologically normal bronchial epithelial cells from smokers have identified transcription-based biomarkers that may prove useful in early diagnosis of lung cancer, including a number of p53-regulated genes. The ability of p53 to regulate transcription is critical for tumor suppression, and this suggests that single-nucleotide polymorphisms (SNPs) in functional p53 binding sites (p53 response elements, or p53REs) that affect gene expression could influence susceptibility to cancer. To connect p53RE SNP genotype with gene expression and cancer risk, we identified a set of 204 SNPs in putative p53REs, and performed cis expression quantitative trait loci (eQTL) analysis, assessing associations between SNP genotypes and mRNA levels of adjacent genes in bronchial epithelial cells obtained from 44 cigarette smokers. To further test and validate these genotype-expression associations, we searched published eQTL studies from independent populations and determined that 53% (39/74) of the bronchial epithelial eQTLs were observed in at least one of other studies. SNPs in p53REs were also evaluated for effects on p53-DNA binding using a quantitative in vitro protein-DNA binding assay. Last, based on linkage disequilibrium, we found 6 p53RE SNPs associated with gene expression were identified as cancer risk SNPs by either genome-wide association studies or candidate gene studies. We provide an approach for identifying and evaluating potentially functional SNPs that may modulate the airway gene expression response to smoking and may influence susceptibility to cancers.

Genetic variation within the Chrna7 gene modulates nicotine reward‐like phenotypes in mice

Mortality from tobacco smoking remains the leading cause of preventable death in the world, yet current cessation therapies are only modestly successful, suggesting new molecular targets are needed. Genetic analysis of gene expression and behavior identified Chrna7 as potentially modulating nicotine place conditioning in the BXD panel of inbred mice. We used gene targeting and pharmacological tools to confirm the role of Chrna7 in nicotine conditioned place preference (CPP). To identify molecular events downstream of Chrna7 that may modulate nicotine preference, we performed microarray analysis of α7 knock-out (KO) and wild-type (WT) nucleus accumbens (NAc) tissue, followed by confirmation with quantitative polymerase chain reaction (PCR) and immunoblotting. In the BXD panel, we found a putative cis expression quantitative trait loci (eQTL) for Chrna7 in NAc that correlated inversely to nicotine CPP. We observed that gain-of-function α7 mice did not display nicotine preference at any dose tested, whereas conversely, α7 KO mice demonstrated nicotine place preference at a dose below that routinely required to produce preference. In B6 mice, the α7 nicotinic acetylcholine receptor (nAChR)-selective agonist, PHA-543613, dose-dependently blocked nicotine CPP, which was restored using the α7 nAChR-selective antagonist, methyllycaconitine citrate (MLA). Our genomic studies implicated a messenger RNA (mRNA) co-expression network regulated by Chrna7 in NAc. Mice lacking Chrna7 demonstrate increased insulin signaling in the NAc, which may modulate nicotine place preference. Our studies provide novel targets for future work on development of more effective therapeutic approaches to counteract the rewarding properties of nicotine for smoking cessation.

Fraction of exhaled nitric oxide values in childhood are associated with 17q11.2-q12 and 17q12-q21 variants

The fraction of exhaled nitric oxide (Feno) value is a biomarker of eosinophilic airway inflammation and is associated with childhood asthma. Identification of common genetic variants associated with childhood Feno values might help to define biological mechanisms related to specific asthma phenotypes. We sought to identify the genetic variants associated with childhood Feno values and their relation with asthma. Feno values were measured in children age 5 to 15 years. In 14 genome-wide association studies (N = 8,858), we examined the associations of approximately 2.5 million single nucleotide polymorphisms (SNPs) with Feno values. Subsequently, we assessed whether significant SNPs were expression quantitative trait loci in genome-wide expression data sets of lymphoblastoid cell lines (n = 1,830) and were related to asthma in a previously published genome-wide association data set (cases, n = 10,365; control subjects: n = 16,110). We identified 3 SNPs associated with Feno values: rs3751972 in LYR motif containing 9 (LYRM9; P = 1.97 × 10(-10)) and rs944722 in inducible nitric oxide synthase 2 (NOS2; P = 1.28 × 10(-9)), both of which are located at 17q11.2-q12, and rs8069176 near gasdermin B (GSDMB; P = 1.88 × 10(-8)) at 17q12-q21. We found a cis expression quantitative trait locus for the transcript soluble galactoside-binding lectin 9 (LGALS9) that is in linkage disequilibrium with rs944722. rs8069176 was associated with GSDMB and ORM1-like 3 (ORMDL3) expression. rs8069176 at 17q12-q21, but not rs3751972 and rs944722 at 17q11.2-q12, were associated with physician-diagnosed asthma. This study identified 3 variants associated with Feno values, explaining 0.95% of the variance. Identification of functional SNPs and haplotypes in these regions might provide novel insight into the regulation of Feno values. This study highlights that both shared and distinct genetic factors affect Feno values and childhood asthma.