Evolution of coronavirus frameshifting elements: Competing stem networks explain conservation and variability.

The frameshifting RNA element (FSE) in coronaviruses (CoVs) regulates the programmed -1 ribosomal frameshift (-1 PRF) mechanism common to many viruses. The FSE is of particular interest as a promising drug candidate. Its associated pseudoknot or stem loop structure is thought to play a large role in frameshifting and thus viral protein production. To investigate the FSE structural evolution, we use our graph theory-based methods for representing RNA secondary structures in the RNA-As-Graphs (RAG) framework to calculate conformational landscapes of viral FSEs with increasing sequence lengths for representative 10 Alpha and 13 Beta-CoVs. By following length-dependent conformational changes, we show that FSE sequences encode many possible competing stems which in turn favor certain FSE topologies, including a variety of pseudoknots, stem loops, and junctions. We explain alternative competing stems and topological FSE changes by recurring patterns of mutations. At the same time, FSE topology robustness can be understood by shifted stems within different sequence contexts and base pair coevolution. We further propose that the topology changes reflected by length-dependent conformations contribute to tuning the frameshifting efficiency. Our work provides tools to analyze virus sequence/structure correlations, explains how sequence and FSE structure have evolved for CoVs, and provides insights into potential mutations for therapeutic applications against a broad spectrum of CoV FSEs by targeting key sequence/structural transitions.

TUT4/7-mediated uridylation of a coronavirus subgenomic RNAs delays viral replication.

Coronaviruses are positive-strand RNA viruses with 3' polyadenylated genomes and subgenomic transcripts. The lengths of the viral poly(A) tails change during infection by mechanisms that remain poorly understood. Here, we use a splint-ligation method to measure the poly(A) tail length and poly(A) terminal uridylation and guanylation of the mouse hepatitis virus (MHV) RNAs. Upon infection of 17-CL1 cells with MHV, a member of the Betacoronavirus genus, we observe two populations of terminally uridylated viral transcripts, one with poly(A) tails ~44 nucleotides long and the other with poly(A) tails shorter than ~22 nucleotides. The mammalian terminal uridylyl-transferase 4 (TUT4) and terminal uridylyl-transferase 7 (TUT7), referred to as TUT4/7, add non-templated uracils to the 3'-end of endogenous transcripts with poly(A) tails shorter than ~30 nucleotides to trigger transcript decay. Here we find that depletion of the host TUT4/7 results in an increased replication capacity of the MHV virus. At late stages of infection, the population of uridylated subgenomic RNAs with tails shorter than ~22 nucleotides is reduced in the absence of TUT4/7 while the viral RNA load increases. Our findings indicate that TUT4/7 uridylation marks the MHV subgenomic RNAs for decay and delays viral replication.

Global Dynamics of Porcine Enteric Coronavirus PEDV Epidemiology, Evolution, and Transmission.

With a possible origin from bats, the alphacoronavirus Porcine epidemic diarrhea virus (PEDV) causes significant hazards and widespread epidemics in the swine population. However, the ecology, evolution, and spread of PEDV are still unclear. Here, from 149,869 fecal and intestinal tissue samples of pigs collected in an 11-year survey, we identified PEDV as the most dominant virus in diarrheal animals. Global whole genomic and evolutionary analyses of 672 PEDV strains revealed the fast-evolving PEDV genotype 2 (G2) strains as the main epidemic viruses worldwide, which seems to correlate with the use of G2-targeting vaccines. The evolving pattern of the G2 viruses presents geographic bias as they evolve tachytely in South Korea but undergo the highest recombination in China. Therefore, we clustered six PEDV haplotypes in China, whereas South Korea held five haplotypes, including a unique haplotype G. In addition, an assessment of the spatiotemporal spread route of PEDV indicates Germany and Japan as the primary hubs for PEDV dissemination in Europe and Asia, respectively. Overall, our findings provide novel insights into the epidemiology, evolution, and transmission of PEDV, and thus may lay a foundation for the prevention and control of PEDV and other coronaviruses.

Prediction and Verification of Curcumin as a Potential Drug for Inhibition of PDCoV Replication in LLC-PK1 Cells.

Porcine deltacoronavirus (PDCoV) is an emerging swine enteropathogenic coronavirus (CoV) that causes lethal watery diarrhea in neonatal pigs and poses economic and public health burdens. Currently, there are no effective antiviral agents against PDCoV. Curcumin is the active ingredient extracted from the rhizome of turmeric, which has a potential pharmacological value because it exhibits antiviral properties against several viruses. Here, we described the antiviral effect of curcumin against PDCoV. At first, the potential relationships between the active ingredients and the diarrhea-related targets were predicted through a network pharmacology analysis. Twenty-three nodes and 38 edges were obtained using a PPI analysis of eight compound-targets. The action target genes were closely related to the inflammatory and immune related signaling pathways, such as the TNF signaling pathway, Jak-STAT signaling pathway, and so on. Moreover, IL-6, NR3C2, BCHE and PTGS2 were identified as the most likely targets of curcumin by binding energy and 3D protein-ligand complex analysis. Furthermore, curcumin inhibited PDCoV replication in LLC-PK1 cells at the time of infection in a dose-dependent way. In poly (I:C) pretreated LLC-PK1 cells, PDCoV reduced IFN-ß production via the RIG-I pathway to evade the host's antiviral innate immune response. Meanwhile, curcumin inhibited PDCoV-induced IFN-ß secretion by inhibiting the RIG-I pathway and reduced inflammation by inhibiting IRF3 or NF-κB protein expression. Our study provides a potential strategy for the use of curcumin in preventing diarrhea caused by PDCoV in piglets.

Three Amino Acid Substitutions in the Spike Protein Enable the Coronavirus Porcine Epidemic Diarrhea Virus To Infect Vero Cells.

Porcine epidemic diarrhea virus (PEDV), a continuously evolving pathogen, causes severe diarrhea in piglets, with high mortality rates. To prevent or mitigate the disease, it is common practice to develop live or inactivated PEDV vaccines based on cell-adapted viral variants. Propagating wild-type PEDV in cultured cells is, however, often challenging due to the lack of knowledge about the requirements for the cell adaptation of PEDV. In the present study, by using the RNA-targeted reverse genetic system for PEDV to apply S protein swapping followed by the rescue of the recombinant viruses, three key amino acid mutations in the S protein, A605E, E633Q, and R891G, were identified, which enable attenuated PEDV strain DR13 (DR13att) to efficiently and productively infect Vero cells, in contrast to the parental DR13 strain (DR13par). The former two key mutations reside inside and in the vicinity of the receptor binding domain (RBD), respectively, while the latter occurs at the N-terminal end of the fusion peptide (FP). Besides the three key mutations, other mutations in the S protein further enhanced the infection efficiency of the recombinant viruses. We hypothesize that the three mutations changed PEDV tropism by altering the S2' cleavage site and the RBD structure. This study provides basic molecular insight into cell adaptation by PEDV, which is also relevant for vaccine design. IMPORTANCE Porcine epidemic diarrhea virus (PEDV) is a lethal pathogen for newborn piglets, and an efficient vaccine is needed urgently. However, propagating wild-type PEDV in cultured cells for vaccine development is still challenging due to the lack of knowledge about the mechanism of the cell adaptation of PEDV. In this study, we found that three amino acid mutations, A605E, E633Q, and R891G, in the spike protein of the Vero cell-adapted PEDV strain DR13att were critical for its cell adaptation. After analyzing the mutation sites in the spike protein, we hypothesize that the cell adaptation of DR13att was achieved by altering the S2' cleavage site and the RBD structure. This study provides new molecular insight into the mechanism of PEDV culture adaptation and new strategies for PEDV vaccine design.

Characterization of RNA G-quadruplexes in porcine epidemic diarrhea virus genome and the antiviral activity of G-quadruplex ligands.

Porcine epidemic diarrhea virus (PEDV), an enteropathogenic coronavirus, has catastrophic impacts on the global pig industry. However, there are still no anti-PEDV drugs with accurate targets. G-quadruplexes (G4s) are non-canonical secondary structures formed within guanine-rich regions of DNA or RNA, and have attracted great attention as potential targets for antiviral strategy. In this study, we reported two putative G4-forming sequences (PQS) in S and Nsp5 genes of PEDV genome based on bioinformatic analysis, and identified that S-PQS and Nsp5-PQS were enabled to fold into G4 structure by using circular dichroism spectroscopy and fluorescence turn-on assay. Furthermore, we verified that both S-PQS and Nsp5-PQS PQS could form G4 structure in live cells by immunofluorescence microscopy. In addition, G4-specific compounds, such as TMPyP4 and PDS, could significantly inhibit transcription, translation and proliferation of PEDV in vitro. Importantly, these compounds exert antiviral activity at the post-entry step of PEDV infection cycle, by inhibiting viral genome replication and protein expression. Lastly, we demonstrated that TMPyP4 can inhibit reporter gene expression by targeting G4 structure in Nsp5. Taken together, these findings not only reinforce the presence of viral G-quadruplex sequences in PEDV genome but also provide new insights into developing novel antiviral drugs targeting PEDV RNA G-quadruplexes.

Seroprevalence of four endemic human coronaviruses and, reactivity and neutralization capability against SARS-CoV-2 among children in the Philippines.

Four endemic human coronaviruses (HCoV), HCoV-229E, HCoV-NL63, HCoV-HKU1, and HCoV-OC43, are closely related to SARS-CoV-2. These coronaviruses are known to infect humans living in temperate areas, including children under 5 years old; however, the seroprevalence of four HCoVs among children in tropical areas, including the Philippines, remains unclear. This study aimed to assess the prevalence of antibodies against four HCoVs and to determine the reactivity and neutralization of these antibodies against SARS-CoV-2 among children in the Philippines. A total of 315 serum samples collected from 2015 to 2018, before the emergence of SARS-CoV-2, in Biliran island, Philippines, were tested for the presence of antibodies against four HCoVs and SARS-CoV-2 using recombinant spike ectodomain proteins by IgG-enzyme-linked immunosorbent assay (ELISA). Reactivity to and neutralization of SARS-CoV-2 were also investigated. The seroprevalence of the four HCoVs was 63.8% for HCoV-229E, 71.4% for HCoV-NL63, 76.5% for HCoV-HKU1, and 83.5% for HCoV-OC43 by ELISA. Age group analysis indicated that seropositivity to all HCoVs reached 80% by 2-3 years of age. While 69/315 (21.9%) of the samples showed reactive to SARS-CoV-2, almost no neutralization against SARS-CoV-2 was detected using neutralization assay. Reactivity of antibodies against SARS-CoV-2 spike protein obtained by ELISA may not correlate with neutralization capability.

Detection of coronaviruses in insectivorous bats of Fore-Caucasus, 2021.

Coronaviruses (CoVs) pose a huge threat to public health as emerging viruses. Bat-borne CoVs are especially unpredictable in their evolution due to some unique features of bat physiology boosting the rate of mutations in CoVs, which is already high by itself compared to other viruses. Among bats, a meta-analysis of overall CoVs epizootiology identified a nucleic acid observed prevalence of 9.8% (95% CI 8.7-10.9%). The main objectives of our study were to conduct a qPCR screening of CoVs' prevalence in the insectivorous bat population of Fore-Caucasus and perform their characterization based on the metagenomic NGS of samples with detected CoV RNA. According to the qPCR screening, CoV RNA was detected in 5 samples, resulting in a 3.33% (95% CI 1.1-7.6%) prevalence of CoVs in bats from these studied locations. BetaCoVs reads were identified in raw metagenomic NGS data, however, detailed characterization was not possible due to relatively low RNA concentration in samples. Our results correspond to other studies, although a lower prevalence in qPCR studies was observed compared to other regions and countries. Further studies should require deeper metagenomic NGS investigation, as a supplementary method, which will allow detailed CoV characterization.

Identification of coronaviruses in bats and rodents in northern and central Argentina.

Due to the present pandemic situation and the many animal species that are epidemiologically involved, there has been a surge of renewed interest in investigating the coronavirus (CoV) population circulating in wildlife, especially bats and rodents, which are potential reservoirs of new human pathogens. In Argentina, information about the viruses present in these mammals is very limited. To investigate the presence of coronaviruses in this country, we obtained 457 samples from hematophagous, insectivorous, and frugivorous bats and rodents from two regions of Argentina. We report here the detection of alphacoronavirus sequences in three groups of bats as well as in rodents. Phylogenetic analysis showed the closest relationships to alphacoronaviruses from Brazil.

Identification of a novel linear B-cell epitope in porcine deltacoronavirus nucleocapsid protein.

Porcine deltacoronavirus (PDCoV) is an emerging swine enteropathogenic coronavirus that caused diarrhea and/or vomiting in neonatal piglets worldwide. Coronaviruses nucleocapsid (N) protein is the most conserved structural protein for viral replication and possesses good antigenicity. In this study, three monoclonal antibodies (mAbs), 3B4, 4D3, and 4E3 identified as subclass IgG2aκ were prepared using the lymphocytic hybridoma technology against PDCoV N protein. Furthermore, the B-cell epitope recognized by mAb 4D3 was mapped by dozens of overlapping truncated recombinant proteins based on the western blotting. The polypeptide 28QFRGNGVPLNSAIKPVE44 (EP-4D3) in the N-terminal of PDCoV N protein was identified as the minimal linear epitope for binding mAb 4D3. And the EP-4D3 epitope's amino acid sequence homology study revealed that PDCoV strains are substantially conserved, with the exception of the Alanine43 substitution Valine43 in the China lineage, the Early China lineage, and the Thailand, Vietnam, and Laos lineage. The epitope sequences shared high similarity (94.1%) with porcine coronavirus HKU15-155 (PorCoV HKU15), Asian leopard cats coronavirus (ALCCoV), sparrow coronavirus HKU17 (SpCoV HKU17), and sparrow deltacoronavirus. In contrast, the epitope sequences shared a very low homology (11.8 to 29.4%) with other porcine CoVs (PEDV, TGEV, PRCV, SADS-CoV, PHEV). Overall, the study will enrich the biological function of PDCoV N protein and provide foundational data for further development of diagnostic applications. KEY POINTS: • Three monoclonal antibodies against PDCoV N protein were prepared. • Discovery of a novel B-cell liner epitope (28QFRGNGVPLNSAIKPVE44) of PDCoV N protein. • The epitope EP-4D3 was conserved among PDCoV strains.

Proteomic screening identifies RPLp2 as a specific regulator for the translation of coronavirus.

Viral mRNA of coronavirus translates in an eIF4E-dependent manner, and the phosphorylation of eIF4E can modulate this process, but the role of p-eIF4E in coronavirus infection is not yet entirely evident. p-eIF4E favors the translation of selected mRNAs, specifically the mRNAs that encode proteins associated with cell proliferation, inflammation, the extracellular matrix, and tumor formation and metastasis. In the present work, two rounds of TMT relative quantitative proteomics were used to screen 77 cellular factors that are upregulated upon infection by coronavirus PEDV and are potentially susceptible to a high level of p-eIF4E. PEDV infection increased the translation level of ribosomal protein lateral stalk subunit RPLp2 (but not subunit RPLp0/1) in a p-eIF4E-dependent manner. The bicistronic dual-reporter assay and polysome profile showed that RPLp2 is essential for translating the viral mRNA of PEDV. RNA binding protein and immunoprecipitation assay showed that RPLp2 interacted with PEDV 5'UTR via association with eIF4E. Moreover, the cap pull-down assay showed that the viral nucleocapsid protein is recruited in m7GTP-precipitated complexes with the assistance of RPLp2. The heterogeneous ribosomes, which are different in composition, regulate the selective translation of specific mRNAs. Our study proves that viral mRNA and protein utilize translation factors and heterogeneous ribosomes for preferential translation initiation. This previously uncharacterized process may be involved in the selective translation of coronavirus.

Effects of immunophilin inhibitors and non-immunosuppressive analogs on coronavirus replication in human infection models.

Rationale: Human coronaviruses (HCoVs) seriously affect human health by causing respiratory diseases ranging from common colds to severe acute respiratory diseases. Immunophilins, including peptidyl-prolyl isomerases of the FK506-binding protein (FKBP) and the cyclophilin family, are promising targets for pharmaceutical inhibition of coronavirus replication, but cell-type specific effects have not been elucidated. FKBPs and cyclophilins bind the immunosuppressive drugs FK506 and cyclosporine A (CsA), respectively. Methods: Primary human bronchial epithelial cells (phBECs) were treated with CsA, Alisporivir (ALV), FK506, and FK506-derived non-immunosuppressive analogs and infected with HCoV-229E. RNA and protein were assessed by RT-qPCR and immunoblot analysis. Treatment with the same compounds was performed in hepatoma cells (Huh-7.5) infected with HCoV-229E expressing Renilla luciferase (HCoV-229E-RLuc) and the kidney cell line HEK293 transfected with a SARS-CoV-1 replicon expressing Renilla luciferase (SARS-CoV-1-RLuc), followed by quantification of luminescence as a measure of viral replication. Results: Both CsA and ALV robustly inhibited viral replication in all models; both compounds decreased HCoV-229E RNA in phBECs and reduced luminescence in HCoV-229E-RLuc-infected Huh7.5 and SARS-CoV-1-RLuc replicon-transfected HEK293. In contrast, FK506 showed inconsistent and less pronounced effects in phBECs while strongly affecting coronavirus replication in Huh-7.5 and HEK293. Two non-immunosuppressive FK506 analogs had no antiviral effect in any infection model. Conclusion: The immunophilin inhibitors CsA and ALV display robust anti-coronaviral properties in multiple infection models, including phBECs, reflecting a primary site of HCoV infection. In contrast, FK506 displayed cell-type specific effects, strongly affecting CoV replication in Huh7.5 and HEK293, but inconsistently and less pronounced in phBECs.

Diversity of gammacoronaviruses and deltacoronaviruses in wild birds and poultry in Russia.

Coronaviruses of the genera Gammacoronavirus and Deltacoronavirus are globally widespread and circulate primarily in wild and domestic birds. Prior studies have established frequently occurring crossover events from avian to mammalian reservoirs. However, there is limited understanding of the diversity and geographical distribution of coronaviruses among birds. In this study, the surveillance of coronaviruses in birds in Russia during 2020 revealed the presence of coronaviruses in 12% of samples from birds. Targeted NGS approach was used for the evaluation of genetic diversity based on RdRp gene. While gammacoronviruses were found in both wild birds and poultry, deltacoronaviruses were found in wild birds only and represent the first detections for Russia. A number of cases with the simultaneous detection of gamma- and deltacoronaviruses in one bird was reported. The results of this study highlight the importance of further research concerning the spread and diversity of coronaviruses among birds within and migrating throughout the territory of Russia across the globe.

Host Kinase CSNK2 is a Target for Inhibition of Pathogenic SARS-like ß-Coronaviruses.

Inhibition of the protein kinase CSNK2 with any of 30 specific and selective inhibitors representing different chemotypes, blocked replication of pathogenic human, bat, and murine ß-coronaviruses. The potency of in-cell CSNK2A target engagement across the set of inhibitors correlated with antiviral activity and genetic knockdown confirmed the essential role of the CSNK2 holoenzyme in ß-coronavirus replication. Spike protein endocytosis was blocked by CSNK2A inhibition, indicating that antiviral activity was due in part to a suppression of viral entry. CSNK2A inhibition may be a viable target for the development of anti-SARS-like ß-coronavirus drugs.

Topoisomerase 3b is dispensable for replication of a positive-sense RNA virus--murine coronavirus.

A recent study demonstrated that a DNA-RNA dual-activity topoisomerase complex, TOP3B-TDRD3, is required for normal replication of positive-sense RNA viruses, including several human flaviviruses and coronaviruses; and the authors proposed that TOP3B is a target of antiviral drugs. Here we examined this hypothesis by investigating whether inactivation of Top3b can inhibit the replication of a mouse coronavirus, MHV, using cell lines and mice that are inactivated of Top3b or Tdrd3. We found that Top3b-KO or Tdrd3-KO cell lines generated by different CRISPR-CAS9 guide RNAs have variable effects on MHV replication. In addition, we did not find significant changes of MHV replication in brains or lungs in Top3B-KO mice. Moreover, immunostaining showed that Top3b proteins are not co-localized with MHV replication complexes but rather, localized in stress granules in the MHV-infected cells. Our results suggest that Top3b does not have a universal role in promoting replication of positive-sense RNA virus, and cautions should be taken when targeting it to develop anti-viral drugs.

Sex-biased expression of the TLR7 gene in severe COVID-19 patients: Insights from transcriptomics and epigenomics.

There is abundant epidemiological data indicating that the incidence of severe cases of coronavirus disease (COVID-19) is significantly higher in males than females worldwide. Moreover, genetic variation at the X-chromosome linked TLR7 gene has been associated with COVID-19 severity. It has been suggested that the sex-biased incidence of COVID-19 might be related to the fact that TLR7 escapes X-chromosome inactivation during early embryogenesis in females, thus encoding a doble dose of its gene product compared to males. We analyzed TLR7 expression in two acute phase cohorts of COVID-19 patients that used two different technological platforms, one of them in a multi-tissue context including saliva, nasal, and blood samples, and a third cohort that included different post-infection timepoints of long-COVID-19 patients. We additionally explored methylation patterns of TLR7 using epigenomic data from an independent cohort of COVID-19 patients stratified by severity and sex. In line with genome-wide association studies, we provide supportive evidence indicating that TLR7 has altered CpG methylation patterns and it is consistently downregulated in males compared to females in the most severe cases of COVID-19.

Myotis fimbriatus Virome, a Window to Virus Diversity and Evolution in the Genus Myotis.

Significant efforts have been made to characterize viral diversity in bats from China. Many of these studies were prospective and focused mainly on Rhinolophus bats that could be related to zoonotic events. However, other species of bats that are part of ecosystems identified as virus diversity hotspots have not been studied in-depth. We analyzed the virome of a group of Myotis fimbriatus bats collected from the Yunnan Province during 2020. The virome of M. fimbriatus revealed the presence of families of pathogenic viruses such as Coronavirus, Astrovirus, Mastadenovirus, and Picornavirus, among others. The viral sequences identified in M. fimbriatus were characterized by significant divergence from other known viral sequences of bat origin. Complex phylogenetic landscapes implying a tendency of co-specificity and relationships with viruses from other mammals characterize these groups. The most prevalent and abundant virus in M. fimbriatus individuals was an alphacoronavirus. The genome of this virus shows evidence of recombination and is likely the product of ancestral host-switch. The close phylogenetic and ecological relationship of some species of the Myotis genus in China may have played an important role in the emergence of this alphacoronavirus.

Detection of subgenomic mRNA from endemic human coronavirus OC43 and NL63 compared to viral genomic loads, single virus detection and clinical manifestations in children with respiratory tract infections.

BACKGROUND: The importance of endemic human coronavirus (HCoV) in children has been insufficiently elucidated upon. Our aims were to develop subgenomic (sg) mRNA tests for HCoV species OC43 and NL63, and to evaluate the relationships to HCoV genomic loads, single HCoV detections and clinical manifestations. METHODS: We have used an 11-yearlong cohort study of children admitted with respiratory tract infection (RTI) and hospital controls. Nasopharyngeal aspirates were analyzed for HCoV subtypes OC43 and NL63 with in-house diagnostic PCR. Positive samples were tested with newly developed real-time PCRs targeting sg mRNA coding for the nucleocapsid protein. RESULTS: OC43 sg mRNA was detected in 86% (105/122) of available OC43-positive samples in the RTI group, and in 63% (12/19) of control samples. NL63 sg mRNA was detected in 72% (71/98) and 71% (12/17) of available NL63-positive patient and control samples, respectively. In RTI samples, sg mRNA detection was strongly associated with a Ct value <32 in both diagnostic PCR tests (OC43: OR = 54, 95% CI [6.8-428]; NL63: OR = 42, 95% CI [9.0-198]) and single NL63 detections (OR = 6.9, 95% CI [1.5-32]). Comparing RTI and controls, only OC43 was associated with RTI when adjusted for age (aOR = 3.2, 95% CI [1.1-9.4]). CONCLUSION: We found strong associations between OC43 and NL63 sg mRNA and high viral genomic loads. sg mRNA for OC43 was associated with RTI. The association between sg mRNA and clinical manifestations needs further evaluation.

Viral Nucleases from Herpesviruses and Coronavirus in Recombination and Proofreading: Potential Targets for Antiviral Drug Discovery.

In this review, we explore recombination in two very different virus families that have become major threats to human health. The Herpesviridae are a large family of pathogenic double-stranded DNA viruses involved in a range of diseases affecting both people and animals. Coronaviridae are positive-strand RNA viruses (CoVs) that have also become major threats to global health and economic stability, especially in the last two decades. Despite many differences, such as the make-up of their genetic material (DNA vs. RNA) and overall mechanisms of genome replication, both human herpes viruses (HHVs) and CoVs have evolved to rely heavily on recombination for viral genome replication, adaptation to new hosts and evasion of host immune regulation. In this review, we will focus on the roles of three viral exonucleases: two HHV exonucleases (alkaline nuclease and PolExo) and one CoV exonuclease (ExoN). We will review the roles of these three nucleases in their respective life cycles and discuss the state of drug discovery efforts against these targets.

Origins and Evolution of Seasonal Human Coronaviruses.

Four seasonal human coronaviruses (sHCoVs) are endemic globally (229E, NL63, OC43, and HKU1), accounting for 5-30% of human respiratory infections. However, the epidemiology and evolution of these CoVs remain understudied due to their association with mild symptomatology. Using a multigene and complete genome analysis approach, we find the evolutionary histories of sHCoVs to be highly complex, owing to frequent recombination of CoVs including within and between sHCoVs, and uncertain, due to the under sampling of non-human viruses. The recombination rate was highest for 229E and OC43 whereas substitutions per recombination event were highest in NL63 and HKU1. Depending on the gene studied, OC43 may have ungulate, canine, or rabbit CoV ancestors. 229E may have origins in a bat, camel, or an unsampled intermediate host. HKU1 had the earliest common ancestor (1809-1899) but fell into two distinct clades (genotypes A and B), possibly representing two independent transmission events from murine-origin CoVs that appear to be a single introduction due to large gaps in the sampling of CoVs in animals. In fact, genotype B was genetically more diverse than all the other sHCoVs. Finally, we found shared amino acid substitutions in multiple proteins along the non-human to sHCoV host-jump branches. The complex evolution of CoVs and their frequent host switches could benefit from continued surveillance of CoVs across non-human hosts.