Original Research ARTICLE
Modulation of Host miRNAs Transcriptome in Lung and Spleen of Peste des Petits Ruminants Virus Infected Sheep and Goats
- 1Computational Biology and Genomics Facility Lab, Division of Veterinary Biotechnology, Indian Council of Agricultural Research-Indian Veterinary Research Institute, Bareilly, India
- 2Division of Biological Products, Indian Council of Agricultural Research-Indian Veterinary Research Institute, Bareilly, India
- 3Division of Biological Standardization, Indian Council of Agricultural Research-Indian Veterinary Research Institute, Bareilly, India
- 4Division of Virology, Indian Council of Agricultural Research-Indian Veterinary Research Institute, Nainital, India
Peste des petits ruminants (PPR) is one of the highly contagious viral disease, characterized by fever, sore mouth, conjunctivitis, gastroenteritis, and pneumonia, primarily affecting sheep and goats. Reports suggested variable host response in goats and sheep and this host response vis-a-vis the expression of microRNAs (miRNAs) has not been investigated. Here, miRNAs were sequenced and proteomics data were generated to identify the role of differentially expressed miRNA (DEmiRNA) in PPR virus (PPRV) infected lung and spleen tissues of sheep and goats. In lungs, 67 and 37 DEmiRNAs have been identified in goats and sheep, respectively. Similarly, in spleen, 50 and 56 DEmiRNAs were identified in goats and sheep, respectively. A total of 20 and 11 miRNAs were found to be common differentially expressed in both the species in PPRV infected spleen and lung, respectively. Six DEmiRNAs—miR-21-3p, miR-1246, miR-27a-5p, miR-760-3p, miR-320a, and miR-363 were selected based on their role in viral infections, apoptosis, and fold change. The target prediction analysis of these six selected DEmiRNAs from the proteome data generated, revealed involvement of more number of genes in lung and spleen of goats than in sheep. On gene ontology analysis of host target genes these DEmiRNAs were found to regulate several immune response signaling pathways. It was observed that the pathways viz. T cell receptor signaling, Rap1 signaling, Toll-like receptor signaling, and B cell receptor signaling governed by DEmiRNAs were more perturbed in goats than in sheep. The data suggests that PPRV-induced miR-21-3p, miR-320a, and miR-363 might act cooperatively to enhance viral pathogenesis in the lung and spleen of sheep by downregulating several immune response genes. The study gives an important insight into the molecular pathogenesis of PPR by identifying that the PPRV—Izatnagar/94 isolate elicits a strong host response in goats than in sheep.
MicroRNAs (miRNAs) are an important part of the host’s regulatory system, involved in post-transcriptional regulation of gene expression in animals, plants, and some DNA viruses (Sevignani et al., 2006). They regulate gene expression by recognizing partial complementary sites, typically within the 3′ untranslated region (3′UTR) of specific mRNAs. Evidence also supports that miRNAs may regulate gene expression by binding to 5′UTR or coding region (Tay et al., 2008; Roberts et al., 2011a). miRNAs are shown to be involved in different biological processes, including reproduction, development, pathogenesis, apoptosis, and signal transduction (Ambros, 2004; Bartel, 2004; Sahu et al., 2015). It has also been suggested that miRNAs may be the effectors in controlling immune regulation, including cellular differentiation and immune response (Rodriguez et al., 2007; Thai et al., 2007; Johnnidis et al., 2008). They are considered as the centered factors in the interaction network between viruses and host. Studies demonstrated that numerous cellular miRNAs (host miRNAs) play a regulatory role in the host–virus interaction network (Scaria et al., 2006; Grassmann and Jeang, 2008).
Cellular miRNAs can greatly influence viral replication and pathogenesis by promoting or inhibiting virus replication (Guo et al., 2013; Li et al., 2014; Mizuguchi et al., 2015). Viral infection also exerts a profound impact on cellular miRNA expression profile, by altering the expression of cellular miRNAs, thereby regulating host or viral RNA targets (Skalsky and Cullen, 2010). It has been observed that miR-142 suppresses replication of Eastern Equine Encephalitis virus (Trobaugh et al., 2014) and miR-122 enhances replication of Hepatitis C virus (Chang et al., 2008). HIV-1, boosts the expression of several host miRNAs, including miR-122, miR-370, miR-373, and miR-297 and suppresses the expression of the miR-17-92 cluster via an unknown mechanism (Roberts et al., 2011b). The emergence of deep sequencing technology has overcome the limitations of miRNA research. Several studies have explored this technology to evaluate global changes in miRNAs expression in response to virus infection (Wang et al., 2009; Cui et al., 2010).
Peste des petits ruminants (PPR) is an acute, highly contagious viral disease of sheep and goats characterized by fever, sore mouth, conjunctivitis, gastroenteritis, and pneumonia. Goats have been found to be more susceptible with severe form of clinical disease than sheep (Lefevre and Diallo, 1990; Nanda et al., 1996; Dhar et al., 2002; Singh et al., 2004a; Delil et al., 2012; Truong et al., 2014). It has also been observed that the rate of recovery is lower in goats than in sheep (Singh et al., 2004a). However, severe outbreaks of PPR in regions having large sheep populations have also been reported (Singh et al., 2004a; Raghavendra et al., 2008; Maganga et al., 2013). Recently, host–virus interaction studies in PPR have uncovered transcription factors modulating immune response to Sungri/96 live attenuated vaccine strain and predicted an immune signaling pathway that induces immune response (Manjunath et al., 2015, 2017). However, the host miRNAome in PPR has not been explored till date. In the present study, miRNAs were sequenced and proteomics data were generated to examine the effect of PPR virus (PPRV) on host miRNAs expression vis-a-vis protein expression in lung and spleen tissues of sheep and goats infected with PPR.
Materials and Methods
Ethics Statement and Animal Experiment
The vaccine potency testing experiment was carried out at ICAR-Indian Veterinary Research Institute Mukteshwar Campus as per the guidelines of Indian Pharmacopoeia-2014. The study was done after obtaining permission from Indian Veterinary Research Institute Animal Ethics Committee (IVRI-IAEC) under the Committee for the Purpose of Control and Supervision of Experiments on Animals (CPCSEA), India. The protocols were approved vide letter no 387/CPCSEA. Animals (ca. 1 year of age) for the experiment were initially tested to be negative for the presence of PPRV antibody by competitive ELISA (Singh et al., 2004b) and serum neutralization test (SNT; Dhinakar Raj et al., 2000). The animals were also found negative for PPRV antigen in nasal, ocular, buccal, and rectal swabs by sandwich ELISA (Singh et al., 2004c). A highly virulent PPRV (Izatnagar/94 - lineage IV) isolate maintained at PPR Laboratory, Division of Virology, Indian Veterinary Research Institute, Mukteshwar was used as a challenge virus (Sreenivasa et al., 2002). The accession number of this isolate is (KR140086.1; Sahu et al., 2017). Splenic suspension (10%) of virulent virus was inoculated subcutaneously (4 ml suspension, 2 ml each at two different sites). The unvaccinated infected group animals were monitored diurnally for, rectal temperature, any secretion from natural orifices, and feeding habit throughout the experimental period. The unvaccinated animals infected with the PPRV, developed symptoms characteristics of PPRV. The infected animals in which the temperature dropped subnormal were euthanized at 10 days post-infection. As PPRV is epitheliotropic and lymphotropic virus, the tissue samples—lung (epithelial) and spleen (lymphoid) were collected from PPRV infected sheep and goats (n = 2 for each of the species). The counterpart healthy tissues (control) were collected from nearby slaughter house from apparently healthy animals that were screened for the absence of PPRV antigen by sandwich ELISA and antibodies by competitive ELISA and SNT.
Confirmation of PPRV Infection
PPRV infection was confirmed in lung and spleen tissues by, RT-PCR, qRT-PCR, and sandwich ELISA.
Small RNA Library Construction
Total RNA from each of the collected samples (lung and spleen) was isolated using the RNeasy Mini kit (Qiagen GmbH, Germany) according to the manufacturer’s protocol. The integrity and quantity of isolated RNA were assessed on a Bioanalyzer (Agilent Technologies, Inc). The RNA integrity number (RIN) value of all the samples was found greater than 8, which is considered suitable for further processing (Kukurba and Montgomery, 2015). The library was prepared using NEBNext Multiplex Small RNA Library Prep Kit (New England Biolabs Inc.) following the manufacturer’s protocol. Hundred nanograms of total RNA from each sample was used for small RNA library preparation. The quality of the libraries was assessed on Bioanalyzer. Libraries were quantified using a Qubit 2.0 Fluorometer (Life Technologies) and by quantitative real-time PCR (qPCR; Robin et al., 2016). The high-throughput sequencing was performed on Illumina – NextSeq500 (75 bp single-end) (manufacturer’s protocol).
miRNAs Prediction and Analysis
The cattle genome sequence was obtained from ftp://ftp.ensembl.org/pub/release-89/fasta/bos_taurus/dna/ The genome was indexed using Bowtie short read aligner program (Langmead et al., 2009). miRNAs are conserved across species (Altuvia et al., 2005). Since there is no complete sheep and goats miRNAs dataset available in miRBase, known mature miRNAs and precursor sequences for cattle were obtained using a Perl script from miRBase database. Data were further processed using miRDeep2 software (Friedlander et al., 2008, 2012). After filtering (read length ≥18 nt) the set of collapsed, non-redundant, clean reads were mapped to the indexed cattle genome using a mapper module. To identify known miRNAs, clean reads were aligned against miRNA precursor sequences reported in the miRBase database using quantifier.pl module. Read counts for each miRNA identified using miRDeep2 were further used for the downstream analysis. The relative expression levels of miRNAs were normalized as TMM (trimmed mean of M-values) using edgeR R/Bioconductor package (Robinson et al., 2010) to identify differentially expressed miRNAs (DEmiRNAs) in lung and spleen.
Proteomics Data Generation and Data Analysis
Protein Extraction and Analysis
Approximately, 1 g of tissue—spleen and lung, from PPRV infected/apparently healthy goats and sheep was taken in 10 ml lysis buffer (50 mM Tris buffer + PIC + PMSF), homogenized on ice, centrifuged at 14,000 rpm and the supernatant was collected into a separate tube (Tris buffer extract). The cell pellet was further added with urea lysis buffer, centrifuged at 14,000 rpm and the resultant supernatant was collected. Supernatants from Tris and urea extractions were run on an SDS-PAGE for quality check (QC) and for further downstream processing. Protein concentration was determined by using Bradford assay and 100 μg of the samples was taken for digestion. Protein samples were treated with 100 mM dithiothreitol at 95°C for 1 h, added with 250 mM iminodiacetic acid and kept at room temperature for 45 min in dark. Samples were then digested with trypsin and incubated overnight at 37°C. Further, 1% of formic acid was added and incubated at 37°C for 45 min. The resulting samples were vacuum dried and dissolved in 10 μl of 0.1% formic acid and centrifuged at 10,000 g. The supernatant was injected on C18 Nano-LC column for separation of peptides followed by analysis on the Waters Synapt G2 Q-TOF instrument for MS and MSMS. The raw data was processed by MassLynx 4.1 WATERS. The individual peptides MSMS spectra were matched to the database sequence for protein identification on PLGS software, WATERS. Based on the m/z values and their probability to match with a specific peptide present in proteins cleaved at arginine (R) or lysine (K) the protein identification was carried with thresholds, minimum number of peptides to be found for a protein—2; minimum number of fragments (MSMS) ions in a peptide—3; minimum number of fragments (MSMS) ions in a protein—7; peptide mass tolerance—30 ppm; and fragment ion mass tolerance—70 ppm. The identified proteins in the three runs of each sample were compared with each other as control (healthy) and infected samples. Expression Analysis package of the PLGS software was then used for quantification. The ion counts matching with the peptides of a specific protein corresponding between the two samples in the three runs, were averaged and the ratio was calculated for the whole protein.
Target Prediction of miRNAs
To better understand the biological function of DEmiRNA, TargetScan tool (Aggarwal et al., 2015) with default parameters was used to predict target genes of the selected DEmiRNAs (six miRNAs selected based on their role). From these predicted genes, the dysregulated genes from the proteomics data were identified for the miRNA selected (downregulated proteins for upregulated miRNA and upregulated proteins for downregulated miRNA). These common target genes from TargetScan and proteomics data were considered for further analysis. The miRNA–protein network was created based on the expression profile of target genes and miRNAs using Cytoscape (ver. 3.1.1; Shannon et al., 2003).
Gene Ontology Enrichment and Pathway Analysis
Functional annotation of the selected DEmiRNAs in each tissue was performed using target genes governed by them in ClueGO (ver. 2.1.4; Bindea et al., 2009) in Cytoscape (ver. 3.1.1; Shannon et al., 2003). Immune system processes and KEGG pathways were selected to generate a functionally organized GO/pathway term networks.
Validation Using qPCR
Total RNA, including small RNA from the lung and spleen of control and infected sheep and goats were isolated using mirVanaTM miRNA isolation kit (Invitrogen). Reverse transcriptase reactions were performed using RT specific primers of miR-363, miR-760-3p, miR-21-3p, and U6snRNA by TaqMan® MicroRNA Reverse Transcription Kit. Real-time PCR was performed using a standard TaqMan PCR kit protocol on an Applied Biosystems 7500 fast Sequence Detection System. The 10 μl PCR included 5 μl of 2× TaqMan Gene Expression Master Mix (Thermo Fisher Scientific Inc., Wilmington, DE, United States, Cat. No. 4369016), 0.5 μl of 20× TaqMan probe, 2 μl (0.134 ng) of RT product and 2.5 μl of NFW. The reactions were incubated in a 96-well plate at 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. All reactions were run in triplicate. The threshold cycle (Aad et al., 2015) is defined as the fractional cycle number at which the fluorescence passes the fixed threshold. The expression of the selected miRNAs described above was assayed taking the expression of U6snRNA as an internal control. The relative expression of each miRNA was calculated using the 2-ΔΔCT method with a control group as calibrator (Schmittgen and Livak, 2008).
Confirmation of PPRV Infection
Viral infection in the lung and spleen of sheep and goats infected with PPRV was confirmed by RT-PCR of 351 bp N gene amplicon in lung and spleen (Supplementary Figure S1). The viral infection was further confirmed by sandwich ELISA and qRT-PCR in both lung and spleen of goats and sheep (data not shown).
miRNAs Prediction and Identification of DEmiRNAs
In goats, miRDeep2 identified 298 and 283 miRNAs in control and PPRV infected lung, and 277 and 274 miRNAs in control and PPRV infected spleen. In sheep, 290 and 298 miRNAs in control and PPRV infected lung, and 274 and 256 miRNAs in control and PPRV infected spleen, respectively, were predicted. DEmiRNAs [FDR ≤ 0.05 and fold change (log2FC) ≥ 1] in lung and spleen of PPRV infected sheep and goats are presented in Table 1. A total of 67 miRNAs (34 downregulated and 33 upregulated) were dysregulated in the lungs of PPRV infected goats. However, a relatively small number of DEmiRNAs—37 miRNAs (16 miRNAs downregulated and 21 miRNAs upregulated) were identified in the lungs of sheep (Table 1). In infected goat’s spleen, 50 miRNAs were dysregulated with 26 of them downregulated and 24 upregulated. In spleen of sheep, 56 miRNAs were differentially expressed after PPRV infection and of these, 26 miRNAs were downregulated and 30 miRNAs were upregulated (Table 1).
TABLE 1. Differentially expressed (log2FC ≥ 1 and FDR ≤ 0.05 up/down regulated miRNAs in goats and sheep spleen and lung.
On comparing tissues across species, 20 and 11 miRNAs were found to be commonly differentially expressed in PPRV infected spleen and lung, respectively, in both sheep and goats. Among these 20 common DEmiRNAs in spleen, 11 DEmiRNAs (miR-199b, miR-1271, miR-217, miR-2887-1, miR-2887-2, miR-6119-3p, miR-221, miR-744, miR-30c, let-7a-5p-2, and miR-211) were downregulated and nine DEmiRNAs (miR-17-3p, miR-486, miR-146b, miR-363, miR-451, miR-193a-3p, miR-760-3p, miR-144, and miR-21-5p) were upregulated in goats, and in sheep, nine DEmiRNAs (miR-199b, miR-1271, miR-217, miR-6119-3p, miR-221, miR-744, miR-30c, let-7a-5p-2, and miR-211) were downregulated and 11 DEmiRNAs (miR-2887-1, miR-2887-2, miR-17-3p, miR-486, miR-146b, miR-363, miR-451, miR-193a-3p, miR-760-3p, miR-144, and miR-21-5p) were upregulated (Figure 1A and Table 2). Of these 20 common DEmiRNAs, miR-21-5p was the most upregulated (log2FC = 2.35) and miR-199b was the most downregulated DEmiRNA (log2FC = -3.03) in the spleen of goats. While in infected spleen of sheep the most upregulated DEmiRNAs were miR-451 (log2FC = 2.75) and miR-144 (log2FC = 2.63) and the most downregulated DEmiRNAs were miR-217 (log2FC = -4.08) and miR-221 (log2FC = -2.60).
FIGURE 1. Commonly differentially expressed miRNAs (log2FC > 1). (A) Commonly differentially expressed miRNAs (log2FC > 1) in spleen tissue of sheep and goats. (B) Commonly differentially expressed miRNAs (log2FC > 1) in lung tissue of sheep and goats.
The expression profile of the 11 common DEmiRNAs in the lung tissue varied in sheep and goats. miR-328 was found downregulated in goats but upregulated in sheep; two miRNAs—miR-2285f-2 and miR-27a-5p were found upregulated in goats but downregulated in sheep; six miRNAs—miR-320a-1, miR-320a-2, miR-1246, miR-363, miR-760-3p, and miR-21-3p were upregulated and two miRNAs—miR-34b and miR-150 were downregulated in both species (Figure 1B and Table 3). Among these 11 common DEmiRNAs, the expression of miR-21-3p, miR-760-3p, and miR-27a-5p was more abundant in lungs of PPRV infected goats with log2 fold change of 5.82, 3.79, and 3.07, respectively, while miR-34b (log2FC = -2.53) and miR-2285f-2 (log2FC = -2.53) were found least abundant in lungs of goats and sheep, respectively.
Among these 31 commonly DEmiRNAs, six DEmiRNAs were selected based on their role in viral infections, apoptosis, and fold change (Table 4). These miRNAs include miR-21-3p, miR-320a, miR-27a-5p, and miR-1246—expressed in lung of both species; miR-760-3p and miR-363—expressed in lung and spleen of both species (Figure 2). The miRNAs—miR-363 and miR-760-3p commonly present in PPRV infected lung and spleen of both species were identified to be upregulated. In lung, the expression of miR-363 and miR-760-3p was higher in goats with log2 fold change of 2.42 and 3.79, respectively, than in sheep (miR-363, log2FC = 1.04 and miR-760-3p, log2FC = 1.24) after PPRV infection. In spleen, the expression profile of miR-363 was higher in goats (log2FC = 1.55) than in sheep (log2FC = 1.28), however, no difference in expression of miR-760-3p was observed in goats (log2FC = 1.63) and sheep (log2FC = 1.67) infected with PPRV.
FIGURE 2. Overview of the miRNA prediction and identification of miRNAs that regulated dysregulated proteins.
Target Prediction and miRNA–Protein Regulatory Network Analysis
A total of 1149 (714 downregulated, 435 upregulated) and 1565 (1041 downregulated, 524 upregulated) differentially expressed proteins were identified in lung of sheep and goats, respectively, and 944 (281 downregulated, 663 upregulated) and 909 (590 downregulated, 319 upregulated) differentially expressed proteins were identified in spleen of sheep and goats, respectively. The number of dysregulated proteins identified through mass spectrophotometry by each of these six miRNAs is shown in Figure 2. The miRNA–protein interactions for each species and tissue are represented in a network (Figure 3). In the miRNA–protein network of lung tissue of goats, three miRNAs—miR-21-3p, miR-363, and miR-320a mutually regulate EGFR (epidermal growth factor receptor), which is involved in immune response. Similarly, IGF1R (insulin like growth factor 1) protein, involved in regulation of immune response was the target of miR-27a-5p, miR-363, miR-320a, and miR-760-3p. The TRIM (tripartite motif family) family members TRIM24, TRIM36, and TRIM45 were identified to be modulated by miR-1246, miR-320a, and miR-21-3p, respectively. The expression level of NF-κB signaling-related molecules IRAK2 and TRAF4 was regulated by miR-1246 and miR-320; and miR-1246 and miR-760-3p, respectively (Figure 3A). In PPRV infected sheep lung, the upregulated miRNAs—miR-21-3p and miR-320a govern immune genes—TRAF6, EGFR, and ERBB4 and the downregulated miR-27a-5p potentially modulates the expression of genes—MAP3K7 and MAPK8IP3, involved in JNK signaling pathways (Figure 3B).
FIGURE 3. miRNA-mRNA regulatory network. (A) Network of upregulated miRNAs (miR-21-3p, miR-27a-5p, miR-1246, miR-320a, miR-363, and miR-760-3p) and their downregulated genes in the lung tissue of goats. (B) Network of upregulated miRNAs (miR-21-3p, miR-1246, miR-320a, miR-363, and miR-760-3p) and their downregulated genes and the downregulated miRNA (miR-27a-5p) and its upregulated genes in the lung tissue of sheep. (C) Network of upregulated miRNAs (miR-363 and miR-760-3p) and their downregulated genes in the spleen tissue of goats. (D) Network of upregulated miRNAs (miR-363 and miR-760-3p) and their downregulated genes in the spleen tissue of sheep. Blue nodes indicate upregulation and red nodes indicates downregulation. Color intensity denoted the level of gene expression.
In PPRV infected spleen of goats, miR-363 and miR-760-3p regulate apoptotic molecules (NFATC2, NOX5, and MYO6) and target genes involved in immunological processes (NFATC2, IGF1R, KLHL21, and NOX5) (Figure 3C). The immune effector molecules IFIT5 and TRIM33 were targeted by miR-363, and CD244 and NFKBIE were regulated by miR-760-3p in infected spleen of sheep (Figure 3D).
Gene Ontology and KEGG Pathway-Based Network Analysis
Functional annotation of a total of 770 and 1226 target proteins governed by selected six miRNAs in infected sheep and goats infected lung, respectively, resulted in higher number of significantly enriched pathways and GO terms in goats than in sheep (Figure 4). The highly enriched common GO terms and pathways targeted by the miRNAs in the lung tissue of sheep and goats includes T cell receptor signaling pathway, Rap1 signaling pathway, Toll-like receptor TLR6:TLR2 signaling pathway, etc. (Figures 4A,B).
FIGURE 4. Gene ontology of immune-related KEGG pathways. (A) Gene ontology of immune-related KEGG pathways in sheep lung visualized in ClueGO (ver. 2.1.4) plugin of Cytoscape (ver. 3.1.1). (B) Gene ontology of immune-related KEGG pathways in goats lung visualized in ClueGO (ver. 2.1.4) plugin of Cytoscape (ver. 3.1.1). (C) Gene ontology of immune-related KEGG pathways in sheep spleen visualized in ClueGO (ver. 2.1.4) plugin of Cytoscape (ver. 3.1.1). (D) Gene ontology of immune-related KEGG pathways in goats spleen visualized in ClueGO (ver. 2.1.4) plugin of Cytoscape (ver. 3.1.1). Color from red to green of the nodes in the network depicts increase in significance. The diameter indicates the percent associated genes to a particular node.
Similarly, functional annotation of 84 and 173 target proteins governed by two out of selected six miRNAs in infected sheep and goats spleen, respectively, identified enrichment of B cell receptor signaling and FC epsilon RI signaling pathways in both the species (Figures 4C,D). Furthermore, targets of DEmiRNAs in sheep spleen were also found enriched in FC-gamma R mediated phagocytosis and myeloid leukocyte mediated immunity. The targets of DEmiRNAs of goat’s spleen were found involved in NF-κB signaling pathway, alpha–beta T cell differentiation, ErbB signaling pathway and regulation of B cell activation. In addition, the number of significantly enriched pathways and GO terms were higher in lung and spleen tissue of goats than in sheep.
Validation of DEmiRNAs by qPCR
To further validate the expression of DEmiRNAs from high-throughput sequencing, qPCR was performed on three DEmiRNAs—miR-21-3p, miR-363, and miR-760-3p. The expression of miR-363 and miR-760-3p in sheep and goats in both the tissues was in concordance with small RNA sequencing results. The expression of miR-21-3p was found to be in concordance with the sequencing results in PPRV infected spleen of both the species. However, in infected lung miR-21-3p was found upregulated on qPCR though not found in small RNA sequencing data in both the species (Figure 5 and Table 5).
FIGURE 5. qPCR validation of sequencing data. The change in the miRNA expression of miR-21-3p, miR-363, and miR-760-3p was calculated with U6snRNA as reference gene for normalization. Fold change is represented as log2FC.
PPR is a major threat to livestock keepers in developing countries, causing a severe disease in goats and sheep. Host encoded miRNAs have been demonstrated to be key regulators of host–virus interactions, and their expression is often affected by viral infection (Hussain and Asgari, 2010; Skalsky and Cullen, 2010). Currently there is no report available suggesting PPRV infection-induced changes in expression of cellular or host miRNAs. In the present study, we investigated the expression pattern of host miRNAs in spleen and lung of sheep and goats infected with PPRV vis-a-vis protein expression.
Detailed analysis revealed many differences in the global expression profile of miRNAs among lung and spleen, suggesting common and unique miRNA transcriptome landscape against PPRV. PPRV infection altered the expression of host miRNAs in lung and spleen. Under PPRV infection, a total of 37 and 67 DEmiRNAs were identified in lung of sheep and goats; and, 56 and 50 DEmiRNAs in the spleen of sheep and goats, respectively. PPRV infection in spleen and lung triggered the expression of many immune-related miRNAs, including, miR-21, miR-150, miR-146b, and let-7 family as reported in Japanese encephalitis virus infection (Cai et al., 2015). Moreover, 20 and 11 common DEmiRNAs expressed in spleen and lung of both species, respectively, suggested variable tissue response to PPRV infection. Among these 31 DEmiRNAs, six DEmiRNAs—miR-21-3p, miR-320a, miR-27a-5p, miR-1246 (expressed in lung of both species), miR-760-3p and miR-363 (expressed in lung and spleen of both species) were selected based on their role in viral infection, apoptosis and fold change. In infected goat’s lung all these six DEmiRNAs were found to upregulated. However, in infected sheep’s lung, miR-27a-5p was found to be downregulated and the rest of the five DEmiRNAs were upregulated.
miR-21-3p induce apoptosis (Lo et al., 2013) and PPRV is also reported to cause apoptosis of host cells (Mondal et al., 2001). The upregulated miR-363 is also known to induce apoptosis (Zhang et al., 2014; Zhou et al., 2014; Hu et al., 2015; Li et al., 2015). The upregulation of miR-21-3p and miR-363 in PPRV infections suggests synergistic effect of these miRs along with the virus in inducing apoptosis. The upregulation of miR-363 has also been reported in Human papillomavirus (HPV)+ and HPV- pharyngeal squamous cell carcinoma and in HPV16+ HNSCC cell lines (Lajer et al., 2011; Wald et al., 2011). Recently, miR-27a-5p was found to be highly expressed in vaccinia virus infection (Buck et al., 2010). Under PPRV infection miR-27a-5p was found to be upregulated in infected lung of goats but downregulated in sheep suggesting a species-specific response. Further, miR-320a is known to inhibit mink enteritis virus infection by downregulating its receptor, transferrin receptor (TfR; Sun et al., 2014). Significant upregulation of miR-320a expression in PPRV infected lung tissue of sheep and goats suggests that miR-320a might serve in triggering antiviral response against PPRV infection. Sheng et al. (2014) reported that the upregulated miR-1246 decreased the expression of cell adhesion target genes and hence promotes the cytotoxicity induced by Ebola virus glycoprotein. Similarly, upregulated miR-1246 was found to promote cell death pathway by reducing the expression levels of DLG3 protein during HEV71 infection in human neuroblastoma cells (Xu et al., 2014). The increased expression of miR-1246 in the PPRV infected lung of goats and sheep in our study, could be a factor contributing to the pathogenesis of PPRV.
For miRNA studies, it is critical to identify targets for understanding its biological function and molecular mechanism (Tang et al., 2015). The miRNA–protein network analysis suggests that one miRNA could participate in several biological processes by targeting different mRNAs, and one biological process could be influenced by multiple miRNAs. The upregulated miR-27a-5p, miR-363, miR-320a, and miR-760-3p were observed to bring about the downregulation of IGF1R in the PPRV infected lung of goats. IGF1R, is a multifunctional receptor that plays an important role in the regulation of immune response, including cell differentiation and proliferation (Smith, 2010). Similarly, TRIM24, tripartite motif containing 24, which is involved in cytokine signaling and secretion (Tisserand et al., 2011) was downregulated by miR-1246 and miR-320a in lung tissue of goats. miRNAs play an important role in regulation of NF-κB signaling pathway during viral infections (Gao et al., 2014) and activation of NF-κB is important for immune defense (Hoesel and Schmid, 2013). NF-κB signaling-related molecules IRAK2 were found modulated by miR-1246 and miR-320a, and TRAF4 was found modulated by miR-1246 and miR-760-3p in PPRV infected lung of goats in the present study. Interferon (IFN)-mediated pathway is a crucial part of the cellular response against viral infection (Wu et al., 2015). TRAF6, a major element in IFN production (Yoshida et al., 2008) was suppressed by PPRV-induced miR-21-3p and miR-320a in the lung of sheep. Similarly, the expression of IFIT5, which is involved in stimulating anti-viral response (Zhang et al., 2013) was suppressed by miR-363 in the spleen of sheep. This suggests that PPRV-induced miR-21-3p, miR-320a, and miR-363 might act cooperatively to enhance viral pathogenesis in the lung and spleen of sheep by downregulating several immune response genes. Further, this could be corroborated by the GO and pathway analysis of the potential targets of all the six DEmiRNAs. It was observed that the pathways governed by DEmiRNAs were more perturbed in goats than in sheep, thereby reflecting on the severity of disease in goats than in sheep.
This study demonstrated for the first time DEmiRNAs in sheep and goats under PPRV infection. The DEmiRNAs identified in this study govern genes involved in immune response processes. It was observed that PPRV elicits a strong host response in goats than in sheep as evident from the number of significantly enriched immune system pathways and genes perturbed. This study revealed that PPRV-induced miR-21-3p, miR-320a, and miR-363 might act cooperatively to enhance viral pathogenesis, which warrants further research.
RS, BPM, AT, and RG conceived and designed the research. ARS, SW, SS, SK, and AK conducted the wet lab work. AP, SW, ARS, and RG analyzed the data. AP, ARS, SW, and RG wrote the manuscript. AP, ARS, SW, RG, VS, APS, KR, BM, and DM helped in manuscript drafting and editing. RS, BPM, AT, and RG proofread the manuscript.
This study was supported in part by Centre for Agricultural Bioinformatics (ICAR-IASRI) and SubDIC (BTISnet), ICAR-IVRI.
Conflict of Interest Statement
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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2017.01146/full#supplementary-material
FIGURE S1 | Confirmation of PPRV infection in spleen and lung tissues of sheep and goats. Amplification of 351 bp N gene by RT-PCR. Lane M, 100 bp ladder; lane 1, NTC; lane 2, sheep spleen; lane 3, goats spleen; lane 4, goats lung; lane 5, sheep lung; lane 6, PPRV positive control.
Aad, G., Abbott, B., Abdallah, J., Abdinov, O., Aben, R., and Abolins, M. (2015). ATLAS Run 1 searches for direct pair production of third-generation squarks at the Large Hadron Collider. Eur. Phys. J C Part Fields 75, 510.
Altuvia, Y., Landgraf, P., Lithwick, G., Elefant, N., Pfeffer, S., Aravin, A., et al. (2005). Clustering and conservation patterns of human microRNAs. Nucleic Acids Res. 33, 2697–2706. doi: 10.1093/nar/gki567
Bindea, G., Mlecnik, B., Hackl, H., Charoentong, P., Tosolini, M., Kirilovsky, A., et al. (2009). ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 25, 1091–1093. doi: 10.1093/bioinformatics/btp101
Buck, A. H., Perot, J., Chisholm, M. A., Kumar, D. S., Tuddenham, L., Cognat, V., et al. (2010). Post-transcriptional regulation of miR-27 in murine cytomegalovirus infection. RNA 16, 307–315. doi: 10.1261/rna.1819210
Cai, Y., Zhu, L., Zhou, Y., Liu, X., Li, X., Lang, Q., et al. (2015). Identification and analysis of differentially-expressed microRNAs in Japanese encephalitis virus-infected PK-15 cells with deep sequencing. Int. J. Mol. Sci. 16, 2204–2219. doi: 10.3390/ijms16012204
Chang, J., Guo, J. T., Jiang, D., Guo, H., Taylor, J. M., and Block, T. M. (2008). Liver-specific microRNA miR-122 enhances the replication of hepatitis C virus in nonhepatic cells. J. Virol. 82, 8215–8223. doi: 10.1128/JVI.02575-07
Cui, L., Guo, X., Qi, Y., Qi, X., Ge, Y., Shi, Z., et al. (2010). Identification of microRNAs involved in the host response to enterovirus 71 infection by a deep sequencing approach. J. Biomed. Biotechnol. 2010:425939. doi: 10.1155/2010/425939
Delil, F., Asfaw, Y., and Gebreegziabher, B. (2012). Prevalence of antibodies to peste des petits ruminants virus before and during outbreaks of the disease in Awash Fentale district, Afar, Ethiopia. Trop. Anim. Health Prod. 44, 1329–1330. doi: 10.1007/s11250-012-0110-8
Dhar, P., Sreenivasa, B. P., Barrett, T., Corteyn, M., Singh, R. P., and Bandyopadhyay, S. K. (2002). Recent epidemiology of peste des petits ruminants virus (PPRV). Vet. Microbiol. 88, 153–159. doi: 10.1016/S0378-1135(02)00102-5
Dhinakar Raj, G., Nachimuthu, K., and Mahalinga Nainar, A. (2000). A simplified objective method for quantification of peste des petits ruminants virus or neutralizing antibody. J. Virol. Methods 89, 89–95. doi: 10.1016/S0166-0934(00)00206-8
Friedlander, M. R., Chen, W., Adamidi, C., Maaskola, J., Einspanier, R., Knespel, S., et al. (2008). Discovering microRNAs from deep sequencing data using miRDeep. Nat. Biotechnol. 26, 407–415. doi: 10.1038/nbt1394
Friedlander, M. R., Mackowiak, S. D., Li, N., Chen, W., and Rajewsky, N. (2012). miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 40, 37–52. doi: 10.1093/nar/gkr688
Guo, X. K., Zhang, Q., Gao, L., Li, N., Chen, X. X., and Feng, W. H. (2013). Increasing expression of microRNA 181 inhibits porcine reproductive and respiratory syndrome virus replication and has implications for controlling virus infection. J. Virol. 87, 1159–1171. doi: 10.1128/JVI.02386-12
Johnnidis, J. B., Harris, M. H., Wheeler, R. T., Stehling-Sun, S., Lam, M. H., Kirak, O., et al. (2008). Regulation of progenitor cell proliferation and granulocyte function by microRNA-223. Nature 451, 1125–1129. doi: 10.1038/nature06607
Lajer, C. B., Nielsen, F. C., Friis-Hansen, L., Norrild, B., Borup, R., Garnaes, E., et al. (2011). Different miRNA signatures of oral and pharyngeal squamous cell carcinomas: a prospective translational study. Br. J. Cancer 104, 830–840. doi: 10.1038/bjc.2011.29
Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10:R25. doi: 10.1186/gb-2009-10-3-r25
Li, S., Duan, X., Li, Y., Liu, B., Mcgilvray, I., and Chen, L. (2014). MicroRNA-130a inhibits HCV replication by restoring the innate immune response. J. Viral. Hepat. 21, 121–128. doi: 10.1111/jvh.12131
Lo, T. F., Tsai, W. C., and Chen, S. T. (2013). MicroRNA-21-3p, a berberine-induced miRNA, directly down-regulates human methionine adenosyltransferases 2A and 2B and inhibits hepatoma cell growth. PLoS ONE 8:e75628. doi: 10.1371/journal.pone.0075628
Maganga, G. D., Verrier, D., Zerbinati, R. M., Drosten, C., Drexler, J. F., and Leroy, E. M. (2013). Molecular typing of PPRV strains detected during an outbreak in sheep and goats in south-eastern Gabon in 2011. Virol. J. 10, 82. doi: 10.1186/1743-422X-10-82
Manjunath, S., Kumar, G. R., Mishra, B. P., Mishra, B., Sahoo, A. P., Joshi, C. G., et al. (2015). Genomic analysis of host - Peste des petits ruminants vaccine viral transcriptome uncovers transcription factors modulating immune regulatory pathways. Vet. Res. 46, 15. doi: 10.1186/s13567-015-0153-8
Manjunath, S., Mishra, B. P., Mishra, B., Sahoo, A. P., Tiwari, A. K., Rajak, K. K., et al. (2017). Comparative and temporal transcriptome analysis of peste des petits ruminants virus infected goat peripheral blood mononuclear cells. Virus Res. 229, 28–40. doi: 10.1016/j.virusres.2016.12.014
Mizuguchi, Y., Takizawa, T., and Uchida, E. (2015). Host cellular microRNA involvement in the control of hepatitis B virus gene expression and replication. World J. Hepatol. 7, 696–702. doi: 10.4254/wjh.v7.i4.696
Mondal, B., Sreenivasa, B. P., Dhar, P., Singh, R. P., and Bandyopadhyay, S. K. (2001). Apoptosis induced by peste des petits ruminants virus in goat peripheral blood mononuclear cells. Virus Res. 73, 113–119. doi: 10.1016/S0168-1702(00)00214-8
Nanda, Y. P., Chatterjee, A., Purohit, A. K., Diallo, A., Innui, K., Sharma, R. N., et al. (1996). The isolation of peste des petits ruminants virus from northern India. Vet. Microbiol. 51, 207–216. doi: 10.1016/0378-1135(96)00025-9
Raghavendra, A. G., Gajendragad, M. R., Sengupta, P. P., Patil, S. S., Tiwari, C. B., Balumahendiran, M., et al. (2008). Seroepidemiology of peste des petits ruminants in sheep and goats of southern peninsular India. Rev. Sci. Tech. 27, 861–867. doi: 10.20506/rst.27.3.1838
Roberts, A. P., Lewis, A. P., and Jopling, C. L. (2011a). miR-122 activates hepatitis C virus translation by a specialized mechanism requiring particular RNA components. Nucleic Acids Res. 39, 7716–7729. doi: 10.1093/nar/gkr426
Robinson, M. D., Mccarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Rodriguez, A., Vigorito, E., Clare, S., Warren, M. V., Couttet, P., Soond, D. R., et al. (2007). Requirement of bic/microRNA-155 for normal immune function. Science 316, 608–611. doi: 10.1126/science.1139253
Sahu, A. R., Wani, S. A., Padhy, A., Kumar, A., Priya, G. B., Sahoo, A. P., et al. (2015). Host-virus interaction: role of miRNA and bioinformatics tools for miRNA target prediction. Adv. Anim. Vet. Sci. 3, 30–36. doi: 10.14737/journal.aavs/2015/3.4s.30.36
Sahu, A. R., Wani, S. A., Saminathan, M., Rajak, K. K., Sahoo, A. P., Pandey, A., et al. (2017). Genome sequencing of an Indian peste des petits ruminants virus isolate, Izatnagar/94, and its implications for virus diversity, divergence and phylogeography. Arch. Virol. 162, 1677–1693. doi: 10.1007/s00705-017-3288-2
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
Sheng, M., Zhong, Y., Chen, Y., Du, J., Ju, X., Zhao, C., et al. (2014). Hsa-miR-1246, hsa-miR-320a and hsa-miR-196b-5p inhibitors can reduce the cytotoxicity of Ebola virus glycoprotein in vitro. Sci. China Life Sci. 57, 959–972. doi: 10.1007/s11427-014-4742-y
Singh, R. P., Saravanan, P., Sreenivasa, B. P., Singh, R. K., and Bandyopadhyay, S. K. (2004a). Prevalence and distribution of peste des petits ruminants virus infection in small ruminants in India. Rev. Sci. Tech. 23, 807–819.
Singh, R. P., Sreenivasa, B. P., Dhar, P., Shah, L. C., and Bandyopadhyay, S. K. (2004b). Development of a monoclonal antibody based competitive-ELISA for detection and titration of antibodies to peste des petits ruminants (PPR) virus. Vet. Microbiol. 98, 3–15.
Singh, R. P., Sreenivasa, B. P., Dhar, P., and Bandyopadhyay, S. K. (2004c). A sandwich-ELISA for the diagnosis of Peste des petits ruminants (PPR) infection in small ruminants using anti-nucleocapsid protein monoclonal antibody. Arch. Virol. 149, 2155–2170.
Sreenivasa, B. P., Dhar, P., Singh, R. P., and Bandyopadhyay, S. K. (2002). “Development of peste des petits ruminants (PPR) challenge virus from a field isolate,” in Proceedings of the XIV Annual Conference and National Seminar on Management of Viral Diseases with Emphasis on Global Trade and WTO Regime of Indian Virological Society, Bangalore.
Sun, J. Z., Wang, J., Wang, S., Yuan, D., Li, Z., Yi, B., et al. (2014). MicroRNA miR-320a and miR-140 inhibit mink enteritis virus infection by repression of its receptor, feline transferrin receptor. Virol. J. 11:210. doi: 10.1186/s12985-014-0210-3
Tang, Z., Yang, Y., Wang, Z., Zhao, S., Mu, Y., and Li, K. (2015). Integrated analysis of miRNA and mRNA paired expression profiling of prenatal skeletal muscle development in three genotype pigs. Sci. Rep. 5:15544. doi: 10.1038/srep15544
Tay, Y., Zhang, J., Thomson, A. M., Lim, B., and Rigoutsos, I. (2008). MicroRNAs to Nanog, Oct4 and Sox2 coding regions modulate embryonic stem cell differentiation. Nature 455, 1124–1128. doi: 10.1038/nature07299
Tisserand, J., Khetchoumian, K., Thibault, C., Dembele, D., Chambon, P., and Losson, R. (2011). Tripartite motif 24 (Trim24/Tif1alpha) tumor suppressor protein is a novel negative regulator of interferon (IFN)/signal transducers and activators of transcription (STAT) signaling pathway acting through retinoic acid receptor alpha (Raralpha) inhibition. J. Biol. Chem. 286, 33369–33379. doi: 10.1074/jbc.M111.225680
Trobaugh, D. W., Gardner, C. L., Sun, C., Haddow, A. D., Wang, E., Chapnik, E., et al. (2014). RNA viruses can hijack vertebrate microRNAs to suppress innate immunity. Nature 506, 245–248. doi: 10.1038/nature12869
Truong, T., Boshra, H., Embury-Hyatt, C., Nfon, C., Gerdts, V., Tikoo, S., et al. (2014). Peste des petits ruminants virus tissue tropism and pathogenesis in sheep and goats following experimental infection. PLoS ONE 9:e87145. doi: 10.1371/journal.pone.0087145
Wald, A. I., Hoskins, E. E., Wells, S. I., Ferris, R. L., and Khan, S. A. (2011). Alteration of microRNA profiles in squamous cell carcinoma of the head and neck cell lines by human papillomavirus. Head Neck 33, 504–512. doi: 10.1002/hed.21475
Wang, Y., Brahmakshatriya, V., Zhu, H., Lupiani, B., Reddy, S. M., Yoon, B. J., et al. (2009). Identification of differentially expressed miRNAs in chicken lung and trachea with avian influenza virus infection by a deep sequencing approach. BMC Genomics 10:512. doi: 10.1186/1471-2164-10-512
Xu, L. J., Jiang, T., Zhao, W., Han, J. F., Liu, J., Deng, Y. Q., et al. (2014). Parallel mRNA and microRNA profiling of HEV71-infected human neuroblastoma cells reveal the up-regulation of miR-1246 in association with DLG3 repression. PLoS ONE 9:e95272. doi: 10.1371/journal.pone.0095272
Yoshida, R., Takaesu, G., Yoshida, H., Okamoto, F., Yoshioka, T., Choi, Y., et al. (2008). TRAF6 and MEKK1 play a pivotal role in the RIG-I-like helicase antiviral pathway. J. Biol. Chem. 283, 36211–36220. doi: 10.1074/jbc.M806576200
Zhang, B., Liu, X., Chen, W., and Chen, L. (2013). IFIT5 potentiates anti-viral response through enhancing innate immune signaling pathways. Acta Biochim. Biophys. Sin. 45, 867–874. doi: 10.1093/abbs/gmt088
Zhou, P., Huang, G., Zhao, Y., Zhong, D., Xu, Z., Zeng, Y., et al. (2014). MicroRNA-363-mediated downregulation of S1PR1 suppresses the proliferation of hepatocellular carcinoma cells. Cell. Signal. 26, 1347–1354. doi: 10.1016/j.cellsig.2014.02.020
Keywords: microRNAs, PPR, sheep, goats, host–pathogen interaction, molecular pathogenesis
Citation: Pandey A, Sahu AR, Wani SA, Saxena S, Kanchan S, Sah V, Rajak KK, Khanduri A, Sahoo AP, Tiwari AK, Mishra B, Muthuchelvan D, Mishra BP, Singh RK and Gandham RK (2017) Modulation of Host miRNAs Transcriptome in Lung and Spleen of Peste des Petits Ruminants Virus Infected Sheep and Goats. Front. Microbiol. 8:1146. doi: 10.3389/fmicb.2017.01146
Received: 23 March 2017; Accepted: 06 June 2017;
Published: 26 June 2017.
Edited by:Lars Kaderali, Technische Universität Dresden, Germany
Reviewed by:Muhammad Zubair Shabbir, University of Veterinary and Animal Sciences, Pakistan
Balamurugan Vinayagamurthy, National Institute of Veterinary Epidemiology and Disease Informatics (ICAR), India
Copyright © 2017 Pandey, Sahu, Wani, Saxena, Kanchan, Sah, Rajak, Khanduri, Sahoo, Tiwari, Mishra, Muthuchelvan, Mishra, Singh and Gandham. 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) or licensor 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: Ravi K. Gandham, firstname.lastname@example.org
†These authors have contributed equally to this work.