Direct Effects of Glyphosate on In Vitro T Helper Cell Differentiation and Cytokine Production

Glyphosate (G) is the active ingredient of the most used herbicides worldwide. Its use is currently very debated, as several studies indicating its hazard and toxicity are emerging. Among them, there is evidence of adverse effects on the immune system. The aim of this work was to investigate if G could directly affect immune cells. Peripheral blood mononuclear cells (PBMC) obtained from healthy donors were used as experimental model. PBMC were expose to G and stimulated with PMA/ionomycin, T helper (Th) cell differentiation and cytokine production were assessed by flow cytometry and enzyme-linked immunosorbent assay, respectively. A reduction of Th1/Th2 ratio, mainly due to a decrease in Th1 cells, was observed following G exposure. Results show an enhancement of IL-4 and IL-17A production, and a reduction of IFN-γ. Based on literature evidence that suggest G being an endocrine disruptor, we investigated the role of nuclear estrogen receptors (ER). ERα/ERβ inhibition by ICI 182,780 abolished the effects of G on IFN-γ and IL-4 release, suggesting a role of ER in the observed effects. To further characterize the mechanism of action of G, miRNAs, both in exosome and intracellular, were investigated. A statistically significant increase in miR-500a-5p was observed following G treatment. The blockage of miR-500a-5p, using a specific antagomir, prevented G-induced reduction of IFN-γ production. Finally a relationship between miR-500a-5p up-regulation and ER was observed. Overall, these results suggest that G can directly act on T cells, altering T cell differentiation and cytokines production.

Introduction

Glyphosate (N-(phosphonomethyl) glycine; G) is the active ingredient of the most used broad-spectrum herbicide worldwide (1, 2). It was commercialized starting from 1974 under the name of Roundup^®, together with adjuvants, like polyoxyethylene tallowamine, that, although considered inert compounds, can increase G toxicity (3).

Its success further increased following the introduction of genetically modified G-resistant crops (4). The herbicidal activity is based on the competitive inhibition of the shikimate pathway, which leads to the blockage of the synthesis of aromatic amino acids in plants (5). Since animals do not have this pathway, G has always been considered safe for humans by regulatory agencies (6, 7).

However, evidence of G toxicity and its alleged threats to the ecosystem and human health are emerging (8, 9). Humans can be exposed to G through different routes; mainly orally, but also through dermal and inhalation exposures (10, 11). Although the majority of studies addressing G toxicity used high doses, there are also studies using environmental relevant doses showing noxious effects in both animals and humans (12–16).

As for general toxicity, there is also evidence of adverse effects against the immune system. Indeed, the latter can be the target of several pesticides (17–19). G-induced immunotoxicity could be the result of a direct action of the herbicide, or can be due to an indirect mechanism, secondary to endocrine, nervous system, or microbiota alterations (20–22). Regarding mammals, evidence indicates the ability of G and G-based herbicides to alter immune responses and to induce inflammation (22). More specifically, intranasal injection of G-rich air samples in mice enhanced eosinophils, neutrophils, and mast cells degranulation in lungs, suggesting an asthma-like pathology (23). In this study, an activation of T helper (Th) 2 response, as reported by the higher production of Th2 type cytokines in lungs, namely IL-4, IL-5, IL-10, IL-13 and IL-33, was observed. Roundup^® was shown to induce pro-inflammatory cytokines (i.e., IL-1β and TNF-α) in liver of rats exposed via feeding at the doses of 100 and 250 mg/kg bw/day (24). Similarly, the exposure via gavage was able to increase the mRNA levels of IL-1β, IL-6, TNF-α, MAPK3, NF-κB, and caspase-3 in the jejunum of rats exposed to the doses of 50 and 500 mg/kg bw/day (25). Authors hypothesized a role of microbiota in G-induced immunotoxicity, as a significant decrease in the relative abundance of Firmicutes and Lactobacillus, while an enrichment in several potentially pathogenic bacteria were observed (25). The same result was obtained by Qiu etal. (26), which treating piglets with Roundup^® in diet, observed increased mRNAs of IL-6 and NF-κB in the jejunum following the exposure to 10, 20 and 40 mg/kg bw/day.

Beside animal studies, there are also few epidemiological studies supporting G-based herbicides affecting human immune system. The Agricultural Health Study demonstrated an association between G occupational exposure and current rhinitis and increased rhinitis episodes (27, 28). This study revealed also an association between allergic and non-allergic wheezes among G sprayers (29).

As microbiota is known to influence the immune system and considering that the shikimate pathway is also found in some microorganisms, the interaction of G with commensal microorganisms represents a plausible explanation of G-induced immunotoxicity. Indeed, recently G-based herbicide ability to inhibit the shikimate pathway in the microbiome was demonstrated (30). Furthermore, Roundup LP Plus^® was demonstrated to alter the metabolism of E. coli, that in turn impacted on the inflammatory immune response, inducing a higher production of TNF-α (31). In addition, associations between G and immune-endocrine disturbances have been described, which also represent a plausible mechanism underlaying G immunotoxicity (reviewed in 22).

The aim of this study was, therefore, to investigate if G could directly affect T helper (Th) cell differentiation and functions. The study focused on Th cells as they have been identified in mice as a possible target of G (23). Investigations were also conducted to get insights on the underlying mechanism of action focusing on the role of the estrogen receptor (ER) and miRNAs in G-induced immunotoxicity. Results obtained demonstrated that G can directly act on T cells, altering T cell differentiation and cytokines production, with a role of estrogen receptor and miR-500-5p in the observed effects.

Materials and methods

Chemicals

G (CAS #1071-83-6, purity ≥ 99%) was purchased from Sigma-Aldrich (St. Louis, MO, USA). It was dissolved in Dulbecco’s Phosphate Buffered Saline. Phorbol 12-myristate 13-acetate (PMA – CAS #16561-29-8, purity ≥ 99%) and ionomycin from Streptomyces conglobatus (CAS #56092-81-0, purity ≥ 98%) were purchased from Sigma-Aldrich and used to stimulate immune cells, as well as phytohaemagglutinin (PHA, CAS #9008-97-3). 17β-estradiol (E2 – CAS #50-28-2, purity 98%) was used as reference compound for ER activity, and it was obtained from Sigma-Aldrich. ICI 182,780 (CAS # 129453-61-8, purity ≥ 99%) was used as ER antagonist and it was obtained from Tocris Bioscience (Bristol, United Kingdom). All the substances, with the exception of G, were dissolved in dimethyl sulfoxide (DMSO, CAS # 67-68-5, purity ≥ 99.5%), with a final DMSO concentration in culture medium ≤ 0.2%. Cell culture medium and all supplements were also purchased from Sigma-Aldrich.

Cells

Peripheral blood mononuclear cells (PBMC) were obtained by Ficoll gradient centrifugation from anonymous buffy coats of 5 male healthy donors, purchased from the Niguarda Hospital in Milan (Italy). After centrifugation, PBMC layers were removed and washed 5 times with Dulbecco’s Phosphate Buffered Saline. Isolated cells were diluted to 10⁶ cells/mL in RPMI 1640 with phenol red containing 2 mM L-glutamine, 0.1 mg/mL streptomycin, 100 IU/mL penicillin, gentamycin 10 µg/mL, 50 µM 2-mercaptoethanol, supplemented with 10% heated-inactivated fetal bovine serum (culture media) and cultured at 37°C in 5% CO₂ incubator. Experiments with E2 and ICI 182,780 were performed using RPMI 1640 without phenol red containing 2 mM L-glutamine, 0.1 mg/mL streptomycin, 100 IU/mL penicillin, gentamycin 10 µg/mL, 50 µM 2-mercaptoethanol, supplemented with 5% heat-inactivated dialyzed fetal bovine serum. Preliminary experiments were conducted to identify non-cytotoxic concentrations (cell viability > 90%). Cytotoxicity was assessed by LDH Cytotoxicity Detection Kit (Takara Bio USA, Inc.) (data not shown).

Th Cells Analysis

Phenotyping of Th cells was conducted using the Human Th1/Th2/Th17 Phenotyping Kit (BD Pharmingen, Inc.). PBMC (10⁶ cells/mL) were incubated overnight (o.n.) and then treated in the presence or absence of G at the concentrations of 0.01, 0.1, 1 and 10 µg/mL. After 1 hour, cells were stimulated with PMA (50 ng/mL), ionomycin (1 µg/mL) and GolgiStop™ Protein Transport Inhibitor for 5 hours, following manufacturer’s instructions. After treatment, cells were centrifuged at 1300 rpm for 5 minutes and stained with anti-human CD4 PerCP-Cy5.5 conjugated, anti-IL-17A-PE, anti-IL-4-APC and anti-IFN-γ-FITC following manufacturer’s instructions. Cells were analyzed using NovoCyte 3000 flow cytometer, and data were quantified using NovoExpress software (Acea Biosciences, Inc). Results are expressed as fold-change of G-treated versus control cells.

Cytokine Production

PBMC (10⁶ cells/mL) were incubated o.n. and then treated with G at the concentrations of 0.01, 0.1, 1 and 10 µg/mL. After 1 hour, cells were stimulated with PMA (10 ng/mL) and ionomycin (100 ng/mL) or PHA (1.2 μg/mL) for 72 hours. Following treatment, cells were centrifuged at 1300 rpm for 5 minutes and supernatants collected for cytokine measurement and stored at -20°C until measurement. Cytokine production was assessed in cell-free supernatants by specific sandwich enzyme-linked immunosorbent assays (ELISA) commercially available. ELISA for IL-4 and IFN-γ were purchased from ImmunoTools (Friesoythe, Germany), while IL-17A ELISA from BioLegend (San Diego, CA, USA). Antibodies dilutions were performed according to manufacturer’s instructions. Limits of detection were 2.3 pg/mL for IL-4, 24 pg/mL for IFN-γ, and 3.9 pg/mL for IL-17A. Results are expressed as fold-change of released cytokines of G-treated versus control cells.

miRNAs Contained in Exosomes Analysis

PBMC (3 x 10⁶) were incubated o.n. and then treated with G at the concentration of 0.1 µg/mL for 72 hours in culture medium without fetal bovine serum. Following treatment, PBMC were centrifuged at 260g at 25°C for 5 minutes and supernatants used for exosome extraction. Cell-free supernatants were centrifuged at 20’000g at 4°C for 30 minutes, and supernatants centrifuged again at 100’000g at 4°C for 90 minutes. The pellet obtained represented the exosomes released by PBMC. Total RNA was extracted using RNeasy Plus Mini Kit (Qiagen, Valencia, CA, USA) following manufacturer’s instructions. Retro-transcription was performed according to manufacturer’s instructions (miScript II RT Kit, Qiagen). Before Real-Time PCR, the selected miRNAs were pre-amplified using the miScript PreAMP PCR kit (Qiagen). For PCR analysis miScript SYBR^® Green PCR Kit (Qiagen) was used. All the primers were purchased from Qiagen and used according to manufacturer’s instructions. The quantification of the miRNAs was performed by the 2^−ΔΔCT method. The fold-changes of 7 miRNAs (hsa_miR-10b, hsa_miR-27a, hsa_miR-100, hsa_miR-136, hsa-miR-424, hsa_miR-500a, hsa_let-7f) were analyzed, and miRNA hsa-RNU6-2 was used for normalization.

Intracellular miR-500a Analysis

Following the above-mentioned treatment and centrifugation at 260g at 25°C for 5 minutes, total RNA was extracted from the pellet, using RNeasy Plus Mini Kit. RNA was retro-transcribed and the fold-change of hsa_miR-500a was analyzed using miScript SYBR^® Green PCR Kit. The quantification was performed by the 2^−ΔΔCT method, using miRNA hsa-RNU6-2 for normalization.

miRNA Silencing

Hsa_miR-500a-5p was silenced through the use of miRCURY LNA™ miRNA inhibitor (Qiagen). 400’000 cells, after o.n. incubation, were silenced with 25 pmol of miRNA inhibitor and 3 µl of HiPerFect Transfection Reagent (Qiagen) for 24 hours. Then, the medium was replaced, and cells treated in the presence or absence of G (0.1 µg/mL) for 1 hour, and then stimulated with PMA and ionomycin for 72 hours. Following treatment, cells were centrifuged at 1300 rpm for 5 minutes and supernatants collected for IFN-γ measurement.

Data Analysis

With exception of PHA-induced cytokine release (n=3), all experiments were conducted using 5 donors. Data are expressed as mean ± standard error (SEM). Statistical analysis was performed using GraphPad Prism version 9.1.1 (GraphPad Software, La Jolla, CA, USA). Data were analyzed by One-way analysis of variance (ANOVA) followed by Dunnett’s multiple comparison test or Tukey’s multiple comparison test or by paired Student’s t test. Differences were considered significant at p ≤ 0.05.

Results

Effects of G on Th Differentiation

Epidemiological studies (27, 29) and mice study (23) suggest that G may be associated with increased risk of allergic reactions with increased Th2 responses. To assess if G was able to directly affect immune cells, its effects on Th differentiation were evaluated in vitro. PBMC were treated with G alone at increasing concentrations for 1 hour and stimulated with PMA and ionomycin for subsequent 5 hours. IFN-γ, IL-4 and IL-17A positive cells were analyzed by flow cytometry (Figure 1). Gating strategies for flow cytometry panels are reported in Supplementary Material (Supplementary Figure 1). A reduction of IFN-γ positive cells was observed following G treatment that reached statistically significance at 1 µg/mL (Figure 1A). No significant changes were observed in IL-4 positive cells (Figure 1B). When analyzing the ratio between IFN-γ and IL-4 positive cells (Figure 1C), a decrease was observed, reaching the statistical significance at the concentration of 10 µg/mL. No statistically significant effects were observed in IL-17A positive cells (Figure 1D). Overall, data indicate that G can directly affect Th cell differentiation, resulting in an imbalance of Th1/Th2 subpopulations, supporting an increase in Th2 responses.

FIGURE 1

Figure 1 Effect of G on PMA plus ionomycin-induced CD4⁺ cell differentiation. PBMC (10⁶/mL) were treated for 1 h with increasing concentrations of G (0.01, 0.1, 1 and 10 µg/mL), and then stimulated with PMA, ionomycin and GolgiStop™ for 5 h, as described in the Materials and Methods section. (A) CD4⁺ IFN-γ⁺ cells. (C) ⁺ IL-4⁺ cells. (E) CD4⁺ IL-17A⁺ cells. (G) CD4⁺ IFN-γ⁺/CD4⁺ IL-4⁺ cell ratio. Results are expressed as fold-change of the number of positive cells in G treated cells compared to control cells. The dotted line reported is set at 1.0 (control). Each value represents the mean ± SEM, n = 5 donors. Each dot represents the value of the single individual. Statistical analysis was performed with Dunnett’s multiple comparison test, with *p < 0.05 vs control cells. Flow cytometric analysis of representative histogram overlay were reported in (B, D, F).

Effects of G on Cytokine Production

Next, the effect of G on cytokine release, namely IFN-γ, IL-4, IL-17A, in response to PMA + ionomycin or PHA was investigated. PBMC were treated with increasing concentrations of G alone for 1 hour and then stimulated with PMA and ionomycin or PHA for 72 hours. Results were overall consistent with data shown in Figure 1. G reduced IFN-γ release (Figure 2A) which resulted statistically significant at 0.1 and 1 µg/mL. The release of IL-4 (Figure 2B) was enhanced, although not linearly, reaching statistically significance at the concentrations of 0.01 and 0.1 µg/mL. This imbalance in IFN-γ and IL-4 release resulted in a reduction in the IFN-γ/IL-4 ratio (Figure 2C), indicating an action in favor of Th2 cells and against Th1 cells. Regarding IL-17A production (Figure 2D), a statistically significant increase was observed at the concentrations of 0.1 and 10 µg/mL, with no clear dose-response. Similar results were also observed using PHA (Table 1), a lectin known to bind to T-cell membranes stimulating activation and proliferation (32), confirming the ability of G to unbalance Th1/Th2 responses, favoring Th2.

FIGURE 2

Figure 2 Effect of G on PMA plus ionomycin-induced IFN-γ, IL-4 and IL-17A release. PBMC (10⁶/mL) were treated for 1 h with increasing concentrations of G (0.01, 0.1, 1 and 10 µg/mL), and then stimulated with PMA and ionomycin for 72 h. Cytokines were measured by ELISA in cell-free supernatants. (A) IFN-γ release. (B) IL-4 release. (C) IL-17A release. (D) IFN-γ/IL-4 ratio. Results are expressed as fold-change of the released cytokines in G treated cells compared to control cells. The dotted line reported is set at 1.0 (control). Each value represents the mean ± SEM, n = 5 donors. Each dot represents the value of the single individual. Statistical analysis was performed with Dunnett’s multiple comparison test, with *p < 0.05 and **p < 0.01 vs control cells.

TABLE 1

Table 1 Effect of G on PHA-induced cytokine production.

Our findings demonstrated a non-monotonic response of G on some endpoints of Th differentiation and cytokine production. For instance, lower concentrations can exert higher effects, with a non-linear trend.

Based on these results, subsequent experiments were conducted using G at the concentration of 0.1 µg/mL. This concentration is biologically relevant as serum level of 0.1891 μg/mL were detected in Thai women involved in agricultural activities (33).

Role of ER in G-Induced Effects

Even if debated, some literature data suggest that G may act through the estrogen receptor (ER) (34). To investigate a possible role of ER in G-induced immunotoxicity the ER inhibitor ICI 182,780 was used (35). E2 was used as positive control. The release of IFN-γ, IL-4 and their ratio was investigated, following treatment with G at the concentration of 0.1 µg/mL or E2 (10 ng/ml), with a pre-treatment of 15 minutes in the presence or absence of ICI (1 µM). G and E2 similarly reduced the production of PMA-ionomycin induced IFN-γ (Figure 3A). ICI pre-treatment combined with G or E2 was able to completely restore IFN-γ production. The same behavior was observed with the release of IL-4 (Figure 3B), which was enhanced by both G and E2 and restored by ICI. The decrease in IFN-γ/IL-4 ratio induced by both G and E2, was prevented by ICI, indicating a role of ER in G-induced Th1/Th2 imbalance.

FIGURE 3

Figure 3 Role of ERα in G-induced cytokine release. PBMC (10⁶/mL) were treated with ICI 182,780 (1 µM) for 15 minutes, and then with G (0.1 µg/mL) or with E2 (10 nM) for 1 h. After that, cells were stimulated with PMA plus ionomycin for 72 h. Cytokines were measured by ELISA in cell-free supernatants. (A) IFN-γ release. (B) IL-4 release. (C) IFN-γ/IL-4 ratio. Results are expressed as fold-change of the cytokine released in G treated cells versus relative control cells (untreated cells for G- or E2-treated cells and cells treated with ICI for G+ICI- or E2+ICI-treated cells). The dotted line reported is set at 1.0 (control). Each column represents the mean ± SEM, n = 5 donors. Each dot represents the value of the single donor. Statistical analysis was performed with Tukey’s multiple comparison test, with *p < 0.05 and **p < 0.01 vs control untreated cells or the corresponding not ICI-treated cells.

Role of miR-500a in G-Induced Effects on Lymphocytes

To further investigate the mechanism underlying G immunotoxicity, starting from a miRNA panel conducted on exosomes obtained from plasma samples of farmers occupationally exposed to G-based herbicide (manuscript in preparation), seven de-regulated miRNAs were selected, namely miR-10b, miR-27a, miR-100, miR-136, miR-424, miR-500a, let-7f. PBMC were treated with G (0.1 µg/mL) for 72 hours, and the presence of the selected miRNAs in exosomes investigated (Figure 4A). Hsa-miR-100 appeared to be slightly down-regulated following G treatment, whereas the other 6 miRNAs investigated resulted to be up-regulated. Even if miR-424 resulted clearly up regulated, the only statistically significant change was observed for hsa-miR-500a-5p, which expression was over 3 times higher following G exposure (Figure 4B). The same over-expression was also observed for the intracellular miRNA (Figure 4C), indicating an effect of G both on intracellular and exosome-released miR-500a-5p.

FIGURE 4

Figure 4 Effect of G on the exosome and intracellular expression of the selected miRNAs. PBMC (10⁶/mL) were treated for 72 h with G (0.1 µg/mL). Exosomes were prepared as described in the Material and Methods section. Total RNA was extracted, retro-transcribed, and target miRNAs were pre-amplified and detected through Real-Time PCR. (A) Fold-changes (2^−ΔΔCt) of the selected miRNAs in exosomes (hsa_miR-10b, hsa_miR-27a, hsa_miR-100, hsa_miR-136, hsa-miR-424, hsa_miR-500a, hsa_let-7f) are reported. (B) Fold-change (2^−ΔCt) of hsa_miR-500a contained in exosomes. (C) Fold-change (2^−ΔCt) of intracellular hsa_miR-500a. hsa-RNU6-2 as negative control. (D) Role of miR-500a in IFN-γ production. PBMC (40’000 cells per well) were silenced with the inhibitor of hsa_miR-500a-5p for 24 h, after medium replacement cells were treated with G (0.1 µg/mL) for 1 h and then with PMA and ionomycin for 72 h. IFN-γ was measured by ELISA in cell-free supernatants. In (B, C) dots represent the individual values of control cells and rhombuses of G-treated cells for each donor. In (D) each dot represents the fold change of IFN-γ release induced by G of each individual, whereas rhombuses represent the fold change of IFN-γ release obtained from antagomir pre-treatment. Each value represents the mean ± SEM, n = 5 donors. Statistical analysis was performed with Paired t test, with *p < 0.05 and **p < 0.01 vs control cells (or G-treated cells for (D).

To investigate if hsa-miR-500a-5p was involved in G-induced inhibition of IFN-γ release, miR-500a-5p was silenced by the use of an antagomir. The pre-incubation with hsa-miR-500a-5p inhibitor was able to restore the production of IFN-γ, in a statistically significant manner (Figure 4D), indicating a role of miR-500a-5p in the observed effects.

Therefore, we highlighted an association between hsa-miR-500a-5p and IFN-γ, which release can be also affected by the pre-treatment with ICI 182,780. To correlate these two phenomena, the expression of hsa-miR-500a-5p induced by G was assessed following ER inhibition. The inhibition of ER was able to statistically significantly reduce the expression of hsa-miR-500a-5p in the intracellular compartment (Figure 5), indicating the implication of ER in G-induced enhancement of hsa-miR-500a-5p expression. Despite the individual variability of miRNA expression, a reduction of hsa-miR-500a-5p expression, following ER inhibition, can be observed.

FIGURE 5

Figure 5 Role of ERα in G-induced intracellular expression of hsa-miR-500a-5p. PBMC (10⁶/mL) were treated with ICI 182,780 (1 µM) for 15 minutes, and subsequently with G (0.1 µg/mL) for 72 h. Total RNA was extracted, retro-transcribed, and hsa-miR-500a-5p detected through Real-Time PCR. Fold-changes (2^−ΔΔCt) are reported, using hsa-RNU6-2 as negative control. Each dot represents the fold change of hsa-miR-500a-5p expression induced by G over control for each individual, whereas rhombuses represent the fold change of hsa-miR-500a-5p expression obtained from cells treated with ICI 182,780 and G over ICI-treated cells. Each value represents the mean ± SEM, n = 5 donors. Statistical analysis was performed with Paired t test, with *p < 0.05 vs G-treated cells.

Discussion

The purpose of this study was to investigate the ability of G to directly affect Th lymphocytes. Results demonstrated that G was able to decrease the Th1/Th2 ratio, inhibiting IFN-γ production. Effects could be reversed by the ER antagonist ICI 182,780, supporting a role of endocrine disturbance, and by miR-500a-5p antagomir. Results are in agreement with the study conducted in mice by Kumar et al. (23) that demonstrated that G-rich air samples promoted Th2 type cytokines.

Even if still controversial, the ability of G to activate ER in vitro has been demonstrated (34, 36), and confirmations also came from experimental animals, showing ERα activation following herbicide exposure in male and female rats (37, 38). The possible role of ER in G-immunotoxicity is strengthened by the ability of estrogens to influence Th1 and Th2 responses, with a promotion toward Th2 activity (39), and to allergic airway inflammation (40). Estrogens can act on immune cells through both ERα and β (41, 42). ERs are involved in shaping the differentiation of Th cells, affecting transcriptional regulation, with possible consequences in inflammation (43). In this study, the role of ER was demonstrated, by the ability of ICI able to restore IFN-γ and IL-4 modulation, and results further support the endocrine disrupting activity of G (44).

G shown a non-monotonic dose-response relationship in its toxic effects against Th. Many endocrine disruptors, and mainly estrogenic substances, display this characteristic (45), strengthening the hypothesis of G considered as an endocrine disruptor. Some possible explanations to this phenomenon could be the cytotoxicity at high concentrations, the different receptor affinity, or receptor desensitization (45, 46).

To further investigate the mechanism underlying G immunotoxicity, a different expression of miRNAs in exosomes by PBMC following G exposure was investigated. MiRNAs are single-stranded RNAs of about 22 nucleotides, that function as post-transcriptional regulators, through the interference with mRNA processes (47). They can affect both mRNA stability and target mRNA for the subsequent degradation (48). The majority of miRNAs are localized in the intracellular compartment, however several miRNAs in various body fluids have been described (49). Exosomes are miRNAs carriers (50), of great relevance in pharmacology and toxicology (51). Notably, circulating miRNAs have been described as potential biomarkers of autoimmune and inflammatory disease (reviewed by 52). Considering the important role of miRNAs in the regulation of gene expression and their possible application as biomarkers and as cell-cell communication mechanism, we examined the differential expression of several miRNAs in farmers occupationally exposed to G (manuscript in preparation). Results presented in this work suggest miR-500a-5p as a possible target of G in immune cells. Only few studies investigated the effect of G on non-coding RNAs, mainly conducted in mice. Ji et al. (53) associated perinatal G exposure in mice with a differential miRNA expression in the prefrontal cortex (55 miRNAs up-regulated and 19 miRNAs down-regulated), mainly involved in neural development. Also, circular RNAs were affected, with 330 up-regulated and 333 down-regulated (54). Supporting this evidence, more recently, perinatal G exposure in mice was associated with aberrant expression of lncRNAs, highlighting a link with impaired neuronal development (55). Also changes in miRNA profile in female rat liver were assessed following G and G-based herbicide treatment, and interestingly these miRNAs were associated with carcinogenesis, like miR-17 and miR-22 (56). Furthermore, an in vitro study performed on carp lymphocytes shown that G-induced lymphocyte apoptosis was mediated by the regulation of miR-203 and PI3K/AKT pathway (57). No studies examined the effect of G on miRNAs in human immune cells, in this regard this study represents the first evidence that G can affect human miRNAs expression linked to the immune system. In particular, a role of miR-500a could be demonstrated, linking miR-500a and lower IFN-γ production. MiR-500a-5p is not a well-characterized miRNA; to date the few available evidence address a role in oxidative stress in breast cancer (58) and in the promotion of breast cancer cell proliferation (59). It has also been described as oncogenic in colorectal cancer (60), as well as responsible for proliferation and metastasis of hepatocarcinoma (61). Based on TargetScan (v7.2; targetscan.org), a useful database to predict miRNA target sites (62), hsa-miR-500a-5p could hypothetically inhibit several mRNAs involved in the differentiation of Th cells towards Th1, preventing it. Among possible targets there are STAT1 (signal transducer and activator of transcription 1), TBX21 (T-box 21; or T-bet), SOCS5 (suppressor of cytokine signaling 5) and ANXA1 (annexin A1), which are all positive regulators of Th1. The up-regulation of miR-500a could be indeed in line with the decrease in Th1 as demonstrated in the current study. Furthermore, we speculate on the existence of a dependance between miR-500a expression and ER, since its inhibition reflected in a lower expression of hsa-miR-500a-5p. Preliminary results indicate indeed a reduced expression of STAT1 following G exposure, which could be reversed by ICI (data not shown), indicating STAT1 as a relevant target to explain G immunotoxicity.

Overall, data obtained indicate that G at biological relevant concentrations can directly affect Th cells, disturbing the balance between Th1 and Th2, through the ER pathway. We also demonstrated a role of miR-500a-5p in G-induced reduced IFN-γ production. Additional studies are required to identify the genes targeted by miR-500a-5p involved in Th1 differentiation and their link with ER.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author Contributions

AM: manuscript writing, reviewing, and editing. MI, VG, CC, SM-R, and EC: manuscript reviewing and editing. All authors contributed to the article and approved the submitted version.

Funding

Research has been supported by Ministero dell’Istruzione, dell’Università e della Ricerca (PRIN 2017, Project number 2017MLC3NF).

Conflict of Interest

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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.854837/full#supplementary-material

References

1. Duke SO, Powles SB. Glyphosate: A Once-in-a-Century Herbicide. Pest Manage Sci (2008) 64(4):319–25. doi: 10.1002/ps.1518

CrossRef Full Text | Google Scholar

2. Benbrook CM. Trends in Glyphosate Herbicide Use in the United States and Globally. Environ Sci Europe (2016) 28(1):3. doi: 10.1186/s12302-016-0070-0

CrossRef Full Text | Google Scholar

3. Mesnage R, Defarge N, Spiroux de Vendômois J, Séralini GE. Potential Toxic Effects of Glyphosate and its Commercial Formulations Below Regulatory Limits. Food Chem Toxicol: an Int J Published Br Ind Biol Res Assoc (2015) 84:133–53. doi: 10.1016/j.fct.2015.08.012

CrossRef Full Text | Google Scholar

4. Acquavella JF, Alexander BH, Mandel JS, Gustin C, Baker B, Chapman P, et al. Glyphosate Biomonitoring for Farmers and Their Families: Results From the Farm Family Exposure Study. Environ Health Perspect (2004) 112(3):321–6. doi: 10.1289/ehp.6667

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Amrhein N, Deus B, Gehrke P, Steinrücken HC. The Site of the Inhibition of the Shikimate Pathway by Glyphosate: Ii. Interference of Glyphosate With Chorismate Formation In Vivo and In Vitro. Plant Physiol (1980) 66(5):830–4. doi: 10.1104/pp.66.5.830

PubMed Abstract | CrossRef Full Text | Google Scholar

6. EFSA. Conclusion on the Peer Review of the Pesticide Risk Assessment of the Active Substance Glyphosate. EFSA J (2015) 13(11):4302. doi: 10.2903/j.efsa.2015.4302

CrossRef Full Text | Google Scholar

7. US EPA. Draft Human Health Risk Assessment in Support of Registration Review, in: Glyphosate. Environmental Protection Agency, Office of Prevention, Pesticides and Toxic Substances, Office of Pesticide Programs (2017). Washington, DC: U.S. Government Printing Office. Available at: https://www.regulations.gov/document?D=EPA-HQ-OPP-2009-0361-0068 (Accessed February 02, 2022).

Google Scholar

8. Myers JP, Antoniou MN, Blumberg B, Carroll L, Colborn T, Everett LG, et al. Concerns Over Use of Glyphosate-Based Herbicides and Risks Associated With Exposures: A Consensus Statement. Environ Health (2016) 15:19. doi: 10.1186/s12940-016-0117-0

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Maggi F, la Cecilia D, Tang F, McBratney A. The Global Environmental Hazard of Glyphosate Use. Sci Total Environ (2020) 717:137167. doi: 10.1016/j.scitotenv.2020.137167

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Martinez TT, Long WC, Hiller R. Comparison of the Toxicology of the Herbicide Roundup by Oral and Pulmonary Routes of Exposure. Proc Western Pharmacol Soc (1990) 33:193–7.

Google Scholar

11. Connolly A, Coggins MA, Galea KS, Jones K, Kenny L, McGowan P, et al. Evaluating Glyphosate Exposure Routes and Their Contribution to Total Body Burden: A Study Among Amenity Horticulturalists. Ann Work Exposures Health (2019) 63(2):133–47. doi: 10.1093/annweh/wxy104

CrossRef Full Text | Google Scholar

12. Eriksson M, Hardell L, Carlberg M, Akerman M. Pesticide Exposure as Risk Factor for non-Hodgkin Lymphoma Including Histopathological Subgroup Analysis. Int J Cancer (2008) 123(7):1657–63. doi: 10.1002/ijc.23589

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Mesnage R, Arno M, Costanzo M, Malatesta M, Séralini GE, Antoniou MN. Transcriptome Profile Analysis Reflects Rat Liver and Kidney Damage Following Chronic Ultra-Low Dose Roundup Exposure. Environ Health (2015) 14:70. doi: 10.1186/s12940-015-0056-1

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Mesnage R, Renney G, Séralini GE, Ward M, Antoniou MN. Multiomics Reveal non-Alcoholic Fatty Liver Disease in Rats Following Chronic Exposure to an Ultra-Low Dose of Roundup Herbicide. Sci Rep (2017) 7:39328. doi: 10.1038/srep39328

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Pham TH, Derian L, Kervarrec C, Kernanec PY, Jégou B, Smagulova F, et al. Perinatal Exposure to Glyphosate and a Glyphosate-Based Herbicide Affect Spermatogenesis in Mice. Toxicol Sci (2019) 169(1):260–71. doi: 10.1093/toxsci/kfz039

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hu J, Lesseur C, Miao Y, Manservisi F, Panzacchi S, Mandrioli D, et al. Low-Dose Exposure of Glyphosate-Based Herbicides Disrupt the Urine Metabolome and its Interaction With Gut Microbiota. Sci Rep (2021) 11(1):3265. doi: 10.1038/s41598-021-82552-2

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Corsini E, Sokooti M, Galli CL, Moretto A, Colosio C. Pesticide Induced Immunotoxicity in Humans: A Comprehensive Review of the Existing Evidence. Toxicology (2013) 307:123–35. doi: 10.1016/j.tox.2012.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Lee GH, Choi KC. Adverse Effects of Pesticides on the Functions of Immune System. Comparative Biochemistry and Physiology. Toxicol Pharmacol: CBP (2020) 235:108789. doi: 10.1016/j.cbpc.2020.108789

CrossRef Full Text | Google Scholar

19. Rajak P, Ganguly A, Sarkar S, Mandi M, Dutta M, Podder S, et al. Immunotoxic Role of Organophosphates: An Unseen Risk Escalating SARS-CoV-2 Pathogenicity. Food Chem Toxicol (2021) 149:112007. doi: 10.1016/j.fct.2021.112007

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Blom J, Ottaviani E. Immune-Neuroendocrine Interactions: Evolution, Ecology, and Susceptibility to Illness. Med Sci Monit Basic Res (2017) 23:362–7. doi: 10.12659/msmbr.907637

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Fung TC, Olson CA, Hsiao EY. Interactions Between the Microbiota, Immune and Nervous Systems in Health and Disease. Nat Neurosci (2017) 20(2):145–55. doi: 10.1038/nn.4476

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Maddalon A, Galbiati V, Colosio C, Mandić-Rajčević S, Corsini E. Glyphosate-Based Herbicides: Evidence of Immune-Endocrine Alteration. Toxicology (2021) 459:152851. doi: 10.1016/j.tox.2021.152851

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Kumar S, Khodoun M, Kettleson EM, McKnight C, Reponen T, Grinshpun SA, et al. Glyphosate-Rich Air Samples Induce IL-33, TSLP and Generate IL-13 Dependent Airway Inflammation. Toxicology (2014) 325:42–51. doi: 10.1016/j.tox.2014.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Pandey A, Dhabade P, Kumarasamy A. Inflammatory Effects of Subacute Exposure of Roundup in Rat Liver and Adipose Tissue. Dose-response (2019) 17(2):1559325819843380. doi: 10.1177/1559325819843380

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Tang Q, Tang J, Ren X, Li C. Glyphosate Exposure Induces Inflammatory Responses in the Small Intestine and Alters Gut Microbial Composition in Rats. Environ Pollut (2020) 261:114129. doi: 10.1016/j.envpol.2020.114129

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Qiu S, Fu H, Zhou R, Yang Z, Bai G, Shi B. Toxic Effects of Glyphosate on Intestinal Morphology, Antioxidant Capacity and Barrier Function in Weaned Piglets. Ecotoxicol Environ Saf (2020) 187:109846. doi: 10.1016/j.ecoenv.2019.109846

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Slager RE, Poole JA, LeVan TD, Sandler DP, Alavanja MC, Hoppin JA. Rhinitis Associated With Pesticide Exposure Among Commercial Pesticide Applicators in the Agricultural Health Study. Occup Environ Med (2009) 66(11):718–24. doi: 10.1136/oem.2008.041798

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Slager RE, Simpson SL, Levan TD, Poole JA, Sandler DP, Hoppin JA. Rhinitis Associated With Pesticide Use Among Private Pesticide Applicators in the Agricultural Health Study. J Toxicol Environ Health Part A (2010) 73(20):1382–93. doi: 10.1080/15287394.2010.497443

CrossRef Full Text | Google Scholar

29. Hoppin JA, Umbach DM, Long S, London SJ, Henneberger PK, Blair A, et al. Pesticides are Associated With Allergic and Non-Allergic Wheeze Among Male Farmers. Environ Health Perspect (2017) 125(4):535–43. doi: 10.1289/EHP315

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mesnage R, Teixeira M, Mandrioli D, Falcioni L, Ducarmon QR, Zwittink RD, et al. Use of Shotgun Metagenomics and Metabolomics to Evaluate the Impact of Glyphosate or Roundup MON 52276 on the Gut Microbiota and Serum Metabolome of Sprague-Dawley Rats. Environ Health Perspect (2021) 129(1):17005. doi: 10.1289/EHP6990

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Mendler A, Geier F, Haange SB, Pierzchalski A, Krause JL, Nijenhuis I, et al. Mucosal-Associated Invariant T-Cell (MAIT) Activation is Altered by Chlorpyrifos- and Glyphosate-Treated Commensal Gut Bacteria. J Immunotoxicol (2020) 17(1):10–20. doi: 10.1080/1547691X.2019.1706672

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Poulton TA, Gallagher A, Potts RC, Beck JS. Changes in Activation Markers and Cell Membrane Receptors on Human Peripheral Blood T Lymphocytes During Cell Cycle Progression After PHA Stimulation. Immunology (1988) 64(3):419–25.

PubMed Abstract | Google Scholar

33. Kongtip P, Nankongnab N, Phupancharoensuk R, Palarach C, Sujirarat D, Sangprasert S, et al. Glyphosate and Paraquat in Maternal and Fetal Serums in Thai Women. J Agromedicine (2017) 22(3):282–9. doi: 10.1080/1059924X.2017.1319315

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Thongprakaisang S, Thiantanawat A, Rangkadilok N, Suriyo T, Satayavivad J. Glyphosate Induces Human Breast Cancer Cells Growth via Estrogen Receptors. Food Chem Toxicol (2013) 59:129–36. doi: 10.1016/j.fct.2013.05.057

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Robertson JF. ICI 182,780 (Fulvestrant)–the First Oestrogen Receptor Down-Regulator–Current Clinical Data. Br J Cancer (2001) 85 Suppl 2(Suppl 2):11–4. doi: 10.1054/bjoc.2001.1982

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Mesnage R, Phedonos A, Biserni M, Arno M, Balu S, Corton JC, et al. Evaluation of Estrogen Receptor Alpha Activation by Glyphosate-Based Herbicide Constituents. Food Chem Toxicol (2017) 108(Pt A):30–42. doi: 10.1016/j.fct.2017.07.025

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Guerrero Schimpf M, Milesi MM, Ingaramo PI, Luque EH, Varayoud J. Neonatal Exposure to a Glyphosate Based Herbicide Alters the Development of the Rat Uterus. Toxicology (2017) 376:2–14. doi: 10.1016/j.tox.2016.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Altamirano GA, Delconte MB, Gomez AL, Ingaramo PI, Bosquiazzo VL, Luque EH, et al. Postnatal Exposure to a Glyphosate-Based Herbicide Modifies Mammary Gland Growth and Development in Wistar Male Rats. Food Chem Toxicol (2018) 118:111–8. doi: 10.1016/j.fct.2018.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Phiel KL, Henderson RA, Adelman SJ, Elloso MM. Differential Estrogen Receptor Gene Expression in Human Peripheral Blood Mononuclear Cell Populations. Immunol Lett (2005) 97(1):107–13. doi: 10.1016/j.imlet.2004.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Watanabe Y, Tajiki-Nishino R, Tajima H, Fukuyama T. Role of Estrogen Receptors α and β in the Development of Allergic Airway Inflammation in Mice: A Possible Involvement of Interleukin 33 and Eosinophils. Toxicology (2019) 411:93–100. doi: 10.1016/j.tox.2018.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Shim GJ, Gherman D, Kim HJ, Omoto Y, Iwase H, Bouton D, et al. Differential Expression of Oestrogen Receptors in Human Secondary Lymphoid Tissues. J Pathol (2006) 208(3):408–14. doi: 10.1002/path.1883

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Straub RH. The Complex Role of Estrogens in Inflammation. Endocr Rev (2007) 28(5):521–74. doi: 10.1210/er.2007-0001

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Mohammad I, Starskaia I, Nagy T, Guo J, Yatkin E, Väänänen K, et al. Estrogen Receptor α Contributes to T Cell-Mediated Autoimmune Inflammation by Promoting T Cell Activation and Proliferation. Sci Signaling (2018) 11(526):eaap9415. doi: 10.1126/scisignal.aap9415

CrossRef Full Text | Google Scholar

44. Muñoz JP, Bleak TC, Calaf GM. Glyphosate and the Key Characteristics of an Endocrine Disruptor: A Review. Chemosphere (2021) 270:128619. doi: 10.1016/j.chemosphere.2020.128619

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Vandenberg LN, Colborn T, Hayes TB, Heindel JJ, Jacobs DRJR, Lee DH, et al. Hormones and Endocrine-Disrupting Chemicals: Low-Dose Effects and Nonmonotonic Dose Responses. Endocr Rev (2012) 33(3):378–455. doi: 10.1210/er.2011-1050

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Lagarde F, Beausoleil C, Belcher SM, Belzunces LP, Emond C, Guerbet M, et al. Non-Monotonic Dose-Response Relationships and Endocrine Disruptors: A Qualitative Method of Assessment. Environ Health: Global Access Sci Source (2015) 14:13. doi: 10.1186/1476-069X-14-13

CrossRef Full Text | Google Scholar

47. Kim VN, Han J, Siomi MC. Biogenesis of Small RNAs in Animals. Nat Rev Mol Cell Biol (2009) 10(2):126–39. doi: 10.1038/nrm2632

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Ghildiyal M, Zamore PD. Small Silencing RNAs: An Expanding Universe. Nat Rev Genet (2009) 10(2):94–108. doi: 10.1038/nrg2504

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Lawrie CH, Gal S, Dunlop HM, Pushkaran B, Liggins AP, Pulford K, et al. Detection of Elevated Levels of Tumour-Associated microRNAs in Serum of Patients With Diffuse Large B-Cell Lymphoma. Br J Haematol (2008) 141(5):672–5. doi: 10.1111/j.1365-2141.2008.07077.x

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are Transported in Plasma and Delivered to Recipient Cells by High-Density Lipoproteins. Nat Cell Biol (2011) 13(4):423–33. doi: 10.1038/ncb2210

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Bowers EC, Hassanin A, Ramos KS. In Vitro Models of Exosome Biology and Toxicology: New Frontiers in Biomedical Research. Toxicology. Vitro (2020) 64:104462. doi: 10.1016/j.tiv.2019.02.016

CrossRef Full Text | Google Scholar

52. Siddeek B, Inoubli L, Lakhdari N, Rachel PB, Fussell KC, Schneider S, et al. MicroRNAs as Potential Biomarkers in Diseases and Toxicology. Mutat Res Genet Toxicol Environ Mutagenesis (2014) 764-765:46–57. doi: 10.1016/j.mrgentox.2014.01.010

CrossRef Full Text | Google Scholar

53. Ji H, Xu L, Wang Z, Fan X, Wu L. Differential microRNA Expression in the Prefrontal Cortex of Mouse Offspring Induced by Glyphosate Exposure During Pregnancy and Lactation. Exp Ther Med (2018) 15(3):2457–67. doi: 10.3892/etm.2017.5669

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Yu N, Tong Y, Zhang D, Zhao S, Fan X, Wu L, et al. Circular RNA Expression Profiles in Hippocampus From Mice With Perinatal Glyphosate Exposure. Biochem Biophys Res Commun (2018) 501(4):838–45. doi: 10.1016/j.bbrc.2018.04.200

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Fan X, Wang D, Shen X, Qiu J, Wu L, Yan J, et al. Identification of lncRNA Expression Profiles and Analysis of ceRNA in the Hippocampus of Perinatal Glyphosate-Exposed Mice. Int J Dev Neurosci (2021) 81(4):312–23. doi: 10.1002/jdn.10102

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Mesnage R, Ibragim M, Mandrioli D, Falcioni L, Tibaldi E, Belpoggi F, et al. Comparative Toxicogenomics of Glyphosate and Roundup Herbicides by Mammalian Stem Cell-Based Genotoxicity Assays and Molecular Profiling in Sprague-Dawley Rats. Toxicol Sci: an Off J Soc Toxicol (2021) 186(1):83–101. doi: 10.1093/toxsci/kfab143

CrossRef Full Text | Google Scholar

57. Wang W, Liu Q, Zhang T, Chen L, Li S, Xu S. Glyphosate Induces Lymphocyte Cell Dysfunction and Apoptosis via Regulation of miR-203 Targeting of PIK3R1 in Common Carp (Cyprinus Carpio L.). Fish Shellfish Immunol (2020) 101:51–7. doi: 10.1016/j.fsi.2020.03.047

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Degli Esposti D, Aushev VN, Lee E, Cros MP, Zhu J, Herceg Z, et al. miR-500a-5p Regulates Oxidative Stress Response Genes in Breast Cancer and Predicts Cancer Survival. Sci Rep (2017) 7(1):15966. doi: 10.1038/s41598-017-16226-3

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Chen B, Sang Y, Song X, Zhang D, Wang L, Zhao W, et al. Exosomal miR-500a-5p Derived From Cancer-Associated Fibroblasts Promotes Breast Cancer Cell Proliferation and Metastasis Through Targeting USP28. Theranostics (2021) 11(8):3932–47. doi: 10.7150/thno.53412

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Li Y, Chen Z. The Oncogenic Role of microRNA-500a in Colorectal Cancer. Oncol Lett (2020) 19(3):1799–805. doi: 10.3892/ol.2020.11291

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Zhao Y, Wang Y, Wang Y. Up-Regulated miR-500a Enhances Hepatocarcinoma Metastasis by Repressing PTEN Expression. Biosci Rep (2017) 37(6):BSR20170837. doi: 10.1042/BSR20170837

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting Effective microRNA Target Sites in Mammalian mRNAs. eLife (2015) 4:e05005. doi: 10.7554/eLife.05005

CrossRef Full Text | Google Scholar