Effect of TCDD exposure in adult female and male mice on the expression of miRNA in the ovaries and testes and associated reproductive functions

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is an environmental contaminant found widely across the world. While animal and human studies have shown that exposure to TCDD may cause significant alterations in the reproductive tract, the effect of TCDD on the expression of miRNA in the reproductive organs has not been previously tested. In the current study, we exposed adult female or male mice to TCDD or vehicle and bred them to study the impact on reproduction. The data showed that while TCDD treatment of females caused no significant change in litter size, it did alter the survival of the pups. Also, TCDD exposure of either the male or female mice led to an increase in the gestational period. While TCDD did not alter the gross morphology of the ovaries and testes, it induced significant alterations in the miRNA expression. The ovaries showed the differential expression of 426 miRNAs, of which 315 miRNAs were upregulated and 111 miRNA that were downregulated after TCDD exposure when compared to the vehicle controls. In the testes, TCDD caused the differential expression of 433 miRNAs, with 247 miRNAs upregulated and 186 miRNAs downregulated. Pathway analysis showed that several of these dysregulated miRNAs targeted reproductive functions. The current study suggests that the reproductive toxicity of TCDD may result from alterations in the miRNA expression in the reproductive organs. Because miRNAs also represent one of the epigenetic pathways of gene expression, our studies suggest that the transgenerational toxicity of TCDD may also result from dysregulation in the miRNAs.


Introduction
, is an environmental pollutant that is considered to be the most persistent and strongest endocrine disruptor (Hites, 2011;Zhang et al., 2021).TCDD is also a historically controversial chemical contaminant that was produced during the high-temperature synthesis of organochloride herbicides, such as Agent Orange, which was heavily utilized during the Vietnam War (Kim et al., 2003; OPEN ACCESS EDITED BY David R. Wallace, Oklahoma State University Center for Health Sciences, United States NASEM, 2018).Veterans exposed to TCDD contamination in Agent Orange during the war have exhibited an increased occurrence of various cancers (Chang et al., 2014;Chang et al., 2017), dementia (Martinez et al., 2021), skin disorders like porphyria and chloracne (Patterson et al., 2016), as well as endocrine and reproductive system disease (Gaspari et al., 2021).
TCDD toxicity is mediated through the aryl hydrocarbon receptor (AhR) which is translocated to the nucleus after ligand binding.Here, it forms a heterodimer with the aryl hydrocarbon receptor nuclear translocator (ARNT) protein which can then directly interact with genes containing Dioxin Response Elements (DREs) (Enan and Matsumara, 1996;Rothhammer and Quintana, 2019).These DREs alter gene expression resulting in changes in physiopathology including both immunosuppressive and inflammatory processes (Clark et al., 1981;Dominques-Acosta et al., 2018;Vasquez-Gomez et al., 2022).TCDD exposure has been observed to disrupt normal immune system development and functions such as increasing thymic atrophy (De Heer et al., 1994;Frazier et al., 1994;Kamath et al., 1997;Camacho et al., 2005), alterations to thymocyte development (Fine et al., 1990) and increased T-cell apoptosis (Singh et al., 2020).Additionally, hormone synthesis and circulation disruption has also been observed, following TCDD exposure, in the thyroid (Kohn, 2000), hypothalamus (Fetissov et al., 2004), pituitary and adrenal glands (Bestervelt et al., 1993).
In addition to modulating gene expression and hormone activity, direct effects on reproductive success have been observed after exposure to dioxin in murine models.Previous research from our laboratory has revealed TCDD induces a loss of mitochondrial membrane potential in epididymal spermatozoa, potentially reducing sperm quality and motility (Fisher et al., 2005) as well as reducing testes' weight and circulating testosterone levels (Choi et al., 2008;Zhang et al., 2021).Paternal exposure to TCDD has been shown to distort the sex ratio of the offspring, increasing the proportion of females per litter (Bircsak et al., 2020).Females exposed to TCDD had increased rates of endometriosis and ovarian cancers (Johnson et al., 1997;Bruner-Tran et al., 2017) as well as altered ovarian function (Tischkau et al., 2011), and reduced estradiol synthesis (Valdez et al., 2009).Epigenetic changes induced by TCDD are heritable and directly observable at the germ line (Patrizi and Siciliani de Cumis, 2018).By studying these alterations at the gonad level, one can better characterize the long-term and transgenerational risk of exposure.
Increased expression of genes containing DREs is not the only transcriptional disruption induced by TCDD.TCDD can also alter the expression of mRNA through dysregulation of miRNA (Singh et al., 2012;AlGhezi et al., 2019;Neamah et al., 2019) which are short, ~20 bp long, non-coding RNAs that regulate gene translation by silencing target mRNA sequences (O'Brien et al., 2018).However, some miRNA express DREs in their 3′ UTR, thereby enabling TCDD to directly alter their expression through DRE-dependent pathways (Disner et al., 2021).miRNAs serve as critical epigenetic modulators, altering the protein levels of the target mRNAs without changing the gene sequences (Yao et al., 2019).
The effect of in vivo exposure to TCDD on miRNA expression in the reproductive organs and its impact on reproductive functions has not been previously studied.However, there are a few reports on the in vitro effects of TCDD on miRNA expression in cells such as Chinese Hamster Ovary (CHO) cells (Sadowska, et al., 2021) and human Sertoli cells (Ribeiro, et al., 2018), the former study showing no significant changes in miRNA expression while the latter showing significant alterations.The present study investigated the miRNA expression in the ovaries and testes of adult mice exposed to TCDD and the genes targeted by those miRNAs.Our data demonstrate that TCDD alters the expression of a large number of miRNAs and several of them target the gene expression involved in the regulation of reproductive functions.

Mice
Male and Female C57BL/6J mice aged 5-6 weeks were purchased from Jackson Laboratories (Bar Harbor, ME, United States 04609) and housed at the AAALAC-accredited DLAR animal facility at the University of South Carolina School of Medicine (Columbia, SC, United States 20208) under specific pathogen-free conditions.Before and throughout the experiment, animals were maintained in temperature and humidity-controlled rooms on 12-h light/dark cycles and provided with access to food and normal rodent chow ad libitum.Female mice were housed together in cages of five without male contact before random pairing and/or treatment to allow for estrus cycles to synchronize.All experiments were performed in accordance with ethical standards approved by the University of South Carolina Institutional Animal Care and Use Committee.

TCDD exposure and mating pairs
TCDD was generously provided by Dr. Steve Safe (Institute of Biosciences & Technology, Texas A&M Health Science Center, College Station, TX, United States) and was dissolved in DMSO to a working stock concentration of 100 μg/mL.The stock was combined with medical-grade corn oil for a total volume of 100 µL injected intraperitoneally (I.P).Mice were dosed at 10 μg/kg body weight similar to other studies which have used similar doses (Besteman, et al., 2007;Li, et al., 2022).Male and female mice were either treated with TCDD or vehicle, the latter consisting of an equal amount of DMSO in corn oil, then immediately paired with one male to one female per cage for initial reproductive success assessment.Treatment and pairing design are described in Figure 1A.For miRNA expression analysis, male and female mice were treated as previously described, kept with their original same-sex cage mates, and sacrificed 72-h post-exposure.

Reproductive parameter assessment
After treatment and pairing, the female mice were examined for the vaginal plug, the sighting of which was considered as possible day 1 of pregnancy.Next, the mated pairs were left undisturbed until females were visibly pregnant at which time males were removed for parturition.Daily monitoring of females was conducted and measurements such as days since pairing, number of pups born in a litter, and number of neonates found deceased in the cage were recorded until weaning age was reached at 4 weeks.

Histology
Mice were euthanized with an overdose of inhaled isoflurane and testes and ovaries were excised from the male and female mice, respectively, rinsed with PBS, and immersed in 4% paraformaldehyde (PFA).Paraffin embedding, cutting, and mounting, and H&E staining was performed on-site through the University of South Carolina School of Medicine Instrument Resource Facility (IRF).Imaging was conducted on a Keyence BZ-X800 Automated High-resolution Fluorescence Microscope (Keyence Corporation of America, Itasca, IL, United States 60143).

Tissue processing and miRNA isolation
Testes and ovaries were excised and placed in blender bags containing RPMI medium containing FBS (10% v/v) and an antibiotic cocktail including penicillin and streptomycin.All tissues were kept on ice during processing.Tissues were homogenized using a Seward ™ Stomacher ™ Model 80 Biomaster Blender 110 V (Seward Inc, Bohemia, NY, United States) and filtered using 100 µm mesh filters.The resulting filtrate was centrifuged at 350 × g for 7 min and the supernatant was discarded.The remaining cell pellet was resuspended in 700 µL of QIAzol ™ Lysis reagent (QIAGEN, Germantown, MD,   United States 20874).miRNA was isolated using a QIAGEN miRNEasy Micro Kit (QIAGEN, Germantown, MD, United States 20874) according to the manufacturer's instructions.miRNA microarray miRNA arrays were performed as described previously (Al-Ghezi et al., 2019).Purified miRNA samples isolated from testes and ovaries harvested from mice 72 h post-exposure to TCDD were labeled using the Affymetrix ® FlashTag ™ Biotin HSR RNA Labeling Kit using the manufacturer's specifications for 100 format miRNA 4.0 arrays.GeneChip ™ Hybridization, Wash, and Stain kit was used in conjunction with the Affymetrix ® Hybridization Oven and Affymetrix ® Fluidics Station on the 450 protocol to prepare Affymetrix ® GeneChip ® miRNA 4.0 arrays.Prepared arrays were then analyzed using the GeneChip ™ Scanner 3,000 7G (All products in this section are sourced from ThermoFisher Scientific, Waltham, MA, United States 02451).Microarray intensity was quantified, and differential expression was analyzed using Applied Biosystems Transcription Analysis Console (TAC) 4.0 software.Expression data were analyzed using Ingenuity Pathway Analysis (IPA) and miRNA targets were determined using the current target prediction database miRDB (Chen and Wang, 2020), and target prediction algorithm,

FIGURE 1
Effect of TCDD exposure on reproductive functions.There were four groups of mice: 1) Adult female mice treated with TCDD and mated with male mice treated with vehicle, 2) Adult female mice treated with vehicle and mated with male mice treated with TCDD, 3) Both male and female mice treated with TCDD, 4) Both male and female mice treated with vehicle.Such mice were then paired and following pregnancy, the pregnant mice were observed for the number of pups delivered, the percentage of pup survival per litter for up to 4 weeks, and the number of days taken from the sighting of vaginal plug to the birth of the pups.Reproductive success parameters by treatment group were measured as follows: (B) litter size at the time of birth, (C) post-natal pup survivorship to weaning age, (D) time from pairing to birth.Significance was determined using Student's t-test with n = 10.The asterisk denotes as follows: * <0.05, **<0.01.
MirTarget as well as web-based TargetScan.MicroRNA raw data have been deposited to Geo database (accession number is GSE241576 and can be accessed by https://www.ncbi.nlm.nih.gov/geo/info/linking.html) of National Institute of Health (NIH).
Real-time quantitative PCR (RT-qPCR) to validate the expression of miRNAs and genes miRNAs and target mRNAs of interest were validated using Real-Time PCR.To synthesize cDNA, iScript ™ cDNA Synthesis Kit from Bio- RAD (BIO-RAD, Herculeus, VCA) was used.RT-qPCR using SYBR ® Green qPCR Supermix from Qiagen (Qiagen, Valencia, CA) was used following the protocol of the company.Real-qPCR was performed on a BIO-RAD CFX Connect RT-PCR Detection System (BIO-RAD, Hercules, CA).RT-qPCR was performed for 40 cycles and using the following conditions: initial activation step (5 min at 95 °C), denaturing temperature (15 s at 94 °C), annealing temperature (30 s at 55 °C), and extension temperature and fluorescence data collection (30 s at 70 °C).Normalized expression (NE) of miRs and genes were calculated using NE ¼ 2_DDCt, where Ct is the threshold cycle to detect fluorescence.The expression of miRNAs was normalized against internal control for miRNA and fold change of miRNAs expression was calculated against the internal control SNORD96A, and the treatment group (TCDD) was compared with the vehicle group.Similarly, the value of genes was normalized against the housekeeping gene (18 S) and fold-change of genes was calculated against 18 S and the TCDD group was compared with the vehicle group.To define significant differences in the expression of miRNAs and genes, ANOVA was performed using GraphPad version 6.0 (GraphPad Software, Inc., San Diego, CA, United States), and differences between the groups were considered significant when p < 0.05.

Data analysis and statistics
We used groups of 10 mice injected with TCDD or the vehicle which were then allowed to mate for studies on reproduction.In other studies, we used groups of 5 mice for studies on the ovaries and testes.For miRNA microarray analysis, we pooled the ovaries or testes from 5 animals while the validation of miRNA using RT-qPCR was performed on individual mice with n = 5.The statistical significance was determined using a one-way ANOVA followed by Tukey's post hoc test.

Results
Reproductive success is modulated by parental exposure to TCDD Wild type C57BL/6J male and female mice were exposed to Vehicle or TCDD as described in Methods, paired, and allowed to mate as described in Figure 1A.There were 4 groups of mice: 1) Female mice exposed to TCDD mated with male mice treated with Vehicle (Maternal).2) Male mice exposed to TCDD mated with female mice treated with Vehicle (Paternal).3) Both male and female mice exposed to TCDD (Dual).4) Both male and female mice exposed to Vehicle only (Vehicle).
The pregnant mice were observed for the number of pups delivered, the percentage of pup survival per litter up to 4 weeks of observation, and the number of days taken from pairing to the birth of the pups.The maternal, paternal, or dual exposure to TCDD did not significantly alter the litter size (# of pups/litter) (Figure 1B).However, the survival of the pups was significantly impacted by TCDD exposure.Specifically, maternal but not paternal exposure to TCDD caused a significant decrease in the survival of pups to weaning age when compared to the Vehicle group (Figure 1C).The non-surviving pups were born alive, but it was unclear if subsequent post-natal mortality was due to health complications or live cannibalism by the mother.Additionally, the gestation period was significantly delayed following maternal exposure to TCDD but not after paternal exposure when compared to the Vehicle controls (Figure 1D).Together, these data suggested that maternal but not paternal exposure to TCDD significantly alters the pregnancy and the survival of the pups.

TCDD induces changes in miRNA expression in the ovaries
Gross and histologic examination of the ovaries of TCDDtreated mice showed no significant differences, and several ovarian follicles at various stages of development were observed in the cortex, that were similar in the TCDD and the Vehicle controls (Figure 2A).Next, we investigated if maternal exposure to TCDD would alter the expression of miRNA in the ovaries by using iGeneChip miRNA microarray.These studies revealed the differential expression of 426 miRNAs in the ovaries, of which 315 miRNAs were upregulated and 111 miRNAs were downregulated after TCDD exposure by ≥ 2-fold when compared to the Vehicle control (Figure 2B).The heatmap of some of the most dysregulated miRNAs has been depicted in Figure 2C, which shows marked differences between the two groups.
Ingenuity Pathway Analysis of some select highly dysregulated miRNA revealed how these miRNAs were predicted to target various upstream regulators and downstream mRNA targets involved with reproductive disease and normal fetal development (Figure 3A).The predicted miRNA/mRNA interactions were also confirmed through alignment using TargetScan (Figure 3B).For example, miR-34c was downregulated in the TCDD-treated group which targeted Tumor protein 53 apoptosis inducing protein 1 (TP53AIP1) and Tumor protein 53 inducible protein 11 (TP53I11), key regulatory molecules involved in the maintenance of homeostasis when the cells are in a state of stress and shown to be involved in preeclampsia (Ali et al., 2021).TCDD also decreased the levels of miR-192 which targeted epiregulin (EREG) known to stimulate oocytes (Romero and Smitz, 2009).TCDD upregulated miR-362 and miR-194 which targeted Alx Homebox Protein 1 (ALX1) and Ras family protein 2 B (RAP2B), respectively.Loss of ALX1 function is associated with birth defects (Iyyanar et al., 2022) while RAP2B is involved in cell proliferation and migration (Di et al., 2015).
To validate the miRNA array data and to confirm the effect of altered miRNA expression on the target genes, we performed RT-qPCR of miRNAs of interest and their potential targets (Figure 3C).The data showed that TCDD-exposed ovaries had decreased levels of miR-34c which correlated with increased expression of TP53AIPI and TP53I11 (Figure 3C).Additionally, TCDDD exposure led to increased expression of miR-362 and miR-194 which correlated with decreased expression of ALX1 and RAP2B (Figure 3C).

TCDD induces changes in miRNA expression in the testes
Next, we investigated the effect of TCDD exposure on miRNA expression in the testes.No significant differences in gross morphology of the testes from male mice treated with TCDD were observed when compared to the Vehicle controls, and histologically, the TCDD-treated groups showed normal spermatogenesis (Figure 4A).
Next, we investigated the effect of TCDD exposure on miRNA expression in the testes.Between males treated with Vehicle or TCDD, no significant differences in gross morphology were observed and normal spermatogenesis was noted in both (Figure 4A).GeneChip miRNA microarray of the testes revealed differential expression of 433 miRNAs, with 247 miRNAs upregulated and 186 miRNAs downregulated after exposure to TCDD (Figure 4B).Some of the highly dysregulated miRNAs following TCDD treatment have been depicted in the form of a heatmap in Figure 4C.
Ingenuity Pathway Analysis of some of the highly dysregulated miRNA showed targeting of genes involved in cell signaling pathways (Figure 5A).These dysregulated functional clusters revealed significant upregulation of certain miRs such as miR-448 and miR-let-7a as well as significant downregulation of miR-219a, and miR-15 and predicted various upstream regulators and downstream mRNA targets (Figure 5A).The predicted miRNA/mRNA interactions were aligned using TargetScan (Figure 5B).RT-qPCR of miRNAs of interest and potential mRNA targets revealed that upregulation of miRNA such as miR-448 and miR-let-7a was associated with decreased expression of specific targets such as DEAD-Box Helicase 20 (DDX20) and Sperm acrosome associated protein 6 (SPACA6), respectively.Moreover, decreased expression of miR-219a, and miR-15 was associated with increased expression of binding motif single-stranded interacting protein 1 (RBMS1), Bcl-like protein L2 (BCL2L2), and Argonaut RISC catalytic component 2 (AGO2).

Comparison of miRNA alterations between testes and ovaries and pathway analysis of miRNA to predict a targeted effect on various disorders
When we compared the miRNA expression changes between the testes and the ovaries, we noted that the TCDD-exposed testes had 247 miRNA upregulated while the ovaries had 315 miRNA that were upregulated (Figure 6A).Also, of these upregulated miRNA, 7.66% shared similar upregulation among these two organs (Figure 6A).Also, the testes showed 186 miRNA to be downregulated following TCDD exposure while the ovaries showed 111 miRNA to be Effect of TCDD on histopathology and miRNA expression in the ovaries.Adult female mice were treated with TCDD or the Vehicle and 3 days later, the ovaries were harvested for analysis as described in Methods.(A) Mouse ovaries sections stained with H&E at ×10 magnification were collected from females 72 h after treatment with either Vehicle or TCDD.(B) GeneChip miRNA 4.0 microarray differential expression analysis of miRNA isolated from mouse ovaries.miRNAs with relative expression ≥2 were considered for downstream analysis.(C) Expression heat map of miRNA that were highly dysregulated in the mouse ovaries.downregulated (Figure 6B).Additionally, these two organs showed 9.19% of such miRNA to be similarly downregulated (Figure 6B).Interestingly, only a small number (1.7%) of the dysregulated miRNA showed reciprocal alterations such as upregulation in the testes and downregulation in the ovaries (Figure 6C).Also, 8.91% of the dysregulated miRNA showed reciprocal alterations such as downregulation in the testes and upregulation in the ovaries (Figure 6D).

Discussion
TCDD or dioxin is a highly prevalent environmental contaminant.It is created by natural events such as volcanic eruptions and forest fires, and waste incineration, paper bleaching using chlorine, and certain plastics manufacture (Zhang et al., 2015).It was also found as a contaminant in Agent Orange, an herbicide used extensively during the Vietnam War which led to significant exposure of the Vietnam War Veterans as well as Vietnam civilians (Scialli et al., 2015).A chemical factory explosion near Seveso, Italy exposed the residents to high levels of TCDD (Eskenazi et al., 2018).TCDD is considered to be a human carcinogen and a potent endocrine disruptor (Bruner-Tran et al., 2017).Because it is highly lipophilic, following exposure, it can stay for a very long period in the body thereby causing chronic toxicity (Huff et al., 1991).The half-life of TCDD in humans is estimated to be around 7-9 years (Pirkle et al., 1989).Reports on human exposures and experimental studies have revealed that TCDD exposure is associated with chloracne of the skin, organ cancers, hepatotoxicity, gonadal and immune changes, pulmonary and other diseases such as diabetes, skewing of the sex ratio, and infertility (Kociba and Schwetz, 1982;Birnbaum, 1994;Kamath et al., 1998;Kamath et al., 1999;Singh et al., 2007).
The ability of TCDD to act as an endocrine disruptor suggests that it can have a significant impact on reproductive health.In fact, animal studies have revealed that exposure to TCDD in adult female mice can lead to transgenerational disorders involving reproductive functions that have been linked to endometriosis and other defects in women (Bruner-Tran et al., 2017).In certain species, exposure to TCDD is known to cause decreased fecundity and reduced ovulatory rate.Similarly, in humans, TCDD exposure was associated with decreased fertility in Seveso mothers and their daughters exposed in utero (Eskenazi et al., 2021).
The current study revealed that exposure of adult female mice to TCDD was significantly linked to poor survivorship of neonates but had no effect on the litter size at birth.Also, TCDD exposure led to an increase in the gestational period.The decreased survival of the pups was seen only when adult female but not adult male mice were exposed to TCDD.However, the increase in the gestational period was seen when either both male or female mice were exposed to TCDD.This suggests that maternal TCDD exposure may interfere with normal pup development or induces postpartum stress which can trigger cannibalism.The precise cause of postnatal mortality in the mouse is unclear and it is suggested that causes other than infanticide play a critical role (Weber et al., 2013).
miRNAs constitute noncoding RNAs that have been shown to play a critical role in the posttranscriptional regulation of gene expression.In the reproductive system, the miRNAs can regulate oocyte maturation, folliculogenesis, corpus luteum function, implantation, and early embryonic development (Eisenberg et al., 2015).It has been shown that deletion of Dicer, the enzyme that cleaves the pre-miRNA to the mature form, leads to postimplantation embryonic lethality in many animal models, thereby suggesting that miRNA plays a critical role in reproduction and development (Meng et al., 2013).
miRNA is expressed in the ovaries and is known to regulate the formation of primordial follicles, oocyte-cumulus cell interactions, granulosal cell function, and luteinization (Maalouf et al., 2016).In the current study, TCDD exposure in female mice upregulated the expression of mmu-miR-362 in the ovaries and decreased the expression of its target, ALX1, a homeobox protein-coding gene necessary for the normal development of the head and face, especially the eyes, nose, and mouth (Iyyanar et al., 2022).Poor development of these critical features may compromise neonate health, thereby prompting the mothers to cannibalize the litter.Additionally, downregulation of miR-34c suppressed the normal post-transcriptional silencing of two genes TP53I11 and TP53AIPI, which play an important role in mediating p53-dependent apoptosis.P53 is a tumor protein that acts as a tumor suppressor and is negatively regulated by proteins encoded by both TP53I11 and TP53AIPI.Both genes were upregulated in the ovaries, possibly as a response to the teratogenic and oncogenic nature of TCDD.TP53I11 and TP53AIPI have also been shown to regulate homeostasis when the cells are in a state of stress and shown to be involved in preeclampsia (Ali et al., 2021).Additionally, miR-192 was upregulated in TCDD-exposed ovaries while downstream target EREG was significantly downregulated.EREG is involved with various biological processes but may also be associated with promoting cancers in various tissues (Maalouf et al., 2016).In the reproductive system, EREG is known to stimulate oocytes (Romero and Smitz, 2009).Likewise, miR-194 was significantly upregulated in TCDD-treated females while downstream target RAP2B was suppressed.RAP2B is a protein-encoding gene in the RAS oncogene family and is involved with the proliferation and migration of various cancer types (Di et al., 2016;Di et al., 2017;Zhang et al., 2017;Li et al., 2018;Miao et al., 2019).

FIGURE 4
Effect of TCDD on histopathology and miRNA expression in the testes.Adult male mice were treated with TCDD or the Vehicle and 3 days later, the testes were harvested for analysis as described in Methods.(A) Mouse testes sections stained with H&E at ×10 magnification were collected from females 72 h after treatment with either Vehicle or TCDD.(B) GeneChip miRNA 4.0 microarray differential expression analysis of miRNA isolated from mouse testes.miRNAs with relative expression ≥2 were considered for analysis.(C) Expression heat map miRNA isolated from mouse testes that were highly dysregulated.Recent studies have shown that miRNAs play a critical role in spermatogenesis, especially during mitotic, meiotic, and postmeiotic stages of spermatogenesis (Wang and Xu, 2015).In the current study, we found that while exposure of male mice to TCDD did not alter the litter size or the survival of the pups, there was an increase in the time from pairing to parturition.Upregulation of two select miRNAs; miR-448 and mir-let-7a and the silencing of their downstream targets, DDX20 and SPACA6, respectively, may be involved in altering the testicular functions.DDX20 is a DEAD box protein-encoding gene that is involved in several cellular activities and based on distribution patterns, may be involved in spermatogenesis and embryogenesis (Lee et al., 2005).SPACA6 is a sperm acrosome membrane-associated protein directly involved in the fusion of sperm to egg plasma membrane (Barbaux et al., 2020).
Additionally, downregulation of two miRNAs; miR-219a and miR-15 in the testes by TCDD correlated with significant upregulation in their respective downstream targets RBMS1, and BCL2L2, as well as AGOT2 involved with cellular proliferation and apoptosis.BCL2L2 has additionally been shown to be critical in adult spermatogenesis and dysregulation of BCL-2 family proteins can induce germ cell apoptosis (Vergara et al., 2011).Lastly, while miR-15 targeted AGO2, it was suggested to act as an upstream regulator of miR-15, and AGO2 was also dysregulated in the ovaries treated with TCDD.This suggested an inverse relationship between AGOT and miR-15.Also, the pathway analysis of all dysregulated miRNA in the ovaries and testes revealed that these miRNA alterations affected to a greater extent, reproductive system diseases, organismal injury and abnormalities, and cancer, consistent with previous studies showing TCDD-mediated systemic toxicities.
AhR activation by TCDD leads to two distinct pathways of the regulation of gene expression: the canonical pathway and the noncanonical pathway (Wright, et al., 2017).In the canonical pathway, the AhR upon ligation, dimerizes with the AhR nuclear translocator (ARNT) and subsequently interacts with DREs leading to the induction of the genes that express DREs.In the non-canonical pathway, the AhR forms a complex with other molecules allowing it to bind to genes lacking DRE sequences leading to either induction or repression of gene expression (Viluksela and Pohjanvirta, 2019;Alhamad, et al., 2022).While the precise mechanisms through which TCDD regulates miRNA in this study is unclear, we and others have shown previously that some miRNAs express dioxin response elements (DRE) in their 3' UTR and the induction of such miRNAs occurs through DRE-dependent pathway (Neamah, et al., 2019;Singh, et al., 2020).In the current study, we investigated the miRNAs for the expression of DREs using in silico analysis and found interestingly that all the miRNA that were altered by TCDD (Figures 3, 5) such as miR-362, -34c, -192 -194, -448, -let-7a, and -15, expressed DREs.Thus, TCDD may use canonical pathway to regulate miRNA expression in the reproductive organs, although the use of non-canonical pathways cannot be ruled out.

Conclusion
The current study demonstrates for the first time that exposure of adult male and female mice to TCDD results in significant alterations in a large number of miRNAs in the ovaries and the testes.These findings suggest that environmental contaminants such as TCDD may alter reproductive functions through alterations in the expression of miRNAs.Future studies will focus on whether such miRNA changes also occur in the germ cells which would explain how TCDD exerts transgenerational effects on the reproductive system by regulating epigenetic pathways involving miRNA.

FIGURE 2
FIGURE 2Effect of TCDD on histopathology and miRNA expression in the ovaries.Adult female mice were treated with TCDD or the Vehicle and 3 days later, the ovaries were harvested for analysis as described in Methods.(A) Mouse ovaries sections stained with H&E at ×10 magnification were collected from females 72 h after treatment with either Vehicle or TCDD.(B) GeneChip miRNA 4.0 microarray differential expression analysis of miRNA isolated from mouse ovaries.miRNAs with relative expression ≥2 were considered for downstream analysis.(C) Expression heat map of miRNA that were highly dysregulated in the mouse ovaries.

FIGURE 3
FIGURE 3 Characterization of TCDD-mediated alterations in miRNA expression in the ovaries and their predicted targets.(A) Ingenuity Pathway Analysis network of specific dysregulated miRNAs, downstream targets, and upstream regulators.(B) Predicted alignments of miRNA and downstream mRNA target sequences sourced from TargetScan.(C) RT-qPCR relative quantification of miRNAs of interest and potential mRNA targets in the ovaries.In Panel C, n = 5 and the Asterix denotes as follows: * <0.05, **<0.01,***<0.001,****<0.0001.

FIGURE 5
FIGURE 5 Characterization of TCDD-mediated alterations in miRNA expression and their predicted targets in the testes.(A) Ingenuity Pathway Analysis network of specific dysregulated miRNAs, downstream targets, and upstream regulators.(B) Predicted alignments of specific miRNA and downstream mRNA target sequences sourced from TargetScan.(C) RT-qPCR relative quantification of miRNAs of interest and potential mRNA targets in the testes.In Panel C, n = 5 and the Asterix denotes as follows: * <0.05, **<0.01,***<0.001,****<0.0001.

FIGURE 6
FIGURE 6 Analysis of miRNA showing similar or dissimilar alterations in the ovaries and testes.Venn diagram showing miRNA that are upregulated (A) or downregulated (B) in ovaries and testes.The data show the numbers of miRNA changes that are unique to each organ as well as shared by the two organs.Panel C shows sharing of miRNA that were upregulated in testes and downregulated in ovaries, and vice versa (D).