Middle East respiratory syndrome coronavirus (MERS-CoV) has evolved to navigate through the sophisticated network of a host's immune system. The immune evasion mechanism including type 1 interferon and protein kinase R-mediated antiviral stress responses has been recently attributed to the involvement of MERS-CoV protein 4a (p4a) that masks the viral dsRNA. However, the structural mechanism of how p4a recognizes and establishes contacts with dsRNA is not well explained. In this study, we report a dynamic mechanism deployed by p4a to engage the viral dsRNA and make it unavailable to the host immune system. Multiple variants of p4a-dsRNA were created and investigated through extensive molecular dynamics procedures to highlight crucial interfacial residues that may be used as potential pharmacophores for future drug development. The structural analysis revealed that p4a exhibits a typical αβββα fold structure, as found in other dsRNA-binding proteins. The α1 helix and the β1-β2 loop play a crucial role in recognizing and establishing contacts with the minor grooves of dsRNA. Further, mutational and binding free energy analyses suggested that in addition to K63 and K67, two other residues, K27 and W45, might also be crucial for p4a-dsRNA stability.

β-defensins are adsorbable on the sperm surface in the male reproductive tract (MRT) and enhance sperm functional characteristics. The beta-defensin 129 (DEFB129) antimicrobial peptide is involved in sperm maturation, motility, and fertilization. However, its role in bovine fertility has not been well investigated. This study examines the relationship between the bovine BBD129 gene and Bos indicus x Bos taurus bull fertility. The complete coding sequence of BBD129 mRNA was identified by RNA Ligase Mediated-Rapid Amplification of cDNA End (RLM-RACE) and Sanger sequencing methodologies. It consisted of 582 nucleotides (nts) including 5' untranslated region (UTR) (46nts) and 3'UTR (23nts). It conserves all beta-defensin-like features. The expression level of BBD129 was checked by RT-qPCR and maximal expression was detected in the corpus-epididymis region compared to other parts of MRT. Polymorphism in BBD129 was also confirmed by Sanger sequencing of 254 clones from 5 high fertile (HF) and 6 low fertile (LF) bulls at two positions, 169 T > G and 329A > G, which change the S57A and N110S in the protein sequence respectively. These two mutations give rise to four types of BBD129 haplotypes. The non-mutated TA-BBD129 (169 T/329A) haplotype was substantially more prevalent among high-fertile bulls (P < 0.005), while the double-site mutated GG-BBD129 (169 T > G/329A > G) haplotype was significantly more prevalent among low-fertile bulls (P < 0.005). The in silico analysis confirmed that the polymorphism in BBD129 results in changes in mRNA secondary structure, protein conformations, protein stability, extracellular-surface availability, post-translational modifications (O-glycosylation and phosphorylation), and affects antibacterial and immunomodulatory capabilities. In conclusion, the mRNA expression of BBD129 in the MRT indicates its region-specific dynamics in sperm maturation. BBD129 polymorphisms were identified as the deciding elements accountable for the changed proteins with impaired functionality, contributing to cross-bred bulls' poor fertility.

The microRNAs (miRs) secreted by the trophectoderm (TE) cells have recently been implicated in the conceptus-endometrial cross talk during implantation and placentation. These miRs modulate various cellular processes during conception and throughout the pregnancy by regulating the gene expression in the foetal and maternal tissues.

Direct inhibition of tumor necrosis factor-alpha (TNF-α) action is considered a promising way to prevent or treat TNF-α-associated diseases. The trimeric form of TNF-α binds to its receptor (TNFR) and activates the downstream signaling pathway. The interaction of TNF-α with molecular-grade dimethyl sulfoxide (DMSO) in an equal volumetric ratio renders TNF-α inert, in this state, TNF-α fails to activate TNFR. Here, we aimed to examine the inhibition of TNF-α function by various concentrations of DMSO. Its higher concentration led to stronger attenuation of TNF-α-induced cytokine secretion by fibroblasts, and of their death. We found that this inhibition was mediated by a perturbation in the formation of the functional TNF-α trimer. Molecular dynamics simulations revealed a transient interaction between DMSO molecules and the central hydrophobic cavity of the TNF-α homodimer, indicating that a brief interaction of DMSO with the TNF-α homodimer may disrupt the formation of the functional homotrimer. We also found that the sensitizing effect of actinomycin D on TNF-α-induced cell death depends upon the timing of these treatments and on the cell type. This study will help to select an appropriate concentration of DMSO as a working solvent for the screening of water-insoluble TNF-α inhibitors.

Low conception rate (CR) despite insemination with morphologically normal spermatozoa is a common reproductive restraint that limits buffalo productivity. This accounts for a significant loss to the farmers and the dairy industry, especially in agriculture-based economies. The immune-related proteins on the sperm surface are known to regulate fertility by assisting the spermatozoa in their survival and performance in the female reproductive tract (FRT). Regardless of their importance, very few studies have specifically catalogued the buffalo sperm surface proteome. The study was designed to determine the identity of sperm surface proteins and to ascertain if the epididymal expressed beta-defensins (BDs), implicated in male fertility, are translated and applied onto buffalo sperm surface along with other immune-related proteins.

Agrichemicals such as organophosphorus pesticides' metabolites (OPPMs) are more hazardous and pervasive than their parent pesticides. Parental germline exposure to such xenobiotics leads to an elevated susceptibility towards reproductive failures e.g. sub- or in-fertility. This study sought to examine the effects of low-dose, acute OPPM exposure on mammalian sperm function using buffalo as the model organism. The buffalo spermatozoa were briefly (2 h) exposed to metabolites of the three most prevalent organophosphorus pesticides (OPPs) viz. Omethoate (from Dimethoate), paraoxon-methyl (from methyl/ethyl parathion) and 3, 5, 6-trichloro-2-pyridinol (from chlorpyrifos). Exposure to OPPMs resulted in compromised structural and functional integrity (dose-dependent) of the buffalo spermatozoa typified by elevated membrane damage, increased lipid peroxidation, precocious capacitation and tyrosine phosphorylation, perturbed mitochondrial activity and function and (P < 0.05). This led to a decline in the in vitro fertilizing ability (P < 0.01) of the exposed spermatozoa, as indicated by reduced cleavage and blastocyst formation rates. Preliminary data indicate that acute exposure to OPPMs, akin to their parent pesticides, induces biomolecular and physiological changes in spermatozoa that compromise their health and function ultimately affecting their fertility. This is the first study demonstrating the in vitro spermatotoxic effects of multiple OPPMs on male gamete functional integrity.

One of the most evolutionary conserved communication systems, the Wnt signaling pathway is a major gene regulatory pathway that affects the developmental competence of oocytes and regulates most embryonic developmental processes. The present study was undertaken to modulate the canonical Wnt (Wingless/integration) signaling pathway in the poor-quality (colorless cytoplasm after Brilliant Cresyl Blue staining, BCB-) buffalo cumulus-oocyte complexes (COCs) to improve their in vitro maturation (IVM) and embryo production (IVEP) rates.

Nucleotide-binding oligomerization domain-containing protein 1 (NOD1) and NOD2 are cytosolic pattern recognition receptors playing pivotal roles in innate immune signaling. NOD1 and NOD2 recognize bacterial peptidoglycan derivatives iE-DAP and MDP, respectively and undergoes conformational alternation and ATP-dependent self-oligomerization of NACHT domain followed by downstream signaling. Lack of structural adequacy of NACHT domain confines our understanding about the NOD-mediated signaling mechanism. Here, we predicted the structure of NACHT domain of both NOD1 and NOD2 from model organism zebrafish (Danio rerio) using computational methods. Our study highlighted the differential ATP binding modes in NOD1 and NOD2. In NOD1, γ-phosphate of ATP faced toward the central nucleotide binding cavity like NLRC4, whereas in NOD2 the cavity was occupied by adenine moiety. The conserved 'Lysine' at Walker A formed hydrogen bonds (H-bonds) and Aspartic acid (Walker B) formed electrostatic interaction with ATP. At Sensor 1, Arg328 of NOD1 exhibited an H-bond with ATP, whereas corresponding Arg404 of NOD2 did not. 'Proline' of GxP motif (Pro386 of NOD1 and Pro464 of NOD2) interacted with adenine moiety and His511 at Sensor 2 of NOD1 interacted with γ-phosphate group of ATP. In contrast, His579 of NOD2 interacted with the adenine moiety having a relatively inverted orientation. Our findings are well supplemented with the molecular interaction of ATP with NLRC4, and consistent with mutagenesis data reported for human, which indicates evolutionary shared NOD signaling mechanism. Together, this study provides novel insights into ATP binding mechanism, and highlights the differential ATP binding modes in zebrafish NOD1 and NOD2.

Nucleotide binding and oligomerization domain (NOD)-like receptors (NLRs) are innate immune receptors that recognize bacterial cell wall components and initiate host immune response. Structure and function of NLRs have been well studied in human and mice, but little information exists on genetic composition and role of these receptors in innate immune system of water buffalo--a species known for its exceptional disease resistance. Here, a comparative study on the functional domains of NOD1 and NOD2 was performed across different species. The NOD mediated in-vitro cellular responses were studied in buffalo peripheral blood mononuclear cells, resident macrophages, mammary epithelial, and fibroblast cells. Buffalo NOD1 (buNOD1) and buNOD2 showed conserved domain architectures as found in other mammals. The domains of buNOD1 and buNOD2 showed analogy in secondary and tertiary conformations. Constitutive expressions of NODs were ubiquitous in different tissues. Following treatment with NOD agonists, peripheral lymphocytes showed an IFN-γ response along-with production of pro-inflammatory cytokines. Alveolar macrophages and mammary epithelial cells showed NOD mediated in-vitro immune response through NF-κB dependent pathway. Fibroblasts showed pro-inflammatory cytokine response following agonist treatment. Our study demonstrates that both immune and non-immune cells could generate NOD-mediated responses to pathogens though the type and magnitude of response depend on the cell types. The structural basis of ligand recognition by buffalo NODs and knowledge of immune response by different cell types could be useful for development of non-infective innate immune modulators and next generation anti-inflammatory compounds.

Toll-like receptors (TLRs) play a fundamental role in the inflammatory response against invading pathogens. However, the dysregulation of TLR-signaling pathways is implicated in several autoimmune/inflammatory diseases. Here, we show that a novel small molecule TLR-inhibitor (TAC5) and its derivatives TAC5-a, TAC5-c, TAC5-d, and TAC5-e predominantly antagonized poly(I:C) (TLR3)-, imiquimod (TLR7)-, TL8-506 (TLR8)-, and CpG-oligodeoxynucleotide (TLR9)-induced signaling pathways. TAC5 and TAC5-a significantly hindered the activation of nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), reduced the phosphorylation of mitogen-activated protein kinases, and inhibited the secretion of tumor necrosis factor-α (TNF-α) and interleukin-6. Besides, TAC5-a prevented the progression of psoriasis and systemic lupus erythematosus (SLE) in mice. Interestingly, TAC5 and TAC5-a did not affect Pam3CSK4 (TLR1/2)-, FSL-1 (TLR2/6)-, or lipopolysaccharide (TLR4)-induced TNF-α secretion, indicating their specificity towards endosomal TLRs (TLR3/7/8/9). Collectively, our data suggest that the TAC5 series of compounds are potential candidates for treating autoimmune diseases such as psoriasis or SLE.

Water buffalo (Bubalus bubalis) are an important animal resource that contributes milk, meat, leather, dairy products, and power for plowing and transport. However, mastitis, a bacterial disease affecting milk production and reproduction efficiency, is most prevalent in populations having intensive selection for higher milk yield, especially where the inbreeding level is also high. Climate change and poor hygiene management practices further complicate the issue. The management of this disease faces major challenges, like antibiotic resistance, maximum residue level, horizontal gene transfer, and limited success in resistance breeding. Bovine mastitis genome wide association studies have had limited success due to breed differences, sample sizes, and minor allele frequency, lowering the power to detect the diseases associated with SNPs. In this work, we focused on the application of targeted gene panels (TGPs) in screening for candidate gene association analysis, and how this approach overcomes the limitation of genome wide association studies. This work will facilitate the targeted sequencing of buffalo genomic regions with high depth coverage required to mine the extremely rare variants potentially associated with buffalo mastitis. Although the whole genome assembly of water buffalo is available, neither mastitis genes are predicted nor TGP in the form of web-genomic resources are available for future variant mining and association studies. Out of the 129 mastitis associated genes of cattle, 101 were completely mapped on the buffalo genome to make TGP. This further helped in identifying rare variants in water buffalo. Eighty-five genes were validated in the buffalo gene expression atlas, with the RNA-Seq data of 50 tissues. The functions of 97 genes were predicted, revealing 225 pathways. The mastitis proteins were used for protein-protein interaction network analysis to obtain additional cross-talking proteins. A total of 1,306 SNPs and 152 indels were identified from 101 genes. Water Buffalo-MSTdb was developed with 3-tier architecture to retrieve mastitis associated genes having genomic coordinates with chromosomal details for TGP sequencing for mining of minor alleles for further association studies. Lastly, a web-genomic resource was made available to mine variants of targeted gene panels in buffalo for mastitis resistance breeding in an endeavor to ensure improved productivity and the reproductive efficiency of water buffalo.

Crossbred bulls produced by crossing Bos taurus and Bos indicus suffer with high incidence of infertility/subfertility problems; however, the etiology remains poorly understood. The uncertain predictability and the inability of semen evaluation techniques to maintain constant correlation with fertility demand for alternate methods for bull fertility prediction. Therefore, in this study, the global differential gene expression between high- and low-fertile crossbred bull sperm was assessed using a high-throughput RNA sequencing technique with the aim to identify transcripts associated with crossbred bull fertility. Crossbred bull sperm contained transcripts for 13,563 genes, in which 2,093 were unique to high-fertile and 5,454 were unique to low-fertile bulls. After normalization of data, a total of 776 transcripts were detected, in which 84 and 168 transcripts were unique to high-fertile and low-fertile bulls, respectively. A total of 176 transcripts were upregulated (fold change > 1) and 209 were downregulated (<1) in low-fertile bulls. Gene ontology analysis identified that the sperm transcripts involved in the oxidative phosphorylation pathway and biological process such as multicellular organism development, spermatogenesis, and in utero embryonic development were downregulated in low-fertile crossbred bull sperm. Sperm transcripts upregulated and unique to low-fertile bulls were majorly involved in translation (biological process) and ribosomal pathway. With the use of RT-qPCR, selected sperm transcripts (n = 12) were validated in crossbred bulls (n = 12) with different fertility ratings and found that the transcriptional abundance of ZNF706, CRISP2, TNP2, and TNP1 genes was significantly (p < 0.05) lower in low-fertile bulls than high-fertile bulls and was positively (p < 0.05) correlated with conception rate. It is inferred that impaired oxidative phosphorylation could be the predominant reason for low fertility in crossbred bulls and that transcriptional abundance of ZNF706, CRISP2, TNP2, and TNP1 genes could serve as potential biomarkers for fertility in crossbred bulls.

Sperm harbours a wide range of proteins regulating their functions and fertility. In the present study, we made an effort to characterize and quantify the proteome of buffalo bull spermatozoa, and to identify fertility associated sperm proteins through comparative proteomics. Using high-throughput mass spectrometry platform, we identified 1305 proteins from buffalo spermatozoa and found that these proteins were mostly enriched in glycolytic process, mitochondrial respiratory chain, tricarboxylic acid cycle, protein folding, spermatogenesis, sperm motility and sperm binding to zona pellucida (p < 7.74E-08) besides metabolic (p = 4.42E-31) and reactive oxygen species (p = 1.81E-30) pathways. Differential proteomic analysis revealed that 844 proteins were commonly expressed in spermatozoa from both the groups while 77 and 52 proteins were exclusively expressed in high- and low-fertile bulls, respectively. In low-fertile bulls, 75 proteins were significantly (p < 0.05) upregulated and 176 proteins were significantly (p < 0.05) downregulated; these proteins were highly enriched in mitochondrial respiratory chain complex I assembly (p = 2.63E-07) and flagellated sperm motility (p = 7.02E-05) processes besides oxidative phosphorylation pathway (p = 6.61E-15). The down regulated proteins in low-fertile bulls were involved in sperm motility, metabolism, sperm-egg recognition and fertilization. These variations in the sperm proteome could be used as potential markers for the selection of buffalo bulls for fertility.

Cathelicidins are an ancient class of antimicrobial peptides (AMPs) with broad spectrum bactericidal activities. In this study, we investigated the diversity and biological activity of cathelicidins of buffalo, a species known for its disease resistance. A series of new homologs of cathelicidin4 (CATHL4), which were structurally diverse in their antimicrobial domain, was identified in buffalo. AMPs of newly identified buffalo CATHL4s (buCATHL4s) displayed potent antimicrobial activity against selected Gram positive (G+) and Gram negative (G-) bacteria. These peptides were prompt to disrupt the membrane integrity of bacteria and induced specific changes such as blebing, budding, and pore like structure formation on bacterial membrane. The peptides assumed different secondary structure conformations in aqueous and membrane-mimicking environments. Simulation studies suggested that the amphipathic design of buCATHL4 was crucial for water permeation following membrane disruption. A great diversity, broad-spectrum antimicrobial action, and ability to induce an inflammatory response indicated the pleiotropic role of cathelicidins in innate immunity of buffalo. This study suggests short buffalo cathelicidin peptides with potent bactericidal properties and low cytotoxicity have potential translational applications for the development of novel antibiotics and antimicrobial peptidomimetics.

The NADPH-dependent HC-toxin reductases (HCTR1 and 2) encoded by enzymatic class of disease resistance homologous genes (Hm1 and Hm2) protect maize by detoxifying a cyclic tetrapeptide, HC-toxin, secreted by the fungus Cochliobolus carbonum race 1(CCR1). Unlike the other classes' resistance (R) genes, HCTR-mediated disease resistance is an inimitable mechanism where the avirulence (Avr) component from CCR1 is not involved in toxin degradation. In this study, we attempted to decipher cofactor (NADPH) recognition and mode of HC-toxin binding to HCTRs through molecular docking, molecular dynamics (MD) simulations and binding free energy calculation methods. The rationality and the stability of docked complexes were validated by 30-ns MD simulation. The binding free energy decomposition of enzyme-cofactor complex was calculated to find the driving force behind cofactor recognition. The overall binding free energies of HCTR1-NADPH and HCTR2-NADPH were found to be -616.989 and -16.9749 kJ mol-1 respectively. The binding free energy decomposition revealed that the binding of NADPH to the HCTR1 is mainly governed by van der Waals and nonpolar interactions, whereas electrostatic terms play dominant role in stabilizing the binding mode between HCTR2 and NADPH. Further, docking analysis of HC-toxin with HCTR-NADPH complexes showed a distinct mode of binding and the complexes were stabilized by a strong network of hydrogen bond and hydrophobic interactions. This study is the first in silico attempt to unravel the biophysical and biochemical basis of cofactor recognition in enzymatic class of R genes in cereal crop maize.

The cDNA sequence of leuteinizing hormone (LH) beta subunit was determined to understand the molecular basis for silent oestrus behavior and poor response to superovulation in buffalo. The LH-beta cDNA contains an open reading frame 426 bp long. The deduced sequence of the LH-beta is 141 amino acids in length. The amino acid sequences of the Indian river buffalo LH-beta subunit showed overall similarity to those of other mammals. The nucleotide sequence variability of LH-beta was studied in more than approximately 350 Indian buffaloes covering five different breeds. The results of the sequence analysis showed that the buffalo LH-beta gene is not highly conserved and non-synonymous mutations are not rare, at least in the samples collected randomly from five different breeds and buffalo populations. A total of seven different variants were obtained. In spite of its crucial role in reproduction, variation of the LH-beta gene was found present in this species. The polymorphisms found were unique in the Indian river buffalo population.

Although, it is known that spermatozoa harbor a variety of RNAs that may influence embryonic development, little is understood about sperm transcriptomic differences in relation to fertility, especially in buffaloes. In the present study, we compared the differences in sperm functional attributes and transcriptomic profile between high- and low-fertile buffalo bulls. Sperm membrane and acrosomal integrity were lower (P < 0.05), while protamine deficiency and lipid peroxidation were higher (P < 0.05) in low- compared to high-fertile bulls. Transcriptomic analysis using mRNA microarray technology detected a total of 51,282 transcripts in buffalo spermatozoa, of which 4,050 transcripts were differentially expressed, and 709 transcripts were found to be significantly dysregulated (P < 0.05 and fold change >1) between high- and low-fertile bulls. Majority of the dysregulated transcripts were related to binding activity, transcription, translation, and metabolic processes with primary localization in the cell nucleus, nucleoplasm, and in cytosol. Pathways related to MAPK signaling, ribosome pathway, and oxidative phosphorylation were dysregulated in low-fertile bull spermatozoa. Using bioinformatics analysis, we observed that several genes related to sperm functional attributes were significantly downregulated in low-fertile bull spermatozoa. Validation of the results of microarray analysis was carried out using real-time qPCR expression analysis of selected genes (YBX1, ORAI3, and TFAP2C). The relative expression of these genes followed the same trend in both the techniques. Collectively, this is the first study to report the transcriptomic profile of buffalo spermatozoa and to demonstrate the dysregulation of functionally relevant transcripts in low-fertile bull spermatozoa. The results of the present study open up new avenues for understanding the etiology for poor fertility in buffalo bulls and to identify fertility biomarkers.

Sperm, which are believed to be transcriptionally and translationally inactive, synthesize RNA and proteins before there is gradual disappearance of the ribosome during chromatin compaction. Sperm transfer several functionally relevant transcripts to the oocyte, controlling maternal-zygotic transition and embryonic development. The present study was undertaken to profile and analyze sperm transcripts comprehensively using Next Generation Ribonucleic acid sequencing technology in Holstein Friesian x Tharparkar crossbred bulls. The results from global transcriptomic profiling revealed transcripts for 13,814 genes; of which 431 transcripts were expressed with >1 FPKM and 13,383 transcripts were expressed with >0 or <1 FPKM. The abundant mRNA transcripts of crossbred bull sperm were PRM1 and HMGB4. Gene ontology of transcripts with>1 FPKM revealed there was a major involvement in the structural constituent of ribosomes and translation. Results from pathway enrichment indicated the connection between ribosome, oxidative phosphorylation and spliceosome pathways and the transcripts of crossbred bull spermatozoa. The transcriptional abundance of selected genes, validated using RT-qPCR, indicated significant variations between bulls. Collectively, it may be inferred that the transcripts in crossbred bull sperm were heavily implicated in functions such as the structural constituent of ribosomes and translation, and pathways such as ribosome, oxidative phosphorylation and spliceosome. Further studies using larger sample sizes are required to understand the possible implications of transcriptomic variations on semen quality and fertility.

RIG1 and MDA5 have emerged as important intracellular innate pattern recognition receptors that recognize viral RNA and mediate cellular signals controlling Type I interferon (IFN-I) response. Buffalo RIG1 and MDA5 genes were investigated to understand the mechanism of receptor induced antiviral response. Sequence analysis revealed that RIG1 and MDA5 maintain a domain arrangement that is common in mammals. Critical binding site residues of the receptors are evolutionary conserved among mammals. Molecular dynamics simulations suggested that RIG1 and MDA5 follow a similar, if not identical, dsRNA binding pattern that has been previously reported in human. Moreover, binding free energy calculation revealed that MDA5 had a greater affinity towards dsRNA compared to RIG1. Constitutive expressions of RLR genes were ubiquitous in different tissues without being specific to immune organs. Poly I:C stimulation induced elevated expressions of IFN-β and IFN-stimulated genes (ISGs) through interferon regulatory factors (IRFs) mediated pathway in buffalo foetal fibroblast cells. The present study provides crucial insights into the structure and function of RIG1 and MDA5 receptors in buffalo.

TLRs are some of the actively pursued drug-targets in immune disorders. Owing to a recent surge in the cognizance of TLR structural biology and signalling pathways, numerous therapeutic modulators, ranging from low-molecular-weight organic compounds to polypeptides and nucleic acid agents have been developed.