Laura B. Ferguson
Lingling Zhang
Shi Wang
Courtney Bridges1
R. Adron Harris
Igor Ponomarev- 1Waggoner Center for Alcohol & Addiction Research, The University of Texas at Austin, Austin, TX, United States
- 2MOE Key Laboratory of Marine Genetics and Breeding, Ocean University of China, Qingdao, China
- 3Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology, Qingdao, China
Peroxisome proliferator activated receptors (PPARs) are nuclear hormone receptors that act as transcription factors in response to endogenous lipid messengers. The fibrates and thiazolidinediones are synthetic PPAR agonists used clinically to treat dyslipidemia and Type 2 Diabetes Mellitus, respectively, but also improve symptoms of several other diseases. Transposable elements (TEs), repetitive sequences in mammalian genomes, are implicated in many of the same conditions for which PPAR agonists are therapeutic, including neurodegeneration, schizophrenia, and drug addiction. We tested the hypothesis that there is a link between actions of PPAR agonists and TE expression. We developed an innovative application of microarray data by mapping Illumina mouse WG-6 microarray probes to areas of the mouse genome that contain TEs. Using this information, we assessed the effects of systemic administration of three PPAR agonists with different PPAR subtype selectivity: fenofibrate, tesaglitazar, and bezafibrate, on TE probe expression in mouse brain [prefrontal cortex (PFC) and amygdala] and liver. We found that fenofibrate, and bezafibrate to a lesser extent, up-regulated probes mapped to retrotransposons: Short-Interspersed Elements (SINEs) and Long-Interspersed Elements (LINEs), in the PFC. Conversely, all PPAR agonists down-regulated LINEs and tesaglitazar and bezafibrate also down-regulated SINEs in liver. We built gene coexpression networks that partitioned the diverse transcriptional response to PPAR agonists into groups of probes with highly correlated expression patterns (modules). Most of the differentially expressed retrotransposons were within the same module, suggesting coordinated regulation of their expression, possibly by PPAR signaling. One TE module was conserved across tissues and was enriched with genes whose products participate in epigenetic regulation, suggesting that PPAR agonists affect TE expression via epigenetic mechanisms. Other enriched functional categories included phenotypes related to embryonic development and learning and memory, suggesting functional links between these biological processes and TE expression. In summary, these findings suggest mechanistic relationships between retrotransposons and PPAR agonists and provide a basis for future exploration of their functional roles in brain and liver.
Introduction
Peroxisome proliferator activated receptors (PPARs), members of the nuclear hormone receptor family, are ligand-activated transcription factors. There are three PPAR isotypes: PPARα, PPARγ and PPARβ (sometimes called PPARδ) with distinct expression patterns, tissue distribution, and ligand specificity (Grygiel-Gorniak, 2014). PPARα is most highly abundant in metabolically active tissues including liver, heart, intestine, and kidney, PPARγ is most abundant in adipose tissue, and PPARβ/δ is ubiquitously expressed, often at relatively higher levels than PPARα or PPARγ (Kliewer et al., 1994; Braissant et al., 1996). All PPAR isotypes are expressed in brain, including the prefrontal cortex (PFC) and amygdala, and appear to be most highly expressed in neurons (Moreno et al., 2004; Warden et al., 2016). Their main function is to regulate fatty acid metabolism but also have roles in inflammation, glucose metabolism, adipogenesis and myelination [for review, see (Tyagi et al., 2011)]. Fibrates are PPARα agonists used to treat dyslipidemia, thiazolidinediones (TZDs) are PPARγ agonists used to treat type 2 diabetes, and the glitazars are dual PPARα and PPARγ agonists used to treat dyslipidemia and type 2 diabetes comorbidities. Although PPARβ /PPARδ agonists are not used clinically, they have been used as performance enhancers by athletes and have been prohibited by the World Anti-Doping Agency.
Transposable elements (TEs) are repeating genomic elements that have the ability to translocate across the genome. First discovered by Barbara McClintock in maize in the 1940s, TEs have been found across different taxonomic groups. They can have both beneficial and deleterious effects on the host. They can drive diversity and evolution through their role as controlling elements that can lead to novel regulation of genes, but they can also drive disease through their role as insertional mutagens [for review, see (Reilly et al., 2013; Chuong et al., 2017)]. Many TEs are conserved across species including fruit flies, rodents and humans, where they are estimated to make up almost half of the genome (Adams et al., 2000; Nellaker et al., 2012). TEs can be divided into two major classes: DNA transposons that translocate via a cut-and-paste method (but are not mobile in mammals), and retrotransposons that translocate via a copy-and-paste method using an RNA intermediate. Retrotransposons can be further subdivided into categories based on whether or not they contain Long Terminal Repeats (LTRs, e.g., endogenous retroviruses - ERVs). Retrotransposons lacking LTRs include Long-Interspersed Elements (LINEs), Short-Interspersed Elements (SINEs), Alu and SVA elements. For more information on different classes of TEs the interested reader is referred to reviews (Muotri et al., 2007; Reilly et al., 2013; Chuong et al., 2017).
Although representing a large portion of the genome in different species, most TE sequences are non-functional and cannot move in the genome. The host uses sophisticated epigenetic mechanisms to create highly condensed heterochromatin to regulate TE expression (Slotkin and Martienssen, 2007; Maze et al., 2011; Ponomarev et al., 2012; Xie et al., 2013; Du et al., 2016). Inhibition of retrotransposon expression occurs through extensive DNA methylation at retrotransposon promoter sequences [reviewed in (Yoder et al., 1997)]. However, retrotransposons can be expressed when the epigenetic silencing is released (Slotkin and Martienssen, 2007). For example, one study linked hypomethylated genomic areas with increased expression of LTRs in human alcoholic frontal cortex compared to controls (Ponomarev et al., 2012), another study linked hypomethylation with increased activation of ERVs in mouse tumors (Howard et al., 2008), and another found that cocaine-mediated reductions of a heterochromatic marker (H3K9me3) were associated with increased LINE-1 expression in mouse nucleus accumbens (Maze et al., 2011).
Transposable elements have been implicated in many pathological conditions, including approximately 100 single-gene diseases (Hancks and Kazazian, 2012), such as hemophilia A (Sukarova et al., 2001; Ganguly et al., 2003; Green et al., 2008), and cystic fibrosis (Chen et al., 2008). Retrotransposons (ERVs and LINE1 elements) have also been linked to more complex diseases including diabetes (Pascual et al., 2001; Dickerson et al., 2008), cancers [for review, see (Burns, 2017)], and several psychiatric disorders (Guffanti et al., 2014), such as alcohol and cocaine addiction (Maze et al., 2011; Ponomarev et al., 2012), autism (Balestrieri et al., 2012), Attention Deficit Hyperactivity Disorder (Balestrieri et al., 2014), depression (Weis et al., 2007), schizophrenia (Yolken et al., 2000; Karlsson et al., 2001, 2004; Frank et al., 2005; Weis et al., 2007; Dickerson et al., 2008; Perron et al., 2008, 2012a,b; Yao et al., 2008; Huang et al., 2011; Lin et al., 2011; Bundo et al., 2014; Slokar and Hasler, 2015), bipolar disorder (Yolken et al., 2000; Weis et al., 2007; Frank et al., 2005; Perron et al., 2012b), post-traumatic stress disorder (Ponomarev et al., 2010; Rusiecki et al., 2012), Rett syndrome (Muotri et al., 2010), and neurodegeneration (Cartault et al., 2012; Li et al., 2012). Besides their clinical uses for dyslipidemia and diabetes, PPAR agonists have shown therapeutic potential (mostly preclinical evidence) for several brain diseases including neurodegeneration (Bordet et al., 2006; Barbiero et al., 2014; Fidaleo et al., 2014; Avagliano et al., 2016; Dickey et al., 2016; Makela et al., 2016; Zhou et al., 2016; Wang Y. et al., 2017), mood disorders (Kemp et al., 2014; Ji et al., 2015; Scheggi et al., 2016; Colle et al., 2017; Jiang et al., 2017; Li et al., 2018; Ni et al., 2018), schizophrenia (Rolland et al., 2012), and several types of substance use disorders (Melis et al., 2010; Stopponi et al., 2011; Panlilio et al., 2012, 2013; Le Foll et al., 2013; Stopponi et al., 2013; Ferguson et al., 2014; Karahanian et al., 2014; Blednov et al., 2015; de Guglielmo et al., 2015; Haile and Kosten, 2017; Rivera-Meza et al., 2017). Therefore, PPAR agonists and TEs are clinically related, i.e., are associated with similar diseases. The rationale for this study is based on this clinical relatedness; we reasoned that because PPAR agonists and TEs are associated with many of the same diseases that their mechanisms of action could be interacting. Here, we tested the hypothesis that activation of PPARs regulates transcriptional activity of TEs. We mapped Illumina microarray probes to genomic locations that correspond to different classes of TEs and used this information to study the effects of PPAR agonists on expression of TEs in the amygdala, PFC, and liver of male C57Bl/6J mice.
Materials and Methods
Animals and Drug Administration
We previously measured the effects of peripherally administered PPAR agonists on brain and liver gene expression (Ferguson et al., 2014). Briefly, we administered (p.o.) saline or one of three PPAR agonists: 150 mg/kg of the PPARα agonist fenofibrate (feno), 75 mg/kg of the pan-PPAR agonist bezafibrate (beza) or 1.5 mg/kg of the dual PPARα-PPARγ agonist tesaglitazar (tesa) to male C57Bl/6J mice (8–12 weeks old, n = 10 per group) once per day for 8 days. See (Ferguson et al., 2014) for details. All experiments were approved by The University of Texas at Austin Institute for Animal Care and Use Committee and were conducted in accordance with NIH guidelines with regard to the use of animals in research.
Gene Expression
Twenty-four hours after the last administration of PPAR agonist, amygdala, PFC, and liver were harvested. We extracted, amplified, and biotin-labeled total RNA and sent aliquots of labeled cRNA to the Yale Center for Genome Analysis, where they were hybridized to Illumina® MouseWG-6 v2 Expression BeadChips. See (Ferguson et al., 2014) for details. Microarray data were deposited in the NCBI Gene Expression Omnibus (GEO) database under accession number GSE67796.
Transposable Element (TE) Annotation
To identify microarray probes targeting potential TE regions, we first obtained probe sequences of all Illumina probes from the probe annotation file MouseWG-6_V2_0_R3_11278593_A1 and mapped them to mouse genome (NCBI37/mm9) downloaded from the University of California, Santa Cruz (UCSC) Genome Bioinformatics website2 using Short Oligonucleotide Analysis Package (Li et al., 2009). Then we compared the probe genomic coordinates with the coordinates of TEs reported in the Repeat Masker table downloaded from UCSC Genome Bioinformatics website. A probe was considered to be targeting a TE region only if it falls completely within an annotated TE region. We constructed four TE data sets covering the major TE classes: DNA transposons, long-term repeat (LTR), long interspersed nuclear element (LINE), and short interspersed nuclear elements (SINE). Supplementary Table S1 contains the list of Illumina probes mapped to each TE class.
Statistical and Gene Network Analyses
Data from each tissue were analyzed separately. Variance stabilization transformation and quantile normalization were used to pre-process the data using the Bioconductor Lumi package (Du et al., 2008). Groups being compared were normalized together, i.e., saline and bezafibrate, saline and tesaglitazar, and saline and fenofibrate. Outlier values for each gene within a group were removed using Grubbs’ test (p < 0.05). To detect genes differentially expressed between treated and control mice, we used the Bioconductor Limma package to fit a linear model for each gene (Ritchie et al., 2015). This analysis generated a t-value for each probe, with a positive t-value indicating up- and a negative t-value indicating a down-regulation after drug treatment.
We conducted weighted gene co-expression network analysis (WGCNA) as previously described (Ferguson et al., 2014). Briefly, all reliably detected genes were included in the network construction for the gene network analysis. Data from all treatments (fenofibrate, tesaglitazar, bezafibrate and saline) were used to detect co-expression patterns. Signed networks were constructed using R and custom functions available at https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/Rpackages/WGCNA/. The power was set to β = 11, β = 14 and β = 7 for the amygdala, PFC and liver network, respectively, and minimum module size of 100, 150, and 100 genes for amygdala, PFC and liver, respectively. A dendrogram cut height of 0.99 was used for amygdala and PFC, and 0.995 was used for liver. Similar modules were merged to reduce the number of modules used for analyses. Gene modules corresponding to branches of the dendrogram were labeled in unique colors. Genes whose profile failed to cluster in the network were labeled in gray.
Bioinformatics Analysis
Hypergeometric test was performed to assess TE enrichment for each set of differentially expressed genes or gene network module and Benjamini-Hochberg’s FDR was calculated to correct for multiple tests (except where noted otherwise). We plotted the t-value distribution of each TE class and used a one-sample, two-tailed t-test to determine if the average t-values for those TE classes were different from 0 (Bonferroni-corrected p-value < 0.05 was the threshold we used for statistical significance).
For functional annotation of the TE-enriched modules, we used the hypergeometric test to assess enrichment of cell types [neurons, oligodendrocytes, microglia, astrocytes, hepatocytes, Kupffer cells, and hepatic stellate cells (HSCs) (Sugino et al., 2006; Takahara et al., 2006; Cahoy et al., 2008; Oldham et al., 2008)] and the PPAR agonist-regulated genesets (i.e., genes differentially expressed at a nominal p-value < 0.05 after PPAR agonist treatment) as previously described (Ferguson et al., 2014). We curated the list of ethanol consumption-related genes from the literature on mutant mouse data (gene knockdown or overexpression). More information can be found in these reviews (Crabbe et al., 2006; Mayfield et al., 2016). Because cells use epigenetic mechanisms to control TEs (See Introduction), we looked for enrichment of 171 genes whose encoded proteins are involved in epigenetic chromatin remodeling and modification in the TE-enriched modules (attained from the SABiosciences website; Supplementary Table S2).
Comparison of Tissue Networks
We built and analyzed co-expression networks for each tissue individually and compared the co-expression patterns of each tissue’s network. The hypergeometric distribution was used to assess the significance of internetwork module overlap. We used Cytoscape 3.2.1 (Shannon et al., 2003) to visualize the network comparisons (Killcoyne et al., 2009). Only conservative modules with overlap p < 10-3 are presented. The 13,623 probes that were common between amygdala, PFC, and liver networks were used in the network comparison. To determine the enriched pathways, phenotypes and ontologies of the genes in the conserved TE module we used WebGestalt3 (Wang J. et al., 2017). We used the criteria that the enriched category must include ≥ 3 genes and pass significance threshold of FDR-adjusted hypergeometric p < 0.05. Phenotype data are derived from literature, knockout data, and other sources curated by Mouse Genome Informatics.
qRT-PCR Validation
We validated Illumina probe ILMN_1244099 (which corresponds to ORF1 of a LINE1 element). We designed the probe in such a way as to detect the LINE1 elements, specifically. Tissue and total RNA isolation were described previously (Ferguson et al., 2014). PFC cDNA was synthesized with the iScript Select cDNA Synthesis Kit (BIORAD, priming strategy: oligo(dT) primers) for 37 samples (due to limited RNA quantity, 3 PFC samples could not be included). Sample information can be found on GEO accession GSE67796. Primer Express was used to design a Taqman assay to detect the LINE1 TE that contained the ILMN_1244099 probe sequence. Gapdh was used as an endogenous control. qBase Plus software (Biogazelle, Zwijnaarde, Belgium) was used for data analysis. A two-sided Mann Whitney test was used to assess the significance of the between group differences.
Taqman Assay [45–108] amp. length of 64
Forward Primer: CCAGAAGAGCCTGGACAGATG
Reverse primer: AGTAGCCTGGGCTGGCATT
Probe: TATACAGACACTAAGAGAACAC
Results
TE Mapping
We previously studied the organization of brain transcriptome in human alcoholics and discovered that TEs, such as LTRs and SINEs, are co-expressed and regulated in alcoholics (Ponomarev et al., 2012). Here, we performed TE enrichment analysis to determine if a TE regulation also exists in mouse, and if PPAR agonist treatment has any effect on TE regulation. Out of 45,281 probes on the microarray, a total of 2,071 probes were found to target TEs, including 104 DNA transposons, 504 LINEs, 576 LTRs and 887 SINEs (Supplementary Tables S1, S4). Overrepresentation analysis indicated that the probes mapped to TEs contained a disproportionate number of unannotated genes, supporting our mapping method (Supplementary Figure S1).
PPAR Agonists Increase TE Expression in the PFC and Decrease TE Expression in the Liver
We previously characterized brain and liver gene expression profiles following an 8-day, treatment with PPAR agonists (Ferguson et al., 2014). Here, we used the aforementioned mapping information to study TE expression in the PPAR agonist-regulated genesets (we define PPAR agonist-regulated genes as those that were differentially expressed between PPAR agonist and saline treatment at p < 0.05). We used the hypergeometric test to determine if more TEs were in the PPAR agonist-regulated genesets than chance level and used Bonferroni corrected p < 0.05 as a statistical threshold. We found that SINEs and LINEs were enriched in the genesets that are up-regulated by fenofibrate in the PFC and SINEs are enriched in the genesets up-regulated by bezafibrate in the PFC (Table 1). Conversely to the overall upregulation of TEs observed in PFC, LINEs were enriched in the genesets down-regulated by bezafibrate and tesaglitazar in liver (Table 1). TEs were not enriched in the genesets regulated by PPAR agonists in amygdala (Supplementary Table S3).
TABLE 1. Transposable elements (TEs) in PPAR regulated genesets.
To examine general tendencies of each TE class to be regulated by PPAR agonists, we used an effect size approach and estimated shifts of the t-value distribution of each TE class in PFC and liver. This analysis confirmed that SINEs and LINEs were up-regulated by fenofibrate and LINEs were up-regulated by bezafibrate in PFC (Figure 1). It also suggested that SINEs and LTRs were generally up-regulated by bezafibrate because the average t-values for those TE classes were statistically greater than 0 (p-value < 0.05), though this was a relatively weak effect. For liver, this method revealed that LINEs were down-regulated by fenofibrate, tesaglitazar, and bezafibrate. Also, SINEs were down-regulated by tesaglitazar and bezafibrate.
FIGURE 1. Global effects of PPAR agonists on transcriptomes of 4 classes of transposable elements (TEs) in mouse prefrontal cortex (PFC) (A) and liver (B). Direction and magnitude of PPAR agonist-induced changes were estimated by plotting t distributions for probes that map to a specific TE type. A t-value represents the magnitude and direction of PPAR agonist effects (t < 0 indicates downregulation; t > 0 indicates upregulation by PPAR agonist). The barplots show the average ± SEM t values (∗p < 0.05; based on one-sample t-test comparing average t-values to zero chance with Bonferroni’s correction). N = 7–10 per group.
We used qPCR to validate some of our microarray findings and specifically measure expression of the LINE1 TE that was originally detected by the ILMN_1244099 probe (Figure 2). This probe was up-regulated in PFC by fenofibrate and we confirmed this up-regulation with the qPCR assay (group difference p = 8.23E-05) (Figure 2). Moreover, qPCR detected a significant difference in the levels of the LINE1 element between the bezafibrate and saline treatment groups that the microarray was not sensitive enough to detect (p = 0.01) (Figure 2). The correlation between the microarray and qPCR data was 0.67 (p < 0.0001).
FIGURE 2. Validation of microarray results. (A) Primer Express software was used to design Taqman primers and probe sequences to detect the LINE-1 element that corresponds to Illumina probe ILMN_1244099 (left). Delta scores were calculated by subtracting the normalized saline group average from each of the normalized gene expression values for each treatment group. (B) Error bars represent the standard error of the mean. The asterisks represent significant differences at p < 0.05 between treatment and saline groups as assessed by the limma moderated t statistic for the microarray group comparisons (Ritchie et al., 2015), and a two-sided Mann Whitney test for the qPCR group comparisons. N = 7–10 per group.
Effects of PPAR Agonists on TE-Enriched Gene Co-expression Networks
We built gene-gene co-expression networks to parse the diverse transcriptional response to PPAR agonist administration in brain and liver. Network construction revealed several groups of highly correlated genes (i.e., modules), and we determined that microarray probes corresponding to TEs were enriched in several of these modules (Figure 3). Interestingly, most of the TE-enriched modules were also enriched with PPAR agonist-regulated genes (Figure 3 and Supplementary Table S5). These findings suggest that PPAR agonists regulate expression of TEs in a coordinated manner. The percentage of PPAR agonist-regulated TEs within TE-enriched modules was higher for PFC and liver than for amygdala (Figure 3). Confirming and extending possible PPAR agonist-specific and tissue-specific effects observed in the differential expression results; fenofibrate and tesaglitazar differentially regulated TEs in a TE-enriched module (tan) in PFC (fenofibrate increased while tesaglitazar decreased TEs) and fenofibrate differentially regulated TEs in the TE-enriched modules in PFC (tan and darkorange) and liver (green) (fenofibrate increased TEs in PFC and decreased in liver).
FIGURE 3. Network analysis of the mouse transcriptome and functional annotation of gene modules after PPAR agonist administration. The dendrogram of the gene network that was constructed for each tissue separately (N = 7–10 per group) (A) The x-axis corresponds to genes detected on the microarray and the y-axis to the co-expression distance between genes determined by the extent of topological overlap. Dynamic tree cutting identified modules, generally dividing them at significant branch points in the dendrogram. Genes in the modules are color-coded. Genes not assigned to a module are labeled gray. The TE-enriched modules are designated by boxes. The adjacent parentheses indicate the percentage of the respective TE class in the module that is differentially expressed (DE). If more PPAR-agonist regulated genes are in the module than expected by chance, it is indicated in the box. The arrows indicate the direction of fold-change induced by the PPAR agonist for the PPAR agonist-regulated genes in that module. F, fenofibrate, T, tesaglitazar, and B, bezafibrate. Heatmap plots of the hypergeometric p-values from the over-representation (enrichment) analysis for the differentially expressed genes (DEGs) and cell type-specific genes (B). Each row in the heatmap corresponds to one module (labeled by color on the left) and each column in the heatmap corresponds to the category being tested for over-representation. Scale bar on the right represents –log(hypergeometric p-values) used to assess statistical significance of over-representation (red, high statistical significance). Rows were arranged by hierarchical clustering. Analysis conducted and graphs made in R.
Because epigenetic modifications are an important mechanism the host uses to control and suppress transposition events, we looked for enrichment of genes involved in epigenetic functions within the modules. One TE-enriched module in each tissue was also enriched with epigenetic genes. In PFC, the tan module was enriched with epigenetic genes (Table 2) and also SINEs, LINEs, LTRs, glutamatergic neuronal genes, fenofibrate up-regulated genes, and tesaglitazar down-regulated genes (Figure 3A). In liver and amygdala, the green modules were enriched with epigenetic genes (Table 2). The green module in the amygdala was also enriched with SINEs, LINEs, LTRs, neuronal genes, and tesaglitazar down-regulated genes (Figure 3A). The green module in the liver was also enriched with SINEs, LINEs, LTRs, and the down-regulated genes of each PPAR agonist (fenofibrate, tesaglitazar, and bezafibrate) (Figure 3A). PPAR-mediated regulation of epigenetic genes in TE-enriched modules points to potential mechanisms underlying the effects of PPAR agonists on TE expression.
TABLE 2. Enrichment of genes whose products are involved in epigenetic functions within network modules.
Preservation of TE-Enriched Modules Across Tissues
To identify conserved patterns of TE expression, we assessed the overlap of gene co-expression modules between the three separate tissue-specific gene networks. As expected, the two brain networks had more conserved modules between them than between the brain and liver networks (Supplementary Figure S2). The TE modules from the various tissues form a distinct cluster within the network and are highly overlapping across tissues (red nodes in Figure 4). Interestingly, several of the TE modules were responsive to PPAR agonists, i.e., enriched with DEGs (node size in Figure 4), suggesting a specific role for PPAR agonists in regulation of TEs. The tan TE module in PFC was more regulated by PPAR agonists (i.e., larger node size) than the green TE modules in liver and amygdala, supporting our finding that TEs are highly enriched in the PPAR-regulated genesets in PFC.
FIGURE 4. A visualization of the conserved regulation of gene expression networks across the tissues analyzed (liver, amygdala, and PFC) created with Cytoscape 3.2.1. The nodes are modules of the individually constructed networks for each tissue (see Figure 3 for details). Only probes common to all tissues were used to assess module overlap. The edges are the –log(hypergeometric p-values) of the number of overlapping probes between 2 nodes (modules). Edge thickness represents the significance of the overlap between the nodes (modules) – the thicker the line, the more significant the overlap. The nodes (modules) were assessed for enrichment of transposable elements (TEs) and genes differentially expressed at p < 0.05 after systemic treatment with one of 3 PPAR agonists (fenofibrate, tesaglitazar, or bezafibrate) (see methods). Node color represents enrichment with TEs (higher intensity = more enrichment). Node size represents enrichment with PPAR agonist-regulated genes (the larger the node, the more enriched with PPAR agonist-regulated genes).
There was a highly significant overlap between the tan module in PFC, green module in liver and green module in amygdala (Figure 4 and Supplementary Figures S2, S3). These modules share an enrichment of LINEs, SINEs, tesaglitazar down-regulated genes, neuronal genes (amygdala and PFC only) and genes known to be involved in epigenetic mechanisms (Figure 3 and Table 2). There were 348 Illumina probes that were common between these 3 modules. We entered these into WebGestalt (Wang J. et al., 2017) to determine the enriched pathways, phenotypes, and ontologies. Of the 348 Illumina probes, 149 were used for analysis, as most probes mapped to TEs and therefore are unannotated and unable to be used for analysis. Two biological process ontologies (histone binding and protein complex scaffold), 2 KEGG pathways (Wnt signaling and adherens junction), and 1 Panther pathway (Cadherin signaling) emerged from this analysis (Table 3 and Supplementary Table S6). Thirteen phenotypes were enriched in the conserved TE module that fell into 3 main categories: learning and memory, craniofacial morphology, and lethality (Table 3 and Supplementary Table S6).
TABLE 3. Enriched functional categories in the conserved TE module.
Cell type is a major contributor to variability in gene expression datasets when whole tissue is analyzed (e.g., Oldham et al., 2006, 2008; Winden et al., 2009; Kang et al., 2011; Grange et al., 2014; Hawrylycz et al., 2015), which is the case here as the samples were from total homogenate. Gene coexpression network analysis, in combination with cell type specific datasets, reliably identifies gene networks specific for individual cell types (e.g., Winden et al., 2009; Grange et al., 2014). We identified that the TE-enriched conserved modules were associated with neurons in PFC and amygdala, and the liver module was not associated with any cell type (Figure 3B). This suggests that the effects of PPAR drugs on TE expression in brain are limited to neurons.
Discussion
Peroxisome proliferator activated receptor agonists have demonstrated therapeutic potential for a range of diseases that include neuropsychiatric states. Many brain disorders are associated with up-regulation of TEs and little is known about TE function in brain, let alone the pharmacological regulation of TE expression patterns. We developed an innovative application of microarray data by mapping array probes to areas of the mouse genome that contain TEs. This enabled us to observe the effects of PPAR agonist administration on the expression of 4 TE classes in mouse brain and liver and make novel hypotheses regarding TE behavior in these tissues. Given the widespread clinical use of PPAR agonists, these findings have potential implications for a significant proportion of the population.
There are a number of ways that TEs can affect biological functions. For example, TEs can modify genomic architecture via insertional mutagenesis or post-insertional rearrangements [reviewed in (Goodier and Kazazian, 2008; Guffanti et al., 2014)]. TEs can also regulate the expression of other genes by introducing transcription start sites, inducing heterochromatization, or interfering with transcription via intron retention, exonization, or polyadenylation [reviewed in (Cowley and Oakey, 2013; Guffanti et al., 2014)]. Although LINEs are more mobile in brain versus other somatic tissue (e.g., skin) (Garcia-Perez et al., 2016), retrotransposition events are relatively rare, estimated to occur in caudate and cortical human neurons (postmortem) or neural progenitor cells (in vitro) at a rate of 1 in 100 to 1 in 10,000 insertion events per neuron or neural progenitor (Coufal et al., 2009; Evrony et al., 2012). Therefore, it is unlikely that the differentially expressed TEs we observed are retrotransposed and inserted into the genome, and it is still unclear if new insertions can occur in postmitotic neurons. Although our findings are observational, and the downstream effects of the changes in TE expression after PPAR agonist treatment and the exact molecular mechanisms by which PPAR agonists modulate TE expression remain unknown, these findings suggest a novel relationship between TEs and PPAR agonists.
Gene-gene coexpression network analysis clusters genes that perform similar biological tasks or act in similar pathways, which allowed us to propose “guilt-by-association” relationships and to hypothesize either regulation or/and functional roles of TEs based on the known functions of other genes within the same module. It is important to note that enrichment analyses are used to generate hypotheses, and future experiments are required to determine how PPAR agonists regulate TE expression and any potential downstream consequences of the observed TE expression changes. In addition to being enriched with probes that are down-regulated by tesaglitazar, the probes within the conserved TE module were also enriched with LINEs, SINEs, and LTRs in both brain areas and LINEs and SINEs in liver, as well as genes related to epigenetic modification. It is known that epigenetic modifications are important for host regulation of TEs [for example, see (Slotkin and Martienssen, 2007; Maze et al., 2011; Ponomarev et al., 2012; Xie et al., 2013; Du et al., 2016)]. We found several members of the Krüppel-associated box zinc finger protein family (Zfp445, Zfp612, Zfp551, Zfp52), in the conserved TE module in both brain areas (amygdala: green, PFC: tan), which are known to play a repressive role in transposition (Wolf et al., 2015). The chromatin remodeler, Chd1, and the lysine methyltransferase, Ehmt1, were also in the conserved TE module, and both have been shown to rather specifically repress retrotransposons due to their importance in heterochromatin formation (Maksakova et al., 2013; Green et al., 2017). Surprisingly, the change in expression after PPAR agonist treatment was the same for these repressive marks and the retrotransposons in these modules, which is opposite of what would be expected. However, it could be that the changes in the heterochromatin marks are a compensatory mechanism in response to changes in TE expression. In other words, if TE expression is increased, the cell will try to inhibit TE expression by increasing expression of repressive epigenetic factors.
The relationship between chromatin modifications and TE expression has also been observed in addiction. In the first study to suggest that TEs play a functional role in AUD, Ponomarev and colleagues used postmortem brains of alcoholic and control cases to show that chronic alcohol abuse resulted in transcriptional activation of LTR-containing transposons accompanied by DNA hypomethylation in their promoters (Ponomarev et al., 2012). It was shown in mice that repeated cocaine exposure decreases heterochromatic histone methylation (H3K9me3) binding and removes epigenetic silencing of LINE-1 elements in the nucleus accumbens (Maze et al., 2011). The role of PPARs in drug dependence is reviewed in (Le Foll et al., 2013), and is discussed below.
Actions of PPAR drugs have also been associated with epigenetic factors. There has been some exploration of the relationship between histone marks and PPARs/PPAR agonists. One study provides evidence that PPAR agonists can regulate epigenetic machinery (Li et al., 2017). In this study, pioglitazone (a PPARγ agonist) suppressed platelet-derived growth factor (PDGF)-stimulated cell proliferation by regulating Hdac1 in vitro (Li et al., 2017). There is also evidence of the inverse relationship, that is, that histone marks regulate PPARs and response to PPAR agonists. Hdac5 is dephosphorylated and translocated to the nucleus by fasting-induced rises in glucagon levels where it interacts with PPARα to promote its transcriptional activity, thereby mediating hepatic fatty acid oxidation by fasting and ER stress (Qiu et al., 2017). The histone H3 lysine 4 methyltransferase, Mll4/Kmt2d, directs overnutrition-induced murine steatosis via its coactivator function for PPARγ2, whereby overnutrition activates Abl1 kinase which phosphorylates PPARγ2, and hence has enhanced interaction with Mll4 (Kim et al., 2016). DNA methylation status of certain CpGs in CD4 + T cells is associated with the inflammatory response in people undergoing a 3-week fenofibrate treatment in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study (Yusuf et al., 2017). Furthermore, Dnmt1, a maintenance enzyme responsible for de novo methylation, regulates the methylation status of the PPARγ promoter in macrophages, where more Dnmt1 correlates with less PPARγ gene expression levels in atherosclerosis patients’ monocytes (Yu et al., 2016). Given the evidence for a relationship between chromatin marks, and both PPARs, and TE expression, the strong, clear relationship we observed from our network analysis between TEs, epigenetic regulators, and PPAR agonist-regulated genes supports our approach and serves as a starting point to guide future experiments exploring the epigenetic control and the biological significance of TEs.
Our findings have two potential implications and interpretations. One possibility is that PPAR agonists are regulating TE expression levels through a PPAR-dependent transcriptional response (i.e., binding of PPARs and associated complexes to DNA in TE promoter regions). Another possibility is that PPAR agonists are regulating TE expression levels in an indirect manner, possibly by modifying the epigenetic landscape. The coexpression of TEs, epigenetic regulators, and PPAR agonist-regulated genes suggests a means by which PPAR agonist administration could affect TE expression. Namely that fenofibrate could increase TE expression by a “passive” release of chromatin-mediated gene silencing, similarly to what was observed for repeated alcohol and cocaine exposure (Maze et al., 2011; Ponomarev et al., 2012). Another line of evidence supports a PPAR-independent effect of PPAR agonist administration on TE expression levels; fenofibrate and bezafibrate, but not tesaglitazar, up-regulated SINEs in PFC. Considering that all of the PPAR agonists in this study target PPARα, it is possible that some of the observed effects on TE expression are off target (meaning they are not mediated directly by PPAR activation). Other off target effects of PPAR agonists have been noted. For example, some of the insulin-sensitizing benefits of PPAR-γ agonists are due in part to their ability to block phosphorylation and not solely to their agonist activity (Choi et al., 2010). It will be important for future experiments to investigate the effects of PPAR agonist administration on TE expression in PPARα knockout mice to determine whether these effects on TE expression are mediated through PPARα.
In our study, the directionality of TE expression indicated that fenofibrate and bezafibrate up-regulated SINEs and LINEs in PFC and down-regulated them in liver. Superficially, the finding in the PFC seems contradictory (i.e., if PPAR agonists are anti-addictive they should decrease, not increase, TEs). We probed the transcriptional response at only a single time point. Therefore, it is possible that we are measuring a compensation for decreased TEs. It could also be that the effect of tesaglitazar was masked when analyzing the full list of differentially expressed genes and TEs, because a downregulation by tesaglitazar was linked to retrotransposons by the network analysis. Another interpretation is that, similarly to TEs, PPAR agonists have both beneficial and detrimental potential. Outside of their clinical uses of improving lipid and glycemic control, PPAR agonists have other therapeutic effects, but not unaccompanied by deleterious side effects. For example, PPAR agonists can be anti-inflammatory, hepatoprotective and caridoprotective (Staels et al., 2013; Youssef and Badr, 2004; Lee and Kim, 2015), but they can also lead to edema, weight gain, and hepatocellular and bladder carcinoma, to name a few (Turner et al., 2014) (Wright et al., 2014). Given the fact that increased expression of TEs can be a mark of genomic instability and lead to a number of diseases via genomic rearrangements and changes in gene expression regulation, perhaps some of PPARs deleterious effects are mediated though their up-regulation of TEs. Thus one consideration emerging from our study is that it might be prudent to test and observe a drug’s effects on TEs as an indication of possible negative effects on the patient, especially in cases where the drug is known to alter chromatin state.
In summary, we found evidence that PPAR agonist administration can modulate retrotransposon expression levels. This is one of only a few studies to demonstrate pharmacological regulation of TEs expression. We provide a basis for further study of the mechanisms and functional consequences of increased or decreased retrotransposon expression in brain and liver after PPAR agonist administration. Understanding whether PPAR agonists are mediating deleterious effects through retrotransposons could lead to the development of treatments with improved safety profiles. More broadly, this information could lead to a better understanding of and novel treatment strategies for the many psychiatric and neurological disorders that are associated with differential expression of retrotransposons.
Author Contributions
LF, LZ, SW, CB, RH, and IP substantially contributed to the design of the work, the acquisition, analysis, and/or interpretation of data. LF drafted the work. RH and IP revised it critically. All authors provided final approval of the version to be published, and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Funding
This work was supported by the National Institutes of Health, National Institute on Alcohol Abuse and Alcoholism grants: F31 AA024332 to LF, R21 AA021462 and K01 AA017234 to IP.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2018.00331/full#supplementary-material
Footnotes
- ^http://support.illumina.com/array/downloads.html
- ^http://genome.ucsc.edu/
- ^http://www.webgestalt.org
References
Adams, M. D., Celniker, S. E., Holt, R. A., Evans, C. A., Gocayne, J. D., Amanatides, P. G., et al. (2000). The genome sequence of Drosophila melanogaster. Science 287, 2185–2195. doi: 10.1126/science.287.5461.2185
Avagliano, C., Russo, R., De Caro, C., Cristiano, C., La Rana, G., Piegari, G., et al. (2016). Palmitoylethanolamide protects mice against 6-OHDA-induced neurotoxicity and endoplasmic reticulum stress: in vivo and in vitro evidence. Pharmacol. Res. 113(Pt A), 276–289. doi: 10.1016/j.phrs.2016.09.004
Balestrieri, E., Arpino, C., Matteucci, C., Sorrentino, R., Pica, F., Alessandrelli, R., et al. (2012). HERVs expression in autism spectrum disorders. PLoS One 7:e48831. doi: 10.1371/journal.pone.0048831
Balestrieri, E., Pitzianti, M., Matteucci, C., D’Agati, E., Sorrentino, R., Baratta, A., et al. (2014). Human endogenous retroviruses and ADHD. World J. Biol. Psychiatry 15, 499–504. doi: 10.3109/15622975.2013.862345
Barbiero, J. K., Santiago, R., Tonin, F. S., Boschen, S., da Silva, L. M., Werner, M. F., et al. (2014). PPAR-alpha agonist fenofibrate protects against the damaging effects of MPTP in a rat model of Parkinson’s disease. Prog. Neuropsychopharmacol. Biol. Psychiatry 53, 35–44. doi: 10.1016/j.pnpbp.2014.02.009
Blednov, Y. A., Benavidez, J. M., Black, M., Ferguson, L. B., Schoenhard, G. L., Goate, A. M., et al. (2015). Peroxisome proliferator-activated receptors alpha and gamma are linked with alcohol consumption in mice and withdrawal and dependence in humans. Alcohol. Clin. Exp. Res. 39, 136–145. doi: 10.1111/acer.12610
Bordet, R., Ouk, T., Petrault, O., Gelé, P., Gautier, S., Laprais, M., et al. (2006). PPAR: a new pharmacological target for neuroprotection in stroke and neurodegenerative diseases. Biochem. Soc. Trans. 34(Pt 6), 1341–1346. doi: 10.1042/BST0341341
Braissant, O., Corriveau, L., Bianchi, A., Keller, J. M., Collet, P., Krémarik-Bouillaud, P., et al. (1996). Differential expression of peroxisome proliferator-activated receptors (PPARs): tissue distribution of PPAR-alpha, -beta, and -gamma in the adult rat. Endocrinology 137, 354–366. doi: 10.1210/endo.137.1.8536636
Bundo, M., Toyoshima, M., Okada, Y., Akamatsu, W., Ueda, J., Nemoto-Miyauchi, T., et al. (2014). Increased l1 retrotransposition in the neuronal genome in schizophrenia. Neuron 81, 306–313. doi: 10.1016/j.neuron.2013.10.053
Burns, K. H. (2017). Transposable elements in cancer. Nat. Rev. Cancer 17, 415–424. doi: 10.1038/nrc.2017.35
Cahoy, J. D., Emery, B., Kaushal, A., Foo, L. C., Zamanian, J. L., Christopherson, K. S., et al. (2008). A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function. J. Neurosci. 28, 264–278. doi: 10.1523/JNEUROSCI.4178-07.2008
Cartault, F., Munier, P., Benko, E., Desguerre, I., Hanein, S., Boddaert, N., et al. (2012). Mutation in a primate-conserved retrotransposon reveals a noncoding RNA as a mediator of infantile encephalopathy. Proc. Natl. Acad. Sci. U.S.A. 109, 4980–4985. doi: 10.1073/pnas.1111596109
Chen, J. M., Dimovski, A., Tchacarova, P., Petkov, G. H., and Efremov, G. D. (2008). Detection of two Alu insertions in the CFTR gene. J. Cyst. Fibros. 7, 37–43. doi: 10.1016/j.jcf.2007.04.001
Choi, J. H., Banks, A. S., Estall, J. L., Kajimura, S., Boström, P., Laznik, D., et al. (2010). Anti-diabetic drugs inhibit obesity-linked phosphorylation of PPARgamma by Cdk5. Nature 466, 451–456. doi: 10.1038/nature09291
Chuong, E. B., Elde, N. C., and Feschotte, C. (2017). Regulatory activities of transposable elements: from conflicts to benefits. Nat. Rev. Genet. 18, 71–86. doi: 10.1038/nrg.2016.139
Colle, R., de Larminat, D., Rotenberg, S., Hozer, F., Hardy, P., Verstuyft, C., et al. (2017). PPAR-gamma agonists for the treatment of major depression: a review. Pharmacopsychiatry 50, 49–55. doi: 10.1055/s-0042-120120
Coufal, N. G., Garcia-Perez, J. L., Peng, G. E., Yeo, G. W., Mu, Y., Lovci, M. T., et al. (2009). L1 retrotransposition in human neural progenitor cells. Nature 460, 1127–1131. doi: 10.1038/nature08248
Cowley, M., and Oakey, R. J. (2013). Transposable elements re-wire and fine-tune the transcriptome. PLoS Genet. 9:e1003234. doi: 10.1371/journal.pgen.1003234
Crabbe, J. C., Phillips, T. J., Harris, R. A., Arends, M. A., and Koob, G. F. (2006). Alcohol-related genes: contributions from studies with genetically engineered mice. Addict. Biol. 11, 195–269. doi: 10.1111/j.1369-1600.2006.00038.x
de Guglielmo, G., Melis, M., De Luca, M. A., Kallupi, M., Li, H. W., Niswender, K., et al. (2015). PPARgamma activation attenuates opioid consumption and modulates mesolimbic dopamine transmission. Neuropsychopharmacology 40, 927–937. doi: 10.1038/npp.2014.268
Dickerson, F., Rubalcaba, E., Viscidi, R., Yang, S., Stallings, C., Sullens, A., et al. (2008). Polymorphisms in human endogenous retrovirus K-18 and risk of type 2 diabetes in individuals with schizophrenia. Schizophr. Res. 104, 121–126. doi: 10.1016/j.schres.2008.05.005
Dickey, A. S., Pineda, V. V., Tsunemi, T., Liu, P. P., Miranda, H. C., Gilmore-Hall, S. K., et al. (2016). PPAR-delta is repressed in Huntington’s disease, is required for normal neuronal function and can be targeted therapeutically. Nat. Med. 22, 37–45. doi: 10.1038/nm.4003
Du, J., Leung, A., Trac, C., Lee, M., Parks, B. W., Lusis, A. J., et al. (2016). Chromatin variation associated with liver metabolism is mediated by transposable elements. Epigenetics Chromatin. 9:28. doi: 10.1186/s13072-016-0078-0
Du, P., Kibbe, W. A., and Lin, S. M. (2008). lumi: a pipeline for processing Illumina microarray. Bioinformatics 24, 1547–1548. doi: 10.1093/bioinformatics/btn224
Evrony, G. D., Cai, X., Lee, E., Hills, L. B., Elhosary, P. C., Lehmann, H. S., et al. (2012). Single-neuron sequencing analysis of L1 retrotransposition and somatic mutation in the human brain. Cell 151, 483–496. doi: 10.1016/j.cell.2012.09.035
Ferguson, L. B., Most, D., Blednov, Y. A., and Harris, R. A. (2014). PPAR agonists regulate brain gene expression: relationship to their effects on ethanol consumption. Neuropharmacology 86, 397–407. doi: 10.1016/j.neuropharm.2014.06.024
Fidaleo, M., Fanelli, F., Ceru, M. P., and Moreno, S. (2014). Neuroprotective properties of peroxisome proliferator-activated receptor alpha (PPARalpha) and its lipid ligands. Curr. Med. Chem. 21, 2803–2821. doi: 10.2174/0929867321666140303143455
Frank, O., Giehl, M., Zheng, C., Hehlmann, R., Leib-Mösch, C., and Seifarth, W. (2005). Human endogenous retrovirus expression profiles in samples from brains of patients with schizophrenia and bipolar disorders. J. Virol. 79, 10890–10901. doi: 10.1128/JVI.79.17.10890-10901.2005
Ganguly, A., Dunbar, T., Chen, P., Godmilow, L., and Ganguly, T. (2003). Exon skipping caused by an intronic insertion of a young Alu Yb9 element leads to severe hemophilia A. Hum. Genet. 113, 348–352. doi: 10.1007/s00439-003-0986-5
Garcia-Perez, J. L., Widmann, T. J., and Adams, I. R. (2016). The impact of transposable elements on mammalian development. Development 143, 4101–4114. doi: 10.1242/dev.132639
Goodier, J. L., and Kazazian, H. H. Jr. (2008). Retrotransposons revisited: the restraint and rehabilitation of parasites. Cell 135, 23–35. doi: 10.1016/j.cell.2008.09.022
Grange, P., Bohland, J. W., Okaty, B. W., Sugino, K., Bokil, H., Nelson, S. B., et al. (2014). Cell-type-based model explaining coexpression patterns of genes in the brain. Proc. Natl. Acad. Sci. U.S.A. 111, 5397–5402. doi: 10.1073/pnas.1312098111
Green, C. D., Huang, Y., Dou, X., Yang, L., Liu, Y., and Han, J. J. (2017). Impact of dietary interventions on noncoding RNA networks and mRNAs encoding chromatin-related factors. Cell Rep. 18, 2957–2968. doi: 10.1016/j.celrep.2017.03.001
Green, P. M., Bagnall, R. D., Waseem, N. H., and Giannelli, F. (2008). Haemophilia A mutations in the UK: results of screening one-third of the population. Br. J. Haematol. 143, 115–128. doi: 10.1111/j.1365-2141.2008.07310.x
Grygiel-Gorniak, B. (2014). Peroxisome proliferator-activated receptors and their ligands: nutritional and clinical implications–a review. Nutr. J. 13:17. doi: 10.1186/1475-2891-13-17
Guffanti, G., Gaudi, S., Fallon, J. H., Sobell, J., Potkin, S. G., Pato, C., et al. (2014). Transposable elements and psychiatric disorders. Am. J. Med. Genet. B Neuropsychiatr. Genet. 165B, 201–216. doi: 10.1002/ajmg.b.32225
Haile, C. N., and Kosten, T. A. (2017). The peroxisome proliferator-activated receptor alpha agonist fenofibrate attenuates alcohol self-administration in rats. Neuropharmacology 116, 364–370. doi: 10.1016/j.neuropharm.2017.01.007
Hancks, D. C., and Kazazian, H. H. Jr. (2012). Active human retrotransposons: variation and disease. Curr. Opin. Genet. Dev. 22, 191–203. doi: 10.1016/j.gde.2012.02.006
Hawrylycz, M., Miller, J. A., Menon, V., Feng, D., Dolbeare, T., Guillozet-Bongaarts, A. L., et al. (2015). Canonical genetic signatures of the adult human brain. Nat. Neurosci. 18, 1832–1844. doi: 10.1038/nn.4171
Howard, G., Eiges, R., Gaudet, F., Jaenisch, R., and Eden, A. (2008). Activation and transposition of endogenous retroviral elements in hypomethylation induced tumors in mice. Oncogene 27, 404–408. doi: 10.1038/sj.onc.1210631
Huang, W., Li, S., Hu, Y., Yu, H., Luo, F., Zhang, Q., et al. (2011). Implication of the env gene of the human endogenous retrovirus W family in the expression of BDNF and DRD3 and development of recent-onset schizophrenia. Schizophr. Bull. 37, 988–1000. doi: 10.1093/schbul/sbp166
Ji, M. J., Yu, X. B., Mei, Z. L., An, Y. Q., Tang, S. S., Hu, M., et al. (2015). Hippocampal PPARdelta overexpression or activation represses stress-induced depressive behaviors and enhances neurogenesis. Int. J. Neuropsychopharmacol. 19, 50–60.
Jiang, B., Wang, Y. J., Wang, H., Song, L., Huang, C., Zhu, Q., et al. (2017). Antidepressant-like effects of fenofibrate in mice via the hippocampal brain-derived neurotrophic factor signalling pathway. Br. J. Pharmacol. 174, 177–194. doi: 10.1111/bph.13668
Kang, H. J., Kawasawa, Y. I., Cheng, F., Zhu, Y., Xu, X., Li, M., et al. (2011). Spatio-temporal transcriptome of the human brain. Nature 478, 483–489. doi: 10.1038/nature10523
Karahanian, E., Quintanilla, M. E., Fernandez, K., and Israel, Y. (2014). Fenofibrate–a lipid-lowering drug–reduces voluntary alcohol drinking in rats. Alcohol 48, 665–670. doi: 10.1016/j.alcohol.2014.08.004
Karlsson, H., Bachmann, S., Schröder, J., McArthur, J., Torrey, E. F., and Yolken, R. H. (2001). Retroviral RNA identified in the cerebrospinal fluids and brains of individuals with schizophrenia. Proc. Natl. Acad. Sci. U.S.A. 98, 4634–4639. doi: 10.1073/pnas.061021998
Karlsson, H., Schröder, J., Bachmann, S., Bottmer, C., and Yolken, R. H. (2004). HERV-W-related RNA detected in plasma from individuals with recent-onset schizophrenia or schizoaffective disorder. Mol. Psychiatry 9, 12–13. doi: 10.1038/sj.mp.4001439
Kemp, D. E., Schinagle, M., Gao, K., Conroy, C., Ganocy, S. J., Ismail-Beigi, F., et al. (2014). PPAR-gamma agonism as a modulator of mood: proof-of-concept for pioglitazone in bipolar depression. CNS Drugs 28, 571–581. doi: 10.1007/s40263-014-0158-2
Killcoyne, S., Carter, G. W., Smith, J., and Boyle, J. (2009). Cytoscape: a community-based framework for network modeling. Methods Mol. Biol. 563, 219–239. doi: 10.1007/978-1-60761-175-2_12
Kim, D. H., Kim, J., Kwon, J. S., Sandhu, J., Tontonoz, P., Lee, S. K., et al. (2016). Critical roles of the histone methyltransferase MLL4/KMT2D in murine hepatic steatosis directed by ABL1 and PPARgamma2. Cell Rep. 17, 1671–1682. doi: 10.1016/j.celrep.2016.10.023
Kliewer, S. A., Forman, B. M., Blumberg, B., Ong, E. S., Borgmeyer, U., Mangelsdorf, D. J., et al. (1994). Differential expression and activation of a family of murine peroxisome proliferator-activated receptors. Proc. Natl. Acad. Sci. U.S.A. 91, 7355–7359. doi: 10.1073/pnas.91.15.7355
Le Foll, B., Di Ciano, P., Panlilio, L. V., Goldberg, S. R., and Ciccocioppo, R. (2013). Peroxisome proliferator-activated receptor (PPAR) agonists as promising new medications for drug addiction: preclinical evidence. Curr. Drug Targets 14, 768–776. doi: 10.2174/1389450111314070006
Lee, W. S., and Kim, J. (2015). Peroxisome proliferator-activated receptors and the heart: lessons from the past and future directions. PPAR Res. 2015:271983. doi: 10.1155/2015/271983
Li, F., Zhu, Y., Wan, Y., Xie, X., Ke, R., Zhai, C., et al. (2017). Activation of PPARgamma inhibits HDAC1-mediated pulmonary arterial smooth muscle cell proliferation and its potential mechanisms. Eur. J. Pharmacol. 814, 324–334. doi: 10.1016/j.ejphar.2017.08.045
Li, J., Xu, B., Chen, Z., Zhou, C., Liao, L., Qin, Y., et al. (2018). PI3K/AKT/JNK/p38 signalling pathway-mediated neural apoptosis in the prefrontal cortex of mice is involved in the antidepressant-like effect of pioglitazone. Clin. Exp. Pharmacol. Physiol. 45, 525–535. doi: 10.1111/1440-1681.12918
Li, R., Yu, C., Li, Y., Lam, T. W., Yiu, S. M., Kristiansen, K., et al. (2009). SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25, 1966–1967. doi: 10.1093/bioinformatics/btp336
Li, W., Jin, Y., Prazak, L., Hammell, M., and Dubnau, J. (2012). Transposable elements in TDP-43-mediated neurodegenerative disorders. PLoS One 7:e44099. doi: 10.1371/journal.pone.0044099
Lin, M., Pedrosa, E., Shah, A., Hrabovsky, A., Maqbool, S., Zheng, D., et al. (2011). RNA-Seq of human neurons derived from iPS cells reveals candidate long non-coding RNAs involved in neurogenesis and neuropsychiatric disorders. PLoS One 6:e23356. doi: 10.1371/journal.pone.0023356
Makela, J., Tselykh, T. V., Kukkonen, J. P., Eriksson, O., Korhonen, L. T., and Lindholm, D. (2016). Peroxisome proliferator-activated receptor-gamma (PPARgamma) agonist is neuroprotective and stimulates PGC-1alpha expression and CREB phosphorylation in human dopaminergic neurons. Neuropharmacology 102, 266–275. doi: 10.1016/j.neuropharm.2015.11.020
Maksakova, I. A., Thompson, P. J., Goyal, P., Jones, S. J., Singh, P. B., Karimi, M. M., et al. (2013). Distinct roles of KAP1, HP1 and G9a/GLP in silencing of the two-cell-specific retrotransposon MERVL in mouse ES cells. Epigenetics Chromatin. 6:15. doi: 10.1186/1756-8935-6-15
Mayfield, J., Arends, M. A., Harris, R. A., and Blednov, Y. A. (2016). Genes and alcohol consumption: studies with mutant mice. Int. Rev. Neurobiol. 126, 293–355. doi: 10.1016/bs.irn.2016.02.014
Maze, I., Feng, J., Wilkinson, M. B., Sun, H., Shen, L., and Nestler, E. J. (2011). Cocaine dynamically regulates heterochromatin and repetitive element unsilencing in nucleus accumbens. Proc. Natl. Acad. Sci. U.S.A. 108, 3035–3040. doi: 10.1073/pnas.1015483108
Melis, M., Carta, S., Fattore, L., Tolu, S., Yasar, S., Goldberg, S. R., et al. (2010). Peroxisome proliferator-activated receptors-alpha modulate dopamine cell activity through nicotinic receptors. Biol. Psychiatry 68, 256–264. doi: 10.1016/j.biopsych.2010.04.016
Moreno, S., Farioli-Vecchioli, S., and Ceru, M. P. (2004). Immunolocalization of peroxisome proliferator-activated receptors and retinoid X receptors in the adult rat CNS. Neuroscience 123, 131–145. doi: 10.1016/j.neuroscience.2003.08.064
Muotri, A. R., Marchetto, M. C., Coufal, N. G., and Gage, F. H. (2007). The necessary junk: new functions for transposable elements. Hum. Mol. Genet. 2, R159–R167. doi: 10.1093/hmg/ddm196
Muotri, A. R., Marchetto, M. C., Coufal, N. G., Oefner, R., Yeo, G., Nakashima, K., et al. (2010). L1 retrotransposition in neurons is modulated by MeCP2. Nature 468, 443–446. doi: 10.1038/nature09544
Nellaker, C., Keane, T. M., Yalcin, B., Wong, K., Agam, A., Belgard, T. G., et al. (2012). The genomic landscape shaped by selection on transposable elements across 18 mouse strains. Genome Biol. 13:R45. doi: 10.1186/gb-2012-13-6-r45
Ni, Y. F., Wang, H., Gu, Q. Y., Wang, F. Y., Wang, Y. J., Wang, J. L., et al. (2018). Gemfibrozil has antidepressant effects in mice: involvement of the hippocampal brain-derived neurotrophic factor system. J. Psychopharmacol. 32, 469–481. doi: 10.1177/0269881118762072
Oldham, M. C., Horvath, S., and Geschwind, D. H. (2006). Conservation and evolution of gene coexpression networks in human and chimpanzee brains. Proc. Natl. Acad. Sci. U.S.A. 103, 17973–17978. doi: 10.1073/pnas.0605938103
Oldham, M. C., Konopka, G., Iwamoto, K., Langfelder, P., Kato, T., Horvath, S., et al. (2008). Functional organization of the transcriptome in human brain. Nat. Neurosci. 11, 1271–1282. doi: 10.1038/nn.2207
Panlilio, L. V., Justinova, Z., and Goldberg, S. R. (2013). Inhibition of FAAH and activation of PPAR: new approaches to the treatment of cognitive dysfunction and drug addiction. Pharmacol. Ther. 138, 84–102. doi: 10.1016/j.pharmthera.2013.01.003
Panlilio, L. V., Justinova, Z., Mascia, P., Pistis, M., Luchicchi, A., Lecca, S., et al. (2012). Novel use of a lipid-lowering fibrate medication to prevent nicotine reward and relapse: preclinical findings. Neuropsychopharmacology 37, 1838–1847. doi: 10.1038/npp.2012.31
Pascual, M., Martin, J., Nieto, A., Giphart, M. J., van der Slik, A. R., de Vries, R. R., et al. (2001). Distribution of HERV-LTR elements in the 5′-flanking region of HLA-DQB1 and association with autoimmunity. Immunogenetics 53, 114–118. doi: 10.1007/s002510100307
Perron, H., Germi, R., Bernard, C., Garcia-Montojo, M., Deluen, C., Farinelli, L., et al. (2012a). Human endogenous retrovirus type W envelope expression in blood and brain cells provides new insights into multiple sclerosis disease. Mult. Scler. 18, 1721–1736. doi: 10.1177/1352458512441381
Perron, H., Hamdani, N., Faucard, R., Lajnef, M., Jamain, S., Daban-Huard, C., et al. (2012b). Molecular characteristics of human endogenous retrovirus type-W in schizophrenia and bipolar disorder. Transl. Psychiatry 2:e201. doi: 10.1038/tp.2012.125
Perron, H., Mekaoui, L., Bernard, C., Veas, F., Stefas, I., and Leboyer, M. (2008). Endogenous retrovirus type W GAG and envelope protein antigenemia in serum of schizophrenic patients. Biol. Psychiatry 64, 1019–1023. doi: 10.1016/j.biopsych.2008.06.028
Ponomarev, I., Rau, V., Eger, E. I., Harris, R. A., and Fanselow, M. S. (2010). Amygdala transcriptome and cellular mechanisms underlying stress-enhanced fear learning in a rat model of posttraumatic stress disorder. Neuropsychopharmacology 35, 1402–1411. doi: 10.1038/npp.2010.10
Ponomarev, I., Wang, S., Zhang, L., Harris, R. A., and Mayfield, R. D. (2012). Gene coexpression networks in human brain identify epigenetic modifications in alcohol dependence. J. Neurosci. 32, 1884–1897. doi: 10.1523/JNEUROSCI.3136-11.2012
Qiu, X., Li, J., Lv, S., Yu, J., Jiang, J., Yao, J., et al. (2017). HDAC5 integrates ER stress and fasting signals to regulate hepatic fatty acid oxidation. J. Lipid Res. 59, 330–338. doi: 10.1194/jlr.M080382
Reilly, M. T., Faulkner, G. J., Dubnau, J., Ponomarev, I., and Gage, F. H. (2013). The role of transposable elements in health and diseases of the central nervous system. J. Neurosci. 33, 17577–17586. doi: 10.1523/JNEUROSCI.3369-13.2013
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43:e47. doi: 10.1093/nar/gkv007
Rivera-Meza, M., Muñoz, D., Jerez, E., Quintanilla, M. E., Salinas-Luypaert, C., Fernandez, K., et al. (2017). Fenofibrate administration reduces alcohol and saccharin intake in rats: possible effects at peripheral and central levels. Front. Behav. Neurosci. 11:133. doi: 10.3389/fnbeh.2017.00133
Rolland, B., Marche, K., Cottencin, O., and Bordet, R. (2012). The PPARalpha agonist fenofibrate reduces prepulse inhibition disruption in a neurodevelopmental model of schizophrenia. Schizophr. Res. Treatment 2012:839853. doi: 10.1155/2012/839853
Rusiecki, J. A., Chen, L., Srikantan, V., Zhang, L., Yan, L., Polin, M. L., et al. (2012). DNA methylation in repetitive elements and post-traumatic stress disorder: a case-control study of US military service members. Epigenomics 4, 29–40. doi: 10.2217/epi.11.116
Scheggi, S., Melis, M., De Felice, M., Aroni, S., Muntoni, A. L., Pelliccia, T., et al. (2016). PPARalpha modulation of mesolimbic dopamine transmission rescues depression-related behaviors. Neuropharmacology 110(Pt A), 251–259. doi: 10.1016/j.neuropharm.2016.07.024
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Slokar, G., and Hasler, G. (2015). Human endogenous retroviruses as pathogenic factors in the development of schizophrenia. Front. Psychiatry 6:183. doi: 10.3389/fpsyt.2015.00183
Slotkin, R. K., and Martienssen, R. (2007). Transposable elements and the epigenetic regulation of the genome. Nat. Rev. Genet. 8, 272–285. doi: 10.1038/nrg2072
Staels, B., Rubenstrunk, A., Noel, B., Rigou, G., Delataille, P., Millatt, L. J., et al. (2013). Hepatoprotective effects of the dual peroxisome proliferator-activated receptor alpha/delta agonist. GFT505, in rodent models of nonalcoholic fatty liver disease/nonalcoholic steatohepatitis. Hepatology 58, 1941–1952. doi: 10.1002/hep.26461
Stopponi, S., de Guglielmo, G., Somaini, L., Cippitelli, A., Cannella, N., Kallupi, M., et al. (2013). Activation of PPARgamma by pioglitazone potentiates the effects of naltrexone on alcohol drinking and relapse in msP rats. Alcohol. Clin. Exp. Res. 37, 1351–1360. doi: 10.1111/acer.12091
Stopponi, S., Somaini, L., Cippitelli, A., Cannella, N., Braconi, S., Kallupi, M., et al. (2011). Activation of nuclear PPARgamma receptors by the antidiabetic agent pioglitazone suppresses alcohol drinking and relapse to alcohol seeking. Biol. Psychiatry 69, 642–649. doi: 10.1016/j.biopsych.2010.12.010
Sugino, K., Hempel, C. M., Miller, M. N., Hattox, A. M., Shapiro, P., Wu, C., et al. (2006). Molecular taxonomy of major neuronal classes in the adult mouse forebrain. Nat. Neurosci. 9, 99–107. doi: 10.1038/nn1618
Sukarova, E., Dimovski, A. J., Tchacarova, P., Petkov, G. H., and Efremov, G. D. (2001). An Alu insert as the cause of a severe form of hemophilia A. Acta Haematol. 106, 126–129. doi: 10.1159/000046602
Takahara, Y., Takahashi, M., Wagatsuma, H., Yokoya, F., Zhang, Q. W., Yamaguchi, M., et al. (2006). Gene expression profiles of hepatic cell-type specific marker genes in progression of liver fibrosis. World J. Gastroenterol. 12, 6473–6499. doi: 10.3748/wjg.v12.i40.6473
Turner, R. M., Kwok, C. S., Chen-Turner, C., Maduakor, C. A., Singh, S., and Loke, Y. K. (2014). Thiazolidinediones and associated risk of bladder cancer: a systematic review and meta-analysis. Br. J. Clin. Pharmacol. 78, 258–273. doi: 10.1111/bcp.12306
Tyagi, S., Gupta, P., Saini, A. S., Kaushal, C., and Sharma, S. (2011). The peroxisome proliferator-activated receptor: a family of nuclear receptors role in various diseases. J. Adv. Pharm. Technol. Res. 2, 236–240. doi: 10.4103/2231-4040.90879
Wang, J., Vasaikar, S., Shi, Z., Greer, M., and Zhang, B. (2017). WebGestalt: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 45, W130–W137. doi: 10.1093/nar/gkx356
Wang, Y., Zhao, W., Li, G., Chen, J., Guan, X., Chen, X., et al. (2017). Neuroprotective effect and mechanism of thiazolidinedione on dopaminergic neurons in vivo and in vitro in parkinson’s disease. PPAR Res. 2017:4089214. doi: 10.1155/2017/4089214
Warden, A., Truitt, J., Merriman, M., Ponomareva, O., Jameson, K., Ferguson, L. B., et al. (2016). Localization of PPAR isotypes in the adult mouse and human brain. Sci. Rep. 6:27618. doi: 10.1038/srep27618
Weis, S., Llenos, I. C., Sabunciyan, S., Dulay, J. R., Isler, L., Yolken, R., et al. (2007). Reduced expression of human endogenous retrovirus (HERV)-W GAG protein in the cingulate gyrus and hippocampus in schizophrenia, bipolar disorder, and depression. J. Neural. Transm. 114, 645–655. doi: 10.1007/s00702-006-0599-y
Winden, K. D., Oldham, M. C., Mirnics, K., Ebert, P. J., Swan, C. H., Levitt, P., et al. (2009). The organization of the transcriptional network in specific neuronal classes. Mol. Syst. Biol. 5:291. doi: 10.1038/msb.2009.46
Wolf, G., Greenberg, D., and Macfarlan, T. S. (2015). Spotting the enemy within: targeted silencing of foreign DNA in mammalian genomes by the Kruppel-associated box zinc finger protein family. Mob. DNA 6:17. doi: 10.1186/s13100-015-0050-8
Wright, M. B., Bortolini, M., Tadayyon, M., and Bopst, M. (2014). Minireview: challenges and opportunities in development of PPAR agonists. Mol. Endocrinol. 28, 1756–1768. doi: 10.1210/me.2013-1427
Xie, M., Hong, C., Zhang, B., Lowdon, R., Xing, X., Li, D., et al. (2013). DNA hypomethylation within specific transposable element families associates with tissue-specific enhancer landscape. Nat. Genet. 45, 836–841. doi: 10.1038/ng.2649
Yao, Y., Schröder, J., Nellåker, C., Bottmer, C., Bachmann, S., Yolken, R. H., et al. (2008). Elevated levels of human endogenous retrovirus-W transcripts in blood cells from patients with first episode schizophrenia. Genes Brain Behav. 7, 103–112.
Yoder, J. A., Walsh, C. P., and Bestor, T. H. (1997). Cytosine methylation and the ecology of intragenomic parasites. Trends Genet. 13, 335–340. doi: 10.1016/S0168-9525(97)01181-5
Yolken, R. H., Karlsson, H., Yee, F., Johnston-Wilson, N. L., and Torrey, E. F. (2000). Endogenous retroviruses and schizophrenia. Brain Res. Brain Res. Rev. 31, 193–199. doi: 10.1016/S0165-0173(99)00037-5
Youssef, J., and Badr, M. (2004). Role of peroxisome proliferator-activated receptors in inflammation control. J. Biomed. Biotechnol. 2004, 156–166. doi: 10.1155/S1110724304308065
Yu, J., Qiu, Y., Yang, J., Bian, S., Chen, G., Deng, M., et al. (2016). DNMT1-PPARgamma pathway in macrophages regulates chronic inflammation and atherosclerosis development in mice. Sci. Rep. 6:30053. doi: 10.1038/srep30053
Yusuf, N., Hidalgo, B., Irvin, M. R., Sha, J., Zhi, D., Tiwari, H. K., et al. (2017). An epigenome-wide association study of inflammatory response to fenofibrate in the genetics of lipid lowering drugs and diet network. Pharmacogenomics 18, 1333–1341. doi: 10.2217/pgs-2017-0037
Keywords: transposable elements, retrotransposons, PPAR, psychiatric, gene expression
Citation: Ferguson LB, Zhang L, Wang S, Bridges C, Harris RA and Ponomarev I (2018) Peroxisome Proliferator Activated Receptor Agonists Modulate Transposable Element Expression in Brain and Liver. Front. Mol. Neurosci. 11:331. doi: 10.3389/fnmol.2018.00331
Received: 16 May 2018; Accepted: 27 August 2018;
Published: 19 September 2018.
Edited by:
Erik Maronde, Goethe-Universität Frankfurt am Main, GermanyReviewed by:
Constance Ciaudo, ETH Zürich, SwitzerlandYukihito Ishizaka, National Center for Global Health and Medicine, Japan
Copyright © 2018 Ferguson, Zhang, Wang, Bridges, Harris and Ponomarev. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Laura B. Ferguson, laurazeavin@gmail.com Igor Ponomarev, ponomarev@utexas.edu