Characterization and implementation of the MarathonRT template-switching reaction to expand the capabilities of RNA-Seq.

End-to-end RNA sequencing methods that capture 5'-sequence content without cumbersome library manipulations are of great interest, particularly for analysis of long RNAs. While template-switching methods have been developed for RNA sequencing by distributive short-read RTs, such as the MMLV RT enzymes used in SMART-Seq methods, they have not been adapted to leverage the power of ultraprocessive RTs, such as those that derive from group II self-splicing introns. To facilitate this transition, we dissected the individual processes that guide the enzymatic specificity and efficiency of the multi-step template switching reaction carried out by RT enzymes, in this case, by a well-characterized enzyme known as MarathonRT. Remarkably, this is the first study of its kind, for any RT. First, we characterized and optimized the enzymatic nontemplated addition (NTA) reaction that occurs when the RT enzyme extends past the RNA 5'-terminus, and we determined the nucleotide specificity of the NTA reaction. We then evaluated the binding specificity of specialized template-switching oligonucleotides, optimizing their sequences and chemical properties to guide efficient template switching reaction. Having dissected and optimized these individual steps, we then unified them into a procedure for performing RNA sequencing with MarathonRT enzymes, using a well-characterized RNA reference set. The resulting reads span a six-log range in transcript concentration and accurately represent the input RNA identities in both length and composition. We also performed RNA-seq starting from total human RNA and poly(A)-enriched RNA, with short and long-read sequencing demonstrating that MarathonRT enhances the discovery of unseen RNA molecules by conventional RT. Altogether, by employing mechanistic enzymology on RT enzymes and using them to modify RNA-seq technologies, we have generated a new pipeline for rapid, accurate sequencing of complex RNA libraries containing mixtures of long RNA transcripts.

Transcriptomic profiling and regulatory pathways of cardiac resident macrophages in aging.

Cardiovascular diseases are an array of age-related disorders, and accumulating evidence suggests a link between cardiac resident macrophages (CRMs) and the age-related disorders. However, how does CRMs alter with aging remains elusive. In the present study, aged mice (20 months old) have been employed to check for their cardiac structural and functional alterations, and the changes in the proportion of CRM subsets as well, followed by sorting of CRMs, including C-C Motif Chemokine Receptor 2 (CCR2)+ and CCR2- CRMs, which were subjected to Smart-Seq. Integrated analysis of the Smart-Seq data with three publicly available single-cell RNA-seq datasets revealed that inflammatory genes were drastic upregulated for both CCR2+ and CCR2- CRMs with aging, but genes germane to wound healing were downregulated for CCR2- CRMs, suggesting the differential functions of these two subsets. More importantly, inflammatory genes involved in damage sensing, complement cascades, and phagocytosis were largely upregulated in CCR2- CRMs, implying the imbalance of inflammatory response upon aging. Our work provides a comprehensive framework and transcriptional resource for assessing the impact of aging on CRMs with a potential for further understanding cardiac aging.

Smart seq2 revealed distinct molecular responses during in vitro porcine oocyte maturation before or after the addition of mogroside V.

Oocyte maturation involves both nuclear and cytoplasmic maturation. Mogroside V (MV) has been shown to enhance nuclear maturation, mitochondrial content, and developmental potential of porcine oocyte during in vitro maturation (IVM). However, the impact of MV on cytoplasmic maturation and its underlying mechanisms are not understood. This study aimed to assess the effect of MV on cytoplasmic maturation. Germinal vesicle (GV) oocytes treated with MV exhibited a noticeable increase in cortical granules (CGs) formation. Additionally, MV enhanced the expression of NNAT and improved glucose uptake in mature oocytes. Further insights were gained through Smart-seq2 analysis of RNA isolated from 100 oocytes. A total of 11,274 and 11,185 transcripts were identified in oocytes treated with and without MV, respectively. Among quantified genes, 438 differentially expressed genes (DEGs) were identified for further analysis. Gene Ontology (GO) enrichment analysis indicated that these DEGs were primarily involved in DNA repair regulation, cellular response to DNA damage, intracellular components, and organelles. Furthermore, the DEGs were significantly enriched in three KEGG pathways: fatty acid synthesis, pyruvate metabolism, and WNT signalling. To validate the results, lipid droplets (LD) and triglyceride (TG) were examined. MV led to an increase in the accumulation of LD and TG production in mature oocytes. These findings suggest that MV enhances cytoplasmic maturation by promoting lipid droplet synthesis. Overall, this study provides valuable insights into the mechanisms through which MV improves oocyte quality during IVM. The results have significant implications for research in livestock reproduction and offer guidance for future studies in this field.

A method to investigate muscle target-specific transcriptional signatures of single motor neurons.

BACKGROUND: Motor neurons in the vertebrate spinal cord have long served as a paradigm to study the transcriptional logic of cell type specification and differentiation. At limb levels, pool-specific transcriptional signatures first restrict innervation to only one particular muscle in the periphery, and get refined, once muscle connection has been established. Accordingly, to study the transcriptional dynamics and specificity of the system, a method for establishing muscle target-specific motor neuron transcriptomes would be required. RESULTS: To investigate target-specific transcriptional signatures of single motor neurons, here we combine ex-ovo retrograde axonal labeling in mid-gestation chicken embryos with manual isolation of individual fluorescent cells and Smart-seq2 single-cell RNA-sequencing. We validate our method by injecting the dorsal extensor metacarpi radialis and ventral flexor digiti quarti wing muscles and harvesting a total of 50 fluorescently labeled cells, in which we detect up to 12,000 transcribed genes. Additionally, we present visual cues and cDNA metrics predictive of sequencing success. CONCLUSIONS: Our method provides a unique approach to study muscle target-specific motor neuron transcriptomes at a single-cell resolution. We anticipate that our method will provide key insights into the transcriptional logic underlying motor neuron pool specialization and proper neuromuscular circuit assembly and refinement.

SMART approaches for genome-wide analyses of skeletal muscle stem and niche cells.

Muscle stem cells (MuSCs) also called satellite cells are the building blocks of skeletal muscle, the largest tissue in the human body which is formed primarily of myofibers. While MuSCs are the principal cells that directly contribute to the formation of the muscle fibers, their ability to do so depends on critical interactions with a vast array of nonmyogenic cells within their niche environment. Therefore, understanding the nature of communication between MuSCs and their niche is of key importance to understand how the skeletal muscle is maintained and regenerated after injury. MuSCs are rare and therefore difficult to study in vivo within the context of their niche environment. The advent of single-cell technologies, such as switching mechanism at 5' end of the RNA template (SMART) and tagmentation based technologies using hyperactive transposase, afford the unprecedented opportunity to perform whole transcriptome and epigenome studies on rare cells within their niche environment. In this review, we will delve into how single-cell technologies can be applied to the study of MuSCs and muscle-resident niche cells and the impact this can have on our understanding of MuSC biology and skeletal muscle regeneration.

Neuronal activity-related transcription is blunted in immature compared to mature dentate granule cells.

Immature dentate granule cells (DGCs) generated in the hippocampus during adulthood are believed to play a unique role in dentate gyrus (DG) function. Although immature DGCs have hyperexcitable membrane properties in vitro, the consequences of this hyperexcitability in vivo remain unclear. In particular, the relationship between experiences that activate the DG, such as exploration of a novel environment (NE), and downstream molecular processes that modify DG circuitry in response to cellular activation is unknown in this cell population. We first performed quantification of immediate early gene (IEG) proteins in immature (5-week-old) and mature (13-week-old) DGCs from mice exposed to a NE. Paradoxically, we observed lower IEG protein expression in hyperexcitable immature DGCs. We then isolated nuclei from active and inactive immature DGCs and performed single-nuclei RNA-Sequencing. Compared to mature nuclei collected from the same animal, immature DGC nuclei showed less activity-induced transcriptional change, even though they were classified as active based on expression of ARC protein. These results demonstrate that the coupling of spatial exploration, cellular activation, and transcriptional change differs between immature and mature DGCs, with blunted activity-induced changes in immature cells.

A survey of RNA editing at single-cell resolution links interneurons to schizophrenia and autism.

Conversion of adenosine to inosine in RNA by ADAR enzymes, termed "RNA editing," is essential for healthy brain development. Editing is dysregulated in neuropsychiatric diseases, but has not yet been investigated at scale at the level of individual neurons. We quantified RNA editing sites in nuclear transcriptomes of 3055 neurons from six cortical regions of a neurotypical female donor, and found 41,930 sites present in at least ten nuclei. Most sites were located within Alu repeats in introns or 3' UTRs, and approximately 80% were cataloged in public RNA editing databases. We identified 9285 putative novel editing sites, 29% of which were also detectable in unrelated donors. Intersection with results from bulk RNA-seq studies provided cell-type and spatial context for 1730 sites that are differentially edited in schizophrenic brain donors, and 910 such sites in autistic donors. Autism-related genes were also enriched with editing sites predicted to modify RNA structure. Inhibitory neurons showed higher overall transcriptome editing than excitatory neurons, and the highest editing rates were observed in the frontal cortex. We used generalized linear models to identify differentially edited sites and genes between cell types. Twenty nine genes were preferentially edited in excitatory neurons, and 43 genes were edited more heavily in inhibitory neurons, including RBFOX1, its target genes, and genes in the autism-associated Prader-Willi locus (15q11). The abundance of SNORD115/116 genes from locus 15q11 was positively associated with editing activity across the transcriptome. We contend that insufficient editing of autism-related genes in inhibitory neurons may contribute to the specific perturbation of those cells in autism.

Transcriptome analyses reveal transcriptional profiles of horse oocytes before and after in vitro maturation.

Oocyte in vitro maturation is necessary for the study and application of animal-assisted reproduction technology in animal reproduction and breeding. The comprehensive transcriptional profile of equine oocyte maturated in vitro has not been fully mined yet, which makes many key transcriptional events still unidentified. Here, Smart-seq2 was performed to analyse the gene expression pattern and the underlying regulatory mechanism of horse germinal vesicle (GV) and in vitro metaphase II (MII) oocytes. The results showed that 6402 genes (2640 up-regulated and 3762 down-regulated in MII samples compared to GV) and 4021 lncRNA transcripts (1210 up-regulated and 2811 down-regulated in MII samples compared to GV) were differentially expressed in GV and MII oocytes. Further, GO and KEGG analysis found that differentially expressed mRNAs and lncRNAs were mainly enriched in the pathways related to energy and lipid metabolism. In addition, LGALS3 was found a key gene in mediating the regulation of oocyte meiosis recovery and fertilization ability. This study provides novel knowledge about gene expression and energy metabolism during equine oocyte maturation and a reference for the further study and application of assisted reproductive technology in horse reproduction and breeding.

Benchmarking full-length transcript single cell mRNA sequencing protocols.

BACKGROUND: Single cell mRNA sequencing technologies have transformed our understanding of cellular heterogeneity and identity. For sensitive discovery or clinical marker estimation where high transcript capture per cell is needed only plate-based techniques currently offer sufficient resolution. RESULTS: Here, we present a performance evaluation of four different plate-based scRNA-seq protocols. Our evaluation is aimed towards applications taxing high gene detection sensitivity, reproducibility between samples, and minimum hands-on time, as is required, for example, in clinical use. We included two commercial kits, NEBNext® Single Cell/ Low Input RNA Library Prep Kit (NEB®), SMART-seq® HT kit (Takara®), and the non-commercial protocols Genome & Transcriptome sequencing (G&T) and SMART-seq3 (SS3). G&T delivered the highest detection of genes per single cell. SS3 presented the highest gene detection per single cell at the lowest price. Takara® kit presented similar high gene detection per single cell, and high reproducibility between samples, but at the absolute highest price. NEB® delivered a lower detection of genes but remains an alternative to more expensive commercial kits. CONCLUSION: For the tested kits we found that ease-of-use came at higher prices. Takara can be selected for its ease-of-use to analyse a few samples, but we recommend the cheaper G&T-seq or SS3 for laboratories where a substantial sample flow can be expected.

Single-cell RNA sequencing of psoriatic skin identifies pathogenic Tc17 cell subsets and reveals distinctions between CD8+ T cells in autoimmunity and cancer.

BACKGROUND: Psoriasis is an inflammatory, IL-17-driven skin disease in which autoantigen-induced CD8+ T cells have been identified as pathogenic drivers. OBJECTIVE: Our study focused on comprehensively characterizing the phenotypic variation of CD8+ T cells in psoriatic lesions. METHODS: We used single-cell RNA sequencing to compare CD8+ T-cell transcriptomic heterogeneity between psoriatic and healthy skin. RESULTS: We identified 11 transcriptionally diverse CD8+ T-cell subsets in psoriatic and healthy skin. Among several inflammatory subsets enriched in psoriatic skin, we observed 2 Tc17 cell subsets that were metabolically divergent, were developmentally related, and expressed CXCL13, which we found to be a biomarker of psoriasis severity and which achieved comparable or greater accuracy than IL17A in a support vector machine classifier of psoriasis and healthy transcriptomes. Despite high coinhibitory receptor expression in the Tc17 cell clusters, a comparison of these cells with melanoma-infiltrating CD8+ T cells revealed upregulated cytokine, cytolytic, and metabolic transcriptional activity in the psoriatic cells that differed from an exhaustion program. CONCLUSION: Using high-resolution single-cell profiling in tissue, we have uncovered the diverse landscape of CD8+ T cells in psoriatic and healthy skin, including 2 nonexhausted Tc17 cell subsets associated with disease severity.

Maturation of 23S rRNA includes removal of helix H1 in many bacteria.

In most bacteria, the three ribosomal RNAs (rRNAs) are encoded together in each of several near-identical operons. As soon as the nascent precursor rRNA emerges from RNA polymerase, ribosome assembly begins. This process entails ribosomal protein binding, rRNA folding, rRNA modification, and rRNA processing. In the model organisms Escherichia coli and Bacillus subtilis, rRNA processing results in similar mature rRNAs, despite substantial differences in the cohort of RNAses involved. A recent study of Flavobacterium johnsoniae, a member of the phylum Bacteroidota (formerly Bacteroidetes), revealed that helix H1 of 23S rRNA is absent from ribosomes, apparently a consequence of rRNA maturation. In this work, we mined RNA-seq data from 19 individual organisms and ocean metatranscriptomic samples to compare rRNA processing across diverse bacterial lineages. We found that mature ribosomes from multiple clades lack H1, and typically these ribosomes also lack an encoded H98. For all groups analysed, H1 is predicted to form in precursor rRNA as part of a longer leader-trailer helix. Hence, we infer that evolutionary loss of H98 sets the stage for H1 removal during 50S subunit maturation.

RNA-seq between asexual archeospores and meiosis-related conchospores in Neopyropia yezoensis using Smart-seq2.

In the life cycle of Neopyropia yezoensis, a potential model system for marine macroalgae, both asexual archeospores and meiosis-related conchospores develop into thalli (gametophyte). To understand this special life phenomenon in macroalgae, we picked out the two kinds of spores (10-30 cells in each sample) and conducted RNA-seq using Smart-seq2. Comparative analysis showed that light capture and carbon fixation associated differentially expressed genes (DEGs) were upregulated in archeospores, thus indicating that archeospores are in a state of rapid vegetative growth. In conchospores, protein synthesis and degradation, especially molecular chaperone, associated DEGs were up-regulated, indicating that complex life activities might be occurring in conchospores. There were 68 genes related to DNA replication and repair expressed in conchospores, showing that active DNA replication might occur in conchospores. Moreover, we found that one conchospore specifically expressed DEG (py04595: DNA helicase) only in diploid stages (conchocelis, sporangial filament) and three archeospores specifically expressed DEGs only in haploid stages (thalli). These molecular level results indicated that conchospores were closer to diploid, and might be the meiotic mother cells of N. yezoensis. In addition, we found that the knotted-like homeobox gene (PyKNOX), which might relate to the transition of gametophyte from sporophyte, was only expressed in sporophyte generation but not expressed in conchospores, archeospores and thalli, indicating the morphogenesis of gametophyte sin N. yezoensis might require the inactivation of PyKNOX.

High-resolution genome-wide expression analysis of single myofibers using SMART-Seq.

Skeletal muscle is a heterogeneous tissue. Individual myofibers that make up muscle tissue exhibit variation in their metabolic and contractile properties. Although biochemical and histological assays are available to study myofiber heterogeneity, efficient methods to analyze the whole transcriptome of individual myofibers are lacking. Here, we report on a single-myofiber RNA-sequencing (smfRNA-Seq) approach to analyze the whole transcriptome of individual myofibers by combining single-fiber isolation with Switching Mechanism at 5' end of RNA Template (SMART) technology. Using smfRNA-Seq, we first determined the genes that are expressed in the whole muscle, including in nonmyogenic cells. We also analyzed the differences in the transcriptome of myofibers from young and old mice to validate the effectiveness of this new method. Our results suggest that aging leads to significant changes in the expression of metabolic genes, such as Nos1, and structural genes, such as Myl1, in myofibers. We conclude that smfRNA-Seq is a powerful tool to study developmental, disease-related, and age-related changes in the gene expression profile of skeletal muscle.

An efficient single-cell transcriptomics workflow for microbial eukaryotes benchmarked on Giardia intestinalis cells.

BACKGROUND: Most diversity in the eukaryotic tree of life is represented by microbial eukaryotes, which is a polyphyletic group also referred to as protists. Among the protists, currently sequenced genomes and transcriptomes give a biased view of the actual diversity. This biased view is partly caused by the scientific community, which has prioritized certain microbes of biomedical and agricultural importance. Additionally, some protists remain difficult to maintain in cultures, which further influences what has been studied. It is now possible to bypass the time-consuming process of cultivation and directly analyze the gene content of single protist cells. Single-cell genomics was used in the first experiments where individual protists cells were genomically explored. Unfortunately, single-cell genomics for protists is often associated with low genome recovery and the assembly process can be complicated because of repetitive intergenic regions. Sequencing repetitive sequences can be avoided if single-cell transcriptomics is used, which only targets the part of the genome that is transcribed. RESULTS: In this study we test different modifications of Smart-seq2, a single-cell RNA sequencing protocol originally developed for mammalian cells, to establish a robust and more cost-efficient workflow for protists. The diplomonad Giardia intestinalis was used in all experiments and the available genome for this species allowed us to benchmark our results. We could observe increased transcript recovery when freeze-thaw cycles were added as an extra step to the Smart-seq2 protocol. Further we reduced the reaction volume and purified the amplified cDNA with alternative beads to test different cost-reducing changes of Smart-seq2. Neither improved the procedure, and reducing the volumes by half led to significantly fewer genes detected. We also added a 5' biotin modification to our primers and reduced the concentration of oligo-dT, to potentially reduce generation of artifacts. Except adding freeze-thaw cycles and reducing the volume, no other modifications lead to a significant change in gene detection. Therefore, we suggest adding freeze-thaw cycles to Smart-seq2 when working with protists and further consider our other modification described to improve cost and time-efficiency. CONCLUSIONS: The presented single-cell RNA sequencing workflow represents an efficient method to explore the diversity and cell biology of individual protist cells.

Comparative transcriptomics implicate mitochondrial and neurodevelopmental impairments in larval zebrafish (Danio rerio) exposed to two selective serotonin reuptake inhibitors (SSRIs).

Pharmaceuticals and personal care products are emerging contaminants that are increasingly detected in the environment worldwide. Certain classes of pharmaceuticals, such as selective serotonin reuptake inhibitors (SSRIs), are a major environmental concern due to their widespread use and the fact that these compounds are designed to have biological effects at low doses. A complication in predicting toxic effects of SSRIs in nontarget organisms is that their mechanism of action is not fully understood. To better understand the potential toxic effects of SSRIs, we employed an ultra-low input RNA-sequencing method to identify potential pathways that are affected by early exposure to two SSRIs (fluoxetine and paroxetine). We exposed wildtype zebrafish (Danio rerio) embryos to 100 µg/L of either fluoxetine or paroxetine for 6 days before extracting and sequencing mRNA from individual larval brains. Differential gene expression analysis identified 1550 genes that were significantly affected by SSRI exposure with a core set of 138 genes altered by both SSRIs. Weighted gene co-expression network analysis identified 7 modules of genes whose expression patterns were significantly correlated with SSRI exposure. Functional enrichment analysis of differentially expressed genes as well as network module genes repeatedly identified various terms associated with mitochondrial and neuronal structures, mitochondrial respiration, and neurodevelopmental processes. The enrichment of these terms indicates that toxic effects of SSRI exposure are likely caused by mitochondrial dysfunction and subsequent neurodevelopmental effects. To our knowledge, this is the first effort to study the tissue-specific transcriptomic effects of SSRIs in developing zebrafish, providing specific, high resolution molecular data regarding the sublethal effects of SSRI exposure.

Single-cell analysis to understand the diversity of immune cell types that drive disease pathogenesis.

Single-cell next-generation sequencing assays are powerful tools to understand the nature of the immune cells that drive disease pathogenesis. In this brief review we explain the value of performing assays at single-cell resolution to better understand the pathogenesis of allergy, asthma, and other lung diseases. We explain the challenges in performing single-cell studies of airways and lung samples from patients with lung diseases. A major limitation comes from the amount of diseased tissue that can be used for research purposes. Finally, we discuss which sequencing strategies can be used to successfully investigate airway and lung diseases at single-cell resolution.

Single-Cell DNA-Seq and RNA-Seq in Cancer Using the C1 System.

Heterogeneous phenotypes of cancer cells enable them to adapt to various environments. The heterogeneity results from diversity of genome, transcriptome, and epigenome at a single-cell level. The C1 system can automatically perform single-cell capture and whole genome amplification (WGA) or whole transcription amplification (WTA) by MDA or Smart-Seq, respectively. Here, we describe the protocols for WGA and WTA from a single cell by using the C1 system and the protocols for sequence library preparation from amplified gDNA and cDNA. We also described about the computational analysis for single-cell data of cancer.

Single-cell RNA-sequencing: The future of genome biology is now.

Genome-wide single-cell analysis represents the ultimate frontier of genomics research. In particular, single-cell RNA-sequencing (scRNA-seq) studies have been boosted in the last few years by an explosion of new technologies enabling the study of the transcriptomic landscape of thousands of single cells in complex multicellular organisms. More sensitive and automated methods are being continuously developed and promise to deliver better data quality and higher throughput with less hands-on time. The outstanding amount of knowledge that is going to be gained from present and future studies will have a profound impact in many aspects of our society, from the introduction of truly tailored cancer treatments, to a better understanding of antibiotic resistance and host-pathogen interactions; from the discovery of the mechanisms regulating stem cell differentiation to the characterization of the early event of human embryogenesis.

The Mechanism of Ovule Abortion in Self-Pollinated 'Hanfu' Apple Fruits and Related Gene Screening.

Apples exhibit S-RNase-mediated self-incompatibility and typically require cross-pollination in nature. 'Hanfu' is a cultivar that produces abundant fruit after self-pollination, although it also shows a high rate of seed abortion afterwards, which greatly reduces fruit quality. In this study, we investigated the ovule development process and the mechanism of ovule abortion in apples after self-pollination. Using a DIC microscope and biomicroscope, we found that the abortion of apple ovules occurs before embryo formation and results from the failure of sperm-egg fusion. Further, we used laser-assisted microdissection (LAM) cutting and sperm and egg cell sequencing at different periods after pollination to obtain the genes related to ovule abortion. The top 40 differentially expressed genes (DEGs) were further verified, and the results were consistent with switching the mechanism at the 5' end of the RNA transcript (SMART-seq). Through this study, we can preliminarily clarify the mechanism of ovule abortion in self-pollinated apple fruits and provide a gene reserve for further study and improvement of 'Hanfu' apple fruit quality.

Screening and Mechanism Study of Three Antagonistic Drugs, Oxysophoridine, Rutin, and Phellodendrine, against Zearalenone-Induced Reproductive Toxicity in Ovine Oocytes.

Zearalenone (ZEN) is a common fungal toxin with reproductive toxicity in various grains. It poses a serious threat to ovine and other animal husbandry industries, as well as human reproductive health. Therefore, investigating the mechanism of toxicity and screening antagonistic drugs are of great importance. In this study, based on the natural compound library and previous Smart-seq2 results, antioxidant and anti-apoptotic drugs were selected for screening as potential antagonistic drugs. Three natural plant compounds (oxysophoridine, rutin, and phellodendrine) were screened for their ability to counteract the reproductive toxicity of ZEN on ovine oocytes in vitro using quantitative polymerase chain reaction (qPCR) and reactive oxygen species detection. The compounds exhibited varying pharmacological effects, notably impacting the expression of antioxidant (GPX, SOD1, and SOD2), autophagic (ATG3, ULK2, and LC3), and apoptotic (CAS3, CAS8, and CAS9) genes. Oxysophoridine promoted GPX, SOD1, ULK2, and LC3 expression, while inhibiting CAS3 and CAS8 expression. Rutin promoted SOD2 and ATG3 expression, and inhibited CAS3 and CAS9 expression. Phellodendrine promoted SOD2 and ATG3 expression, and inhibited CAS9 expression. However, all compounds promoted the expression of genes related to cell cycle, spindle checkpoint, oocyte maturation, and cumulus expansion factors. Although the three drugs had different regulatory mechanisms in enhancing antioxidant capacity, enhancing autophagy, and inhibiting cell apoptosis, they all maintained a stable intracellular environment and a normal cell cycle, promoted oocyte maturation and release of cumulus expansion factors, and, ultimately, counteracted ZEN reproductive toxicity to promote the in vitro maturation of ovine oocytes. This study identified three drugs that antagonize the reproductive toxicity of ZEN on ovine oocytes, and compared their mechanisms of action, providing data support and a theoretical basis for their subsequent application in the ovine breeding industry, reducing losses in the breeding industry, screening of ZEN reproductive toxicity antagonists and various toxin antagonists, improving the study of ZEN reproductive toxicity mechanisms, and even protection of human reproductive health.