Impact Factor 2.245 | CiteScore 2.6
More on impact ›

Original Research ARTICLE

Front. Vet. Sci., 12 November 2020 | https://doi.org/10.3389/fvets.2020.587273

Toxocara canis Differentially Affects Hepatic MicroRNA Expression in Beagle Dogs at Different Stages of Infection

  • 1Heilongjiang Key Laboratory for Zoonosis, College of Veterinary Medicine, Northeast Agricultural University, Harbin, China
  • 2State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Lanzhou, China
  • 3Faculty of Medicine and Health Sciences, School of Veterinary Medicine and Science, University of Nottingham, Loughborough, United Kingdom
  • 4College of Veterinary Medicine, Shanxi Agricultural University, Taigu, China

Toxocara canis is a neglected zoonotic parasite, which threatens the health of dogs and humans worldwide. The molecular mechanisms that underlie the progression of T. canis infection remain mostly unknown. MicroRNAs (miRNAs) are small non-coding RNAs that have been identified in T. canis; however, the regulation and role of miRNAs in the host during infection remain incompletely understood. In this study, we determined hepatic miRNA expression at different stages of T. canis infection in beagle dogs. Individual dogs were infected by 300 embryonated T. canis eggs, and their livers were collected at 12 hpi (hours post-infection), 24 hpi, and 36 dpi (days post-infection). The expression profiles of liver miRNAs were determined using RNA-sequencing. Compared to the control groups, 9, 16, and 34 differentially expressed miRNAs (DEmiRNAs) were detected in the livers of infected dogs at the three infection stages, respectively. Among those DEmiRNAs, the novel-294 and cfa-miR-885 were predicted to regulate inflammation-related genes at the initial stage of infection (12 hpi). The cfa-miR-1839 was predicted to regulate the target gene TRIM71, which may influence the development of T. canis larvae at 24 hpi. Moreover, cfa-miR-370 and cfa-miR-133c were associated with immune response at the final stage of infection (36 dpi). Some immunity-related Gene Ontology terms were enriched particularly at 24 hpi. Likewise, Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that many significantly enriched pathways were involved in inflammation and immune responses. The expression level of nine DEmiRNAs was validated using quantitative real-time PCR (qRT-PCR). These results show that miRNAs play critical roles in the pathogenesis of T. canis during the hepatic phase of parasite development. Our data provide fundamental information for further investigation of the roles of miRNAs in the innate/adaptive immune response of dogs infected by T. canis.

Introduction

Toxocariasis, mainly caused by Toxocara canis infection in dogs (1), is a neglected zoonosis worldwide. Dogs, as the definitive hosts of T. canis, can excrete eggs with feces, leading to environment contamination (2, 3). Humans can be infected by T. canis via ingestion of food contaminated with embryonated eggs or larvae (2, 4). The larvae hatch inside the intestinal tract of the host and then migrate to the other parts of the body, leading to serious health problems, such as visceral larva migrans, ocular larva migrans, and neurotoxocariasis (5, 6). Owing to non-specific symptoms in human infection and the diagnostic challenges (7), the public health impact caused by T. canis infection may be ignored. More efforts are needed to achieve better understanding of the pathogenesis of toxocariasis (1).

In recent years, multiple approaches, such as genomics, transcriptomics, and proteomics, have been applied to characterize the biological and molecular features of T. canis (811). However, there is a dearth of information regarding the interaction between T. canis and its host. MicroRNAs (miRNAs) are endogenous, small, non-coding RNAs that have received significant attention from the scientific community due to their involvement in many and diverse biological processes, such as apoptosis, proliferation, metabolism, and immune response (12). miRNAs together with transcription factors have been considered key regulators that modulate the expression levels of almost all genes that mediate various pathophysiological processes (13, 14). Also, parasites can alter the expression level of host miRNAs in order to regulate gene expression of the target tissues (1517). miRNAs can repress mRNA expression through binding to the 3′ untranslated regions of target genes or by enhancing mRNA degradation (18). A previous study showed that T. canis miRNAs, Tc-let-7-5p, Tc-miR-34, and Tc-miR-100, play roles in host–parasite interactions (8). However, the regulation of miRNAs and their roles in the pathogenesis of T. canis during the hepatic phase of infection remains mostly unknown.

Therefore, the present study aimed to investigate the alteration of miRNA expression profiles in the livers of beagle dogs infected by T. canis at different stages of infection using small RNA transcriptome sequencing and bioinformatics analysis. The study findings revealed that some differentially expressed miRNAs (DEmiRNAs) play roles in the regulation of inflammatory and immune responses of puppies against T. canis infection.

Materials and Methods

Experimental Infection of Dogs

The adult T. canis worms were collected from naturally infected dogs in Rongchang District, Chongqing Municipality, China, and female adult T. canis were identified based on gross morphology. The eggs of T. canis were obtained from the uteri of female T. canis worms. The unembryonated eggs were incubated with 0.5% formalin solution at 28°C (85–95% humidity) for 28 days, and then the eggs were harvested and filtered through a 200-mesh screen. The finally embryonated eggs were stored in 1% formalin solution at 4°C. Eighteen beagle puppies (6–7 weeks old) were provided by the National Canine Laboratory Animal Resource Center and housed following Good Laboratory Practice (GLP) in an animal facility, according to the GB standard (GB 14922.2-2011) of China. All puppies had tested negative for T. canis infection by indirect ELISA, using larval ES antigen, before the start of the experiment. Feces of individual puppies were collected daily and examined using a standard sugar floatation method. Puppies were equally allocated into three experimental and three control groups (three puppies per group). The puppy groups were housed separately and were provided with similar conditions including access to the same water source and food supply. Each puppy was orally infected with 300 embryonated eggs in 1 ml normal saline solution, while the control puppies were inoculated with the same amount of saline but without any eggs.

Detection of T. canis Infection

Blood samples of each puppy were collected from the jugular vein into sterile tubes containing EDTA-K2 and tubes without anticoagulant. The T. canis IgG antibody was detected using indirect ELISA as previously described (19). Light microscopy was also used to observe whether T. canis larvae are present in the livers of puppies. The genomic DNA of infectious eggs (used to infect puppies), larvae (isolated from the liver of infected puppies), and adult T. canis (recovered from the small intestine of infected puppies) were isolated using a DNA extraction kit (TianGen™, Beijing, China) according to the manufacturer's instruction. The isolated DNA samples were analyzed using PCR with primer sequences specific to the internal transcribed spacer (ITS) region (partial sequence of ITS-1 and ITS-2) (20). The positive PCR products were sequenced by Sangon Biotech (Shanghai, China), and the obtained sequences were searched against similar sequences available in the GenBank database using Nucleotide BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi).

Sample Collection, RNA Extraction, and Quantification

Infected puppies (n = 9) and naive (control) puppies (n = 9) were killed by potassium chloride (KCl) under a general anesthetic [50 mg/kg tiletamine-zolazepam (Zoletil®), Virbac, France]. The liver samples were collected from all puppies at 12 h post-infection (hpi), 24 hpi, and 36 days post-infection (dpi) according to methods described in our previous study (21). Three biological replicates were examined in each group at each time point post-infection. The collected liver samples were quickly stored in liquid nitrogen until used for RNA extraction. The larvae were recovered from the remaining fresh liver sample of each puppy using the modified Baermann funnel method as previously described (22). The total RNA was extracted from the liver samples of puppies using TRIZOL (Life Technologies, CA, USA). The genomic DNA was removed using DNase I (NEB, Ipswich, USA). RNA concentration was measured using the Qubit® RNA Assay Kit and Qubit® 2.0 Fluorometer (Life Technologies, CA, USA). The integrity of RNA was assessed by the RNA Nano 6000 Assay Kit and the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA). The purity of RNA was examined using the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA). High-quality RNA samples with RNA integrity numbers (RINs) >8.0 were used to construct the sequencing library.

Small RNA Library Preparation and Sequencing

A total of 3 μg RNA of each sample was used for the construction of the small RNA library by using the NEBNext® Multiplex Small RNA Library Prep Set for Illumina® (NEB, USA). The constructed libraries were sequenced on an Illumina Hiseq 2500 platform. The reads that contain poly N, 5′ adapter contaminants, without 3′ adapter or the insert tag (including poly G, C, A, or T and low-quality reads) were filtered from raw data using custom perl and python scripts. The clean reads that ranged from 18 to 35 nt were mapped against the reference sequence by Bowtie (23). The mapped sRNAs were searched against miRBase20.0 data to identify known miRNAs. In addition, mirdeep2 (24) and miREvo (25) were used to identify potential novel miRNAs. The novel miRNA was predicted by the characteristics of the hairpin structure of the miRNA precursor. The expression levels of miRNA were estimated by TPM (transcript per million) (26). Differential expression analysis was performed using the DESeq R package (1.8.3). A P-value < 0.05 was used as the significance threshold value of differential expression.

Bioinformatics Analysis of DEmiRNA

The target genes of miRNAs were predicted by RNAhybrid, PITA tools, and miRanda (27). To predict the function of the target gene of DEmiRNAs, GOseq R package (28), and KOBAS software (29, 30) were used for Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, respectively. A P-value < 0.05 was considered as significantly enriched.

Verification of miRNA Expression by Quantitative Real-Time PCR

Quantitative real-time PCR (qRT-PCR) was performed to confirm the upregulation and downregulation of nine miRNAs randomly selected for verification of the RNA-seq results. qRT-PCR was performed using an miRcute enhanced miRNA qRT-PCR Kit (TianGen, Beijing, China) on a LightCycler480 (Roche, Basle, Switzerland). The cDNA of miRNA was synthesized using an miRcute enhanced miRNA cDNA first chain synthesis Kit (TianGen, Beijing, China). Then, 1 μl cDNA was used for qRT-PCR as follows: initial denaturation at 95°C for 15 min, followed by 40 cycles of 94°C for 20 s and 60°C for 34 s. All primers are listed in Table 1. The U6 small nuclear RNA (snRNA) was used as an internal control gene. Melting curve analysis (95°C for 10 s, 65°C for 1 min, and progressive increase from 65°C to 95°C) was performed to ensure specific amplification in each reaction. Each reaction included a blank control to rule out the presence of contamination. The relative expression quantity was analyzed by the 2−ΔΔCt method (31, 32).

TABLE 1
www.frontiersin.org

Table 1. Primers used in microRNA (miRNA)–specific quantitative real-time PCR (qRT-PCR) analysis.

Results

Detection of T. canis Infection in Beagle Dogs

A blood sample of each puppy was collected and tested for IgG antibody against T. canis. The IgG antibodies of T. canis were detected in infected puppies at 36 dpi. At 12 hpi, T. canis larvae were recovered from the livers of three infected puppies, and at 24 hpi, T. canis larvae were recovered from the livers of all infected puppies. At 36 dpi, T. canis larvae were detected in the liver of one infected puppy, and adult T. canis were recovered in the small intestine of all infected puppies. Furthermore, no T. canis larvae and anti–T. canis IgG antibodies were found in the control puppies. The sequences obtained from the embryonated eggs, larvae, and adult worms were found to match the sequence of T. canis (GenBank Accession No. JF837169.1).

Characteristics of the Sequenced Data

In each miRNA library, 142,605,411 raw reads and 7.131 Gb raw data were obtained from infected puppy groups, whereas 136,764,149 raw reads and 6.898 Gb raw data were obtained from the control groups. More than 99% of reads had sequencing quality >Q20 (Table 2). The clean reads with appropriate 18–35 nt lengths were selected for further analysis. Moreover, 52.72–65.17 and 0.17–0.24% unique reads were confirmed as known and novel miRNAs, respectively. In addition, 22.16–29.70% non-annotated reads were found in this study (Supplementary Table 1), suggesting that these non- annotated reads could be involved in the pathogenesis and progression of T. canis infection.

TABLE 2
www.frontiersin.org

Table 2. Summary of the quality control parameters of the reads.

Differentially Expressed Hepatic miRNAs (DEmiRNAs) at Different Infection stages

A total of 59 DEmiRNAs were identified at three infection stages, including 9, 16, and 34 DEmiRNAs at 12 hpi, 24 hpi, and 36 dpi, respectively. Details of DEmiRNAs are shown in Supplementary Table 2. Among these DEmiRNAs, 23 DEmiRNAs were upregulated, whereas 36 DEmiRNAs were downregulated (Figure 1). However, no common miRNA was found at the three infection stages (Figure 2).

FIGURE 1
www.frontiersin.org

Figure 1. The numbers of differentially expressed microRNAs (DEmiRNAs) at three infection stages [12 h post-infection (hpi), 24 hpi, and 36 days post-infection (dpi)]. The red and blue colors represent the upregulated and downregulated miRNAs, respectively.

FIGURE 2
www.frontiersin.org

Figure 2. Venn diagram showing the shared and unique DEmiRNAs at different time points after infection (12 hpi, 24 hpi, and 36 dpi).

Target Gene Prediction and Functional Analysis

Candidate target genes of DEmiRNAs of beagle dog livers were predicted by RNAhybrid, miRanda, and PITA tools. The prediction results showed that a total of 384, 518, and 1,225 hepatic genes were targeted by DEmiRNAs at 12 hpi, 24 hpi, and 36 dpi, respectively. The immunoglobulin superfamily (IgSF) is a molecular superfamily with an immunoglobulin-like domain, and most members of IgSF are found on the surface of lymphocytes and participate in various immune activities. In this study, the IgSF genes were predicted to be regulated by some miRNAs, such as novel-294, cfa-miR-25, cfa-miR-15b, cfa-miR-145, cfa-miR-150, cfa-miR-497, cfa-miR-1839 and cfa-miR-151, cfa-miR-133c, novel-337, cfa-miR-127, cfa-miR-205, and cfa-miR-194 (Supplementary Table 3). The GO enrichment analysis showed that 9 DEmiRNAs were significantly enriched in 797 GO terms at 12 hpi; 16 DEmiRNAs were significantly enriched in 768 GO terms at 24 hpi; and 34 DEmiRNAs were significantly enriched in 665 GO terms at 36 dpi (P < 0.05) (Supplementary Table 4). The top 30 most significant GO terms (P < 0.05) belonging to biological process and molecular function at each time point are shown in Figure 3. Among of these GO terms, the immune-related GO terms were found in infected livers (Supplementary Table 5). At 12 hpi, one DEmiRNA was related to three immune-related GO terms, including natural killer (NK) cell differentiation involved in immune response (GO:0002325), regulation of NK cell differentiation involved in immune response (GO:0032826), and negative regulation of NK cell differentiation involved in immune response (GO:0032827). At 24 hpi, six DEmiRNAs were significantly related to 32 immune-related GO terms, including positive regulation of B-cell–mediated immunity (GO:0002714), positive regulation of immunoglobulin-mediated immune response (GO:0002891), positive regulation of leukocyte-mediated immunity (GO:0002705), positive regulation of lymphocyte-mediated immunity (GO:0002708), and regulation of innate immune response (GO:0045088). At 36 dpi, 12 DEmiRNAs were significantly related to four immune-related GO terms, including immune system development (GO:0002520), negative regulation of innate immune response (GO:0045824), positive regulation of immune system process (GO:0002684), and negative regulation of immune response (GO:0050777) (Supplementary Table 5).

FIGURE 3
www.frontiersin.org

Figure 3. The top 30 Gene Ontology (GO) terms of the DEmiRNAs. (A) The significantly enriched biological process (BP) and molecular function (MF) terms of the target genes of DEmiRNAs at 12 hpi. (B) The significantly enriched BP terms of target genes of DEmiRNAs at 24 hpi. (C) The significantly enriched BP and MF terms of target genes of DEmiRNAs at 36 dpi.

KEGG Pathway Analysis

KEGG pathway analysis revealed that 5, 10, and 6 pathways were significantly enriched (P < 0.05) at 12 hpi, 24 hpi, and 36 dpi, respectively (Supplementary Table 6). The top 20 targeted pathways of DEmiRNAs are shown in Figure 4. At 12 hpi, some immune-related pathways were significantly enriched, including the mitogen-activated protein kinase (MAPK) signaling pathway (cfa04010) and Fc epsilon RI signaling pathway (cfa04664) (Figure 4A). The highly enriched pathways at 24 hpi (Figure 4B) included glycerolipid metabolism (cfa00561), the renin–angiotensin system pathway (cfa04614), and other types of O-glycan biosynthesis (cfa00514). Moreover, significantly inflammation-related pathways were found at 36 dpi, including cell adhesion molecules (CAMs) (cfa04514), vitamin B6 metabolism (cfa00750), cytokine–cytokine receptor interaction (cfa04060), and glycosaminoglycan biosynthesis–keratan sulfate (cfa00533) (Figure 4C).

FIGURE 4
www.frontiersin.org

Figure 4. The top 20 signaling pathways enriched by the target genes of differently expressed microRNAs (miRNAs) at (A) 12 hpi, (B) 24 hpi, and (C) 36 dpi. The X-axis labels denote the gene numbers, and the Y-axis labels represent the names of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The pathways enclosed within the blue rectangle represent the significantly enriched pathways.

qRT-PCR Validation of RNA-Sequencing Data

Nine miRNAs were randomly selected for qRT-PCR verification (Figure 5), including cfa-miR-381, cfa-miR-10b, cfa-miR-146, cfa-miR-10a, cfa-miR-194, cfa-miR-30a, cfa-miR-125a, cfa-miR-371, and cfa-miR-16. Although the expression levels of the nine miRNAs obtained by qRT-PCR were slightly higher those obtained by RNA-sequencing, the expression trends obtained by both methods were consistent, which suggested an increased expression of the examined miRNAs.

FIGURE 5
www.frontiersin.org

Figure 5. Quantitative real-time PCR. The Y-axis shows the relative change of miRNA levels expressed as fold increase compared with the control U6 small nuclear RNA (snRNA). The X-axis shows the name of the nine miRNAs used in the cross-validation of expression analysis. *P < 0.05; **P < 0.01; ***P < 0.001.

Discussion

In this study, hepatic miRNA expression patterns of beagle dogs during T. canis infection were investigated using RNA-sequencing analysis. A total of 59 DEmiRNAs were identified, including 9, 16, and 34 DEmiRNAs at 12 hpi, 24 hpi, and 36 dpi, respectively. Most of the downregulated DEmiRNAs were found in the infected group at 36 dpi (Figure 1). This finding suggested that miRNAs were negatively regulated to alter the expression of the target genes during the later stage of infection in infected puppies' livers. Moreover, among the dysregulated miRNAs, only one and three were common between 12 hpi and 36 dpi and between 24 hpi and 36 dpi, respectively (Figure 2). This finding indicates that the expression of hepatic miRNAs of beagle dogs evolves during the T. canis infection. Two important miRNAs, novel-294 and cfa-miR-885, were differentially expressed at the initial stage of infection (12 hpi) (Supplementary Table 2). Furthermore, the IgSF gene FGFR1 was predicted to be regulated by miRNA novel-294 at 12 hpi (Supplementary Table 3). A previous study found that FGFR1 expression improved beta-cell survival in cytokine-induced inflammation (33). Hence, we speculated that the upregulation of the novel-294 in puppies' livers promotes the immune response of puppies to T. canis infection via increasing the expression of FGFR1. miRNA cfa-miR-885 was predicted to regulate IGSF3 gene of the IgSF (Supplementary Table 3). The IgSF and leukocyte integrins play an important role in the regulation of leukocyte recruitment to the inflammation sites (34). These findings suggest that cfa-miR-885 plays a protective role during T. canis infection.

At 24 hpi, 16 miRNAs were found differentially expressed in infected livers (e.g., cfa-miR-1839). The cfa-miR-1839 was downregulated at 24 hpi, and TRIM71, the predicted target gene of cfa-miR-1839 (Supplementary Table 3), regulates juvenile-to-adult transition events in nematodes and mammals (35). The majority of T. canis larvae reach the liver at 24 hpi (36); subsequently, the larvae migrate to lung, muscle, and brain tissue via the circulation (36). However, some larvae cannot continue to migrate and are trapped in the hepatic capillaries (36). Thus, we assumed that the downregulation of hepatic cfa-miR-1839 possibly regulates the juvenile-to-adult development of T. canis larvae in the puppies' livers. Whether cfa-miR-1839 regulates juvenile-to-adult transition of T. canis larvae remains to be further investigated.

cfa-miR-370 and cfa-miR-133c were significantly downregulated at 36 dpi (Supplementary Table 2). The target gene prediction showed that cfa-miR-370 targets the CD3E gene (Supplementary Table 3). CD3E, the T-cell antigen receptor epsilon subunit gene, is essential for TCR signaling and T-cell differentiation (37, 38). In our study, cfa-miR-370 was downregulated in the infected livers at 36 dpi, indicating that the downregulation of cfa-miR-370 may represent an immunoreaction to resist T. canis infection. In addition, two differently expressed miRNAs, cfa-miR-370 and cfa-miR-133c, were predicted to regulate the ARRDC1 gene (Supplementary Table 3). Additionally, cfa-miR-133c was predicted to regulate the IL18R1 gene. Previous research found that IL-18R/MyD88 plays a crucial role in the development of a robust Th1 response during Trypanosoma cruzi infection (39). Therefore, the abnormal expression of cfa-miR-133c seems to be a potential candidate for further study of the role of cfa-miR-133c in promoting a protective Th1 immune response to T. canis infection.

GO enrichment analysis showed that 797, 769, and 665 GO terms were significantly enriched at 12 hpi, 24 hpi, and 36 dpi, respectively (Supplementary Table 4). Based on the analysis of the 797 significantly enriched GO terms at 12 hpi, nine DEmiRNAs were mainly involved in the signal-related biological process, including instance single-organism processes, single-organism cellular process, and single-organism metabolic process (Figure 3A). These processes are involved in acute liver injury (40). At 24 hpi, 16 DEmiRNAs were mainly associated with biological processes, such as protein phosphorylation, intracellular signal transduction, regulation of phosphorylation, and regulation of protein phosphorylation (Figure 3B). Some phosphorylated proteins were related to TNF-α signaling (41), which regulate immune response (42). At 36 dpi, 10 upregulated and 24 downregulated miRNAs were involved in response to stress, positive regulation of the immune system process, and regulation of nervous system development (Figure 3C). Moreover, a total of 39 immunization-related GO terms were significantly enriched at three time points (Supplementary Table 5). These results show that T. canis larvae can elicit a significant immune response after infecting the liver.

The top 20 enriched pathways are shown in Figure 4. The target genes of the DEmiRNAs were significantly enriched in the MAPK signaling pathway, galactose metabolism, and Fc epsilon RI signaling pathway at 12 hpi (Figure 4A). The MAPK signaling pathway participated in diverse cellular responses, such as inflammation, differentiation, proliferation, and apoptosis (43). Additionally, the MAPK signaling pathway is involved in the development, reproduction, and survival of Schistosoma mansoni (44). The Fc epsilon RI is the primary receptor in mast cells, which can influence the synthesis of proinflammatory cytokines and other molecules that are involved in inflammatory responses (45) (Figure 4A). These findings suggested that T. canis can cause inflammatory responses in puppy livers at the initial stage. At 24 hpi, the glycerolipid metabolism pathway, renin–angiotensin system pathway, and other types of O-glycan biosynthesis pathways were significantly enriched (Figure 4B). The glycerolipid metabolism pathway could perturb the host immune system and metabolism following anisakid infection (46). Further, the renin–angiotensin system pathway influences a range of processes from inflammation and immune responses to longevity (47). A previous study found that the O-glycoprotein biosynthesis pathway was associated with the immune response of dendritic cells (48). According to these findings, we assumed that T. canis larvae can trigger the immune responses in the infected host liver at 24 hpi. Some inflammation-related pathways were significantly enriched at 36 dpi, such as the cell adhesion molecules (CAMs) pathway, vitamin B6 metabolism pathway, cytokine–cytokine receptor interaction pathway, and glycosaminoglycan biosynthesis–keratan sulfate pathway (Figure 4C). The CAMs can direct mediate leukocyte migration, which is essential for generating effective inflammatory responses (49). Previous research found that vitamin B6 as a co-factor was involved in the inflammation response (50). It was reported that a high level of cytokine–cytokine motif chemokine ligand 1 could enhance and prolong the inflammatory response (51). The glycosaminoglycan pathway interacts with multiple ligands, which play an essential role in the inflammatory reaction (52). These results indicated that the puppies' livers may trigger inflammatory responses at the late stage of T. canis infection.

Conclusion

The present study, for the first time, revealed the hepatic miRNA expression patterns of beagle dogs at three T. canis infection stages. A total of 59 DEmiRNAs were identified in the infected livers. The functional enrichment analysis of predicted target genes showed that miRNA cfa-miR-1839 could be related to the juvenile-to-adult transition of T. canis larvae. KEGG pathway analysis found that some significantly enriched pathways were related to the inflammatory response at 12 hpi and 36 dpi. Several miRNAs, such as novel-294, cfa-miR-88, cfa-miR-370, and cfa-miR-133c, were associated with immune responses. These findings should enrich our understanding of the interactions between T. canis and its definitive host. Further studies to elucidate the detailed molecular mechanisms and the physiological functions of the DEmiRNAs in the pathogenesis of T. canis infection are warranted.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA630302.

Ethics Statement

The animal study was reviewed and approved by the Animal Administration and Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences (Approval No. 2018-015). The dogs used in the study were handled in accordance with good animal practices required by the Animal Ethics Procedures and Guidelines of the People's Republic of China.

Author Contributions

Y-XL and X-QZ conceived and designed the experiments. YZ and W-BZ performed the experiments. YZ analyzed the data and wrote the paper. J-JH and HME participated in improving the English of the manuscript. HME, X-QZ, and Y-XL critically revised the manuscript. All authors have read and approved the final version of the manuscript. All authors contributed to the preparation of the manuscript.

Funding

Project support was provided by the International Science and Technology Cooperation Project of Gansu Provincial Key Research and Development Program (Grant No. 17JR7WA031) and the Agricultural Science and Technology Innovation Program (ASTIP) (Grant No. CAAS-ASTIP-2016-LVRI-03).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thanked Novogene Bioinformatics Technology Co., Ltd. (Beijing, China) for technical support.

Supplementary Material

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

References

1. Chen J, Zhou DH, Nisbet AJ, Xu MJ, Huang SY, Li MW, et al. Advances in molecular identification, taxonomy, genetic variation and diagnosis of Toxocara spp. Infect Genet Evol. (2012) 12:1344–8. doi: 10.1016/j.meegid.2012.04.019

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Overgaauw PA, van Knapen F. Veterinary and public health aspects of Toxocara spp. Vet Parasitol. (2013) 193:398–403. doi: 10.1016/j.vetpar.2012.12.035

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Holland CV. Knowledge gaps in the epidemiology of Toxocara: the enigma remains. Parasitology. (2017) 144:81–94. doi: 10.1017/S0031182015001407

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Macpherson CN. The epidemiology and public health importance of toxocariasis: a zoonosis of global importance. Int J Parasitol. (2013) 43:999–1008. doi: 10.1016/j.ijpara.2013.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Ma G, Holland CV, Wang T, Hofmann A, Fan CK, Maizels RM, et al. Human toxocariasis. Lancet Infect Dis. (2018) 18:14–24. doi: 10.1016/S1473-3099(17)30331-6

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Fan CK, Liao CW, Cheng YC. Factors affecting disease manifestation of toxocarosis in humans: genetics and environment. Vet Parasitol. (2013) 193:342–52. doi: 10.1016/j.vetpar.2012.12.030

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Rubinsky-Elefant G, Hirata CE, Yamamoto JH, Ferreira MU. Human toxocariasis: diagnosis, worldwide seroprevalences and clinical expression of the systemic and ocular forms. Ann Trop Med Parasitol. (2010) 104:3–23. doi: 10.1179/136485910X12607012373957

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Ma G, Luo Y, Zhu H, Luo Y, Korhonen PK, Young ND, et al. MicroRNAs of Toxocara canis and their predicted functional roles. Parasit Vectors. (2016) 9:229. doi: 10.1186/s13071-016-1508-3

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Sperotto RL, Kremer FS, Aires Berne ME, Costa de Avila LF, da Silva Pinto L, Monteiro KM, et al. Proteomic analysis of Toxocara canis excretory and secretory (TES) proteins. Mol Biochem Parasitol. (2017) 211:39–47. doi: 10.1016/j.molbiopara.2016.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zheng WB, Zou Y, Zhu XQ, Liu GH. Toxocara “omics” and the promises it holds for medicine and veterinary medicine. Adv Parasitol. (2020) 109:89–108. doi: 10.1016/bs.apar.2020.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Zhu XQ, Korhonen PK, Cai H, Young ND, Nejsum P, von Samson-Himmelstjerna G, et al. Genetic blueprint of the zoonotic pathogen Toxocara canis. Nat Commun. (2015) 6:6145. doi: 10.1038/ncomms7145

PubMed Abstract | CrossRef Full Text

12. Cai Y, Yu X, Hu S, Yu J. A brief review on the mechanisms of miRNA regulation. Genom Proteom Bioinform. (2009) 7:147–54. doi: 10.1016/S1672-0229(08)60044-3

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chen K, Rajewsky N. The evolution of gene regulation by transcription factors and microRNAs. Nat Rev Genet. (2007) 8:93–103. doi: 10.1038/nrg1990

CrossRef Full Text | Google Scholar

14. Hobert O. Gene regulation by transcription factors and microRNAs. Science. (2008) 319:1785–6. doi: 10.1126/science.1151651

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Deng M, Lancto CA, Abrahamsen MS. Cryptosporidium parvum regulation of human epithelial cell gene expression. Int J Parasitol. (2004) 34:73–82. doi: 10.1016/j.ijpara.2003.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Sturm A, Amino R, van de Sand C, Regen T, Retzlaff S, Rennenberg A, et al. Manipulation of host hepatocytes by the malaria parasite for delivery into liver sinusoids. Science. (2006) 313:1287–90. doi: 10.1126/science.1129720

PubMed Abstract | CrossRef Full Text | Google Scholar

17. McDonald V, Korbel DS, Barakat FM, Choudhry N, Petry F. Innate immune responses against Cryptosporidium parvum infection. Parasite Immunol. (2013) 35:55–64. doi: 10.1111/pim.12020

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ambros V. The functions of animal microRNAs. Nature. (2004) 431:350–5. doi: 10.1038/nature02871

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Regis SC, Mendonça LR, Silva Ndos S, Dattoli VC, Alcântara-Neves NM, Barrouin-Melo SM. Seroprevalence and risk factors for canine toxocariasis by detection of specific IgG as a marker of infection in dogs from Salvador, Brazil. Acta Trop. (2011) 120:46–51. doi: 10.1016/j.actatropica.2011.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Mikaeili F, Mathis A, Deplazes P, Mirhendi H, Barazesh A, Ebrahimi S, et al. Differentiation of Toxocara canis and Toxocara cati based on PCR-RFLP analyses of rDNA-ITS and mitochondrial cox1 and nad1 regions. Acta Parasitol. (2017) 62:549–56. doi: 10.1515/ap-2017-0066

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Zheng WB, Zou Y, Elsheikha HM, Liu GH, Hu MH, Wang SL, et al. Serum metabolomic alterations in Beagle dogs experimentally infected with Toxocara canis. Parasit Vectors. (2019) 12:447. doi: 10.1186/s13071-019-3703-5

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Takamiya S, Mita T. Large-scale purification of active liquid-cultured Caenorhabditis elegans using a modified Baermann apparatus. Parasitol Int. (2016) 65:580–3. doi: 10.1016/j.parint.2016.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. (2009) 10:R25. doi: 10.1186/gb-2009-10-3-r25

PubMed Abstract | CrossRef Full Text | Google Scholar

24. 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. doi: 10.1093/nar/gkr688

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Wen M, Shen Y, Shi S, Tang T. miREvo: an integrative microRNA evolutionary analysis platform for next-generation sequencing experiments. BMC Bioinformatics. (2012) 13:140. doi: 10.1186/1471-2105-13-140

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Zhou L, Chen J, Li Z, Li X, Hu X, Huang Y, et al. Integrated profiling of microRNAs and mRNAs: microRNAs located on Xq27.3 associate with clear cell renal cell carcinoma. PLoS ONE. (2010) 5:e15224. doi: 10.1371/journal.pone.0015224

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. (2003) 5:R1. doi: 10.1186/gb-2003-5-1-r1

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Chen CM, Lu YL, Sio CP, Wu GC, Tzou WS, Pai TW. Gene Ontology based housekeeping gene selection for RNA-seq normalization. Methods. (2014) 67:354–63. doi: 10.1016/j.ymeth.2014.01.019

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. (2008) 36:D480–4. doi: 10.1093/nar/gkm882

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. (2005) 21:3787–93. doi: 10.1093/bioinformatics/bti430

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. (2001) 25:402–8. doi: 10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Adnan M, Morton G, Hadi S. Analysis of rpoS and bolA gene expression under various stress-induced environments in planktonic and biofilm phase using 2−ΔΔCT method. Mol Cell Biochem. (2011) 357:275–82. doi: 10.1007/s11010-011-0898-y

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Regeenes R, Silva PN, Chang HH, Arany EJ, Shukalyuk AI, Audet J, et al. Fibroblast growth factor receptor 5 (FGFR5) is a co-receptor for FGFR1 that is up-regulated in beta-cells by cytokine-induced inflammation. J Biol Chem. (2018) 293:17218–28. doi: 10.1074/jbc.RA118.003036

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Wetzel A, Chavakis T, Preissner KT, Sticherling M, Haustein UF, Anderegg U, et al. Human Thy-1 (CD90) on activated endothelial cells is a counterreceptor for the leukocyte integrin Mac-1 (CD11b/CD18). J Immunol. (2004) 172:3850–9. doi: 10.4049/jimmunol.172.6.3850

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Aeschimann F, Neagu A, Rausch M, Großhans H. let-7 coordinates the transition to adulthood through a single primary and four secondary targets. Life Sci Alliance. (2019) 2:e201900335. doi: 10.26508/lsa.201900335

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Schnieder T, Laabs EM, Welz C. Larval development of Toxocara canis in dogs. Vet Parasitol. (2011) 175:193–206. doi: 10.1016/j.vetpar.2010.10.027

CrossRef Full Text | Google Scholar

37. Le Deist F, Hivroz C, Partiseti M, Rieux-Laucat F, Debatin KM, Choquet D, et al. T cell activation deficiencies. Clin Immunol Immunopathol. (1995) 76:163–4. doi: 10.1016/S0090-1229(95)90072-1

CrossRef Full Text | Google Scholar

38. Soudais C, De Villartay JP, Le Deist F, Fischer A, Lisowska-Grospierre B. Genetic analysis of the human CD3-epsilon gene in a T cell receptor/CD3 immunodeficiency. Immunodeficiency. (1993) 4:117–9. doi: 10.1038/ng0193-77

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Oliveira AC, Gomes-Neto JF, Barbosa CD, Granato A, Reis BS, Santos BM, et al. Crucial role for T cell-intrinsic IL-18R-MyD88 signaling in cognate immune response to intracellular parasite infection. Elife. (2017) 6:e30883. doi: 10.7554/eLife.30883

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Wang G, Jin S, Ling X, Li Y, Hu Y, Zhang Y, et al. Proteomic profiling of LPS-induced macrophage-derived exosomes indicates their involvement in acute liver injury. Proteomics. (2019) 19:e1800274. doi: 10.1002/pmic.201800274

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Oliveira KC, Carvalho ML, Bonatto JM, Schechtman D, Verjovski-Almeida S. Human TNF-α induces differential protein phosphorylation in Schistosoma mansoni adult male worms. Parasitol Res. (2016) 115:817–28. doi: 10.1007/s00436-015-4812-5

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Pereira-Fonseca DCM, Oliveira-Rovai FM, Rodas LAC, Beloti CAC, Torrecilha RBP, Ito PKRK, et al. Dog skin parasite load, TLR-2, IL-10 and TNF-α expression and infectiousness. Parasite Immunol. (2017) 39:e12493. doi: 10.1111/pim.12493

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Kim L, Butcher BA, Denkers EY. Toxoplasma gondii interferes with lipopolysaccharide-induced mitogen-activated protein kinase activation by mechanisms distinct from endotoxin tolerance. J Immunol. (2004) 172:3003–10. doi: 10.4049/jimmunol.172.5.3003

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Gava SG, Tavares NC, Falcone FH, Oliveira G, Mourão MM. Profiling Transcriptional regulation and functional roles of schistosoma mansoni c-Jun N-Terminal Kinase. Front Genet. (2019) 10:1036. doi: 10.3389/fgene.2019.01036

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Kraft S, Bieber T. Fc epsilon ri-mediated activation of transcription factors in antigen-presenting cells. Int Arch Allergy Immunol. (2001) 125:9–15. doi: 10.1159/000053791

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Liu K, Yin D, Shu Y, Dai P, Yang Y, Wu H. Transcriptome and metabolome analyses of Coilia nasus in response to Anisakidae parasite infection. Fish Shellfish Immunol. (2019) 87:235–42. doi: 10.1016/j.fsi.2018.12.077

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Sparks MA, Crowley SD, Gurley SB, Mirotsou M, Coffman TM. Classical Renin-Angiotensin system in kidney physiology. Compr Physiol. (2014) 4:1201–28. doi: 10.1002/cphy.c130040

PubMed Abstract | CrossRef Full Text

48. Tomek MB, Maresch D, Windwarder M, Friedrich V, Janesch B, Fuchs K, et al. A general protein O-glycosylation gene cluster encodes the species-specific Glycan of the oral pathogen Tannerella forsythia: O-Glycan biosynthesis and immunological implications. Front Microbiol. (2018) 9:2008. doi: 10.3389/fmicb.2018.02008

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Golias C, Tsoutsi E, Matziridis A, Makridis P, Batistatou A, Charalabopoulos K. Review. Leukocyte and endothelial cell adhesion molecules in inflammation focusing on inflammatory heart disease. In Vivo. (2007) 21:757–69.

PubMed Abstract | Google Scholar

50. Ueland PM, McCann A, Midttun Ø, Ulvik A. Inflammation, vitamin B6 and related pathways. Mol Aspects Med. (2017) 53:10–27. doi: 10.1016/j.mam.2016.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Zhou Y, Qiao H, Yin N, Chen L, Xie Y, Wu J, et al. Immune and cytokine/chemokine responses of PBMCs in rotavirus-infected rhesus infants and their significance in viral pathogenesis. J Med Virol. (2019) 91:1448–69. doi: 10.1002/jmv.25460

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Morla S. Glycosaminoglycans and Glycosaminoglycan mimetics in cancer and inflammation. Int J Mol Sci. (2019) 20:1963. doi: 10.3390/ijms20081963

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Toxocara canis, miRNAs, beagle dogs, liver, RNA-seq

Citation: Zou Y, Zheng W-B, He J-J, Elsheikha HM, Zhu X-Q and Lu Y-X (2020) Toxocara canis Differentially Affects Hepatic MicroRNA Expression in Beagle Dogs at Different Stages of Infection. Front. Vet. Sci. 7:587273. doi: 10.3389/fvets.2020.587273

Received: 25 July 2020; Accepted: 01 October 2020;
Published: 12 November 2020.

Edited by:

Alessia Libera Gazzonis, University of Milan, Italy

Reviewed by:

Quan Liu, Foshan University, China
Guo-Hua Liu, Hunan Agricultural University, China

Copyright © 2020 Zou, Zheng, He, Elsheikha, Zhu and Lu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yi-Xin Lu, luyixin@neau.edu.cn; Xing-Quan Zhu, xingquanzhu1@hotmail.com