Positive selection, genetic recombination, and intra-host evolution in novel equine coronavirus genomes and other members of the Embecovirus subgenus.

There are several examples of coronaviruses in the Betacoronavirus subgenus Embecovirus that have jumped from an animal to the human host. Studying how evolutionary factors shape coronaviruses in non-human hosts may provide insight into the coronavirus host-switching potential. Equids, such as horses and donkeys, are susceptible to equine coronaviruses (ECoVs). With increased testing prevalence, several ECoV genome sequences have become available for molecular evolutionary analyses, especially those from the United States of America (USA). To date, no analyses have been performed to characterize evolution within coding regions of the ECoV genome. Here, we obtain and describe four new ECoV genome sequences from infected equines from across the USA presenting clinical symptoms of ECoV, and infer ECoV-specific and Embecovirus-wide patterns of molecular evolution. Within two of the four data sets analyzed, we find evidence of intra-host evolution within the nucleocapsid (N) gene, suggestive of quasispecies development. We also identify 12 putative genetic recombination events within the ECoV genome, 11 of which fall in ORF1ab. Finally, we infer and compare sites subject to positive selection on the ancestral branch of each major Embecovirus member clade. Specifically, for the two currently identified human coronavirus (HCoV) embecoviruses that have spilled from animals to humans (HCoV-OC43 and HCoV-HKU1), we find that there are 42 and 2 such sites, respectively, perhaps reflective of the more complex ancestral evolutionary history of HCoV-OC43, which involves several different animal hosts.IMPORTANCEThe Betacoronavirus subgenus Embecovirus contains coronaviruses that not only pose a health threat to animals and humans, but also have jumped from animal to human host. Equids, such as horses and donkeys are susceptible to equine coronavirus (ECoV) infections. No studies have systematically examined evolutionary patterns within ECoV genomes. Our study addresses this gap and provides insight into intra-host ECoV evolution from infected horses. Further, we identify and report natural selection pattern differences between two embecoviruses that have jumped from animals to humans [human coronavirus OC43 and HKU1 (HCoV-OC43 and HCoV-HKU1, respectively)], and hypothesize that the differences observed may be due to the different animal host(s) that each virus circulated in prior to its jump into humans. Finally, we contribute four novel, high-quality ECoV genomes to the scientific community.

Molecular detection using hybridization capture and next-generation sequencing reveals cross-species transmission of feline coronavirus type-1 between a domestic cat and a captive wild felid.

Feline coronavirus (FCoV) infection normally causes mild or subclinical signs and is common in domestic cats. However, in some cats, FCoV infection can also lead to the development of feline infectious peritonitis (FIP)-a typically lethal disease. FCoV has two serotypes or genotypes, FCoV-1 and FCoV-2, both of which can cause FIP. The main difference between the genotypes is the viral spike (S) protein that determines tropism and pathogenicity, crucial mechanisms in the development of FIP. Subclinical infection and FIP have both been reported in wild felids, including in threatened species. Due to the high genetic variability of the S gene and the technical challenges to sequencing it, detection and characterization of FCoV in wild felids have mainly centered on other more conserved genes. Therefore, the genotype causing FIP in most wild felids remains unknown. Here, we report a retrospective molecular epidemiological investigation of FCoV in a zoological institution in the U.Ss. In 2008, a domestic cat (Felis catus) and a Pallas' cat (Otocolobus manul) sharing the same room succumbed to FIP. Using in situ hybridization, we detected FCoV RNA in different tissues of both felids. Using hybridization capture and next-generation sequencing, we detected, sequenced, and characterized the whole genome of the FCoV infecting both felids. Our data show for the first time that FCoV-1 can be transmitted between domestic and wild felids and extends the known host range of FCoV-1. Our findings highlight the importance of identifying the genotype causing FIP, to develop effective control measures. IMPORTANCE: Feline coronavirus (FCoV) is highly prevalent in domestic cats worldwide and has also been reported in wild felids, including endangered species, in which it has caused substantial population declines. Characterizing the genetic diversity of FCoV is crucial due to recent reports of novel pathogenic recombinant variants causing high mortality in feral cats in Cyprus. In this retrospective molecular epidemiology study, we used archived samples collected in a zoological institution in the U.S. in which a domestic and a wild felid succumbed to FCoV. Using hybridization capture (HC) and next-generation sequencing, we show for the first time that FCoV can be naturally transmitted between domestic and wild felids. We demonstrate the efficacy of HC for detecting and sequencing the whole genome of FCoV, which is essential to characterize its different genotypes.

Intranasal Vaccination with Recombinant TLR2-Active Outer Membrane Vesicles Containing Sequential M2e Epitopes Protects against Lethal Influenza a Challenge.

Influenza is a highly contagious respiratory disease, resulting in an estimated 3 to 5 million cases of severe illness annually. While most influenza vaccines are administered parenterally via injection, one shortcoming is that they do not generate a strong immune response at the site of infection, which can become important in a pandemic. Intranasal vaccines can generate both local and systemic protective immune responses, can reduce costs, and enhance ease of administration. Previous studies showed that parenterally administered outer membrane vesicles (OMVs) that carry sequences of the M2e protein (OMV-M2e) protect against influenza A/PR8 challenge in mice and ferrets. In the current study, we measured the effectiveness of the intranasal route of the OMV-M2e vaccine against the influenza A/PR8 strain in mice. We observed high anti-M2e IgG and IgA titers post-challenge in mice vaccinated intranasally with OMV-M2e. In addition, we observed a Th1/Tc1 bias in the vaccinated mice, and an increased Th17/Tc17 response, both of which correlated with survival to A/PR8 challenge and significantly lower lung viral titers. We conclude that the intranasal-route administration of the OMV-M2e vaccine is a promising approach toward generating protection against influenza A as it leads to an increased proinflammatory immune response correlating with survival to viral challenge.

Host Cell Proteases Involved in Human Respiratory Viral Infections and Their Inhibitors: A Review.

Viral tropism is most commonly linked to receptor use, but host cell protease use can be a notable factor in susceptibility to infection. Here we review the use of host cell proteases by human viruses, focusing on those with primarily respiratory tropism, particularly SARS-CoV-2. We first describe the various classes of proteases present in the respiratory tract, as well as elsewhere in the body, and incorporate the targeting of these proteases as therapeutic drugs for use in humans. Host cell proteases are also linked to the systemic spread of viruses and play important roles outside of the respiratory tract; therefore, we address how proteases affect viruses across the spectrum of infections that can occur in humans, intending to understand the extrapulmonary spread of SARS-CoV-2.

SARS-COV-2 SEROSURVEY OF HEALTHY, PRIVATELY OWNED CATS PRESENTING TO A NEW YORK CITY ANIMAL HOSPITAL IN THE EARLY PHASE OF THE COVID-19 PANDEMIC (2020-2021).

SARS-CoV-2, the cause of the ongoing COVID-19 pandemic, not only infects humans but is also known to infect various species, including domestic and wild animals. While many species have been identified as susceptible to SARS-CoV-2, there are limited studies on the prevalence of SARS-CoV-2 in animals. Both domestic and non-domestic cats are now established to be susceptible to infection by SARS-CoV-2. While serious disease in cats may occur in some instances, the majority of infections appear to be subclinical. Differing prevalence data for SARS-CoV-2 infection of cats have been reported, and are highly context-dependent. Here, we report a retrospective serological survey of cats presented to an animal practice in New York City, located in close proximity to a large medical center that treated the first wave of COVID-19 patients in the U.S. in the Spring of 2020. We sampled 79, mostly indoor, cats between June 2020 to May 2021, the early part of which time the community was under a strict public health "lock-down". Using a highly sensitive and specific fluorescent bead-based multiplex assay, we found an overall prevalence of 13/79 (16%) serologically-positive animals for the study period; however, cats sampled in the Fall of 2020 had a confirmed positive prevalence of 44%. For SARS-CoV-2 seropositive cats, we performed viral neutralization test with live SARS-CoV-2 to additionally confirm presence of SARS-CoV-2 specific antibodies. Of the thirteen seropositive cats, 7/13 (54%) were also positive by virus neutralization, and two of seropositive cats had previously documented respiratory signs, with high neutralization titers of 1/1024 and 1/4096; overall however, there was no statistically significant association of SARS-CoV-2 seropositivity with respiratory signs, or with breed, sex or age of the animals. Follow up sampling of cats showed that positive serological titers were maintained over time. In comparison, we found an overall confirmed positive prevalence of 51% for feline coronavirus (FCoV), an endemic virus of cats, with 30% confirmed negative for FCoV. We demonstrate the impact of SARS-CoV-2 in a defined feline population during the first wave of SARS-CoV-2 infection of humans, and suggest that human-cat transmission was substantial in our study group. Our study provide a new context for SARS-CoV-2 transmission events across species.

A Mechanistic Understanding of the Modes of Ca2+ Ion Binding to the SARS-CoV-1 Fusion Peptide and Their Role in the Dynamics of Host Membrane Penetration.

The SARS-CoV-1 spike glycoprotein contains a fusion peptide (FP) segment that mediates the fusion of the viral and host cell membranes. Calcium ions are thought to position the FP optimally for membrane insertion by interacting with negatively charged residues in this segment (E801, D802, D812, E821, D825, and D830); however, which residues bind to calcium and in what combinations supportive of membrane insertion are unknown. Using biological assays and molecular dynamics studies, we have determined the functional configurations of FP-Ca2+ binding that likely promote membrane insertion. We first individually mutated the negatively charged residues in the SARS CoV-1 FP to assay their roles in cell entry and syncytia formation, finding that charge loss in the D802A or D830A mutants greatly reduced syncytia formation and pseudoparticle transduction of VeroE6 cells. Interestingly, one mutation (D812A) led to a modest increase in cell transduction, further indicating that FP function likely depends on calcium binding at specific residues and in specific combinations. To interpret these results mechanistically and identify specific modes of FP-Ca2+ binding that modulate membrane insertion, we performed molecular dynamics simulations of the SARS-CoV-1 FP and Ca2+ions. The preferred residue pairs for Ca2+ binding we identified (E801/D802, E801/D830, and D812/E821) include the two residues found to be essential for S function in our biological studies (D802 and D830). The three preferred Ca2+ binding pairs were also predicted to promote FP membrane insertion. We also identified a Ca2+ binding pair (E821/D825) predicted to inhibit FP membrane insertion. We then carried out simulations in the presence of membranes and found that binding of Ca2+ to SARS-CoV-1 FP residue pairs E801/D802 and D812/E821 facilitates membrane insertion by enabling the peptide to adopt conformations that shield the negative charges of the FP to reduce repulsion by the membrane phospholipid headgroups. This calcium binding mode also optimally positions the hydrophobic LLF region of the FP for membrane penetration. Conversely, Ca2+ binding to the FP E801/D802 and D821/D825 pairs eliminates the negative charge screening and instead creates a repulsive negative charge that hinders membrane penetration of the LLF motif. These computational results, taken together with our biological studies, provide an improved and nuanced mechanistic understanding of the dymanics of SARS-CoV-1 calcium binding and their potential effects on host cell entry.

Analysis of the molecular determinants for furin cleavage of the spike protein S1/S2 site in defined strains of the prototype coronavirus murine hepatitis virus (MHV).

We analyzed the spike protein S1/S2 cleavage of selected strains of a prototype coronavirus, mouse hepatitis virus (MHV) by the cellular protease furin, in order to understand the structural requirements underlying the sequence selectivity of the scissile segment. The probability of cleavage of selected MHV strains was first evaluated from furin cleavage scores predicted by the ProP computer software, and then cleavage was measured experimentally with a fluorogenic peptide cleavage assay consisting of S1/S2 peptide mimics and purified furin. We found that in vitro cleavability varied across MHV strains in line with predicted results-but with the notable exception of MHV-A59, which was not cleaved despite a high score predicted for its sequence. Using the known X-Ray structure of furin in complex with a substrate-like inhibitor as an initial structural reference, we carried out molecular dynamics (MD) simulations to learn the modes of binding of the peptides in the furin active site, and the suitability of the complex for initiation of the enzymatic cleavage. We identified the 3D structural requirements of the furin active site configuration that enable bound peptides to undergo cleavage, and the way in which the various strains tested experimentally are fulfilling these requirements. We find that despite some flexibility in the organization of the peptide bound to the active site of the enzyme, the presence of a histidine at P2 of MHV-A59 fails to properly orient the sidechain of His194 of the furin catalytic triad and therefore produces a distortion that renders the peptide/complex structural configuration in the active site incompatible with requirements for cleavage initiation. The Ser/Thr in P1 of MHV-2 and MHV-S has a similar effect of distorting the conformation of the furin active site residues produced by the elimination of the canonical salt-bridge formed by arginine in P1 position. This work informs a study of coronavirus infection and pathogenesis with respect to the function of the viral spike protein, and suggests an important process of viral adaptation and evolution within the spike S1/S2 structural loop.

Backyard zoonoses: The roles of companion animals and peri-domestic wildlife.

The spillover of human infectious diseases from animal reservoirs is now well appreciated. However, societal and climate-related changes are affecting the dynamics of such interfaces. In addition to the disruption of traditional wildlife habitats, in part because of climate change and human demographics and behavior, there is an increasing zoonotic disease risk from companion animals. This includes such factors as the awareness of animals kept as domestic pets and increasing populations of free-ranging animals in peri-domestic environments. This review presents background and commentary focusing on companion and peri-domestic animals as disease risk for humans, taking into account the human-animal interface and population dynamics between the animals themselves.

Inhibiting Glutamine Metabolism Blocks Coronavirus Replication in Mammalian Cells.

Developing therapeutic strategies against COVID-19 has gained widespread interest given the likelihood that new viral variants will continue to emerge. Here we describe one potential therapeutic strategy which involves targeting members of the glutaminase family of mitochondrial metabolic enzymes (GLS and GLS2), which catalyze the first step in glutamine metabolism, the hydrolysis of glutamine to glutamate. We show three examples where GLS expression increases during coronavirus infection of host cells, and another in which GLS2 is upregulated. The viruses hijack the metabolic machinery responsible for glutamine metabolism to generate the building blocks for biosynthetic processes and satisfy the bioenergetic requirements demanded by the 'glutamine addiction' of virus-infected host cells. We demonstrate how genetic silencing of glutaminase enzymes reduces coronavirus infection and that newer members of two classes of small molecule allosteric inhibitors targeting these enzymes, designated as SU1, a pan-GLS/GLS2 inhibitor, and UP4, which is specific for GLS, block viral replication in mammalian epithelial cells. Overall, these findings highlight the importance of glutamine metabolism for coronavirus replication in human cells and show that glutaminase inhibitors can block coronavirus infection and thereby may represent a novel class of anti-viral drug candidates. Teaser: Inhibitors targeting glutaminase enzymes block coronavirus replication and may represent a new class of anti-viral drugs.

Age-dependent acquisition of pathogenicity by SARS-CoV-2 Omicron BA.5.

Pathology studies of SARS-CoV-2 Omicron variants of concern (VOC) are challenged by the lack of pathogenic animal models. While Omicron BA.1 and BA.2 replicate in K18-hACE2 transgenic mice, they cause minimal to negligible morbidity and mortality, and less is known about more recent Omicron VOC. Here, we show that in contrast to Omicron BA.1, BA.5-infected mice exhibited high levels of morbidity and mortality, correlating with higher early viral loads. Neither Omicron BA.1 nor BA.5 replicated in brains, unlike most prior VOC. Only Omicron BA.5-infected mice exhibited substantial weight loss, high pathology scores in lungs, and high levels of inflammatory cells and cytokines in bronchoalveolar lavage fluid, and 5- to 8-month-old mice exhibited 100% fatality. These results identify a rodent model for pathogenesis or antiviral countermeasure studies for circulating SARS-CoV-2 Omicron BA.5. Further, differences in morbidity and mortality between Omicron BA.1 and BA.5 provide a model for understanding viral determinants of pathogenicity.

Intra-host variation in the spike S1/S2 region of a feline coronavirus type-1 in a cat with persistent infection.

Feline coronavirus type 1 (FCoV-1) is widely known for causing feline infectious peritonitis (FIP), a systemic infection that is often fatal, with the virus known as the FIPV biotype. However, subclinical disease also occurs, in which cats may not show signs and intermittently shed the virus, including in feces, possibly for long periods of time. This virus is known as the FECV biotype. Progression of FECV to FIPV has been linked to several genomic changes, however a specific region of the viral spike protein at the interface of the spike S1 and S2 domains has been especially implicated. In this study, we followed a cat (#576) for six years from 2017, at which time FCoV-1 was detected in feces and conjunctival swabs, until 2022, when the animal was euthanized based on a diagnosis of alimentary small cell lymphoma. Over this time period, the cat was clinically diagnosed with inflammatory bowel disease and chronic rhinitis, and cardiac problems were also suspected. Using hybridization capture targeting the spike (S) gene of FCoV followed by next-generation sequencing, we screened 27 clinical samples. We detected FCoV-1 in 4 samples taken in 2017 (intestine and nasal tissue, feces, and conjunctiva), and 3 samples taken in 2022 (feces, and intestinal and heart tissue), but not in fecal samples taken in 2019 and 2020. Next, we focused on the S1/S2 region within S, which contains the furin cleavage site (FCS), a key regulator of viral transmission and pathogenesis. We show that the FCoV-1 variants obtained from feces in 2017 and 2022 were identical, while the ones from conjunctiva (2017), heart (2022), and intestine (2017 and 2022) were distinct. Sequence comparison of all the variants obtained showed that most of the non-synonymous changes in the S1/S2 region occur within the FCS. In the heart, we found two variants that differed by a single nucleotide, resulting in distinct FCS motifs that differ in one amino acid. It is predicted that one of these FCS motifs will down-regulate spike cleavability. The variant from the conjunctiva (2017) had a 6-nucleotide in-frame insertion that resulted in a longer and more exposed S1/S2 loop, which is predicted to be more accessible to the furin protease. Our studies indicate that FCoV-1 can independently persist in the gastrointestinal tract and heart of a cat over a long period of time without evidence of typical FIP signs, with intermittent viral shedding from the gastrointestinal and respiratory tracts.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic feline coronavirus phenotypes.

Feline coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed feline enteric coronavirus [FECV]), with around 12 per cent developing into deadly feline infectious peritonitis (FIP; feline infectious peritonitis virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV- and FECV-specific signals of positive selection. We analyzed full-length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site (FCS) and the other within the fusion domain of Spike. We also found fifteen sites subject to positive selection associated with FIPV within Spike, eleven of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were fourteen sites (twelve novel sites) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 FCS and adjacent to C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype and included twenty-four positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV-wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that it is unlikely to be one singular 'switch' mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Interactions of SARS-CoV-2 and MERS-CoV fusion peptides measured using single-molecule force methods.

We address the challenge of understanding how hydrophobic interactions are encoded by fusion peptide (FP) sequences within coronavirus (CoV) spike proteins. Within the FPs of severe acute respiratory syndrome CoV 2 and Middle East respiratory syndrome CoV (MERS-CoV), a largely conserved peptide sequence called FP1 (SFIEDLLFNK and SAIEDLLFDK in SARS-2 and MERS, respectively) has been proposed to play a key role in encoding hydrophobic interactions that drive viral-host cell membrane fusion. Although a non-polar triad (Leu-Leu-Phe (LLF)) is common to both FP1 sequences, and thought to dominate the encoding of hydrophobic interactions, FP1 from SARS-2 and MERS differ in two residues (Phe 2 versus Ala 2 and Asn 9 versus Asp 9, respectively). Here we explore whether single-molecule force measurements can quantify hydrophobic interactions encoded by FP1 sequences, and then ask whether sequence variations between FP1 from SARS-2 and MERS lead to significant differences in hydrophobic interactions. We find that both SARS-2 and MERS wild-type FP1 generate measurable hydrophobic interactions at the single-molecule level, but that SARS-2 FP1 encodes a substantially stronger hydrophobic interaction than its MERS counterpart (1.91 ± 0.03 nN versus 0.68 ± 0.03 nN, respectively). By performing force measurements with FP1 sequences with single amino acid substitutions, we determine that a single-residue mutation (Phe 2 versus Ala 2) causes the almost threefold difference in the hydrophobic interaction strength generated by the FP1 of SARS-2 versus MERS, despite the presence of LLF in both sequences. Infrared spectroscopy and circular dichroism measurements support the proposal that the outsized influence of Phe 2 versus Ala 2 on the hydrophobic interaction arises from variation in the secondary structure adopted by FP1. Overall, these insights reveal how single-residue diversity in viral FPs, including FP1 of SARS-CoV-2 and MERS-CoV, can lead to substantial changes in intermolecular interactions proposed to play a key role in viral fusion, and hint at strategies for regulating hydrophobic interactions of peptides in a range of contexts.

Analysis of the molecular determinants for furin cleavage of the spike protein S1/S2 site in defined strains of the prototype coronavirus murine hepatitis virus (MHV).

We have analyzed the spike protein S1/S2 cleavage site of selected strains of MHV by the cellular protease furin, in order to understand the structural requirements underlying the sequence selectivity of the scissile segment. The probability of cleavage of the various MHV strains was first evaluated from furin cleavage scores predicted by the ProP computer software, and then cleavage was measured experimentally with a fluorogenic peptide cleavage assay consisting of S1/S2 peptide mimics and purified furin. We found that in vitro cleavability varied across MHV strains in line with predicted results-but with the notable exception of MHV-A59, which was not cleaved despite a high score predicted for its sequence. Using the known X-Ray structure of furin in complex with a substrate-like inhibitor as an initial structural reference, we carried out molecular dynamics (MD) simulations to learn the modes of binding of the peptides in the furin active site, and the suitability of the complex for initiation of the enzymatic cleavage. We thus identified the 3D structural requirements of the furin active site configuration that enable bound peptides to undergo cleavage, and the way in which the various strains tested experimentally are fulfilling these requirements. We find that despite some flexibility in the organization of the peptide bound to the active site of the enzyme, the presence of a histidine at P2 of MHV-A59 fails to properly orient the sidechain of His194 of the furin catalytic triad and therefore produces a distortion that renders the peptide/complex structural configuration in the active site incompatible with requirements for cleavage initiation. The Ser/Thr in P1 of MHV-2 and MHV-S has a similar effect of distorting the conformation of the furin active site residues produced by the elimination of the canonical salt-bridge formed by arginine in P1 position. This work informs a study of coronavirus infection and pathogenesis with respect to the function of the viral spike protein, and suggests an important process of viral adaptation and evolution within the spike S1/S2 structural loop.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic Feline Coronavirus phenotypes.

Feline Coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed Feline Enteric Coronavirus [FECV]), with around 12% developing into deadly Feline Infectious Peritonitis (FIP; Feline Infectious Peritonitis Virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV and FECV specific signals of positive selection. We analyzed full length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site, and the other within the fusion domain of Spike. We also found 15 sites subject to positive selection associated with FIPV within Spike, 11 of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were 14 sites (12 novel) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 furin cleavage site and adjacent C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype, and included 24 positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that is unlikely to be one singular "switch" mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

The Omicron Variant BA.1.1 Presents a Lower Pathogenicity than B.1 D614G and Delta Variants in a Feline Model of SARS-CoV-2 Infection.

Omicron (B.1.1.529) is the most recent SARS-CoV-2 variant of concern, which emerged in late 2021 and rapidly achieved global predominance by early 2022. In this study, we compared the infection dynamics, tissue tropism, and pathogenesis and pathogenicity of SARS-CoV-2 D614G (B.1), Delta (B.1.617.2), and Omicron BA.1.1 (B.1.1.529) variants in a highly susceptible feline model of infection. Although D614G- and Delta-inoculated cats became lethargic and showed increased body temperatures between days 1 and 3 postinfection (pi), Omicron-inoculated cats remained subclinical and, similar to control animals, gained weight throughout the 14-day experimental period. Intranasal inoculation of cats with D614G- and the Delta variants resulted in high infectious virus shedding in nasal secretions (up to 6.3 log10 TCID50.Ml-1), whereas strikingly lower level of viruses shedding (<3.1 log10 TCID50.Ml-1) was observed in Omicron-inoculated animals. In addition, tissue distribution of the Omicron variant was markedly reduced in comparison to the D614G and Delta variants, as evidenced by lower in situ viral RNA detection, in situ viral immunofluorescence staining, and viral loads in tissues on days 3, 5, and 14 pi. Nasal turbinate, trachea, and lung were the main-but not the only-sites of replication for all three viral variants. However, only scarce virus staining and lower viral titers suggest lower levels of viral replication in tissues from Omicron-infected animals. Notably, while D614G- and Delta-inoculated cats presented pneumonia, histologic examination of the lungs from Omicron-infected cats revealed mild to modest inflammation. Together, these results demonstrate that the Omicron variant BA.1.1 is less pathogenic than D614G and Delta variants in a highly susceptible feline model. IMPORTANCE The SARS-CoV-2 Omicron (B.1.1.529) variant of concern emerged in South Africa late in 2021 and rapidly spread across the world causing a significant increase in the number of infections. Importantly, this variant was also associated with an increased risk of reinfections. However, the number of hospitalizations and deaths due to COVID-19 did not follow the same trends. These early observations suggested effective protection conferred by immunizations and/or overall lower virulence of the highly mutated variant virus. In this study we present novel evidence demonstrating that the Omicron BA.1.1 variant of concern presents a lower pathogenicity when compared to D614G- or Delta variants in cats. Clinical, virological, and pathological evaluations revealed lower disease severity, viral replication, and lung pathology in Omicron-infected cats when compared with D614G and Delta variant inoculated animals, confirming that Omicron BA.1.1 is less pathogenic in a highly susceptible feline model of infection.

Intrinsic furin-mediated cleavability of the spike S1/S2 site from SARS-CoV-2 variant B.1.1.529 (Omicron).

The ability of SARS-CoV-2 to be primed for viral entry by the host cell protease furin has become one of the most investigated of the numerous transmission and pathogenicity features of the virus. SARS-CoV-2 The variant B.1.1.529 (Omicron) emerged in late 2020 and has continued to evolve and is now present in several distinct sub-variants. Here, we analyzed the "furin cleavage site" of the spike protein of SARS-CoV-2 B.1.1.529 (Omicron variant) in vitro, to assess the role of two key mutations (spike, N679K and P681H) that are common across all subvariants compared to the ancestral B.1 virus and other notable lineages. We observed significantly increased intrinsic cleavability with furin compared to an original B lineage virus (Wuhan-Hu1), as well as to two variants, B.1.1.7 (Alpha) and B.1.617 (Delta) that subsequently had wide circulation. Increased furin-mediated cleavage was attributed to the N679K mutation, which lies outside the conventional furin binding pocket. Our findings suggest that B.1.1.529 (Omicron variant) has gained genetic features linked to intrinsic furin cleavability, in line with its evolution within the population as the COVID-19 pandemic has proceeded.

The Omicron variant BA.1.1 presents a lower pathogenicity than B.1 D614G and Delta variants in a feline model of SARS-CoV-2 infection.

Omicron (B.1.1.529) is the most recent SARS-CoV-2 variant of concern (VOC), which emerged in late 2021 and rapidly achieved global predominance in early 2022. In this study, we compared the infection dynamics, tissue tropism and pathogenesis and pathogenicity of SARS-CoV-2 D614G (B.1), Delta (B.1.617.2) and Omicron BA.1.1 sublineage (B.1.1.529) variants in a highly susceptible feline model of infection. While D614G- and Delta-inoculated cats became lethargic, and showed increased body temperatures between days 1 and 3 post-infection (pi), Omicron-inoculated cats remained subclinical and, similar to control animals, gained weight throughout the 14-day experimental period. Intranasal inoculation of cats with D614G- and the Delta variants resulted in high infectious virus shedding in nasal secretions (up to 6.3 log10 TCID 50 .ml -1 ), whereas strikingly lower level of viruses shedding (<3.1 log10 TCID 50 .ml -1 ) was observed in Omicron-inoculated animals. In addition, tissue distribution of the Omicron variant was markedly reduced in comparison to the D614G and Delta variants, as evidenced by in situ viral RNA detection, in situ immunofluorescence, and quantification of viral loads in tissues on days 3, 5, and 14 pi. Nasal turbinate, trachea, and lung were the main - but not the only - sites of replication for all three viral variants. However, only scarce virus staining and lower viral titers suggest lower levels of viral replication in tissues from Omicron-infected animals. Notably, while D614G- and Delta-inoculated cats had severe pneumonia, histologic examination of the lungs from Omicron-infected cats revealed mild to modest inflammation. Together, these results demonstrate that the Omicron variant BA.1.1 is less pathogenic than D614G and Delta variants in a highly susceptible feline model. Author Summary: The SARS-CoV-2 Omicron (B.1.1.529) variant of concern (VOC) emerged in South Africa late in 2021 and rapidly spread across the world causing a significant increase in the number of infections. Importantly, this variant was also associated with an increased risk of reinfections. However, the number of hospitalizations and deaths due to COVID-19 did not follow the same trends. These early observations, suggested effective protection conferred by immunizations and/or overall lower virulence of the highly mutated variant virus. In this study we present novel evidence demonstrating that the Omicron BA.1.1 variant of concern (VOC) presents a lower pathogenicity when compared to D614G- or Delta variants in cats. Clinical, virological and pathological evaluations revealed lower disease severity, viral replication and lung pathology in Omicron-infected cats when compared to D614G and Delta variant inoculated animals, confirming that Omicron BA.1.1 is less pathogenic in a highly susceptible feline model of infection.

Spike Protein Cleavage-Activation in the Context of the SARS-CoV-2 P681R Mutation: an Analysis from Its First Appearance in Lineage A.23.1 Identified in Uganda.

Based on its predicted ability to affect transmissibility and pathogenesis, surveillance studies have highlighted the role of a specific mutation (P681R) in the S1/S2 furin cleavage site of the SARS-CoV-2 spike protein. Here we analyzed A.23.1, first identified in Uganda, as a P681R-containing virus several months prior to the emergence of B.1.617.2 (Delta variant). We performed assays using peptides mimicking the S1/S2 from A.23.1 and B.1.617 and observed significantly increased cleavability with furin compared to both an original B lineage (Wuhan-Hu1) and B.1.1.7 (Alpha variant). We also performed cell-cell fusion and functional infectivity assays using pseudotyped particles and observed an increase in activity for A.23.1 compared to an original B lineage spike. However, these changes in activity were not reproduced in the B lineage spike bearing only the P681R substitution. Our findings suggest that while A.23.1 has increased furin-mediated cleavage linked to the P681R substitution, this substitution needs to occur on the background of other spike protein changes to enable its functional consequences. IMPORTANCE During the course of the SARS-CoV-2 pandemic, viral variants have emerged that often contain notable mutations in the spike gene. Mutations that encode changes in the spike S1/S2 (furin) activation site have been considered especially impactful. The S1/S2 change from proline to arginine at position 681 (P681R) first emerged in the A.23.1 variant in Uganda, and subsequently occurred in the more widely transmitted Delta variant. We show that the A.23.1 spike is more readily activated by the host cell protease furin, but that this is not reproduced in an original SARS-CoV-2 spike containing the P681R mutation. Changes to the S1/S2 (furin) activation site play a role in SARS-CoV-2 infection and spread, but successful viruses combine these mutations with other less well identified changes, occurring as part of natural selection.