The virus-host interface: Molecular interactions of Alphacoronavirus-1 variants from wild and domestic hosts with mammalian aminopeptidase N

Epigenetic signatures of social status in wild female spotted hyenas (Crocuta crocuta)

In mammalian societies, dominance hierarchies translate into inequalities in health, reproductive performance and survival. DNA methylation is thought to mediate the effects of social status on gene expression and phenotypic outcomes, yet a study of social status-specific DNA methylation profiles in different age classes in a wild social mammal is missing. We tested for social status signatures in DNA methylation profiles in wild female spotted hyenas (Crocuta crocuta), cubs and adults, using non-invasively collected gut epithelium samples. In spotted hyena clans, female social status influences access to resources, foraging behavior, health, reproductive performance and survival. We identified 149 differentially methylated regions between 42 high- and low-ranking female spotted hyenas (cubs and adults). Differentially methylated genes were associated with energy conversion, immune function, glutamate receptor signalling and ion transport. Our results provide evidence that socio-environmental inequalities are reflected at the molecular level in cubs and adults in a wild social mammal.

Detection and characterization of novel luchacoviruses, genus Alphacoronavirus, in saliva and feces of meso-carnivores in the northeastern United States

IMPORTANCE

Several coronaviruses (CoVs) have been detected in domesticated, farmed, and wild meso-carnivores, causing a wide range of diseases and infecting diverse species, highlighting their important but understudied role in the epidemiology of these viruses. Assessing the viral diversity hosted in wildlife species is essential to understand their significance in the cross-species transmission of CoVs. Our focus here was on CoV discovery in meso-carnivores in the Northeast United States as a potential "hotspot" area with high density of humans and urban wildlife. This study identifies novel alphacoronaviruses circulating in multiple free-ranging wild and domestic species in this area and explores their potential epidemiological importance based on regions of the Spike gene, which are relevant for virus-host interactions.

Characterization of CCoV-HuPn-2018 spike protein-mediated viral entry

Canine coronavirus-human pneumonia-2018 (CCoV-HuPn-2018) was recently isolated from a child with pneumonia. This novel human pathogen resulted from cross-species transmission of a canine coronavirus. It has been known that CCoV-HuPn-2018 uses aminopeptidase N (APN) from canines, felines, and porcines, but not humans, as functional receptors for cell entry. The molecular mechanism of cell entry in CCoV-HuPn-2018 remains poorly understood. In this study, we demonstrated that among the nine APN orthologs tested, the APN of the Mexican free-tailed bat could also efficiently support CCoV-HuPn-2018 spike (S) protein-mediated entry, raising the possibility that bats may also be an alternative host epidemiologically important for the transmission of this virus. The glycosylation at residue N747 of canine APN is critical for its receptor activity. The gain of glycosylation at the corresponding residues in human and rabbit APNs converted them to functional receptors for CCoV-HuPn-2018. Interestingly, the CCoV-HuPn-2018 spike protein pseudotyped virus infected multiple human cancer cell lines in a human APN-independent manner, whereas sialic acid appeared to facilitate the entry of the pseudotyped virus into human cancer cells. Moreover, while host cell surface proteases trypsin and TMPRSS2 did not promote the entry of CCoV-HuPn-2018, endosomal proteases cathepsin L and B are required for the entry of CCoV-HuPn-2018 in a pH-dependent manner. IFITMs and LY6E are host restriction factors for the CCoV-HuPn-2018 entry. Our results thus suggest that CCoV-HuPn-2018 has not yet evolved to be an efficient human pathogen. Collectively, this study helps us understand the cell tropism, receptor usage, cross-species transmission, natural reservoir, and pathogenesis of this potential human coronavirus. IMPORTANCE Viral entry is driven by the interaction between the viral spike protein and its specific cellular receptor, which determines cell tropism and host range and is the major constraint to interspecies transmission of coronaviruses. Aminopeptidase N (APN; also called CD13) is a cellular receptor for HCoV-229E, the newly discovered canine coronavirus-human pneumonia-2018 (CCoV-HuPn-2018), and many other animal alphacoronaviruses. We examined the receptor activity of nine APN orthologs and found that CCoV-HuPn-2018 utilizes APN from a broad range of animal species, including bats but not humans, to enter host cells. To our surprise, we found that CCoV-HuPn-2018 spike protein pseudotyped viral particles successfully infected multiple human hepatoma-derived cell lines and a lung cancer cell line, which is independent of the expression of human APN. Our findings thus provide mechanistic insight into the natural hosts and interspecies transmission of CCoV-HuPn-2018-like coronaviruses.

Detection and characterization of novel luchacoviruses, genus Alphacoronavirus, shed in saliva and feces of meso-carnivores in the northeastern United States

Small to mid-sized carnivores, or meso-carnivores, comprise a group of diverse mammals, many of which can adapt to anthropogenically disturbed environments. Wild meso-carnivores living in urban areas may get exposed to or spread pathogens to other species, including stray/feral domestic animals. Several coronaviruses (CoVs) have been detected in domesticated and farmed meso-carnivores, but knowledge of CoVs circulating in free-ranging wild meso-carnivores remains limited. In this study, we analyzed 321 samples collected between 2016 and 2022 from 9 species of free-ranging wild meso-carnivores and stray/feral domestic cats in the northeastern United States. Using a pan-CoV PCR, we screened tissues, feces, and saliva, nasal, and rectal swabs. We detected CoV RNA in fecal and saliva samples of animals in four species: fisher (Pekania pennanti), bobcat (Lynx rufus), red fox (Vulpes vulpes), and domestic cat (Felis catus). Next-generation sequencing revealed that all these viruses belonged to the Luchacovirus subgenus (Alphacoronavirus genus), previously reported only in rodents and lagomorphs (i.e., rabbits). Genetic comparison of the 3'-end of the genome (~12,000bp) revealed that although the viruses detected group with, and have a genetic organization similar to other luchacoviruses, they are genetically distinct from those from rodents and lagomorphs. Genetic characterization of the spike protein revealed that the meso-carnivore luchacoviruses do not have an S1/S2 cleavage motif but do have highly variable structural loops containing cleavage motifs similar to those identified in certain pathogenic CoVs. This study highlights the importance of characterizing the spike protein of CoVs in wild species for further targeted epidemiologic monitoring.

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.

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.

LncRNA-mRNA integrated profiling analysis in response to white spot syndrome virus in hepatopancreas in Penaeus japonicus

Penaeus japonicas is an important shrimp species, which is exposed to stressors including a variety of epidemic diseases. To date, little is known about the mechanisms involved in the response to white spot syndrome virus (WSSV) mediated by long non-coding RNAs (lncRNAs). A total of 6544 putative lncRNAs were identified in the hepatopancreas in P. japonicas, which provides a useful lncRNA reference resource for use in future studies. In addition, a total of 444 differentially expressed mRNAs and 457 differentially expressed lncRNAs were identified at 6, 12, and 24 h after WSSV infection in the hepatopancreas of P. japonicas. Functional enrichment analysis showed that the differentially expressed mRNAs were enriched in terms related to immune response and viral infectivity such as defense response, aminopeptidase activity, whereas the differentially expressed lncRNA partner genes were enriched in ubiquitin-dependent protein catabolic process, lipoprotein metabolic process, and antigen processing and presentation. Moreover, several lncRNAs were induced by WSSV infection, indicating these lncRNAs might participate in regulating many immune processes referring to their partner genes. Co-expression analysis of the lncRNAs and their partner genes identified some high lncRNA-mRNA correlations. These results suggest that WSSV stimulates the immune response in the hepatopancreas potentially through an important coding and non-coding gene network, thereby providing valuable information regarding non-coding responses to WSSV in Penaeus species.

Analysis of the docking property of host variants of hACE2 for SARS-CoV-2 in a large cohort

The recent novel coronavirus disease (COVID-19) outbreak, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is threatening global health. However, an understanding of the interaction of SARS-CoV-2 with human cells, including the physical docking property influenced by the host's genetic diversity, is still lacking. Here, based on germline variants in the UK Biobank covering 502,543 individuals, we revealed the molecular interactions between human angiotensin-converting enzyme 2 (hACE2), which is the representative receptor for SARS-CoV-2 entry, and COVID-19 infection. We identified six nonsense and missense variants of hACE2 from 2585 subjects in the UK Biobank covering 500000 individuals. Using our molecular dynamics simulations, three hACE2 variants from 2585 individuals we selected showed higher levels of binding free energy for docking in the range of 1.44-3.69 kcal/mol. Although there are diverse contributors to SARS-CoV-2 infections, including the mobility of individuals, we analyzed the diagnosis records of individuals with these three variants of hACE2. Our molecular dynamics simulations combined with population-based genomic data provided an atomistic understanding of the interaction between hACE2 and the spike protein of SARS-CoV-2.