Impact Factor 4.235 | CiteScore 6.4
More on impact ›

Original Research ARTICLE

Front. Microbiol., 19 October 2020 | https://doi.org/10.3389/fmicb.2020.570903

Transcriptomic Profiling of Mouse Brain During Acute and Chronic Infections by Toxoplasma gondii Oocysts

Rui-Si Hu1,2, Jun-Jun He1*, Hany M. Elsheikha3, Yang Zou1, Muhammad Ehsan1, Qiao-Ni Ma1, Xing-Quan Zhu1,4 and Wei Cong5*
  • 1State 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
  • 2College of Animal Science and Technology, Jilin Agricultural University, Changchun, 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
  • 5Marine College, Shandong University, Weihai, China

Infection by the protozoan Toxoplasma gondii can have a devastating impact on the structure and function of the brain of the infected individuals, particularly immunocompromised patients. A systems biology view of the brain transcriptome can identify key molecular targets and pathways that mediate the neuropathogenesis of cerebral toxoplasmosis. Here, we performed transcriptomic analysis of the brain of mice infected by T. gondii Pru strain oocysts at 11 and 33 days post-infection (dpi) compared to uninfected (control) mice using RNA sequencing (RNA-seq). T. gondii altered the expression of 936 and 2,081 transcripts at 11 and 33 dpi, respectively, and most of these were upregulated in the infected brains. Gene Ontology (GO) enrichment and pathway analysis showed that immune response, such as interferon-gamma (IFN-γ) responsive genes were strongly affected at 11dpi. Likewise, differentially expressed transcripts (DETs) related to T cell activation, cytokine production and immune cell proliferation were significantly altered at 33 dpi. Host-parasite interactome analysis showed that some DETs were involved in immune signaling, metabolism, biosynthesis-related processes and interspecies interaction. These findings should increase knowledge of the mouse brain transcriptome and the changes in transcriptional regulation and downstream signaling pathways during acute and chronic T. gondii infections.

Introduction

Toxoplasma gondii is a widely prevalent intracellular apicomplexan protozoan parasite which chronically infects approximately one-third of the world’s population (Tenter et al., 2000). Although infection by T. gondii does not normally cause clinical illness in humans, individuals whose immune system is compromised can develop serious health problems. T. gondii has the ability to infect all tissue types in the affected host, however, this parasite has preference for the brain and muscles, where it remains for a long time in these tissues as dormant cysts. Chronic T. gondii infection can lead to neurological disorders, such as bipolar disease and schizophrenia (Elsheikha et al., 2016; Elsheikha and Zhu, 2016; Shiojiri et al., 2019). A correlation between dysregulated immune-inflammation and brain dysfunctions during T. gondii infection has been discussed previously (Parlog et al., 2015; Elsheikha and Zhu, 2016; Wohlfert et al., 2017). Additionally, the dormant T. gondii cysts can reactivate and cause severe encephalitis, leading neurological manifestations and even fatal consequences (Wohlfert et al., 2017). T. gondii infections represent one of the major health challenges that compromise the life of people with HIV and remain a prominent cause of morbidity and mortalities in immunocompromised patients in general (Mayor et al., 2011). One study reported that 65% of HIV patients died within the first year of diagnosis with reactivated latent toxoplasmosis (Mayor et al., 2011).

How T. gondii reaches and colonizes the brain and the molecular mechanisms that underpin these processes continue to challenge scientists. The neuro-pathogenesis of T. gondii infection is driven by a complicated network of molecular processes and pathways (Ngo et al., 2017). Data supporting this view are emerging from various “-omics”-based profiling of the brains of mice infected by T. gondii (Jia et al., 2013; Tanaka et al., 2013; Pittman et al., 2014; Zhou et al., 2015; Hu et al., 2018; Garfoot et al., 2019a; Ma et al., 2019). The magnitude of this challenge has increased the interest in understanding how T. gondii infection affects murine host at the transcriptional level in various tissues other than the brain, including the liver (He et al., 2016a), spleen (He et al., 2016b,c), peripheral lymphocytes (Jia et al., 2013), and the uterus (Zhou et al., 2020). Increased understanding of the brain transcriptomic signature associated with T. gondii infection can provide new testable hypotheses that may ultimately facilitate the development of new treatment interventions to control cerebral toxoplasmosis. However, the effect of infection by T. gondii oocysts on the brain transcriptome remains obscure, although consumption of oocysts with food or water contaminated with cat feces represents a major route of infection of the intermediate hosts (Dubey, 2009; Blader et al., 2015).

In the present study, the transcriptomic response of mouse brain to experimental T. gondii infection by oocysts of type II Pru strain was examined using RNA sequencing (RNA-seq). Our computational analysis revealed significant transcriptomic changes in the infected and uninfected (control) mice at 11 and 33 days post-infection (dpi). We identified 936 differently expressed transcripts (DETs) during acute (11 dpi) infection and 2,081 DETs during chronic (33 dpi) infection. The transcriptional changes correlated with T. gondii load in infected brain tissues. The identified DETs were involved in metabolism, immune-related, biosynthesis-related processes, and interspecies interaction between organisms. Our data not only show how T. gondii oocyst infection alters the transcriptome of the brain of mice at early and late stage of infections, but also provide insight for future research efforts into the development of therapeutic interventions to combat T. gondii infection of the brain.

Materials and Methods

Collection and Preparation of Sporulated Oocysts

One, 10-week-old, specific-pathogen-free (SPF), kitten was fed 100 cysts of Pru strain freshly prepared from brain homogenates obtained from experimentally infected mice. Cat feces were examined microscopically on a daily basis. Oocyst collection was performed as previously described (Zhou et al., 2016), and purification of oocysts was performed using sucrose flotation and cesium chloride gradient as previously described (Staggs et al., 2009). To obtain sporozoites-containing (i.e., infectious) oocysts, the purified oocysts were centrifuged at 360 × g, and the oocyst pellet was suspended in 2% sulfuric acid and was aerated on a shaker for 7 days at room temperature. Sporulated oocysts were washed twice with 0.85% normal saline and suspended in 2% sulfuric acid. Finally, oocysts were counted using a hemocytometer and their number was adjusted to 100 oocyst/ml in phosphate buffered saline (PBS), and stored at 4°C until use.

BALB/c Mice and Parasite Challenge

SPF female BALB/c mice, 7-week-old, were purchased from Lanzhou University Laboratory Animal Centre (Lanzhou, Gansu Province, China). Twelve T. gondii-seronegative mice were randomly divided into four groups (3 mice/group): two mouse groups were infected by T. gondii oocysts for 11 and 33 days, representing the acutely and chronically infected groups, respectively. Two mouse groups were left uninfected for the same length of times as the matched control groups. The mice were handled following the protocol approved by the Animal Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences. T. gondii infection was performed in each mouse by oral gavage with 100 oocysts in 1 ml of PBS. The uninfected (control) mice were sham inoculated with 1 ml of PBS only without any oocysts. Mice were maintained in an Animal Biosafety Level 2 (ABSL-2) containment laboratory, with room temperature (22 ± 0.5°C) under 12/12 h light/dark cycle. All mice were fed commercial diet and had access to purified water ad libitum. Following the infection, mice were monitored for the development of clinical signs, such as ruffled hair, and physical inactivity. At 11 dpi and 33 dpi, mice were anesthetized by intraperitoneal injection with 100 μl xylazine (20 mg/ml) and ketamine (1 mg/ml) in PBS, and brain tissues were harvested, quickly washed in PBS, flash frozen in liquid nitrogen, and preserved at −80°C until use.

Detection of T. gondii in the Brain Tissue

The presence of T. gondii cysts in infected brain tissues has been confirmed in our previous study of microRNA profiling using a microscope (Hu et al., 2018). Additionally, genomic DNA was extracted from the brain tissues using TIANamp Genomic DNA kit (TianGenTM, Beijing, China) and DNA samples were stored at −20°C. PCR amplification of T. gondii B1 gene was used to further confirm the presence of T. gondii in the brain tissue (Cong et al., 2015).

RNA Extraction and Transcriptome Sequencing (RNA-Seq) Analysis

The RNA of each sample were extracted separately using TRIzol Reagent (Invitrogen, CA, United States) according to the manufacturer’s instructions. The quality of the extracted RNA was checked in order to ensure the high quality of RNA library. The purity (OD260/280) of the RNA preparation was examined using a NanoPhotometer® spectrophotometer (IMPLEN, CA, United States). All RNA preparations had an absorbance ratio OD260/280 > 1.8. RNA samples were further treated with PQI DNase (Promega, MI, United States) to remove any genomic DNA contaminates. The integrity and concentration of RNA were examined using RNA Nano 6000 Assay Kit (Agilent Technologies, CA, United States) and Qubit® RNA Assay Kit (Life Technologies, CA, United States), respectively. About one microgram of the total RNA per sample (of 3 pooled replicates) was used for the construction of mRNA libraries using NEBNext® UltraTM RNA Library Prep Kit (Illumina, NEB, United States). Each library was sequenced on an Illumina HiSeqTM 2500 (Illumina, San Diego, CA, United States). Adaptor sequences and low-quality reads were removed using Trimmomatic (v.0.39) (Bolger et al., 2014). Genomic data of Mus musculus (Mouse; GRCm38.p6) was downloaded from Ensembl database (Release-100). Transcript abundance was estimated using free-alignment method by Salmon (v.1.2.0) because of its high sensitivity and accuracy for bias-aware quantification (Patro et al., 2017). Expression profile of each transcript was normalized to TPM (transcripts per million) (Dillies et al., 2013). Differential expression analysis was performed by using count-based technology-DEseq2 package (v.3.6.1) (Love et al., 2014). The DETs were identified according to P-value < 0.05 and | log2 fold-change| ≥ 1.5.

Verification of DETs by Quantitative Real Time qRT-PCR

Total RNA was treated with 1 μl gDNA Eraser, and 2 μl 5 × gDNA Eraser Buffer at 42°C for 2 min in order to remove residual genomic DNA contaminates. Single cDNA was synthesized using a PrimeScriptTM RT Reagent Kit (Takara, Japan). qRT-PCR reaction was carried out using the TB GreenTM Premix Ex TaqTM II (Takara, Japan). The amplification protocol of qRT-PCR was performed as follows: 95°C for 30 s, 40 cycles of 95°C for 5 s, 60°C for 30 s. Melting curve analysis ranging from 65 to 95°C was performed to ensure the specificity of amplification. The measurements of mRNA were normalized to glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and the relative abundances of gene expression were calculated using CT method (2–ΔΔCT method) (Schmittgen and Livak, 2008). Details of the oligonucleotide primers and results of qRT-PCR analysis used to verify the sequencing data are shown in Supplementary Table 2.

Bioinformatic Analysis

ClusterProfiler (v.3.6.2) in RStudio1 was used for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the DETs. The clusterProfiler is an ontology-based tool that not only automates the process of biological-term classification and the enrichment analysis of cluster genes, but also enables the analysis and visualization of the modules (Yu et al., 2012). FDR (false discovery rate) was calculated using Benjamini and Hochberg method and adjusted P-value < 0.05 was set as the cut-off threshold for the significant enrichment. We used TRRUST (v2) (Han et al., 2018) to analyze the transcriptional regulatory networks of mouse. TRRUST is a manually curated database that provides large number of known (activation and repression) and unknown interactions. In this study, only interactions between the transcription factors (TFs) and their targets belonging to the DETs were selected, and external candidates were excluded. For host-parasite interaction analysis, we retrieved mouse and their relevant T. gondii interactions from OTHOHPI database (Cuesta-Astroz et al., 2019). Only host proteins with dysregulation at transcript levels were considered. Interaction network from OTHOHPI was visualized in Cytoscape (v.3.7.0) (Shannon et al., 2003) and functional analysis was performed using Gene Ontology database (Ashburner et al., 2000).

Results

Characterization of the Transcriptomic Dataset

The presence of T. gondii in the mouse brain at 11 and 33 dpi was confirmed using specific semi-nested PCR primers that target the B1 gene (Figure 1). The integrity of each RNA template used in the study was >8.0. The high quality RNAs were used to construct mRNA libraries from infected and uninfected brains of mice. In this study, 94% of the reads showed good quality values >Q20, 90% of clean reads was up to Q30, and in each sample >60 million clean reads were generated (Table 1). The results showed clear separation between infected and control transcriptomes and low variations between the biological replicates for all DETs (adjusted P-value < 0.05; | log2FC| > 1.5) in T. gondii-infected and uninfected (control) brains (Figure 2A). In total, 936 DETs (789 upregulated and 147 downregulated) were detected in the brain of acutely infected mice, and 2,081 DETs (1,908 upregulated and 173 downregulated) were detected in the brain of chronically infected mice (Supplementary Table 1). Two pair-wise comparisons of early and late DETs as shown in the Venn diagram, where the upregulated DETs during early and late infections in the brain tissue show that a total of 618 transcripts were shared between acute and chronic infections (Figure 2B). However, in regard to the downregulated DETs, only 6 transcripts were shared (Figure 2B), suggesting a modest commonality in the downregulated DETs during early and late T. gondii infections. Ten randomly selected transcripts were confirmed by qRT-PCR, and the results were consistent with those obtained by RNA-seq, supporting the validity of the obtained transcriptomic RNA-seq data (Supplementary Table 2).

FIGURE 1
www.frontiersin.org

Figure 1. Detection of T. gondii B1 gene amplified by specific semi-nested PCR in infected and uninfected mouse brain tissues. The amplified band of 98 bp was detected in infected mouse brain (11 and 33 dpi). However, no product was detected in the brain tissues of uninfected mice at the same time points. Lane M: 500 bp DNA molecular ladder; Lane P: positive control; Lane N: negative control; Lanes 1–3: infected brain tissues at 11 dpi; Lanes 4–6: infected brain tissues at 33 dpi; Lanes 7–9: uninfected brain tissues at 11 dpi; Lanes 10–12: uninfected brain tissues at 33 dpi.

TABLE 1
www.frontiersin.org

Table 1. Characteristics of the RNA-seq results obtained in the present study.

FIGURE 2
www.frontiersin.org

Figure 2. Results of RNA sequencing (RNA-seq) and differential expression analysis of transcripts. (A) Heat maps showing the hierarchical clustering of 936 and 2,081 differentially expressed transcripts (DETs) at 11 and 33 dpi, respectively. Each mouse group included three biological replicates. (B) Venn diagram showing the overlap of DETs in the brain of acutely- and chronically-infected mice, including upregulated and downregulated transcripts.

GO Enrichment and KEGG Analysis

We analyzed the GO functional enrichment and significant pathways associated with upregulated and downregulated transcripts in acutely and chronically infected groups (Supplementary Figure 1). GO terms were classified into three categories: Biological Process (BP), Cellular Component (CC) and Molecular Function (MF).

The upregulated transcripts in the brain tissues of acutely infected mice were significantly enriched in infection- or immune-related processes, such as responses to virus, interferon-mediated host response, antigen processing and presentation, and regulation of innate immune responses. For the upregulated transcripts in the brain tissues of chronically infected mice, most of the transcripts enriched in BP terms were involved in activation of immune cells, such as T-cell activation, lymphocyte proliferation, leukocyte proliferation, and mononuclear cell proliferation. Several significantly enriched pathways in both infection stages exhibited similar patterns, which were involved in disease-related physiological pathways or immune responses, for instance, Epstein-Barr virus infection, viral myocarditis, graft-versus-host disease, Type I diabetes mellitus, phagosome, cell adhesion molecules (CAMs), and antigen processing and presentation pathway.

The downregulated transcripts for the two stages of infection were not enriched in any pathway. However, most significantly enriched CC terms of the downregulated transcripts in early infected mice were involved in cell leading edge, presynaptic and synaptic membrane or their components, as well as cell-cell adherens junctions. The most significantly enriched MF terms of the downregulated transcripts in acutely infected mice were involved in phosphatase binding and extracellular matrix structural constituent. BP terms were not significantly enriched in acutely infected mice. However, during chronic infection, the downregulated transcripts were significantly enriched in BP terms involved in the response to metal and calcium ion. However, CC and MF terms were not significantly enriched during chronic infection.

We also performed GO functional enrichment and KEGG analyses of all DETs. According to the adjusted P-value, the top 10 significantly enriched GO terms and pathways with respect to reciprocal comparisons of early and late infections showed that genes were mainly involved in biological processed related to immune responses (Figure 3A). There were 6 significantly enriched GO terms that had cross relationships with 28 transcripts (Figure 3B). These included Acod1, Ccl2, Ccl4, Ccl5, Ccl19, Gbp2, Gbp2b, Gbp3, Gbp5, Gbp6, H2-Ab1, Ifi204, Ifi205, Ifi211, Ifng, Irgm1, Irgm2, IL12rb1, Irf1, Irf8, Mndal, Slc11a1, Socs1, Stat1, Stx11, Tgtp1, Tmem173, and Tlr2. Notably, all of these transcripts were upregulated in the mouse brain during early and late infections.

FIGURE 3
www.frontiersin.org

Figure 3. Functional analysis of the differentially expressed transcripts (DETs). (A) The top 10 Gene Ontology (GO) terms and KEGG pathways of the DETs detected in the brain of acutely- and chronically-infected mice are shown. (B) Circos plot showing the relationship between GO terms and the corresponding genes. A number of genes related to immune response were enriched for biological process such as response to interferon-beta and interferon-gamma, T-cell activation, and positive regulation of cytokine production.

Differentially Expressed Transcription Factors (TFs) and Their Target Genes

According to TRRUST database, we detected 38 unique differentially expressed TFs (32 upregulated and 6 downregulated) in infected brain tissues (Table 2). TRRUST provides information of known interactions in term of regulatory mode for activation or suppression of TFs interacting with their target genes, each with a variety of evidence types and reference sources. We identified 102 unique DETs targeted by the detected TFs that were involved in 185 unique interactions (Figure 4 and Supplementary Table 3). There were 11 upregulated TFs that were co-activated in both acute and chronic infections (Table 2), which functioned as regulators for shared target genes (interaction lines with purple color); of which, Irf1 and Stat1 also can regulate two exclusive target genes in early infection (interaction lines with orange color); and as shown in Figure 4, 35 TFs except Eif2ak2 and Irf9 were involved in more interactions together with the corresponding target genes in late stage of infection (interaction lines with blue color).

TABLE 2
www.frontiersin.org

Table 2. Description of the differentially expressed transcription factors (TFs) detected in the present study.

FIGURE 4
www.frontiersin.org

Figure 4. Co-expression regulatory network of the differentially expressed TFs and their target genes. TFs are represented by orange circles and their gene targets are indicated by green circles. The connected lines are shown in pink (shared early and late infections), orange (early infection only), and blue (late infection only). The regulatory mode of the known interactions (i.e., activation and repression) are shown in Supplementary Table 3.

Mouse Genes Targeted by T. gondii Factors

Based on the identified DETs in infected mouse brain, we detected the host-parasite protein-protein interactions and retrieved 354 computationally predicted interactions (Supplementary Table 4). We found that 159 T. gondii proteins related to the secretome and/or the membrane had relevant roles in mouse-T. gondii interactome networks in which 151 mouse genes were regulated (Supplementary Table 4). As shown in Figure 5, there were 32 genes (upregulated) that were detected in early and chronic stages of infection (i.e., highly targeted by T. gondii), and were involved in GO terms associated with biological processes relevant to the host-parasite interaction such as interspecies interaction between organisms; 34 genes (20 upregulated and 14 downregulated) exclusive to early infection were targeted by T. gondii, and the functional relevance of these interactions as shown by the enriched BPs were related to catabolism, biosynthesis, and metabolism-related processes; and 85 genes (64 upregulated and 21 downregulated) specific to chronic infection were more significantly enriched in specific molecular biosynthesis or localization, and metabolism process, such as small molecule biosynthetic process, macromolecule localization, and metabolic process.

FIGURE 5
www.frontiersin.org

Figure 5. Host-parasite interaction networks and function enrichment analysis of mouse brain genes. DETs were distinguished by acute infection only, shared acute and chronic infections, and chronic infection only. The top 10 biological processes GO terms of the interaction networks are shown.

Discussion

In the present study, we compared the transcriptome of brain tissues of BALB/C mice infected by oocysts of the type II T. gondii Pru strain with that of uninfected mice at 11 and 33 dpi using RNA-seq profiling. The expression of 936 and 2,081 transcripts were altered during acute and chronic infection, respectively. This increased alterations of transcripts during chronic infection is in line with previous results obtained using tachyzoites of type II ME49 strain at acute (10 dpi) and chronic infection (28 dpi) (Pittman et al., 2014). A correlation between expressed host genes and developing cysts has been previously shown (Pittman et al., 2014; Garfoot et al., 2019a). Also, more enriched transcriptome profiles were found in the cysts (i.e., bradyzoite stage) compared to tachyzoite stage (Garfoot et al., 2019b). Long-term infection caused by T. gondii may produce greater changes of transcript abundance, reflecting more host response to the infection (Garfoot et al., 2019b).

In our study, we examined the transcriptomic changes for transcripts at two time points. As the infection progressed within the brain from early to late phases, the upregulated transcripts increased and so as the transcription of immune-related IFN-γ in response to combat the infection (Figure 3A). Also, functional enrichment analysis suggested that chronic infection leads to T cell activation and proliferation of immune cells, such as leukocytes and mononuclear cells (Figure 3A). IFN-γ, CD8+ T cells, NFB1 have been shown to mediate cell mediated immunity against T. gondii infection and contribute to the long-term resistance to T. gondii in the brain (Harris et al., 2010). Interestingly, the downregulation of mice response to calcium was significantly enriched in late infection (Supplementary Figure 1). In T. gondii, calcium-dependent protein kinases exhibited stage-specific expression, for example CDPK5 was more abundant in bradyzoite stage (Pittman et al., 2014). Thus, it is possible that Ca2+ inside parasite cysts play a role in the regulation of host cell environment (Stewart et al., 2017; Zhang et al., 2018).

IFN-γ is a key cytokine in the control and elimination of T. gondii particularly in mice (Suzuki et al., 1988). IFN-γ killing of T. gondii is mediated by upregulation of the immunoregulatory enzyme indoleamine 2,3-dioxygenase (IDO) to deplete cellular tryptophan; inducible nitric oxide synthase (iNOS) to limit cellular arginine; IFN-γ inducible immunity-related GTPase (IRGs) and guanylate binding proteins (GBPs) to induce autophagy-dependent degradation of the parasitophorous vacuole membrane (PVM) (Collazo et al., 2001; Yarovinsky, 2014; Blader et al., 2015). In the present study, we revealed that IFN-γ (Ifng), iNOS (Nos2), six IRGs (Irgm1, Irgm2, Igtp, Iigp1, Tgtp1, and Tgtp2), and 11 GBPs (Gbp2, Gbp2b, Gbp3, Gbp4, Gbp5, Gbp6, Gbp7, Gbp8, Gbp9, Gbp10, and Gbp11) were commonly upregulated during acute and chronic infections (adjusted P-value < 0.05; | log2FC| > 1.5). Consistent with the previous finding reported in T. gondii-infected mouse liver (He et al., 2016a), many IRGs and GBPs-related genes were upregulated. Although increased expression of IDO in the brain as anti-toxoplasma mechanism to control the parasite replication has been reported in previous studies (Nagineni et al., 1996; Zhu et al., 2019), the expression level of IDO did not change in response to T. gondii infection in this study.

Inflammasome-dependent pyroptosis is an important host cell death pathway for controlling T. gondii (Blader et al., 2015; Zamboni and Lima-Junior, 2015). In the present study, we detected nine pyroptosis-related DETs in the infected brain, including three (Casp4, Gsdmd, and Pycard) that were upregulated in both early and late infections; five (Casp1, Naip2, Naip5, Naip6, and Nlrp1b) were upregulated and one (Gsdme) was downregulated in chronic infection. Among them, Casp1 (Caspase-1, an inflammatory caspase molecule), Pycard (ASC, an adapter molecule), and Nlrp1b (a tandem Nrlp1 paralog, a sensor molecule) are components of the canonical inflammasome (Zamboni and Lima-Junior, 2015). Consistent with this, inflammasome in response to T. gondii requires Nrlp1b-mediated Casp1 activation during late infection (Ewald et al., 2014). Our results also showed that two pro-inflammatory cytokines, IL-1β and IL-18 related to Casp1 activation (Rathinam et al., 2012), were upregulated during chronic infection (Supplementary Table 1), suggesting that these cytokines may help mice to mount an inflammatory response that prevents reactivation of latent T. gondii infection. In addition, the specificity of Nrlp1b activation depends on the polymorphism of the sensor molecule (Zamboni and Lima-Junior, 2015), whether this is a response specific to T. gondii or a generic response to any infection remains to be determined.

The present study detected 38 differentially expressed transcription factors (TFs) and 102 differentially expressed gene targets in T. gondii-infected brain tissues. There were 11 TFs commonly activated during early and late infections (Figure 4 and Supplementary Table 3). For example, the Irf8 within the infected dendritic cells (DCs), especially the CD8+ DCs, has a critical role in detecting Toll-like receptor-mediated T. gondii profilin and in the subsequent induction of IL-12 by Myd88-mediated signaling pathway (Yarovinsky, 2014). The Irf1 and Stat1 contribute to the control of T. gondii by regulating the expression of factors essential for the host to resist infection, such as Socs1 (Zimmermann et al., 2006), Ciita (Rosowski et al., 2014), and NOS2 (Gavrilescu et al., 2004). We also detected 35 TFs exclusive to chronic infection that were involved in more interactions within the mouse brain (Table 2 and Figure 4 blue line color).

Previous data suggested correlated gene expressions between the host and the parasite, which can be attributed to the molecular interactions that mediate infection (Reid and Berriman, 2013). T. gondii can utilize its secretory effector proteins to manipulate the transcriptome and various signaling pathways in the host cells to its own benefits (Blader et al., 2015; Hakimi et al., 2017). To study host-pathogen interactomes, homology-based prediction is the one of most common approaches that has been applied to study a variety of intra-species interactions (Lee et al., 2008; Reid and Berriman, 2013; Cuesta-Astroz et al., 2019). More recently, bioinformatic analysis (Cuesta-Astroz et al., 2019) employed this strategy, and predicted effector molecules of T. gondii and their interactions with the host. There are certain ways by which T. gondii can evade host defense mechanisms, such as directly inactivating the IFN-γ-inducible GTPase by using effector molecules and avoiding IFN-dependent killing by targeting host signaling and gene transcription (Blader et al., 2015; Hakimi et al., 2017). By comparing the effect of T. gondii-related factors on mouse brain over the course of infection and distinguishing functional differences of the related DETs, we found that dysregulated transcripts expressed exclusively in early and late infection stages were significantly enriched in biochemical processes-related regulations (Figure 5). These results show the ability of T. gondii factors to alter the brain homeostasis, leading to neuropathology and disruption of neurotransmitter (David et al., 2016). Our analysis also revealed shared DETs between early and late infections which were significantly enriched in immune-related terms (Figure 5). Collectively, these results suggest that T. gondii has significant capabilities to manipulate host cell and can dysregulate many transcripts and processes in the mouse brain.

Conclusion

RNA sequencing was used to characterize the transcriptional changes of the mouse brain during infection by T. gondii oocysts at 11 and 33 dpi. The differentially expressed transcripts (DETs) during late infection were nearly double the number of DETs during early infection. Most of the upregulated transcripts are shared between early and late infection stages, and the greatest difference detected in the downregulated transcripts. Acute infection-induced DETs were mainly associated with IFN-γ inducible immune response and host defense, and IFN-γ was also necessary to maintain chronic infection. Notably, chronic infection resulted in the activation of T-cells and the proliferation of other immune cells such as leukocytes and lymphocytes. T. gondii infection altered the expression of TFs and their target genes, particularly at late stage of infection. Significant interactions between T. gondii and the mouse brain DETs were also detected. These data improve our understanding of the transcriptomic changes that occur during T. gondii oocyst infection in mice and provide new insight into the neuropathogenesis of T. gondii infection.

Data Availability Statement

The datasets supporting the findings of this article are included within the article and its Supplementary Material. The RNA-seq raw data described in the present study has been submitted to the NCBI Short Read Archive database (https://www.ncbi.nlm.nih.gov/sra) under the bio-project number PRJNA483261.

Ethics Statement

The animal study was reviewed and approved by the experimental protocol of this study was reviewed and approved by Animal Research Ethics Committee of Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences. The animals were handled in compliance with the animal ethics requirements of the People’s Republic of China. All efforts were made to minimize animal suffering during the experiment.

Author Contributions

WC, J-JH, and X-QZ conceived and designed the study. R-SH, YZ, ME, Q-NM, and WC performed the experiments. R-SH, J-JH, HE, and WC contributed reagents, materials and analysis tools. R-SH and WC analyzed the data and wrote the manuscript. HE, J-JH, WC, and X-QZ critically revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

Project support was kindly provided by the National Natural Science Foundation of China (Grant No. 31902291), 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 would like to thank Beijing Novogene Biotechnology Company Limited for technical assistance.

Supplementary Material

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

Supplementary Figure 1 | (Top panel) the top 10 GO terms in the categories: biological process (BP), cellular component (CC) and molecular function (MF), and KEGG analysis of the enriched upregulated transcripts in the mouse brain during acute and chronic T. gondii infections. Bottom panel the top 10 GO terms of the enriched down-regulated transcripts. BP terms were not significantly enriched in acute infection, and CC and MF terms were not significantly enriched in chronic infection.

Supplementary Table 1 | Differentially expressed transcripts identified in the present study during acute and chronic T. gondii infections.

Supplementary Table 2 | Oligonucleotide primers and results of qRT-PCR analysis used to verify the sequencing data.

Supplementary Table 3 | Details of the correlation between TFs and their target genes in the co-expression regulatory network.

Supplementary Table 4 | Details of host-parasite interaction analysis, including T. gondii-related factors and their functions.

Footnotes

  1. ^ https://rstudio.com/

References

Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat. Genet. 25, 25–29. doi: 10.1038/75556

PubMed Abstract | CrossRef Full Text | Google Scholar

Blader, I. J., Coleman, B. I., Chen, C. T., and Gubbels, M. J. (2015). Lytic cycle of Toxoplasma gondii: 15 years later. Annu. Rev. Microbiol. 69, 463–485. doi: 10.1146/annurev-micro-091014-104100

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Collazo, C. M., Yap, G. S., Sempowski, G. D., Lusby, K. C., Tessarollo, L., Vande Woude, G. F., et al. (2001). Inactivation of LRG-47 and IRG-47 reveals a family of interferon gamma-inducible genes with essential, pathogen-specific roles in resistance to infection. J. Exp. Med. 194, 181–188. doi: 10.1084/jem.194.2.181

PubMed Abstract | CrossRef Full Text | Google Scholar

Cong, W., Liu, G. H., Meng, Q. F., Dong, W., Qin, S. Y., Zhang, F. K., et al. (2015). Toxoplasma gondii infection in cancer patients: prevalence, risk factors, genotypes and association with clinical diagnosis. Cancer Lett. 359, 307–313. doi: 10.1016/j.canlet.2015.01.036

CrossRef Full Text | Google Scholar

Cuesta-Astroz, Y., Santos, A., Oliveira, G., and Jensen, L. J. (2019). Analysis of predicted host-parasite interactomes reveals commonalities and specificities related to parasitic lifestyle and tissues tropism. Front. Immunol. 10:212. doi: 10.3389/fimmu.2019.00212

PubMed Abstract | CrossRef Full Text | Google Scholar

David, C. N., Frias, E. S., Szu, J. I., Vieira, P. A., Hubbard, J. A., Lovelace, J., et al. (2016). GLT-1-dependent disruption of CNS glutamate homeostasis and neuronal function by the protozoan parasite Toxoplasma gondii. PLoS Pathog. 12:e1005643. doi: 10.1371/journal.ppat.1005643

PubMed Abstract | CrossRef Full Text | Google Scholar

Dillies, M. A., Rau, A., Aubert, J., Hennequet-Antier, C., Jeanmougin, M., Servant, N., et al. (2013). A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis. Brief. Bioinform. 14, 671–683. doi: 10.1093/bib/bbs046

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubey, J. P. (2009). History of the discovery of the life cycle of Toxoplasma gondii. Int. J. Parasitol. 39, 877–882. doi: 10.1016/j.ijpara.2009.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Elsheikha, H. M., Büsselberg, D., and Zhu, X. Q. (2016). The known and missing links between Toxoplasma gondii and schizophrenia. Metab. Brain Dis. 31, 749–759. doi: 10.1007/s11011-016-9822-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Elsheikha, H. M., and Zhu, X. Q. (2016). Toxoplasma gondii infection and schizophrenia: an inter-kingdom communication perspective. Curr. Opin. Infect. Dis. 29, 311–318. doi: 10.1097/QCO.0000000000000265

CrossRef Full Text | Google Scholar

Ewald, S. E., Chavarria-Smith, J., and Boothroyd, J. C. (2014). NLRP1 is an inflammasome sensor for Toxoplasma gondii. Infect. Immun. 82, 460–468. doi: 10.1128/IAI.01170-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Garfoot, A. L., Cervantes, P. W., and Knoll, L. J. (2019a). Transcriptional analysis shows a robust host response to Toxoplasma gondii during early and late chronic infection in both male and female mice. Infect. Immun. 87:e00024-19. doi: 10.1128/IAI.00024-19

CrossRef Full Text | Google Scholar

Garfoot, A. L., Wilson, G. M., Coon, J. J., and Knoll, L. J. (2019b). Proteomic and transcriptomic analyses of early and late-chronic Toxoplasma gondii infection shows novel and stage specific transcripts. BMC Genomics 20:859. doi: 10.1186/s12864-019-6213-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Gavrilescu, L. C., Butcher, B. A., Del Rio, L., Taylor, G. A., and Denkers, E. Y. (2004). STAT1 is essential for antimicrobial effector function but dispensable for gamma interferon production during Toxoplasma gondii infection. Infect. Immun. 72, 1257–1264. doi: 10.1128/iai.72.3.1257-1264.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Hakimi, M. A., Olias, P., and Sibley, L. D. (2017). Toxoplasma effectors targeting host signaling and transcription. Clin. Microbiol. Rev. 30, 615–645. doi: 10.1128/CMR.00005-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, H., Cho, J. W., Lee, S., Yun, A., Kim, H., Bae, D., et al. (2018). TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucleic Acids Res. 46, D380–D386. doi: 10.1093/nar/gkx1013

PubMed Abstract | CrossRef Full Text | Google Scholar

Harris, T. H., Wilson, E. H., Tait, E. D., Buckley, M., Shapira, S., Caamano, J., et al. (2010). NF-κB1 contributes to T cell-mediated control of Toxoplasma gondii in the CNS. J. Neuroimmunol. 222, 19–28. doi: 10.1016/j.jneuroim.2009.12.009

PubMed Abstract | CrossRef Full Text | Google Scholar

He, J. J., Ma, J., Elsheikha, H. M., Song, H. Q., Huang, S. Y., and Zhu, X. Q. (2016a). Transcriptomic analysis of mouse liver reveals a potential hepato-enteric pathogenic mechanism in acute Toxoplasma gondii infection. Parasit. Vectors 9:427. doi: 10.1186/s13071-016-1716-x

CrossRef Full Text | Google Scholar

He, J. J., Ma, J., Li, F. C., Song, H. Q., Xu, M. J., and Zhu, X. Q. (2016b). Transcriptional changes of mouse splenocyte organelle components following acute infection with Toxoplasma gondii. Exp. Parasitol. 167, 7–16. doi: 10.1016/j.exppara.2016.04.019

PubMed Abstract | CrossRef Full Text | Google Scholar

He, J. J., Ma, J., Song, H. Q., Zhou, D. H., Wang, J. L., Huang, S. Y., et al. (2016c). Transcriptomic analysis of global changes in cytokine expression in mouse spleens following acute Toxoplasma gondii infection. Parasitol. Res. 115, 703–712. doi: 10.1007/s00436-015-4792-5

CrossRef Full Text | Google Scholar

Hu, R. S., He, J. J., Elsheikha, H. M., Zhang, F. K., Zou, Y., Zhao, G. H., et al. (2018). Differential brain microRNA expression profiles after acute and chronic infection of mice with Toxoplasma gondii oocysts. Front. Microbiol. 9:2316. doi: 10.3389/fmicb.2018.02316

CrossRef Full Text | Google Scholar

Jia, B., Lu, H., Liu, Q., Yin, J., Jiang, N., and Chen, Q. (2013). Genome-wide comparative analysis revealed significant transcriptome changes in mice after Toxoplasma gondii infection. Parasit. Vectors 6:161. doi: 10.1186/1756-3305-6-161

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S. A., Chan, C. H., Tsai, C. H., Lai, J. M., Wang, F. S., Kao, C. Y., et al. (2008). Ortholog-based protein-protein interaction prediction and its application to inter-species interactions. BMC Bioinformatics 9(Suppl. 12):S11. doi: 10.1186/1471-2105-9-S12-S11

PubMed Abstract | CrossRef Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8

CrossRef Full Text | Google Scholar

Ma, J., He, J. J., Hou, J. L., Zhou, C. X., Zhang, F. K., Elsheikha, H. M., et al. (2019). Metabolomic signature of mouse cerebral cortex following Toxoplasma gondii infection. Parasit. Vectors 12:373. doi: 10.1186/s13071-019-3623-4

CrossRef Full Text | Google Scholar

Mayor, A. M., Fernandez Santos, D. M., Dworkin, M. S., Rios-Olivares, E., and Hunter-Mellado, R. F. (2011). Toxoplasmic encephalitis in an AIDS cohort at Puerto Rico before and after highly active antiretroviral therapy (HAART). Am. J. Trop. Med. Hyg. 84, 838–841. doi: 10.4269/ajtmh.2011.10-0718

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagineni, C. N., Pardhasaradhi, K., Martins, M. C., Detrick, B., and Hooks, J. J. (1996). Mechanisms of interferon-induced inhibition of Toxoplasma gondii replication in human retinal pigment epithelial cells. Infect. Immun. 64, 4188–4196. doi: 10.1128/iai.64.10.4188-4196.1996

PubMed Abstract | CrossRef Full Text | Google Scholar

Ngo, H. M., Zhou, Y., Lorenzi, H., Wang, K., Kim, T. K., Zhou, Y., et al. (2017). Toxoplasma modulates signature pathways of human epilepsy, neurodegeneration & cancer. Sci. Rep. 7:11496. doi: 10.1038/s41598-017-10675-6

CrossRef Full Text | Google Scholar

Parlog, A., Schluter, D., and Dunay, I. R. (2015). Toxoplasma gondii-induced neuronal alterations. Parasite Immunol. 37, 159–170. doi: 10.1111/pim.12157

CrossRef Full Text | Google Scholar

Patro, R., Duggal, G., Love, M. I., Irizarry, R. A., and Kingsford, C. (2017). Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419. doi: 10.1038/nmeth.4197

CrossRef Full Text | Google Scholar

Pittman, K. J., Aliota, M. T., and Knoll, L. J. (2014). Dual transcriptional profiling of mice and Toxoplasma gondii during acute and chronic infection. BMC Genomics 15:806. doi: 10.1186/1471-2164-15-806

PubMed Abstract | CrossRef Full Text | Google Scholar

Rathinam, V. A., Vanaja, S. K., and Fitzgerald, K. A. (2012). Regulation of inflammasome signaling. Nat. Immunol. 13, 333–342. doi: 10.1038/ni.2237

CrossRef Full Text | Google Scholar

Reid, A. J., and Berriman, M. (2013). Genes involved in host-parasite interactions can be revealed by their correlated expression. Nucleic Acids Res. 41, 1508–1518. doi: 10.1093/nar/gks1340

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosowski, E. E., Nguyen, Q. P., Camejo, A., Spooner, E., and Saeij, J. P. (2014). Toxoplasma gondii inhibits gamma interferon (IFN-gamma)- and IFN-beta-induced host cell STAT1 transcriptional activity by increasing the association of STAT1 with DNA. Infect. Immun. 82, 706–719. doi: 10.1128/IAI.01291-13

CrossRef Full Text | Google Scholar

Schmittgen, T. D., and Livak, K. J. (2008). Analyzing real-time PCR data by the comparative C(T) method. Nat. Protoc. 3, 1101–1108. doi: 10.1038/nprot.2008.73

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Shiojiri, D., Kinai, E., Teruya, K., Kikuchi, Y., and Oka, S. (2019). Combination of clindamycin and azithromycin as alternative treatment for Toxoplasma gondii encephalitis. Emerg. Infect. Dis. 25, 841–843. doi: 10.3201/eid2504.181689

PubMed Abstract | CrossRef Full Text | Google Scholar

Staggs, S. E., See, M. J., Dubey, J. P., and Villegas, E. N. (2009). Obtaining highly purified Toxoplasma gondii oocysts by a discontinuous cesium chloride gradient. J. Vis. Exp. 33:1420. doi: 10.3791/1420

PubMed Abstract | CrossRef Full Text | Google Scholar

Stewart, R. J., Whitehead, L., Nijagal, B., Sleebs, B. E., Lessene, G., McConville, M. J., et al. (2017). Analysis of Ca(2)(+) mediated signaling regulating Toxoplasma infectivity reveals complex relationships between key molecules. Cell. Microbiol. 19:e12685. doi: 10.1111/cmi.12685

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, Y., Orellana, M. A., Schreiber, R. D., and Remington, J. S. (1988). Interferon-gamma: the major mediator of resistance against Toxoplasma gondii. Science 240, 516–518. doi: 10.1126/science.3128869

CrossRef Full Text | Google Scholar

Tanaka, S., Nishimura, M., Ihara, F., Yamagishi, J., Suzuki, Y., and Nishikawa, Y. (2013). Transcriptome analysis of mouse brain infected with Toxoplasma gondii. Infect. Immun. 81, 3609–3619. doi: 10.1128/IAI.00439-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Tenter, A. M., Heckeroth, A. R., and Weiss, L. M. (2000). Toxoplasma gondii: from animals to humans. Int. J. Parasitol. 30, 1217–1258. doi: 10.1016/s0020-7519(00)00124-7

CrossRef Full Text | Google Scholar

Wohlfert, E. A., Blader, I. J., and Wilson, E. H. (2017). Brains and brawn: Toxoplasma infections of the central nervous system and skeletal muscle. Trends Parasitol. 33, 519–531. doi: 10.1016/j.pt.2017.04.001

CrossRef Full Text | Google Scholar

Yarovinsky, F. (2014). Innate immunity to Toxoplasma gondii infection. Nat. Rev. Immunol. 14, 109–121. doi: 10.1038/nri3598

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–287. doi: 10.1089/omi.2011.0118

CrossRef Full Text | Google Scholar

Zamboni, D. S., and Lima-Junior, D. S. (2015). Inflammasomes in host response to protozoan parasites. Immunol. Rev. 265, 156–171. doi: 10.1111/imr.12291

CrossRef Full Text | Google Scholar

Zhang, N. Z., Gao, Q., Wang, M., Hou, J. L., Zhang, F. K., Hu, L. Y., et al. (2018). Protective efficacy against acute and chronic Toxoplasma gondii infection induced by immunization with the DNA vaccine TgDOC2C. Front. Microbiol. 9:2965. doi: 10.3389/fmicb.2018.02965

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, C. X., Zhou, D. H., Elsheikha, H. M., Liu, G. X., Suo, X., and Zhu, X. Q. (2015). Global metabolomic profiling of mice brains following experimental infection with the cyst-forming Toxoplasma gondii. PLoS One 10:e0139635. doi: 10.1371/journal.pone.0139635

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, C. X., Zhu, X. Q., Elsheikha, H. M., He, S., Li, Q., Zhou, D. H., et al. (2016). Global iTRAQ-based proteomic profiling of Toxoplasma gondii oocysts during sporulation. J. Proteomics 148, 12–19. doi: 10.1016/j.jprot.2016.07.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Zhang, X. X., Mahmmod, Y. S., Hernandez, J. A., Li, G. F., Huang, W. Y., et al. (2020). A transcriptome analysis: various reasons of adverse pregnancy outcomes caused by acute Toxoplasma gondii infection. Front. Physiol. 11:115. doi: 10.3389/fphys.2020.00115

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, W., Li, J., Pappoe, F., Shen, J., and Yu, L. (2019). Strategies developed by Toxoplasma gondii to survive in the host. Front. Microbiol. 10:899. doi: 10.3389/fmicb.2019.00899

PubMed Abstract | CrossRef Full Text | Google Scholar

Zimmermann, S., Murray, P. J., Heeg, K., and Dalpke, A. H. (2006). Induction of suppressor of cytokine signaling-1 by Toxoplasma gondii contributes to immune evasion in macrophages by blocking IFN-gamma signaling. J. Immunol. 176, 1840–1847. doi: 10.4049/jimmunol.176.3.1840

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Toxoplasma gondii, oocysts, RNA-seq, transcriptome, cerebral toxoplasmosis, host-parasite interaction

Citation: Hu R-S, He J-J, Elsheikha HM, Zou Y, Ehsan M, Ma Q-N, Zhu X-Q and Cong W (2020) Transcriptomic Profiling of Mouse Brain During Acute and Chronic Infections by Toxoplasma gondii Oocysts. Front. Microbiol. 11:570903. doi: 10.3389/fmicb.2020.570903

Received: 09 June 2020; Accepted: 21 September 2020;
Published: 19 October 2020.

Edited by:

Lihua Xiao, South China Agricultural University, China

Reviewed by:

Si-Yang Huang, Yangzhou University, China
Fangli Lu, Sun Yat-sen University, China
Rui Fang, Huazhong Agricultural University, China

Copyright © 2020 Hu, He, Elsheikha, Zou, Ehsan, Ma, Zhu and Cong. 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: Jun-Jun He, hejunjun617@163.com; Wei Cong, messicw@163.com