Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 30 September 2020
Sec. Computational Genomics
This article is part of the Research Topic Bioinformatics Analysis of Omics Data for Biomarker Identification in Clinical Research View all 45 articles

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

  • 1Department of Biology, Hamilton College, Clinton, NY, United States
  • 2Department of Biomedical Sciences, College of Health Sciences and Technology, Rochester Institute of Technology, Rochester, NY, United States

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.

Introduction

Coronaviruses are pathogenic RNA viruses of the subfamily Coronovirinae, engendering severe disease in many bird and mammalian species, such diseases becoming fatal in chickens, pigs, cattle, dogs, and mice (Fehr and Perlman, 2015; Amer, 2018; Ellis, 2019). They are responsible for the ongoing COVID-19 pandemic and the attendant significant morbidity and mortality (Gussow et al., 2020; Lucas et al., 2020). Bovine coronavirus (BCoV) is endemic worldwide and affects both wild and domesticated cattle (Saif, 2010), causing disease affecting the respiratory system and gastrointestinal tract, with pathogenesis ranging from mild to fatal (Underwood et al., 2015; Oma et al., 2016). BCoV has high mortality rates in calves, inducing wasting and flattening of the villi in the GI and a substantial cause of calf diarrhea. It is indicated in winter dysentery, and drives the consequent decrease in milk production in adult cattle. It results in significant economic losses in both the dairy and beef industries, and an outbreak in meatpacking industries can alter the supply chain and meat distribution system of any country (Amer, 2018). To date, there is no effective vaccine for BCoV, which is of significant concern, since the virus has been demonstrated to show low fidelity to its host, infecting other animal in the process (Fehr and Perlman, 2018; Saif and Jung, 2020). Due to the paucity of information on bovine coronaviruses, the possible consequences arising from the current pandemic and its potential to become zoonotic, more research is needed to understand and disentangle host-virus interaction in the context of immune response and disease pathogenesis.

Bovine coronavirus is transmitted via the respiratory or oral-fecal route, and it has been documented that recurrent nasal shedding occurs in infected cattle. Infection is detected using fecal or nasal swabs analyzed via real-time reverse transcriptase-polymerase chain reactions (RT-PCR) that detect nucleocapsid (N) or transmembrane (M) protein-coding genes of the virus (Takiuchi et al., 2006; Decaro et al., 2008; Oma et al., 2016). BCoV has a highly glycosylated spike (S) protein as well as a hemagglutinin esterase (HE), both major participants in the attachment of the virus to the host cell (Saif and Jung, 2020). Viral entry induces activation of both signaling and inflammatory receptors that nurture differential expression of genes, mediating host immune response in the process. Due to the close evolutionary relationship between BCoV and severe acute respiratory syndrome (SARS-CoV), there is a potential for overlap in immune response pathway that benefit our current efforts. It is expected that bovine coronavirus will instigate the response of particular immune-related genes in cattle, with toll-like receptor 9 (TLR9), interleukin-6 (IL-6), and IL-8 implicated in this process (Breuer et al., 2012). Significantly, these genes in the presence of exogenous stimuli, have been demonstrated to mediate the inflammatory response and initiation of innate immunity, regulation of the MAPK cascade, and regulation of T cell proliferation, a major requirement for adaptive and long-lived immune response (Breuer et al., 2012). MicroRNAs have been proposed as players in the regulation process and now occupy a central role in the host immune response to BCoV.

MicroRNAs (miRNAs) are small non-coding regions of RNA which act as post-transcriptional regulatory factors and lessen the expression of their target genes, by either degrading the target mRNAs or inhibit their translation (Lu and Rothenberg, 2018). Regions of their sequences are complementary to that of their target mRNAs, and only 7−8 base-pairs are adequate for target-binding (Mohr and Mott, 2015). A single miRNA can have many target genes, and miRNAs often influence the expression of numerous genes within a functioning pathway (Simonson and Das, 2015; Lu and Rothenberg, 2018). It is thought that more than half of protein-coding mammalian genes are regulated by miRNAs. Bovine miRNAs such as miR-125b and miR-141 have been found to have effects on lactation, miR-122 regulates cholesterol in the liver, and miR-133 controls proliferation of skeletal muscle (Krützfeldt et al., 2005; Chen et al., 2006; Li et al., 2012). In a disease system, the expression of miRNAs can be reduced or amplified (Simonson and Das, 2015), with bovine miRNAs especially found to be differentially expressed upon viral infection (Glazov et al., 2009). miRNAs have shown potential as diagnostic markers, and have been used to distinguish between types of cancer (Mishra et al., 2016), and have been identified as potential biomarkers and promising therapeutic agents (Lu and Rothenberg, 2018). Such therapeutic potential is due to their dynamic ability to target multiple proteins and pathways (Simonson and Das, 2015). They have been studied as a treatment for atherosclerosis, hepatitis C virus, and cancers (Scheel et al., 2017). Additionally, miRNA have been demonstrated to have significant implications in bovine disease elucidation, including immune response to viral infection (Vidigal and Ventura, 2015; Scheel et al., 2017; Morenikeji et al., 2019).

Numerous existing algorithms have proven useful as tools for predicting miRNA targets using seed matches or other compensating factors (Griffiths-Jones et al., 2008; Agarwal et al., 2015; Fan et al., 2016; Hanif et al., 2018). Databases are available to compile associations and interactions between proteins for purposes such as identifying disease genes or biomarkers (Hayashida and Akutsu, 2016; Khurana et al., 2017). Software to assemble functional annotation and significantly associated pathways for genes enhances our understanding of their biological processes (Gene Ontology Consortium, 2004; Mi et al., 2019). Computational tools for predicting the effects of variants have been used to recognize the need for individual preventative treatments and to determine which variants play a role in disease susceptibility (Weil and Chen, 2011; Saint Pierre and Génin, 2014; McLaren et al., 2016).

In this study, miRNAs are predicted to play a significant regulatory role relating to immune response to BCoV infection. Notably, we found no record concerning miRNAs as regulatory factors in bovine coronavirus infection. We utilized a combination of different bioinformatic software to identify miRNAs that target BCoV-associated genes, their functional roles, and existing variants. Increased understanding of the regulatory function of miRNA in relation to BCoV could be valuable for gene therapy and drug development. This study aims to elucidate the role of miRNAs in the host immune response to BCoV, demonstrate miRNA-target interactions, identify potential markers for disease diagnosis, and possible therapeutic drug targets.

Materials and Methods

Determination of Genes Associated With Bovine Coronavirus Infection

Figure 1A shows the pipeline we utilized for this study, including a step by step flowchart starting from literature search to miRNA variant identification, as previously described (Morenikeji et al., 2019). To begin, we searched for recently published, publicly available literature, pertaining to gene expression, and particularly bovine coronavirus from biomedical science repository PubMed and Google Scholar. Using “BCoV” as the keyword returns 166 results, out of which we identified 77 genes as significant and were included for further analysis (Aich et al., 2007). These 77 genes were used as seed to query Gene Ontology database1 with the “Bos taurus” and “biological process” as reference. The biological pathways associated with inputted genes and published papers were retrieved. The results were sorted by fold enrichment above 70 and p-Values less than 0.005 which eventually catalog a total of 31 significant pathways for further analysis. In all, seventeen genes were reported in these pathways and cataloged for literature search to determine which of these genes have been reported in relation to the greatest number of other coronavirus species (murine, avian, SARS-CoV, feline, porcine, and canine were searched for in particular). Seven genes returned zero supporting articles and were removed from further analysis. The gene NANS (N-acethylneuraminate synthase) was included in the final list of 11 genes due to its relation to the receptor for bovine coronavirus. The eleven genes were then used as seed for further computational analyses.

FIGURE 1
www.frontiersin.org

Figure 1. Flowchart showing steps in the project design and analyses process (A), while the pie chart indicates how the number of genes associated with bovine coronavirus were sorted according to significance. (B) The blue represents genes that were only from the article by Aich et al. (2007). The orange are genes from the article that also appeared in significant pathways from GO. The green are genes from the article, from significant pathways, and found in supporting literature. The green also includes the gene NANS that was added.

Search for miRNAs That Target Identified Significant Genes

Using NCBI2, the complete sequences were collected for the 11 genes determined to be significant in bovine coronavirus infection. In order to predict which miRNAs could target the 11 genes, multiple databases were searched for each gene. miRBase (Kozomara and Griffiths-Jones, 2014; Kozomara et al., 2019) was searched using the mRNA sequence of a gene with the “Bos taurus” as reference and an E-value cutoff of 12. Also, miRNet (Fan et al., 2016) was searched with Bos taurus as the reference species, using the “Official gene symbol,” and “targeted by: miRNA.” Additionally, TargetScanHuman (Agarwal et al., 2015) was searched using the gene symbol, and only the “highly conserved” resulting miRNAs were recorded. The returned miRNAs from the databases miRBase, miRNet, and TargetScanHuman were compiled for a total of 69 miRNAs. To determine which miRNAs could target multiple of the 11 genes, an online Venn diagram creator3, which determines overlapping elements in the entered lists and creates Venn diagrams to represent them was utilized.

Creation of Gene−Gene Interaction and miRNA Network

Using STRING4, a network for the final 11 genes was created. STRING is a database containing both experimentally determined and predicted interactions between proteins. It creates a visual representation of these interactions. In order to tie the identified miRNAs, target genes and pathways in to a critical network and infer significant biological process/mechanism for gene regulation, Cytoscape version 3.7.2 was utilized to create a network that incorporated the miRNAs, their target genes, and the pathways that the target genes are associated with as described (Oliveira et al., 2018; Lucas et al., 2020).

Evolutionary Analysis of Concomitant miRNAs Across Gene Set

Three miRNAs that were identified in the previous analysis to have overlapping and multiple targets among the 11 genes included miRNAs bta-22-3p, bta-miR-2325c, and bta-miR-11976 were regarded most significant. The sequences of these concomitant miRNAs were retrieved from multiple sources including miRBase and GenBank. To conduct evolutionary trace analysis, sequences for the concomitant miRNAs from other species were retrieved using Ensembl BLAST5. Sequences were retrieved from fifteen species, including gorilla, chimpanzee, human, pig, wild yak, horse, cow, goat, sheep, cat, dog, chicken, mouse, rat, and bat (Table 3). For bta-miR-2325c and bta-miR-11976, the longer NCBI sequence was necessary to return results. One sequence from each species returned was used. Only the region of the returned sequence that aligned with the input sequence, and brief flanking regions, were included for the creation of the phylogeny. MEGA version X was used to align the sequences using MUSCLE alignment for the individual miRNAs and all of them combined (Kumar et al., 2018; Stecher et al., 2020). For the individual miRNAs, a “neighbor-joining” tree was created using the bootstrap method in MEGA version X. A tree containing all three miRNAs was created using iTOL6. The display mode was set to “circular.” Phylogenetic analysis was performed in order to determine evolutionary conservation of miRNAs bta-22-3p, bta-miR-2325c, and bta-miR-11976 across species.

Functional Enrichment of Identified Gene Targets by Concomitant miRNAs

Enrichment analysis was performed using Gene Ontology’s PANTHER tool for the seven genes that were targeted by concomitant miRNAs. These genes were CEBPB, SRC, TLR9, CEBPD, IRF1, NANS, and RHOA. They were searched under “Bos taurus” and “biological process.” The results were sorted by the p-Values and FDR as described (Oliveira et al., 2018). The significant pathways that related to immune response, and the genes indicated in each, were recorded. For more specific pathways, GeneAnalytics7 was utilized. As the species cow is not supported, the seven genes were searched under “human symbols.” The immune-related pathways, their score, and the genes involved in the pathway were recorded. Only the pathways of “high” quality according to GeneAnalytics were considered. miRBase, GenBank UniProt, miRCarta, miRGator, Ensembl, and NON-CODE were searched for functional annotation for the concomitant miRNAs.

Ascertainment of Variants Within Concomitant miRNAs

The sequences of the four concomitant miRNAs were searched against the cow genome via Ensembl BLAST. The variant feature of Ensembl was used in order to examine the occurrence of variants as they could affect the functionality of the miRNAs. The flanking regions were not considered. For each variant returned, the SNP ID, class, most severe consequence, change tolerance, and location was recorded. The change tolerance score is the Genomic Evolutionary Rate Profiling (GERP) score. The definition of a GERP score is “the reduction in the number of substitutions in the multi-species sequence alignment compared to the neutral expectation” (Huber et al., 2020). GERP scores reflect “rejected substitutions” and quantify negative selection experienced at a particular location (Cooper et al., 2005; Pollard et al., 2010). Previously, and in this study, a cutoff of an absolute value of greater than or equal to 2 was used to determine important variants (Jha et al., 2015).

Results

Identification of Significant Genes in BCoV Infection

The article by Aich et al. (2007) was used as the initial source from which to catalog the differentially expressed genes in bovine coronavirus infection (Supplementary Table 1). Since literature pertaining to gene expression in BCoV is limited, further stipulations were placed on the 77 seed genes to ascertain if the gene was significant in the disease. These included association with significant immune-related pathways and supporting literature from other coronaviruses (Table 1). The genes TLR9, RELA and TLR7 appeared in 11, 9, and 8 out of the 31 returned pathways respectively. Interestingly, TLR7 featured in the greatest number of sources pertaining to other coronaviruses. It has been indicated in porcine, avian, human, murine, and feline coronavirus infections (Calzada-Nova et al., 2010; Kameka et al., 2014; Li et al., 2016; Athmer et al., 2018; Malbon et al., 2019). A total of 11 genes comprised the final list: TLR7, TLR9, IRF1, IL-6, CEBPD, TRAF2, RHOA, TP53, CEBPB, SRC, and NANS (Table 1). Figure 1B summarizes the diminution in the number of genes considered for further analysis.

TABLE 1
www.frontiersin.org

Table 1. Lists of the genes indicated to be significant in bovine coronavirus infection.

Discovery of miRNAs That Target Significant Genes

Many miRNAs were predicted for the list of 11 genes as depicted in Table 2. The gene predicted to be targeted by the greatest number of miRNAs was IRF1, which is targeted by 21 miRNAs such as bta-miR-240, bta-miR-22-3p, bta-miR-214, and others. It was followed by TP53 which had 14 predicted miRNAs. Of all the genes, NANS and CEBPD are targeted by the fewest miRNAs with only one predicted miRNA for each. In order to determine the significant miRNAs that could potentially serve as drug targets, concomitant miRNAs were evaluated, with the assumption that miRNAs targeting multiple genes would be choice therapeutic candidates. Both bta-miR-11976 and bta-miR-11975 were predicted to target the genes CEBPB, SRC, and TLR9 (Figure 2A). The genes CEBPD and IRF1 are both targeted by bta-miR-22-3p (Figure 2B) while bta-miR-2325c was predicted to target NANS and RHOA (Figure 2C). These concomitant miRNAs (bta-miR-11976, bta-miR-11975, bta-miR-22-3p, and bta-miR-2325c) were used for further analysis.

TABLE 2
www.frontiersin.org

Table 2. List of miRNAs predicted by miRBase, miRNet, or TargetScanHuman for the 11 significant genes.

FIGURE 2
www.frontiersin.org

Figure 2. (A) Venn diagram to show that CEBPB, SRC, and TLR9 can all be targeted by the miRNAs bta-miR-11976 and bta-miR-11975; (B) show that bta-miR-22-3p can target both CEBPD and IRF1; (C) showing that NANS and RHOA are both targeted by bta-miR-2325c.

Network Analysis of Genes Significant in BCoV

To understand the gene-gene interaction or miRNA-target association that could guide the selection of possible therapeutic agents or discovery of biomarkers for disease diagnosis, we examined the gene relationships in BCoV infection. The online application STRING version 11.0, was used to create a network for the 11 final genes as shown; two of the genes, CEBPB and CEBPD were not found giving a total of nine genes in the network. Of the nine genes that were found in STRING network, only TRAF2 and NANS were not connected with others in the network (Figure 3). IL-6, IRF1, and TP53 are depicted as key drivers of the network. IL-6 is connected to five other genes including TP53, SRC, IRF1, TLR9, and TLR7. IRF1 and TP53 are both connected to four other genes. Many of the interactions, represented in light green, were determined via text-mining. The black lines represent co-occurrence. Two of the interactions come from curated databases, and three of them have been experimentally determined. Two of the interactions have been identified from four different tools. These interactions are between RHOA and SRC and between TLR9 and TLR7.

FIGURE 3
www.frontiersin.org

Figure 3. Network for nine of the final genes created using STRING. The Ensembl ID represents cytokine IL-6.

Figure 4 depicts the interactions between all of the predicted miRNAs, the 11 final genes, the interactions between the genes and the pathways that the genes are associated with. The diamonds represent the miRNAs. miRNAs that were only predicted for one gene appear in green while those predicted to target two or three genes appear in red. Pathways are represented by blue rectangles and other genes are shown in pink. The genes CEBPB, TLR9, IRF1, and RHOA appear as major hubs in the network having several interactions with other elements in the network, and likely having more influence. The genes TLR7, NANS, and CEBPD are connected to the network by only a few interactions. The genes IL-6 and TRAF2 are isolated with connections only to predicted non-concomitant miRNAs.

FIGURE 4
www.frontiersin.org

Figure 4. A network incorporating the identified genes, the miRNAs that target them, the interactions between them, and the pathways that they are indicated in.

Elucidating Evolutionary Relationships of Concomitant miRNAs Across Species

Our hypothesis is that miRNAs with lower rates of evolution are conserved across many species and might be candidates to study immune gene regulation in BCoV. Therefore, we performed an evolutionary analysis of the selected candidates from our previous analysis. Information pertaining to the sequences that were used to create the phylogenies for bta-miR-22-3p, bta-miR-2325c, and bta-miR-11975 is shown in Table 3. For bta-miR-11976, no results were found for any of the species searched and was removed from further analysis while miRNA-22 was found in 14 species. The alignment analysis revealed that miR-22 was perfectly conserved across all species that it appeared in Figure 5. From the phylogenetic tree in Figure 6A, it is interesting to note that Capra hircus and Sus scrofa are closely related. The other ruminants, Ovis aries, Bos taurus, and Bos mutus, are clustered together with Ovis aries as the outgroup. While Capra hircus and Sus scrofa are not closely related species, their miR-22s are indicated to be. Surprisingly, we found Rattus norvegicus to be part of the same clade as Gallus gallus rather than Mus musculus. Lastly, Gorilla gorilla gorilla and Homo sapiens most recent common ancestor is at the node with a bootstrap value of 67 and are not sister taxa, as expected.

TABLE 3
www.frontiersin.org

Table 3. Characteristics of sequences used to create the phylogenies for bta-miR-22-3p, bta-miR-2325c, and bta-miR-11975.

FIGURE 5
www.frontiersin.org

Figure 5. The alignment of miR-22. The yellow box encompasses the sequence of miR-22 for Bos taurus.

FIGURE 6
www.frontiersin.org

Figure 6. The phylogeny of individual miRNA; (A) miRNA-22; (B) miR-2325c; (C) miR-11975 (Fewer species returned results for miR-11975 than compared to miR-22 and miR-2325c).

miR-11975 and miR-2325c did not display as much conservation as miR-22. For miR-2325c, 14 species such as sheep, mouse, rat, gorilla, and humans were returned from BLAST search results. It is interesting that Gorilla gorilla gorilla is the outgroup for Homo sapiens and Pteropus vampyrus (Figure 6B). Pan troglodytes is multiple nodes away from this grouping, and in the same clade as Equus caballus. All four ruminants cluster together, but Capra hircus and Ovis aries are sister taxa with Bos mutus and Bos taurus as the outgroups. Once again, Rattus norvegicus and Mus musculus are not indicated to be closely related.

miR-11975 was only found in 9 species including pig, horse, goat, cow, wild yak, chicken, megabat, cat, and dog. It is interesting to note that Bos taurus and Bos mutus are multiple nodes apart, whereas they were indicated to be closely related for miR-22 and miR-2325c (Figure 6C). Bos mutus is indicated as the outgroup to Gallus gallus and Pteropus vampyrus, which is surprising.

When plotting all of the miRNAs together, the three different miRNAs cluster together as seen in Figure 7. Sus scrofa miR-22 is the only sequence that is not clustered with other sequences of the same miRNA. It is interesting to note that Gallus gallus and Equus caballus are sister taxa for miR-22. Bos taurus is in the same clade as Bos mutus for miR-22 and Capra hircus for miR-11975. Homo sapiens and Gorilla gorilla gorilla are not sister groups in either miR-22 nor miR-2325c. It is especially interesting to note that Canis lupus familiaris miR-22 shares a node with branches that terminate in miR-2325c sequences.

FIGURE 7
www.frontiersin.org

Figure 7. The phylogenetic tree created using iTOL for the three miRNAs that returned homologous sequences from other species on Ensembl. miR-2325c is represented in blue, mir-22 in green, and miR-11975 in black.

Determining Function of Identified Genes Through Pathways

Specific pathways in which the genes targeted by concomitant miRNAs are indicated in Table 4. Gene Ontology returned the pathways numbered 1-5, while GeneAnalytics returned pathways 6-9. The greatest number of genes, five, was indicated in the “IL-2 pathway” (Figure 8). This pathway includes the genes CEBPB, CEBPD, TLR9, RHOA, and SRC (Table 4). The pathways with the fewest indicated genes, three, were “MAPK-Erk pathway,” “immune response IFN gamma signaling pathway,” “IL-6 signaling pathway,” and “positive regulation of leukocyte differentiation” (Figure 8). Particular combinations of genes appeared in multiple pathways. For example, the genes TLR9, IRF1, CEBPB, and RHOA appeared in four different pathways. bta-miR-22-3p can target CEBPD and IRF1. CEBPD was only present in the IL-6 and IL-2 pathways. IRF1 was indicated in all of the pathways except the “MAPK-Erk pathway” and the “IL-2 pathway.” bta-miR-22-3p will lessen all of the pathways of Table 4. CEBPB, SRC, and TLR9 were predicted to be targeted by bta-miR-11975 and bta-miR. CEBPB was present in all of the pathways including “IL-6 signaling pathway,” “regulation of lymphocyte activation,” and “immune response IFN gamma signaling pathway.” TLR9 was present in six of the pathways such as “regulation of lymphocyte activation” and “positive regulation of cytokine production.” SRC is present in “IL-2 pathway,” “immune response IFN gamma signaling pathway,” and “MAPK-Erk pathway.” As bta-miR-2325c targets RHOA, it will affect the regulation of lymphocyte and leukocyte activation, the regulation of leukocyte differentiation, and the IL-2 pathway. There is considerable overlap in the pathways that will be affected by the concomitant miRNAs. bta-miR-22-3p, bta-miR-11975, and bta-miR-11976 affect all of the pathways in Table 4.

TABLE 4
www.frontiersin.org

Table 4. Immune-related and specific pathways in which significant genes are indicated from Gene Ontology and GeneAnalytics.

FIGURE 8
www.frontiersin.org

Figure 8. Bar chart representing the number of genes indicated in each of the significant and specific immune-related pathways found using GO and GeneAnalytics.

Examination of Variants Within Concomitant miRNAs

Two of the miRNAs, miR-22, and miR-2325c, were found to have polymorphic variants (Table 5). miR-22-3p revealed only one variant, although it could result in a difference in the mature miRNA. The change tolerance for this variant was 3.61. bta-miR-2325c revealed eight variants, all of which could cause intron variance. The variants of bta-miR-2325c had change tolerances ranging from −2.4 to 1.15.

TABLE 5
www.frontiersin.org

Table 5. Information pertaining to the variants found for miR-22-3p and miR-2325c when compared to the cow genome.

Discussion

MicroRNAs are short RNA non-coding regions that repress the expression of their target genes post-transcriptionally, and have been implicated in multiple disease conditions, including as potential biomarkers for disease diagnosis and targets for development of therapeutics (Scheel et al., 2017; Lu and Rothenberg, 2018). They have been found to have numerous functions in Bovidae (Krützfeldt et al., 2005). Likewise, the suppression of a such as miRNA-122 has been shown to inhibit human hepatitis C propagation (van Rooij and Kauppinen, 2014) and miRNA-34 aids tumor suppressor in human liver cancer (van Rooij and Kauppinen, 2014; Vidigal and Ventura, 2015). Increasing work is being done to enhance our understanding of their role in bovine immune response to viruses (Lawless et al., 2014; Li et al., 2018; Morenikeji et al., 2019). To the best of our knowledge, the functionality of miRNAs in bovine disease, especially related to bovine coronavirus (BCoV) infection, remains poorly understood. BCoV is an economically costly RNA virus which affects the cattle industry through calf mortality and decreased milk production (Saif, 2010; Amer, 2018). No work has been conducted to determine the role of miRNAs in immune response to BCoV.

In this study, computational tools were employed to predict important miRNAs that target significant immune-related genes indicated in the host response to BCoV. Stringent criteria were placed on the genes whose expression reportedly differed upon infection with BCoV. Conditions included significant association with immune-related pathways and citation in host responses to other coronaviruses. Literature concerning other coronaviruses was considered due to the paucity of published reports regarding gene expression in BCoV. As coronaviruses are genetically similar to each other, it is likely that pathogenesis is similar as well (Siddell et al., 1983; Fehr and Perlman, 2015; McGruder and Leibowitz, 2015; Mousavizadeh and Ghasemi, 2020). The gene NANS, which codes for N-acetyl-9-0-acetylneuraminic acid (van Karnebeek et al., 2016), was included in our study due to its particular involvement in BCoV, despite not satisfying all the stipulations. This sialic acid has been determined to be the receptor onto which the S protein of BCoV binds (Schultze and Herrler, 1993; Belouzard et al., 2012; Park, 2019). Interestingly, NANS is not involved in immune response unlike the remainder of the examined genes. Instead, it only serves as the exclusive entry receptor to BCoV, and therefore not expected to interact with the other genes analyzed in this study. It is also expected that NANS would not appear in literature pertaining to other coronaviruses as the receptor is specific only to BCoV and the human strain hCoV-OC43 (Reguera et al., 2014). bta-miR-2325c was predicted to target NANS, and would therefore repress the expression of the receptor. This miRNA could mitigate viral entry into the host.

This study identified four miRNAs that targets multiple BCoV-associated genes. The concomitant miRNAs were miR-22, miR-2325c, miR-11975, and miR-11976. It is known that an individual miRNA can have multiple phenotypic effects, and by repressing distinct targets, it can be pleiotropic (Valastyan et al., 2009; Chen et al., 2019; Jiang et al., 2019). miRNAs with multiple targets have been shown to have several effects on the pathogenesis of a single disease including human cancer (Ventura et al., 2008; Fazi and Blandino, 2013; Baranwal et al., 2018). Select miRNAs, such as miRNA-155, control diverse functions of the immune system including regulation of inflammatory response, activation of macrophages, and immunological memory, through regulation of numerous immune response genes (Vigorito et al., 2013; Alivernini et al., 2018). In this study, we argue that a concomitant miRNA would have a greater effect on BCoV pathogenesis than a miRNA that targets only a single disease-related gene. These miRNAs could have more than one effect on the host response BCoV. The concomitant miRNAs predicted in this study may serve as effective drug targets as they repress the expression of multiple genes relating to BCoV.

miRNA prediction and functional annotation revealed the potential consequences of gene repression by the concomitant miRNAs. miR-22 targets the genes CEBPD and IRF1. CEBPD is known to enhance inflammatory and immune response genes (Ko et al., 2015). Increased expression of CEBPD also elevates the expression of IL-6 (Hungness et al., 2002). IRF1 mediates both innate and adaptive responses to viral infection such as inducing apoptosis and restricting viral replication (Kanazawa et al., 2004; Kung et al., 2014; Nair et al., 2017; Panda et al., 2019). Both target genes are highly indicated in the immune related pathways found in this study. As these genes partake in the host immune response to the virus, the expression of these genes is beneficial for the host. As miR-22 represses the expression of these advantageous genes, a drug would aim to decrease the expression of miR-22. Previously, miR-22-3p has been identified as a potential drug target and regulator in autoimmune diseases (Wang et al., 2017) and has also been proposed as a biomarker for pancreatic cancer and bovine pulmonary hypertension (Hussein et al., 2017; Grunig et al., 2018). Altering the expression of miR-22-3p has been shown to relieve symptoms of Alzheimer’s Disease, and conserved miRNAs, such as bta-mir-288 and bta-mir-193a, have been identified as drug targets for bovine disease (Ji et al., 2019; Morenikeji et al., 2019). The current study suggests that miR-22-3p may be an effective drug target or diagnostic marker for treating BCoV infection.

Our study also indicated that bta-miR-11976 and bta-miR-11975 target the genes CEBPB, SRC, and TLR9 which are all related to host immune response. Specifically, TLR9 recognizes viral nucleic acids, increases T cell and CD8 + production, induces cytokines, and enhances antibody response (Griebel et al., 2005; Lester and Li, 2014). TLR9 has been found to confer immunity against the respiratory pathogen Legionella pneumophila in mice (Bhan et al., 2008). bta-miR-11976 and bta-miR-11975 repress TLR9, and therefore their expression may decrease the host’s capacity to response to a BCoV. A drug repressing bta-miR-11976 and bta-miR-11975 may decrease the severity of BCoV infection. miRNA-2325c has not been functionally annotated nor studied to the best of our knowledge but through the in silico tools, we are able to infer the regulatory function of this miRNA. miRNA-2325c targets the genes NANS and RHOA. NANS codes for BCoV’s receptor. As bta-miR-2325c targets RHOA, it will decrease the regulation of lymphocyte and leukocyte activation, decrease the regulation of leukocyte differentiation, and affect the IL-2 pathway. RHOA has been determined to regulate the cell cycle, cytoskeleton formation, and cell morphology (Oishi et al., 2012; Satoh et al., 2012). We are hereby proposing that miRNA-2325c therefore may regulate structural components of bovine cells which might have a protective function against the virus invasion by repression of NANS that codes for BCoV receptor.

Recognition of miRNA targets in conjunction with network analysis creates an extensive visualization of the identified genes’ biological processes. Gene interaction networks have previously been used to identify important genes associated with diseases such as bovine mastitis, trypanosomosis, tuberculosis, and human cancer (Aranday-Cortes et al., 2012; Al-Aamri et al., 2019; Li et al., 2019; Morenikeji et al., 2019). Certain genes were revealed to be central in the gene interaction network of this study. One such gene, IL-6, is responsible for numerous host viral responses. For example, it affects inflammation, activates B cells, controls the differentiation of T cells, and induces antibody synthesis (Tanaka et al., 2014). While none of the concomitant miRNAs target IL-6, the miRNA miR-22-3p would affect the expression of IL-6 through its regulation of CEBPB and IRF1. IL-6 was found to be connected to IRF1 via co-occurrence and text-mining. Functional annotation of CEBPB indicates that it is associated with IL-6 expression (Ko et al., 2015; Delgobo et al., 2019). Both these genes were also found to be associated with the IL-6 signaling pathway. Expression of miR-22-3p would therefore effect more than simply its target genes as it might affect several pathways or biological processes.

Furthermore, evolutionary analysis revealed that miR-22 was perfectly conserved across all of the species it was compared to including humans, sheep, bat, and mouse indicating its functional role has been selected and preserved during evolution. Flanking regions were also perfectly conserved. This evolutionary conservation reveals the fact that mutations were not tolerated in this region, suggesting that miR-22 has an important functional role that is shared across the species it appears in.

Moreover, miR-22-3p was found to have a variant that results in a difference in the mature miRNA. This variant had a change tolerance of 3.61. The change tolerance score, which is the Genomic Evolutionary Rate Profiling (GERP) score, is defined as “the reduction in the number of substitutions in the multi-species sequence alignment compared to the neutral expectation” (Huber et al., 2020). GERP scores reflect “rejected substitutions” and quantify negative selection experienced at a particular location (Cooper et al., 2005; Pollard et al., 2010). Commonly, a cutoff of an absolute value of greater than or equal to 2 is used to determine important variants (Jha et al., 2015). With a GERP score of 3.61, the variant of miR-22-3p experienced 3.61 fewer substitutions at that location compared to the neutral rate of evolution. This high preclusion of substitutions shows that the variant is deleterious. miR-2325c revealed eight variants. Only three of them had a GERP score with an absolute value greater than two. The GERP scores were negative, indicating that substitutions occur at a rate that is greater than neutral, suggesting that substitutions are tolerated in these locations, and it is likely that they are not detrimental. The variants may even be beneficial, as previously described in which similar variation in mature miRNAs targeted CD14 gene, thus regulating its expression during bovine trypanosomosis (Morenikeji et al., 2020). Further work is needed to validate these predictions with in vitro analyses examining miRNA expression, binding, and their effect on the host immune response to bovine coronaviruses.

Conclusion

This study reveals computational evidence that particular miRNAs could regulate the host immune response to bovine coronavirus. Eleven genes were identified as predominant players in the host response, and an integrated pipeline was developed in order to determine the importance of the miRNAs that targeted them. The concomitant miRNAs were assumed to play a larger role in disease pathogenesis, as they could affect the expression of multiple BCoV-associated genes such as TLR9, IRF1, and RHOA. Functional capacities of the predicted miRNAs were inferred using annotations of the target genes, their network interactions, as well as their associated pathways. Some important pathways identified include the IL-6 signaling pathway, regulation of lymphocyte activation, and the IL-2 pathway. A critical function of miR-22-3p is suggested as the miRNA was found to be completely conserved across numerous species. miR-22-3p was also found to contain a potentially deleterious variant. Overall, this study expands our understanding of the host immune response to BCoV. Four miRNAs, bta-miR-11975, bta-miR-11976, bta-miR-22-3p, and bta-miR-2325c that were identified could serve as potential biomarkers with which to diagnose BCoV and as drug targets for therapeutic treatments. However, further work is recommended to validate these predictions with in vitro analyses examining miRNA expression, binding, and their effect on the host immune response to BCoV.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Author Contributions

OBM and BNT conceptualized and designed the experiments, revised the manuscript, and contributed to the discussion and scientific content. OBM, MW, ES, KB, and EY carried out the experiments, analyzed the data, and drafted the manuscript. All authors read and approved the final version of the manuscript.

Funding

This work was supported by the Ralph E. Hansmann, Science Students Support Fund and the Sergei S. Zlinkoff Student Medical Research Fund, Hamilton College SSRF. Page charges were graciously supported by the College of Health Sciences and Technology, Rochester Institute of Technology. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.584392/full#supplementary-material

Footnotes

  1. ^ http://geneontology.org/
  2. ^ https://www.ncbi.nlm.nih.gov/
  3. ^ http://bioinformatics.psb.ugent.be/webtools/Venn/
  4. ^ https://string-db.org/
  5. ^ https://useast.ensembl.org/Multi/Tools/Blast
  6. ^ https://itol.embl.de/tree/
  7. ^ https://ga.genecards.org/

References

Agarwal, V., Bell, G. W., Nam, J. W., and Bartel, D. P. (2015). Predicting effective microRNA target sites in mammalian mRNAs. eLife 4:e05005. doi: 10.7554/eLife.05005

PubMed Abstract | CrossRef Full Text | Google Scholar

Aich, P., Wilson, H. L., Kaushik, R. S., Potter, A. A., Babiuk, L. A., and Griebel, P. (2007). Comparative analysis of innate immune responses following infection of newborn calves with bovine rotavirus and bovine coronavirus. J. Gen. Virol. 88(Pt 10), 2749–2761. doi: 10.1099/vir.0.82861-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Aamri, A., Taha, K., Al-Hammadi, Y., Maalouf, M., and Homouz, D. (2019). Analyzing a co-occurrence gene-interaction network to identify disease-gene association. BMC Bioinformatics 20:70. doi: 10.1186/s12859-019-2634-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Alivernini, S., Gremese, E., McSharry, C., Tolusso, B., Ferraccioli, G., McInnes, I. B., et al. (2018). MicroRNA-155-at the critical interface of innate and adaptive immunity in arthritis. Front. Immunol. 8:1932. doi: 10.3389/fimmu.2017.01932

PubMed Abstract | CrossRef Full Text | Google Scholar

Amer, H. M. (2018). Bovine-like coronaviruses in domestic and wild ruminants. Anim. Health Res. Rev. 19, 113–124. doi: 10.1017/S1466252318000117

PubMed Abstract | CrossRef Full Text | Google Scholar

Aranday-Cortes, E., Hogarth, P. J., Kaveh, D. A., Whelan, A. O., Villarreal-Ramos, B., Lalvani, A., et al. (2012). Transcriptional profiling of disease-induced host responses in bovine tuberculosis and the identification of potential diagnostic biomarkers. PLoS One 7:e30626. doi: 10.1371/journal.pone.0030626

PubMed Abstract | CrossRef Full Text | Google Scholar

Athmer, J., Fehr, A. R., Grunewald, M. E., Qu, W., Wheeler, D. L., Graepel, K. W., et al. (2018). Selective packaging in murine coronavirus promotes virulence by limiting type I interferon responses. mBio 9:e00272-18. doi: 10.1128/mBio.00272-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Baranwal, S., Rawat, S. G., and Gupta, P. (2018). miR-301, pleiotropic MicroRNA in regulation of inflammatory bowel disease and colitis-associated cancer. Front. Immunol. 9:522. doi: 10.3389/fimmu.2018.00522

PubMed Abstract | CrossRef Full Text | Google Scholar

Belouzard, S., Millet, J. K., Licitra, B. N., and Whittaker, G. R. (2012). Mechanisms of coronavirus cell entry mediated by the viral spike protein. Viruses 4, 1011–1033. doi: 10.3390/v4061011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhan, U., Trujillo, G., Lyn-Kew, K., Newstead, M. W., Zeng, X., Hogaboam, C. M., et al. (2008). Toll-like receptor 9 regulates the lung macrophage phenotype and host immunity in murine pneumonia caused by Legionella pneumophila. Infect. Immun. 76, 2895–2904. doi: 10.1128/IAI.01489-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Breuer, K., Foroushani, A. K., Laird, M. R., Chen, C., Sribnaia, A., Lo, R., et al. (2012). InnateDB: systems biology of innate immunity and beyond - recent updates and continuing curation. Nucleic Acids Res. 41, D1228–D1233.

Google Scholar

Calzada-Nova, G., Schnitzlein, W., Husmann, R., and Zuckermann, F. A. (2010). Characterization of the cytokine and maturation responses of pure populations of porcine plasmacytoid dendritic cells to porcine viruses and toll-like receptor agonists. Vet. Immunol. Immunopathol. 135, 20–33. doi: 10.1016/j.vetimm.2009.10.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H., Zhang, Z., and Feng, D. (2019). Prediction and interpretation of miRNA-disease associations based on miRNA target genes using canonical correlation analysis. BMC Bioinformatics 20:404. doi: 10.1186/s12859-019-2998-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Mandel, E., Thomson, J., Wu, Q., Callis, T. E., Hammond, S. M., et al. (2006). The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat. Genet. 38, 228–233. doi: 10.1038/ng1725

PubMed Abstract | CrossRef Full Text | Google Scholar

Cooper, G. M., Stone, E. A., Asimenos, G. Nisc Comparative Sequencing Program, Green, E. D., Batzoglou, S., et al. (2005). Distribution and intensity of constraint in mammalian genomic sequence. Genome Res. 15, 901–913. doi: 10.1101/gr.3577405

PubMed Abstract | CrossRef Full Text | Google Scholar

Decaro, N., Elia, G., Campolo, M., Desario, C., Mari, V., Radogna, A., et al. (2008). Detection of bovine coronavirus using a TaqMan-based real-time RT-PCR assay. J. Virol. Methods 151, 167–171. doi: 10.1016/j.jviromet.2008.05.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Delgobo, M., Mendes, D. A., Kozlova, E., Rocha, E. L., Rodrigues-Luiz, G. F., Mascarin, L., et al. (2019). An evolutionary recent IFN/IL-6/CEBP axis is linked to monocyte expansion and tuberculosis severity in humans. eLife 8:e47013. doi: 10.7554/eLife.47013

PubMed Abstract | CrossRef Full Text | Google Scholar

Ellis, J. (2019). What is the evidence that bovine coronavirus is a biologically significant respiratory pathogen in cattle? Can. Vet. J. 60, 147–152.

Google Scholar

Fan, Y., Siklenka, K., Arora, S. K., Ribeiro, P., Kimmins, S., and Xia, J. (2016). miRNet - dissecting miRNA-target interactions and functional associations through network-based visual analysis. Nucleic Acids Res. 44, W135–W141. doi: 10.1093/nar/gkw288

PubMed Abstract | CrossRef Full Text | Google Scholar

Fazi, F., and Blandino, G. (2013). MicroRNAs: non-coding pleiotropic factors in development, cancer prevention and treatment. Microrna 2:81. doi: 10.2174/2211536611302020001

PubMed Abstract | CrossRef Full Text | Google Scholar

Fehr, A. R., and Perlman, S. (2015). Coronaviruses: an overview of their replication and pathogenesis. Meth. Mol. Biol. 1282, 1–23. doi: 10.1007/978-1-4939-2438-7_1

PubMed Abstract | CrossRef Full Text | Google Scholar

Gene Ontology Consortium (2004). The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 32(Suppl._1), D258–D261. doi: 10.1093/nar/gkh036

PubMed Abstract | CrossRef Full Text | Google Scholar

Glazov, E. A., Kongsuwan, K., Assavalapsakul, W., Horwood, P. F., Mitter, N., and Mahony, T. J. (2009). Repertoire of bovine miRNA and miRNA-like small regulatory RNAs expressed upon viral infection. PLoS One 4:e6349. doi: 10.1371/journal.pone.0006349

PubMed Abstract | CrossRef Full Text | Google Scholar

Griebel, P. J., Brownlie, R., Manuja, A., Nichani, A., Mookherjee, N., Popowych, Y., et al. (2005). Bovine toll-like receptor 9: a comparative analysis of molecular structure, function and expression. Vet. Immunol. Immunopathol. 108, 11–16. doi: 10.1016/j.vetimm.2005.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffiths-Jones, S., Kaur Saini, H. K., van Dongen, S., and Enright, A. J. (2008). miRBase: tools for microRNA genomics. Nucleic. Acids Res. 36, D154–D158. doi: 10.1093/nar/gkm952

PubMed Abstract | CrossRef Full Text | Google Scholar

Grunig, G., Eichstaedt, C. A., Verweyen, J., Durmus, N., Saxer, S., Krafsur, G., et al. (2018). Circulating MicroRNA markers for pulmonary hypertension in supervised exercise intervention and nightly oxygen intervention. Front. Phys. 9:955. doi: 10.3389/fphys.2018.00955

PubMed Abstract | CrossRef Full Text | Google Scholar

Gussow, A. B., Auslander, N., Faure, G., Wolf, Y. I., Zhang, F., and Koonin, E. V. (2020). Genomic determinants of pathogenicity in SARS-CoV-2 and other human coronaviruses. Proc. Natl. Acad. Sci. U.S.A. 117, 15193–15199. doi: 10.1073/pnas.2008176117

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanif, Q., Farooq, M., Amin, I., Mansoor, S., Zhang, Y., and Khan, Q. M. (2018). In silico identification of conserved miRNAs and their selective target gene prediction in indicine (Bos indicus) cattle. PLoS One 13:e0206154. doi: 10.1371/journal.pone.0206154

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayashida, M., and Akutsu, T. (2016). Complex network-based approaches to biomarker discovery. Biomark. Med. 10, 621–632. doi: 10.2217/bmm-2015-0047

PubMed Abstract | CrossRef Full Text | Google Scholar

Huber, C. D., Kim, B. Y., and Lohmueller, K. E. (2020). Population genetic models of GERP scores suggest pervasive turnover of constrained sites across mammalian evolution. PLoS Genet. 16:e1008827. doi: 10.1371/journal.pgen.1008827

PubMed Abstract | CrossRef Full Text | Google Scholar

Hungness, E. S., Luo, G. J., Pritts, T. A., Sun, X., Robb, B. W., Hershko, D., et al. (2002). Transcription factors C/EBP-beta and -delta regulate IL-6 production in IL-1beta-stimulated human enterocytes. J. Cell. Physiol. 192, 64–70. doi: 10.1002/jcp.10116

PubMed Abstract | CrossRef Full Text | Google Scholar

Hussein, N. A., Kholy, Z. A., Anwar, M. M., Ahmad, M. A., and Ahmad, S. M. (2017). Plasma miR-22-3p, miR-642b-3p and miR-885-5p as diagnostic biomarkers for pancreatic cancer. J. Cancer Res. Clin. Oncol. 143, 83–93. doi: 10.1007/s00432-016-2248-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Jha, P., Lu, D., and Xu, S. (2015). Natural selection and functional potentials of human noncoding elements revealed by analysis of next generation sequencing data. PLoS One 10:e0129023. doi: 10.1371/journal.pone.0129023

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, Q., Wang, X., Cai, J., Du, X., Sun, H., and Zhang, N. (2019). MiR-22-3p regulates amyloid β deposit in mice model of Alzheimer’s disease by targeting mitogen-activated protein kinase 14. Curr. Neurovasc. Res. 16, 473–480. doi: 10.2174/1567202616666191111124516

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, C., Guo, Y., Yu, H., Lu, S., and Meng, L. (2019). Pleiotropic microRNA-21 in pulmonary remodeling: novel insights for molecular mechanism and present advancements. Allergy Asthma Clin. Immunol. 15:33. doi: 10.1186/s13223-019-0345-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kameka, A. M., Haddadi, S., Kim, D. S., Cork, S. C., and Abdul-Careem, M. F. (2014). Induction of innate immune response following infectious bronchitis corona virus infection in the respiratory tract of chickens. Virology 450–451, 114–121. doi: 10.1016/j.virol.2013.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanazawa, N., Kurosaki, M., Sakamoto, N., Enomoto, N., Itsui, Y., Yamashiro, T., et al. (2004). Regulation of hepatitis C virus replication by interferon regulatory factor 1. J. Virol. 78, 9713–9720. doi: 10.1128/JVI.78.18.9713-9720.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Khurana, V., Peng, J., Chung, C. Y., Auluck, P. K., Fanning, S., Tardiff, D. F., et al. (2017). Genome-scale networks link neurodegenerative disease genes to α-synuclein through specific molecular pathways. Cell Syst. 4, 157–170.e14. doi: 10.1016/j.cels.2016.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Ko, C., Chang, W., and Wang, J. (2015). Biological roles of CCAAT/Enhancer-binding protein delta during inflammation. J. Biomed. Sci. 22:6. doi: 10.1186/s12929-014-0110-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozomara, A., and Griffiths-Jones, S. (2014). miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 42, D68–D73. doi: 10.1093/nar/gkt1181

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozomara, A., Birgaoanu, M., and Griffiths-Jones, S. (2019). miRBase: from microRNA sequences to function. Nucleic Acids Res. 47, D155–D162. doi: 10.1093/nar/gky1141

PubMed Abstract | CrossRef Full Text | Google Scholar

Krützfeldt, J., Rajewsky, N., Braich, R., Rajeev, K. G., Tuschl, T., Manoharan, M., et al. (2005). Silencing of microRNAs in vivo with ‘antagomirs’. Nature 438, 685–689. doi: 10.1038/nature04303

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096

PubMed Abstract | CrossRef Full Text | Google Scholar

Kung, H. C., Evensen, Ø., Hong, J. R., Kuo, C. Y., Tso, C. H., Ngou, F. H., et al. (2014). Interferon regulatory factor-1 (IRF-1) is involved in the induction of phosphatidylserine receptor (PSR) in response to dsRNA virus infection and contributes to apoptotic cell clearance in CHSE-214 cell. Int. J. Mol. Sci. 15, 19281–19306. doi: 10.3390/ijms151019281

PubMed Abstract | CrossRef Full Text | Google Scholar

Lawless, N., Vegh, P., O’Farrelly, C., and Lynn, D. J. (2014). The role of microRNAs in bovine infection and immunity. Front. Immunol. 5:611. doi: 10.3389/fimmu.2014.00611

PubMed Abstract | CrossRef Full Text | Google Scholar

Lester, S. N., and Li, K. (2014). Toll-like receptors in antiviral innate immunity. J. Mol. Biol. 426, 1246–1264. doi: 10.1016/j.jmb.2013.11.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Mao, L., Li, W., Hao, F., Zhong, C., Zhu, X., et al. (2018). Analysis of microRNAs expression profiles in Madin-Darby bovine kidney cells infected with caprine parainfluenza virus type 3. Front. Cell. Infect. Microbiol. 8:93. doi: 10.3389/fcimb.2018.00093

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Chen, X., and Chen, Z. (2019). Identification of key candidate genes in dairy cow in response to Escherichia coli mastitis by bioinformatical analysis. Front. Genet. 10:1251. doi: 10.3389/fgene.2019.01251

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S. W., Wang, C. Y., Jou, Y. J., Huang, S. H., Hsiao, L. H., Wan, L., et al. (2016). SARS coronavirus papain-like protease inhibits the TLR7 signaling pathway through removing Lys63-linked polyubiquitination of TRAF3 and TRAF6. Int. J. Mol. Sci. 17:678. doi: 10.3390/ijms17050678

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Liu, H., Jin, X., Lo, L., and Liu, J. (2012). Expression profiles of microRNAs from lactating and non-lactating bovine mammary glands and identification of miRNA related to lactation. BMC Genomics 13:731. doi: 10.1186/1471-2164-13-731

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, T. X., and Rothenberg, M. E. (2018). MicroRNA. J. Allergy Clin. Immunol. 141, 1202–1207. doi: 10.1016/j.jaci.2017.08.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucas, C., Wong, P., Klein, J., Castro, T. B. R., Silva, J., Sundaram, M., et al. (2020). Longitudinal analyses reveal immunological misfiring in severe COVID-19. Nature 584, 463–469. doi: 10.1038/s41586-020-2588-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Malbon, A., Meli, M. L., Barker, E. N., Davidson, A. D., Tasker, S., and Kipar, A. (2019). Inflammatory mediators in the mesenteric lymph nodes, site of a possible intermediate phase in the immune response to feline coronavirus and the pathogenesis of feline infectious peritonitis? J. Comp. Pathol. 166, 69–86. doi: 10.1016/j.jcpa.2018.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

McGruder, B., and Leibowitz, J. L. (2015). A review of genetic methods and models for analysis of coronavirus-induced severe pneumonitis. J. Gen. Virol. 96(Pt 3), 494–506. doi: 10.1099/vir.0.069732-0

PubMed Abstract | CrossRef Full Text | Google Scholar

McLaren, W., Gil, L., Hunt, S. E., Ritchie, G. R., Thormann, A., Flicek, P., et al. (2016). The ensembl variant effect predictor. Genome Biol. 17:122. doi: 10.1186/s13059-016-0974-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Mi, H., Muruganujan, A., Ebert, D., Huang, X., and Thomas, P. D. (2019). PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 47, D419–D426. doi: 10.1093/nar/gky1038

PubMed Abstract | CrossRef Full Text | Google Scholar

Mishra, S., Yadav, T., and Rani, V. (2016). Exploring miRNA-based approaches in cancer diagnostics and therapeutics. Crit. Rev. Oncol. Hematol. 98, 12–23. doi: 10.1016/j.critrevonc.2015.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Mohr, A. M., and Mott, J. L. (2015). Overview of MicroRNA biology. Semin. Liver Dis. 35, 3–11. doi: 10.1055/s-0034-1397344

PubMed Abstract | CrossRef Full Text | Google Scholar

Morenikeji, O. B., Capria, A. L., Ojurongbe, O., and Thomas, B. N. (2020). 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. Genes 11:112. doi: 10.3390/genes11010112

PubMed Abstract | CrossRef Full Text | Google Scholar

Morenikeji, O. B., Hawkes, M. E., Hudson, A. O., and Thomas, B. N. (2019). Computational network analysis identifies evolutionarily conserved miRNA gene interactions potentially regulating immune response in bovine trypanosomosis. Front. Microbiol. 10:2010. doi: 10.3389/fmicb.2019.02010

PubMed Abstract | CrossRef Full Text | Google Scholar

Mousavizadeh, L., and Ghasemi, S. (2020). Genotype and phenotype of COVID-19: their roles in pathogenesis. J. Microbiol. Immunol. Infect. doi: 10.1016/j.jmii.2020.03.022 [Epub ahead of print].

CrossRef Full Text | PubMed Abstract | Google Scholar

Nair, S., Poddar, S., Shimak, R. M., and Diamond, M. S. (2017). Interferon regulatory factor 1 protects against Chikungunya virus-induced immunopathology by restricting infection in muscle cells. J. Virol. 91:e01419-17. doi: 10.1128/JVI.01419-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Oliveira, G. B., Regitano, L., Cesar, A., Reecy, J. M., Degaki, K. Y., Poleti, M. D., et al. (2018). Integrative analysis of microRNAs and mRNAs revealed regulation of composition and metabolism in Nelore cattle. BMC Genom. 19:126. doi: 10.1186/s12864-018-4514-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Oishi, A., Makita, N., Sato, J., and Iiri, T. (2012). Regulation of RhoA signaling by the cAMP-dependent phosphorylation of RhoGDIα. J. Biol. Chem. 287, 38705–38715. doi: 10.1074/jbc.M112.401547

PubMed Abstract | CrossRef Full Text | Google Scholar

Oma, V. S., Tråvén, M., Alenius, S., Myrmel, M., and Stokstad, M. (2016). Bovine coronavirus in naturally and experimentally exposed calves; viral shedding and the potential for transmission. Virol. J. 13:100. doi: 10.1186/s12985-016-0555-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Panda, D., Gjinaj, E., Bachu, M., Squire, E., Novatt, H., Ozato, K., et al. (2019). IRF1 maintains optimal constitutive expression of antiviral genes and regulates the early antiviral response. Front. Immunol. 10:1019. doi: 10.3389/fimmu.2019.01019

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, S. S. (2019). Post-glycosylation modification of sialic acid and its role in virus pathogenesis. Vaccines 7:171. doi: 10.3390/vaccines7040171

PubMed Abstract | CrossRef Full Text | Google Scholar

Pollard, K. S., Hubisz, M. J., Rosenbloom, K. R., and Siepel, A. (2010). Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res. 20, 110–121. doi: 10.1101/gr.097857.109

PubMed Abstract | CrossRef Full Text | Google Scholar

Reguera, J., Mudgal, G., Santiago, C., and Casasnovas, J. M. (2014). A structural view of coronavirus-receptor interactions. Virus Res. 194, 3–15. doi: 10.1016/j.virusres.2014.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Saif, L. J. (2010). Bovine respiratory coronavirus. The veterinary clinics of North America. Food Anim. Pract. 26, 349–364. doi: 10.1016/j.cvfa.2010.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Saif, L. J., and Jung, K. (2020). Comparative pathogenesis of bovine and porcine respiratory coronaviruses in the animal host species and SARS-CoV-2 in humans. J. Clin. Microbiol. 58:e01355-20. doi: 10.1128/JCM.01355-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Saint Pierre, A., and Génin, E. (2014). How important are rare variants in common disease? Brief. Funct. Genomics 13, 353–361. doi: 10.1093/bfgp/elu025

PubMed Abstract | CrossRef Full Text | Google Scholar

Satoh, N., Toyohira, Y., Itoh, H., Zhang, H., Ueno, S., Tsutsui, M., et al. (2012). Stimulation of norepinephrine transporter function by fasudil, a Rho kinase inhibitor, in cultured bovine adrenal medullary cells. Naunyn Schmiedebergs Arch. Pharamcol. 385, 921–931. doi: 10.1007/s00210-012-0773-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Scheel, T. K. H., Moore, M. J., Luna, J. M., Nishiuchi, E., Fak, J., Darnell, R. B., et al. (2017). Global mapping of miRNA-target interactions in cattle (Bos taurus). Sci. Rep. 7:8190. doi: 10.1038/s41598-017-07880-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Schultze, B., and Herrler, G. (1993). Recognition of N-acetyl-9-O-acetylneuraminic acid by bovine coronavirus and hemagglutinating encephalomyelitis virus. Adv. Exp. Med. Biol. 342, 299–304. doi: 10.1007/978-1-4615-2996-5_46

CrossRef Full Text | Google Scholar

Siddell, S. G., Anderson, R., Cavanagh, D., Fujiwara, K., Klenk, H. D., Macnaughton, M. R., et al. (1983). Coronaviridae. Intervirology 20, 181–189. doi: 10.1159/000149390

PubMed Abstract | CrossRef Full Text | Google Scholar

Simonson, B., and Das, S. (2015). MicroRNA therapeutics: the next magic bullet? Mini Rev. Med. Chem. 15, 467–474. doi: 10.2174/1389557515666150324123208

PubMed Abstract | CrossRef Full Text | Google Scholar

Stecher, G., Tamura, K., and Kumar, S. (2020). Molecular Evolutionary Genetics Analysis (MEGA) for macOS. Mol. Biol. Evol. 37, 1237–1239. doi: 10.1093/molbev/msz312

PubMed Abstract | CrossRef Full Text | Google Scholar

Takiuchi, E., Stipp, D. T., Alfieri, A. F., and Alfieri, A. A. (2006). Improved detection of bovine coronavirus N gene in faeces of calves infected naturally by a semi-nested PCR assay and an internal control. J. Virol. Methods 131, 148–154. doi: 10.1016/j.jviromet.2005.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka, T., Narazaki, M., and Kishimoto, T. (2014). IL-6 in inflammation, immunity, and disease. Cold Spring Harb. Prespect. Biol. 6:a016295. doi: 10.1101/cshperspect.a016295

PubMed Abstract | CrossRef Full Text | Google Scholar

Underwood, W. J., Blauwiekel, R., Delano, M. L., Gillesby, R., Mischler, S. A., and Schoell, A. (2015). “Biology and diseases of ruminants (sheep, goats, and cattle),” in Laboratory Animal Medicine, 3rd Edn, eds J. G. Fox, L. C. Anderson, F. M. Loew, and F. W. Quimby (San Diego, CA: Academic Press), 623–625.

Google Scholar

Valastyan, S., Benaich, N., Chang, A., Reinhardt, F., and Weinberg, R. A. (2009). Concomitant suppression of three target genes can explain the impact of a microRNA on metastasis. Genes Dev. 23, 2592–2597. doi: 10.1101/gad.1832709

PubMed Abstract | CrossRef Full Text | Google Scholar

van Karnebeek, C. D., Bonafé, L., Wen, X. Y., Tarailo-Graovac, M., Balzano, S., Royer-Bertrand, B., et al. (2016). NANS-mediated synthesis of sialic acid is required for brain and skeletal development. Nat. Genet. 48, 777–784. doi: 10.1038/ng.3578

PubMed Abstract | CrossRef Full Text | Google Scholar

van Rooij, E., and Kauppinen, S. (2014). Development of microRNA therapeutics is coming of age. EMBO Mol. Med. 6, 851–864. doi: 10.15252/emmm.201100899

PubMed Abstract | CrossRef Full Text | Google Scholar

Ventura, A., Young, A. G., Winslow, M. M., Lintault, L., Meissner, A., Erkeland, S. J., et al. (2008). Targeted deletion reveals essential and overlapping functions of the miR-17 through 92 family of miRNA clusters. Cell 132, 875–886. doi: 10.1016/j.cell.2008.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Vidigal, J. A., and Ventura, A. (2015). The biological functions of miRNAs: lessons from in vivo studies. Trends Cell Biol. 25, 137–147. doi: 10.1016/j.tcb.2014.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Vigorito, E., Kohlhaas, S., Lu, D., and Leyland, R. (2013). miR-155: an ancient regulator of the immune system. Immunol. Rev. 253, 146–157. doi: 10.1111/imr.12057

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, B., Yao, Q., Xu, D., and Zhang, J. A. (2017). MicroRNA-22-3p as a novel regulator and therapeutic target for autoimmune diseases. Int. Rev. Immunol. 36, 176–181. doi: 10.1080/08830185.2017.1281272

PubMed Abstract | CrossRef Full Text | Google Scholar

Weil, M. K., and Chen, A. P. (2011). PARP inhibitor treatment in ovarian and breast cancer. Curr. Probl. Cancer 35, 7–50. doi: 10.1016/j.currproblcancer.2010.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: miRNA, bovine, coronavirus, prediction, gene regulation, immune response

Citation: Morenikeji OB, Wallace M, Strutton E, Bernard K, Yip E and Thomas BN (2020) Integrative Network Analysis of Predicted miRNA-Targets Regulating Expression of Immune Response Genes in Bovine Coronavirus Infection. Front. Genet. 11:584392. doi: 10.3389/fgene.2020.584392

Received: 17 July 2020; Accepted: 04 September 2020;
Published: 30 September 2020.

Edited by:

Shibiao Wan, St. Jude Children’s Research Hospital, United States

Reviewed by:

Aline Silva Mello Cesar, University of São Paulo, Brazil
Li Guo, Nanjing University of Posts and Telecommunications, China

Copyright © 2020 Morenikeji, Wallace, Strutton, Bernard, Yip and Thomas. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Olanrewaju B. Morenikeji, omorenik@hamilton.edu; Bolaji N. Thomas, bntsbi@rit.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.