Abstract
Introduction:
Flaviviridae comprise a group of enveloped, positive-stranded RNA viruses that are mainly transmitted through either mosquitoes or tick bites and/or contaminated blood, blood products, or other body secretions. These viruses cause diseases ranging from mild to severe and are considered important human pathogens. MicroRNAs (miRNAs) are non-coding molecules involved in growth, development, cell proliferation, protein synthesis, apoptosis, and pathogenesis. These small molecules are even being used as gene suppressors in antiviral therapeutics, inhibiting viral replication. In the current study, we used bioinformatic tools to predict a possible miRNA sequence that could be complementary to the nucleocapsid (NP) and/or capsid (CP) gene of the Flaviviridae family and provide an inhibitory solution.
Methods:
Bioinformatics is a field of science that includes tremendous computational analysis, logarithms, and sequence alignments. To predict the right alignments between miRNA and viral mRNA genomes, we used computational databases such as miRBase, NCBI, and Basic Alignment Search Tool–nucleotides (BLAST-n).
Results:
Of the 2,600 mature miRNAs, hsa-miR-548d-3p revealed complementary sequences with the flavivirus capsid gene and bovine viral diarrhea virus (BVDV) capsid gene and was selected as a possible candidate to inhibit flaviviruses.
Conclusion:
Although more detailed in vitro and in vivo studies are required to test the possible inhibitory effects of hsa-miR-548d-3p against flaviviruses, this computational study may be the first step to study further, developing a novel therapeutic for lethal viruses within the Flaviviridae family using suggested candidate miRNAs.
Introduction
Arboviruses (arthropod-borne viruses) are a group of viruses that are classified into different taxonomic families such as Flaviviridae, Bunyaviridae, Togaviridae, Rhabdoviridae, Reoviridae, and Asfarviridae, with Flaviviridae, Togaviridae, and Bunyaviridae being the families that cause disease in humans (). In this study, we primarily focused on the Flaviviridae family, which includes four species of pestiviruses, namely, bovine viral diarrhea viruses 1 and 2 (BVDV 1 and BVDV2), classical swine fever virus (CSFV), and border disease virus (BDV) (Mari et al., 2016; Maurer et al., 2004; Warrener and Collett, 1995; Schweize and Peterhans, 2001). The Flavivirus genus also includes global human pathogens such as Zika virus (ZIKV), West Nile virus (WNV), Japanese encephalitis virus (JENV), dengue virus (DENV), yellow fever virus (YFV), and tick-borne encephalitis virus (TBEV), which all pose a threat to global public health (; Reed et al., 1998). Hepacivirus, another member of the flavivirus family, or hepatitis C virus or simply HCV, is responsible for non-A and non-B hepatitis among humans (; Merwaiss et al., 2019; Suzich et al., 1993). This genus also includes additional 50 arthropod-borne viruses, which are mainly transmitted via mosquito bites and tick bites () (Figure 1). Although the primary vectors are mainly mosquitos and ticks, these viruses have also been detected and isolated from bats and rodents ().
FIGURE 1
The incubation period of flavivirus infections in humans can range from 3 to 6 days (), with presentation of acute flavivirus diseases ranging from being mild to severe and can be life-threatening (Pierson and Diamond, 2020). Pierson and Diamond mentioned that the symptoms of mild illness are mostly similar to flu-like symptoms, which includes asymptomatic infection and/or self-limiting febrile episodes, while severe illness includes hemorrhagic fever, shock syndrome, encephalitis, paralysis, congenital defects, hepatitis, and hepatic failure (2020; ). Although there are vaccines currently available for most of these viruses, which have also been successful, however, due to re-establishment of vectors, globalization, and urbanizations, epidemics continue to occur, which restricts effectiveness of these vaccines (; ; van Leur et al., 2021). People infected with DENV can develop more severe manifestations like dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS), which includes vascular leakage or hypovolemic shock and coagulopathy, followed by bleeding, organ impartment, and death (Simmons et al., 2012; ).
Of these flaviviruses, WNV and JENV are known as neurotropic viruses and cause acute encephalopathy, causing severe neuroinflammation of the central nervous system (CNS) and the blood–brain barrier (). In WNV, symptoms include flaccid paralysis, convulsions, cranial neuropathies, optic neuritis, ataxia, stiffness, rigidity spasms, and tremors that might cause long-term neurological changes (World Health Organization, 2019). JENV shows symptoms similar to that of WNV but is rare and has a much higher fatality rate of 30% (World Health Organization, 2019). TBEV, is also a member of the encephalitis virus family like WNV and JENV, but on the contrary, it is not transmitted by mosquitos like the other members of the arbovirus family, rather it is transmitted by infected tick (Ixodes ricinus) bites that can spread from animals to humans (Turtle et al., 2012).
HCV, another member of flaviviruses, attacks the hepatocytes (liver cells) in humans (Song et al., 2001), causing inflammation of the liver. HCV is a blood-borne virus and is transmitted primarily through infected blood and/or blood products or contaminated body fluids. One example of possible HCV transmission is the sharing of needles among drug abusers who use needle injections. In the mid 2000s, HCV transmission has also occurred among men who have sexual encounters with other men, also known as Men sex with other men (MSM) (Nijmeijer et al., 2019) (Figure 2). Despite treatments currently available, there is no vaccine for HCV (). When left untreated, HCV can lead to liver cirrhosis and chronic hepatitis C infection, leading to liver carcinoma (). Progression of HCV is rather slow and can remain unnoticed (asymptomatic) for decades until the patient develops liver disease, which results in delay in diagnosis and treatment ().
FIGURE 2
YFV infections occur in 12% individuals with a 95% confidence interval and in 5 to 26% individuals with manifestations of jaundice, hemorrhage, and organ failure (Waggoneret al., 2018). Mosquitoes are primary carriers in areas of endemicity and are mainly recorded in Africa and South America, and despite successful vaccinations, outbreaks continue and lead to significant high morbidity and mortality rates (). Julander et al. stated that despite vaccinations, there is a great need for more therapies as there is no antiviral agent available for YFV (2009).
Similar to all flaviviruses, ZIKV is another member, which is transmitted to humans by Aedes (Stegomyia subgenus) mosquitoes (), and the disease caused by ZIKV can range from mild to severe, with a 3–12-day incubation period (; Musso and Gubler, 2016). Basarab et al. stated that ZIKV symptoms can include fever, conjunctivitis, arthralgia, myalgia, and itchy rashes (Musso and Gubler, 2016; Musso and Nhan, 2015; ). However, Basarab et al. claimed that symptoms also include headache, retro-orbital pain, peripheral edema, joint pain, and even gastrointestinal disturbances (2016).
Overall, flaviviruses are small positive-sense, single-stranded RNA viruses that harbor structural proteins such as the capsid (C), which is responsible for protecting the viral genome; the pre-membrane protein (prM); the envelope protein (E); and non-structural (NS) proteins that are categorized as NS2A, NS2B, NS3, NS4A, 2K, NA4B, and NS5 (Mutebi et al., 2004) with a genome of approximately 11 kb (). During viral entry, replication occurs in the endoplasmic reticulum, where ribosomes are present (Figure 3). Because the genome of these viruses can act as a messenger RNA (mRNA), the genome is readily translated into proteins, making more virus particles (van den Elsen et al., 2021). Viral attachment is accomplished by the E protein attachment to the cognate receptors (). Laureti et al. conferred that the E protein binds to receptors such as glycosaminoglycans that increase the viral density on the host cell surface, allowing for more effective receptor binding (2018; Perera-Lecoin et al., 2013). On the surface of the E protein, the ectodomain harbors three domains, namely, E-D1, E-2, and E-D3, where E-D3 interacts with attachment factors and receptors and is mainly the target of neutralizing antibodies (; Pierson Kielian, 2013).
FIGURE 3
BVDV is a causative agent of bovine diarrhea and mucosal disease and hemorrhagic syndrome with high mortality among cattle (). The virion size ranges between 40 and 60 nm (), and the genome is approximately 13.3 kb in size (Murrayet al., 2008). The viral proteins of BVDV are organized in the following order: NH2-Npro-C-Erns-E1-E2-p7-NS2- NS3-NS4A-NS4B-NS5A-NS5B-COOH (Tellinghuisen et al., 2006; Neill, 2013; ; ), which is very similar to that of flavivirus polyprotein.
Transmission of BVDV among cattle includes fomites, such as contaminated feed, water, and equipment, and among other surfaces such as the nose; tongue; milk bottle nipples; needles; palpitations; secretions; and excretion of urine feces, mucus, milk, and other contaminated materials (Niskanen et al., 2000) (Figure 4). When cattle are exposed, they usually recover over time and shed the virus temporarily; however, pregnant cattle are more susceptible, and the outcome depends on the gestational stage of the fetus (. Although cows are the main host, BVDV infects various cattle, including bisons, and can cause immune dysfunction and result in asymptomatic infections and seroconversion, including fatal mucosal disease (Hause et al., 2021).
FIGURE 4
Diseases associated with BVDV can range from clinically inappropriate to severe, even with the availability of vaccines (Xue et al., 2009). Acute and persistent BVDV infections among pregnant cows are often accompanied by transmission into the fetuses, resulting in abortions, teratogenic changes, or delivery of persistently infected, immunotolerant calves, depending on the gestation period (). In the transmission process, if a cow is pregnant and is infected with the virus, the virus is transmitted to the fetus (). The virus has the ability to cause transplacental infection, resulting in different outcomes depending on the stage, which includes fetal death, malformation, acute syndromes of the neonate, immune tolerance, and lifelong viral persistence (Peterhans et al., 2003).
In the 90s, two small RNAs were discovered in Caenorhabditis elegans (C. elegans); it was later identified that the longer RNA, about 70 nucleotides, was the precursor of shorter RNAs that were about 22 nucleotides, which were classified as microRNAs (miRNAs) due to their short length (). miRNAs are small non-coding segments of RNAs that, unlike mRNAs, which encodes proteins, control various levels of important roles such as animal and human growth regulation, development, gene expression, cell proliferation, apoptosis, and even serves as an initiator for protein synthesis (; Ranganathan and Sivasankar, 2014; ; ). Most miRNAs are transcribed from DNA sequences into primary miRNAs (pri-miRNAs), then processed into precursor miRNAs (pre-miRNAs), and then into mature miRNAs (; O'Brien et al., 2018).
There are three distinct types of miRNAs: small interference RNA (siRNA), RNA interferences (RNAi), and miRNAs (Qian et al., 2022). These molecules not only regulate gene expression or growth and development but can also suppress viral replication by targeting specific genes, resulting in inhibiting viral growth in its host. In the process of suppressing viral replication, mature miRNAs bind to complementary sequences on the 3′ end of the target mRNAs (Skalsky ans Cullen, 2010). Perfect complementarity miRNAs generally lead to potential cleavage of the mRNA genome, while imperfect complementarity results in repression and destabilization or degradation (Skalsky and Cullen, 2010; ; Selbach et al., 2008). In this study, we utilize advance bioinformatics tools to identify a possible complementary miRNA sequence to the nucleocapsid (NP) and/or capsid (CP) gene sequences of the flavivirus family.
Methods
Collection of Flavivirus genome sequences from NCBI
Complementary alignments were carried out using viral genome sequences that are responsible for the nucleocapsid and capsid protein synthesis of BVDV and all flaviviruses and were obtained from the National Center for Biotechnology Information (NCBI) database (Table 1). Figure 5 shows the flowchart of the computational analysis and multiple sequence alignments using miRBase, NCBI, and Basic Local Alignment Search Tool–nucleotides (BLAST-n).
TABLE 1
| Virus name | Accession number | Gene name | Sequence |
|---|---|---|---|
| BVDV | AJ715397.1 | NP | TCCGACACAAAAGATAAAGGGGCGGTGAGAAAGGAGCAGCAAAAGCCAGATAGGTTGGAAAAGGGGAGAATGAAGATAACACCTAAAGAGTCAGAGAAAGACAGTAAGACCAGGCCACCAGATGCCACAATAGTGGTCGATGGAGTCAAATATCAGGTAAAGAAAAAAGGGAAAGTCAAGAGCAAGAACACCCAGGACGGCTTGTACCACAACAAAAATAAACCTCAAGAGTCTCGCAAGAAACTAGAGAAAGCCCTACTGGCCTGGGCAGTAATAGCCCTGGTTTTGTTTCAAGTC |
| HCV 1a Strain THCM-NR1/03 capsid protein gene | GQ913857.1 | Capsid/core protein | ATGAGCACGAATCCTAAACCTCAAAGAAAAACCAAACGTAACACCAACCGTCGCCCACAGGACGTTAAGTTCCCGGGTGGCGGTCAGATCGTTGGTGGAGTTTACTTGTTGCCGCGCAGGGGCCCTAGATTGGGTGTGCGCGCGACGAGGAAGACTTCCGAGCGGTCGCAACCTCGAGGTAGACGTCAGCCTATCCCCAAGGCGCGTCGGCCCGAGGGCAGGACCTGGGCTCAGCCCGGGTACCCTTGGCCCCTCTATGGCAATGAGGGCTGCGGGTGGGCGGGATGGCTCCTGTCCCCCCGTGGCTCTCGGCCTAGCTGGGGCCCCACAGACCCCCGGCGTAGGTCGCGCAATTTGGGTAAGGTCATCGATACCCTCACGTGCGGCTTCGCCGACCTCATGGGGTACATCCCGCTCGTCGGCGCCCCTCTTGGAGGCGCCGCCAGGGCCCTGGCGCATGGCGTCCGGGTTCTGGAAGACGGCGTGAACTATGCAACAGGGAACCTTCCTGGTTGCTCTTTTTCTATCTTCCTTCTAGCCCTGCTCTCTTGCCTGACTGTGCCCGCGTCAGCC |
| HCV Genotype 2 isolate MOR34 | JN055424.1 | Capsid/core protein | ATGAGCACGAATCCTAAACCTCAAAGACAAACCAAAAGAAACACCAACCGCCGCCCAAAGGACGTTAAGTTCCCGGGCGGTGGTCAGATCGTTGGCGGGGTGTACTTGTTGCCGCGCAGGGGCCCTAGATTGGGTGTGCGCGCGACGAGGAAGACCTCGGAGCGATCCCAGCCGCGTGAAAGGCGCCAACCCATCCCCAGGGCTGGGCGCACCACCGGCAGGTCCTGGCAGCAGCCGGGATATCCTTGGCCCCTTTATGGGAACGA |
| HCV subtype 3a isolate THCM-L3/03 | HM042020.1 | Capsid/core protein | ATGAGCACACTTCCTAAACCTCAAAGAAAAACCAAAAGAAACACCATCCGTCGCCCACAGGACGTCAAGTTCCCGGGTGGCGGACAGATCGTTGGTGGAGTATACGTGTTGCCGCGCAGGGGCCCACGATTGGGTGTGCGCGCGACGCGTAAAACTTCTGAACGGTCACAGCCTCGCGGACGACGACAGCCTATCCCCAAGGCACGTCGGAGCGAAGGCCGGTCCTGGGCTCAGCCTGGGTACCCTTGGCCCCTCTATGGTAACGAGGGCTGCGGGTGGGCAGGGTGGCTCCTGTCCCCACGCGGCTCCCGTCCATCTTGGGGCCCAAACGACCCCCGGCGACGGTCCCGCAATTTGGGTAAAGTCATCGATACCCTCACATGCGGATTCGCCGACCTCATGGGGTACATCCCGCTCGTCGGCGCTCCCGTAGGGGGCGTCGCAAGGGCCCTCGCACATGGCGTGAGGGCCCTTGAAGACGGGATAAATTTCGCAACAGGGAACTTGCCCGGTTGCTCCTTTTCTATCTTTCTTCTTGCTCTACTCTCTTGCTTAATCCATCCAGCAGCTAGC |
| HCV type 4 isolate QC27 | U33436.1 | Capsid/core protein | ATGAGCACGAATCCTAAACCTCAAAGAAAAACCAAACATAACACCAACCGCCGCCCCATGGACGTCAAGTTCCCGGGTGGTGGTCAGATCGTTGGCGGAGTTTACTTGTTGCCGCGCAGGGGCCCTCGTTTGGGTGTGCGCGCGACTCGGAAGACTTCGGAGCGGTCGCAACCTCGTGGGAGACGCCAGCCTATCCCCAAGGCGCGTCGATCCGAGGGAAGGTCCTGGGCACAGCCAGGATACCCATGGCCTCTTTACGGTAATGAGGGTTGCGGGTGGGCAGGATGGCTCCTGTCCCCCCGTGGTTCTCGACCGTCTTGGGGTCCAAATGATCCCCGGCGGAGGTCCCGCAACTTGGGTAAGGTCATCGATACCCTAACCTGCGGCTTCGCCGACCTCATGGGATACATCCCGCTCGTGGGCGCCCCCGTTGGTGGCGTCGCCAGGGCCCTGGCACATGGTGTCAGGGCCGTGGAGGACGGGATTAATTACGCAACAGGGAACCTTCCCGGTTGCTCCTTTTCTATCTTCCTCCTAGCACTCTTTTCGTGCCTGACTGTCCCCGCTTCGGCC |
| HCV type 5 isolates QC21 | U33434.1 | Capsid/core protein | ATGAGCACGAATCCTAAACCTCAAAGAAAAACCAAAAGAAACACCAACCGCCGCCCACAGGACGTCAAGTTCCCGGGCGGTGGTCAGATCGTTGGTGGAGTTTACTTGTTGCCGCGCAGGGGCCCTAGGTTGGGTGTGCGCGCGACTCGGAAGACTTCAGAACGGTCGCAACCCCGCGGACGGCGTCAGCCTATTCCCAAGGCGCGCCAATCCGCGGGCCGGTCCTGGGGTCAACCCGGGTACCCTTGGCCCCTTTATGGCAATGAGGGCCTCGGATGGGCAGGGTGGTTGCTCTCCCCCCGGGGTTCTCGGCCTAGTTGGGGCCCCAATGACCCCCGGCGAAAGTCACGTAATTTGGGTAAGGTCATCGATACCCTAACGTGCGGATTCGCCGACCTCATGGGGTATATCCCGCTCGTAGGCGGCCCCGTAGGGGGCGTCGCAAGGGCTCTCGCGCATGGTGTGAGGGTTCTTGAGGACGGGGTAAACTATGCGACAGGGAATTTGCCCGGTTGCTCTTTCTCTATCTTCCTCCTTGCACTTCTCTCGTGCTTGACTGTCCCGGCCTCTGCA |
| HCV type 6 isolate QC26 | U33435.1 | Capsid/core protein | ATGAGCACACTTCCAAAACCCCAAAGAAAAACCAAAAGAAACACCAACCGTCGCCCAATGGACGTCAAGTTCCCGGGTGGCGGTCAGATCGTTGGCGGAGTTTACTTGTTGCCGCGCAGGGGCCCACGGTTGGGTGTGCGCGCGACGAGAAAGACTTCCGAGCGATCCCAGCCCAGAGGTAGGCGTCAACCTATACCAAAAGCACGCCAGCCTCAGGGCAGGCACTGGGCTCAGCCCGGATACCCTTGGCCTCTTTATGGAAACGAGGGCTGCGGGTGGGCGGGATGGCTCTTGTCCCCCCGCGGTTCCCGGCCACATTGGGGCCCCAATGACCCCCGGCGTCGATCCCGCAATTTGGGTAAGGTCATCGATACCCTAACGTGTGGATTCGCCGATCTCATGGGGTACATTCCCGTCGTGGGCGCGCCTCTAGGCGGCGTCGCGGCTGCGCTCGCACATGGTGTGAGGGCAATCGAGGACGGGATCAATTATGCAACAGGGAATCTTCCCGGTTGCTCTTTCTCTATCTTCCTTTTGGCACTATTCTCGTGCCTCACGACGCCAGCCTCGGCC |
| DENV | KM519590.1 | Capsid protein | TTCTCAACCGGACTTTTTTCTGGGAAAGGACCCTTACGGATGGTGCTAGCATTCATCACGTTTTTGCGAGTCCTTTCCATCCCACCAACAGCAGGGATTCTGAAAAGATGGGGACAGTTGAAGAAAAATAAGGCCATCAGGATACTGATTGGATTCAGGAAGGAGATAGGCCGCATGCTGAACATCTTGAACGGGAGAAAAAGGTCAACGATAACATTGCTGTGCTTGATTCCCACCGTAATGGCGTTTCACTT |
| DENV 1 | KY346993.1 | Capsid protein | ATGAACAACCAACGGAAAAAGACGGGTCGACCGTCTTTCAATATGCTGAAACGCGCGAGAAACCGCGTGTCAACTGGTTCACAGTTGGCGAAGAGATTCTCAAAAGGATTGCTTTCAGGCCAAGGACCCATGAAATTGGTGATGGCTTTCATAGCATTTCTAAGATTTCTAGCCATACCCCCAACAGCAGGAATTTTGGCTAGATGGAGCTCATTCAAGAAGAATGGAGCGATCAAAGTGTTACGGGGTTTCAAAAAAGAGATCTCAAGCATGTTGAACATAATGGATAGAAGGAAAAGA |
| DENV 2 | JQ846016.1 | Capsid protein | ATGAATAACCAACGGAAAAAGGCGAAAAACACGCCTTTCAATATGCTGAAACGCGAGAGAAACCGCGTGTCGACTGTGCAACAGCTGACAAAGAGATTCTCACTTGGAATGCTGCAGGGACGAGGACCATTAAAACTGTTCATGGCCCTGGTGGCGTTCCTTCGTTTCCTAACAATCCCACCAACAGCAGGGATATTGAAGAGATGGGGAACAATTAAAAAATCAAAAGCTATTAATGTTTTGAGAGGGTTCAGGAAAGAGATTGGAAGGATGCTGAACATCTTGAATAGGAGACGCAGATCTGCAGGCATGATCATTATGCTGATTCCAACAGTGATGGCG |
| DENV 3 | HQ223036.1 | Capsid protein | ATGAACAACCAACGAAAAAAGACGGGAAAACCGTATATCAATATGCTGAAACGCGTGAGAAACCCTGTGTCAACTGGATCACAGTTGGCGAAGAGATTCTCAAGAGGATTGCTCAACGGCCAAGGACCAATGAAATTGGTTATGGCGTTCATAGCTTTCAGATTTCTAGCCATTCCACCAACAGCAGGAGTCTTGGCTAGATGGGGAACCTTCAAGAAGTCAGGGGCTATTAAGGTCCTAAAAGGCTTCAAGAAGGAGATTTCAAACATGCTGAGTATTATCAACAAAAGGAAAAAGACATCGCTCTGTCTCATGATGATGTTACCAGCAACACTTGCT |
| DENV 4 | GQ890685.1 | Capsid protein | TTGGTGAAGAGATTCTCAACCGGACTTTTCTCTGGGAAAGGAACCTTACGGATGGTGCTAGCATTCATCACGTTTTTGCGAGTCCTTTCCATCCCGCCAACAGCAGGGATTTTGAAAAGATGGGGACAGTTGAAAAAGAATAAGGCCATCAAGATACTGATTGGATTCAGGAAGGAGATAGGTCGCATGTTAAACATCTTAAATAGGAGAAGAAGGTCAACAATGACATTGCTGTGTTTGATTCCCACCGTAATGGCATTTCACCTGTCAACAAGAGACGGCGAACCCCTCATGATAGTGGCAAAACACGAAAGGGGGAGACCTCTCTTGTTTAAGACAACAGAAGA |
| JENV | KJ420596.1 | Capsid protein | ATCAATATGCTGAAACGCGGCATACCCCGCGTATCCCCACTTGTGGGGGTGAAGAGGGTAATTATGAACTTGCTCGACGGCAGAGGGCCAATACGATTCGTTTTGGCTCTCTTGGCGTTTTTCAAGTTCACAGCACTAGCCCCGACCAAGGCACTCGTTAGCCGATGGAAGGCAGTAGAGAAGAGCGTTGCAATGAAACATCTCACCAGTTTCAAACGAGAACTTGG |
| YFV | L06480.1 | Capsid protein | TCTGGTCGTAAAGCTCAGGGAAAAACCCTGGGCGTCAATATGGTACGACGAGGAGTTCGCTCCTTGTCAAACAAAATAAAACAAAAAACAAAACAAATTGGAAACAGACCTGGACCTTCAAGAGGTGTTCAAGGATTTATCTTTTTCTTTTTGTTCAACATTTTGACTGGAAAAAAGATCACAGCCCACCTAAAGAGGTTGTGGAAAATGCTGGACCCAAGACAAGGCTTGGCTGTTCTAAGGAAAGTTAAGAGAGTGGTGGCCAGTTTAATGAGAGGATTGTCCTCAAGGAAACGCCGTTCCCATGATGTTCTGACTGTGCAATTCCTAATTTTGGGAATGCTGTTGATGACGGGTGGA |
| WENV | FJ425728.1 | Capsid protein | TAACAACAATTAACACAGTGCGAGCTGTTTCTTAGCACGAAGATCTCGATGTCTAAGAAACCAGGAGGGCCCGGTAAAAACCGGGCTGTCAATATGCTAAAACGCGGTATGCCCCGCGGATTGTCCTTGATAGGACTAAAGAGGGCTATGCTGAGTCTGATTGACGGGAAGGGCCCAATACGTTTCGTGTTGGCTCTTTTGGCGTTT |
| ZIKV | KX443145.1 | Capsid protein | TGACAGTTCGAGTTTGAAGCGAAAGCTAGCAACAGTATCAACAGGTTTTATTTTGGATTTGGAAACGAGAGTTTCTGGTCATGAAAAACCCAAAAAAGAAATCCGGAGGATTCCGGATTGTCAATATGCTAAAACGCGGAGTAGCCCGTGTGAGCCCCTTTGGGGGCTTGAAGAGGCTGCCAGCCGGACTTCTGCTGGGTCATGGGCCCATCAGGATGGTCTTGGCGATTCTAGCCTTTTTGAGATTCACGGCAATCAAGCCATCACTGGGTCTCATCAATAGATGGGGTTCAGTGGGGAAAAAAGAGGCTATGGAAATAATAAAGAAGTTCAAGAAAGATCTGGCTGCCATACTGAGAATAATCAATGCTAGGAAGGAGAAGAAGAGACGAGGCGCAGATACTAGTGTCGGAATTGTTGGCCTCCTGCTGACCACAGCTATGGCAGCGGAGGTCACTAGACGTGGGAGTGCATACTATATGTACTTGGACAGAAACGATGCTGGGGAGGCCATATCTTTTCCAACCACATTGGGGATGAATAAGTGTTATATA |
| TBEV | EU715176.1 | Capsid protein | AAATTTTATTACACGCCCAGGGGTTTGCTTCAGACACCAACAGGAGGGCCAGGTTCGGAAGAAACAATCTTTGGTTACTACTAGTCGTGAACGTGTTGAGAAAAAGACAGCTTAGGAGAACAAGAGCTGGGGATGGCCAGGAAGGCCATTCTGAAAGGAAAGGGGGGGGGTCCCCCTCGACAGGTGTCGAAAGGACCCAAAAAGC |
Name of the virus genome, with the accession number, name of genome, and the sequence that were retrieved from the NCBI database.
FIGURE 5
Collection of miRNAs from miRBase and sequence alignments
Figures 8, 6 show the method of predicting the right miRNA sequence using miRBase (https://mirbase.org/) and NCBI BLAST–n) (https://blast.ncbi.nlm.nih.gov/Blast.cgi) (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastn&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome). For launching alignments, the miRNA sequences were entered in the Query section of BLAST and the viral mRNA sequence was entered in the Subject section.
FIGURE 6
Results
After running series of alignments, our results revealed that hsa-miR-548d-3p (MI0003668) showed complementary sequence structures with the viral genome sequences that are responsible for the nucleocapsid gene of the BVDV and flavivirus capsid protein synthesis. Table 2 shows the details about miR-548d-3p, which includes, name, species, accession number, tissues, sequence, and website. According to BLAST-n, hsa-miR-548d-3p exhibited 100% similarities and showed the highest numbers of alignment positions on the YFV capsid gene (7 locations), as compared to the BVDV NP gene and ZIKV CP gene (Figure 7; 5 locations), and 4 alignment positions on the WNV capsid gene. Figures 8A, B shows the number of alignment positions of hsa-miR-548d-3p on our virus' genome. Table 3 shows all alignment data exhibited by hsa-miR-548d-3p on viral genomes of DENV (Figures 9), HCV, and the other flavivirus members, and Figures 10–16 show the number alignment locations of hsa-miR-548d-3p on the flavivirus family's genome.
TABLE 2
| Name of miRNA | Species | Accession # | Tissue | Sequence | Website |
|---|---|---|---|---|---|
| hsa-miR-548-3p | Homo sapiens (Human) | MIMAT0003323 | Melanoblast | CAAAAACCACAGUUUCUUUUGC | https://mirbase.org/mature/MIMAT0003323?mature_acc=MIMAT0003323 |
Details about the hsa-miR-548-3p sequence.
FIGURE 7
FIGURE 8
TABLE 3
| Virus genome | miRNA name and sequence | # Of matches | Score | E-Value | Identities | Query location | Subject location |
|---|---|---|---|---|---|---|---|
| BVDV Nucleocapsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 5 | 14.4 bits (7) | 0.22 | 7/7 (100%) | Query 11–18 | Subject 245–238 |
| 14.4 bits (7) | 0.22 | 7/7 (100%) | Query 15–21 | Subject 14–8 | |||
| 14.4 bits (7) | 0.22 | 7/7 (100%) | Query 16–22 | Subject 45–39 | |||
| 14.4 bits (7) | 0.22 | 7/7 (100%) | Query 13–19 | Subject 165–159 | |||
| 14.4 bits (7) | 0.22 | 7/7 (100%) | Query 3–9 | Subject 287–281 | |||
| HCV Genotype 1a Capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 2 | 16.3 bits (8) | 0.11 | 8/8 (100%) | Query 2–9 | Subject 27–34 |
| 14.4 bits (7) | 0.43 | 7/7 (100%) | Query 13–19 | Subject 29–23 | |||
| HCV Genotype 2 capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 1 | 20.3 bits (10) | 0.003 | 10/10 (100%) | Query 12–21 | Subject 42–33 |
| HCV Genotype 3 capsid protein | has-miR-548d-3p CAAAAACCACAGUUUCUUUUGC0.43 | 3 | 20.3 bits (10) | 0.007 | 10/10 (100%) | Query 12–21 | Subject 42–33 |
| 16.4 bits (8) | 0.11 | 8/8 (100%) | Query 2–9 | Subject 27–34 | |||
| 14.4 bits (7) | 0.43 | 7/7 (100%) | Query 13–19 | Subject 29–23 | |||
| HCV Genotype 4 capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 2 | 16.4 bits (8) | 0.11 | 8/8 (100%) | Query 2–9 | Subject 27–34 |
| 14.4 bits (7) | 0.43 | 7/7 (100%) | Query 13–9 | Subject 29–23 | |||
| HCV Genotype 5 capsid Protein | has-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 3 | 20.3 bits (10) | 0.007 | 10/10 (100%) | Query 12–21 | Subject 42–33 |
| 16.4 bits (8) | 0.11 | 8/8 (100%) | Query 2–9 | Subject 27–34 | |||
| 14.4 bits (7) | 0.43 | 7/7 (100%) | Query 13–19 | Subject 29–23 | |||
| HCV Genotype 6 capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 2 | 16.4 bits (8) | 0.11 | 8/8 (100%) | Query 2–9 | Subject 27–34 |
| 14.4 bits (7) | 0.43 | 7/7 (100%) | Query 13–19 | Subject 29–23 | |||
| Dengue virus capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 1 | 14.4 bits (7) | 0.19 | 7/7 (100%) | Query 1–7 | Subject 66–60 |
| Dengue virus 1 capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 1 | 14.4 bits (7) | 0.22 | 7/7 (100%) | Query 8–14 | Subject 80–86 |
| Dengue Virus 2 capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 1 | 14.4 bits (7) | 0.25 | 7/7 (100%) | Query 9–15 | Subject 139–133 |
| Dengue Virus 3 capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 1 | 14.4 bits (7) | 0.25 | 7/7 (100%) | Query 8–14 | Subject 80–86 |
| Dengue Virus 4 capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 2 | 14.4 bits (7) | 0.26 | 7/7 (100%) | Query 1–7 | Subject 78–72 |
| 14.4 bits (7) | 0.26 | 7/7 (100%) | Query 14–20 | Subject 140–134 | |||
| Japanese encephalitis virus capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 1 | 14.4 bits (7) | 0.18 | 7/7 (100%) | Query 10–16 | Subject 207–213 |
| Yellow fever virus capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 7 | 16.4 bits (8) | 0.42 | 8/8 (100%) | Query 13–20 | Subject 180–187 |
| 14.4 bits (7) | 1.7 | 7/7 (100%) | Query 2–8 | Subject 57–63 | |||
| 14.4 bits (7) | 1.7 | 7/7 (100%) | Query 3–9 | Subject 686–692 | |||
| 14.4 bits (7) | 1.7 | 7/7 (100%) | Query 15–21 | Subject 757–751 | |||
| 14.4 bits (7) | 1.7 | 7/7 (100%) | Query 12–18 | Subject 1729–1735 | |||
| 14.4 bits (7) | 1.7 | 7/7 (100%) | Query 12–18 | Subject 2,298–2,292 | |||
| 14.4 bits (7) | 1.7 | 7/7 (100%) | Query 1–7 | Subject 2,349–2,343 | |||
| West Nile virus capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 4 | 14.4 bits (7) | 0.16 | 7/7 (100%) | Query 12–18 | Subject 27–33 |
| 14.4 bits (7) | 0.16 | 7/7 (100%) | Query 12–18 | Subject 61–55 | |||
| 14.4 bits (7) | 0.16 | 7/7 (100%) | Query 2–8 | Subject 76–82 | |||
| 14.4 bits (7) | 0.16 | 7/7 (100%) | Query 15–21 | Subject 195–201 | |||
| Zika virus capsid protein | hsa-miR-548d-3p CAAAAACCACAGUUUCUUUUGC | 5 | 16.4 bits (8) | 0.10 | 8/8 (100%) | Query 13–20 | Subject 101–94 |
| 14.4 bits (7) | 0.41 | 7/7 (100%) | Query 11–17 | Subject 70–76 | |||
| 14.4 bits (7) | 0.41 | 7/7 (100%) | Query 2–8 | Subject 84–90 | |||
| 14.4 bits (7) | 0.41 | 7/7 (100%) | Query 6–12 | Subject 432–439 | |||
| 14.4 bits (7) | 0.41 | 7/7 (100%) | Query 5–11 | Subject 524–530 |
All the alignment data on hsa-miR-548-3p to the capsid genome of all members of the flaviviruses retrieved from NCBI BLASTn, including the name of the genome and the name of the miRNA hsa-miR-548d-3p.
This Table also includes the number (#) of complementary matches between the miRNA and the viral mRNA genome sequence. The alignment score shows a significantly high score, which indicates a high degree of similar alignments between our miRNA sequence and our virus genome sequence. This table includes the E-value that measures the number of alignments similar found by hsa-miR-548d-3p by chance. The 100% of identity states that the similarity of our sequence alignments has a perfect match. These results included the query locations, which is the location on our miRNA alignment that has the perfect math, while the subject location, which is our mRNA alignment location on our viral genome sequence. On our subject, we see that our miRNA candidate aligns at multiple locations ranging from 14 to 2,349 locations.
FIGURE 9
FIGURE 10
FIGURE 11
FIGURE 12
FIGURE 13
FIGURE 14
FIGURE 15
FIGURE 16
hsa-miR-548d-3p aligned with dengue virus
Our results showed that the hsa-miR-548d-3p sequence is identical to the capsid sequence of dengue virus (DENV) and its four serotypes (Figures 9, 10; Table 1). We found that miR-548-3p has a 100% perfect match on DENV virus genome sequences, as displayed in Table 3. Figures 12A, B show that out of the four serotypes of DENV, 1 out of 5 (20%) showed 2 alignment locations and aligns at 2 nucleotides on the DENV 4 genome at subject locations 78–72 and 140–134.
hsa-miR-548d-3p alignment with HCV
The alignment analysis also revealed multiple sequence alignments between hsa-miR-548d-3p and the capsid sequence of the HCV virus and its genotypes. Figure 15 shows that out of the six genotypes of HCV, the highest number of alignment locations is observed between our miRNA and HCV genotype 3 and genotype 5 (33.33%) and the second highest number of locations is observed on genotype 1, genotype 4, and genotype 6 (50%).
hsa-miR-548d-3p aligned with Flavivirus capsid gene
Additional alignments were carried out, and it showcased that miR-548d-3p harbors sequence similarities with the capsid genome of the remaining Flavivirus members [ZIKV (Figure 7), YFV, WNV, and TBEV]. Figure 15 shows that our candidate miRNA shows similarities on various sections on the flavivirus capsid genome, and Figures 16A, B demonstrate that hsa-miR-548d-3p has the highest number of alignments on the YFV capsid genome (7 alignment locations), 4 alignment positions on the WNV capsid genome, and 5 alignment positions on the ZIKV capsid genome, depicted by this sequence alignment analysis.
Discussion
miRNAs are a class of small non-coding RNA segments, ranging up to 22 nt long, and serve as possible inhibiting regulators against viral mRNA expression during virus replication (). A perfect complementary sequence between miRNA and mRNA regions is believed to be sufficient for successful cleavage or degradation of the mRNA sequence, but imperfect alignments may block viral translation (; ). We studied the sequence homology of our miRNA sequence and the flavivirus genome sequences, and we found that after numerous sequence alignments, this study confirms the significant complementary sequence of our candidate miRNA sequence on the flavivirus capsid genome. After careful analysis, we have analyzed that hsa-miR-548d-3p showed identical alignment locations on the capsid gene of DENV 1, 3, and 4 viruses with some minor differences (Figure 11; Table 3). Figure 13 also confirms that hsa-miR-548d-3p also has identical alignment locations on the HCV virus and its genotypes. Hence, those miRNAs are used as antiviral therapeutics; these findings suggest that hsa-miR-548d-3p may be a possible candidate as a universal antiviral therapeutic agent against infections caused by the flavivirus family.
Alignment data of NCBI BLAST-n
To understand the idea of a good alignment between two sequences, it is necessary to understand the score, E-value (expected value), percentage of identity, and gaps. The bits score indicates how significant the alignment is; the higher the score on the alignment, the better. Observing the expected or, simply, the E-value indicates the significance of an alignment; the lower the E-value signifies, the better the alignment between two sequences. According to Tom Madden, if an alignment has an E-value of 0.05, then the similarities have a 5 in 100 possibility of occurring by chance (). The percentage (%) of identity signifies how perfect the alignment is between two sequences. Table 3 shows that all the alignments demonstrated are 7/7 or 8/8, which indicates two sequences are 100% similar. According to , the percentage of identity of 100% means that the nucleotides of the subject sequence are identical to the reference sequence or the query at every position of the alignment.
Using miRBase for the prediction of miRNA targets
Many bioinformatics tools were developed for biogenesis and to help biologists investigate miRNA biology. Among these tools, miRBase (https://mirbase.org/) is the most widely used software program (), which was developed in the year 2002 (). Later, the name changed from “the microRNA Registry” provided molecular researchers with stable and unique gene names for their novel miRNA discoveries and storages of miRNA sequences (). miRBase is a primary repository database for retrieving data on miRNA and has three main functions: 1) provides confidential services for independent assignment of miRNA genes, 2) sequences provide miRNA data, annotation, references, and links to other published miRNAs, and 3) provides miRNA target pipelines for the prediction of the target (; ). miRBase also allows searching published pre-miRNA and mature miRNA sequences, in addition to readily available annotation and sequence data that are available for download (). Overall, miRBase provides scientists a variety of data on miRNAs when obtaining sequences that include the accession number, symbols, description, and gene family.
Using miRBase and BLASTn
One of the important function of miRbase is to provide a microRNA target pipeline for the prediction of targets for all published animal miRNAs (; ), and all miRNA sequences from this database revealed interactions against 3′-untranslated regions, which are predicted from all available species ().
We used BLAST to perform alignments for its potential to detect similar regions within parts of long sequences. It is a fast, sensitive, and accurate tool for analyzing sequences for alignments (; ). The reason that miRBase does not run alignment analysis was BLAST was used when searching for the right miRNA sequences over 2,600 miRNAs. Additionally, BLAST was used because it is the most widely used software package in bioinformatics research due to its main function of comparing sequence(s) of interest (Stover and Cavalcanti, 2017).
Roles of hsa-miR-548d-3p in humans
hsa-miR-548d-3p is a mature miRNA that is found in primates and comprises over 69 identified miRNAs that are presented in all human chromosomes, but also as a more poorly conserved miRNA (Ramos-Sanchez et al., 2022), it demonstrates to enhance cell proliferation and inhibit apoptosis in breast cancer (Souza et al., 2016; Souza et al., 2021). Functions can include many biological processes, such as signaling pathways like MAPK, phosphatidylinositol (P13K), p53, B-cell receptor, T-cell receptor, TGF-beta, PPAR, calcium, and insulin signaling pathways, and in human tumorigenesis, such as colorectal cancer, glioma, and non-small cell lung cancer (Ramox-Sanchez, 2022; ). hsa-miR-548d-3p is proven to be involved in homeostasis of stress damage, and metabolic and survival pathways for cell proliferation (; ). In an experiment done by Rooda L. et al., their results indicated that hsa-miR-548d-3p and its family may play additional roles in humans, such as in ovarian follicle activation, development, granulosa cell differentiation, and proliferation (Rooda et al., 2021).
Bovine viral diarrhea virus as a model for flaviviruses
Prestiviruses are more closely related to HCV than the classical flaviviruses and have been used as surrogate models for HCV (Tellinghuisen et al., 2006; ) to test in vitro infectivity (). According to , as one of the most characterized members of the Flaviviridae family, BVDV serves as a good model system to study flaviviruses and has primarily been used as a surrogate model for HCV in identification and characterization of antiviral agents (). This approach leverages the similarities between BVDV and HCV to develop and test potential treatments for HCV more effectively. stated that both viruses BVDV and HCV utilize the IRES within the 5′ Untranslated Region (UTR) necessary for translation of viral polyprotein, while NS3 proteases of both viruses require NS4A as a cofactor for polyprotein processing.
Limitations
This was a pure computer-based study using bioinformatic tools to showcase possible miRNA–mRNA sequence similarities. Due to pestiviruses like the BVDV, which is used as a surrogate model for studying HCV virus, we hypothesized that if we can utilize the BVDV genome sequence as our test subject, then we could find a possible universal miRNA-based antiviral therapeutic for the family of flaviviruses. Based on our results and Table 3 and Figures 10–16, we found hsa-miR-548d-3p as a possible candidate due to its perfect match with the genome of all our viruses rather than just one. Again, this is a full bioinformatic-based analytical study, where in vivo lab equipment was not used. Hence, the results are not considered final until proved using in vivo experimentation.
Conclusion
After performing a series of sequence alignments, we predicted hsa-miR-548d-3p, a mature miRNA sequence, as a potential candidate to target flaviviruses showing perfect alignments with BVDV; HCV genotype 1a, 2, 3, 4, 5, and 6; DENV serotype 1, 2, 3, and 4; JENV; WNV; ZIKA; and TBEV. Although more detailed in vitro and in vivo studies are required to utilize hsa-miR-548d-3p as an antiviral therapeutic, this study may be considered a first step to develop a new type of miRNA treatment against a range of viruses within the Flaviviridae family. This study also recognizes that the BVDV may not be the surrogate model for only HCV virus but can also prove to be a good model system for antiviral therapeutic studies against other members of the Flaviviridae family.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.
Author contributions
HC: writing–original draft and writing–review and editing. SH: writing–original draft and writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The authors gratefully acknowledge the moral support from School of STEM, Diné College, and financial support from MS Biology Grant # NSF-TCUP-ICE-TI-2202023 (PI: Donald K. Robinson).
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.
Publisher’s note
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.
References
1
AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool. J. Mol. Biol.215 (3), 403–410. 10.1006/jmbi.1990.9999
2
ArdekaniA. M.NaeiniM. M. (2010). The role of MicroRNAs in human diseases. Avicenna J. Med. Biotechnol.2 (4), 161–179.
3
BabikerA.JeudyJ.KligermanS.KhambatyM.ShahA.BagchiS. (2017). Risk of cardiovascular disease due to chronic hepatitis C infection: a review. J. Clin. Transl. hepatology5 (4), 1–20. 10.14218/jcth.2017.00021
4
BaekD.VillénJ.ShinC.CamargoF. D.GygiS. P.BartelD. P. (2008). The impact of microRNAs on protein output. Nature455 (7209), 64–71. 10.1038/nature07242
5
BarrowsN. J.CamposR. K.LiaoK. C.PrasanthK. R.Soto-AcostaR.YehS. C.et al (2018). Biochemistry and molecular biology of flaviviruses. Chem. Rev.118 (8), 4448–4482. 10.1021/acs.chemrev.7b00719
6
BasarabM.BowmanC.AaronsE. J.CropleyI. (2016). Zika virus. Bmj352, i1049. 10.1136/bmj.i1049
7
BecherP.OrlichM.ThielH. J. (1998). Complete genomic sequence of border disease virus, a pestivirus from sheep. J. virology72 (6), 5165–5173. 10.1128/jvi.72.6.5165-5173.1998
8
BenzartiE.LindenA.DesmechtD.GariglianyM. (2019). Mosquito-borne epornitic flaviviruses: an update and review. J. General Virology100 (2), 119–132. 10.1099/jgv.0.001203
9
CannataroR.CioneE. (2019). Antioxidant and microRNAs: an applied overview. J. Microbiol. Biotechnol. Rep.
10
CasalM. L.DambachD. M.MeisterT.JezykP. F.PattersonD. F.HenthornP. S. (2004). Familial glomerulonephropathy in the bullmastiff. Veterinary pathol.41 (4), 319–325. 10.1354/vp.41-4-319
11
ChenL.HeikkinenL.WangC.YangY.SunH.WongG. (2019). Trends in the development of miRNA bioinformatics tools. Briefings Bioinforma.20 (5), 1836–1852. 10.1093/bib/bby054
12
ChenN.LiuY.BaiT.ChenJ.ZhaoZ.LiJ.et al (2022). Quercetin inhibits Hsp70 blocking of bovine viral diarrhea virus infection and replication in the early stage of virus infection. Viruses14 (11), 2365. 10.3390/v14112365
13
ChiS.ChenS.JiaW.HeY.RenL.WangX. (2022). Non-structural proteins of bovine viral diarrhea virus. Virus genes58 (6), 491–500. 10.1007/s11262-022-01914-8
14
CondeJ. N.SilvaE. M.BarbosaA. S.Mohana-BorgesR. (2017). The complement system in flavivirus infections. Front. Microbiol.8, 213. 10.3389/fmicb.2017.00213
15
de Oliveira FigueiredoP.Stoffella-DutraA. G.Barbosa CostaG.Silva de OliveiraJ.Dourado AmaralC.Duarte SantosJ.et al (2020). Re-emergence of yellow fever in Brazil during 2016-2019: challenges, lessons learned, and perspectives. Viruses12 (11), 1233. 10.3390/v12111233
16
DuncanJ. D.UrbanowiczR. A.TarrA. W.BallJ. K. (2020). Hepatitis C virus vaccine: challenges and prospects. Vaccines8 (1), 90. 10.3390/vaccines8010090
17
DurantelD.Carrouée-DurantelS.Branza-NichitaN.DwekR. A.ZitzmannN. (2004). Effects of interferon, ribavirin, and iminosugar derivatives on cells persistently infected with noncytopathic bovine viral diarrhea virus. Antimicrob. agents Chemother.48 (2), 497–504. 10.1128/aac.48.2.497-504.2004
18
FasslerJ.CooperP. (2008). “BLAST glossary. 2011 jul 14,” in BLAST® help (Bethesda (MD): National Center for Biotechnology Information US). Available at: https://www.ncbi.nlm.nih.gov/books/NBK62051/.
19
FinkielszteinL. M.MoltrasioG.CaputtoM.CastroE.CavallaroL.MoglioniA. (2010). What is known about the antiviral agents active against bovine viral diarrhea virus (BVDV)?Curr. Med. Chem.17 (26), 2933–2955. 10.2174/092986710792065036
20
FinneganE. F.PasquinelliA. E. (2013). MicroRNA biogenesis: regulating the regulators. Crit. Rev. Biochem. Mol. Biol.48 (1), 51–68. 10.3109/10409238.2012.738643
21
FuG.BrkićJ.HayderH.PengC. (2013). MicroRNAs in human placental development and pregnancy complications. Int. J. Mol. Sci.14 (3), 5519–5544. 10.3390/ijms14035519
22
FultonR. W.PurdyC. W.ConferA. W.SalikiJ. T.LoanR. W.BriggsR. E.et al (2000). Bovine viral diarrhea viral infections in feeder calves with respiratory disease: interactions with Pasteurella spp., parainfluenza-3 virus, and bovine respiratory syncytial virus. Can. J. veterinary Res. = Revue Can. de recherche veterinaire64 (3), 151–159.
23
Giménez-RicharteÁ.Ortiz de SalazarM. I.Giménez-RicharteM. P.ColladoM.FernándezP. L.ClavijoC.et al (2022). Transfusion-transmitted arboviruses: update and systematic review. PLOS Neglected Trop. Dis.16 (10), e0010843. 10.1371/journal.pntd.0010843
24
Griffiths-JonesS.et al (2006). miRBase: microRNA sequences, targets and gene nomenclature. Nucleic acids Res.34 (Database issue), D140–D144. 10.1093/nar/gkj112
25
Griffiths-JonesS.SainiH. K.van DongenS.EnrightA. J. (2007). miRBase: tools for microRNA genomics. Nucleic acids Res.36 (Suppl. l_1), D154–D158. 10.1093/nar/gkm952
26
HaM.KimV. N. (2014). Regulation of microRNA biogenesis. Nat. Rev. Mol. cell Biol.15 (8), 509–524. 10.1038/nrm3838
27
HamelR.DejarnacO.WichitS.EkchariyawatP.NeyretA.LuplertlopN.et al (2015). Biology of Zika virus infection in human skin cells. J. virology89 (17), 8880–8896. 10.1128/jvi.00354-15
28
HaradaT.TautzN.ThielH. J. (2000). E2-p7 region of the bovine viral diarrhea virus polyprotein: processing and functional studies. J. virology74 (20), 9498–9506. 10.1128/jvi.74.20.9498-9506.2000
29
HasanM. M.AkterR.UllahM. S.AbedinM. J.UllahG. M. A.HossainM. Z. (2014). A computational approach for predicting role of human microRNAs in MERS‐CoV genome. Adv. Bioinforma.2014 (1), 1–8. 10.1155/2014/967946
30
HillsS. L.FischerM.PetersenL. R. (2017). Epidemiology of Zika virus infection. J. Infect. Dis.216 (Suppl. l_10), S868–S874. 10.1093/infdis/jix434
31
HuT.WuZ.WuS.ChenS.ChengA. (2021). The key amino acids of E protein involved in early flavivirus infection: viral entry. Virology J.18 (1), 136. 10.1186/s12985-021-01611-2
32
HauseB. MPillatzkiA.ClementT.BraggT.RidpathJ.ChaseC. C. L. (2021). Persistent infection of American bison (Bison bison) with bovine viral diarrhea virus and bosavirus. Vet. Microbiol.252, 108949. 10.1016/j.vetmic.2020.108949
33
IskenO.BarothM.GrassmannC. W.WeinlichS.OstareckD. H.Ostareck-LedererA.et al (2007). Nuclear factors are involved in hepatitis C virus RNA replication. Rna13 (10), 1675–1692. 10.1261/rna.594207
34
JackovaA.NovackovaM.PelletierC.AudevalC.GueneauE.HaffarA.et al (2008). The extended genetic diversity of BVDV-1: typing of BVDV isolates from France. Veterinary Res. Commun.32 (1), 7–11. 10.1007/s11259-007-9012-z
35
JulanderJ. G.ShaferK.SmeeD. F.MorreyJ. D.FurutaY. (2009). Activity of T-705 in a hamster model of yellow fever virus infection in comparison with that of a chemically related compound, T-1106. Antimicrob. agents Chemother.53 (1), 202–209. 10.1128/aac.01074-08
36
JunglenS.KoppA.KurthA.PauliG.EllerbrokH.LeendertzF. H. (2009). A new flavivirus and a new vector: characterization of a novel flavivirus isolated from uranotaenia mosquitoes from a tropical rain forest. J. virology83 (9), 4462–4468. 10.1128/jvi.00014-09
37
Khodakaram-TaftiA.FarjanikishG. H. (2017). Persistent bovine viral diarrhea virus (BVDV) infection in cattle herds. Iran. J. veterinary Res.18 (3), 154–163.
38
KosinovaE.PsikalI.RobesovaB.KovarcikK. (2007). Real-time PCR for quantitation of bovine viral diarrhea virus RNA using SYBR Green I fluorimetry. VETERINARNI MEDICINA-PRAHA-52 (6), 253–261. 10.17221/1882-vetmed
39
KozomaraA.Griffiths-JonesS. (2010). miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic acids Res.39 (Suppl. l_1), D152–D157. 10.1093/nar/gkq1027
40
LacknerT.MüllerA.PankrazA.BecherP.ThielH. J.GorbalenyaA. E.et al (2004). Temporal modulation of an auto protease is crucial for replication and pathogenicity of an RNA virus. J. virology78 (19), 10765–10775. 10.1128/jvi.78.19.10765-10775.2004
41
LaiV. C.ZhongW.SkeltonA.IngravalloP.VassilevV.DonisR. O.et al (2000). Generation and characterization of a hepatitis C virus NS3 protease-dependent bovine viral diarrhea virus. J. virology74 (14), 6339–6347. 10.1128/jvi.74.14.6339-6347.2000
42
LauretiM.NarayananD.Rodriguez-AndresJ.FazakerleyJ. K.KedzierskiL. (2018). Flavivirus receptors: diversity, identity, and cell entry. Front. Immunol.9, 2180. 10.3389/fimmu.2018.02180
43
LiF.WangY.YuL.CaoS.WangK.YuanJ.et al (2015). Viral infection of the central nervous system and neuroinflammation precede blood-brain barrier disruption during Japanese encephalitis virus infection. J. virology89 (10), 5602–5614. 10.1128/jvi.00143-15
44
LiY.WangJ.KanaiR.ModisY. (2013). Crystal structure of glycoprotein E2 from bovine viral diarrhea virus. Proc. Natl. Acad. Sci.110 (17), 6805–6810. 10.1073/pnas.1300524110
45
LiangT.GuoL.LiuC. (2012). Genome‐wide analysis of mir‐548 gene family reveals evolutionary and functional implications. BioMed Res. Int.2012 (1), 1–8. 10.1155/2012/679563
46
LoboI. (2008). Basic local alignment search tool (BLAST). Nat. Educ.1 (1).
47
Luna BuitragoD.LoveringR. C.CaporaliA. (2023). Insights into online microRNA bioinformatics tools. Non-coding RNA9 (2), 18. 10.3390/ncrna9020018
48
MaddenT. (2013). The BLAST sequence analysis tool. NCBI Handb.2 (5), 425–436.
49
MaiorinoM.Bosello-TravainV.CozzaG.MiottoG.RoveriA.ToppoS.et al (2015). Understanding mammalian glutathione peroxidase 7 in the light of its homologs. Free Radic. Biol. Med.83, 352–360. 10.1016/j.freeradbiomed.2015.02.017
50
MariV. L.LosurdoM.LucenteM. S.LorussoE.EliaG.MartellaV.et al (2016). Multiplex real-time RT-PCR assay for bovine viral diarrhea virus type 1, type 2 and HoBi-like pestivirus. J. virological methods229, 1–7. 10.1016/j.jviromet.2015.12.003
51
MaurerK.KreyT.MoennigV.ThielH. J.RümenapfT. (2004). CD46 is a cellular receptor for bovine viral diarrhea virus. J. virology78 (4), 1792–1799. 10.1128/jvi.78.4.1792-1799.2004
52
MerwaissF.CzibenerC.AlvarezD. E. (2019). Cell-to-cell transmission is the main mechanism supporting bovine viral diarrhea virus spread in cell culture. J. Virology93 (3), 017766-18–e11128. 10.1128/jvi.01776-18
53
MurrayC. L.MarcotrigianoJ.RiceC. M. (2008). Bovine viral diarrhea virus core is an intrinsically disordered protein that binds RNA. J. virology82 (3), 1294–1304. 10.1128/jvi.01815-07
54
MussoD.GublerD. J. (2016). Zika virus. Clin. Microbiol. Rev.29 (3), 487–524. 10.1128/cmr.00072-15
55
MussoD.NhanT. X. (2015). Emergence of Zika virus. Clin. Microbiol.4, 222. 10.4172/2327-5073.1000222
56
MutebiJ. P.RijnbrandR. C. A.WangH.RymanK. D.WangE.FulopL. D.et al (2004). Genetic relationships and evolution of genotypes of yellow fever virus and other members of the yellow fever virus group within the Flavivirus genus based on the 3′ noncoding region. J. virology78 (18), 9652–9665. 10.1128/jvi.78.18.9652-9665.2004
57
NeillJ. D. (2013). Molecular biology of bovine viral diarrhea virus. Biologicals41 (1), 2–7. 10.1016/j.biologicals.2012.07.002
58
NijmeijerB. M.KoopsenJ.SchinkelJ.PrinsM.GeijtenbeekT. B. (2019). Sexually transmitted hepatitis C virus infections: current trends, and recent advances in understanding the spread in men who have sex with men. J. Int. AIDS Soc.22 (l), e25348. 10.1002/jia2.25348
59
NiskanenR.LindbergA.LarssonB.AleniusS. (2000). Lack of virus transmission from bovine viral diarrhoea virus infected calves to susceptible peers. Acta Veterinaria Scand.41 (1), 93–99. 10.1186/bf03549659
60
O'BrienJ.HayderH.ZayedY.PengC. (2018). Overview of microRNA biogenesis, mechanisms of actions, and circulation. Front. Endocrinol.9, 402. 10.3389/fendo.2018.00402
61
Perera-LecoinM.MeertensL.CarnecX.AmaraA. (2013). Flavivirus entry receptors: an update. Viruses6 (1), 69–88. 10.3390/v6010069
62
PeterhansE.JungiT. W.SchweizerM. (2003). BVDV and innate immunity. Biologicals31 (2), 107–112. 10.1016/s1045-1056(03)00024-1
63
PiersonT. C.DiamondM. S. (2020). The continued threat of emerging flaviviruses. Nat. Microbiol.5 (6), 796–812. 10.1038/s41564-020-0714-0
64
PiersonT. C.KielianM. (2013). Flaviviruses: braking the entering. Curr. Opin. virology3 (1), 3–12. 10.1016/j.coviro.2012.12.001
65
QianQ.XuR.WangY., (2022). “The NS4A protein of classical swine fever virus suppresses RNA silencing in mammalian cells.” J. Virology, 96(15), 018744-e1921. 10.1128/jvi.01874-21
66
Ramos-SanchezE. M.ReisL. C.SouzaM. A.MuxelS. M.SantosK. R.LagosD.et al (2022). miR-548d-3p is up-regulated in human visceral leishmaniasis and suppresses parasite growth in macrophages. Front. Cell. Infect. Microbiol.12, 826039. 10.3389/fcimb.2022.826039
67
RanganathanK.SivasankarV. (2014). MicroRNAs-Biology and clinical applications. J. Oral Maxillofac. Pathology18 (2), 229–234. 10.4103/0973-029x.140762
68
ReedK. E.GorbalenyaA. E.RiceC. M. (1998). The NS5A/NS5 proteins of viruses from three genera of the family flaviviridae are phosphorylated by associated serine/threonine kinases. J. virology72 (7), 6199–6206. 10.1128/jvi.72.7.6199-6206.1998
69
RoodaI.KaseltB.LiivrandM.SmolanderO. P.SalumetsA.Velthut-MeikasA. (2021). Hsa-mir-548 family expression in human reproductive tissues. BMC Genomic Data22, 40–13. 10.1186/s12863-021-00997-w
70
SchweizerM.PeterhansE. (2001). Noncytopathic bovine viral diarrhea virus inhibits double-stranded RNA-induced apoptosis and interferon synthesis. J. Virology75 (10), 4692–4698. 10.1128/jvi.75.10.4692-4698.2001
71
SelbachM.SchwanhäusserB.ThierfelderN.FangZ.KhaninR.RajewskyN. (2008). Widespread changes in protein synthesis induced by microRNAs. nature455 (7209), 58–63. 10.1038/nature07228
72
SimmonsC. P.FarrarJ. J.van Vinh ChauN.WillsB. (2012). Dengue. N. Engl. J. Med.366 (15), 1423–1432. 10.1056/nejmra1110265
73
SkalskyR. L.CullenB. R. (2010). Viruses, microRNAs, and host interactions. Annu. Rev. Microbiol.64, 123–141. 10.1146/annurev.micro.112408.134243
74
SongZ. Q.HaoF.MinF.MaQ. Y.LiuG. D. (2001). Hepatitis C virus infection of human hepatoma cell line 7721 in vitro. World J. Gastroenterology7 (5), 685. 10.3748/wjg.7.685
75
SouzaM. A.AlmeidaT. M.CastroM. C. A.Oliveira-MendesA. P.AlmeidaA. F.OliveiraB. C.et al (2016). American tegumentary leishmaniasis: mRNA expression for Th1 and Treg mediators are predominant in patients with recent active disease. Immunobiology221 (2), 253–259. 10.1016/j.imbio.2015.08.009
76
SouzaM. D. A.Ramos-SanchezE. M.MuxelS. M.LagosD.ReisL. C.PereiraV. R. A.et al (2021). miR-548d-3p alters parasite growth and inflammation in leishmania (Viannia) braziliensis infection. Front. Cell. Infect. Microbiol.11, 687647. 10.3389/fcimb.2021.687647
77
StoverN. A.CavalcantiA. R. (2017). Using NCBI BLAST. Curr. Protoc. Essent. Lab. Tech.14 (1), 11–1. 10.1002/cpet.8
78
SuzichJ. A.TamuraJ. K.Palmer-HillF.WarrenerP.GrakouiA.RiceC. M.et al (1993). Hepatitis C virus NS3 protein polynucleotide-stimulated nucleoside triphosphatase and comparison with the related pestivirus and flavivirus enzymes. J. virology67 (10), 6152–6158. 10.1128/jvi.67.10.6152-6158.1993
79
TellinghuisenT. L.PaulsonM. S.RiceC. M. (2006). The NS5A protein of bovine viral diarrhea virus contains an essential zinc-binding site similar to that of the hepatitis C virus NS5A protein. J. virology80 (15), 7450–7458. 10.1128/jvi.00358-06
80
TurtleL.GriffithsM. J.SolomonT. (2012). Encephalitis caused by flaviviruses. QJM Mon. J. Assoc. Physicians105 (3), 219–223. 10.1093/qjmed/hcs013
81
van den ElsenK.QuekJ. P.LuoD. (2021). Molecular insights into the flavivirus replication complex. Viruses13 (6), 956. 10.3390/v13060956
82
van LeurS. W.HeunisT.MunnurD.SanyalS. (2021). Pathogenesis and virulence of flavivirus infections. Virulence12 (1), 2814–2838. 10.1080/21505594.2021.1996059
83
WaggonerJ. J.RojasA.PinskyB. A. (2018). Yellow fever virus: diagnostics for a persistent arboviral threat. J. Clin. Microbiol.56 (10), e00827-18–e01128. 10.1128/jcm.00827-18
84
WarrenerP.CollettM. S. (1995). Pestivirus NS3 (p80) protein possesses RNA helicase activity. J. virology69 (3), 1720–1726. 10.1128/jvi.69.3.1720-1726.1995
85
World Health Organization (2019). Japanese encephalitis. World Health Organization. Available at: https://www.who.int/news-room/fact-sheets/detail/japanese-encephalitis.
86
XueW.MattickD.SmithL.MaxwellJ. (2009). Fetal protection against bovine viral diarrhea virus types 1 and 2 after the use of a modified-live virus vaccine. Can. J. Veterinary Res.73 (4), 292–297.
Summary
Keywords
Flavivirus, miRNA, BLAST, NCBI, alignments, antiviral
Citation
Cayatineto HW and Hakim ST (2025) hsa-miR-548d-3p: a potential microRNA to target nucleocapsid and/or capsid genes in multiple members of the Flaviviridae family . Front. Bioinform. 4:1487292. doi: 10.3389/fbinf.2024.1487292
Received
27 August 2024
Accepted
11 December 2024
Published
14 January 2025
Volume
4 - 2024
Edited by
Deepak Agarwal, Tamil Nadu Fisheries University, India
Reviewed by
Nan Zheng, Nanjing University, China
Muhammad Aleem Ashraf, Emerson University Multan, Pakistan
Carolina Polonio, Brigham and Women’s Hospital and Harvard Medical School, United States
Updates
Copyright
© 2025 Cayatineto and Hakim.
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: S. T. Hakim, stabassum@dinecollege.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.