Transcriptome profiling of buffalo endometrium reveals molecular signature distinct to early pregnancy.

Buffalo reproduction struggles with a high incidence of early embryonic mortality. Effective treatment and prevention strategies for this condition are not available due to lack of understanding of molecular pathways in early pregnancy of this species. In the present study, we have attempted to understand these molecular pathways by characterizing the endometrial transcriptomic profiles of pregnant buffalos during early pregnancy. For the transcriptome profiling, buffalo endometrial tissues of 29-36 days of pregnancy and of nonpregnant luteal phase were collected from the local slaughterhouse. We confirmed the status of pregnancy based on the crown vertebral length of the foetus. Total RNA was isolated and sequencing was performed using the Illumina nextseq platform. The raw reads were filtered and mapped to the Bos taurus UMD 3.1 reference genome assembly. An average of 24,597 genes was investigated for differential expression between the two groups. Transcriptome data identified a total of 450 differentially expressed genes (using a cut off value of log2 fold changes >2 and <-2) in early pregnancy in comparison to the nonpregnant group (Padj < 0.05). Among these, 270 genes were significantly upregulated and 180 genes were downregulated. The most impacted pathways were related to secretion, transport, ionic homeostasis, mitosis and negative regulation of viral processes. In conclusion, our study characterized a unique set of DEGs, during the early pregnancy of buffalo, which potentially modulate the endometrial environment to establish and maintain a successful pregnancy.

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.

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.

RNA Seq analysis for transcriptome profiling in response to classical swine fever vaccination in indigenous and crossbred pigs.

In present investigation, differential expression of transcriptome after classical swine fever (CSF) vaccination has been explored at the cellular level in crossbred and indigenous (desi) piglets. RNA Sequencing by Expectation-Maximization (RSEM) package was used to quantify gene expression from RNA Sequencing data, and differentially expressed genes (DEGs) were identified using EBSeq, DESeq2, and edgeR softwares. After analysis, 5222, 6037, and 6210 common DEGs were identified in indigenous post-vaccinated verses pre-vaccinated, crossbred post-vaccinated verses pre-vaccinated, and post-vaccinated crossbred verses indigenous pigs, respectively. Functional annotation of these DEGs showed enrichment of antigen processing-cross presentation, B cell receptor signaling, T cell receptor signaling, NF-κB signaling, and TNF signaling pathways. The interaction network among the immune genes included more number of genes with greater connectivity in vaccinated crossbred than the indigenous piglets. Higher expression of IRF3, IL1ß, TAP1, CSK, SLA2, SLADM, and NF-kB in crossbred piglets in comparison to indigenous explains the better humoral response observed in crossbred piglets. Here, we predicted that the processed CSFV antigen through the T cell receptor signaling cascade triggers the B cell receptor-signaling pathway to finally activate MAPK kinase and NF-κB signaling pathways in B cell. This activation results in expression of genes/transcription factors that lead to B cell ontogeny, auto immunity and immune response through antibody production. Further, immunologically important genes were validated by qRT-PCR.

Assessment of expected breeding values for fertility traits of Murrah buffaloes under subtropical climate.

AIM: The aim of the present study was to assess the influence of temperature and humidity prevalent under subtropical climate on the breeding values for fertility traits viz. service period (SP), pregnancy rate (PR) and conception rate (CR) of Murrah buffaloes in National Dairy Research Institute (NDRI) herd. MATERIALS AND METHODS: Fertility data on 1379 records of 581 Murrah buffaloes spread over four lactations and climatic parameters viz. dry bulb temperature and relative humidity (RH) spanned over 20 years (1993-2012) were collected from NDRI and Central Soil and Salinity Research Institute, Karnal, India. Monthly average temperature humidity index (THI) values were estimated. Threshold THI value affecting fertility traits was identified by fixed least-squares model analysis. Three zones of non-heat stress, heat stress and critical heat stress zones were developed in a year. The genetic parameters heritability (h(2)) and repeatability (r) of each fertility trait were estimated. Genetic evaluation of Murrah buffaloes was performed in each zone with respect to their expected breeding values (EBV) for fertility traits. RESULTS: Effect of THI was found significant (p<0.001) on all fertility traits with threshold THI value identified as 75. Based on THI values, a year was classified into three zones: Non heat stress zone(THI 56.71-73.21), HSZ (THI 75.39-81.60) and critical HSZ (THI 80.27-81.60). The EBVfor SP, PR, CR were estimated as 138.57 days, 0.362 and 69.02% in non-HSZ while in HSZ EBV were found as 139.62 days, 0.358 and 68.81%, respectively. EBV for SP was increased to 140.92 days and for PR and CR, it was declined to 0.357 and 68.71% in critical HSZ. CONCLUSION: The negative effect of THI was observed on EBV of fertility traits under the non-HSZ and critical HSZ Thus, the influence of THI should be adjusted before estimating the breeding values for fertility traits in Murrah buffaloes.