Next Article in Journal
“Pseudo”-Secondary Treatment Failure Explained via Disease Progression and Effective Botulinum Toxin Therapy: A Pilot Simulation Study
Previous Article in Journal
Crotoxin Modulates Macrophage Phenotypic Reprogramming
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 15
Issue 10
10.3390/toxins15100617
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Smart-seq2 Technology Reveals a Novel Mechanism That Zearalenone Inhibits the In Vitro Maturation of Ovine Oocytes by Influencing TNFAIP6 Expression

by
Zongshuai Li
1,2,
Yali Liu
3,
Tian Ma
2,
Chen Lv
2,
Yina Li
2,
Hongwei Duan
2,
Xingxu Zhao
2,
Jianlin Wang
1,* and
Yong Zhang
2,*
1
State Key Laboratory of Grassland Agro–Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Grassland Agriculture Engineering Center, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730020, China
2
Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Gansu Agricultural University, Lanzhou 730070, China
3
Lanzhou University Second Hospital, Lanzhou 730030, China
*
Authors to whom correspondence should be addressed.
Toxins 2023, 15(10), 617; https://doi.org/10.3390/toxins15100617
Submission received: 24 September 2023 / Revised: 5 October 2023 / Accepted: 13 October 2023 / Published: 17 October 2023
(This article belongs to the Section Mycotoxins)
Browse Figures
Review Reports Versions Notes

Abstract

:
Zearalenone (ZEN), a non-steroidal estrogenic fungal toxin widely present in forage, food, and their ingredients, poses a serious threat to animal and human reproductive health. ZEN also threatens ovine, a major source of human food and breeding stock. However, the mechanisms underlying the impact of ZEN on the in vitro maturation (IVM) of ovine oocytes remain unclear. This study aimed to elucidate these mechanisms using the Smart-seq2 technology. A total of 146 differentially expressed genes were obtained, using Smart-seq2, from sheep oocytes cultured in vitro after ZEN treatment. ZEN treatment inhibited RUNX2 and SPP1 expression in the PI3K signaling pathway, leading to the downregulation of THBS1 and ultimately the downregulation of TNFAIP6; ZEN can also decrease TNFAIP6 by reducing PTPRC and ITGAM. Both inhibit in vitro maturation of ovine oocytes and proliferation of cumulus cells by downregulating TNFAIP6. These findings provide data and a theoretical basis for elucidating ZEN’s toxicity mechanisms, screening therapeutic drugs, and reducing ZEN-related losses in the ovine industry.
Keywords:
Zearalenone; ovine oocyte; Smart-seq2; TNFAIP6
Key Contribution: This research discovered a new mechanism of ZEN inhibiting in vitro maturation of ovine oocytes by downregulating the TNFAIP6 gene through Smart–seq2 technology, providing data support for elucidating ZEN’s toxicity mechanisms and screening therapeutic drugs.

1. Introduction

The non-steroidal estrogenic fungal toxin zearalenone (ZEN), also known as F-2 toxin, is mainly produced by Fusarium species such as F. graminearum [1,2]. It was first isolated from moldy corn contaminated with F. graminearum and is one of nature’s most common fungal toxins [3]. ZEN (C18H22O5) has an estrogen-like structure; therefore, it has reproductive toxicity and primarily harms female animals [4,5]. ZEN can be produced at any point during crop development, maturity, harvest, or storage. It is found in various forage, food, and raw materials, including corn, rice, soybean, rye, and wheat, with corn having the highest detection rate and contents [6]. ZEN has a melting point of 161–163 °C, and it is not particularly sensitive to environmental changes or heat treatment. Therefore, it can stably exist during the storage and processing of feed, food, and byproducts [4,7]. Ovines are a primary source of human food and a major breeding stock in the livestock industry. Their coarse feed primarily comprises corn silage, corn straw, and alfalfa, all of which are sensitive to ZEN [5,7,8]. ZEN therefore poses a significant threat to ovine farming. ZEN can remain in meat products and enter the human body, and, as with other food contaminants, can pose a risk to human reproductive health [8].
In mice, ZEN can reduce the expression of cyclin B1 and CDK1, inducing G2/M phase arrest via the IPCHK1/2–Wee1–MPF pathway, and thus interfering with oocyte maturation [9]. Moreover, in mice, single-cell sequencing revealed that ZEN altered the status of germ cell and granulosa cells via the Hippo signaling pathway, and blocked the assembly of mouse primordial follicles [10]. In pigs, ZEN can significantly increase the phosphatidylcholine or phosphatidylethanolamine content and consume hemolysin phosphatidylcholine, thereby impairing follicular cavity formation and oocyte maturation [11]. In mice and pigs, ZEN can increase oxidative stress and disrupt spindle assembly and chromosome separation, leading to premature condensation of chromatin and interfering with the expression of genes related to oocyte activity, spindle assembly, redox potential, and apoptosis, thereby inhibiting the discharge of polar bodies and cumulus expansion and affecting oocyte maturation [12,13]; it also causes autophagy, apoptosis, and oxidative stress and disrupts embryo development by affecting the expression of epigenetic histones and organelle structure and function (mitochondria, endoplasmic reticulum, Golgi apparatus, and lysosomes), in turn affecting oocyte maturation and ovarian follicle formation [14,15,16]. It has a damaging effect on preantral follicles in ovine ovaries and participates in oocyte autophagy, although the specific mechanism of toxicity is unknown [17].
Smart-seq2 sequencing is particularly suitable for studying early embryonic cell populations, making it possible to sequence limited oocyte samples without requiring the ten-times more expensive single-cell sequencing technology [18]. However, there have been no reports to date on the mechanisms underlying the impacts of ZEN on the in vitro maturation (IVM) of sheep oocytes using Smart-seq2 technology.
In summary, the toxic effects of ZEN on both animal and human reproductive health remain unclear. This research on its mechanisms of influence on sheep oocyte IVM, via Smart-seq2, is expected to reveal new pathways and genes that will improve research into its reproductive toxicity mechanisms. Referring to the functions exercised by the core pathways or genes in Smart-seq2 (such as oxidative stress, autophagy, etc.), this research is also expected to screen antagonistic drugs for ZEN, as well as provide a theoretical basis and data for reducing or even eliminating its threat to reproductive health.

2. Results

2.1. ZEN Concentration Screening

The increase in yolk gap in the three experimental groups indicates that ZEN has a significant toxic effect on sheep oocytes (Figure 1A). Based on qPCR, SOD1 and CCNB1 expression was the lowest in the 20 μmol/L treatment (Figure 1B). Although GPX1 did not show significant changes compared to the other groups, the 20 μmol/L treatment had the lowest value compared to the other two treatment groups (Figure 1B). CAT has a dose effect, but it has been significantly upregulated in the 20 μmol/L treatment (p < 0.05) (Figure 1B). CDK1 and SOD2 were both inhibited by ZEN, but their expression was highest in the 20 μmol/L treatment (Figure 1B). These results indicate that 20 μmol/L was an effective working concentration for ZEN.

2.2. Data Validation and Basic Omics Analysis

The number of oocytes used for sequencing is listed in Table 1. After obtaining omics data, GPX1, CAT, BMP15, SOD1, CDC20, GDF9, CDK1, and CCNB1 were randomly selected for qRT-PCR (Figure 1C). The qPCR and transcriptome results were consistent, indicating that the transcriptome data were accurate and reliable (Figure 1D). Based on the correlation heat map (Figure 2A), the correlation coefficients between samples were ≥0.70. The violin plot of the sample shows that there is no significant degradation outlier in the sample, and the obtained data are valid and reliable (Figure 2B). Screening via edgeR (criteria: q Value ≤ 0.1 and multiple difference ≥ 1.2) identified 102 upregulated genes and 44 downregulated genes; these were selected for subsequent analysis (Figure 2C).

2.3. Enrichment Pathway Analyses

GO analysis was performed on the differentially expressed genes (DEGs). The 20 most enriched GO terms belonged mainly to biological processes and were related to intercellular interactions, cell development, and immune responses (Figure 2D). KEGG analysis of DEGs revealed that ZEN mainly affected signaling pathways related to human diseases and environmental information processes such as MicroRNAs in cancer, the PI3K-Akt and JAK-STAT signaling pathways, systemic lupus erythematosus, the TGF-beta signaling pathway, and extracellular matrix-receptor interactions (Figure 2E).
GSEA analysis was conducted on the sequencing results: the P450, autophagy, and cancer-related pathways were activated (Figure 3A), whereas signaling pathways related to primary immunodeficiency, autoimmune deficiency, and various diseases were suppressed (Figure 3B). ITGAM and MHC2 were suppressed, indicating that ZEN inhibited the ovine oocyte immune system.

2.4. Key Gene Mining

Genes related to oocyte maturation, cumulus expansion factor, oxidative stress, apoptosis, autophagy, and spindle assembly detection points in the transcriptome data were extracted and visualized as histograms (Figure 4A,C–F). The expression of most genes related to oxidative stress, oocyte apoptosis, and autophagy was promoted, while that of most genes related to spindle, and all (tested) genes related to oocyte nuclear maturation, and cumulus expansion factor was suppressed. Fluorescence staining results showed that ZEN treatment promoted the production and accumulation of ROS in sheep oocytes cultured in vitro (Figure 4B).
Using String software (https://string-db.org/, accessed on 12 October 2023) to analyze the 146 DEGs, a network diagram of their protein–protein interactions (PPIs) was obtained, with ovine as the reference species. A total of 68 predicted proteins was obtained (Figure 5A). By labeling gene expression, we discovered the distribution of the downregulated and upregulated proteins (Figure 5A). Through Markov clustering in String, we obtained four sets with proteins number ≥ 4 (a gene set with a proteins number of 15 was suppressed) (Figure 5B).
Through a literature review, we classified 68 genes into five main categories: cancer (Figure 5C), oocyte development (Figure 5D), proliferation and aging (Figure 5E), oxidative stress (Figure 5F), and immunity (Figure 5G). Their respective heat maps indicated that except for a large number of cancer-related genes being activated, the other four categories of genes had varying degrees of activation and inhibition.

2.5. Diagram of the Mechanism of the Effects of ZEN on Ovine Oocyte IVM

We constructed a PPI network from transcriptome data, including DEG, estrogen receptors, and six related gene sets. The results showed that DEGs were most closely related to oocyte maturation, cumulus expansion factor, and cell apoptosis (Figure 6A,C). After removing the gene sets of oxidative stress, spindle assembly detection points, and autophagy from the network, we screened proteins with connection numbers ≥3 to obtain a closely related proteins set network (Figure 6C) and conduct a joint analysis with the directed network of DEG corresponding proteins (Figure 6B).
Finally, we inferred a conceptual diagram of the mechanism whereby ZEN affects ovine oocyte IVM (Figure 7).

3. Discussion

ZEN toxicity is known to be dose-dependent [12]. This was also confirmed early in the experiment; however, the maturation rates of the 20 and 30 μmol/L groups were very similar. qPCR detection showed that genes related to oxidative stress and cell cycle obstruction had peaks in the 20 μmol/L treatment group. We speculated that when the concentration of ZEN was 30 μmol/L, oocytes activated a limited correction mechanism for survival, so it was no longer conducive to research. Therefore, we selected the 20 μmol/L treatment group with a peak for research.
The samples are prepared and transported in batches. There are batch differences between the samples; however, the sample repeatability is strongly correlated. GO analysis revealed that cell adhesion, cell development, and immune-related terms were mainly affected; KEGG analysis revealed that the PI3K signaling pathway was at the core position; GSEA revealed that various immune mechanisms, such as primary immunity, were suppressed, primarily due to downregulation of ITGAM and MHC. These results indicate that ZEN treatment mainly affects ovine oocyte IVM by affecting cell development, cancer, and immunity; however, the specific mechanisms require further analysis.
Previous studies have shown that ZEN affects animal oocyte maturation mainly via the expression of genes involved in oocyte maturation, cumulus expansion factor, oxidative stress, cell apoptosis, autophagy, and spindle assembly checkpoint [9,10,11,12,13,14,15,16,17]. The overall transcriptome data analysis revealed that ZEN induced oxidative stress, promoted apoptosis and autophagy, and inhibited oocyte development and spindle assembly, thus impeding ovine oocyte IVM. This is consistent with existing research results; hence, we further classified DEG. The results showed that the genes are mainly enriched in the PI3K signaling pathway and have an impact on IVM of ovine oocytes through the following five aspects: (1) Oocyte maturation, the overexpression of genes such as CCND2, PDGFC, IFRD1, H2AFZ, INHBA, CD200, and RPS26 indicated that both groups of oocytes had undergone IVM [19,20,21,22,23,24,25]. However, in the experimental group, the expression levels of TNFAIP6 related to cumulus cell proliferation [13], SPP1 related to in vitro development of oocytes [26], RUNX2 related to follicular development [27], and C1QB related to embryonic development were downregulated [28], indicating that ZEN may inhibit the process of oocyte maturation by influencing these genes. (2) Overexpression of cancer genes, such as FAM20C, PDGFRB, LAMB2, PSME2, CAD, ANXA3, MBNL1, EWSR1, RPL23A, NPNT, and SERPINH1, is consistent with the reports that ZEN promotes cancer occurrence, which will lead to disruption of the oocyte metabolic system [29,30,31,32,33,34,35,36,37,38,39]. (3) Cell proliferation and aging, as well as the expression of genes such as CCN1, A1CF, IL7R, DHRS3, and FBLN1, indicated that ZEN promotes cell aging and apoptosis [40,41,42,43,44]. (4) Oxidative stress, overexpression of CPE and NDRG1, indicated the generation of ROS in cells [45,46]. Downregulation of DAPP1 and APOD leads to a decrease in ROS production and antioxidant capacity [47,48], ultimately leading to an increase in ROS in oocytes. (5) Immunity, downregulation of LCP1, LRRN3, and ERAP2 indicated that ZEN can cause damage to the immune system, leading to cell apoptosis [49,50,51].
Next, we further analyzed the relationship among estrogen receptor genes, enriched differential genes, and the six selected gene sets. We found that for IVM of ovine oocytes, the DEGs enriched in the ZEN treatment group were mainly related to estrogen receptors, cumulus expansion factors, oocyte maturation, and apoptosis-related genes. Although ZEN treatment caused oxidative stress, promoted autophagy, and interfered with spindle assembly, these processes are closely related to apoptosis-related genes but not to DEGs. The results of the directed network graph indicate that in ovines, ZEN–estrogen receptor binding inhibits the expression of RUNX2 and SPP1 in the PI3K signaling pathway, thereby downregulating THBS1 and ultimately TNFAIP6, which in turn inhibits ovine oocyte IVM and cumulus cell proliferation. Furthermore, reducing PTPRC and ITGAM expression downregulates the expression of genes related to PI3K signaling, and ultimately TNFAIP6, thereby promoting apoptosis and inhibiting cumulus cell proliferation.

4. Conclusions

This study used Smart-seq2 single-cell transcriptome technology to elucidate the mechanism in which ZEN affects ovine oocyte IVM. We hypothesize two mechanisms for its action. First, inhibiting the expression of RUNX2 and SPP1 in the PI3K signaling pathway causes the downregulation of THBS1, and ultimately that of TNFAIP6, thereby inhibiting ovine oocyte IVM and cumulus cell proliferation. Second, reducing PTPRC and ITGAM expression downregulates the expression of genes related to PI3K signaling, and ultimately TNFAIP6, thereby promoting apoptosis and inhibiting cumulus cell proliferation. Further, these findings reveal that ZEN can cause autophagy-related gene overexpression and ROS accumulation in ovine oocytes, although further research is needed to clarify the related mechanisms. This study provides data and a theoretical basis for improving the study of the mechanisms of ZEN toxicity, screening therapeutic drugs, and reducing ZEN-related losses in the ovine industry.

5. Materials and Methods

5.1. In Vitro Culture of Ovine Oocytes

Washing solution (40 mL) was prepared using 39 mL M199 (Biological Industries, Kibbutz Beit Haemek, Israel), 1 mL HEPES (1 M) (Sigma Aldrich, St. Louis, MO, USA), and 160 mg of bovine serum albumin (BSA; Sigma Aldrich). Then, 10 mL of maturation solution was prepared using 10 mL M199 (Biological Industries), 10 µL follicle stimulating hormone (FSH, 5 μg/mL; Sigma Aldrich), 10 µL luteinizing hormone (LH, 5 μg/mL; Sigma Aldrich), 10 µL estradiol (1 μg/mL; Sigma Aldrich), 100 µL double antibody (10,000 units/mL Penicillin G Sodium Salt, 10 mg/mL Streptomycin Sulfate; Biological Industries), and 60 mg BSA (Sigma Aldrich). Both culture media were disinfected using the 0.22 μm filter.
The sheep ovaries were removed immediately after slaughter and placed in sterile phosphate-buffered saline (PBS) containing double antibody (final concentration 10.0 units/mL Penicillin G Sodium Salt, 100 µg/mL Streptomycin Sulfate) at 37 °C before being transported back to the laboratory within 2–3 h. The ovaries were placed in a sterile beaker, rinsed 2–3 times with PBS, and immersed in PBS containing double antibodies (final concentration 10.0 units/mL Penicillin G Sodium Salt, 100 µg/mL Streptomycin Sulfate). The beaker with ovaries was placed in a 37 °C water bath for later use. A 5 mL syringe was used to extract the follicular fluid, which was injected into a centrifuge tube containing 10 mL wash solution. At the same time, a 35 mm cell culture dish was prepared on an ultra-clean table, a 3 mm × 3 mm square was drawn on the underside of the dish, washing solution was added, and the dish was placed in the incubator for later use.
The centrifuge tube containing the follicular fluid was rested for 5 min, disinfected, and placed on an ultra-clean table. The lower layer of liquid was gently removed and transferred into a spare 35 mm culture dish. Under a microscope, cumulus–oocyte complexes (COCs) were extracted and placed in a new 35 mm culture dish. The COCs were then moved to another new 35 mm culture dish. Finally, the COCs were moved into a four-well plate (Thermo Fisher Scientific, Waltham, MA, USA) with droplets of maturation solution, and transferred to an incubator (5% CO2, 37 °C) for cultivation. The criteria for selecting COCs were as follows: oocyte cytoplasm uniformly filled within the zona pellucida, with three or more layers of dense cumulus cells (Grade A), and one to three layers of cumulus cells (Grade B) surrounding the COCs.

5.2. Screening of ZEN Working Concentration

The control and experimental groups were designated as CK (no ZEN) and T1, T2, and T3 (10, 20, and 30 μmol/L, respectively, in the maturation solution). After 26 h, hyaluronidase was used to remove granulosa cells. The naked eggs were transferred into a centrifuge tube containing PBS, centrifuged at 106 g for 5 min, and the supernatant was discarded. Next, 350 µL of TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, USA) was added, total RNA was extracted from cells according to the manufacturer’s protocol, and the final sample was stored at −80 °C. cDNA was prepared using the RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA) following the manufacturer’s protocol. A QuantiNova SYBR Green PCR Kit (Qiagen, Dusseldorf, Germany) was used for quantitative real-time PCR (qRT-PCR) assays.
ZEN can cause oxidative stress and cell cycle obstruction in oocytes. Therefore, qPCR was performed to quantify SOD1, SOD2, CAT, GPX, CDK1, and CCNB1 expression (Table 1). Working concentrations were determined based on the expression of the six genes at various concentrations.

5.3. ZEN Sample Preparation and Delivery

We established control and experimental groups using the optimal working ZEN concentration for the experimental groups. Because a limited number of oocytes was obtained each time, the samples were prepared six times, with the control and experimental groups set up each time, and any two sets of replicated samples were randomly combined. Finally, the control and experimental groups each had three biological replicates, for a total of six samples (Table 2). The cell samples were sent to Genedenovo Biotechnology Co., Ltd. (Guangzhou, China) on dry ice to establish a cDNA library.

5.4. Immunofluorescence

The reactive oxygen specifications (ROS) detection assay kit (BioVision, San Francisco Bay, CA, USA) was used for detection. Cumulus cells were digested, and the naked eggs were rinsed twice with PBS before being collected in centrifuge tubes. A 100 μL reactive oxygen species (ROS) assay buffer was added, and the cells were gently blown and centrifuged at 106× g for 5 min. The supernatant was discarded and rinsed with PBS once. A 100 uL 1× ROS Label Buffer was added, and the resulting mixture was incubated at 37 °C in the dark for 45 min and centrifuged at 106× g for 5 min. The supernatant was discarded, and 100 μL ROS assay buffer was added. The cells were gently blown and centrifuged at 106× g for 5 min to discard the supernatant. A 100 uL ROS Inducer Buffer was added, the mixture was incubated in the dark for 1 h, and the results were observed under an inverted fluorescence microscope.

5.5. Smart-seq2 Data Validation, Visualization, and Analysis

Eight genes (GPX1, CAT, BMP15, SOD1, CDC20, GDF9, CDK1, and CCNB1) were randomly selected to validate the sequencing results (Tabel 1). After verifying the sequencing results, the enriched Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were subjected to omics analysis, supplemented with Gene Set Enrichment Analysis (GSEA), in order to select the important signaling pathways and genes and to make preliminary inferences about the mechanisms of action of ZEN.

Author Contributions

Z.L.: Conceptualization, formal analysis, investigation, methodology, validation, visualization, writing—original draft; Y.L. (Yali Liu) and T.M.: visualization; C.L., Y.L. (Yina Li) and H.D.: methodology; X.Z.: writing—review and editing; J.W. and Y.Z.: conceptualization, funding acquisition, project administration, resources. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the national natural science fund joint fund key project: National Key R&D Program of China (2021YFD1300902) and 13th Five-Year National Key Research and development plan food safety technology research and development major project (2019YFC1605705).

Institutional Review Board Statement

This work was performed with the Animal Care and Use Committee of the College of Veterinary Medicine at Gansu Agricultural University Approval GSAU-Eth-LST-2021-003.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors have no competing interests to declare.

References

  1. Al-Jaal, B.A.; Jaganjac, M.; Barcaru, A.; Horvatovich, P.; Latiff, A. Aflatoxin, fumonisin, ochratoxin, zearalenone and deoxynivalenol biomarkers in human biological fluids: A systematic literature review, 2001–2018. Food Chem. Toxicol. 2019, 129, 211–228. [Google Scholar] [CrossRef] [PubMed]
  2. Caglayan, M.O.; Şahin, S.; Üstündağ, Z. Detection Strategies of Zearalenone for Food Safety: A Review. Crit. Rev. Anal. Chem. 2022, 52, 294–313. [Google Scholar] [CrossRef]
  3. Stob, M.; Baldwin, R.S.; Tuite, J.; Andrews, F.N.; Gillette, K.G. Isolation of an anabolic, uterotrophic compound from corn infected with Gibberella zeae. Nature 1962, 196, 1318. [Google Scholar] [CrossRef] [PubMed]
  4. Han, X.; Huangfu, B.; Xu, T.; Xu, W.; Asakiya, C.; Huang, K.; He, X. Research Progress of Safety of Zearalenone: A Review. Toxins 2022, 14, 386. [Google Scholar] [CrossRef] [PubMed]
  5. Ropejko, K.; Twarużek, M. Zearalenone and Its Metabolites-General Overview, Occurrence, and Toxicity. Toxins 2021, 13, 35. [Google Scholar] [CrossRef] [PubMed]
  6. Golge, O.; Bulent, K. Occurrence of deoxynivalenol and zearalenone in cereals and cereal products from Turkey. Food Control 2020, 110, 106982. [Google Scholar] [CrossRef]
  7. EFSA Panel on Contaminants in the Food Chain (CONTAM); Knutsen, H.K.; Alexander, J.; Barregård, L.; Bignami, M.; Brüschweiler, B.; Ceccatelli, S.; Cottrill, B.; Dinovi, M.; Edler, L.; et al. Risks for animal health related to the presence of zearalenone and its modified forms in feed. EFSA J. 2017, 15, e04851. [Google Scholar]
  8. Liu, J.; Applegate, T. Zearalenone (ZEN) in Livestock and Poultry: Dose, Toxicokinetics, Toxicity and Estrogenicity. Toxins 2020, 12, 377. [Google Scholar] [CrossRef]
  9. Ji, Y.M.; Zhang, K.H.; Pan, Z.N.; Ju, J.Q.; Zhang, H.L.; Liu, J.C.; Wang, Y.; Sun, S.C. High-dose zearalenone exposure disturbs G2/M transition during mouse oocyte maturation. Reprod. Toxicol. 2022, 110, 172–179. [Google Scholar] [CrossRef]
  10. Tian, Y.; Zhang, M.Y.; Zhao, A.H.; Kong, L.; Wang, J.J.; Shen, W.; Li, L. Single-cell transcriptomic profiling provides insights into the toxic effects of Zearalenone exposure on primordial follicle assembly. Theranostics 2021, 11, 5197–5213. [Google Scholar] [CrossRef]
  11. Lai, F.N.; Liu, X.L.; Li, N.; Zhang, R.Q.; Zhao, Y.; Feng, Y.Z.; Nyachoti, C.M.; Shen, W.; Li, L. Phosphatidylcholine could protect the defect of zearalenone exposure on follicular development and oocyte maturation. Aging 2018, 10, 3486–3506. [Google Scholar] [CrossRef] [PubMed]
  12. Pan, L.Z.; Cheng, H.; Zhang, J.; Gong, S.; Tian, X.D.; Pan, C.J.; Luo, M.J.; Tan, J.H. Invivo zearalenone exposure dose-dependently compromises mouse oocyte competence by impairing chromatin configuration and gene transcription. Reprod. Fertil. Dev. 2021, 33, 229–238. [Google Scholar] [CrossRef] [PubMed]
  13. Lai, F.N.; Ma, J.Y.; Liu, J.C.; Wang, J.J.; Cheng, S.F.; Sun, X.F.; Li, L.; Li, B.; Nyachoti, C.M.; Shen, W. The influence of N-acetyl-l-cysteine on damage of porcine oocyte exposed to zearalenone in vitro. Toxicol. Appl. Pharmacol. 2015, 289, 341–348. [Google Scholar] [CrossRef] [PubMed]
  14. Feng, Y.Q.; Wang, J.J.; Li, M.H.; Tian, Y.; Zhao, A.H.; Li, L.; De Felici, M.; Shen, W. Impaired primordial follicle assembly in offspring ovaries from zearalenone-exposed mothers involves reduced mitochondrial activity and altered epigenetics in oocytes. Cell Mol. Life Sci. 2022, 79, 258. [Google Scholar] [CrossRef] [PubMed]
  15. Xu, J.; Sun, L.; He, M.; Zhang, S.; Gao, J.; Wu, C.; Zhang, D.; Dai, J. Resveratrol Protects against Zearalenone-Induced Mitochondrial Defects during Porcine Oocyte Maturation via PINK1/Parkin-Mediated Mitophagy. Toxins 2022, 14, 641. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, Y.; Xing, C.H.; Chen, S.; Sun, S.C. Zearalenone exposure impairs organelle function during porcine oocyte meiotic maturation. Theriogenology 2022, 177, 22–28. [Google Scholar] [CrossRef] [PubMed]
  17. Silva, I.P.; Brito, D.C.C.; Silva, T.E.S.; Silva, R.F.; Guedes, M.I.F.; Silva, J.Y.G.; Rodrigues, A.P.R.; Santos, R.R.; Figueiredo, J.R. In vitro exposure of sheep ovarian tissue to the xenoestrogens zearalenone and enterolactone: Effects on preantral follicles. Theriogenology 2021, 174, 124–130. [Google Scholar] [CrossRef] [PubMed]
  18. Björklund, Å.K.; Forkel, M.; Picelli, S.; Konya, V.; Theorell, J.; Friberg, D.; Sandberg, R.; Mjösberg, J. The heterogeneity of human CD127(+) innate lymphoid cells revealed by single-cell RNA sequencing. Nat. Immunol. 2016, 17, 451–460. [Google Scholar] [CrossRef]
  19. Chermuła, B.; Jeseta, M.; Sujka-Kordowska, P.; Konwerska, A.; Jankowski, M.; Kranc, W.; Kocherova, I.; Celichowski, P.; Antosik, P.; Bukowska, D.; et al. Genes regulating hormone stimulus and response to protein signaling revealed differential expression pattern during porcine oocyte in vitro maturation, confirmed by lipid concentration. Histochem. Cell Biol. 2020, 154, 77–95. [Google Scholar] [CrossRef]
  20. Wang, Q.; Guo, X.; Yao, D.; Wang, B.; Li, Y.; Zhang, J.; Zhang, X. Insights into Transcriptomic Differences in Ovaries between Lambs and Adult Sheep after Superovulation Treatment. Animals 2023, 13, 665. [Google Scholar] [CrossRef]
  21. Li, F.X.; Liu, Y.; Miao, X.P.; Fu, G.Q.; Curry, T.E., Jr. Expression and regulation of the differentiation regulators ERBB Receptor Feedback Inhibitor 1 (ERRFI1) and Interferon-related Developmental Regulator 1 (IFRD1) during the periovulatory period in the rat ovary. Mol. Reprod. Dev. 2016, 83, 714–723. [Google Scholar] [CrossRef] [PubMed]
  22. Chediek Dall’Acqua, P.; Barros Nunes, G.; Rodrigues da Silva, C.; Fontes, P.K.; Fábio Gouveia Nogueira, M.; Lombardi Lopes, F.; Marinho, M.; Zoccal Mingoti, G. Differences in embryonic gene expression and quality indicate the benefit of epidermal growth factor receptor inhibitor during prematuration to improve competence in bovine oocytes. Reprod. Domest. Anim. 2019, 54, 666–677. [Google Scholar] [CrossRef] [PubMed]
  23. Brązert, M.; Kranc, W.; Nawrocki, M.J.; Sujka-Kordowska, P.; Konwerska, A.; Jankowski, M.; Kocherova, I.; Celichowski, P.; Jeseta, M.; Ożegowska, K.; et al. New markers for regulation of transcription and macromolecule metabolic process in porcine oocytes during in vitro maturation. Mol. Med. Rep. 2020, 21, 1537–1551. [Google Scholar] [CrossRef]
  24. Assou, S.; Anahory, T.; Pantesco, V.; Le Carrour, T.; Pellestor, F.; Klein, B.; Reyftmann, L.; Dechaud, H.; De Vos, J.; Hamamah, S. The human cumulus--oocyte complex gene-expression profile. Hum. Reprod. 2006, 21, 1705–1719. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, X.M.; Yan, M.Q.; Ji, S.Y.; Sha, Q.Q.; Huang, T.; Zhao, H.; Liu, H.B.; Fan, H.Y.; Chen, Z.J. Loss of oocyte Rps26 in mice arrests oocyte growth and causes premature ovarian failure. Cell Death Dis. 2018, 9, 1144. [Google Scholar] [CrossRef] [PubMed]
  26. Munakata, Y.; Kawahara-Miki, R.; Shiratsuki, S.; Tasaki, H.; Itami, N.; Shirasuna, K.; Kuwayama, T.; Iwata, H. Gene expression patterns in granulosa cells and oocytes at various stages of follicle development as well as in in vitro grown oocyte-and-granulosa cell complexes. J. Reprod. Dev. 2016, 62, 359–366. [Google Scholar] [CrossRef] [PubMed]
  27. Gao, K.; Wang, P.; Peng, J.; Xue, J.; Chen, K.; Song, Y.; Wang, J.; Li, G.; An, X.; Cao, B. Regulation and function of runt-related transcription factors (RUNX1 and RUNX2) in goat granulosa cells. J. Steroid Biochem. Mol. Biol. 2018, 181, 98–108. [Google Scholar] [CrossRef]
  28. Ortega, M.S.; Kurian, J.J.; McKenna, R.; Hansen, P.J. Characteristics of candidate genes associated with embryonic development in the cow: Evidence for a role for WBP1 in development to the blastocyst stage. PLoS ONE 2017, 12, e0178041. [Google Scholar] [CrossRef]
  29. Liu, X.; Zhan, Y.; Xu, W.; Liu, X.; Geng, Y.; Liu, L.; Da, J.; Wang, J.; Zhang, X.; Jin, H.; et al. Prognostic and immunological role of Fam20C in pan-cancer. Biosci. Rep. 2021, 41, BSR20201920. [Google Scholar] [CrossRef]
  30. Sun, T.; Bi, F.; Liu, Z.; Yang, Q. TMEM119 facilitates ovarian cancer cell proliferation, invasion, and migration via the PDGFRB/PI3K/AKT signaling pathway. J. Transl. Med. 2021, 19, 111. [Google Scholar] [CrossRef]
  31. Pingili, A.K.; Chaib, M.; Sipe, L.M.; Miller, E.J.; Teng, B.; Sharma, R.; Yarbro, J.R.; Asemota, S.; Al Abdallah, Q.; Mims, T.S.; et al. Immune checkpoint blockade reprograms systemic immune landscape and tumor microenvironment in obesity-associated breast cancer. Cell Rep. 2021, 35, 109285. [Google Scholar] [CrossRef] [PubMed]
  32. Wang, X.; Wu, F.; Deng, Y.; Chai, J.; Zhang, Y.; He, G.; Li, X. Increased expression of PSME2 is associated with clear cell renal cell carcinoma invasion by regulating BNIP3-mediated autophagy. Int. J. Oncol. 2021, 59, 106. [Google Scholar] [CrossRef] [PubMed]
  33. Li, G.; Li, D.; Wang, T.; He, S. Pyrimidine Biosynthetic Enzyme CAD: Its Function, Regulation, and Diagnostic Potential. Int. J. Mol. Sci. 2021, 22, 10253. [Google Scholar] [CrossRef] [PubMed]
  34. Kim, J.Y.; Jung, E.J.; Kim, J.M.; Son, Y.; Lee, H.S.; Kwag, S.J.; Park, J.H.; Cho, J.K.; Kim, H.G.; Park, T.; et al. MiR-221 and miR-222 regulate cell cycle progression and affect chemosensitivity in breast cancer by targeting ANXA3. Exp. Ther. Med. 2023, 25, 127. [Google Scholar] [CrossRef] [PubMed]
  35. Zhang, Q.; Wu, Y.; Chen, J.; Tan, F.; Mou, J.; Du, Z.; Cai, Y.; Wang, B.; Yuan, C. The Regulatory Role of Both MBNL1 and MBNL1-AS1 in Several Common Cancers. Curr. Pharm. Des. 2022, 28, 581–585. [Google Scholar] [CrossRef] [PubMed]
  36. Lee, J.; Nguyen, P.T.; Shim, H.S.; Hyeon, S.J.; Im, H.; Choi, M.H.; Chung, S.; Kowall, N.W.; Lee, S.B.; Ryu, H. EWSR1, a multifunctional protein, regulates cellular function and aging via genetic and epigenetic pathways. Biochim. Biophys. Acta Mol. Basis Dis. 2019, 1865, 1938–1945. [Google Scholar] [CrossRef] [PubMed]
  37. Zhang, Z.; Wu, Q.; Fang, M.; Liu, Y.; Jiang, J.; Feng, Q.; Hu, R.; Xu, J. HERC3 directly targets RPL23A for ubiquitination degradation and further regulates Colorectal Cancer proliferation and the cell cycle. Int. J. Biol. Sci. 2022, 18, 3282–3297. [Google Scholar] [CrossRef] [PubMed]
  38. Huang, L.; Guan, S.; Feng, L.; Wei, J.; Wu, L. Integrated analysis identified NPNT as a potential key regulator in tumor metastasis of hepatocellular carcinoma. Gene 2022, 825, 146436. [Google Scholar] [CrossRef]
  39. Tian, S.; Peng, P.; Li, J.; Deng, H.; Zhan, N.; Zeng, Z.; Dong, W. SERPINH1 regulates EMT and gastric cancer metastasis via the Wnt/β-catenin signaling pathway. Aging 2020, 12, 3574–3593. [Google Scholar] [CrossRef]
  40. Jun, J.I.; Lau, L.F. The matricellular protein CCN1 induces fibroblast senescence and restricts fibrosis in cutaneous wound healing. Nat. Cell Biol. 2010, 12, 676–685. [Google Scholar] [CrossRef]
  41. Ni, D.; Yi, Q.; Liu, J.; Hu, Y.; Lv, T.; Tan, G.; Liu, Y.; Xu, L.; Xia, H.; Zhou, Q.; et al. A1CF-promoted colony formation and proliferation of RCC depends on DKK1-MEK/ERK signal axis. Gene 2020, 730, 144299. [Google Scholar] [CrossRef] [PubMed]
  42. Lee, K.Y.; Ho, S.C.; Sun, W.L.; Feng, P.H.; Lin, C.W.; Chen, K.Y.; Chuang, H.C.; Tseng, C.H.; Chen, T.T.; Wu, S.M. Lnc-IL7R alleviates PM2.5-mediated cellular senescence and apoptosis through EZH2 recruitment in chronic obstructive pulmonary disease. Cell Biol. Toxicol. 2022, 38, 1097–1120. [Google Scholar] [CrossRef] [PubMed]
  43. Sumei, S.; Xiangyun, K.; Fenrong, C.; Xueguang, S.; Sijun, H.; Bin, B.; Xiaolei, S.; Yongjiu, T.; Kaichun, W.; Qingchuan, Z.; et al. Hypermethylation of DHRS3 as a Novel Tumor Suppressor Involved in Tumor Growth and Prognosis in Gastric Cancer. Front. Cell Dev. Biol. 2021, 9, 624871. [Google Scholar] [CrossRef] [PubMed]
  44. Xu, G.; Geng, X.; Yang, F.; Zhang, H. FBLN1 promotes chondrocyte proliferation by increasing phosphorylation of Smad2. J. Orthop. Sci. 2022, 27, 242–248. [Google Scholar] [CrossRef]
  45. Li, J.; Dong, Z.; Pan, Y.; Wang, L.; Zhao, W.; Zhang, J. CPE Regulates Proliferation and Apoptosis of Primary Myocardial Cells Mediated by Ischemia and Hypoxia Injury. J. Healthc. Eng. 2022, 2022, 3155171. [Google Scholar] [CrossRef] [PubMed]
  46. Nishigaki, A.; Tsubokura, H.; Ishida, M.; Hashimoto, Y.; Yoshida, A.; Hisamatsu, Y.; Tsuzuki-Nakao, T.; Murata, H.; Okada, H. NDRG1 is expressed in human granulosa cells: An implicative role of NDRG1 in the ovary. Reprod. Med. Biol. 2022, 21, e12437. [Google Scholar] [CrossRef] [PubMed]
  47. Hao, L.; Marshall, A.J.; Liu, L. Bam32/DAPP1-Dependent Neutrophil Reactive Oxygen Species in WKYMVm-Induced Microvascular Hyperpermeability. Front. Immunol. 2020, 11, 1028. [Google Scholar] [CrossRef] [PubMed]
  48. Zhang, Y.; Cong, Y.; Wang, S.; Zhang, S. Antioxidant activities of recombinant amphioxus (Branchiostoma belcheri) apolipoprotein D. Mol. Biol. Rep. 2011, 38, 1847–1851. [Google Scholar] [CrossRef]
  49. Linehan, J.B.; Lucas Zepeda, J.; Mitchell, T.A.; LeClair, E.E. Follow that cell: Leukocyte migration in L-plastin mutant zebrafish. Cytoskeleton 2022, 79, 26–37. [Google Scholar] [CrossRef]
  50. Chou, J.P.; Ramirez, C.M.; Wu, J.E.; Effros, R.B. Accelerated aging in HIV/AIDS: Novel biomarkers of senescent human CD8+ T cells. PLoS ONE 2013, 8, e64702. [Google Scholar] [CrossRef]
  51. Yang, Z.; Tian, H.; Bie, F.; Xu, J.; Zhou, Z.; Yang, J.; Li, R.; Peng, Y.; Bai, G.; Tian, Y.; et al. ERAP2 Is Associated With Immune Infiltration and Predicts Favorable Prognosis in SqCLC. Front. Immunol. 2021, 12, 788985. [Google Scholar] [CrossRef]
Toxins 15 00617 g001
Figure 1. ZEN concentration screening and transcriptome data validation. (A) The growth status of sheep oocytes cultured in vitro after 36 h of exposure to different concentrations of ZEN(CK,10 μM, 20 μM, and 30 μM) (10×, 100 μm), (B) qRT-PCR results of SOD1, CCNB1, GPX1, CAT, CDK1, and SOD2 in the four groups, (C) qRT-PCR results of randomly selected GPX1, CAT, BMP15, SOD1, CDC20, GDF9, CDK1, and CCNB1, (D) histogram of log2 values for GPX1, CAT, BMP15, SOD1, CDC20, GDF9, CDK1, and CCNB in transcriptome data.
Figure 1. ZEN concentration screening and transcriptome data validation. (A) The growth status of sheep oocytes cultured in vitro after 36 h of exposure to different concentrations of ZEN(CK,10 μM, 20 μM, and 30 μM) (10×, 100 μm), (B) qRT-PCR results of SOD1, CCNB1, GPX1, CAT, CDK1, and SOD2 in the four groups, (C) qRT-PCR results of randomly selected GPX1, CAT, BMP15, SOD1, CDC20, GDF9, CDK1, and CCNB1, (D) histogram of log2 values for GPX1, CAT, BMP15, SOD1, CDC20, GDF9, CDK1, and CCNB in transcriptome data.
Toxins 15 00617 g001
Toxins 15 00617 g002
Figure 2. Visualization of transcriptome data. (A) Correlation heat map, (B) violin diagrams of six samples, (C) differentially expressed genes, (D) the 20 most enriched GO terms, and (E) KEGG analysis of DEGs.
Figure 2. Visualization of transcriptome data. (A) Correlation heat map, (B) violin diagrams of six samples, (C) differentially expressed genes, (D) the 20 most enriched GO terms, and (E) KEGG analysis of DEGs.
Toxins 15 00617 g002
Toxins 15 00617 g003
Figure 3. GSEA analysis. (A) Signal pathways upregulated by genes in GSEA analysis, and (B) signal pathways downregulated by genes in GSEA analysis.
Figure 3. GSEA analysis. (A) Signal pathways upregulated by genes in GSEA analysis, and (B) signal pathways downregulated by genes in GSEA analysis.
Toxins 15 00617 g003
Toxins 15 00617 g004
Figure 4. Visualization of transcriptome data. (A) Histogram of log2 values of oxidative stress related genes in transcriptome data, (B) fluorescence staining of reactive oxygen species in the control and experimental groups (fluorescent and white light) (20×, 50 μm), (C) histogram of log2 values of apoptosis related genes in transcriptome data, (D) histogram of log2 values of autophagy-related genes in transcriptome data, (E) histogram of log2 values of oocyte maturation and cumulus expansion factor related genes in transcriptome data, and (F) histogram of log2 values of spindle assembly detection point related genes in transcriptome data.
Figure 4. Visualization of transcriptome data. (A) Histogram of log2 values of oxidative stress related genes in transcriptome data, (B) fluorescence staining of reactive oxygen species in the control and experimental groups (fluorescent and white light) (20×, 50 μm), (C) histogram of log2 values of apoptosis related genes in transcriptome data, (D) histogram of log2 values of autophagy-related genes in transcriptome data, (E) histogram of log2 values of oocyte maturation and cumulus expansion factor related genes in transcriptome data, and (F) histogram of log2 values of spindle assembly detection point related genes in transcriptome data.
Toxins 15 00617 g004
Toxins 15 00617 g005
Figure 5. Gene classification analysis in DEGs. (A) PPI network diagram of DEGs, (B) four gene sets with gene number ≥4, (C) heat map of cancer related genes, (D) heat map of oocyte development related genes, (E) heat map of proliferation and aging related genes, (F) heat map of oxidative stress related genes, and (G) heat map of immunity related genes.
Figure 5. Gene classification analysis in DEGs. (A) PPI network diagram of DEGs, (B) four gene sets with gene number ≥4, (C) heat map of cancer related genes, (D) heat map of oocyte development related genes, (E) heat map of proliferation and aging related genes, (F) heat map of oxidative stress related genes, and (G) heat map of immunity related genes.
Toxins 15 00617 g005
Toxins 15 00617 g006
Figure 6. PPI analysis of DEGs. (A) PPI analysis of DEGs (within the red box) and its relationship with seven other selected genes, (B) the directed network graph of DEGs, (C) relationship diagram of proteins with connection number ≥3 in PPI.
Figure 6. PPI analysis of DEGs. (A) PPI analysis of DEGs (within the red box) and its relationship with seven other selected genes, (B) the directed network graph of DEGs, (C) relationship diagram of proteins with connection number ≥3 in PPI.
Toxins 15 00617 g006
Toxins 15 00617 g007
Figure 7. Conceptual diagram of the mechanism of ZEN on in vitro maturation of sheep oocytes.
Figure 7. Conceptual diagram of the mechanism of ZEN on in vitro maturation of sheep oocytes.
Toxins 15 00617 g007
Table 1. Primers for qRT-PCR.
Table 1. Primers for qRT-PCR.
PrimersPrimer SequencesProduct Length (bp)
SOD1GGCAATGTGAAGGCTGACAA130
TGCCCAAGTCATCTGGTCTT
SOD2GGACAAATCTGAGCCCCAAC180
CAATCTGTAAGCGTCCCTGC
CATCCAGCCCTGACAAAATGCTT242
AAAGCGGGTCCTATGTTCCA
GPXCAGTTTGGGCATCAGGAAAAC100
CGAAGAGCATGAAATTGGGC
CDK1ATGGCTTGGATCTGCTCTCGAA154
TGCTCTTGACACAACACAGGA
CCNB1GCTTGGAGACATCGGTAACA129
GGAGCCTTTTCCAGAGGTTTTG
BMP15GGACACCCTAGGGAAAACCG101
TGTATGTGCCAGGAGCCTCT
CDC20GGCTGAGCTGAAAGGTCACA214
AACACCGTGAGGAGTTGGTC
GDF9TGACAGAGCTTTGCGCTACA166
TGATGGAAAGGTTCCTGCCG
Table 2. The number of oocytes in Smart-seq2 sequencing samples.
Table 2. The number of oocytes in Smart-seq2 sequencing samples.
GroupCK1CK2CK3T1T2T3
Number of oocytes (pcs)251238222223235217
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Z.; Liu, Y.; Ma, T.; Lv, C.; Li, Y.; Duan, H.; Zhao, X.; Wang, J.; Zhang, Y. Smart-seq2 Technology Reveals a Novel Mechanism That Zearalenone Inhibits the In Vitro Maturation of Ovine Oocytes by Influencing TNFAIP6 Expression. Toxins 2023, 15, 617. https://doi.org/10.3390/toxins15100617

AMA Style

Li Z, Liu Y, Ma T, Lv C, Li Y, Duan H, Zhao X, Wang J, Zhang Y. Smart-seq2 Technology Reveals a Novel Mechanism That Zearalenone Inhibits the In Vitro Maturation of Ovine Oocytes by Influencing TNFAIP6 Expression. Toxins. 2023; 15(10):617. https://doi.org/10.3390/toxins15100617

Chicago/Turabian Style

Li, Zongshuai, Yali Liu, Tian Ma, Chen Lv, Yina Li, Hongwei Duan, Xingxu Zhao, Jianlin Wang, and Yong Zhang. 2023. "Smart-seq2 Technology Reveals a Novel Mechanism That Zearalenone Inhibits the In Vitro Maturation of Ovine Oocytes by Influencing TNFAIP6 Expression" Toxins 15, no. 10: 617. https://doi.org/10.3390/toxins15100617

APA Style

Li, Z., Liu, Y., Ma, T., Lv, C., Li, Y., Duan, H., Zhao, X., Wang, J., & Zhang, Y. (2023). Smart-seq2 Technology Reveals a Novel Mechanism That Zearalenone Inhibits the In Vitro Maturation of Ovine Oocytes by Influencing TNFAIP6 Expression. Toxins, 15(10), 617. https://doi.org/10.3390/toxins15100617

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop