Temporal transcriptomic changes in microRNAs involved in the host immune response and metabolism during Neospora caninum infection
Parasites & Vectors volume 16, Article number: 28 (2023)
Neospora caninum infection is a major cause of abortion in cattle, which results in serious economic losses to the cattle industry. However, there are no effective drugs or vaccines for the control of N. caninum infections. There is increasing evidence that microRNAs (miRNAs) are involved in many physiological and pathological processes, and dysregulated expression of host miRNAs and the biological implications of this have been reported for infections by various protozoan parasites. However, to our knowledge, there is presently no published information on host miRNA expression during N. caninum infection.
The expression profiles of miRNAs were investigated by RNA sequencing (RNA-seq) in caprine endometrial epithelial cells (EECs) infected with N. caninum at 24 h post infection (pi) and 48 hpi, and the functions of differentially expressed (DE) miRNAs were predicted by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The transcriptome data were validated by using quantitative real-time polymerase chain reaction. One of the upregulated DEmiRNAs, namely chi-miR-146a, was selected to study the effect of DEmiRNAs on the propagation of N. caninum tachyzoites in caprine EECs.
RNA-seq showed 18 (17 up- and one downregulated) and 79 (54 up- and 25 downregulated) DEmiRNAs at 24 hpi and 48 hpi, respectively. Quantitative real-time polymerase chain reaction analysis of 13 randomly selected DEmiRNAs (10 up- and three downregulated miRNAs) confirmed the validity of the RNA-seq data. A total of 7835 messenger RNAs were predicted to be potential targets for 66 DEmiRNAs, and GO and KEGG enrichment analysis of these predicted targets revealed that DEmiRNAs altered by N. caninum infection may be involved in host immune responses (e.g. Fc gamma R-mediated phagocytosis, Toll-like receptor signaling pathway, tumor necrosis factor signaling pathway, transforming growth factor-β signaling pathway, mitogen-activated protein kinase signaling pathway) and metabolic pathways (e.g. lysine degradation, insulin signaling pathway, AMP-activated protein kinase signaling pathway, Rap1 signaling pathway, calcium signaling pathway). Upregulated chi-miR-146a was found to promote N. caninum propagation in caprine EECs.
This is, to our knowledge, the first report on the expression profiles of host miRNAs during infection with N. caninum, and shows that chi-miR-146a may promote N. caninum propagation in host cells. The novel findings of the present study should help to elucidate the interactions between host cells and N. caninum.
Neospora caninum is an important intracellular protozoan parasite that infects a variety of animals, including domesticated ruminants (e.g. cattle, sheep, goats) and wildlife [1, 2]. Infection with N. caninum has been reported to cause reproductive disorders in pregnant animals and neuromuscular disorders in newborn animals [3, 4]. Notably, neosporosis, which is caused by N. caninum infection, was found to be one of the main causes of abortion in cattle [5, 6]. A meta-analysis also found that the combined seroprevalence of N. caninum infection in goats was 5.99%, and seropositive goats were more likely to abort than seronegative ones . Although there is insufficient evidence to indicate that N. caninum is of zoonotic significance, antibodies against N. caninum have been detected in humans . However, no effective drugs or vaccines have yet been developed to control neosporosis or infection with N. caninum [9,10,11,12]. At present, eliminating cattle positive for antibodies against N. caninum is the most common strategy used in cattle production against this parasite [13, 14]. This results in huge economic losses, with a total annual loss of more than $1 billion to the cattle industry of 10 countries alone .
MicroRNAs (miRNAs) are a type of endogenous non-coding RNA (ncRNA), approximately 19–25 nucleotides (nt) in length . In studies by Lei  and Bartel , miRNA expression was generally dysregulated, and miRNAs were found to function as key elements in the regulation of target mRNAs through their combination with the 3' untranslated region of the latter. They have also been found to play a role in translational inhibition and/or mRNA degradation at post-transcriptional levels  in a great number of biological processes (e.g. cell proliferation , apoptosis [21, 22], autophagy  and pyroptosis ), and during the progression of various diseases [25, 26], including infections with protozoan parasites. For example, a total of 81, 126, 82 and 131 miRNAs were differentially expressed (DE) in porcine alveolar macrophages , pig splenocytes , livers of domestic cats , and mouse spleen , respectively, during infections with Toxoplasma gondii, an obligate intracellular zoonotic protozoan parasite similar to N. caninum in its morphological and biological features but differing in its life cycle and biological effects on host cells [31,32,33]. Previously, our group found that N. caninum infection significantly altered the expression patterns of host cell long ncRNAs, which are some of the most important regulatory ncRNAs . In the present study, we investigated dynamic miRNA profiles in caprine endometrial epithelial cells (EECs) during their infections with N. caninum.
African green monkey kidney cells (Vero cells) and caprine EECs were provided by Prof. Xuefeng Qi and Prof. Yaping Jin from Northwest A&F University (Shaanxi, China), respectively. Vero cells were cultured in DMEM medium supplemented with 10% fetal bovine serum (FBS; ExCell Bio, Shanghai, China), and caprine EECs were cultured in DMEM/F12 medium supplemented with 10% FBS. All cells were cultured with 5% CO2 at 37 ℃ in a CO2 incubator (Thermo Fisher Scientific, Wilmington, NC).
Parasite purification and infection
NC-1 N. caninum tachyzoites were gifted by Prof. Qun Liu from the China Agricultural University (Beijing, China), and passaged in Vero cells in DMEM medium supplemented with 2% FBS, 100 U/ml penicillin, and 100 μg/ml streptomycin. Before infection, Vero cells infected with tachyzoites were scraped, pipetted using a 27-gauge needle, filtered through a 5.0-µm filter, and centrifuged at 716 g for 10 min. After the supernatant had been discarded, the pellets were resuspended in fresh medium for further study, and the number of tachyzoites was counted with a hemocytometer.
The in vitro infection model in caprine EECs was established in accordance with our previous study  at a multiplicity of infection of 3:1 (tachyzoites:cells).
Sample preparation and RNA extraction
Caprine EECs infected (experimental group) or not infected (control group) with N. caninum tachyzoites were collected at 24 h (experimental groups, TZ1_24h to TZ3_24h; control groups, C1_24h to C3_24h) and 48 h (experimental groups, TZ1_48h to TZ3_48h; control groups, C1_48h to C3_48h) post-infection [pi; hours pi (hpi)], treated with Trizol reagent (Invitrogen, Carlsbad, CA), and stored at —80 °C until RNA extraction. All the analyses were performed with three biological replicates.
The total RNA for each sample was extracted by using a mirVana miRNA isolation kit (Ambion, Austin, TX) following the manufacturer’s protocol. The extracted RNA samples were quantified by using Nanodrop 2000 (Thermo Fisher Scientific). The RNA integrity of each RNA sample was checked using an Agilent 2100 Bioanalyzer (Agilent Technologies, CA), and RNA samples with 28S/18S ≥ 0.7 and an RNA integrity number ≥ 7 were used for further analysis.
Small RNA library construction and RNA sequencing
A total of 1 μg RNA from each sample was used to construct the small RNA library by using TruSeq Small RNA Sample Prep Kits (Illumina, USA) following the manufacturer’s recommendations. Briefly, the total RNA samples were ligated to adapters at each end, and then the adapter-ligated RNA specimens were reverse transcribed into complementary DNA (cDNA) samples for polymerase chain reaction (PCR) amplification. PCR products ranging from 140 to 160 base pairs in length were purified as a small RNA library. The quality of the library was assessed on the Agilent Bioanalyzer 2100 system by using DNA high-sensitivity chips. These small RNA libraries were sequenced using the Illumina HiSeq X Ten platform to generate a 150-base pairs paired-end. All analyses were performed by OE Biotechnology (Shanghai, China).
The basic reads obtained by RNA sequencing (RNA-seq) were converted into raw reads by base calling. Adapter sequences of raw reads were removed by using cutadapt (version 1.14) , and sequences shorter than 15 nt and longer than 41 nt were filtered out. The obtained sequences were subjected to Q20 quality control by using fastx_toolkit (version 0.0.13) software (http://hannonlab.cshl.edu/fastx_toolkit), and the reads containing N bases were filtered out by using NGSQCToolkit (version 2.3.2)  to obtain high-quality clean reads for subsequent analysis.
The obtained clean reads were mapped to the reference genome from Capra hircus (ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/001/704/415/GCF_001704415.1_ ARS1/gcf_00-1,704,415.1_ars1_genome.fna.gz) and the percentage of these reads that were aligned to the genome counted. Using bowtie software , the clean reads were aligned with the Rfam v.10.0 database (http://www.sanger.ac.uk/software/Rfam)  to annotate ribosomal RNA, small nuclear RNA, small nucleolar RNA, and transfer RNA, and then the sequences that were annotated to the Rfam database were filtered and removed. The sequences aligned to the transcripts less than 15 nt and more than 26 nt in length were removed. Using Repeat Masker software , the filtered sequences were aligned to the repeat database to identify possible repetitive sequences and filtered. The identified miRNAs were aligned with the miRbase v.22.0 database (http://www.mirbase.org/)  to analyze the expression levels of the known miRNAs, and the unannotated small RNA sequences were analyzed by Mirdeep2  to predict novel miRNAs, and the secondary structure of these novel miRNAs were predicted by using RNAfold software.
To analyze the expression levels of DEmiRNAs between the experimental group and the control group at 24 hpi and 48 hpi after N. caninum infection, the expression level of each sample was calculated using transcripts per million , and the differential expression was calculated by using DESeq2 1.16.1, with a q-value < 0.05 and | log2 fold change |> 0 considered as significantly different.
Verification of miRNA expression by quantitative real-time PCR
Thirteen (10 up- and three downregulated) DEmiRNAs were randomly selected for quantitative real-time PCR (qRT-PCR) analysis to verify the sequencing data. A total of 12 samples from the experimental (six samples) and control (six samples) groups were collected at 24 hpi and 48 hpi, and the total RNA of each sample was extracted using Trizol reagent (Invitrogen, Carlsbad, CA). The quality of each RNA sample was evaluated by using a Nano-100 spectrophotometer (Hangzhou, China). The cDNA samples were synthesized using 0.8-μg RNA samples in a Mir-XTM miRNA First-Strand Synthesis Kit, in accordance with the manufacturer’s instructions. qRT-PCR reactions were performed in a 25-μl mixture containing 2 μl cDNA, 9 μl ddH2O, 0.5 μl ROX reference dye, 0.5 μl miRNA-specific forward primer, 0.5 μl mRQ 3' primer, and 12.5 μl TB Green Premix Ex TaqTM II (Tli RNaseH Plus) under the following conditions: 95 ℃ for 10 min, 40 cycles of 95 ℃ for 15 s, 55–58 ℃ (Additional file 1: Data S1) for 30 s, in a Four-channel Real-time Fluorescence Quantitative PCR system (Tianlong TL988, Shaanxi, China). Three biological replicates were performed for each reaction, and the u6 small nuclear RNA gene was used as the expression level standard. The relative expression of each gene was calculated by using the 2−ΔΔCt method, with P < 0.05 indicating statistically significant difference.
Target prediction and function analysis of DEmiRNAs
Targets of all DEmiRNAs were predicted by using miranda software  with S ≥ 150, ΔG ≤ − 30 kcal/mol and strict 5’ seed pairing. Functions of the DEmiRNAs were predicted by submitting their targets to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases for enrichment analysis, using R based on the hypergeometric distribution. GO enrichment analysis included three elements: biological process (BP), cellular composition (CC) and molecular function (MF). KEGG enrichment analysis was used to identify the pathways and predict biological functions. Elements with a q-value < 0.05 were considered to be significantly enriched.
Propagation of N. caninum tachyzoites affected by DEmiRNAs
To study the effect of DEmiRNAs on the propagation of N. caninum tachyzoites in caprine EECs, one upregulated miRNA, namely chi-miR-146a, was selected. The mimics and inhibitor of chi-miR-146a together with their negative controls were obtained from GenePharma (Shanghai, China). First, a total of 100 pmol of mimics or inhibitor was transfected into caprine EECs by using Lipofectamine 2000 reagent (Invitrogen, Gaithersburg, MD). Then, N. caninum tachyzoites were infected at 24 h post-transfection with a multiplicity of infection of 3:1 (tachyzoites:cells), and the number of parasites per vacuole was determined at 30 hpi or 42 hpi in 100 parasitophorous vacuoles under fluorescence microscopy (Olympus, Tokyo, Japan).
The differences between the control and the experimental groups were analyzed by using GraphPad Prism 5.0 software (http://www.graphpad.com), and P < 0.05 was considered to indicate a statistically significant difference in the two-tailed t-test, which was used as the parametric test.
Identification of caprine miRNAs
A total of 343.84 M raw reads were generated from 12 samples by RNA-seq. Through QC quality control (adapter sequences, N base sequence and Q20), a total of 294.44 M clean reads and 15.77 M unique reads were obtained (Table 1), and 75.10% (221112855) of these clean reads (294429204) were matched to the reference genome of C. hircus (Additional file 2: Data S2). Through Rfam alignment, transcript sequence alignment, repeat sequence alignment and miRNA (miRBase database) alignment annotation, 348 known and 891 novel predicted miRNAs were identified (Additional file 3: Data S3).
Differential expression profiles of miRNAs in caprine EECs infected with N. caninum
To analyze the differential expression profiles of miRNAs in caprine EECs during N. caninum infection, two categories were used for comparison, namely TZ_24h-vs-C_24h and TZ_48h-vs-C_48h. When using the criteria of a q-value < 0.05 and | log2fold change |> 0, a total of 84 (59 up- and 25 downregulated) DEmiRNAs were found to be dysregulated, with 13 (12 up- and one downregulated) miRNAs differentially expressed in both of the categories (TZ_24h-vs-C_24h and TZ_48h-vs-C_48h). Notably, five of the upregulated miRNAs were only found in category TZ_24h-vs-C_24h, while 66 (42 up- and 24 downregulated) DEmiRNAs were detected in category TZ_48h-vs-C_48h (Fig. 1). The volcano maps (Fig. 2a, b) show the overall distribution of the DEmiRNAs, and the hierarchical clustering heatmaps (Fig. 2c, d) clearly differentiate the control and experimental groups. Detailed information on the DEmiRNAs is given in Additional file 4: Data S4.
Validation of DEmiRNAs by qRT-PCR
To validate the accuracy of the RNA-seq data, five upregulated (chi-miR-146b-5p, chi-miR-146a, chi-miR-200a, chi-miR-218, chi-miR-381) DEmiRNAs were selected for qRT-PCR analysis from category TZ_24h-vs-C_24h (Fig. 3a). Five up- (chi-miR-383, chi-miR-155-3p, chi-miR-125b-3p, chi-miR-30a-3p, chi-miR-147-5p) and three downregulated (chi-miR-877-3p, chi-miR-214-3p, chi-miR-20a-5p) DEmiRNAs were selected for qRT-PCR analysis from category TZ_48h-vs-C_48h (Fig. 3b). The expression of all the determined DEmiRNAs was consistent with the RNA-seq data, indicating high credibility of the latter.
Target prediction and functional annotation of DEmiRNAs
A total of 7835 mRNAs were predicted to be potential targets for 66 DEmiRNAs by miranda software, and comprised 2094 and 7712 mRNAs for 12 and 63 DEmiRNAs at 24 hpi and 48 hpi, respectively (Additional file 5: Data S5). GO analysis of these targets indicated that 883 and 1752 terms were significantly enriched for the categories TZ_24h-vs-C_24h and TZ_48h-vs-C_48h, respectively (Additional file 6: Data S6). The top 30 significantly enriched terms are shown in Fig. 4. Interestingly, most of the significantly enriched terms in BP, CC and MF were the same for the two categories, while only one, one and two terms in BP, CC and MF, respectively, differed between the two categories.
KEGG enrichment analysis of targets for DEmiRNAs showed that 42 and 90 terms were significantly enriched for the categories TZ_24h-vs-C_24h and TZ_48h-vs-C_48h, respectively (Additional file 7: Data S7). The top 20 significantly enriched terms are shown in Fig. 5. Surprisingly, in the KEGG enrichment analysis, contrary to the results of the GO analysis, most of the terms that were significantly enriched differed between the two categories. For example, several pathways involved in signal transduction [e.g. Hippo signaling pathway—fly, tumor necrosis factor (TNF) signaling pathway, transforming growth factor (TGF)-β signaling pathway, mTOR signaling pathway, Jak-STAT signaling pathway, AMP-activated protein kinase (AMPK) signaling pathway], amino acid metabolism (e.g. lysine degradation), and signaling molecules and interaction (e.g. neuroactive ligand-receptor interaction) were significantly enriched in the category TZ_24h-vs-C_24h, while some pathways involved in other signal transduction pathways [e.g. mitogen-activated protein kinase (MAPK) signaling pathway, Rap1 signaling pathway, phospholipase D signaling pathway, ErbB signaling pathway], and development and regeneration (e.g. axon guidance, dorso-ventral axis formation) were significantly enriched in the category TZ_48h-vs-C_48h. Of these, most of the targets were predicted to be involved in immune (e.g. Hippo signaling pathway—fly, Fc gamma R-mediated phagocytosis, Toll-like receptor signaling pathway, TNF signaling pathway, TGF-β signaling pathway, mTOR signaling pathway, Jak-STAT signaling pathway, MAPK signaling pathway, ErbB signaling pathway, inflammatory mediator regulation of TRP channels) and metabolic (e.g. lysine degradation, insulin signaling pathway, AMPK signaling pathway, Rap1 signaling pathway, calcium signaling pathway) pathways.
Propagation of N. caninum tachyzoites in caprine EECs affected by chi-miR-146a
The average number of N. caninum tachyzoites per vacuole in caprine EECs was significantly increased by transfection with chi-miR-146a mimics at both 30 hpi (Fig. 6a, b) and 42 hpi (Fig. 6c, d), while the opposite effect was found at these two time points when a chi-miR-146a inhibitor was used during transfection (Fig. 6). These findings indicated that the propagation of N. caninum tachyzoites in caprine EECs was promoted by chi-miR-146a.
Dysregulation of miRNA expression has been reported for several tissues/organs (e.g. ovary , liver , pituitary , lung ) and types of cells (e.g. skeletal muscle satellite cells , skin fibroblast cells , Leydig cells , blood leukocytes and milk somatic cells , intramuscular preadipocytes , endometrial epithelium cells ) in physiological and pathological processes in goats, including in response to infectious diseases. For example, the expression of 316 DEmiRNAs was found in peripheral blood mononuclear cells of goats infected with peste des petits ruminants virus, and most of the predicted targets of these DEmiRNAs were found to be involved in immune escape . Brucella melitensis M5-90 infection induced 777 DEmiRNAs in goat fibroblasts, and targets for these DEmiRNAs were predicted to function in immune responses (e.g. cytokine-cytokine receptor interaction, natural killer cell-mediated cytotoxicity and Toll-like receptor signaling pathway, TNF signaling pathway, MAPK signaling pathway and JAK/STAT signaling pathway) . Additionally, Fasciola gigantica-derived excretory-secretory products were found to induce 30 DEmiRNAs of goat peripheral blood mononuclear cells, and predicted targets for these DEmiRNAs were significantly enriched in biological processes (e.g. cell differentiation, cell development and regulation of nervous system development) and the TGF-β signaling pathway . These findings suggested that miRNAs have significant biological roles in pathophysiological processes.
The uterus is the reproductive organ in which the fetus develops in placental mammals . Using RNA-seq, 578 and 464 miRNAs were identified in the endometrium of Xinong Saanen dairy goats and Chuanzhong black goats, respectively [59, 60]. A total of 143 miRNAs were differentially expressed during embryo implantation in Xinong Saanen dairy goats , and 33 DEmiRNAs were identified in goats on day 16 of pregnancy that were not present in non-pregnant goats on day 16 of the estrous cycle . Additionally, 434 miRNAs were identified in the extracellular vesicles of Chuanzhong black goat uterine fluid, with 106 DEmiRNAs identified during peri-implantation . To elucidate the pathophysiological events that take place in the uterus during infection with N. caninum, in the present study we determined the expression profiles of miRNAs in caprine EECs infected with N. caninum at 24 hpi and 48 hpi. A total of 1239 miRNAs were identified in the caprine EECs, of which 891 are novel miRNAs, and thus our findings greatly enhance the miRNA database with respect to information on goats.
In the present study, N. caninum infection induced 18 and 79 DEmiRNAs at 24 hpi and 48 hpi, respectively. Functional analysis showed that most of the targets of these DEmiRNAs were significantly enriched in immune-related signaling pathways (e.g. Fc gamma R-mediated phagocytosis, Toll-like receptor signaling pathway, TNF signaling pathway, TGF-β signaling pathway, MAPK signaling pathway). Of these DEmiRNAs, downregulated miR-193b-5p has been found to target HMGA2 to inhibit 5-fluorouracil-induced apoptosis through the MAPK signaling pathway . Downregulated miRNA-455-3p promoted TGF-β signaling and inhibited the development of osteoarthritis by targeting PAK2 . Furthermore, downregulated miR-361-3p was predicted to regulate the expression of MAP3K8 through the TNF signaling pathway in our study, and MAP3K8 has been reported to promote angiogenesis and inhibit inflammation with the participation of the tumor necrosis factor-α (TNF-α) signaling pathway . Interestingly, C57BL/6 mice infected with N. caninum exhibited higher mortality associated with inflammatory lesions and increased parasite burden in the brain through TNF-TNFR1 signaling . In addition, in our study, downregulated miR-128-5p was predicted to activate MKNK2 through the MAPK signaling pathway. Previous studies showed that N. caninum evaded antigen presentation from bone marrow-derived macrophages by activating p38 MAPK , and extracellular vesicles secreted by N. caninum regulated the cytokine expression of bone marrow-derived macrophages through Toll-like receptor 2 and MAPK signaling pathways . Furthermore, a 14-3-3 protein of N. caninum could induce effective immune responses and stimulate cytokine expression by activation of the MAPK signaling pathway . These findings suggest that these DEmiRNAs play important roles in immune defense or immunopathogenesis during N. caninum infection.
Additionally, the functional analysis also showed that targets of the DEmiRNAs were significantly enriched in metabolism-related pathways (e.g. lysine degradation, insulin signaling pathway, AMPK signaling pathway, Rap1 signaling pathway, calcium signaling pathway). Among these DEmiRNAs, downregulated miR-128-5p was predicted in our study to regulate the expression of inhibitor of DNA-binding 1 through the Rap1 signaling pathway; inhibitor of DNA-binding 1 has been reported to play a key role in cell adhesion in neural stem cells through the preservation of Rap1 signaling . Downregulated miR-877-3p was predicted in our study to regulate the expression of sirtuin 1 (SIRT1) through the AMPK signaling pathway. A previous study showed that SIRT1 could suppress lipogenesis through the AMPK signaling pathway , and our group found that SIRT1 promoted cell autophagy and intracellular proliferation of N. caninum tachyzoites in caprine EECs through inducing mitochondrial dysfunction . These findings indicate that these DEmiRNAs play a significant role in the interaction between the host and N. caninum through metabolic regulation.
Interestingly, five of the DEmiRNAs (miR-146a, miR-155-5p, miR-155-3p, miR-17-5p, miR-20a-5p) in the present study were also found to be dysregulated during T. gondii infection in previous studies [29, 72,73,74]. Of these, miR-146a, miR-155-5p and miR-155-3p were upregulated during infections with two different protozoan parasites [72,73,74], but the opposite was found for miR-17-5p  and miR-20a-5p expression . Notably, miR-146a has been identified as a microRNA fingerprint associated with Toxoplasma persistence in the host brain . Previous studies showed that miR-146a is a negative regulator of innate immune response , and is associated with inflammatory immune responses (e.g. systemic lupus erythematosus , osteoarthritis , rheumatoid arthritis ), tumors (e.g. colorectal cancer , breast cancer , gastric cancer , lung cancer ), angiogenesis [83, 84], apoptosis , and autophagy . In addition, miR-146a was reported to affect gene expression through various signaling pathways, e.g. TNF-α, NF-κB, MEK-1/2 and JNK-1/2 . In our study, chi-miR-146a was upregulated and promoted the propagation of N. caninum tachyzoites in caprine EECs, suggesting that this miRNA plays a role in the intracellular survival of N. caninum tachyzoites in host cells. Thus, the regulatory mechanism of chi-miR-146a should be further studied in future work.
N. caninum induced significant differential expression of miRNA in caprine EECs at 24 hpi and 48 hpi, and functional analysis showed that the DEmiRNAs may play significant roles in the interaction between N. caninum and caprine EECs. Upregulated chi-miR-146a promoted the propagation of N. caninum tachyzoites in caprine EECs. These findings provide fundamental data for further elucidation of the significance of ncRNAs in host defense and pathogenesis during N. caninum infection.
Availability of data and materials
The datasets supporting the findings of this article are included within the article and its additional files.
AMP-activated protein kinase
Endometrial epithelial cells
Fetal bovine serum
Kyoto Encyclopedia of Genes and Genomes
Mitogen-activated protein kinase
Quantitative real-time polymerase chain reaction
Transforming growth factor
Tumor necrosis factor
Lindsay DS, Dubey JP. Neosporosis, toxoplasmosis, and sarcocystosis in ruminants: an update. Vet Clin North Am Food Anim Pract. 2020;36:205–22.
Almería S. Neospora caninum and wildlife. ISRN Parasitol. 2013;2013:947347.
Ying Z, Zhu ZF, Yang X, Liu J, Liu Q. Prevalence and associated risk factors of Neospora caninum infection among cattle in mainland China: a systematic review and meta-analysis. Prev Vet Med. 2022;201:105593.
Dubey JP. Review of Neospora caninum and neosporosis in animals. Korean J Parasitol. 2003;41:1–16.
Dubey JP. Recent advances in Neospora and neosporosis. Vet Parasitol. 1999;84:349–67.
Dubey JP, Schares G, Ortega-Mora LM. Epidemiology and control of neosporosis and Neospora caninum. Clin Microbiol Rev. 2007;20:323–67.
Rodrigues AA, Reis SS, Sousa ML, Moraes EDS, Garcia JL, Nascimento TVC, et al. A systematic literature review and meta-analysis of risk factors for Neospora caninum seroprevalence in goats. Prev Vet Med. 2020;185:105176.
Oshiro LM, Motta-Castro AR, Freitas SZ, Cunha RC, Dittrich RL, Meirelles AC, et al. Neospora caninum and Toxoplasma gondii serodiagnosis in human immunodeficiency virus carriers. Rev Soc Bras Med Trop. 2015;48:568–72.
Ojo KK, Reid MC, Kallur Siddaramaiah L, Müller J, Winzer P, Zhang Z, et al. Neospora caninum calcium-dependent protein kinase 1 is an effective drug target for neosporosis therapy. PLoS One. 2014;9:e92929.
Nishikawa Y. Towards a preventive strategy for neosporosis: challenges and future perspectives for vaccine development against infection with Neospora caninum. J Vet Med Sci. 2017;79:1374–80.
Monney T, Hemphill A. Vaccines against neosporosis: what can we learn from the past studies? Exp Parasitol. 2014;140:52–70.
Marugan-Hernandez V. Neospora caninum and bovine neosporosis: current vaccine research. J Comp Pathol. 2017;157:193–200.
Dubey JP, Schares G. Neosporosis in animals–the last five years. Vet Parasitol. 2011;180:90–108.
Hall CA, Reichel MP, Ellis JT. Neospora abortions in dairy cattle: diagnosis, mode of transmission and control. Vet Parasitol. 2005;128:231–41.
Reichel MP, Alejandra Ayanegui-Alcérreca M, Gondim LF, Ellis JT. What is the global economic impact of Neospora caninum in cattle—the billion dollar question. Int J Parasitol. 2013;43:133–42.
Lu TX, Rothenberg ME. MicroRNA. J Allergy Clin Immunol. 2018;141:1202–7.
Lai EC. Micro RNAs are complementary to 3’ UTR sequence motifs that mediate negative post-transcriptional regulation. Nat Genet. 2002;30:363–4.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.
O’Brien J, Hayder H, Zayed Y, Peng C. Overview of microRNA biogenesis, mechanisms of actions, and circulation. Front Endocrinol (Lausanne). 2018;9:402.
Sirotkin AV, Lauková M, Ovcharenko D, Brenaut P, Mlyncek M. Identification of microRNAs controlling human ovarian cell proliferation and apoptosis. J Cell Physiol. 2010;223:49–56.
Subramanian S, Steer CJ. MicroRNAs as gatekeepers of apoptosis. J Cell Physiol. 2010;223:289–98.
Babashah S, Soleimani M. The oncogenic and tumour suppressive roles of microRNAs in cancer and apoptosis. Eur J Cancer. 2011;47:1127–37.
Fu LL, Wen X, Bao JK, Liu B. MicroRNA-modulated autophagic signaling networks in cancer. Int J Biochem Cell Biol. 2012;44:733–6.
Zha X, Xi X, Fan X, Ma M, Zhang Y, Yang Y. Overexpression of METTL3 attenuates high-glucose induced RPE cell pyroptosis by regulating miR-25-3p/PTEN/Akt signaling cascade through DGCR8. Aging (Albany NY). 2020;12:8137–50.
Lee YS, Dutta A. MicroRNAs in cancer. Annu Rev Pathol. 2009;4:199–227.
Bhaskaran M, Mohan M. MicroRNAs: history, biogenesis, and their evolving role in animal development and disease. Vet Pathol. 2014;51:759–74.
Li S, Yang J, Wang L, Du F, Zhao J, Fang R. Expression profile of microRNAs in porcine alveolar macrophages after Toxoplasma gondii infection. Parasit Vectors. 2019;12:65.
Hou Z, Liu D, Su S, Wang L, Zhao Z, Ma Y, et al. Comparison of splenocyte microRNA expression profiles of pigs during acute and chronic toxoplasmosis. BMC Genomics. 2019;20:97.
Cong W, Zhang XX, He JJ, Li FC, Elsheikha HM, Zhu XQ. Global miRNA expression profiling of domestic cat livers following acute Toxoplasma gondii infection. Oncotarget. 2017;8:25599–611.
He JJ, Ma J, Wang JL, Xu MJ, Zhu XQ. Analysis of miRNA expression profiling in mouse spleen affected by acute Toxoplasma gondii infection. Infect Genet Evol. 2016;37:137–42.
Dubey JP, Barr BC, Barta JR, Bjerkås I, Björkman C, Blagburn BL, et al. Redescription of Neospora caninum and its differentiation from related coccidia. Int J Parasitol. 2002;32:929–46.
Howe DK, Sibley LD. Comparison of the major antigens of Neospora caninum and Toxoplasma gondii. Int J Parasitol. 1999;29:1489–96.
Al-Bajalan MMM, Xia D, Armstrong S, Randle N, Wastling JM. Toxoplasma gondii and Neospora caninum induce different host cell responses at proteome-wide phosphorylation events; a step forward for uncovering the biological differences between these closely related parasites. Parasitol Res. 2017;116:2707–19.
Zhao SS, Tao DL, Chen JM, Wu JP, Yang X, Song JK, et al. RNA sequencing reveals dynamic expression of lncRNAs and mRNAs in caprine endometrial epithelial cells induced by Neospora caninum infection. Parasit Vectors. 2022;15:297.
Zhao SS, Tao DL, Chen JM, Chen X, Geng XL, Wang JW, et al. Neospora caninum infection activated autophagy of caprine endometrial epithelial cells via mTOR signaling. Vet Parasitol. 2022;304:109685.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. Embnet J. 2011;17:10–2.
Patel RK, Jain M. NGS QC Toolkit: a toolkit for quality control of next generation sequencing data. PLoS One. 2012;7:e30619.
Langmead B. Aligning short sequencing reads with bowtie. Curr Protoc Bioinform. 2010. https://doi.org/10.1002/0471250953.bi1107s32.
Griffiths-Jones S, Bateman A, Marshall M, Khanna A, Eddy SR. Rfam: an RNA family database. Nucleic Acids Res. 2003;31:439–41.
Chen N. Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinform. 2004;5:4–10.
Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ. miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008;36:D154–8.
Friedländer MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40:37–52.
Sun J, Wang S, Li C, Ren Y, Wang J. Novel expression profiles of microRNAs suggest that specific miRNAs regulate gene expression for the sexual maturation of female Schistosoma japonicum after pairing. Parasit Vectors. 2014;7:177.
Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. 2003;5:R1.
Liu Y, Zhou Z, He X, Tao L, Jiang Y, Lan R, et al. Integrated analyses of miRNA-mRNA expression profiles of ovaries reveal the crucial interaction networks that regulate the prolificacy of goats in the follicular phase. BMC Genomics. 2021;22:812.
Zhao X, Ji Z, Xuan R, Wang A, Li Q, Zhao Y, et al. Characterization of the microRNA expression profiles in the goat kid liver. Front Genet. 2022;12:794157.
Ye J, Yao Z, Si W, Gao X, Yang C, Liu Y, et al. Identification and characterization of microRNAs in the pituitary of pubescent goats. Reprod Biol Endocrinol. 2018;16:51.
Zhang W, Jiao Z, Huang H, Wu Y, Wu H, Liu Z, et al. Effects of Pasteurella multocida on histopathology, miRNA and mRNA expression dynamics in lung of goats. Animals (Basel). 2022;12:1529.
Liao R, Lv Y, Dai J, Zhang D, Zhu L, Lin Y. chi-miR-99b-3p regulates the proliferation of goat skeletal muscle satellite cells in vitro by targeting caspase-3 and NCOR1. Animals (Basel). 2022;12:2368.
Pang F, Wang X, Chen Z, Zhang Z, Zhang M, Wang C, et al. Integrated analysis of differentially expressed miRNAs and mRNAs in goat skin fibroblast cells in response to orf virus infection reveals that cfa-let-7a regulates thrombospondin 1 expression. Viruses. 2020;12:118.
An SY, Zhang GM, Liu ZF, Zhou C, Yang PC, Wang F. MiR-1197-3p regulates testosterone secretion in goat Leydig cells via targeting PPARGC1A. Gene. 2019;710:131–9.
Urbańska DM, Jarczak J, Czopowicz M, Kaba J, Horbańczuk K, Bagnicka E. miRNA expression patterns in blood leukocytes and milk somatic cells of goats infected with small ruminant lentivirus (SRLV). Sci Rep. 2022;12:13239.
Du Y, Wang Y, Li Y, Emu Q, Zhu J, Lin Y. miR-214-5p regulating differentiation of intramuscular preadipocytes in goats via targeting KLF12. Front Genet. 2021;12:748629.
Zhang L, Liu XR, Liu JZ, Song YX, Zhou ZQ, Cao BY. miR-182 selectively targets HOXA10 in goat endometrial epithelium cells in vitro. Reprod Domest Anim. 2017;52:1081–92.
Qi X, Wang T, Xue Q, Li Z, Yang B, Wang J. MicroRNA expression profiling of goat peripheral blood mononuclear cells in response to peste des petits ruminants virus infection. Vet Res. 2018;49:62.
Li B, Chen S, Wang C, Chen Q, Man C, An Q, et al. Integrated mRNA-seq and miRNA-seq analysis of goat fibroblasts response to Brucella Melitensis strain M5–90. PeerJ. 2021;9:e11679.
Wang SS, Chen D, He JJ, Zheng WB, Tian AL, Zhao GH, et al. Fasciola gigantica-derived excretory-secretory products alter the expression of mRNAs, miRNAs, lncRNAs, and circRNAs involved in the immune response and metabolism in goat peripheral blood mononuclear cells. Front Immunol. 2021;12:653755.
Turner ML, Healey GD, Sheldon IM. Immunity and inflammation in the uterus. Reprod Domest Anim. 2012;47:402–9.
Song Y, An X, Zhang L, Fu M, Peng J, Han P, et al. Identification and profiling of microRNAs in goat endometrium during embryo implantation. PLoS ONE. 2015;10:e0122202.
Zang X, Zhou C, Wang W, Gan J, Li Y, Liu D, et al. Differential microRNA expression involved in endometrial receptivity of goats. Biomolecules. 2021;11:472.
Xie Y, Liu G, Zang X, Hu Q, Zhou C, Li Y, et al. Differential expression pattern of goat uterine fluids extracellular vesicles miRNAs during peri-implantation. Cells. 2021;10:2308.
Wang J, Zhang X, Zhang J, Chen S, Zhu J, Wang X. Long noncoding RNA CRART16 confers 5-FU resistance in colorectal cancer cells by sponging miR-193b-5p. Cancer Cell Int. 2021;21:638.
Hu S, Zhao X, Mao G, Zhang Z, Wen X, Zhang C, et al. MicroRNA-455-3p promotes TGF-β signaling and inhibits osteoarthritis development by directly targeting PAK2. Exp Mol Med. 2019;51:1–13.
Li J, Lv H, Che Y. microRNA-381-3p confers protection against ischemic stroke through promoting angiogenesis and inhibiting inflammation by suppressing cebpb and Map3k8. Cell Mol Neurobiol. 2020;40:1307–19.
Ferreira França FB, Silva MV, Silva MF, Ramos ELP, Miranda VDS, Mota CM, et al. TNF-TNFR1 signaling enhances the protection against Neospora caninum infection. Front Cell Infect Microbiol. 2022;11:789398.
Mota CM, Oliveira AC, Davoli-Ferreira M, Silva MV, Santiago FM, Nadipuram SM, et al. Neospora caninum activates p38 MAPK as an evasion mechanism against innate immunity. Front Microbiol. 2016;7:1456.
Li S, Gong P, Tai L, Li X, Wang X, Zhao C, et al. Extracellular vesicles secreted by Neospora caninum are recognized by Toll-Like receptor 2 and modulate host cell innate immunity through the MAPK signaling pathway. Front Immunol. 2018;9:1633.
Li S, Gong P, Zhang N, Li X, Tai L, Wang X, et al. 14-3-3 protein of Neospora caninum modulates host cell innate immunity through the activation of MAPK and NF-κB pathways. Front Microbiol. 2019;10:37.
Tan R, Lee YJ, Chen X. Id-1 plays a key role in cell adhesion in neural stem cells through the preservation of RAP1 signaling. Cell Adh Migr. 2012;6:1–3.
Long JK, Dai W, Zheng YW, Zhao SP. miR-122 promotes hepatic lipogenesis via inhibiting the LKB1/AMPK pathway by targeting Sirt1 in non-alcoholic fatty liver disease. Mol Med. 2019;25:26.
Tao DL, Zhao SS, Chen JM, Chen X, Yang X, Song JK, et al. Neospora caninum infection induced mitochondrial dysfunction in caprine endometrial epithelial cells via downregulating SIRT1. Parasit Vectors. 2022;15:274.
Cannella D, Brenier-Pinchart MP, Braun L, van Rooyen JM, Bougdour A, Bastien O, et al. miR-146a and miR-155 delineate a microRNA fingerprint associated with Toxoplasma persistence in the host brain. Cell Rep. 2014;6:928–37.
Jiang D, Wu S, Xu L, Xie G, Li D, Peng H. Anti-infection roles of miR-155-5p packaged in exosomes secreted by dendritic cells infected with Toxoplasma gondii. Parasit Vectors. 2022;15:3.
Zhu S, Lu J, Lin Z, Abuzeid AMI, Chen X, Zhuang T, et al. Anti-tumoral effect and action mechanism of exosomes derived from Toxoplasma gondii-infected dendritic cells in mice colorectal cancer. Front Oncol. 2022;12:870528.
Taganov KD, Boldin MP, Chang KJ, Baltimore D. NF-kappaB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. Proc Natl Acad Sci USA. 2006;103:12481–6.
El-Akhras BA, Ali YBM, El-Masry SA, Bassyouni IH, El-Sayed IH, Talaat RM. mir-146a genetic polymorphisms in systemic lupus erythematosus patients: correlation with disease manifestations. Noncoding RNA Res. 2022;7:142–9.
Liu JN, Lu S, Fu CM. MiR-146a expression profiles in osteoarthritis in different tissue sources: a meta-analysis of observational studies. J Orthop Surg Res. 2022;17:148.
Pauley KM, Satoh M, Chan AL, Bubb MR, Reeves WH, Chan EK. Upregulated miR-146a expression in peripheral blood mononuclear cells from rheumatoid arthritis patients. Arthritis Res Ther. 2008;10:R101.
Garo LP, Ajay AK, Fujiwara M, Gabriely G, Raheja R, Kuhn C, et al. MicroRNA-146a limits tumorigenic inflammation in colorectal cancer. Nat Commun. 2021;12:2419.
Li Y, Li W, Lin J, Lv C, Qiao G. miR-146a enhances the sensitivity of breast cancer cells to paclitaxel by downregulating IRAK1. Cancer Biother Radiopharm. 2022;37:624–35.
Shomali N, Mansoori B, Mohammadi A, Shirafkan N, Ghasabi M, Baradaran B. MiR-146a functions as a small silent player in gastric cancer. Biomed Pharmacother. 2017;96:238–45.
Wani JA, Majid S, Khan A, Arafah A, Ahmad A, Jan BL, et al. Clinico-pathological importance of miR-146a in lung cancer. Diagnostics (Basel). 2021;11:274.
Rau CS, Yang JC, Chen YC, Wu CJ, Lu TH, Tzeng SL, et al. Lipopolysaccharide-induced microRNA-146a targets CARD10 and regulates angiogenesis in human umbilical vein endothelial cells. Toxicol Sci. 2014;140:315–26.
Li Y, Zhu H, Wei X, Li H, Yu Z, Zhang H, et al. LPS induces HUVEC angiogenesis in vitro through miR-146a-mediated TGF-β1 inhibition. Am J Transl Res. 2017;9:591–600.
Hsieh YT, Chou YC, Kuo PY, Tsai HW, Yen YT, Shiau AL, et al. Down-regulated miR-146a expression with increased neutrophil extracellular traps and apoptosis formation in autoimmune-mediated diffuse alveolar hemorrhage. J Biomed Sci. 2022;29:62.
Zhang F, Wang J, Chu J, Yang C, Xiao H, Zhao C, et al. MicroRNA-146a induced by hypoxia promotes chondrocyte autophagy through Bcl-2. Cell Physiol Biochem. 2015;37:1442–53.
The authors are grateful to Prof. Yaping Jin and Prof. Xuefeng Qi from Northwest A&F University for providing the passage cells, and to Prof. Qun Liu from China Agricultural University for providing the NC-1 N. caninum tachyzoites.
This work was supported by grants from the Key Research and Development Program of Shaanxi Province (no. 2022NY-097) and the Innovation Support Plan of Shaanxi Province (no. 2021TD-31).
Ethics approval and consent to participate
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The sequences and annealing temperatures of the primers for the quantitative real-time polymerase chain reaction of selected differentially expressed (DE) microRNAs (miRNAs).
The statistics of the reference genomes matched to Capra hircus.
The sequences and lengths of the identified miRNAs.
All the DEmiRNAs.
The predicted target mRNAs for the DEmiRNAs.
Gene Ontology (GO) enrichment analysis of targets for all the DEmiRNAs.
Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of targets for all the DEmiRNAs.
About this article
Cite this article
Chen, JM., Zhao, SS., Tao, DL. et al. Temporal transcriptomic changes in microRNAs involved in the host immune response and metabolism during Neospora caninum infection. Parasites Vectors 16, 28 (2023). https://doi.org/10.1186/s13071-023-05665-9