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.

Systems biology under heat stress in Indian cattle.

Transcriptome profiling of Vrindavani and Tharparkar cattle (n = 5 each) revealed that more numbers of genes were dysregulated in Vrindavani than in Tharparkar. A contrast in gene expression was observed with 18.9 % of upregulated genes in Vrindavani downregulated in Tharparkar and 17.8% upregulated genes in Tharparkar downregulated in Vrindavani. Functional annotation of genes differentially expressed in Tharparkar and Vrindavani revealed that the systems biology in Tharparkar is moving towards counteracting the effects due to heat stress. Unlike Vrindavani, Tharparkar is not only endowed with higher expression of the scavengers (UBE2G1, UBE2S, and UBE2H) of misfolded proteins but also with protectors (VCP, Serp1, and CALR) of naïve unfolded proteins. Further, higher expression of the antioxidants in Tharparkar enables it to cope up with higher levels of free radicals generated as a result of heat stress. In this study, we found relevant genes dysregulated in Tharparkar in the direction that can counter heat stress.

Signature of genome wide gene expression in classical swine fever virus infected macrophages and PBMCs of indigenous vis-a-vis crossbred pigs.

The genetic basis of differential host immune response vis-à-vis transcriptome profile was explored in PBMCs of indigenous (Ghurrah) and crossbred pigs after classical swine fever vaccination and in monocyte derived macrophages (MDMs) challenged with virulent classical swine fever (CSF) virus. The humoral immune response (E2 antibody) was higher (74.87%) in crossbred than indigenous pigs (58.20%) at 21st days post vaccination (21dpv). The rate of reduction of ratio of CD4+/CD8+ was higher in crossbred pigs than indigenous pigs at 7th days post vaccination (7dpv). The immune genes IFIT1, IFIT5, RELA, NFKB2, TNF and LAT2 were up regulated at 7dpv in RNA seq data set and was in concordance during qRT-PCR validation. The Laminin Subunit Beta 1 (LAMB1) was significantly (p ≤ 0.05) down-regulated in MDMs of indigenous pigs and consequently a significantly (p ≤ 0.01) higher copy number of virulent CSF virus was evidenced in macrophages of crossbred pigs than indigenous pigs. Activation of LXR:RXR pathway at 60 h post infection (60hpi) in MDMs of indigenous versus crossbred pigs inhibited nuclear translocation of NF-κB, resulted into transrepression of proinflammatory genes. But it helped in maintenance of HDL level by lowering down cholesterol/LDL level in MDMs of indigenous pigs. The key immune genes (TLR2, TLR4, IL10, IL8, CD86, CD54, CASP1) of TREM1 signaling pathway were upregulated at 7dpv in PBMCs but those genes were downregulated at 60hpi in MDMs indigenous pigs. Using qRT-PCR, the validation of differentially expressed, immunologically important genes (LAMB1, OAS1, TLR 4, TLR8 and CD86) in MDMs revealed that expression of these genes were in concordance with RNA-seq data.

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.

Genome-wide integrated analysis of miRNA and mRNA expression profiles to identify differentially expressed miR-22-5p and miR-27b-5p in response to classical swine fever vaccine virus.

The present study was conducted to identify the differentially expressed miRNAs (DE miRNAs) in the peripheral blood mononuclear cells of crossbred pigs in response to CSF vaccination on 7 and 21 days of post vaccination as compared to unvaccinated control (0 dpv). Simultaneously, set of miRNA was predicted using mRNA seq data at same time point. The proportion of CD4-CD8+ and CD4+CD8+ increased after vaccination, and the mean percentage inhibition was 86.89% at 21 dpv. It was observed that 22 miRNAs were commonly expressed on both the time points. Out of predicted DE miRNAs, it was found that 40 and 35 DE miRNAs were common, obtained from miRNA seq analysis and predicted using mRNA seq data on 7 dpv versus 0 dpv and 21 dpv versus 0 dpv respectively. Two DE miRNAs, ssc-miR-22-5p and ssc-miR-27b-5p, were selected based on their log2 fold change and functions of their target genes in immune process/pathway of viral infections. The validations of DE miRNAs using qRT-PCR were in concordance with miRNA seq analysis. Two set of target genes, CD40 and SWAP70 (target gene of ssc-miR-22-5p) and TLR4 and Lyn (target gene of ssc-miR-27b-5p), were validated and were in concordance with results of RNA seq analysis at a particular time point (except TLR4). The first report of genome-wide identification of differentially expressed miRNA in response to live attenuated vaccine virus of classical swine fever revealed miR-22-5p and miR-27b-5p were differentially expressed at 7 dpv and 21 dpv.

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.

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.

Draft genome sequence of field isolate Brucella melitensis strain 2007BM/1 from India.

OBJECTIVES: Brucellosis is among one of the most widespread important global zoonotic diseases that is endemic in many parts of India. Brucella melitensis is supposed to be the most pathogenic species for humans. Here we report the draft genome sequence of B. melitensis strain 2007BM/1 isolated from a human in India. METHODS: Genomic DNA was extracted from Brucella culture and was sequenced using an Illumina MiSeq platform. The generated reads were assembled using three de novo assemblers and the draft genome was annotated. RESULTS: This monoisolate, with a genome length of 3268756bp, was found to be resistant to azithromycin and trimethoprim/sulfamethoxazole but susceptible to tetracycline, ofloxacin, rifampicin, ciprofloxacin and doxycycline. The presence of virulence genes in the strain was identified. CONCLUSIONS: The results obtained will help in understanding drug resistance mechanisms and virulence factors in highly zoonotic B. melitensis and suggest the need for judicious use of antibiotics in livestock health and management practices.

Expression kinetics of ISG15, IRF3, IFNγ, IL10, IL2 and IL4 genes vis-a-vis virus shedding, tissue tropism and antibody dynamics in PPRV vaccinated, challenged, infected sheep and goats.

Here, we studied the in vivo expression of Th1 (IL2 and IFN gamma) and Th2 (IL4 and IL10) - cytokines and antiviral molecules - IRF3 and ISG15 in peripheral blood mononuclear cells in relation to antigen and antibody dynamics under Peste des petits ruminants virus (PPRV) vaccination, infection and challenge in both sheep and goats. Vaccinated goats were seropositive by 9 days post vaccination (dpv) while in sheep idiosyncratic response was observed between 9 and 14 dpv for different animals. Expression of PPRV N gene was not detected in PBMCs of vaccinated and vaccinated challenged groups of both species, but was detected in unvaccinated infected PBMCs at 9 and 14 days post infection. The higher viral load at 9 dpi coincided with the peak clinical signs of the disease. The peak in viral replication at 9 dpi correlated with significant expression of antiviral molecules IRF3, ISG15 and IFN gamma in both the species. With the progression of disease, the decrease in N gene expression also correlated with the decrease in expression of IRF3, ISG15 and IFN gamma. In the unvaccinated infected animals ISG15, IRF3, IFN gamma and IL10 expression was higher than vaccinated animals. The IFN gamma expression predominated over IL4 in both vaccinated and infected animals with the infected exhibiting a stronger Th1 response. The persistent upregulation of this antiviral molecular signature - ISG15 and IRF3 even after 2 weeks post vaccination, presumably reflects the ongoing stimulation of innate immune cells.

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.

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

Peste des petits ruminants is an important transboundary disease infecting small ruminants. Genome or gene sequence analysis enriches our knowledge about the evolution and transboundary nature of the causative agent of this disease, peste des petits ruminants virus (PPRV). Although analysis using whole genome sequences of pathogens leads to more precise phylogenetic relationships, when compared to individual genes or partial sequences, there is still a need to identify specific genes/genomic regions that can provide evolutionary assessments consistent with those predicted with full-length genome sequences. Here the virulent Izatnagar/94 PPRV isolate was assembled and compared to all available complete genome sequences (currently in the NCBI database) to estimate nucleotide diversity and to deduce evolutionary relationships between genes/genomic regions and the full length genomes. Our aim was to identify the preferred candidate gene for use as a phylogenetic marker, as well as to predict divergence time and explore PPRV phylogeography. Among all the PPRV genes, the H gene was identified to be the most diverse with the highest evolutionary relationship with the full genome sequences. Hence it is considered as the most preferred candidate gene for phylogenetic study with 93% identity set as a nucleotide cutoff. A whole genome nucleotide sequence cutoff value of 94% permitted specific differentiation of PPRV lineages. All the isolates examined in the study were found to have a most recent common ancestor in the late 19th or in the early 20th century with high posterior probability values. The Bayesian skyline plot revealed a decrease in genetic diversity among lineage IV isolates since the start of the vaccination program and the network analysis localized the ancestry of PPRV to Africa.

Comparative and temporal transcriptome analysis of peste des petits ruminants virus infected goat peripheral blood mononuclear cells.

Peste des petits ruminanats virus (PPRV), a morbillivirus causes an acute, highly contagious disease - peste des petits ruminants (PPR), affecting goats and sheep. Sungri/96 vaccine strain is widely used for mass vaccination programs in India against PPR and is considered the most potent vaccine providing long-term immunity. However, occurrence of outbreaks due to emerging PPR viruses may be a challenge. In this study, the temporal dynamics of immune response in goat peripheral blood mononuclear cells (PBMCs) infected with Sungri/96 vaccine virus was investigated by transcriptome analysis. Infected goat PBMCs at 48h and 120h post infection revealed 2540 and 2000 differentially expressed genes (DEGs), respectively, on comparison with respective controls. Comparison of the infected samples revealed 1416 DEGs to be altered across time points. Functional analysis of DEGs reflected enrichment of TLR signaling pathways, innate immune response, inflammatory response, positive regulation of signal transduction and cytokine production. The upregulation of innate immune genes during early phase (between 2-5 days) viz. interferon regulatory factors (IRFs), tripartite motifs (TRIM) and several interferon stimulated genes (ISGs) in infected PBMCs and interactome analysis indicated induction of broad-spectrum anti-viral state. Several Transcription factors - IRF3, FOXO3 and SP1 that govern immune regulatory pathways were identified to co-regulate the DEGs. The results from this study, highlighted the involvement of both innate and adaptive immune systems with the enrichment of complement cascade observed at 120h p.i., suggestive of a link between innate and adaptive immune response. Based on the transcriptome analysis and qRT-PCR validation, an in vitro mechanism for the induction of ISGs by IRFs in an interferon independent manner to trigger a robust immune response was predicted in PPRV infection.

Systems biology approach: Panacea for unravelling host-virus interactions and dynamics of vaccine induced immune response.

Systems biology is an interdisciplinary research field in life sciences, which involves a comprehensive and quantitative analysis of the interactions between all of the components of biological systems over time. For the past 50 years the discipline of virology has overly focused on the pathogen itself. However, we now know that the host response is equally or more important in defining the eventual pathological outcome of infection. Systems biology has in recent years been increasingly recognised for its importance to infectious disease research. Host-virus interactions can be better understood by taking into account the dynamical molecular networks that constitute a biological system. To decipher the pathobiological mechanisms of any disease requires a deep knowledge of how multiple and concurrent signal-transduction pathways operate and are deregulated. Hence the intricacies of signalling pathways can be dissected only by system level approaches.