2,3,7,8-Tetrachlorodibenzo-<i>p</i>-dioxin (TCDD) and Polychlorinated Biphenyl Coexposure Alters the Expression Profile of MicroRNAs in the Liver Associated with Atherosclerosis

Shan, Qiuli; Qu, Fan; Chen, Ningning

doi:https://doi.org/10.1155/2020/2652756

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 2652756 | https://doi.org/10.1155/2020/2652756

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) and Polychlorinated Biphenyl Coexposure Alters the Expression Profile of MicroRNAs in the Liver Associated with Atherosclerosis

Qiuli Shan,^1,2Fan Qu,¹and Ningning Chen¹

Academic Editor: Davor Zeljezic

Received06 Dec 2019

Revised14 Jul 2020

Accepted20 Jul 2020

Published13 Aug 2020

Abstract

MicroRNAs (miRNAs) are a class of small RNAs that regulate gene expression. 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) and polychlorinated biphenyls (PCBs) are persistent organic pollutants that exist as complex mixtures in vivo. When humans are simultaneously exposed to these compounds, the development of atherosclerosis is known to be enhanced. However, the roles of miRNA in TCDD- and PCB-induced atherosclerosis are largely unknown. Therefore, the present study is aimed at elucidating the possible dysregulation of miRNAs in atherogenesis induced by coexposure to TCDD and PCBs. Eight-week-old male ApoE^-/- mice were coexposed to TCDD (15 μg/kg) and Aroclor1254 (55 mg/kg, a representative mixture of PCBs) by intraperitoneal injection four times over a 6-week period. Microarray analysis of miRNAs and mRNAs in the liver of ApoE^-/- mice with or without TCDD and Aroclor1254 coexposure was performed. We discovered that 68 miRNAs and 1312 mRNAs exhibited significant expression changes in response to TCDD and PCB coexposure and revealed that both changed miRNAs and mRNAs are involved in cardiovascular disease processes. An integrated miRNA-mRNA approach indicated that miRNA-26a-5p, miRNA-193a-3p, and miRNA-30c-5p participated in specific TCDD and Aroclor1254 coresponsive networks which are relevant to the cardiovascular system development and function network. Furthermore, our results also indicated that miRNA-130a-3p and miRNA-376a-3p were novel players in the regulation of TCDD- and Aroclor1254-induced atherosclerosis pathways. In summary, our finding provided new insights into the mechanism of atherosclerosis in response to TCDD and PCB coexposure.

1. Introduction

Cardiovascular disease (CVD) has become a major cause of morbidity and mortality in many countries [1]. Atherosclerosis is regarded as the leading cause of CVD, which is present in all humans of advanced age and continues to develop throughout life [2, 3]. Atherosclerosis is a chronic disease of the arterial wall that is characterized by the formation of atherosclerosis plaques on the inner walls of the arterial vessel [4]. Generally, the plaques are composed of the accumulated foam cells containing lipid particles, proliferated vascular smooth muscle cells (VSMCs), and extracellular matrix [5]. Moreover, monolayer cells, macrophages, and lymphocytes are also critical players in the development of plaques [6]. In this context, it is important to mention that lipids have a central role in the pathogenesis of plaques, and the liver is the main tissue which plays important roles in the regulation of lipid metabolism [3].

MicroRNAs (miRNAs) are a family of noncoding, single-stranded, 17–25-nucleotide long RNAs. miRNAs can bind to the complementary site of the target mRNA through its 3-UTR, thereby regulating protein-coding gene expression at posttranscriptional and translational levels [7]. MicroRNAs regulated diverse biological processes, including cancer and metabolic diseases [8–10]. Recently, several research reports showed that miRNAs also played a vital role in atherosclerosis [11–13]. For example, Chen and Juo [5] showed that miR-21, miR-126, miR-365, and miR-200 are related to endothelial dysfunction, while miR-26a and Let-7 are involved in the regulation of the VSMC phenotype. Moreover, Fernandez-Hernando and Baldan [14] found that miR-758 and miR-33a/b are involved in regulating the functions of macrophages. In addition, because miRNA promoter regions contain toxic response regions, miRNAs also play an important role in toxicity responses in vivo [15]. Fukushima et al. [16] showed that exposure to liver toxins such as acetaminophen or carbon tetrachloride in rat caused alteration in the expression of various miRNAs. Furthermore, Izzotti et al. [17] have demonstrated that cigarette smoking can cause changes in the miRNA expression profile in the lungs of rats. However, the relationships among toxic insult, CVD processes, and miRNAs were reported less.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) and polychlorinated biphenyls (PCBs) are persistent organic pollutants that exist as complex mixtures in vivo [18, 19]. In recent years, persistent organic pollutant- (POP-) induced CVD development has been reported in many studies [20–22]. For example, we have previously shown that in vivo coexposure of PCBs and TCDD can enhance atherosclerosis [23], but no reports have given attention to miRNA expression profiles in response to CVD induced by TCDD and PCBs. Moreover, coherent miRNA-mRNA expression changes can occur when miRNAs decrease target mRNA levels through destabilization [15, 24]. Thus, it is of great significance to identify miRNAs with differential expression from the miRNA-mRNA CVD network.

The liver not only plays a central role in the interlinked diseases of diabetes, metabolic syndrome, and CVD but is also considered as the common target organ for environmental toxins [25, 26]. Compared with other tissues, a large number of miRNAs are located within liver cells [27]. Several studies have revealed that liver miRNAs play a pivotal role in controlling the gene expression of CVD [28, 29]. For example, Goedeke et al. found miR-148a to be a key regulator of hepatic LDL-C clearance, through the direct modulation of LDLR expression. Furthermore, liver miRNAs could be considered as hepatocyte-derived extracellular vesicles that promote endothelial inflammation and atherogenesis [30]. In our previous study, we have shown that in vivo coexposure of PCBs (55 mg/kg body weight) and TCDD (15 μg/kg body weight) can enhance atherosclerosis [23]. Thus, in the present study, we sought to determine whether miRNA changes corresponding to the altered liver function in Apolipoprotein E-deficient (ApoE^-/-) mice could be detected after exposure to TCDD and Aroclor1254. MicroRNA screen and miRNA-mRNA interaction analysis were therefore conducted to show whether these miRNAs play an important role in CVD processes induced by a combination of TCDD and PCBs.

2. Materials and Methods

2.1. Chemicals

TCDD (lot number: ER052609) was obtained from Cerilliant Corporation (Cerilliant, Round Rock, TX, USA), and Aroclor1254 was obtained from AccuStandard (AccuStandard, New Haven, CT, USA). TCDD and Aroclor1254 were dissolved in a minimum amount of DMSO (Amresco, Solon, OH, USA) (0.1%) and diluted to the required concentration in corn oil (Fulinmen, China). The final concentration of DMSO in TCDD and Aroclor1254 solutions administered to mice was <0.01%. All other chemicals were obtained from Sigma-Aldrich (St. Louis, MO, USA) and were of analytical grade or of the highest grade available, unless otherwise mentioned.

2.2. Mice

ApoE^-/- mice are a good model for simulating similar atherosclerotic pathological processes in the arteries of humans [31]. Thus, in the present study, twelve male ApoE^-/- mice (8 weeks old and 20-25 g in weight) were purchased from Weitong Lihua Biocompany (Beijing, China). The mice were housed in cages with a 12-12 h light-dark cycle, 22°C ambient temperature, and a relative humidity of 40-60%. All mice had ad libitum access to water and standard chow obtained from Peking University. The research presented was performed in compliance with the national standards established by the National Council of Animal’s Use of China and was approved by the Institutional Animal Care and Use Committee of Peking University.

2.3. In Vivo TCDD and Aroclor1254 Coexposure

To determine the adverse effects of Aroclor1254 and TCDD on miRNA and mRNA profiles in the liver of ApoE^-/- mice, per the previously study [32], by using the complete randomization method, 12 animals were randomly assigned to the following two groups: (1) control group (corn oil, 0.4 ml i.p.); (2) TCDD and Aroclor1254 coexposure group. Mice were injected with the combined TCDD (15 μg/kg b.w.) and Aroclor1254 (55 mg/kg b.w.) four times, with the first two injections in the 1^st week (3 days apart) and the final two injections in the 4^th week (3 days apart) of the 6-week study. According to previous studies, TCDD (15 μg/kg b.w.) and PCB77 (49 mg/kg b.w.) led to form atherosclerotic plaques in ApoE^-/- mice [20, 33]. In our previous studies, we saw adhesion molecule expression in endothelial cells treated by 3.4 μM PCB77 [34]. Thus, 15 μg/kg b.w. of TCDD and 55 mg/kg b.w. of Aroclor1254 were chosen for the present study. This amount of POPs was based on levels that are usually found in humans after acute exposure. After 6 weeks, ApoE^-/- mice were fasted overnight and then anesthetized for tissue and plasma collection. The livers were excised, rinsed with phosphate-buffered saline (PBS, pH 7.4), and weighed. A portion of each liver sample was stored at -80°C until analysis. The entire aorta from the heart was removed and rinsed with PBS.

2.4. Oil Red O Staining

The heart tissues from ApoE^-/- mice in each group were frozen in optimal cutting temperature (OCT) compound, and then, Oil Red O staining was performed. Next, the atherosclerotic lesions were observed under a dissecting Zeiss AxioSkop 40 microscope (Zeiss, Germany). Finally, using ImageJ 1.44 software, mean lesion area was calculated from the digital pictures which are acquired by 9 to 12 serial sections from each mouse and quantified as described previously [23].

2.5. RNA Isolation

Total RNA or miRNA was isolated from the same liver and heart samples of each animal (control: , TCDD+Aroclor1254: ). Total RNA was isolated using Trizol reagent (Invitrogen, Carlsbad, CA, USA), and miRNA was isolated using a miRNA Isolation Kit (Ambion, Austin, TX, USA) according to the manufacturer’s instructions, respectively. The mRNA was purified with RNeasy Kit (Qiagen, Hilden, Germany), and the quantities and qualities of the mRNA were determined by 1% formaldehyde denaturing gel electrophoresis. The samples with bright bands of ribosomal 28S and 18S RNA in a were used for microarray analysis. The quality of miRNA was determined using an Agilent Small RNA Kit according to the manufacturer’s protocol. Both total RNA and miRNA samples were also quantified on a NanoDrop ND-2000 (Thermo Fisher Scientific, Waltham, MA, USA).

2.6. Microarray Experiments

mRNA microarrays were conducted as previously described [35]. The miRNA microarrays were performed using the Agilent Mouse miRNA array ( format, V18.0, Agilent Technologies, Inc.) at CapitalBio Corporation (Beijing, China) with each array containing probes for 1079 mature mouse miRNAs and 47 virus miRNAs. Information for the entire probe set is based on Sanger miRBase R18.0. Each miRNA was detected by probes repeated 30 times. miRNA microarray experiments were then conducted in accordance with the following description. Briefly, the miRNAs were labeled using the Agilent miRNA labeling reagent. Total RNA was dephosphorylated and ligated with pCp-Cy3, and the labeled RNA was purified and hybridized to miRNA arrays. Images were scanned with the Agilent microarray scanner (Agilent) and then analyzed using Agilent feature extraction software version 10.10.

The miRNA array data were analyzed for data summarization, normalization, and quality control using GeneSpring software V12 (Agilent). The default 90th percentile normalization method was performed for data preprocessing. To select the differentially expressed genes, we used threshold values of ≥1.5 and ≤−1.5-fold change and a Benjamini-Hochberg corrected value of 0.05.

2.7. Real-Time Quantitative PCR

Two micrograms of total RNA was used to create single-stranded cDNA using the miRNA First-Strand cDNA synthesis kit (Tiangen, Beijing, China). The forward miRNA primer was synthesized according to the sequence generated by miRBase database (http://www.mirbase.org/). The universal reverse primer was purchased from Tiangen (Beijing, China). PCR was performed using GoTaq® qPCR Master Mix (Vazyme Biotech Co., Ltd., Nanjing, China) according to the manufacturer’s instructions. Real-time PCR was performed using CFX96 Touch System (Bio-Rad, Hercules, CA, USA). PCR thermocycling conditions were 95°C for 10 min, followed by 40 cycles of 95°C for 30 s, 60°C for 30 s, and 72°C for 30 s. A melting curve program was used, and continuous fluorescence measurements were taken. The results are reported as the relative expression after normalization by β-actin, using the 2^−ΔΔCt method. Significant differences were determined using Student’s -test and considered significant when . mRNA and miRNA primers are listed in Tables 1 and 2, respectively.

2.8. Analysis of miRNAs and Their Association with Various Diseases and Networks

For heat map generation and analysis of miRNA expression, we selected miRNAs that were up- or downregulated more than 1.5-fold in ApoE^-/- mice exposed to TCDD and Aroclor1254 by using MultiExperiment Viewer software (version 4.6.0, Agilent Technologies). Next, the selected miRNAs were analyzed for their role in regulating the expression of various genes and networks using Ingenuity Pathway Analysis (IPA) software and database (version 15, Ingenuity Systems Inc., CA, USA).

2.9. miRNA-mRNA Target Interactions after TCDD and Aroclor1254 Exposure

This process is based on two principles. The first principle, which is widely accepted and experimentally validated, is the assumption that the expression of a given miRNA is inversely correlated with the expression of its mRNA targets [24, 36]. The second principle is that we select enriched gene sets from the mRNA microarray experiment using IPA. Note that other proposed pathway-centric methodologies can be substituted for IPA in this step [37]. The rationale for this approach is that biological processes result from coordinated activation of coherent sets of genes known as modules [36, 38]. The final step of this procedure involves systematic linkage of differentially expressed miRNAs to their corresponding pathways based on the presence of predicted target genes within enriched modules. A genetic pathway was constructed for the enriched “atherosclerosis signaling pathway,” and a network was constructed for the enriched “cardiovascular system development and function” module based on experimentally verified gene product interactions derived from IPA knowledge bases.

2.10. Statistical Analysis

Data are depicted as . Statistical analyses were performed using GraphPad Prism software (GraphPad Software, San Diego, CA, USA). Single comparisons were performed using unpaired Student’s -test with considered as significant.

3. Results

3.1. Increased Atherosclerosis and Altered Gene Expression Profile in ApoE^-/- Mice after Exposure to TCDD and Aroclor1254

After 6 weeks of administration, the results of Oil Red O showed that TCDD and Aroclor1254 coexposure significantly increased ( in the TCDD and Aroclor1254 coexposure group compared with in the control group, 2.31-fold, ) the atherosclerotic lesion compared to the control group (Figures 1(a)–1(c)). In addition, we investigated the changes in atherogenesis-related genes in response to coexposure to TCDD and Aroclor1254 in the heart tissue. As shown in Figure 1(d), the mRNA levels of matrix metalloproteinase-12 (Mmp12), matrix metalloproteinase-13 (Mmp13), CD36, and intercellular cell adhesion molecule-1 (Icam-1) were significantly upregulated in the TCDD and Aroclor1254 coexposure group, relative to the control (increased by 11.45-fold, 63.83-fold, 5.25-fold, and 2.25-fold, respectively, ).

(a)

(b)

(c)

(d)

Figure 1

TCDD and Aroclor1254 coexposure promoted the formation of atherosclerotic lesions. Lesions of the aortic valve lipids in the (a) control and (b) TCDD+Aroclor1254 groups were visualized by Oil Red O staining (representative cross sections, magnification ×40). (c) Values are the means of the atherosclerotic lesion from individual mouse. (d) Effect of TCDD and Aroclor1254 coexposure on the mRNA expression of genes regulating atherogenesis in ApoE^-/- mice by RT-PCR. Data are expressed as (). and , when compared between different groups.

Based on the microarrays, a total of 68 miRNAs with () and 46 miRNAs with () were found in ApoE^-/- mice exposed to TCDD and Aroclor1254 compared with mice in the control group (Figure 2(a)). Additionally, compared with the control group, 1312 mRNAs were found to be differentially expressed (; ) in the Aroclor1254 and TCDD coexposure group (data not shown) [35]. The 68 differentially expressed miRNAs were then analyzed with cluster software. A similarity measure of miRNAs in the two groups was done using the Half Square Euclidean Distance method as previously described [37]. The visualization of the miRNA cluster analysis based on similarity between groups is shown as a dendrogram (heat map, Figure 2(b)), and miRNA expression is reflected in the range of −2 to +2 (Figure 2(b)). These data demonstrated that the Aroclor1254 and TCDD coexposure significantly altered the expression profile of genes in the liver of ApoE^-/- mice.

(a)

(b)

(c)

(d)

Figure 2

Differentially expressed miRNAs in the liver of ApoE^-/- mice after exposure to TCDD and Aroclor1254 and validation by RT-PCR. (a) Venn diagram showed the downregulated miRNAs (green circle) and upregulated miRNAs (red circle) after TCDD and Aroclor1254 coexposure when compared to vehicle. (b) Heat map showed changes in miRNA expression profile in the liver after coexposure to TCDD and Aroclor1254. The expression pattern (green to red) represents the spectrum of downregulated to upregulated expression pattern of miRNAs. (c) The expression profile of the miRNAs in the liver was determined by real-time PCR. Data are depicted as from mice/group; as compared between different groups. (d) The comparison between qPCR results and microarray data is depicted.

3.2. Validation of miRNA Expression by Real-Time PCR

To validate the miRNA array data, 4 differentially expressed miRNAs (upregulated miRs: miR-34a and miR-223; downregulated miRs: miR-122 and miR-130a) were chosen for real-time PCR analysis randomly. In Figures 2(c) and 2(d), the -axis represents the miRNAs that were selected for validation, while the -axis represents the fold change in miRNA expression. Within this figure, the upward bar chart indicates upregulated miRNA expression, while the downward bar chart indicates downregulated miRNA expression. Results showed upregulated expressions of miR-34a and miR-223 in the livers treated with TCDD and Aroclor1254 when compared with control-treated livers (Figures 2(c) and 2(d)). Similarly, we observed the downregulation of miR-122 and miR-130a in TCDD- and Aroclor1254-cotreated ApoE^-/- mice when compared with controls (Figures 2(c) and 2(d)). Thus, the real-time PCR data validated the expression profiles obtained from the arrays.

3.3. Dysregulation of miRNAs and mRNAs Involved in Cardiovascular Disease Pathways after TCDD and Aroclor1254 Exposure

To further analyze the functions of differentially expressed miRNAs found in the liver of ApoE^-/- mice after TCDD and Aroclor1254 coexposure treatment, we used the IPA software and database to analyze the selected 68 miRNAs. The analysis revealed that there were as many as 20 pathways that may be affected by the dysregulated miRNAs (Figure 3 and Supplementary Table 1). There were 35 miRNAs involved in cancer, 24 in metabolic disease, 20 in cellular growth and proliferation, and 17 in hematological disease (Figure 3). Similarly, as shown in Figure 3, several differentially expressed miRNAs were involved in other pathways such as inflammatory disease (16 miRNAs), immunological disease (15 miRNAs), inflammatory response (13 miRNAs), cell death and survival (12 miRNAs), cellular movement (10 miRNAs), and cell cycle (8 miRNAs). Remarkably, there were five miRNAs involved in cardiovascular disease, including miR-16-5p, miR-30c-5p, miR-181a-5p, miR-26a-5p, and miR-130a-3p (data not shown).

Similarly, 1312 differentially expressed mRNAs were analyzed using IPA software. As shown in Figure 4, there were more than 20 pathways that may be affected by the dysregulation of mRNAs. Figure 4 presented the top 22 pathways according to the canonical pathways in IPA. The prominent passage regulation occurred in LPS/IL-1-mediated inhibition of RXR function, IL-10 signaling, NRF2-mediated oxidative stress response, and atherosclerosis signaling.

Figure 4

Significant canonical pathways that responded to TCDD and Aroclor1254 coexposure according to IPA software analysis. Up- and downregulated (1312 genes, , ) mRNAs were analyzed using IPA software. The graph showed the various pathways regulated by coexposure to TCDD and Aroclor1254. The value for each pathway was indicated by the line and is expressed as −1 times the log of the value. The bar represented the ratio of the number of genes differentially expressed in a given pathway divided by the total number of genes that make up that pathway.

3.4. miRNA-mRNA Interactions Participate in Cardiovascular Disease

mRNA network analysis suggested that TCDD and Aroclor1254 coexposure could change cardiovascular system development and function network and atherosclerosis signaling pathway (Supplementary Table 2). To evaluate the miRNAs which participated in cardiovascular disease induced by TCDD and Aroclor1254 coexposure, a miRNA-mRNA network approach was applied. We used the differentially expressed miRNAs and mRNAs that participated in TCDD- and Aroclor1254-responsive networks and pathways as input data, and the IPA software output is shown in Figures 5 and 6.

Figure 5

MicroRNA-mRNA interaction related to the cardiovascular system development and function network. miRNAs were used to link mRNA network related with cardiovascular system development and function to demonstrate the regulation role through IPA. Green and blue icons meant the downregulated mRNAs and miRNAs, respectively. Red and yellow icons meant the upregulated mRNAs and miRNAs, respectively. The solid arrow between miRNA and mRNAs denoted the interaction between miRNAs and mRNAs. Significantly upregulated or downregulated mRNAs: ( value < 0.05) or ( value < 0.05). Significantly upregulated genes or downregulated miRNAs: ( value < 0.05) or ( value < 0.05), respectively.

Figure 6

MicroRNA-mRNA networks in the atherosclerosis signaling pathway after TCDD and Aroclor1254 coexposure in ApoE^-/- mice. The atherosclerosis signaling pathway was significantly altered after treatment with TCDD and Aroclor1254. The mRNAs with red meant altered genes after treatment. The solid arrow between miRNA and mRNAs denoted the interaction between miRNAs and mRNAs. Significantly upregulated or downregulated mRNAs: ( value < 0.05) or ( value < 0.05). Significantly upregulated genes or downregulated miRNAs: ( value < 0.05) or ( value < 0.05), respectively.

As shown in Figure 5, the network for cardiovascular system development and function revealed the participation of miRNAs, such as miR-26a-5p, miR-17-5p, miR-30c-5p, miR-18-5P, miR-532-3P, miR-532-5p, miR-30a-3p, miR-193a-3p, miR-34a-5p, and miR-100-5p, and mRNAs including cyclin D1 (CCND1), transcription factor 4 (TCF4), connective tissue growth factor (CTGF), and transforming growth factor, beta receptor II (TGFBR2). Figure 6 presents the effects of Aroclor1254 and TCDD coexposure on miRNAs and mRNAs associated with the atherosclerosis signaling pathway. We can see from this figure that there are six kinds of cell (monocyte, endothelial cell, macrophage, foam cell, T lymphocyte, and smooth muscle cell) in ApoE^-/- mice that are associated with the atherosclerosis induced by Aroclor1254 and TCDD coexposure. miR-376a-3p, miR-484, miR-7a-5p, miR-532-5p, miR-130a-3p, miR-718, miR-122-5p, and miR-192-5p showed interactions with mRNAs that participated in atherosclerosis signaling pathways.

4. Discussion

The mechanism of TCDD- and PCB-induced atherosclerosis has been well characterized by many studies [21, 23, 39, 40]. Emerging data continue to demonstrate the participation of miRNAs in atherosclerosis [41–43] and in biological responses induced by environmental pollutants [44–46]. However, there have been no previous reports on miRNA expression in TCDD- and PCB-induced atherosclerosis. In our previous study, we found that TCDD and PCB coexposure in ApoE^-/- mice enhanced atherosclerosis development [23]. Moreover, a growing body of evidence has shown that individual exposure to TCDD or PCBs in vivo or in vitro models caused liver disease while enhancing atherosclerosis [33, 47]. Importantly, liver disease activates the expression of inflammatory factors and oxidative stress in the development of the course, further stimulating the formation of atherosclerosis [48]. Therefore, exploring the disease-specific miRNA-mRNA networks in the liver can serve as a novel mechanism for TCDD- and PCB-induced atherosclerosis in an in vivo animal model.

The cluster analysis data of miRNAs showed that TCDD and Aroclor1254 treatment caused significant changes in miRNA expression profiles in the livers of ApoE^-/- mice compared with livers of control-treated mice. Of the miRNAs screened, the expressions of 68 miRNAs were changed by more than 1.5-fold and those of 46 miRNAs were altered twofold or more. We further validated the expression profile of selected miRNAs by performing real-time PCR. All the miRNAs that we analyzed by real-time PCR corroborated the data obtained from the miRNA array analysis. The IPA pathway analysis revealed that there are five miRNAs, miR-16-5p, miR-30c-5p, miR-181a-5p, miR-26a-5p, and miR-130a-3p, that play a role in the cardiovascular disease induced by TCDD and Aroclor1254 coexposure. Moreover, the mRNA pathway analysis indicated that the atherosclerosis signaling pathway was significantly activated by TCDD and Aroclor1254 exposure.

It is widely accepted that expression changes in a small number of miRNAs have strong effects on mRNA levels. Therefore, to further determine which miRNAs regulate the atherosclerosis process, we investigated miRNAs and mRNAs that may be involved in miRNA-mRNA interactions and associated with the cardiovascular system development and function network and the atherosclerosis pathway.

4.1. Cardiovascular System Development and Function Network

Some of the dysregulated miRNAs could interact with multiple CVD-related mRNAs, thus playing a prominent role in the cardiovascular system development and function network. In the present study, we found that the downregulation of miRNA-26a-5p was associated with the upregulation of TGFBR2 in the liver of ApoE^-/^- mice after coexposure to TCDD and Aroclor1254. TGFBR2 is one of CVD development effectors, which plays an essential role in vascular development [49]. Previous study showed that selective deletion of TGFBR2 in smooth muscle caused the vascular deformities and embryonic lethality [50]. Thus, we hypothesized that miRNA-26a-5p may play an important role in the regulation of CVD induced by TCDD and Aroclor1254 coexposure treatments by regulating the expression of its target gene TGFBR2. Additionally, TGFBR2 can also regulate cyclin-dependent kinase (CDK) activity by interacting with CCND1, thus playing an important role in cardiovascular disease [51, 52]. In this study, the cardiovascular system development and function network showed that the expression of CCND1 was upregulated by miRNA-193a-3p [53]. Moreover, the interactions between miRNA-30c-5p and CTGF were also detected [54]. It has been reported that the mRNA level of CTGF was elevated in the aorta of diabetic ApoE^-/- mice, which was an important factor in the vascular ECM accumulation [55]. These gene expression changes together suggested that miRNA-26a-5p, miRNA-193a-3p, and miRNA-30c-5p may play important roles in the promotion of the CVD response following TCDD and Aroclor1254 exposure.

4.2. Atherosclerosis Pathway

Atherosclerosis is a major risk factor for CVD and is characterized by the accumulation of lipids and fibrous elements in the large arteries [56]. In this study, atherosclerosis pathway analysis showed that VSMCs, lymphocytes, endothelial cells, monolayer cells, macrophage, and foam cells are critical players in the atherosclerosis pathway induced by TCDD and Aroclor1254 coexposure in ApoE^-/- mice. There are several miRNAs and important genes together with validated and/or predicted microRNA-mRNA interactions involved in the atherosclerosis pathway. An important endothelial cell mRNA is ICAM-1, which is a possible target of miR-130a-3p. ICAM-1, as an adhesion molecule of the immunoglobulin superfamily, is critically involved in both cancer metastasis and atherogenesis. There are several pieces of evidence demonstrating that PCBs induce ICAM-1 expression. For example, PCB104 exposure in human vessel endothelial cells could upregulate the level of ICAM-1 [57], while exposure to PCB153 not only increased the ICAM-1 expression but also elevated adhesion of leukocytes to brain endothelial cells [58]. In the present study, our gene chip analysis showed that TCDD and PCB coexposure could upregulate the mRNA level of ICAM-1 and that ICAM-1 was regulated by miR-130a-3p. Notably, in a previous study, miR-130a was suggested to play an important role in maintaining normal endothelial progenitor cell (EPC) function, and dysfunction of EPCs contributes to diabetic vascular disease [59]. Other atherosclerosis risk factors, such as Mmp13, M-CSF, sdf-1α, and Integrinβ2, are also targets of miR-130a-3p possibly; however, only Mmp13 and ICAM-1 were upregulated in accordance with downregulated miR-130a-3p. Another important player involved in the macrophage cell functional regulation is CD36, and its mRNA is a possible target of miR-376a-3p. CD36 is a lipid scavenger receptor located on the surface of macrophages that can bind oxidized low-density lipid from plasma. CD36 can be induced by TCDD [60] and PCBs [61]. It has been reported that CD36 is a trigger for the inflammation response and influences the development of foam cells to activate atherosclerosis [62]. Therefore, miR-376a-3p interaction with CD36 may be an important atherosclerotic mechanism induced by TCDD and PCBs. Although other miRNA-mRNA interactions were displayed in the atherosclerosis pathway, the changes in mRNA levels were not consistent with miRNA regulation. Overall, those results strengthened the claim that miRNAs participate in atherosclerosis pathway regulation and potentially have roles in CVD processes induced by TCDD and PCBs.

5. Conclusion

In summary, we have detected and validated the expression levels of mRNAs and miRNAs after TCDD and PCB coexposure in ApoE^-/- mice. Our results suggested that miRNAs played an important role in CVD processes, such as the CVD development and function network and the atherosclerosis pathway. Furthermore, this study presented a possible miRNA-mRNA network involved in these pathways, thereby providing new insights into the mechanism of CVD in response to TCDD and PCB coexposure. Additionally, our findings established a basis for the development of combined miRNA- and mRNA-based approaches for predicting the development of CVD in response to various compounds.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (NSFC) (grant number 21507040) and the Natural Science Foundation of Shandong Province (grant number BS2015SW023).

Supplementary Materials

Two tables were contained in our supplementary material file; the contents of the two tables were TCDD and Aroclor1254 co-exposure-regulated miRNAs and their association with diseases and modulation of component expression of the atherosclerosis signaling pathway in ApoE^-/- mice coexposed to TCDD and PCBs. (Supplementary Materials)

References

G. A. Roth, C. Johnson, A. Abajobir et al., “Global, regional, and national burden of cardiovascular diseases for 10 causes, 1990 to 2015,” Journal of the American College of Cardiology, vol. 70, no. 1, pp. 1–25, 2017.
View at: Publisher Site | Google Scholar
P. Abizanda, P. Atienzar, L. Casado et al., “Cardiovascular risk factors are associated with subclinical atherosclerosis in high functioning older adults,” Maturitas, vol. 67, no. 1, pp. 54–59, 2010.
View at: Publisher Site | Google Scholar
X. H. Yu, D. W. Zhang, X. L. Zheng, and C. K. Tang, “Cholesterol transport system: an integrated cholesterol transport model involved in atherosclerosis,” Progress in Lipid Research, vol. 73, pp. 65–91, 2019.
View at: Publisher Site | Google Scholar
E. Falk, “Pathogenesis of atherosclerosis,” Journal of the American College of Cardiology, vol. 47, 8 Suppl, pp. C7–12, 2006.
View at: Publisher Site | Google Scholar
K.-C. Chen and S.-H. H. Juo, “MicroRNAs in atherosclerosis,” The Kaohsiung Journal of Medical Sciences, vol. 28, no. 12, pp. 631–640, 2012.
View at: Publisher Site | Google Scholar
P. Libby, “Inflammation in atherosclerosis,” Nature, vol. 420, no. 6917, pp. 868–874, 2002.
View at: Publisher Site | Google Scholar
V. Ambros, “The functions of animal microRNAs,” Nature, vol. 431, no. 7006, pp. 350–355, 2004.
View at: Publisher Site | Google Scholar
P. Ahmad, J. Sana, M. Slavik, P. Slampa, P. Smilek, and O. Slaby, “MicroRNAs involvement in radioresistance of head and neck cancer,” Disease Markers, vol. 2017, Article ID 8245345, 8 pages, 2017.
View at: Publisher Site | Google Scholar
J. Xu, J. Li, and Z.-j. Liu, “MicroRNAs in the occurrence and development of primary hepatocellular carcinoma,” Advances in Clinical and Experimental Medicine, vol. 25, no. 5, pp. 971–975, 2016.
View at: Publisher Site | Google Scholar
A. Matsuda, I. K. Yan, C. Foye, M. Parasramka, and T. Patel, “MicroRNAs as paracrine signaling mediators in cancers and metabolic diseases,” Best Practice & Research Clinical Endocrinology & Metabolism, vol. 30, no. 5, pp. 577–590, 2016.
View at: Publisher Site | Google Scholar
H. Giral, A. Kratzer, and U. Landmesser, “MicroRNAs in lipid metabolism and atherosclerosis,” Best Practice & Research Clinical Endocrinology & Metabolism, vol. 30, no. 5, pp. 665–676, 2016.
View at: Publisher Site | Google Scholar
B. Laffont and K. J. Rayner, “MicroRNAs in the pathobiology and therapy of atherosclerosis,” Canadian Journal of Cardiology, vol. 33, no. 3, pp. 313–324, 2017.
View at: Publisher Site | Google Scholar
A. Churov, V. Summerhill, A. Grechko, V. Orekhova, and A. Orekhov, “MicroRNAs as potential biomarkers in atherosclerosis,” International Journal of Molecular Sciences, vol. 20, no. 22, p. 5547, 2019.
View at: Publisher Site | Google Scholar
C. Fernández-Hernando and Á. Baldán, “MicroRNAs and cardiovascular disease,” Current Genetic Medicine Reports, vol. 1, no. 1, pp. 30–38, 2013.
View at: Publisher Site | Google Scholar
D. Lizarraga, S. Gaj, K. J. Brauers, L. Timmermans, J. C. Kleinjans, and J. H. M. van Delft, “Benzo[a]pyrene-induced changes in microRNA-mRNA networks,” Chemical Research in Toxicology, vol. 25, no. 4, pp. 838–849, 2012.
View at: Publisher Site | Google Scholar
T. Fukushima, Y. Hamada, H. Yamada, and I. Horii, “Changes of micro-RNA expression in rat liver treated by acetaminophen or carbon tetrachloride--regulating role of micro-RNA for RNA expression,” The Journal of Toxicological Sciences, vol. 32, no. 4, pp. 401–409, 2007.
View at: Publisher Site | Google Scholar
A. Izzotti, G. A. Calin, P. Arrigo, V. E. Steele, C. M. Croce, and S. De Flora, “Downregulation of microRNA expression in the lungs of rats exposed to cigarette smoke,” vol. 23, no. 3, pp. 806–812, 806.
View at: Publisher Site | Google Scholar
K. T. Brunner, A. Pedley, J. M. Massaro, U. Hoffmann, E. J. Benjamin, and M. T. Long, “Increasing Liver Fat Is Associated With Incident Cardiovascular Risk Factors,” Clinical Gastroenterology and Hepatology, vol. 18, no. 8, pp. 1884–1886, 2020.
View at: Publisher Site | Google Scholar
C. Chirollo, M. Ceruso, T. Pepe et al., “Levels and congeners distribution of dioxins, furans and dioxin-like PCBs in buffaloes adipose tissues sampled in vivo and milk,” CyTA - Journal of Food, vol. 16, no. 1, pp. 1109–1114, 2018.
View at: Publisher Site | Google Scholar
D. Wu, N. Nishimura, V. Kuo et al., “Activation of aryl hydrocarbon receptor induces vascular inflammation and promotes atherosclerosis in apolipoprotein E-/- mice,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 31, no. 6, pp. 1260–1267, 2011.
View at: Publisher Site | Google Scholar
M. C. Petriello, J. A. Brandon, J. Hoffman et al., “Dioxin-like PCB 126 inflammation and atherosclerosis in LDL Receptor-DeficientMice,” Toxicological Sciences, vol. 162, no. 2, pp. 548–558, 2018.
View at: Publisher Site | Google Scholar
B. Yang, Y. Wang, Q. Qin et al., “Polychlorinated quinone Macrophage-Derivedfoam Formation,” Chemical Research in Toxicology, vol. 32, no. 12, pp. 2422–2432, 2019.
View at: Publisher Site | Google Scholar
Q. Shan, J. Wang, F. Huang, X. Lv, M. Ma, and Y. Du, “Augmented atherogenesis in ApoE-null mice co-exposed to polychlorinated biphenyls and 2,3,7,8-tetrachlorodibenzo-p-dioxin,” Toxicology and Applied Pharmacology, vol. 276, no. 2, pp. 136–146, 2014.
View at: Publisher Site | Google Scholar
H. Guo, N. T. Ingolia, J. S. Weissman, and D. P. Bartel, “Mammalian microRNAs predominantly act to decrease target mRNA levels,” Nature, vol. 466, no. 7308, pp. 835–840, 2010.
View at: Publisher Site | Google Scholar
J. Luo, L. Xu, J. Li, and S. Zhao, “Nonalcoholic fatty liver disease as a potential risk factor of cardiovascular disease,” European Journal of Gastroenterology & Hepatology, vol. 27, no. 3, pp. 193–199, 2015.
View at: Publisher Site | Google Scholar
P. A. Kern, R. B. Fishman, W. Song, A. D. Brown, and V. Fonseca, “The effect of 2,3,7,8-tetrachlorodibenzo- _p_ -dioxin (TCDD) on oxidative enzymes in adipocytes and liver,” Toxicology, vol. 171, no. 2-3, pp. 117–125, 2002.
View at: Publisher Site | Google Scholar
X. M. Chen, “MicroRNA signatures in liver diseases,” World Journal of Gastroenterology, vol. 15, no. 14, pp. 1665–1672, 2009.
View at: Publisher Site | Google Scholar
P. Willeit, P. Skroblin, S. Kiechl, C. Fernandez-Hernando, and M. Mayr, “Liver microRNAs: potential mediators and biomarkers for metabolic and cardiovascular disease?” European Heart Journal, vol. 37, no. 43, pp. 3260–3266, 2016.
View at: Publisher Site | Google Scholar
L. Goedeke, N. Rotllan, A. Canfrán-Duque et al., “MicroRNA-148a regulates LDL receptor and ABCA1 expression to control circulating lipoprotein levels,” Nature Medicine, vol. 21, no. 11, pp. 1280–1289, 2015.
View at: Publisher Site | Google Scholar
F. Jiang, Q. Chen, W. Wang, Y. Ling, Y. Yan, and P. Xia, “Hepatocyte-derived extracellular vesicles promote endothelial inflammation and atherogenesis via microRNA-1,” Journal of Hepatology, vol. 72, no. 1, pp. 156–166, 2020.
View at: Publisher Site | Google Scholar
Y. Nakashima, A. S. Plump, E. W. Raines, J. L. Breslow, and R. Ross, “ApoE-deficient mice develop lesions of all phases of atherosclerosis throughout the arterial tree,” Arteriosclerosis and Thrombosis: A Journal of Vascular Biology, vol. 14, no. 1, pp. 133–140, 1994.
View at: Publisher Site | Google Scholar
M. F. W. Festing, “Randomized block experimental designs can increase the power and reproducibility of laboratory animal experiments,” ILAR Journal, vol. 55, no. 3, pp. 472–476, 2014.
View at: Publisher Site | Google Scholar
V. Arsenescu, R. I. Arsenescu, V. King, H. Swanson, and L. A. Cassis, “Polychlorinated biphenyl-77 induces adipocyte differentiation and proinflammatory adipokines and promotes obesity and atherosclerosis,” Environmental Health Perspectives, vol. 116, no. 6, pp. 761–768, 2008.
View at: Publisher Site | Google Scholar
J. Wang, X. Lv, and Y. Du, “Inflammatory response and insulin signaling alteration induced by PCB77,” Journal of Environmental Sciences (China), vol. 22, no. 7, pp. 1086–1090, 2010.
View at: Publisher Site | Google Scholar
Q. Shan, F. Huang, J. Wang, and Y. Du, “Effects of co-exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin and polychlorinated biphenyls on nonalcoholic fatty liver disease in mice,” Environmental Toxicology, vol. 30, no. 12, pp. 1364–1374, 2015.
View at: Publisher Site | Google Scholar
S. A. Gharib, A. Khalyfa, A. Abdelkarim, B. Bhushan, and D. Gozal, “Integrative miRNA-mRNA profiling of adipose tissue unravels transcriptional circuits induced by sleep fragmentation,” PLoS One, vol. 7, no. 5, p. e37669, 2012.
View at: Publisher Site | Google Scholar
N. P. Singh, U. P. Singh, H. Guan, P. Nagarkatti, and M. Nagarkatti, “Prenatal exposure to TCDD triggers significant modulation of microRNA expression profile in the thymus that affects consequent gene expression,” PLoS One, vol. 7, no. 9, p. e45054, 2012.
View at: Publisher Site | Google Scholar
M. Vidal, M. E. Cusick, and A. L. Barabasi, “Interactome networks and human disease,” Cell, vol. 144, no. 6, pp. 986–998, 2011.
View at: Publisher Site | Google Scholar
Y. Yu, Q. Liu, S. Guo et al., “2, 3, 7, 8-Tetrachlorodibenzo-p-dioxin promotes endothelial cell apoptosis through activation of EP3/p38MAPK/Bcl-2 pathway,” Journal of Cellular and Molecular Medicine, vol. 21, no. 12, pp. 3540–3551, 2017.
View at: Publisher Site | Google Scholar
C. Wang, M. C. Petriello, B. Zhu, and B. Hennig, “PCB 126 induces monocyte/macrophage polarization and inflammation through AhR and NF-κB pathways,” Toxicology and Applied Pharmacology, vol. 367, pp. 71–81, 2019.
View at: Publisher Site | Google Scholar
M. C. Blaser and E. Aikawa, “Differential miRNA underpins harmful and roles for vesicles in Atherogenesis,” Circulation Research, vol. 124, no. 4, pp. 467–469, 2019.
View at: Publisher Site | Google Scholar
M. Fang, Y. Li, Y. Wu, Z. Ning, X. Wang, and X. Li, “miR-185 silencing promotes the progression of atherosclerosis via targeting stromal interaction molecule 1,” Cell Cycle, vol. 18, no. 6-7, pp. 682–695, 2019.
View at: Publisher Site | Google Scholar
X. Shi and X. Chen, “Effect of microRNA-370 on coronary atherosclerosis and its underlying mechanism,” Experimental and Therapeutic Medicine, vol. 17, no. 1, pp. 115–122, 2019.
View at: Publisher Site | Google Scholar
J. Krauskopf, T. M. de Kok, D. G. Hebels et al., “MicroRNA profile for health risk assessment: Environmental exposure to persistent organic pollutants strongly affects the human blood microRNA machinery,” Scientific Reports, vol. 7, no. 1, p. 9262, 2017.
View at: Publisher Site | Google Scholar
M. A. Ribeiro, M. S. Estill, G. J. Fernandez, L. N. Moraes, S. A. Krawetz, and W. R. Scarano, “Integrative transcriptome and microRNome analysis identifies dysregulated pathways in human Sertoli cells exposed to TCDD,” Toxicology, vol. 409, pp. 112–118, 2018.
View at: Publisher Site | Google Scholar
C. A. Waugh, A. Arukwe, and V. L. B. Jaspers, “Deregulation of microRNA-155 and its transcription factor NF-kB by polychlorinated biphenyls during viral infections,” Apmis, vol. 126, no. 3, pp. 234–240, 2018.
View at: Publisher Site | Google Scholar
S. G. Han, S. S. Han, M. Toborek, and B. Hennig, “EGCG protects endothelial cells against PCB 126-induced inflammation through inhibition of AhR and induction of Nrf2-regulated genes,” Toxicology and Applied Pharmacology, vol. 261, no. 2, pp. 181–188, 2012.
View at: Publisher Site | Google Scholar
N. Alkhouri, C. Carter-Kent, M. Elias, and A. E. Feldstein, “Atherogenic dyslipidemia and cardiovascular risk in children with nonalcoholic fatty liver disease,” Clinical Lipidology, vol. 6, no. 3, pp. 305–314, 2011.
View at: Publisher Site | Google Scholar
W. Li, Q. Li, Y. Jiao et al., “Tgfbr2 disruption in postnatal smooth muscle impairs aortic wall homeostasis,” vol. 124, no. 2, pp. 755–767, 755.
View at: Publisher Site | Google Scholar
R. L. C. Carvalho, F. Itoh, M.-J. Goumans et al., “Compensatory signalling induced in the yolk sac vasculature by deletion of TGF receptors in mice,” Journal of Cell Science, vol. 120, no. 24, pp. 4269–4277, 2007.
View at: Publisher Site | Google Scholar
J. Massague, “TGFβ signalling in context,” Nature Reviews Molecular Cell Biology, vol. 13, no. 10, pp. 616–630, 2012.
View at: Publisher Site | Google Scholar
C. Ihling, K. Technau, V. Gross, J. Schulte-Monting, A. M. Zeiher, and H. E. Schaefer, “Concordant upregulation of type II-TGF-β-receptor, the cyclin-dependent kinases inhibitor P27^Kip1 and cyclin E in human atherosclerotic tissue: implications for lesion cellularity,” Atherosclerosis, vol. 144, no. 1, pp. 7–14, 1999.
View at: Publisher Site | Google Scholar
H. Nakano, Y. Yamada, T. Miyazawa, and T. Yoshida, “Gain-of-function microRNA screens identify miR-193a regulating proliferation and apoptosis in epithelial ovarian cancer cells,” International Journal of Oncology, vol. 42, no. 6, pp. 1875–1882, 2013.
View at: Publisher Site | Google Scholar
R. F. Duisters, A. J. Tijsen, B. Schroen et al., “miR-133 and miR-30 regulate connective tissue growth factor: implications for a role of microRNAs in myocardial matrix remodeling,” Circulation Research, vol. 104, no. 2, pp. 170–178, 2009.
View at: Publisher Site | Google Scholar
Y. Pham, Y. Tu, T. Wu et al., “Cell division autoantigen 1 plays a profibrotic role by modulating downstream signalling of TGF- β in a murine diabetic model of atherosclerosis,” Diabetologia, vol. 53, no. 1, pp. 170–179, 2010.
View at: Publisher Site | Google Scholar
P. Libby, P. M. Ridker, and G. K. Hansson, “Progress and challenges in translating the biology of atherosclerosis,” Nature, vol. 473, no. 7347, pp. 317–325, 2011.
View at: Publisher Site | Google Scholar
W. Choi, S. Y. Eum, Y. W. Lee, B. Hennig, L. W. Robertson, and M. Toborek, “PCB 104-induced proinflammatory reactions in human vascular endothelial cells: relationship to cancer metastasis and atherogenesis,” Toxicological Sciences, vol. 75, no. 1, pp. 47–56, 2003.
View at: Publisher Site | Google Scholar
S. Y. Eum, I. Andras, B. Hennig, and M. Toborek, “NADPH oxidase and lipid raft-associated redox signaling are required for PCB153-induced upregulation of cell adhesion molecules in human brain endothelial cells,” Toxicology and Applied Pharmacology, vol. 240, no. 2, pp. 299–305, 2009.
View at: Publisher Site | Google Scholar
S. Meng, J. Cao, X. Zhang et al., “Downregulation of microRNA-130a contributes to endothelial progenitor cell dysfunction in diabetic patients via its target Runx3,” PLoS One, vol. 8, no. 7, p. e68611, 2013.
View at: Publisher Site | Google Scholar
M. M. Angrish, B. D. Mets, A. D. Jones, and T. R. Zacharewski, “Dietary fat is a lipid source in 2,3,7,8-Tetrachlorodibenzo-ρ-Dioxin (TCDD)-elicited hepatic steatosis in C57BL/6 mice,” Toxicological Sciences, vol. 128, no. 2, pp. 377–386, 2012.
View at: Publisher Site | Google Scholar
I. K. Lai, K. Dhakal, G. S. Gadupudi et al., “N-acetylcysteine (NAC) diminishes the severity of PCB 126-induced fatty liver in male rodents,” Toxicology, vol. 302, no. 1, pp. 25–33, 2012.
View at: Publisher Site | Google Scholar
L. Chavez-Sanchez, M. G. Garza-Reyes, J. E. Espinosa-Luna, K. Chavez-Rueda, M. V. Legorreta-Haquet, and F. Blanco-Favela, “The role of TLR2, TLR4 and CD36 in macrophage activation and foam cell formation in response to oxLDL in humans,” Human Immunology, vol. 75, no. 4, pp. 322–329, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Qiuli Shan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

504

Downloads

1011

Citations