Comprehensive selection signature analyses in dairy cattle exploiting purebred and crossbred genomic data

The main objective of the current research was to locate, annotate, and highlight specific areas of the bovine genome that are undergoing intense positive selection. Here, we are analyzing selection signatures in crossbred (Bos taurus X Bos indicus), taurine (Bos taurus), and indicine (Bos indicus) cattle breeds. Indicine cattle breeds found throughout India are known for their higher heat tolerance and disease resilience. More breeds and more methods can provide a better understanding of the selection signature. So, we have worked on nine distinct cattle breeds utilizing seven different summary statistics, which is a fairly extensive approach. In this study, we carried out a thorough genome-wide investigation of selection signatures using bovine 50K SNP data. We have included the genotyped data of two taurine, two crossbreds, and five indicine cattle breeds, for a total of 320 animals. During the 1950s, these indicine (cebuine) cattle breeds were exported with the aim of enhancing the resilience of taurine breeds in Western countries. For this study, we employed seven summary statistics, including intra-population, i.e., Tajima's D, CLR, iHS, and ROH and inter-population statistics, i.e., FST, XP-EHH, and Rsb. The NCBI database, PANTHER 17.0, and CattleQTL database were used for annotation after finding the important areas under selection. Some genes, including EPHA6, CTNNA2, NPFFR2, HS6ST3, NPR3, KCNIP4, LIPK, SDCBP, CYP7A1, NSMAF, UBXN2B, UGDH, UBE2K, and DAB1, were shown to be shared by three or more different approaches. Therefore, it gives evidence of the most intense selection in these areas. These genes are mostly linked to milk production and adaptability traits. This study also reveals selection regions that contain genes which are crucial to numerous biological functions, including those associated with milk production, coat color, glucose metabolism, oxidative stress response, immunity and circadian rhythms.

Identification of novel polymorphism in mammary-derived growth inhibitor gene of water buffalo and its expression analysis in the mammary gland

Mammary-derived growth inhibitor (MDGI), a member of the lipophilic family of fatty acid-binding proteins, plays an important role in the development, regulation, and differentiation of the mammary gland. The aim of the study was to identify polymorphism in the MDGI gene and its expression analysis in the mammary gland at various stages of lactation, in Indian buffalo. Nucleotide sequence analysis of MDGI gene in different breeds of riverine and swamp buffaloes revealed a total of 16 polymorphic sites and one Indel. Different transcription factor binding sites were predicted for buffalo MDGI gene promoter sequence, using online tools and in-silico analysis indicating that the SNPs in this region can impact the gene expression regulation. Phylogenetic analysis exhibited the MDGI of buffalo being closer to other ruminants like cattle, yak, sheep, and goats. Further, the expression analysis revealed that buffalo MDGI being highly expressed in well-developed mammary glands of lactating buffalo as compared to involution/non-lactating and before functional development to start the milk production stage in heifers. Stage-specific variation in expression levels signifies the important functional role of the MDGI gene in mammary gland development and milk production in buffalo, an important dairy species in Southeast Asia.

Landmarks in the history of selective sweeps

Half a century ago, a seminal article on the hitchhiking effect by Smith and Haigh inaugurated the concept of the selection signature. Selective sweeps are characterised by the rapid spread of an advantageous genetic variant through a population and hence play an important role in shaping evolution and research on genetic diversity. The process by which a beneficial allele arises and becomes fixed in a population, leading to a increase in the frequency of other linked alleles, is known as genetic hitchhiking or genetic draft. Kimura's neutral theory and hitchhiking theory are complementary, with Kimura's neutral evolution as the 'null model' and positive selection as the 'signal'. Both are widely accepted in evolution, especially with genomics enabling precise measurements. Significant advances in genomic technologies, such as next-generation sequencing, high-density SNP arrays and powerful bioinformatics tools, have made it possible to systematically investigate selection signatures in a variety of species. Although the history of selection signatures is relatively recent, progress has been made in the last two decades, owing to the increasing availability of large-scale genomic data and the development of computational methods. In this review, we embark on a journey through the history of research on selective sweeps, ranging from early theoretical work to recent empirical studies that utilise genomic data.

Uncovering genes underlying coat color variation in indigenous cattle breeds through genome-wide positive selection

The identification of candidate genes related to pigmentation and under selective sweep provides insights into the genetic basis of pigmentation and the evolutionary forces that have shaped this variation. The selective sweep events in the genes responsible for normal coat color in Indian cattle groups are still unknown. To find coat color genes displaying signs of selective sweeps in the indigenous cattle, we compiled a list of candidate genes previously investigated for their association with coat color and pigmentation. After that, we performed a genome-wide scan of positive selection signatures using the BovineSNP50K Bead Chip in 187 individuals of seven indigenous breeds. We applied a wide range of methods to find evidence of selection, such as Tajima's D, CLR, iHS, varLD, ROH, and FST. We found a total of sixteen genes under selective sweep, that were involved in coat color and pigmentation physiology. These genes are CRIM1 in Gir, MC1R in Sahiwal, MYO5A, PMEL and POMC in Tharparkar, TYRP1, ERBB2, and ASIP in Red Sindhi, MITF, LOC789175, PAX3 and TYR in Ongole, and IRF2, SDR165 and, KIT in Nelore, ADAMTS19 in Hariana. These genes are related to melanin synthesis, the biology of melanocytes and melanosomes, and the migration and survival of melanocytes during development.

Biomarker discovery and authentication of cold-slaughtered chicken through classical analytical procedures and mass spectrometry based proteomic approaches

1. This study evaluated the suitability of routine analytical procedures and used mass spectrometry-based proteomic approaches to distinguish meat from dead chicken/ cold-slaughtered birds (CS), electrically stunned and slaughtered birds, as per standard protocols (ES), and birds slaughtered according to halal guidelines (HS).2. Meat from CS birds had lower (P < 0.05) pH, water-holding capacity and higher (P < 0.05) lipid oxidation, haem iron content, residual blood and total viable counts relative to ES and HS meat indicating poor quality.3. The results demonstrated the presence of unique protein bands on SDS-PAGE only in CS meat that can be used for routine screening.4. Protein analysis using MALDI-TOF mass spectrometry identified haemoglobin subunit alpha-A and alpha-D; Adenylate kinase isoenzyme 1 as reliable and stable marker proteins for authentication of dead chicken meat under raw and cooked conditions and halal slaughtered chicken, respectively.5. The methods used may be employed by the food safety and regulatory agencies for regular screening of meat quality and to authenticate CS or HS chicken.

Genome-wide expression analysis reveals different heat shock responses in indigenous (Bos indicus) and crossbred (Bos indicus X Bos taurus) cattle

Environmental heat stress in dairy cattle leads to poor health, reduced milk production and decreased reproductive efficiency. Multiple genes interact and coordinate the response to overcome the impact of heat stress. The present study identified heat shock regulated genes in the peripheral blood mononuclear cells (PBMC). Genome-wide expression patterns for cellular stress response were compared between two genetically distinct groups of cattle viz., Hariana (B. indicus) and Vrindavani (B. indicus X B. taurus). In addition to major heat shock response genes, oxidative stress and immune response genes were also found to be affected by heat stress. Heat shock proteins such as HSPH1, HSPB8, FKB4, DNAJ4 and SERPINH1 were up-regulated at higher fold change in Vrindavani compared to Hariana cattle. The oxidative stress response genes (HMOX1, BNIP3, RHOB and VEGFA) and immune response genes (FSOB, GADD45B and JUN) were up-regulated in Vrindavani whereas the same were down-regulated in Hariana cattle. The enrichment analysis of dysregulated genes revealed the biological functions and signaling pathways that were affected by heat stress. Overall, these results show distinct cellular responses to heat stress in two different genetic groups of cattle. This also highlight the long-term adaptation of B. indicus (Hariana) to tropical climate as compared to the crossbred (Vrindavani) with mixed genetic makeup (B. indicus X B. taurus).

N6-Methyladenosine RNA Modification in Host Cells Regulates Peste des Petits Ruminants Virus Replication

N6-methyladenosine (m6A) modification is a major RNA epigenetic regulatory mechanism. The dynamics of m6A levels in viral genomic RNA and their mRNAs have been shown to have either pro- or antiviral functions, and therefore, m6A modifications influence virus-host interactions. Currently, no reports are available on the effect of m6A modifications in the genome of Peste des petits ruminants virus (PPRV). In the present study, we took PPRV as a model for nonsegmented negative-sense single-stranded RNA viruses and elucidate the role of m6A modification on viral replication. We detected m6A-modified sites in the mRNA of the virus and host cells, as well as the PPRV RNA genome. Further, it was found that the level of m6A modification in host cells alters the viral gene expression. Knockdown of the METTL3 and FTO genes (encoding the m6A RNA modification writer and eraser proteins, respectively) results in alterations of the levels of m6A RNA modifications in the host cells. Experiments using these genetically modified clones of host cells infected with PPRV revealed that both higher and lower m6A RNA modification in the host cells negatively affect PPRV replication. We found that m6A-modified viral transcripts had better stability and translation efficiency compared to the unmodified mRNA. Altogether, from these data, we conclude that the m6A modification of RNA regulates PPRV replication. These findings contribute toward a way forward for developing novel antiviral strategies against PPRV by modulating the dynamics of host m6A RNA modification. IMPORTANCE Peste des petits ruminants virus (PPRV) causes a severe disease in sheep and goats. PPRV infection is a major problem, causing significant economic losses to small ruminant farmers in regions of endemicity. N6-methyladenosine (m6A) is an important RNA modification involved in various functions, including virus-host interactions. In the present study, we used stable clones of Vero cells, having knocked down the genes encoding proteins involved in dynamic changes of the levels of m6A modification. We also used small-molecule compounds that interfere with m6A methylation. This resulted in a platform of host cells with various degrees of m6A RNA modification. The host cells with these different microenvironments were useful for studying the effect of m6A RNA modification on the expression of viral genes and viral replication. The results pinpoint the level of m6A modifications that facilitate the maximum replication of PPRV. These findings will be useful in increasing the virus titers in cultured cells needed for the economical development of the vaccine. Furthermore, the findings have guiding significance for the development of novel antiviral strategies for limiting PPRV replication in infected animals.

A single amino acid substitution in the VP2 protein of Indian foot-and-mouth disease virus serotype O vaccine strain confers thermostability and protective immunity in cattle

Foot-and-mouth disease (FMD) is a significant threat to animal health globally. Prophylactic vaccination using inactivated FMD virus (FMDV) antigen is being practised for the control in endemic countries. A major limitation of the current vaccine is its susceptibility to high environmental temperature causing loss of immunogenicity, thus necessitating the cold chain for maintenance of its efficacy. Hence, the FMD vaccine with thermostable virus particles will be highly useful in sustaining the integrity of whole virus particle (146S) during storage at 4°C. In this study, 12 recombinant mutants of Indian vaccine strain of FMDV serotype O (O/IND/R2/1975) were generated through reverse genetics approach and evaluated for thermostability. One of the mutant viruses, VP2_Y98F was more thermostable than other mutants and the parent FMDV. The oil-adjuvanted vaccine formulated with the inactivated VP2_Y98F mutant FMDV was stable up to 8 months when stored at 4°C and induced protective antibody response till dpv 180 after primary vaccination. It is concluded that the VP2_Y98F mutant FMDV was thermostable and has the potential to replace the parent vaccine strain.

Temporal dysregulation of genes in Lamb testis cell during sheeppox virus infection

The present study was aimed to elucidate the host-virus interactions using RNA-Seq analysis at 1h and 8h of post-infection of SPPV in LTC. The differentially expressed genes (DEGs) and the underlying mechanisms linked to the host immune responses were obtained. The protein-protein interaction (PPI) network analysis and Ingenuity pathway analysis (IPA) illustrated the interaction between the DEGs and their involvement in cell signalling responses. Highly connected hubs viz. AURKA, CHEK1, CCNB2, CDC6, and MAPK14 were identified through PPI network analysis. IPA analysis showed that IL-6 and ERK5 mediated signalling pathways were highly enriched at both time points. The TP53 gene was identified to be the leading upstream regulator that directly responded to SPPV infection, resulting in downregulation at both time points. The study provides an overview of how the lamb testis genes and their underlying mechanisms link to growth and immune response during SPPV infection.

Effect of acute heat shock on stress gene expression and DNA methylation in zebu (Bos indicus) and crossbred (Bos indicus × Bos taurus) dairy cattle

Environmental temperature is one of the major factors to affect health and productivity of dairy cattle. Gene expression networks within the cells and tissues coordinate stress response, metabolism, and milk production in dairy cattle. Epigenetic DNA methylations were found to mediate the effect of environment by regulating gene expression patterns. In the present study, we compared three Indian native zebu cattle, Bos indicus (Sahiwal, Tharparkar, and Hariana) and one crossbred Bos indicus × Bos taurus (Vrindavani) for stress gene expression and differences in the DNA methylation patterns. The results indicated acute heat shock to cultured PBMC affected their proliferation, stress gene expression, and DNA methylation. Interestingly, expressions of HSP70, HSP90, and STIP1 were found more pronounced in zebu cattle than the crossbred cattle. However, no significant changes were observed in global DNA methylation due to acute heat shock, even though variations were observed in the expression patterns of DNA methyltransferases (DNMT1, DNMT3a) and demethylases (TET1, TET2, and TET3) genes. The treatment 5-AzaC (5-azacitidine) that inhibit DNA methylation in proliferating PBMC caused significant increase in heat shock-induced HSP70 and STIP1 expression indicating that hypomethylation facilitated stress gene expression. Further targeted analysis DNA methylation in the promoter regions revealed no significant differences for HSP70, HSP90, and STIP1. However, there was a significant hypomethylation for BDNF in both zebu and crossbred cattle. Similarly, NR3C1 promoter region showed hypomethylation alone in crossbred cattle. Overall, the results indicated that tropically adapted zebu cattle had comparatively higher expression of stress genes than the crossbred cattle. Furthermore, DNA methylation may play a role in regulating expression of certain genes involved in stress response pathways.

Selection and validation of suitable reference gene for qPCR gene expression analysis in lamb testis cells under Sheep pox virus infection

Virus infection alters host gene expression, therefore ideal and stable reference housekeeping genes are required to normalise the expression of other expressed host genes in quantitative real-time PCR (qRT-PCR). The suitable reference gene may vary in response to different viral infections in different hosts or cells. In the present study, we cultured primary lamb testis cells (LTC) and assessed the expression stability of seven widely used housekeeping genes (B2M, HMBS, HPRT1, HSP-90, POLR2A, 18s_RNA, GAPDH) as reference genes in Sheeppox virus (SPPV) infected and control (uninfected-0h) LTC at 0.5h, 4.0h, 8.0h, and 12.0h post-infection) using NormFinder, Bestkeeper, geNorm, and the comparative ΔCT method in RefFinder based on their expression levels. Analysis revealed that HSP90, 18s_RNA, HPRT, POLR2A, and B2M were the most stable genes from the panel in the individual analysis group in 0h, 0.5h, 4.0h, 8.0h, and 12.0h, respectively. Furthermore, B2M was shown to be the most stable reference gene in the combined control with the respective and overall infected groups, except the control group of 4.0hpi of SPPV infection. In this study, we selected the most suitable reference genes in LTC for particular time points of SPPV infection. The identified most suitable housekeeping gene can be used during normalization of expression of other targeted genes at aspecific time point of SPPV infection.

Identification of important genomic footprints using eight different selection signature statistics in domestic cattle breeds

In the present study, the population genomic data of different cattle breeds were explored to decipher the genomic regions affected due to selective events and reflected in the productive, reproductive, thermo-tolerance, and health-related traits. To find out these genomic deviations due to selective sweeps, we used eight different statistical tools (Tajima's D, Fu & Li's D*, CLR, ROH, iHS, F_ST, FLK, and hapFLK) on seven indigenous and five exotic cattle breeds. We further performed composite analysis by comparing their covariance matrix. Several candidate genes were found to be related to milk production (ADARB, WDR70, and CA8), reproductive (PARN, FAM134B2, and ZBTB20), and health-related traits (SP110, CXCL2, CLXCL3, CXCL5, IRF8, and MYOM1). The outcome of this investigation provides a basis for detecting selective sweeps that explain the genetic variation of traits. They may possess functional importance for multiple cattle breeds in different subcontinents. However, further studies are required to improve the findings using high-density arrays or whole-genome sequencing with higher resolution and greater sample sizes.

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.

Characterisation and comparison of immune response mechanisms in an indigenous and a commercial pig breed after classical swine fever vaccination

The live attenuated classical swine fever (CSF) vaccine has been successfully used to prevent and control CSF outbreaks for 6 decades. However, the immune response mechanisms against the vaccine remain poorly understood. Moreover, very few reports exist regarding the breed differences in the response to CSF vaccine. In this study, we generated the peripheral blood mononuclear cell transcriptomes of indigenous Ghurrah and commercial Landrace pig breeds, before and 7 days after CSF vaccination. Subsequently, between and within-breed differential gene expression analyses were carried out. Results revealed large differences in pre-vaccination peripheral blood mononuclear cell transcriptome profiles of the two breeds, which were homogenised 7 days after vaccination. Before vaccination, gene set enrichment analysis showed that pathways related to antigen sensing and innate immune response were enriched in Ghurrah, while pathways related to adaptive immunity were enriched in Landrace. Ghurrah exhibited greater immunomodulation compared to Landrace following the vaccination. In Ghurrah, cell-cycle processes and T-cell response pathways were upregulated after vaccination. However, no pathways were upregulated in Landrace after vaccination. Pathways related to inflammation were downregulated in both the breeds after vaccination. Key regulators of inflammation such as IL1A, IL1B, NFKBIA and TNF genes were strongly downregulated in both the breeds after vaccination. Overall, our results have elucidated the mechanisms of host immune response against CSF vaccination in two distinct breeds and revealed common key genes instrumental in the global immune response to the vaccine.

Proteome Modulation in Peripheral Blood Mononuclear Cells of Peste des Petits Ruminants Vaccinated Goats and Sheep

In the present study, healthy goats and sheep (n = 5) that were confirmed negative for peste des petits ruminants virus (PPRV) antibodies by monoclonal antibody-based competitive ELISA and by serum neutralization test and for PPRV antigen by s-ELISA were vaccinated with Sungri/96. A quantitative study was carried out to compare the proteome of peripheral blood mononuclear cells (PBMCs) of vaccinated goat and sheep [5 days post-vaccination (dpv) and 14 dpv] vs. unvaccinated (0 day) to divulge the alteration in protein expression following vaccination. A total of 232 and 915 proteins were differentially expressed at 5 and 14 dpv, respectively, in goats. Similarly, 167 and 207 proteins were differentially expressed at 5 and 14 dpv, respectively, in sheep. Network generated by Ingenuity Pathway Analysis was “infectious diseases, antimicrobial response, and inflammatory response,” which includes the highest number of focus molecules. The bio functions, cell-mediated immune response, and humoral immune response were highly enriched in goats at 5 dpv and at 14 dpv. At the molecular level, the immune response produced by the PPRV vaccine virus in goats is effectively coordinated and stronger than that in sheep, though the vaccine provides protection from virulent virus challenge in both. The altered expression of certain PBMC proteins especially ISG15 and IRF7 induces marked changes in cellular signaling pathways to coordinate host immune responses.

Genome-wide detection of copy number variations in Tharparkar cattle

Copy number variations (CNVs) are major forms of genetic variation with an increasing importance in animal genomics. This study used the Illumina BovineSNP 50 K BeadChip to detect the genome-wide CNVs in the Tharparkar cattle. With the aid of PennCNV software, we noticed a total of 447 copy number variation regions (CNVRs) across the autosomal genome, occupying nearly 2.17% of the bovine genome. The average size of detected CNVRs was found to be 122.2 kb, the smallest CNVR being 50.02 kb in size, to the largest being 1,232.87 Kb. Enrichment analyses of the genes in these CNVRs gave significant associations with molecular adaptation-related Gene Ontology (GO) terms. Most CNVR genes were significantly enriched for specific biological functions; signaling pathways, sensory responses to stimuli, and various cellular processes. In addition, QTL analysis of CNVRs described them to be linked with economically essential traits in cattle. The findings here provide crucial information for constructing a more comprehensive CNVR map for the indigenous cattle genome.

Identification of genes affecting milk fat and fatty acid composition in Vrindavani crossbred cattle using 50 K SNP-Chip

The aim of this study was to identify candidate genes associated with milk fat per cent and fatty acid (FA) composition in Vrindavani cattle using the Illumina 50 K single-nucleotide polymorphism (SNP) array. After quality control, a total of 41,427 informative and high-quality SNPs were used for a genome-wide association study (GWAS) for milk fat percentage and 16 different types of fatty acids. Lactation stage, parity, test day milk yield, and proportion of exotic inheritance were included as fixed effects in the GWAS model. A total of 67 genome-wide significant (P < 1.20 × 10^-06) SNPs and 176 suggestive significant (P < 2.41 × 10^-05) SNPs were identified. Out of these, 15 SNPs were associated with more than one trait. The strongest associations were found on BTA14 for milk fat percentage and on BTA2 and BTA16 for polyunsaturated fatty acids. Several significant SNPs were identified close to or within the genes ELOVL6, FABP4, PMP2, PLIN1, MFGE8, GHRL2, and LDLRAD3 which are known to be associated with fat percentage and FA composition in dairy cattle breeds. This study is a step forward to better characterize the molecular mechanisms of phenotypic variation in milk fatty acids in a taurine-indicine composite cattle breed reared in tropical environments.

Genomic scans for selection signatures revealed candidate genes for adaptation and production traits in a variety of cattle breeds

Domestication and selection are the major driving forces responsible for the determinative genetic variability in livestock. These selection patterns create unique genetic signatures within the genome. BovineSNP50 chip data from 236 animals (seven indicine and five taurine cattle breeds) were analyzed in the present study. We implemented three complementary approaches viz. iHS (Integrated haplotype score), ROH (Runs of homozygosity), and F_ST, to detect selection signatures. A total of 179, 56, and 231 regions revealed 518, 277, and 267 candidate genes identified by iHS, ROH, and F_ST methods, respectively. We found several candidate genes (e.g., NCR3, ARID5A, HIST1H2BN, DEFB4, DEFB7, HSPA1L, HSPA1B, and DNAJB4) related to production traits and the adaptation of indigenous breeds to local environmental constraints such as heat stress and disease susceptibility. However, further studies are warranted to refine the findings using a larger sample size, whole-genome sequencing, and/or high density genotyping.

SNPs with intermediate minor allele frequencies facilitate accurate breed assignment of Indian Tharparkar cattle

Tharparkar cattle breed is widely known for its superior milch quality and hardiness attributes. This study aimed to develop an ultra-low density breed-specific single nucleotide polymorphism (SNP) genotype panel to accurately quantify Tharparkar populations in biological samples. In this study, we selected and genotyped 72 Tharparkar animals randomly from Cattle & Buffalo Farm of IVRI, India. This Bovine SNP50 BeadChip genotypic datum was merged with the online data from six indigenous cattle breeds and five taurine breeds. Here, we used a combination of pre-selection statistics and the MAF-LD method developed in our laboratory to analyze the genotypic data obtained from 317 individuals of 12 distinct breeds to identify breed-informative SNPs for the selection of Tharparkar cattle. This methodology identified 63 unique Tharparkar-specific SNPs near intermediate gene frequencies. We report several informative SNPs in genes/QTL regions affecting phenotypes or production traits that might differentiate the Tharparkar breed.

Signatures of Selection in Composite Vrindavani Cattle of India

Vrindavani is an Indian composite cattle breed developed by crossbreeding taurine dairy breeds with native indicine cattle. The constituent breeds were selected for higher milk production and adaptation to the tropical climate. However, the selection response for production and adaptation traits in the Vrindavani genome is not explored. In this study, we provide the first overview of the selection signatures in the Vrindavani genome. A total of 96 Vrindavani cattle were genotyped using the BovineSNP50 BeadChip and the SNP genotype data of its constituent breeds were collected from a public database. Within-breed selection signatures in Vrindavani were investigated using the integrated haplotype score (iHS). The Vrindavani breed was also compared to each of its parental breeds to discover between-population signatures of selection using two approaches, cross-population extended haplotype homozygosity (XP-EHH) and fixation index (F_ST). We identified 11 common regions detected by more than one method harboring genes such as LRP1B, TNNI3K, APOB, CACNA2D1, FAM110B, and SPATA17 associated with production and adaptation. Overall, our results suggested stronger selective pressure on regions responsible for adaptation compared to milk yield.

Selection and validation of differentially expressed metabolic and immune genes in weaned Ghurrah versus crossbred piglets

In the present investigation, differentially expressed genes (DEGs) were studied using RNA sequencing (RNA-seq) technique in porcine peripheral blood mononuclear cells (PBMC) of weaned Ghurrah and crossbred piglets at 3-month age. Transcriptomic analysis was done using three different packages, namely, EBSeq, DESeq2, and edgeR, to identify the DEGs between Ghurrah and crossbred piglets. Total 7717 DEGs were commonly identified by all three packages, out of which 4151 genes found to be up-regulated, and 3566 genes were down-regulated. Functional annotation of these DEGs indicated metabolism as the most commonly enriched category followed by the immune response. Genes related to metabolism and growth were up-regulated in crossbred piglets as compared with Ghurrah piglets, whereas immunity-related genes were up-regulated in Ghurrah piglets elucidating the disease resistance nature of this indigenous breed over crossbred counterparts. Further, eight DEGs, namely, LRP-1, ADCY4, ERRFI1, LDHD, ARG1, OASL, MGARP, and S100A8, were validated by qRT-PCR in a separate set of biological samples and found to be in concordance with RNA-seq results. Finding in the present study provides insight into genes and their molecular mechanisms governing difference in growth performance between Ghurrah and crossbred pigs.

Detection of recent infection of Japanese encephalitis virus in swine population using IgM ELISA: A suitable sentinel to predict infection in humans

Japanese encephalitis (JE) is a mosquito-borne flaviviral zoonotic disease and is one of the major causes of encephalitis in children. Swine, being an amplifier host of Japanese encephalitis virus (JEV), play an important role in its epidemiology. Therefore, early detection of either JEV or antibodies against JEV in swine is a feasible alternative for initiating necessary measures to prevent the spread of infection to humans. Since IgM antibodies appear early in swine sera, recombinant NS1 protein based indirect IgM ELISA was developed in the present study with the objective to know the recent infection of swine population with JEV. The relative diagnostic sensitivity and specificity of the developed ELISA was 95.34% and 98.6%, respectively. The developed ELISA was found to have excellent reproducibility on inter-laboratory and inter-institutional validation studies. A total of 3,027 field swine sera samples were screened using the developed ELISA and 488 samples were found positive for IgM against JEV with an overall sero-positivity of 16.12% in swine population of India. The highest sero-positivity was observed in swine population of Eastern zone of India which coincided with the maximum number of human JE cases reported from this zone during the same period. Further, antibody kinetics study revealed that the IgM antibodies against NS1 protein of JEV started appearing in swine sera at day 5 and disappeared completely by day 40. The IgG antibodies started appearing at day 7, and remained for more than 365 days indicating the suitability of IgM ELISA to know the recent infection of JEV. The developed IgM ELISA can be readily incorporated into surveillance programs for detection of JEV activity in swine population so that outbreaks in humans can be prevented by taking suitable preventive measures.

Genome-wide assessment of genetic diversity, linkage disequilibrium and haplotype block structure in Tharparkar cattle breed of India

Knowledge about genetic diversity is very essential for the management and sustainable utilization of livestock genetic resources. In this study, we presented a comprehensive genome-wide analysis of genetic diversity, ROH, inbreeding, linkage disequilibrium, effective population size and haplotype block structure in Tharparkar cattle of India. A total of 24 Tharparkar animals used in this study were genotyped with Illumina BovineSNP50 array. After quality control, 22,825 biallelic SNPs were retained, which were in HWE, MAF > 0.05 and genotyping rate >90%. The overall mean observed (H_O) and expected heterozygosity (H_E) were 0.339 ± 0.156 and 0.325 ± 0.129, respectively. The average minor allele frequency was 0.234 with a standard deviation of ± 0.131. We identified a total of 1832 ROH segments and the highest autosomal coverage of 13.87% was observed on chromosome 23. The genomic inbreeding coefficients estimates by F_ROH, F_HOM, F_GRM and F_UNI were 0.0589, 0.0215, 0.0532 and 0.0160 respectively. The overall mean linkage disequilibrium (LD) for a total of 133,532 pairwise SNPs measured by D' and r² was 0.6452 and 0.1339, respectively. In addition, we observed a gradual decline in effective population size over the past generations.

Foot-and-mouth disease virus induces PERK-mediated autophagy to suppress the antiviral interferon response

Foot-and-mouth disease virus (FMDV) is a picornavirus that causes contagious acute infection in cloven-hoofed animals. FMDV replication-associated viral protein expression induces endoplasmic reticulum (ER) stress and the unfolded protein response (UPR), in turn inducing autophagy to restore cellular homeostasis. We observed that inhibition of BiP (also known as HSPA5 and GRP78), a master regulator of ER stress and UPR, decreased FMDV infection confirming their involvement. Further, we show that the FMDV infection induces UPR mainly through the PKR-like ER kinase (PERK; also known as EIF2AK3)-mediated pathway. Knockdown of PERK and chemical inhibition of PERK activation resulted in decreased expression of FMDV proteins along with the reduction of autophagy marker protein LC3B-II [the lipidated form of LC3B (also known as MAP1LC3B)]. There are conflicting reports on the role of autophagy in FMDV multiplication. Our study systematically demonstrates that during FMDV infection, PERK-mediated UPR stimulated an increased level of endogenous LC3B-II and turnover of SQSTM1, thus confirming the activation of functional autophagy. Modulation of the UPR and autophagy by pharmacological and genetic approaches resulted in reduced numbers of viral progeny, by enhancing the antiviral interferon response. Taken together, this study underscores the prospect of exploring PERK-mediated autophagy as an antiviral target.

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.

Ancestry informative markers derived from discriminant analysis of principal components provide important insights into the composition of crossbred cattle

The cost of SNP genotyping to screen different breeds and to estimate the exact proportion of ancestry level is quite high, which can be compensated through deriving a small panel of ancestry informative markers (AIMs). Hence, we carried out the present study to provide an insight into ancestry level inferred from a panel of informative markers in the crossbred Vrindavani population developed at ICAR-IVRI, India. We have performed a new method i.e., discriminant analysis of principal components (DAPC) for the first time on the dataset of Vrindavani cattle. To confirm our method, we had performed DAPC on two other well-known crossbred cattle, i.e., Frieswal and Beefmaster. Three sets of panels (500, 1000 and 2000 markers) were tested for clustering of individuals. Among all the panels, we found the panel (1000 markers) with DAPC based contribution method was of the smallest size and comparatively of the highest accuracy.

Revelation of genomic breed composition in a crossbred cattle of India with the help of Bovine50K BeadChip

The aim of the present study was to assess the population structure and admixture levels in the Vrindavani composite population in India by using Bovine50KSNP BeadChip data. Genotypic data were generated for randomly selected animals (n = 72) of Vrindavani population and the data for parental breeds i.e., Hariana (n = 10), Holstein-Friesian (n = 63), Jersey (n = 28) and Brown Swiss (n = 22) were retrieved from a public repository. The indices of population structure were calculated using PLINK software and R-program. The merged dataset was analysed for assessing admixture levels and population stratification using three different approaches i.e., principal component analysis (PCA), multi-dimensional scaling (MDS) approach and the model-based approach in STRUCTURE software. The average minor allele frequency (MAF) value for Vrindavani population was estimated to be 0.235. Vrindavani population was found to possess an average ancestry of 39.5, 22.9, 26.9, and 10.7% inheritance levels from Holstein Friesian, Jersey, Hariana and Brown Swiss cattle breeds, respectively.

Comparative analysis of five different methods to design a breed-specific SNP panel for cattle

Single nucleotide polymorphisms (SNPs) have now replaced microsatellite markers in several species for various genetic investigations like parentage assignment, genetic breed composition, assessment for individuality and, most popularly, as a useful tool in genomic selection. However, such a resource, which can offer to assist breed identification in a cost-effective manner is still not explored in cattle breeding programs. In our study, we have tried to describe methods for reducing the number of SNPs to develop a breed-specific panel. We have used SNP data from Dryad open public access repository. Starting from a global dataset of 178 animals belonging to 10 different breeds, we selected five panels each comprising of similar number of SNPs using different methods i.e., Delta, Pairwise Wright's F_ST, informativeness for assignment, frequent item feature selection (FIFS) and minor allele frequency-linkage disequilibrium (MAF-LD) based method. MAF-LD based method has been recently developed by us for construction of breed-specific SNP panels. The STRUCTURE software analysis of MAF-LD based method showed appropriate clustering in comparison to other panels. Later, the panel of 591 breed-specific SNPs was called to their respective breeds using Venny 2.1.0 and UGent web tools software. Breed-specific SNPs were later annotated by using various Bioinformatics softwares.

Development of recombinant nonstructural 1 protein based indirect enzyme linked immunosorbent assay for sero-surveillance of Japanese encephalitis in swine

Japanese encephalitis virus (JEV) causes severe neurological disease in humans, especially among children. The disease is endemic in several South Asian countries including India. Swine play a major role as amplifier host for JEV and act as a source of infection to humans through mosquito bite. Early detection of either virus or antibodies in swine will aid to undertake control measures to prevent virus spread to humans. Swine seldom show symptoms of JEV infection and the viraemic phase lasts for a short period of 3 to 4 days indicating the potential of detection of antibodies, which remain for relatively longer period, as a suitable alternative. Cost effective and sensitive assays for the detection of JEV antibodies in swine are not available indigenously. Hence, we have developed a recombinant nonstructural protein 1 (rNS1) based enzyme linked immunosorbent assay for the detection of IgG antibodies against JEV in swine. The test is robust, highly sensitive (91%), specific (97%), reproducible and affordable. Field validation of the assay was done by screening 3628 swine Serum samples collected from different parts of India. The overall sero-positivity was found to be 32.22%. The developed ELISA can be readily incorporated into surveillance programs for detection of Japanese encephalitis virus activity in swine population thereby aiding in prediction of outbreaks in humans.

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 log₂ 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.

Early transcriptome profile of goat peripheral blood mononuclear cells (PBMCs) infected with peste des petits ruminant's vaccine virus (Sungri/96) revealed induction of antiviral response in an interferon independent manner

Sungri/96 vaccine strain is considered the most potent vaccine providing long-term immunity against peste des petits ruminants (PPR) in India. Previous studies in our laboratory highlighted induction of robust antiviral response in an interferon independent manner at 48 h and 120 h post infection (p.i.). However, immune response at the earliest time point 6 h p.i. (time taken to complete one PPRV life cycle), in PBMCs infected with Sungri/96 vaccine virus has not been investigated. This study was taken up to understand the global gene expression profiling of goat PBMCs after Sungri/96 PPRV vaccine strain infection at 6 h post infection (p.i.). A total of 1926 differentially expressed genes (DEGs) were identified with 616 - upregulated and 1310 - downregulated. TLR7/TLR3, IRF7/IRF1, ISG20, IFIT1/IFIT2, IFITM3, IL27 and TREX1 were identified as key immune sensors and antiviral candidate genes. Interestingly, type I interferons (IFNα/β) were not differentially expressed at this time point as well. TREX1, an exonuclease which inhibits type I interferons at the early stage of virus infection was found to be highly upregulated. IL27, an important antiviral host immune factor was significantly upregulated. ISG20, an antiviral interferon induced gene with exonuclease activity specific to ssRNA viruses was highly expressed. Functional profiling of DEGs showed significant enrichment of immune system processes with 233 genes indicating initiation of immune defense response in host cells. Protein interaction network showed important innate immune molecules in the immune network with high connectivity. The study highlights important immune and antiviral genes at the earliest time point.

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.

Multiplexed Autoantibody Signature for Serological Detection of Canine Mammary Tumours

Spontaneously occurring canine mammary tumours (CMTs) are the most common neoplasms of female unspayed dogs and are of potential importance as models for human breast cancer as well. Mortality rates are thrice higher in dogs as compared to humans with breast cancer, which can partly be attributed to lack of diagnostic techniques for their early detection. Human breast cancer studies reveal role of autoantibodies in early cancer diagnosis and also the usefulness of autoantibody panels in increasing the sensitivity, as well as, specificity of diagnostic assays. Therefore, in this study, we took advantage of high-throughput Luminex technique for developing a multiplex assay to detect autoantibody signatures against 5 canine mammary tumour-associated autoantigens (TAAs). These TAAs were expressed separately as fusion proteins with halo tag at the N-terminus, which allows easy and specific covalent coupling with magnetic microspheres. The multiplex assay, comprising a panel of candidate autoantigens (TPI, PGAM1, MNSOD, CMYC & MUC1) was used for screening circulating autoantibodies in 125 dog sera samples, including 75 mammary tumour sera and 50 healthy dog sera. The area under curve (AUC) of the combined panel of biomarkers is 0.931 (p < 0.0001), which validates the discriminative potential of the panel in differentiating tumour patients from healthy controls. The assay could be conducted in 3hrs using only 1ul of serum sample and could detect clinical cases of canine mammary tumour with sensitivity and specificity of 78.6% and 90%, respectively. In this study, we report for the first time a multiplexed assay for detection of autoantibodies in canine tumours, utilizing luminex technology and halo-tag coupling strategy. Further to the best of our knowledge, autoantibodies to CMYC and MUC1 have been reported for the first time in canines in this study.

In vitro study on role of σB protein in avian reovirus pathogenesis

Avian reoviruses, members of Orthoreovirus genus was known to cause diseases like tenosynovitis, runting-stunting syndrome in chickens. Among eight structural proteins, the proteins of S-class are mainly associated with viral arthritis but the significance of σB protein in arthritis is not established till date. In this infection pathological condition together with infection of joints often leads to arthritis because joints consists of cartilage which forms lubricating surface between two bones, and has limited metabolic, replicative and repair capacity. To establish the role of σB protein in arthritis, an in-vitro microarray study was conducted consisting four groups viz. virus infected and control; pDsRed-Express-N1-σB and empty pDs-Red transfected, CEF cells. With cut-off value as FC ≥2, p value <0.05, 6709 and 4026 numbers of DEGs in virus and σB, respectively were identified. The Ingenuity Pathway Analysis gave an idea about the involvement of σB protein in "osteoarthritis pathway", which was activated with z-score with 3.151. The pathway "Role of IL-17A in arthritis pathway" was also enriched with -log (p-value) 1.64. Among total 122 genes involved in osteoarthritis pathway, 28 upregulated and 11 downregulated DEGs were common to both virus and σB treated cells. Moreover, 14 upregulated and 7 downregulated were unique in σB transfected cells. Using qRT-PCR for IL-1B, BMP2, SMAD1, SPP1 genes, the microarray data was validated. We concluded that during ARV infection σB protein, if not fully partially leads to molecular alteration of various genes of host orchestrating the different molecular pattern in joints, leading to tenosynovitis syndrome.

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.

PCR-SSCP analysis of MDGI gene and its association with milk production traits in river buffalo (Bubalus bubalis)

In this study, we investigated the genetic variation within 3'UTR of Mammary-Derived Growth Inhibitor (MDGI) gene of buffalo using PCR-SSCP and sequencing; and also analyzed association of polymorphism with the milk production traits. The study revealed two conformational patterns, 'A' and 'B' among 234 Mehsana buffaloes maintained with their records in the field and at farm. The frequency of SSCP variant 'A' was found to be invariably high in the buffalo population under study. Further, association analysis of SSCP variants with various milk production and milk quality traits indicated no significant effect on any of the traits investigated. Sequencing of SSCP variant 'A' showed homozygous G/G and A/A and 'B' had heterozygous G/C and A/G at positions +124 and +140 respectively, in the 3'UTR of buffalo MDGI. The preliminary results showed the substantial variations in the distribution of SSCP variants' frequencies within Mehsana buffaloes, however these variants had non-significant association with milk yield, fat yield and fat percentage in Mehsana buffaloes.

Hemato-biochemical alterations and urinalysis in dogs suffering from benign prostatic hyperplasia

Aim:

The study was designed to evaluate the hemato-biochemical alterations, urinalysis along with histomorphological and histological changes of prostate glands in dogs affected with benign prostatic hyperplasia (BPH) in and around Bhubaneswar, Odisha, India.

Materials and Methods:

In toto, 445 dogs presented to the Teaching Veterinary Clinical Complex of the College of Veterinary Sciences and Animal Husbandry, one Government Veterinary Hospital and two pet clinics in and around Bhubaneswar screened for the presence of BPH. Most of the 57 dogs were 6 years and above as reported by the owners. Only 57 dogs found positive for BPH basing on the presence of typical clinical signs subjected for a detailed hemato-biochemical study. Most of the 57 dogs were 6 years and above as reported by the owners. Routine and microscopic urinalyses were done as per the routine procedure. Histomorphological evaluations of prostate glands were done through manual rectal palpation. Histological examinations of prostate tissue sections of two dead dogs were conducted with routine hematoxylin and eosin stain.

Results:

The study revealed about 12.8% (57/445) of dogs was suffering from BPH. Typical clinical signs - such as passing small thin tape-shaped feces, holding tail away from backward, tenesmus, and straining during urination and defecation - were seen in most of the cases. Urine samples of affected dogs were positive for glucose, occult blood, and protein. A significant decrease in lymphocytes and increase in eosinophil counts in dogs with BPH was recorded. Serum biochemical analysis showed a nonsignificant increase in creatinine and blood urea nitrogen with a significant decrease in total protein, albumin, globulin, A:G ratio. Histology of prostate glands collected during postmortem was characterized by fibrosis of prostate gland, and hyperplasia of the acinar epithelium.

Conclusions:

High rate of the prevalence of BPH in dogs poses an alarming condition which if diagnosed at an early stage can certainly prolong the longevity of the dogs.

Identification of genetic variation in NOD-like receptor 2 gene and influence of polymorphism on gene structure and function in buffalo (Bubalus bubalis)

Nucleotide-binding oligomerization domain (NOD)-like receptor 2 is one of the important mediators of innate as well as adaptive immune response to microbial infections. In this study, NOD-like receptor-2 was characterized by determining the full gene sequence and analyzing genetic diversity in Indian buffaloes. Sequence analysis of buffalo NOD2 revealed 3042 nucleotides long ORF, encoding 1013 amino acids from 12 exons. Domain structure analysis indicated existence of 8 leucine-rich repeat (LRR) domains in buffalo, cattle, sheep and mouse, along with central NACHT/NOD domain and two N-terminal CARD domains. Comparative sequence analysis among different buffalo breeds identified 46 polymorphic sites in NOD2 gene. Among coding region SNPs, 10 were non-synonymous, 7 synonymous and 3 were present in 5'UTR. Genotyping of two nsSNPs, revealed significant differences in the allele frequencies, distinguishing swamp and riverine buffaloes, having different utilities. Association analysis with mastitis in dairy buffaloes indicated significant variation in allelic frequencies at G1135A locus, between mastitis affected and non-affected animals. Further, NOD2 gene expression was quantified in different riverine buffalo tissues, using real-time PCR and lymph node displayed highest expression, compared to others organs included in the study. Overall, the study revealed buffalo NOD2 gene attributes, important to understand species specific immune response in ruminants.

Expression Analysis of Genes Associated with Prolificacy in FecB Carrier and Noncarrier Indian Sheep

The effect of FecB mutation on the gene expression in FecB carrier and noncarrier estrous synchronized ewes, has been analyzed. For this study the whole ovarian tissues and Graafian follicles were collected from estrus synchronized FecB carrier Garole, and non-carrier Deccani Indian sheep, showing remarkable differences in the numbers of preovulatory follicles among two groups. Eleven potential candidate genes (BMP15, GDF9, BMP4, BMP7, BMPR1B, BMPR1A, SMAD9, LHCGR, FSHR, IGF1R, and STAT5) were selected for their expression analysis by SybrGreen based real-time PCR, across ovaries and Graafian follicles of different fecundity groups, for having better insights into the effect of FecB genotypes on follicular development. Variable expression was observed for almost all the genes included in the present study among high and low fecundity groups that was most significant for the BMP7, BMP4, LHCGR, and FSHR transcripts in the ovarian follicles of high and low fecundity ewes, indicating their importance in governing the fecundity in FecB carrier, Indian Garole sheep. BMP4 expression among the genes studied was significantly higher in FecB carrier Garole sheep. This study confirms the changes in mRNA expression of the genes implicated in follicular development in FecB carrier and noncarrier Indian sheep breeds.

EvaGreen-based Multiplex Real-time PCR Assay for Rapid Differentiation of Wild-Type and Glycoprotein E-Deleted Bovine Herpesvirus-1 Strains

Bovine herpesvirus-1 (BoHV-1) is an important viral pathogen causing significant economic losses to the cattle industry. Glycoprotein E-deleted marker vaccines form the basis for BoHV-1 control programs widely, wherein detection and differentiation of wild-type and gE-deleted vaccine strains is of crucial importance for proper disease management. In the present study, we report an EvaGreen-based multiplex real-time polymerase chain reaction (EGRT-PCR) assay for rapid differentiation of wild-type and glycoprotein E-deleted strains of BoHV-1. The EGRT-PCR assay could simultaneously detect two viral genes (glycoprotein B and E) and an internal positive control gene (bovine growth hormone- bGH), in a single-tube reaction. The analytical sensitivity of the EGRT-PCR assay was as little as 10 copies of the BoHV-1 DNA per reaction. The modified real-time PCR assay could successfully differentiate wild-type and gE-deleted BoHV-1 strains based on gene specific melting temperatures (T_m) peaks. Our results have shown that the EGRT-PCR developed in this study might prove to be a promising tool in disease management by enabling rapid differentiation of wild-type and gE-deleted strains of BoHV-1.

A double-stranded probe coupled with isothermal amplification for qualitative and quantitative detection of avian reovirus

We applied a probe-based real-time loop-mediated isothermal amplification (Cy5-RTqLAMP) technique targeting the avian reovirus (ARV) S3 gene to develop a rapid, sensitive, and specific method for virus detection and quantification. This test specifically detected the presence of ARV, but not other viruses or bacteria present in clinical or artificially spiked samples, including Newcastle disease virus, infectious bursal disease virus, fowl adenovirus, Marek's disease virus, Escherichia coli, and Salmonella spp. This test can detect ARV in less than one hour with an analytical sensitivity of 10 viral gene copies and 1 fg of total cDNA. The Cy5-RTqLAMP does not yield false positive results and is 100 times more sensitive than conventional PCR. This test was shown to be able to detect the presence of ARV in clinical samples. A similar strategy may be used for detection of other important human and animal viral pathogens.

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.

•

Deciphering the host virus interactions through system biology approach reviewed

•

High throughput techniques to understand the host pathogen interactions examined

•

Shift from virus-centric perspective to spectrum of virus-host interactions

•

Modeling of host-virus cross talk

MARKER ASSISTED EVALUATION OF MORPHOLOGICAL AND GENETIC ATTRIBUTES OF SUB-POPULATIONS OF NILI-RAVI BUFFALO: A VULNERABLE DAIRY TYPE RIVERINE BREED OF INDIA

In the present study, we report the distribution of true to type and atypical Nili-Ravi buffalo, a vulnerable dairy type riverine breed of North India and its underlying genetic structure. Out of total investigated buffaloes 73.5% had bilateral wall eyes while 5.4% had unilateral wall eyes and 21.1% had no wall eyes. 41.15% of Nili-Ravi buffaloes maintained in the breeding farm were having typical true to the type characteristics (both eyes walled, white markings in forehead, muzzle/chin, all the four legs and tail) while only 28.5% of Nili-Ravi buffaloes were true to the type under field conditions. Genotypic data were generated in four groups of Nili-Ravi buffalo (FMTNR--Typical Nili-Ravi from farm; FMANR--Atypical Nili-Ravi from farm; FDTNR--Typical Nili-Ravi from field; FDANR--Atypical Nili-Ravi from field) at 16 microsatellite loci. Comparative genetic analysis of various groups of Nili-Ravi buffaloes with Murrah revealed significant between group differences with an estimated global F(ST) of 0.063. Pair-wise F(ST) values ranged from 0.003 (between FDTNR and FDANR) to 0.112 (between FMTNR and FDTNR). Phylogenetic analysis and multi-dimensional scaling revealed clustering of FDTNR and FDANR together while FMTNR and FMANR clustered separately with Murrah in between farm and field Nili-Ravi buffaloes. Based on the results, the paper also proposes three pronged strategy for conservation and sustainable genetic improvement of Nili-Ravi buffalo in India.

Genetic analysis of river, swamp and hybrid buffaloes of north-east India throw new light on phylogeography of water buffalo (Bubalus bubalis)

This study analysed buffaloes from north-east India and compared their nuclear and mitochondrial DNA variations with buffaloes of mainland India, China, Mediterranean and South-East Asia. Microsatellite genotypes of 338 buffaloes including 210 from six north-east Indian buffalo populations and three mainland Indian breeds were analysed to evaluate their genetic structure and evolutionary relationships. Phylogenetic analysis and multidimensional scaling plot of pairwise FST revealed the clustering of all swamp-type buffaloes of north-east India with Lower Assamese (significantly hybrid type) buffaloes in one plane and all the mainland river buffaloes in another plane while the upper Assamese buffaloes being distinct from both these clusters. Analysis of mtDNA D-loop region of 530-bp length was performed on 345 sequences belonging to 23 buffalo populations from various geographical regions to establish the phylogeography of Indian water buffalo. The swamp buffaloes of north-east India clustered with both the lineages of Chinese swamp buffalo. Multidimensional scaling display of pairwise FST derived from mitochondrial DNA data showed clustering of upper Assamese, Chilika and Mediterranean buffaloes distinctly from all the other Indian buffalo populations. Median-joining network analysis further confirmed the distinctness and ancestral nature of these buffaloes. The study revealed north-east region of India forming part of the wider hybrid zone of water buffalo that may probably extend from north-east India to South-East Asia.

Genomic analysis of host - Peste des petits ruminants vaccine viral transcriptome uncovers transcription factors modulating immune regulatory pathways

Peste des petits ruminants (PPR), is an acute transboundary viral disease of economic importance, affecting goats and sheep. Mass vaccination programs around the world resulted in the decline of PPR outbreaks. Sungri 96 is a live attenuated vaccine, widely used in Northern India against PPR. This vaccine virus, isolated from goat works efficiently both in sheep and goat. Global gene expression changes under PPR vaccine virus infection are not yet well defined. Therefore, in this study we investigated the host-vaccine virus interactions by infecting the peripheral blood mononuclear cells isolated from goat with PPRV (Sungri 96 vaccine virus), to quantify the global changes in the transcriptomic signature by RNA-sequencing. Viral genome of Sungri 96 vaccine virus was assembled from the PPRV infected transcriptome confirming the infection and demonstrating the feasibility of building a complete non-host genome from the blood transcriptome. Comparison of infected transcriptome with control transcriptome revealed 985 differentially expressed genes. Functional analysis showed enrichment of immune regulatory pathways under PPRV infection. Key genes involved in immune system regulation, spliceosomal and apoptotic pathways were identified to be dysregulated. Network analysis revealed that the protein - protein interaction network among differentially expressed genes is significantly disrupted in infected state. Several genes encoding TFs that govern immune regulatory pathways were identified to co-regulate the differentially expressed genes. These data provide insights into the host - PPRV vaccine virus interactome for the first time. Our findings suggested dysregulation of immune regulatory pathways and genes encoding Transcription Factors (TFs) that govern these pathways in response to viral infection.