Characterization of rumen microbiome and immune genes expression of crossbred beef steers with divergent residual feed intake phenotypes.

We investigated whole blood and hepatic mRNA expressions of immune genes and rumen microbiome of crossbred beef steers with divergent residual feed intake phenotype to identify relevant biological processes underpinning feed efficiency in beef cattle. Low-RFI beef steers (n = 20; RFI = - 1.83 kg/d) and high-RFI beef steers (n = 20; RFI = + 2.12 kg/d) were identified from a group of 108 growing crossbred beef steers (average BW = 282 ± 30.4 kg) fed a high-forage total mixed ration after a 70-d performance testing period. At the end of the 70-d testing period, liver biopsies and blood samples were collected for total RNA extraction and cDNA synthesis. Rumen fluid samples were also collected for analysis of the rumen microbial community. The mRNA expression of 84 genes related to innate and adaptive immunity was analyzed using pathway-focused PCR-based arrays. Differentially expressed genes were determined using P-value ≤ 0.05 and fold change (FC) ≥ 1.5 (in whole blood) or ≥ 2.0 (in the liver). Gene ontology analysis of the differentially expressed genes revealed that pathways related to pattern recognition receptor activity, positive regulation of phagocytosis, positive regulation of vitamin metabolic process, vascular endothelial growth factor production, positive regulation of epithelial tube formation and T-helper cell differentiation were significantly enriched (FDR < 0.05) in low-RFI steers. In the rumen, the relative abundance of PeH15, Arthrobacter, Moryella, Weissella, and Muribaculaceae was enriched in low-RFI steers, while Methanobrevibacter, Bacteroidales_BS11_gut_group, Bacteroides and Clostridium_sensu_stricto_1 were reduced. In conclusion, our study found that low-RFI beef steers exhibit increased mRNA expression of genes related to immune cell functions in whole blood and liver tissues, specifically those involved in pathogen recognition and phagocytosis regulation. Additionally, these low-RFI steers showed differences in the relative abundance of some microbial taxa which may partially account for their improved feed efficiency compared to high-RFI steers.

The functions of long non-coding RNA (lncRNA)-MALAT-1 in the pathogenesis of renal cell carcinoma.

Renal cell carcinoma (RCC), a prevalent form of renal malignancy, is distinguished by its proclivity for robust tumor proliferation and metastatic dissemination. Long non-coding RNAs (lncRNAs) have emerged as pivotal modulators of gene expression, exerting substantial influence over diverse biological processes, encompassing the intricate landscape of cancer development. Metastasis-associated lung adenocarcinoma transcript 1 (MALAT-1), an exemplar among lncRNAs, has been discovered to assume functional responsibilities within the context of RCC. The conspicuous expression of MALAT-1 in RCC cells has been closely linked to the advancement of tumors and an unfavorable prognosis. Experimental evidence has demonstrated the pronounced ability of MALAT-1 to stimulate RCC cell proliferation, migration, and invasion, thereby underscoring its active participation in facilitating the metastatic cascade. Furthermore, MALAT-1 has been implicated in orchestrating angiogenesis, an indispensable process for tumor expansion and metastatic dissemination, through its regulatory influence on pro-angiogenic factor expression. MALAT-1 has also been linked to the evasion of immune surveillance in RCC, as it can regulate the expression of immune checkpoint molecules and modulate the tumor microenvironment. Hence, the potential utility of MALAT-1 as a diagnostic and prognostic biomarker in RCC emerges, warranting further investigation and validation of its clinical significance. This comprehensive review provides an overview of the diverse functional roles exhibited by MALAT-1 in RCC.

Population Genomics Reveals Incipient Speciation, Introgression, and Adaptation in the African Mona Monkey (Cercopithecus mona).

Guenons (tribe Cercopithecini) are the most widely distributed nonhuman primate in the tropical forest belt of Africa and show considerable phenotypic, taxonomic, and ecological diversity. However, genomic information for most species within this group is still lacking. Here, we present a high-quality de novo genome (total 2.90 Gb, contig N50 equal to 22.7 Mb) of the mona monkey (Cercopithecus mona), together with genome resequencing data of 13 individuals sampled across Nigeria. Our results showed differentiation between populations from East and West of the Niger River ∼84 ka and potential ancient introgression in the East population from other mona group species. The PTPRK, FRAS1, BNC2, and EDN3 genes related to pigmentation displayed signals of introgression in the East population. Genomic scans suggest that immunity genes such as AKT3 and IL13 (possibly involved in simian immunodeficiency virus defense), and G6PD, a gene involved in malaria resistance, are under positive natural selection. Our study gives insights into differentiation, natural selection, and introgression in guenons.

Transcriptome Profiling Reveals CD73 and Age-Driven Changes in Neutrophil Responses against Streptococcus pneumoniae.

Neutrophils are required for host resistance against Streptococcus pneumoniae, but their function declines with age. We previously found that CD73, an enzyme required for antimicrobial activity, is downregulated in neutrophils (also known as polymorphonuclear leukocytes [PMNs]) from aged mice. This study explored transcriptional changes in neutrophils induced by S. pneumoniae to identify pathways controlled by CD73 and dysregulated with age. Pure bone marrow-derived neutrophils isolated from wild-type (WT) young and old and CD73 knockout (CD73KO) young mice were mock challenged or infected with S. pneumoniae ex vivo. RNA sequencing (RNA-Seq) was performed to identify differentially expressed genes (DEGs). We found that infection triggered distinct global transcriptional changes across hosts that were strongest in CD73KO neutrophils. Surprisingly, there were more downregulated than upregulated genes in all groups upon infection. Downregulated DEGs indicated a dampening of immune responses in old and CD73KO hosts. Further analysis revealed that CD73KO neutrophils expressed higher numbers of long noncoding RNAs (lncRNAs) than those in WT controls. Predicted network analysis indicated that CD73KO-specific lncRNAs control several signaling pathways. We found that genes in the c-Jun N-terminal kinase (JNK)-mitogen-activated protein kinase (MAPK) pathway were upregulated upon infection in CD73KO mice and in WT old mice, but not in WT young mice. This corresponded to functional differences, as phosphorylation of the downstream AP-1 transcription factor component c-Jun was significantly higher in neutrophils from infected CD73KO mice and old mice. Importantly, inhibition of JNK/AP-1 rescued the ability of these neutrophils to kill S. pneumoniae. Together, our findings revealed that the ability of neutrophils to modify their gene expression to better adapt to bacterial infection is in part regulated by CD73 and declines with age.

The Protective Role of miR-125b in Hepatocellular Carcinoma: Unraveling Tumor-Suppressive Mechanisms.

MicroRNAs (miRNAs) have emerged as crucial regulators of gene expression, playing pivotal roles in various biological processes, including cancer development and progression. Among them, miR-125b has garnered significant attention due to its multifaceted functional roles in human hepatocellular carcinoma (HCC). Extensive research has revealed that miR-125b plays a dual role in HCC, acting as both a tumor suppressor and an oncogene depending on the context. As a tumor suppressor, miR-125b exerts its inhibitory effects on HCC by targeting key oncogenic pathways and genes involved in cell proliferation, migration, invasion, and angiogenesis. Its downregulation in HCC is frequently observed and correlates with aggressive tumor characteristics and poor prognosis. Conversely, miR-125b can also function as an oncogene in specific HCC subtypes or under certain conditions. It has been shown to promote HCC growth, metastasis, and therapeutic resistance by targeting tumor suppressor genes, modulating the epithelial-mesenchymal transition (EMT) process, and enhancing cancer stem cell-like properties. The upregulation of miR-125b in HCC has been associated with advanced disease stages and unfavorable clinical outcomes. Furthermore, the dysregulation of miR-125b expression in HCC is influenced by a complex network of regulatory mechanisms. Understanding these regulatory mechanisms is crucial for deciphering the precise functional roles of miR-125b in HCC and exploring its potential as a diagnostic biomarker or therapeutic target. In the current review study, we comprehensively elucidated the diverse functional roles of miR-125b in HCC, providing a comprehensive overview of its regulatory mechanisms and impact on key cellular processes involved in HCC progression.

Deciphering sperm functions using biological networks.

The global human population is exponentially increasing, which requires the production of quality food through efficient reproduction as well as sustainable production of livestock. Lack of knowledge and technology for assessing semen quality and predicting bull fertility is hindering advances in animal science and food animal production and causing millions of dollars of economic losses annually. The intent of this systemic review is to summarize methods from computational biology for analysis of gene, metabolite, and protein networks to identify potential markers that can be applied to improve livestock reproduction, with a focus on bull fertility. We provide examples of available gene, metabolic, and protein networks and computational biology methods to show how the interactions between genes, proteins, and metabolites together drive the complex process of spermatogenesis and regulate fertility in animals. We demonstrate the use of the National Center for Biotechnology Information (NCBI) and Ensembl for finding gene sequences, and then use them to create and understand gene, protein and metabolite networks for sperm associated factors to elucidate global cellular processes in sperm. This study highlights the value of mapping complex biological pathways among livestock and potential for conducting studies on promoting livestock improvement for global food security.

Molecular genotyping reveals mixed bovine and human trypanosomiasis in cattle from West Africa.

Background and Aim: Animal trypanosomiasis is a major contributor to agricultural and economic losses, especially in sub-Saharan Africa. We have shown that some animal species expressed genes that are significant players in immune response to bovine trypanosomosis, impeding signs and symptoms of the disease. We hypothesize that such animals are contributors to disease transmission dynamics and severe outcomes. Therefore, this study aims to ascertain trypanosome species diversity in cattle and their potential role as reservoirs for the transmission of human disease. Materials and Methods: We performed a molecular genotyping of trypanosome internal transcribed spacer 1 (ITS-1) and 18S ribosomal RNA genes on genomic DNA extracts from randomly sampled N'Dama cattle from slaughterhouses in Nigeria. We identified trypanosome species circulating among the animals through polymerase chain reaction and genomic sequencing. We performed multiple sequence alignments as well as conducted a phylogenetic relationship between identified species. Results: In all, 9 of 127 (7.1%) samples were positively amplified (band sizes ranging from 250 bp to 710 bp), including an isolate with two distinct bands (700 and 710 bp), indicating two trypanosome types. Sequence similarity and homology analysis identified four species, namely: Trypanosoma vivax, Trypanosoma congolense forest type, T. congolense savannah type, and Trypanosoma brucei. Interestingly, one of the bands, additionally verified by nucleotide sequencing, was identified as a human trypanosome (Trypanosoma brucei gambiense), confirming our hypothesis that cattle are potential reservoir hosts for human trypanosomes. Conclusion: Overall, we observed different trypanosome species in our study area, with animals on the same farm infected with multiple species, which could complicate treatment and disease control strategies. Finding human trypanosome species strengthens the argument that disease transmission dynamics are modulated by other vertebrates, further complicating control programs.

The differential plasma and ruminal metabolic pathways and ruminal bacterial taxa associated with divergent residual body weight gain phenotype in crossbred beef steers.

We applied ruminal and plasma metabolomics and ruminal 16S rRNA gene sequencing to determine the metabolic pathways and ruminal bacterial taxa associated with divergent residual body weight gain phenotype in crossbred beef steers. A group of 108 crossbred growing beef steers (average BW = 282.87 ± 30 kg) were fed a forage-based diet for a period of 56 d in a confinement dry lot equipped with GrowSafe intake nodes to determine their residual body weight gain (RADG) phenotype. After RADG identification, blood and rumen fluid samples were collected from beef steers with the highest RADG (most efficient; n = 16; 0.76 kg/d) and lowest RADG (least efficient; n = 16; -0.65 kg/d). Quantitative untargeted metabolome analysis of the plasma and rumen fluid samples were conducted using chemical isotope labelling/liquid chromatography-mass spectrometry. Differentially abundant metabolites in each of the plasma and rumen fluid samples between the two groups of beef steers were determined using a false discovery rate (FDR)-adjusted P-values ≤ 0.05 and area under the curve (AUC) > 0.80. Rumen and plasma metabolic pathways that were differentially enriched or depleted (P ≤ 0.05) in beef steers with positive RADG compared to those with negative RADG were determined by the quantitative pathway enrichment analysis. A total of 1,629 metabolites were detected and identified in the plasma of the beef steers; eight metabolites including alanyl-phenylalanine, 8-hydroxyguanosine, and slaframine were differentially abundant (FDR ≤ 0.05; AUC > 0.80) in beef steers with divergent RADG; five metabolic pathways including steroid hormone biosynthesis, thiamine metabolism, propanoate metabolism, pentose phosphate pathway, and butanoate metabolism were enriched (P ≤ 0.05) in beef steers with positive RADG, relative to negative RADG steers. A total of 1,908 metabolites were detected and identified in the rumen of the beef steers; results of the pathway enrichment analysis of all the metabolites revealed no metabolic pathways in the rumen were altered (P > 0.05). The rumen fluid samples were also analyzed using 16S rRNA gene sequencing to assess the bacterial community composition. We compared the rumen bacterial community composition at the genus level using a linear discriminant analysis effect size (LEfSe) to identify the differentially abundant taxa between the two groups of beef steers. The LEfSe results showed greater relative abundance of Bacteroidetes_vadinHA17 and Anaerovibrio in steers with positive RADG compared to the negative RADG group, while steers in the negative RADG group had greater relative abundance of Candidatus_Amoebophilus, Clostridium_sensu_stricto_1, Pseudomonas, Empedobacter, Enterobacter, and Klebsiella compared to the positive RADG group. Our results demonstrate that beef steers with positive or negative RADG exhibit differences in plasma metabolic profiles and some ruminal bacterial taxa which probably explain their divergent feed efficiency phenotypes.

Evolutionary conservation of motifs within vanA and vanB of vancomycin-resistant enterococci.

Background and Aim: Global Health is threatened by the rapid emergence of multidrug-resistant bacteria. Antibiotic resistomes rapidly evolve, yet conserved motifs elucidated in our study have the potential for future drug targets for precision medicine. This study aimed to identify conserved genetic sequences and their evolutionary pathways among vancomycin-resistant Enterococcus species such as Enterococcus faecium and Enterococcus faecalis. Materials and Methods: We retrieved a total of 26 complete amino acid and nucleotide sequences of resistance determinant genes against vancomycin (vanA and vanB), streptomycin (aac-aah), and penicillin (pbp5) from the publicly available genetic sequence database, GenBank. The sequences were comprised of bacteria classified under the genera of Enterococcus, Staphylococcus, Amycolatopsis, Ruminococcus, and Clostridium. Sequences were aligned with Clustal Omega Multiple Sequence Alignment program and Percent Identity Matrices were derived. Phylogenetic analyses to elucidate evolutionary relationships between sequences were conducted with the neighbor-end joining method through the Molecular Evolutionary Genetics Analysis (MEGAX) software, developed by the Institute of Molecular Evolutionary Genetics at Pennsylvania State University. Subsequent network analyses of the resistance gene, vanB, within E. faecium were derived from ScanProsite and InterPro. Results: We observed the highest nucleotide sequence similarity of vanA regions within strains of E. faecium (100%) and E. faecalis (100%). Between Enterococcus genera, we continued to observe high sequence conservation for vanA and vanB, up to 99.9% similarity. Phylogenetic tree analyses suggest rapid acquisition of these determinants between strains within vanA and vanB, particularly between strains of Enterococcus genera, which may be indicative of horizontal gene transfer. Within E. faecium, Adenosine 5'-Triphosphate (ATP)-Grasp and D-ala-D-ala ligase (Ddl) were found as conserved domains of vanA and vanB. We additionally found that there is notable sequence conservation, up to 66.67%, between resistomes against vancomycin and streptomycin among E. faecium. Conclusion: Resistance genes against vancomycin have highly conserved sequences between strains of Enterococcus bacteria. These conserved sequences within vanA and vanB encode for ATP-Grasp and Ddl motifs, which have functional properties for maintaining cell wall integrity. High sequence conservation is also observed among resistance genes against penicillin and streptomycin, which can inform future drug targets for broader spectrum therapies.

Population genetic analysis of Plasmodium falciparum cell-traversal protein for ookinetes and sporozoite among malaria patients from southern Nigeria.

Plasmodium falciparum immune escape mechanisms affect antigens being prioritized for vaccine design. As a result of the multiple surface antigens the parasite exhibits at different life cycle stages, designing a vaccine that would efficiently boost the immune system in clearing infections has been challenging. The P. falciparum cell-traversal protein for ookinetes and sporozoite (Pfceltos) is instrumental for ookinete traversal of the mosquito midgut and sporozoites invasion of the human liver cells. Pfceltos elicits both humoral and cellular immune response but has been reported with multiple single nucleotide polymorphisms in global isolates. A cross-sectional survey, conducted in southern Nigeria, between January-March 2021 recruited 283 individuals. Of this, 166 demonstrated P. falciparum infections (86 from Cross River and 80 from Edo), 48 (55.8%) while only 36 (45%) were amplified for Pfceltos gene from both sites respectively. Fifty amplified samples were sequenced and analysed for their diversity, polymorphisms and population structure of the gene. The number of segregating sites in Edo State was higher (34) than that of Cross River State. Though nucleotide diversity was higher for Edo compared to Cross River State (θw = 0.02505; π = 0.03993 versus θw = 0.00930; π = 0.01033 respectively), the reverse was the case for haplotype diversity (0.757 versus 0.890 for Edo and Cross River respectively). Of the twelve haplotypes observed from both states, only two (KASLPVEK and NAFLSFEK) were shared, with haplotype prevalence higher in Edo (16% and 36%) than Cross River (8% and 4%). The Tajima's D test was positive for both states, with Fst value showing a strong genetic differentiation (Fst = 0.25599), indicating the occurrence of balancing selection favoring haplotype circulation at a low frequency. The shared haplotypes, low Hst and Fst values presents a challenge to predict the extent of gene flow. High LD values present a grim public health consequence should a Pfceltos-conjugated vaccine be considered for prophylaxis in Nigeria.

Deciphering inhibitory mechanism of coronavirus replication through host miRNAs-RNA-dependent RNA polymerase interactome.

Despite what we know so far, Covid-19, caused by SARS-CoV-2 virus, remains a pandemic that still require urgent healthcare intervention. The frequent mutations of the SARS-CoV-2 virus has rendered disease control with vaccines and antiviral drugs quite challenging, with newer variants surfacing constantly. There is therefore the need for newer, effective and efficacious drugs against coronaviruses. Considering the central role of RNA dependent, RNA polymerase (RdRp) as an enzyme necessary for the virus life cycle and its conservation among coronaviruses, we investigated potential host miRNAs that can be employed as broad-range antiviral drugs averse to coronaviruses, with particular emphasis on BCoV, MERS-CoV, SARS-CoV and SARS-CoV-2. miRNAs are small molecules capable of binding mRNA and regulate expression at transcriptional or translational levels. Our hypothesis is that host miRNAs have the potential of blocking coronavirus replication through miRNA-RdRp mRNA interaction. To investigate this, we retrieved the open reading frame (ORF1ab) nucleotide sequences and used them to interrogate miRNA databases for miRNAs that can bind them. We employed various bioinformatics tools to predict and identify the most effective host miRNAs. In all, we found 27 miRNAs that target RdRp mRNA sequence of multiple coronaviruses, of which three - hsa-miR-1283, hsa-miR-579-3p, and hsa-miR-664b-3p target BCoV, SARS-CoV and SARS-CoV-2. Additionally, hsa-miR-374a-5p has three bovine miRNA homologs viz bta-miR-374a, bta-miR-374b, and bta-miR-374c. Inhibiting the expression of RdRp enzyme via non-coding RNA is novel and of great therapeutic importance in the control of coronavirus replication, and could serve as a broad-spectrum antiviral, with hsa-miR-1283, hsa-miR-579-3p, and hsa-miR-664b-3p as highly promising.

Molecular Evolution and Inheritance Pattern of Sox Gene Family among Bovidae.

Sox genes are an evolutionarily conserved family of transcription factors that play important roles in cellular differentiation and numerous complex developmental processes. In vertebrates, Sox proteins are required for cell fate decisions, morphogenesis, and the control of self-renewal in embryonic and adult stem cells. The Sox gene family has been well-studied in multiple species including humans but there has been scanty or no research into Bovidae. In this study, we conducted a detailed evolutionary analysis of this gene family in Bovidae, including their physicochemical properties, biological functions, and patterns of inheritance. We performed a genome-wide cataloguing procedure to explore the Sox gene family using multiple bioinformatics tools. Our analysis revealed a significant inheritance pattern including conserved motifs that are critical to the ability of Sox proteins to interact with the regulatory regions of target genes and orchestrate multiple developmental and physiological processes. Importantly, we report an important conserved motif, EFDQYL/ELDQYL, found in the SoxE and SoxF groups but not in other Sox groups. Further analysis revealed that this motif sequence accounts for the binding and transactivation potential of Sox proteins. The degree of protein-protein interaction showed significant interactions among Sox genes and related genes implicated in embryonic development and the regulation of cell differentiation. We conclude that the Sox gene family uniquely evolved in Bovidae, with a few exhibiting important motifs that drive several developmental and physiological processes.

Regulatory network of miRNA, lncRNA, transcription factor and target immune response genes in bovine mastitis.

Pre- and post-transcriptional modifications of gene expression are emerging as foci of disease studies, with some studies revealing the importance of non-coding transcripts, like long non-coding RNAs (lncRNAs) and microRNAs (miRNAs). We hypothesize that transcription factors (TFs), lncRNAs and miRNAs modulate immune response in bovine mastitis and could potentially serve as disease biomarkers and/or drug targets. With computational analyses, we identified candidate genes potentially regulated by miRNAs and lncRNAs base pair complementation and thermodynamic stability of binding regions. Remarkably, we found six miRNAs, two being bta-miR-223 and bta-miR-24-3p, to bind to several targets. LncRNAs NONBTAT027932.1 and XR_003029725.1, were identified to target several genes. Functional and pathway analyses revealed lipopolysaccharide-mediated signaling pathway, regulation of chemokine (C-X-C motif) ligand 2 production and regulation of IL-23 production among others. The overarching interactome deserves further in vitro/in vivo explication for specific molecular regulatory mechanisms during bovine mastitis immune response and could lay the foundation for development of disease markers and therapeutic intervention.

Dissecting Transcription Factor-Target Interaction in Bovine Coronavirus Infection.

Coronaviruses are RNA viruses that cause significant disease within many species, including cattle. Bovine coronavirus (BCoV) infects cattle and wild ruminants, both as a respiratory and enteric pathogen, and possesses a significant economic threat to the cattle industry. Transcription factors are proteins that activate or inhibit transcription through DNA binding and have become new targets for disease therapies. This study utilized in silico tools to identify potential transcription factors that can serve as biomarkers for regulation of BCoV pathogenesis in cattle, both for testing and treatment. A total of 11 genes were identified as significantly expressed during BCoV infection through literature searches and functional analyses. Eleven transcription factors were predicted to target those genes (AREB6, YY1, LMO2, C-Rel, NKX2-5, E47, RORAlpha1, HLF, E4BP4, ARNT, CREB). Function, network, and phylogenetic analyses established the significance of many transcription factors within the immune response. This study establishes new information on the transcription factors and genes related to host-pathogen interactome in BCoV infection, particularly transcription factors YY1, AREB6, LMO2, and NKX2, which appear to have strong potential as diagnostic markers, and YY1 as a potential target for drug therapies.

Evolutionarily Conserved Long Non-coding RNA Regulates Gene Expression in Cytokine Storm During COVID-19.

Coronavirus is a family of viruses including alpha-, beta-, gamma-, delta-coronaviruses. Only alpha- and betacoronaviruses have been observed to infect humans. Past outbreaks of SARS-CoV and MERS-CoV, both betacoronavirus, are the result of a spillover from animals. Recently, a new strain termed SARS-CoV-2 emerged in December 2019 in Wuhan, China. Severe cases of COVID-19, the disease caused by SARS-CoV-2, lead to acute respiratory distress syndrome (ARDS). One contributor to the development of ARDS is cytokine storm, an overwhelming inflammatory immune response. Long non-coding RNAs (lncRNAs) are genetic regulatory elements that, among many functions, alter gene expression and cellular processes. lncRNAs identified to be pertinent in COVID-19 cytokine storm have the potential to serve as disease markers or drug targets. This project aims to computationally identify conserved lncRNAs potentially regulating gene expression in cytokine storm during COVID-19. We found 22 lncRNAs that can target 10 cytokines overexpressed in COVID-19 cytokine storm, 8 of which targeted two or more cytokine storm cytokines. In particular, the lncRNA non-coding RNA activated by DNA damage (NORAD), targeted five out of the ten identified cytokine storm cytokines, and is evolutionarily conserved across multiple species. These lncRNAs are ideal candidates for further in vitro and in vivo analysis.

Integrative Network Analysis of Predicted miRNA-Targets Regulating Expression of Immune Response Genes in Bovine Coronavirus Infection.

Bovine coronavirus (BCoV) infection that causes disease outbreaks among farm animals, resulting in significant economic losses particularly in the cattle industry, has the potential to become zoonotic. miRNAs, which are short non-coding segments of RNA that inhibits the expression of their target genes, have been identified as potential biomarkers and drug targets, though this potential in BCoV remains largely unknown. We hypothesize that certain miRNAs could simultaneously target multiple genes, are significantly conserved across many species, thereby demonstrating the potential to serve as diagnostic or therapeutic tools for bovine coronavirus infection. To this end, we utilized different existing and publicly available computational tools to conduct system analysis predicting important miRNAs that could affect BCoV pathogenesis. Eleven genes including CEBPD, IRF1, TLR9, SRC, and RHOA, significantly indicated in immune-related pathways, were identified to be associated with BCoV, and implicated in other coronaviruses. Of the 70 miRNAs predicted to target the identified genes, four concomitant miRNAs (bta-miR-11975, bta-miR-11976, bta-miR-22-3p, and bta-miR-2325c) were found. Examining the gene interaction network suggests IL-6, IRF1, and TP53 as key drivers. Phylogenetic analysis revealed that miR-22 was completely conserved across all 14 species it was searched against, suggesting a shared and important functional role. Functional annotation and associated pathways of target genes, such as positive regulation of cytokine production, IL-6 signaling pathway, and regulation of leukocyte differentiation, indicate the miRNAs are major participants in multiple aspects of both innate and adaptive immune response. Examination of variants evinced a potentially deleterious SNP in bta-miR-22-3p and an advantageous SNP in bta-miR-2325c. Conclusively, this study provides new insight into miRNAs regulating genes responding to BCoV infection, with bta-miR-22-3p particularly indicated as a potential drug target or diagnostic marker for bovine coronavirus.

SNP Diversity in CD14 Gene Promoter Suggests Adaptation Footprints in Trypanosome Tolerant N'Dama (Bos taurus) but not in Susceptible White Fulani (Bos indicus) Cattle.

Immune response to infections has been shown to be mediated by genetic diversity in pattern recognition receptors, leading to disease tolerance or susceptibility. We elucidated naturally occurring variations within the bovine CD14 gene promoter in trypanosome-tolerant (N'Dama) and susceptible (White Fulani) cattle, with genomic and computational approaches. Blood samples were collected from White Fulani and N'Dama cattle, genomic DNA extracted and the entire promoter region of the CD14 gene amplified by PCR. We sequenced this region and performed in silico computation to identify SNP variants, transcription factor binding sites, as well as micro RNAs in the region. CD14 promoter sequences were compared with the reference bovine genome from the Ensembl database to identify various SNPs. Furthermore, we validated three selected N'Dama specific SNPs using custom Taqman SNP genotyping assay for genetic diversity. In all, we identified a total of 54 and 41 SNPs at the CD14 promoter for N'Dama and White Fulani respectively, including 13 unique SNPs present in N'Dama only. The significantly higher SNP density at the CD14 gene promoter region in N'Dama may be responsible for disease tolerance, possibly an evolutionary adaptation. Our genotype analysis of the three loci selected for validation show that mutant alleles (A/A, C/C, and A/A) were adaptation profiles within disease tolerant N'Dama. A similar observation was made for our haplotype analysis revealing that haplotypes H1 (ACA) and H2 (ACG) were significant combinations within the population. The SNP effect prediction revealed 101 and 89 new transcription factor binding sites in N'Dama and White Fulani, respectively. We conclude that disease tolerant N'Dama possessing higher SNP density at the CD14 gene promoter and the preponderance of mutant alleles potentially confirms the significance of this promoter in immune response, which is lacking in susceptible White Fulani. We, therefore, recommend further in vitro and in vivo study of this observation in infected animals, as the next step for understanding genetic diversity relating to varying disease phenotypes in both breeds.

Differential Evolution of α-Glucan Water Dikinase (GWD) in Plants.

The alpha-glucan water dikinase (GWD) enzyme catalyzes starch phosphorylation, an integral step in transitory starch degradation. The high phosphate content in stored starch has great industrial value, due to its physio-chemical properties making it more versatile, although the phosphate content of stored starch varies depending on the botanical source. In this study, we used various computational approaches to gain insights into the evolution of the GWD protein in 48 plant species with possible roles in enzyme function and alteration of phosphate content in their stored starch. Our analyses identified deleterious mutations, particularly in the highly conserved 5 aromatic amino acid residues in the dual tandem carbohydrate binding modules (CBM-45) of GWD protein in C. zofingiensis, G. hirsutum, A. protothecoides, P. miliaceum, and C. reinhardtii. These findings will inform experimental designs for simultaneous repression of genes coding for GWD and the predicted interacting proteins to elucidate the role this enzyme plays in starch degradation. Our results reveal significant diversity in the evolution of GWD enzyme across plant species, which may be evolutionarily advantageous according to the varying needs for phosphorylated stored starch between plants and environments.

RNA-seq profiling of skin in temperate and tropical cattle.

Skin is a major thermoregulatory organ in the body controlling homeothermy, a critical function for climate adaptation. We compared genes expressed between tropical- and temperate-adapted cattle to better understand genes involved in climate adaptation and hence thermoregulation. We profiled the skin of representative tropical and temperate cattle using RNA-seq. A total of 214,754,759 reads were generated and assembled into 72,993,478 reads and were mapped to unique regions in the bovine genome. Gene coverage of unique regions of the reference genome showed that of 24,616 genes, only 13,130 genes (53.34%) displayed more than one count per million reads for at least two libraries and were considered suitable for downstream analyses. Our results revealed that of 255 genes expressed differentially, 98 genes were upregulated in tropically-adapted White Fulani (WF; Bos indicus) and 157 genes were down regulated in WF compared to Angus, AG (Bos taurus). Fifteen pathways were identified from the differential gene sets through gene ontology and pathway analyses. These include the significantly enriched melanin metabolic process, proteinaceous extracellular matrix, inflammatory response, defense response, calcium ion binding and response to wounding. Quantitative PCR was used to validate six representative genes which are associated with skin thermoregulation and epithelia dysfunction (mean correlation 0.92; p < 0.001). Our results contribute to identifying genes and understanding molecular mechanisms of skin thermoregulation that may influence strategic genomic selection in cattle to withstand climate adaptation, microbial invasion and mechanical damage.

Evolutionary Pattern of Interferon Alpha Genes in Bovidae and Genetic Diversity of IFNAA in the Bovine Genome.

Interferons are secretory proteins induced in response to specific extracellular stimuli which stimulate intra- and intercellular networks for regulating innate and acquired immunity, resistance to viral infections, and normal and tumor cell survival and death. Type 1 interferons plays a major role in the CD8 T-cell response to viral infection. The genomic analysis carried out here for type I interferons within Bovidae family shows that cattle, bison, water buffalo, goat, and sheep (all Bovidae), have different number of genes of the different subtypes, with a large increase in the numbers, compared to human and mouse genomes. A phylogenetic analysis of the interferon alpha (IFNA) proteins in this group shows that the genes do not follow the evolutionary pattern of the species, but rather a cycle of duplications and deletions in the different species. In this study we also studied the genetic diversity of the bovine interferon alpha A (IFNAA), as an example of the IFNA genes in cattle, sequencing a fragment of the coding sequence in 18 breeds of cattle from Pakistan, Nigeria and USA. Similarity analysis allowed the allocation of sequences into 22 haplotypes. Bhagnari, Brangus, Sokoto Gudali, and White Fulani, had the highest number of haplotypes, while Angus, Hereford and Nari Master had the least. However, when analyzed by the average haplotype count, Angus, Bhagnari, Hereford, Holstein, Muturu showed the highest values, while Cholistani, Lohani, and Nari Master showed the lowest values. Haplotype 4 was found in the highest number of individuals (74), and in 15 breeds. Sequences for yak, bison, and water buffalo, were included within the bovine haplotypes. Medium Joining network showed that the sequences could be divided into 4 groups: one with highly similar haplotypes containing mostly Asian and African breeds, one with almost all of the Bos taurus American breeds, one mid-diverse group with mostly Asian and African sequences, and one group with highly divergent haplotypes with five N'Dama sequences and one from each of White Fulani, Dhanni, Tharparkar, and Bhagnari. The large genetic diversity found in IFNAA could be a very good indication of the genetic variation among the different genes of IFNA and could be an adaptation for these species in response to viral challenges they face.