Molecular Vision 2018; 24:633-646
Received 26 February 2018 | Accepted 21 September 2018 | Published 24 September 2018
Zeljka Smit-McBride,1 Johnny Nguyen,1 Garrett W. Elliott,1 Zhe Wang,2 Ryan A. McBride,3 Anthony T. Nguyen,1 Sharon L. Oltjen,1 Glenn Yiu,1 Sara M. Thomasy,1,4 Kent E. Pinkerton,5,6 Eugene S. Lee,7 David Cunefare,8 Sina Farsiu,8,9 Lawrence S. Morse1
1Department of Ophthalmology & Vision Science, School of Medicine, University of California, Davis, Davis, CA; 2School of Veterinary Medicine, University of California, Davis, Davis, CA; 3Department of Cell Biology and Human Anatomy, School of Medicine, University of California, Davis, Davis, CA; 4Department of Surgical and Radiological Sciences, School of Veterinary Medicine, University of California, Davis, Davis, CA; 5California National Primate Research Center, University of California, Davis, Davis, CA; 6Center for Health & the Environment (CHE), University of California, Davis, Davis, CA; 7Sacramento VA Medical Center, Mather, CA; 8Department of Biomedical Engineering, Duke University, Durham, NC; 9Department of Ophthalmology, Duke University, Durham, NC
Correspondence to: Zeljka Smit-McBride, Department of Ophthalmology & Vision Science, Vitreoretinal Research Lab, School of Medicine, University of California, Davis, Davis, CA 95616, One Shields Ave, Tupper Hall #2403, Phone: (530)320-8474, FAX: (530) 752-2270; email: email@example.com
Purpose: To identify changes induced by environmental tobacco smoke (ETS) in circulatory microRNA (miRNA) in plasma and ocular fluids of the Rhesus macaque and compare these changes to normal age-related changes. Tobacco smoke has been identified as the leading environmental risk factor for age-related macular degeneration (AMD).
Methods: All Rhesus macaques were housed at the California National Primate Research Center (CNPRC), University of California, Davis. Four groups of animals were used: Group 1 (1–3 years old), Group 2 (19–28 years old), Group 3 (10–16 years old), and Group 4 (middle aged, 9–14 years old). Group 4 was exposed to smoke for 1 month. Ocular fluids and plasma samples were collected, miRNAs isolated, and expression data obtained using Affymetrix miRNA GeneTitan Array Plates 4.0. Bioinformatics analysis was done on the Affymetrix Expression Console (EC), Transcriptome Analysis Software (TAS) using ANOVA for candidate miRNA selection, followed by Ingenuity Pathway Analysis (IPA).
Results: The expression of circulatory miRNAs showed statistically significant changes with age and ETS. In the plasma samples, 45 miRNAs were strongly upregulated (fold change >±1.5, p<0.05) upon ETS exposure. In the vitreous, three miRNAs were statistically significantly downregulated with ETS, and two of them (miR-6794 and miR-6790) were also statistically significantly downregulated with age. Some retinal layers exhibited a thinning trend measured with optical coherence tomography (OCT) imaging. The pathways activated were IL-17A, VEGF, and recruitment of eosinophils, Th2 lymphocytes, and macrophages.
Conclusions: ETS exposure of Rhesus macaques resulted in statistically significant changes in the expression of the circulatory miRNAs, distinct from those affected by aging. The pathways activated appear to be common for ETS and AMD pathogenesis. These data will be used to develop an animal model of early dry AMD.
Age-related macular degeneration (AMD) is the most common form of visual impairment in elderly individuals in the Western world. AMD is a complex blinding disease, affected by genetic and environmental factors, as well as epigenetic regulators. Cigarette smoke is the leading environmental risk factor, second only to aging. Recently, a central regulatory role for a novel type of epigenetic regulator, miRNA, has been discovered in many human diseases. MiRNAs are considered master regulatory switches, playing roles in numerous health and disease processes, as well as in aging and epigenetics. Some miRNAs can be induced in response to environmental challenges, thus allowing the organism to modify its underlying genetic programming in response to a changing environment.
Although knowledge about miRNAs and their function is accumulating rapidly, knowledge specific to their role in healthy ocular aging and age-related pathologies of the eye is lacking. In this project, we used tobacco smoke exposure as a model of a common environmental insult that affects ocular tissue. Smoke-induced damage is mediated in part through direct oxidation, depletion of antioxidants, complement activation, and inflammation . Interestingly, similar processes have been shown to be critical in the pathology of AMD. Although research in recent years has identified some mechanisms of AMD pathogenesis, this research is limited by the lack of appropriate animal models for AMD. The most commonly used model is the rodent which lacks macula and never truly develops hallmark lesions of AMD, such as drusen. Although primates develop an age-related maculopathy that is equivalent to human AMD , this process takes decades, which is impractical for an experimental approach.
Our long-term goal is to develop a primate model of retinal degeneration, which would resemble AMD. Given that cigarette smoke, or more encompassing, environmental tobacco smoke (ETS) exposure, is the leading environmental risk factor for AMD, second only to aging, we hypothesize that there is a common mechanism between ETS-induced damage and AMD pathogenesis. We hypothesize that miRNAs regulate normal age-related homeostatic processes in the retina. Therefore, the purpose of this study was to identify the initial changes that occur in circulatory miRNA in the ocular fluids and plasma during normal aging and compare the changes to the changes that occur with ETS exposure. We identified differentially expressed miRNAs in plasma and ocular fluids during normal, healthy aging and determined that process was considerably different compared to pathologic changes induced by an environmental insult, modeled by short-term exposure to tobacco smoke. Interestingly, some of the pathways activated with ETS were the same as those involved in AMD pathogenesis.
Rhesus macaques used for this study were housed at the California National Primate Research Center (CNPRC) at University of California, Davis. There were four experimental groups of female animals (a total of 16 animals): Group 1 (young, 1–3 years; n=4), Group 2 (old, 19–28 years; n=5), Group 3 (middle aged, 10–16; n=3), and Group 4 (middle aged, 9–14; n=4). The Group 4 was exposed to smoke for 1 month at the CNPRC facility (Appendix 1). All procedures conformed to the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research. These procedures were also authorized by the Institutional Animal Care and Use Committee (#18988) at the University of California, Davis.
Rhesus macaques selected for the ETS exposure experiment were clinically assessed and found healthy in the absence of (n=6) or following (n=4) ETS exposure. Fundus photography using the CF-1 retinal camera (Canon, Tokyo, Japan) and spectral domain optical coherence tomography (SD-OCT) imaging were performed using the Spectralis device (Heidelberg Engineering, Heidelberg, Germany) as detailed below.
SD-OCT imaging was performed using the Spectralis SD-OCT device (Heidelberg Engineering), which had been modified with a flat chin-rest to allow the Rhesus monkeys’ heads to be positioned. OCT images were captured with a 30° × 5° SD-OCT raster scan with 1536 A-scans per B-scan and 234 μm of spacing between B-scans. This was performed in the high-resolution mode with 25 scans averaged for each B-scan using the device’s eye-tracking capability. Images were semiautomatically segmented using the Duke Optical Coherence Tomography Retinal Analysis Program (DOCTRAP, version 62.0), custom image analysis software designed using MATLAB (Mathworks) . Segmentation boundaries were automatically calculated by DOCTRAP using graph-based determination of reflectivity profiles and then manually refined by the two graders. Thickness measurements were averaged across the 3 mm segment centered on the fovea in each horizontal line-scan, and included the nerve fiber layer (NFL), ganglion cell layer (GCL), inner plexiform layer (IPL), inner nuclear layer (INL), outer plexiform layer (OPL), outer nuclear layer (ONL), photoreceptor inner segments (IS), photoreceptor outer segments (OS), retinal pigment epithelium (RPE), choriocapillaris (CC), and choroid (C). Total retinal thickness was measured from the internal limiting membrane (ILM) to Bruch’s membrane (BM).
All measurements reported in this study were calculated using the average measurements from two graders. The mixed-effect ANOVA was performed to test the pre/post exposure effects on all retinal layers. A mixed-design ANOVA model (also known as a split-plot ANOVA) is used to test for differences between two or more independent groups while subjecting participants to repeated measures. Differences in chorioretinal layer thicknesses obtained from the SD-OCT pre/post exposure were compared using a mixed-effect ANOVA, with two independent variables: The within-subject factor variable was pre/post exposure and eyes [oculus sinister (OS) / oculus dextrus (OD)] while the between-subject factor variable was an animal ID. The dependent variable was the thickness of each layer. All statistical analysis was performed using SPSS software (version 22, IBM). The segmentation analysis and quantification of the retinal layers are shown for the total retinal thickness measurement along the central 3 mm region for all four animals. The mean value of retinal layer thickness was calculated for each layer in four monkeys pre/post exposure. No statistical significance was found for pre/post exposure eyes.
Groups 1 through 3 were housed in standard indoor animal housing, while Group 4 was housed in standard caging in a modified animal holding room specifically designed for ETS. Tobacco smoke exposure occurred over 4 weeks with ETS delivered for approximately 6 h per day, 5 days a week. The exposure room was ventilated with filtered air during non-ETS exposure periods or with ETS diluted with filtered air during exposure periods. Necropsies were performed on the day immediately following the last exposure day.
Carbon monoxide (CO)-- A Teledyne-API Model 300E Carbon Monoxide (CO) Analyzer (Teledyne Instruments, San Diego, CA) was used to continuously monitor the CO concentration in the exposure room for the duration of the exposure. Data were acquired in parts per million (ppm) every 4 min.
Nicotine-- Two samples were collected per day, for 1 h each. Tobacco smoke–exposed air was drawn through a stainless steel sample tube at a flow rate of 2,000 ml/min for 1 h through an XAD-4 sorbent tube, which was sealed after sample collection and stored in the refrigerator. Nicotine samples were analyzed by the UC Davis Center for Health and the Environment (CHE) as previously described .
Total mass concentration-- Tobacco smoke–exposed air was drawn through a stainless steel sample tube at a flow rate of 2,000 ml/min for 1 h. The sample was collected on a Pallflex® 25 mm Fiberfilm filter (Pall Corp., East Hills, NY), which was weighed on a Sartorius Model MC5 microbalance before and after sampling.
Cotinine assays-- Blood samples were collected in 7 ml EDTA vacutainer tubes for the control animals and at 15- and 30-day time points in the ETS-exposed animals. A cotinine ELISA kit was purchased from Calbiotech (Spring Valley, CA), and samples were assayed according to the manufacturer’s protocol.
Ocular fluid samples were collected when the animals were euthanized. Blood samples were collected at the beginning of the experiment and when the animals were euthanized. Animals were euthanized following the standard CNPRC protocol: the animal was sedated with Ketamine at a dosage of 10 mg/kg. After that sodium pentobarbital greater or equal to 120 mg/kg Fatal Plus or Beuthanasia Plus 1 ml/4.5 kg was administered intravenously to the animal for euthanasia. After each animal was euthanized, 100–200 μl samples of aqueous and vitreous humor and 10 ml of venous blood for plasma were collected. Plasma was isolated using ACCUSPIN tubes (MilliporeSigma, St. Louis, MO) following the manufacturer’s protocol. Plasma samples were stored at −80 °C.
Enucleated eyes were immersion fixed overnight at room temperature in 2% phosphate-buffered paraformaldehyde containing 0.5% glutaraldehyde, pH 7.4 (0.1 M phosphate buffer). After overnight fixation, the eyes were bisected, the lens removed, and the bisected eye was placed in 4% phosphate-buffered formalin for approximately 1 week to reverse glutaraldehyde fixation before embedding in paraffin. The eyes were embedded in paraffin as follows: washed in PBS (1X; 154 mM NaCl, 1.1 mM KH2PO4, 5.6 mM Na2HP04, pH 7.4), dehydrated through graded ethanol washes (50%, 75%, 95%, and 100% ethanol), and gradually cleared in xylene using a step-wise procedure (2:1 ethanol to xylene; 1:1 ethanol to xylene; 1:2 ethanol to xylene), followed by xylene before embedding in paraffin. Sections were cut (6 μm) using a Leica RM2125RT microtome (Leica, Nussloch, Germany), placed on SuperFrost Plus microscope slides, and dried overnight at room temperature. Standard hematoxylin-eosin staining was performed on the sectioned eyes.
All the samples were collected within a 1 h after the animals were euthanized, processed immediately, and stored in aliquots at −80 °C until RNA isolation. The total RNA, which included miRNAs, was isolated from each of the individual samples using miRCURY™ RNA Isolation Kits -Biofluids (Qiagen, Germantown, MD). A small aliquot of the sample was frozen separately for quantification and quality control on BioAnalyzer, so the main sample tube was thawed only once during the process of labeling probe for GeneChip microarray hybridization. Total isolated RNA (5–10 ng/μl) from fluid samples was analyzed using the BioAnalyzer using Small RNA Analysis Kit (Agilent Technologies, Santa Clara, CA) following the manufacturer’s protocol.
MiRNA samples (20 ng) were labeled as probes, using the FlashTag kit (ThermoFisher Scientific, Waltham, MA) following the manufacturer’s protocol and hybridized to microRNA microarrays. Two Affymetrix miRNA microarrays GeneTitan Plates 4.0, each with 24 samples, for a total of 48 samples, were run at the UCSF Institute for Human Genetics (IHG). Raw data in the form of probe cell intensity data (.cel) files were retrieved and subjected to bioinformatics analysis.
Probe cell intensity data (.cel) from Affymetrix GeneChip® miRNA Arrays were analyzed with the Affymetrix® Expression Console (EC) software (ThermoFisher Scientific). The application uses the RMA + DABG analysis for doing background subtraction, normalization, and summarization of the probe sets from Affymetrix expression microarrays, and the final output is a .chp file. Spearman rank order correlation coefficients between biologic replicates were computed by EC to assess the between samples correlations (Appendix 2). Generated .chp files for each sample were imported into Affymetrix Transcriptome Analysis Console (TAC; ThermoFisher Scientific) for subsequent analysis of differentially expressed genes.
Statistical analysis of microarray data-- The criteria used to determine a statistically significant change in miRNA abundance was −1.5 ≥ fold change (FC) ≥1.5, p≤0.05. ANOVA was used for calculating the statistical significance and identifying the p values for each microRNA. The data were filtered for Macaca mulatta (Rhesus macaque).
Ingenuity pathway analysis-- Ingenuity Pathway Analysis software (Qiagen) was used for gene pathway analysis and the identification of potential miRNA gene targets. Raw data were deposited at NCBI Gene Expression Omnibus (GEO), Accession Number GSE108951.
Real-time qPCR assays (TaqMan assays, ThermoFisher Scientific) were run on miRNA isolated from plasma samples for ETS versus control on three replicas for each group and in triplicate for each replica. Assays were run at the Real-time PCR Research and Diagnostics Core Facility at UC Davis. Each miRNA sample was amplified using a TaqMan preamplification procedure and TaqMan preamp primers (Applied Biosystems, Life Technologies, Grand Island, NY) following manufacturers protocols. Reverse transcription thermal-cycling conditions were: 16 °C/30 min, 42 °C /30 min, 85 °C /5 min, hold 4 °C /∞. Pre-amplification thermal-cycling conditions were: 95 °C /10 min, 55 °C /2 min, 72 °C /2 min, 95 °C /15 s, 12 cycles of 95 °C /15 s, 60 °C /4 min, hold at 99.9 °C /10 min, hold 4 °C /∞. The samples were run in 384 well plates. Real-time thermal-cycling conditions: hold 95 °C /10 min, 40 cycles of 95 °C /15 sec, 60 °C /60 s, hold 4°C /∞. The TaqMan assays performed were miR-221 (ABI No. 000524, hsa-miR-221), miR-26a (ABI No. 000405, hsa-miR-26a), let-7c-5p (ABI No. 000379, hsa-let-7c), and three control genes, miR-486-5p (ABI No. 001278, hsa-miR-486), miR-320b (ABI No. 002844, hsa-miR-320b), and miR-638 (ABI No. 001582, hsa-miR-638). MiRNA assays were conducted at a dilution of 1:100s, and each run was done in triplicate. A control with no enzyme in the reverse transcriptase (RT) step was included as a negative control.
The adult female Rhesus macaques in Group 4 (exposure) received four 7-day cycles of 1.42–2.15 mg/m3 of tobacco smoke for an average of 6 h/day on days 1–5, followed by 2 days in the filtered air (Table 1). During exposure, filtered air ventilation was lowered to 3.5 room air exchanges per hour. The daily average nicotine level of exposure measured through XAD-4 sorbent tubes was (0.610 mg/m3 ±0.090). In animals that underwent ETS exposure, ocular health was evaluated by veterinary and physician ophthalmologists (SMT and GY) pre- and post-exposure. The control group (Group 3) was housed in standard indoor animal housing.
Cotinine is an alkaloid found in tobacco and is the predominant metabolite of nicotine . Cotinine concentrations in the control and ETS-exposed animals at day 30 are displayed in Table 2. On average, ETS-exposed animals contained statistically significantly higher (p<0.01), concentrations of cotinine compared to control animals at the 15-day time point (48.6±13.1 ng/ml versus 0 ng/ml) and the 30-day time point (54.9±5.0 ng/ml versus 0.06±.001 ng/ml).
Rhesus macaques from Groups 3 and 4 were clinically assessed before and following ETS exposure. Ophthalmic examinations did not identify any visible ocular changes pre- and post-exposure (Figure 1). Retinal imaging, performed using OCT, followed by segmentation analysis and quantification of the retinal layers is shown for the total retinal thickness measurement along the central 3 mm region for all four animals (Table 3). There were trends noted in the thinning of the nuclear fiber layer (NFL; 17.36±0.860 um versus 16.01±0.470 um (p<0.17) and in the thickening of the choriocapillaris layer (CC; 18.60±2.050 versus 19.46±2.240 (p<0.34), pre- and post-ETS exposure, although these trends were not statistically significant.
Representative profiles from the BioAnalyzer small RNA kit are shown in Figure 2. There was an abundance of circulatory miRNA molecules in each of the fluids analyzed. Representation of miRNA in each of the fluid compartments isolated from 200 μl of fluid was within the range of 2.30–14.6 ng/μl in the aqueous, 0.8–2.2 ng/μl in the vitreous, and 0.16–0.32 ng/μl in the plasma. The small RNA bands were abundant (Figure 2, red arrow) and the range of the sizes of the small RNA molecules was comparable to the range in the profile of the high-quality small RNA fraction in the solid tissue RNA prep (data not shown). The total concentration of miRNA in nanograms/microliters did not change statistically significantly in the plasma samples from the same animals pre- and post-exposure.
In the plasma samples, 40 microRNAs were ubiquitously present (p<0.05) in all samples, and up to 204 variably. In vitreous, ten miRNAs were ubiquitously present in all samples and up to 88 variably, and in aqueous, seven were ubiquitously present and up to 148 variably. The presence of microRNAs was established according to the Affymetrix P/A/M detection algorithm based on the Wilcox signed-rank test introduced by Affymetrix .
In the plasma samples, 45 miRNAs were statistically significantly upregulated after ETS exposure (Figure 3A). In the vitreous samples, three miRNAs were statistically significantly downregulated after ETS exposure (Figure 3B). There were no statistically significantly dysregulated miRNAs observed in the aqueous samples. All miRNAs with their fold change and p-values (−1.5 ≥ FC ≥ 1.5; p<0.05) are listed in Appendix 3, Appendix 4, Appendix 5, Appendix 6, Appendix 7.
In the plasma samples from the aging monkeys, two miRNAs were statistically significantly downregulated between old (Group 2) versus young (Group 1) macaques (Figure 4). In the vitreous and aqueous samples, six miRNAs were statistically significantly downregulated in each compartment (−1.5 ≥ FC ≥ 1.5; p<0.05). All dysregulated miRNAs differed between the vitreous and aqueous samples.
Independent quantitative PCR (qPCR) validation of representative differentially expressed genes was performed using commercial gene expression assays (TaqMan; Applied Biosystems), following manufacturers protocol (for detailed reaction conditions see Materials and Methods section on Quantitative polymerase chain reaction (qPCR) assays). The fold changes and their direction were confirmed for the chosen sets of genes. Data are presented in Table 4.
Ingenuity Pathway Analysis identified two top networks that are activated during ETS exposure: IL-17A and VEGF. One of the top IPA pathways that incorporates the activation of both networks is the inflammatory pathway of recruitment and activation of eosinophils, Th2 lymphocytes, and macrophages (Figure 5).
In addition to a role in cardiovascular and cerebrovascular diseases, as well as a host of cancers, tobacco smoking is the principal risk factor for AMD in the industrialized world. The genetics of AMD has been well studied with several risk alleles identified, such as those found in CFH (Gene ID: 3075, OMIM 134370) and HTRA1 (Gene ID: 5654, OMIM 610149) [7,8]. However, environmental factors, which are easily modifiable and potentially could prevent or mitigate AMD, are understudied. Thus, there is an increasing urgency to understand the mechanisms involved in ETS pathogenesis, particularly given the proliferation of new forms of nicotine delivery [9-11].
Curtis’s and other groups have shown that 1 month of high ETS exposure produces chronic obstructive pulmonary disease (COPD) in primate animal models [12,13]. Our working hypothesis was that the retina and the choroid, highly vascular tissues that are sensitive to oxidative stress, might be similarly impacted by short-term ETS exposure. Furthermore, environmental second-hand smoke has been shown to have serious health effects [14-16]. We measured animal exposure to ETS using several criteria: particulate matter, nicotine, and cotinine. Cotinine is an alkaloid found in tobacco and the predominant metabolite of nicotine . Cotinine has an in vivo half-life of approximately 20 h and is typically detectable up to 1 week after the use of tobacco. The amounts of cotinine in blood, saliva, and urine are proportionate to the amount of exposure to tobacco smoke; thus, cotinine is a valuable indicator of tobacco smoke exposure, including secondary (passive) smoke . The nicotine concentrations measured in the current study were consistent with second-hand smoke exposure, and the cotinine assays confirmed marked ETS exposure with a concentration of about 50 ng/ml. In humans, cotinine concentrations of about 30 ng/ml indicate marked exposure .
All macaques were clinically assessed before and following ETS exposure. Fundus photography of the back of the eye did not show any visible changes pre- or post-exposure. Retinal imaging using OCT, followed by segmentation analysis and quantification of the retinal layers, showed no statistical significance in the total retinal thickness or any of the retinal layers. This result was consistent with findings in human smokers, where the central macular thickness shows no statistically significant differences between smokers and non-smokers [19,20]. Although not reaching a statistical significance level, perhaps because of the small animal numbers and only 1 month of exposure, a trend in the thinning of the NFL was observed. This observation is consistent with the literature showing that smoking tobacco may cause considerable thinning of the NFL in heavy smokers [20-22]. In addition, we observed the trend of choriocapillaris thickening in primates exposed to ETS. This result is consistent with the findings that exposure to ETS results in the formation of sub-RPE deposits, thickening of Bruch’s membrane, and the accumulation of deposits within Bruch’s membrane. All these processes lead to an injury stimulus of the choriocapillaris and RPE, which may explain the association between tobacco smoking and early AMD . One of the limitations of this study is that the OCT-based measurements are limited by the inclusion of Henle’s fiber layer within the ONL rather than the OPL to which the layer belongs . Thus, true longitudinal changes in the OPL at the macula may not be accurate.
Upon ETS exposure, the greatest changes in circulatory miRNAs were observed in the plasma (45 miRNAs were strongly upregulated), followed by the vitreous (three miRNAs were downregulated). However, there were no statistically significantly dysregulated miRNAs detected in the aqueous humor. These results differ markedly from what occurs during normal aging where only two miRNAs were statistically significantly downregulated in the plasma, while in the vitreous and aqueous six different miRNAs were statistically significantly downregulated (-1.5>FC>1.5; p<0.05). The only exceptions, two out of three novel miRNAs identified as downregulated following ETS exposure in the vitreous (miR-6794 and miR-6790) were the same ones observed to be downregulated in the vitreous during normal aging.
Circulating miRNAs are excellent biomarkers. They can be easily sampled from plasma and bodily fluids, including ocular fluids, thus gaining the name of the liquid biomarker or liquid biopsy . They are stable for hours, even at room temperature, because they travel protected associated with proteins or encapsulated by small vesicles called exosomes .
In addition to being accessible biomarkers, miRNAs are biologically powerful molecules that regulate post-transcriptional gene expression. MiRNA expression is necessary for cellular function, such as cell proliferation, apoptosis, inflammation, and a wide variety of cellular processes. MiRNAs function through base-pairing with mRNA, which leads to degradation or impaired translation of that particular mRNA . Therefore, they can be referred to as “dimmer switches” for gene expression, because of their ability to repress gene expression without completely silencing it. Six of the miRNAs identified in the present study in primates exposed to ETS (let-7a, let-7c, miR-106a-5p, miR-106b-5p, miR-17-5p, and miR-361-5p) were also identified as circulatory candidate miRNAs in AMD [28-31]. The let-7 family and miR-361-5p are upregulated in retinal tissue in advanced AMD, targeting TGFBR1 [29,32]. MiR-17-5p has been shown to be upregulated in the plasma of patients with wet AMD . MiR-17-5p negatively regulates TGFBR2 expression by directly binding to the 3′ untranslated region (UTR) of TGFBR2 mRNA, thus promoting cell growth and migration in cancer . MiR-106 has been identified in a cohort of miRNAs that respond to oxidative stress exposure in ARPE-19 cells . MiR-106a has been shown to target PTEN mRNA and lncRNA LINC00657, thus playing an important role in the regulation of PTEN expression in cancer .
Ingenuity Pathway Analysis of candidate miRNAs and their target genes identified two top networks that were activated during ETS exposure: IL-17A and VEGF. Upregulation of IL-17A has been previously observed with ETS exposure. Specifically, Siew and colleagues demonstrated that ETS increases bronchial mucosal IL-17A expression, and neutrophilic inflammation in asthmatics . Levanen and colleagues implied that long-term exposure to tobacco smoke increases extracellular cytokine signaling via IL-17A in the peripheral airways of smokers with COPD and that this signaling may involve cytotoxic T-cells . Interestingly, IL-17A also appears to be important in the pathogenesis of AMD. For example, genetic variants of IL-17A have been functionally associated with an increased risk of AMD . Zhang and colleagues suggested that IL-17A triggers a key inflammatory mediator, IL-1β, in RPE cells, via NLRP3 inflammasome activation . Therefore, the IL-17A pathway may be an important therapeutic target for AMD. Another study demonstrated that a functional single nucleotide polymorphism in the IL-17A 3′ UTR is targeted by miRNA in vitro and that miRNA-mediated gene dysregulation may play a role in the pathogenesis of AMD .
ETS exposure causes extensive tissue remodeling in the upper airways, and VEGF plays a critical role in this process. It has been shown that ETS extract has a stimulatory effect on VEGF expression through the TLR4/ROS/MAPKs/NF-kappaB signaling pathway in nasal fibroblasts . In addition, a human RPE cell line (ARPE-19) exposed to an ETS concentrate resulted in premature senescence, thus supporting the role of oxidative damage in ETS-induced senescence activation. Furthermore, the senescent ARPE-19 cells upregulated VEGF and simultaneously downregulated complement factor H expression. As these phenomena are involved in AMD pathogenesis, these studies coupled with the present results support the hypothesis that ETS-induced senescence activation is important in the induction and progression of AMD. Moreover, this hypothesis would also explain the striking association of AMD pathogenesis with tobacco smoking .
One of the top IPA pathways that incorporates activation of both networks is the inflammatory recruitment and activation of eosinophils, Th2 lymphocytes, and macrophages pathway. This result fits with studies performed on infant monkeys after perinatal ETS exposure, which showed a statistically significant increase in mast cell, eosinophil, monocyte, and lymphocyte cell numbers in the lungs of infant monkeys . It has been demonstrated that ETS exposure during early life enhances local Th2 immunity by impairing normal pulmonary Th1 immune maturation. This effect was greater in animals that began ETS exposure in utero compared to those that were exposed postnatal . However, an association of the inflammatory response and two distinct types of immune responses—innate (complement system, phagocytes, mast cells, eosinophils, macrophages) and adaptive (T- and B-lymphocytes)—have been strongly implicated in AMD. An overactive complement system is associated with AMD pathogenesis, and serum proinflammatory cytokines, including IL-17, are elevated in patients with AMD. IL-17 is produced by complement C5a-receptor-expressing T-cells, such as gamma-delta T-cells, members of the innate immune response [46-48]. VEGF plays an important role in pathological angiogenesis, and anti-VEGF agents have been used to treat ocular diseases that are driven by pathological angiogenesis [49-51]. Based on the present results from the ETS-exposed monkeys and literature published on AMD, it appears likely that activation of IL-17A and VEGF networks, as well as inflammatory pathway of recruitment and activation of eosinophils, Th2 lymphocytes, and macrophages, are common for tobacco smoke and AMD pathogenesis.
In conclusion, this study identified changes in circulatory miRNAs in the plasma and ocular fluids of the Rhesus macaque following ETS exposure, which differed markedly in comparison to age-related changes. The results of this study suggest that ETS exposure may be an appropriate animal model for early dry AMD. The identification and characterization of eye-specific miRNA in ocular tissue, as well as their target molecular pathways, will generate opportunities to identify novel strategies for therapeutic and environmental modulation or intervention for blinding ocular diseases.
Appendix 1. Supplementary table 1. Rhesus macaques in the study.
Appendix 2. Supplementary figure 1
Appendix 3. Supplementary table 2
Appendix 4. Supplementary table 3
Appendix 5. Supplementary table 4
Appendix 6. Supplementary table 5
Appendix 7. Supplementary table 6
The authors are grateful to the CNPRC, the UC Davis Inhalation Exposure Core operators Louise Olsen, Chris Royer and Ross Allen, Dale Uyeminami from Environmental Health Center for nicotine exposure calculations, Kevin NH Nguyen for assisting with OCT scan readings, and Dr. Philip H. Kass, Professor Population Health & Reproduction, UC Davis, for mentoring ZW in statistical analysis. This project was funded by Barr Family Retina Research Foundation generous gift to the UC Davis Department of Ophthalmology (LSM, ZSM) and a pilot grant from the Environmental Health Center National Institute of Environmental Health Sciences P30ES023513 (ZSM) and NIH P30EY005722 (SF). The NEI Eye Core grant P30 EY012576 provided support for the tissue fixation, paraffin embedding and H&E staining. The CNPRC NIH P51OD011107 grant provided support for the colony and infrastructure, including the cores. The Affymetrix microRNA TitanChips 4.0 were processed through genomics services offered at the Institute for Human Genetics at UCSF.