Genome sequencing of an Indian peste des petits ruminants virus isolate, Izatnagar/94, and its implications for virus diversity, divergence and phylogeography

Changes in m6A RNA methylation of goat lung following PPRV infection

Peste des petits ruminants (PPR) is an acute, highly contagious viral disease of goats and sheep, caused by the Peste des petits ruminants virus (PPRV). Earlier studies suggest the involvement of diverse regulatory mechanisms in PPRV infection. Methylation at N6 of Adenosine called m6A is a type RNA modification that influences various physiological and pathological phenomena. As the lung tissue represents the primary target organ of PPRV, the present study explored the m6A changes and their functional significance in PPRV disease pathogenesis. m6A-seq analysis revealed 1289 m6A peaks to be significantly altered in PPRV infected lung in comparison to normal lung, out of which 975 m6A peaks were hypomethylated and 314 peaks were hypermethylated. Importantly, hypomethylated genes were enriched in Interleukin-4 and Interleukin-13 signaling and various processes associated with extracellular matrix organization. Further, of the 843 differentially m6A-containing cellular transcripts, 282 transcripts were also found to be differentially expressed. Functional analysis revealed that these 282 transcripts are significantly enriched in signaling by Interleukins, extracellular matrix organization, cytokine signaling in the immune system, signaling by receptor tyrosine kinases, and Toll-like Receptor Cascades. We also found m6A reader HNRNPC and the core component of methyltransferase complex METTL14 to be highly upregulated than the m6A readers - HNRNPA2B1 and YTHDF1 at the transcriptome level. These findings suggest that alteration in the m6A landscape following PPRV is implicated in diverse processes including Interleukin signaling.

Differential expression of long non-coding RNAs under Peste des petits ruminants virus (PPRV) infection in goats

Peste des petits ruminants (PPR) characterized by fever, sore mouth, conjunctivitis, gastroenteritis, and pneumonia, is an acute, highly contagious viral disease of sheep and goats. The role of long non-coding RNAs (lncRNAs) in PPRV infection has not been explored to date. In this study, the transcriptome profiles of virulent Peste des petits ruminants virus (PPRV) infected goat tissues – lung and spleen were analyzed to identify the role of lncRNAs in PPRV infection. A total of 13,928 lncRNA transcripts were identified, out of which 170 were known lncRNAs. Intergenic lncRNAs (7625) formed the major chunk of the novel lncRNA transcripts. Differential expression analysis revealed that 15 lncRNAs (11 downregulated and 4 upregulated) in the PPRV infected spleen samples and 16 lncRNAs (13 downregulated and 3 upregulated) in PPRV infected lung samples were differentially expressed as compared to control. The differentially expressed lncRNAs (DElncRNAs) possibly regulate various immunological processes related to natural killer cell activation, antigen processing and presentation, and B cell activity, by regulating the expression of mRNAs through the cis- or trans-regulatory mechanism. Functional enrichment analysis of differentially expressed mRNAs (DEmRNAs) revealed enrichment of immune pathways and biological processes in concordance with the pathways in which correlated lncRNA-neighboring genes were enriched. The results suggest that a coordinated immune response is raised in both lung and spleen tissues of the goat through mRNA-lncRNA crosstalk.

Ongoing Assessment of the Molecular Evolution of Peste Des Petits Ruminants Virus Continues to Question Viral Origins

Understanding the evolution of viral pathogens is critical to being able to define how viruses emerge within different landscapes. Host susceptibility, which is spread between different species and is a contributing factor to the subsequent epidemiology of a disease, is defined by virus detection and subsequent characterization. Peste des petits ruminants virus is a plague of small ruminant species that is a considerable burden to the development of sustainable agriculture across Africa and much of Asia. The virus has also had a significant impact on populations of endangered species in recent years, highlighting its significance as a pathogen of high concern across different regions of the globe. Here, we have re-evaluated the molecular evolution of this virus using novel genetic data to try and further resolve the molecular epidemiology of this disease. Viral isolates are genetically characterized into four lineages (I-IV), and the historic origin of these lineages is of considerable interest to the molecular evolution of the virus. Our re-evaluation of viral emergence using novel genome sequences has demonstrated that lineages I, II and IV likely originated in West Africa, in Senegal (I) and Nigeria (II and IV). Lineage III sequences predicted emergence in either East Africa (Ethiopia) or in the Arabian Peninsula (Oman and/or the United Arab Emirates), with a paucity of data precluding a more refined interpretation. Continual refinements of evolutionary emergence, following the generation of new data, is key to both understanding viral evolution from a historic perspective and informing on the ongoing genetic emergence of this virus.

A comprehensive global perspective on phylogenomics and evolutionary dynamics of Small ruminant morbillivirus

A string of complete genome sequences of Small ruminant morbillivirus (SRMV) have been reported from different parts of the globe including Asia, Africa and the Middle East. Despite individual genome sequence-based analysis, there is a paucity of comparative genomic and evolutionary analysis to provide overarching and comprehensive evolutionary insights. Therefore, we first enriched the existing database of complete genome sequences of SRMVs with Pakistan-originated strains and then explored overall nucleotide diversity, genomic and residue characteristics, and deduced an evolutionary relationship among strains representing a diverse geographical region worldwide. The average number of pairwise nucleotide differences among the whole genomes was found to be 788.690 with a diversity in nucleotide sequences (0.04889 ± S.D. 0.00468) and haplotype variance (0.00001). The RNA-dependent-RNA polymerase (L) gene revealed phylogenetic relationship among SRMVs in a pattern similar to those of complete genome and the nucleoprotein (N) gene. Therefore, we propose another useful molecular marker that may be employed for future epidemiological investigations. Based on evolutionary analysis, the mean evolution rate for the complete genome, N, P, M, F, H and L genes of SRMV was estimated to be 9.953 × 10^–4, 1.1 × 10^–3, 1.23 × 10^–3, 2.56 × 10^–3, 2.01 × 10^–3, 1.47 × 10^–3 and 9.75 × 10^–4 substitutions per site per year, respectively. A recombinant event was observed in a Pakistan-originated strain (KY967608) revealing Indian strains as major (98.1%, KR140086) and minor parents (99.8%, KT860064). Taken together, outcomes of the study augment our knowledge and current understanding towards ongoing phylogenomic and evolutionary dynamics for better comprehensions of SRMVs and effective disease control interventions.

Epidemic and evolutionary characteristics of peste des petits ruminants virus infecting Procapra przewalskii in Western China

Due to the migration or transboundary spread of domestic and wild animals, peste des petits ruminants virus posed a high potential threat to them. In this study, we initially detected that a class of animal named Procapra przewalskii was infected with peste des petits ruminants virus (PPRV ChinaGS2018) in Gansu province. According to phylogenetic relationships analysis, we found that ChinaGS2018 comprised of 15,954 nucleotides and was classified into IV genotypes. In addition, indirect immunofluorescence assay (IFA) showed that ChinaGS2018 could infect isolated primary goat tracheal epithelium cells (GTC). Comparing with full-length genome sequences revealed that ChinaGS2018 strain has high identity to the reference complete genomes (87.16-99.55%) at the nucleotide level. Multiple sequence alignment showed that F protein has the highest identity of 99.8%, and H protein has the highest nucleotide substitution ratio. Our study also suggested this strain may be transmitted from Xinjiang, China. Along with the migratory of Procapraprzewalskii, this wild ruminant infected with PPRV can pose a huge threat to other wild ruminants and domestic ones. This is the first report describing infected with PPRV which will provide insights into the epidemiology and pathogenesis of this important virus.

A comparative phylogenomic analysis of peste des petits ruminants virus isolated from wild and unusual hosts

Peste des petits ruminants virus (PPRV) infects a wide range of domestic and wild ruminants, and occasionally unusual hosts such as camel, cattle and pig. Given their broad host-spectrum and disease endemicity in several developing countries, it is imperative to elucidate the viral evolutionary insights for their dynamic pathobiology and differential host-selection. For this purpose, a dataset of all available (n = 37) PPRV sequences originating from wild and unusual hosts was composed and in silico analysed. Compared to domestic small ruminant strains of same geographical region, phylogenomic and residue analysis of PPRV sequences originating from wild and unusual hosts revealed a close relationship between strains. A lack of obvious difference among the studied sequences and deduced residues suggests that these are the host factors that may play a role in their susceptibility to PPRV infection, immune response, pathogenesis, excretion patterns and potential clinical signs or resistance to clinical disease. Summarizing together, the comparative analysis enhances our understanding towards molecular epidemiology of the PPRV in wild and unusual hosts for appropriate intervention strategies particularly at livestock-wildlife interface.

Contrasting Gene Expression Profiles of Monocytes and Lymphocytes From Peste-Des-Petits-Ruminants Virus Infected Goats

In this study, transcriptome analysis of PPRV infected PBMC subsets-T helper cells, T cytotoxic cells, monocytes, and B lymphocytes was done to delineate their role in host response. PPRV was found to infect lymphocytes and not monocytes. The established receptor for PPRV-SLAM was found downregulated in lymphocytes and non-differentially expressed in monocytes. A profound deviation in the global gene expression profile with a large number of unique upregulated genes (851) and downregulated genes (605) was observed in monocytes in comparison to lymphocytes. ISGs-ISG15, Mx1, Mx2, RSAD2, IFIT3, and IFIT5 that play a role in antiviral response and the genes for viral sensors-MDA5, LGP2, and RIG1, were found to be upregulated in lymphocytes and downregulated in monocytes. The transcription factors-IRF-7 and STAT-1 that regulate expression of most of the ISGs were found activated in lymphocytes and not in monocytes. Interferon signaling pathway and RIG1 like receptor signaling pathway were found activated in lymphocytes and not in monocytes. This contrast in gene expression profiles and signaling pathways indicated the predominant role of lymphocytes in generating the antiviral response against PPRV in goats, thus, giving us new insights into host response to PPRV.

Dysregulated miRNAome and Proteome of PPRV Infected Goat PBMCs Reveal a Coordinated Immune Response

In this study, the miRNAome and proteome of virulent Peste des petits ruminants virus (PPRV) infected goat peripheral blood mononuclear cells (PBMCs) were analyzed. The identified differentially expressed miRNAs (DEmiRNAs) were found to govern genes that modulate immune response based on the proteome data. The top 10 significantly enriched immune response processes were found to be governed by 98 genes. The top 10 DEmiRNAs governing these 98 genes were identified based on the number of genes governed by them. Out of these 10 DEmiRNAs, 7 were upregulated, and 3 were downregulated. These include miR-664, miR-2311, miR-2897, miR-484, miR-2440, miR-3533, miR-574, miR-210, miR-21-5p, and miR-30. miR-664 and miR-484 with proviral and antiviral activities, respectively, were upregulated in PPRV infected PBMCs. miR-210 that inhibits apoptosis was downregulated. miR-21-5p that decreases the sensitivity of cells to the antiviral activity of IFNs and miR-30b that inhibits antigen processing and presentation by primary macrophages were downregulated, indicative of a strong host response to PPRV infection. miR-21-5p was found to be inhibited on IPA upstream regulatory analysis of RNA-sequencing data. This miRNA that was also highly downregulated and was found to govern 16 immune response genes in the proteome data was selected for functional validation vis-a-vis TGFBR2 (TGF-beta receptor type-2). TGFBR2 that regulates cell differentiation and is involved in several immune response pathways was found to be governed by most of the identified immune modulating DEmiRNAs. The decreased luciferase activity in Dual Luciferase Reporter Assay indicated specific binding of miR-21-5p and miR-484 to their target thus establishing specific binding of the miRNAs to their targets.This is the first report on the miRNAome and proteome of virulent PPRV infected goat PBMCs.

VARV B22R homologue as phylogenetic marker gene for Capripoxvirus classification and divergence time dating

Sheeppox disease is associated with significant losses in sheep production world over. The sheep pox virus, the goatpox virus, and the lumpy skin disease virus cannot be distinguished by conventional serological tests. Identification of these pathogens needs molecular methods. In this study, seven genes viz. EEV maturation protein-F12L, Virion protein-D3R, RNA polymerase subunit-A5R, Virion core protein-A10L, EEV glycoprotein-A33R, VARV B22R homologue, and Kelch like protein-A55R that cover the start, middle, and end of the genome were selected. These genes were amplified from Roumanian-Fanar vaccine strain and Jaipur virulent strain, cloned, and sequenced. On analysis with the available database sequences, VARV B22R homologue was identified as a marker for phylogenetic reconstruction for classifying the sheeppox viruses of the ungulates. Further, divergence time dating with VARV B22R gene accurately predicted the sheeppox disease outbreak involving Jaipur virulent strain.

Selection and validation of suitable reference genes for qPCR gene expression analysis in goats and sheep under Peste des petits ruminants virus (PPRV), lineage IV infection

Identification of suitable candidate reference genes is an important prerequisite for validating the gene expression data obtained from downstream analysis of RNA sequencing using quantitative real time PCR (qRT-PCR). Though existence of a universal reference gene is myth, commonly used reference genes can be assessed for expression stability to confer their suitability to be used as candidate reference genes in gene expression studies. In this study, we evaluated the expression stability of ten most commonly used reference genes (GAPDH, ACTB, HSP90, HMBS, 18S rRNA, B2M, POLR2A, HPRT1, ACAC, YWHAZ) in fourteen different Peste des petits ruminants virus (PPRV) infected tissues of goats and sheep. RefFinder and RankAggreg software were used to deduce comprehensive ranking of reference genes. Our results suggested HMBS and B2M in goats and HMBS and HPRT1 in sheep can be used as suitable endogenous controls in gene expression studies of PPRV infection irrespective of tissues and condition as a whole, thus eliminating the use of tissue specific/ condition specific endogenous controls. We report for the first time suitable reference genes for gene expression studies in PPRV infected tissues. The reference genes determined here can be useful for future studies on gene expression in sheep and goat infected with PPRV, thus saving extra efforts and time of repeating the reference gene determination and validation.

Modulation of Host miRNAs Transcriptome in Lung and Spleen of Peste des Petits Ruminants Virus Infected Sheep and Goats

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.