VThunter: a database for single-cell screening of virus target cells in the animal kingdom.

Viral infectious diseases are a devastating and continuing threat to human and animal health. Receptor binding is the key step for viral entry into host cells. Therefore, recognizing viral receptors is fundamental for understanding the potential tissue tropism or host range of these pathogens. The rapid advancement of single-cell RNA sequencing (scRNA-seq) technology has paved the way for studying the expression of viral receptors in different tissues of animal species at single-cell resolution, resulting in huge scRNA-seq datasets. However, effectively integrating or sharing these datasets among the research community is challenging, especially for laboratory scientists. In this study, we manually curated up-to-date datasets generated in animal scRNA-seq studies, analyzed them using a unified processing pipeline, and comprehensively annotated 107 viral receptors in 142 viruses and obtained accurate expression signatures in 2 100 962 cells from 47 animal species. Thus, the VThunter database provides a user-friendly interface for the research community to explore the expression signatures of viral receptors. VThunter offers an informative and convenient resource for scientists to better understand the interactions between viral receptors and animal viruses and to assess viral pathogenesis and transmission in species. Database URL: https://db.cngb.org/VThunter/.

Genetic characteristics of human papillomavirus type 16, 18, 52 and 58 in southern China.

Persistent infections of high-risk human papillomaviruses (HPVs) are the leading cause of cervical cancers. We collected cervical exfoliated cell samples from females in Changsha city, Hunan Province and obtained 338 viral genomes of four major HPV types, including HPV 16 (n = 82), 18 (n = 35), 52 (n = 121) and 58 (n = 100). The lineage/sublineage distribution of the four HPVs confirmed previous epidemiological reports, with the predominant prevailing sublineage as A4 (50%), A1 (37%) and A3 (13%) for HPV16, A1 (83%) for HPV18, B2 (86%) for HPV52 and A1 (65%), A3 (19%) and A2 (12%) for HPV58. We also identified two potentially novel HPV18 sublineages, i.e. A6 and A7. Virus mutation analysis further revealed the presence of HPV16 and HPV58 sublineages associated with potentially high oncogenicity. These findings expanded our knowledge of the HPV genetic diversity in China, providing valuable evidence to facilitate HPV DNA screening, vaccine effectiveness evaluation and control strategy development.

Survival of viable but nonculturable Cronobacter sakazakii in macrophages contributes to infections.

Cronobacter sakazakii (C. sakazakii), a pathogen that exists in dry and low-moisture environments, such as powder infant formula (PIF), can enter a viable but nonculturable (VBNC) state under harsh conditions, which enables it to escape traditional detection methods and thus poses a potential public health risk. This study aimed at assessing the virulent nature of VBNC C. sakazakii. Our results showed that VBNC C. sakazakii induced intestinal inflammation in neonatal rats. However, the degree of inflammation was significantly lower than that of culturable bacteria due to decreasing endotoxin production, motility, adhesion, and invasion ability in the VBNC state. From the perspective of bacterial translocation, the numbers of C. sakazakii in the blood, liver, and spleen of rats treated with VBNC cells were in the same order of magnitude as those treated with its culturable counterpart and may lead to the same degree of bacteremia. According to the macrophage survival assays, the survival rate of VBNC C. sakazakii within macrophages was 4.7 times higher than that of culturable cells. Based on these findings, we hypothesize that VBNC C. sakazakii evaded the host immune defense system, penetrated the tissue barrier, and translocated to the bloodstream, liver, and spleen through macrophages. Thus, our study reveals that VBNC C. sakazakii could be a potential risk for infants' health.

The genesis and source of the H7N9 influenza viruses causing human infections in China.

A novel H7N9 influenza A virus first detected in March 2013 has since caused more than 130 human infections in China, resulting in 40 deaths. Preliminary analyses suggest that the virus is a reassortant of H7, N9 and H9N2 avian influenza viruses, and carries some amino acids associated with mammalian receptor binding, raising concerns of a new pandemic. However, neither the source populations of the H7N9 outbreak lineage nor the conditions for its genesis are fully known. Using a combination of active surveillance, screening of virus archives, and evolutionary analyses, here we show that H7 viruses probably transferred from domestic duck to chicken populations in China on at least two independent occasions. We show that the H7 viruses subsequently reassorted with enzootic H9N2 viruses to generate the H7N9 outbreak lineage, and a related previously unrecognized H7N7 lineage. The H7N9 outbreak lineage has spread over a large geographic region and is prevalent in chickens at live poultry markets, which are thought to be the immediate source of human infections. Whether the H7N9 outbreak lineage has, or will, become enzootic in China and neighbouring regions requires further investigation. The discovery here of a related H7N7 influenza virus in chickens that has the ability to infect mammals experimentally, suggests that H7 viruses may pose threats beyond the current outbreak. The continuing prevalence of H7 viruses in poultry could lead to the generation of highly pathogenic variants and further sporadic human infections, with a continued risk of the virus acquiring human-to-human transmissibility.

COVID-19 mRNA vaccine protects against SARS-CoV-2 Omicron BA.1 infection in diet-induced obese mice through boosting host innate antiviral responses.

BACKGROUND: Obesity is a worldwide epidemic and is considered a risk factor of severe manifestation of Coronavirus Disease 2019 (COVID-19). The pathogenicity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and host responses to infection, re-infection, and vaccination in individuals with obesity remain incompletely understood. METHODS: Using the diet-induced obese (DIO) mouse model, we studied SARS-CoV-2 Alpha- and Omicron BA.1-induced disease manifestations and host immune responses to infection, re-infection, and COVID-19 mRNA vaccination. FINDINGS: Unlike in lean mice, Omicron BA.1 and Alpha replicated to comparable levels in the lungs of DIO mice and resulted in similar degree of tissue damages. Importantly, both T cell and B cell mediated adaptive immune responses to SARS-CoV-2 infection or COVID-19 mRNA vaccination are impaired in DIO mice, leading to higher propensity of re-infection and lower vaccine efficacy. However, despite the absence of neutralizing antibody, vaccinated DIO mice are protected from lung damage upon Omicron challenge, accompanied with significantly more IFN-α and IFN-ß production in the lung tissue. Lung RNAseq and subsequent experiments indicated that COVID-19 mRNA vaccination in DIO mice boosted antiviral innate immune response, including the expression of IFN-α, when compared to the nonvaccinated controls. INTERPRETATION: Our findings suggested that COVID-19 mRNA vaccination enhances host innate antiviral responses in obesity which protect the DIO mice to a certain degree when adaptive immunity is suboptimal. FUNDING: A full list of funding bodies that contributed to this study can be found in the Acknowledgements section.

Single-Nucleus RNA Sequencing and Spatial Transcriptomics Reveal the Immunological Microenvironment of Cervical Squamous Cell Carcinoma.

The effective treatment of advanced cervical cancer remains challenging. Herein, single-nucleus RNA sequencing (snRNA-seq) and SpaTial enhanced resolution omics-sequencing (Stereo-seq) are used to investigate the immunological microenvironment of cervical squamous cell carcinoma (CSCC). The expression levels of most immune suppressive genes in the tumor and inflammation areas of CSCC are not significantly higher than those in the non-cancer samples, except for LGALS9 and IDO1. Stronger signals of CD56+ NK cells and immature dendritic cells are found in the hypermetabolic tumor areas, whereas more eosinophils, immature B cells, and Treg cells are found in the hypometabolic tumor areas. Moreover, a cluster of pro-tumorigenic cancer-associated myofibroblasts (myCAFs) are identified. The myCAFs may support the growth and metastasis of tumors by inhibiting lymphocyte infiltration and remodeling of the tumor extracellular matrix. Furthermore, these myCAFs are associated with poorer survival probability in patients with CSCC, predict resistance to immunotherapy, and might be present in a small fraction (< 30%) of patients with advanced cancer. Immunohistochemistry and multiplex immunofluorescence staining are conducted to validate the spatial distribution and potential function of myCAFs. Collectively, these findings enhance the understanding of the immunological microenvironment of CSCC and shed light on the treatment of advanced CSCC.

Screening of cell-virus, cell-cell, gene-gene crosstalk among animal kingdom at single cell resolution.

BACKGROUND: The exact animal origin of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains obscure and understanding its host range is vital for preventing interspecies transmission. METHODS: Herein, we applied single-cell sequencing to multiple tissues of 20 species (30 data sets) and integrated them with public resources (45 data sets covering 26 species) to expand the virus receptor distribution investigation. While the binding affinity between virus and receptor is essential for viral infectivity, understanding the receptor distribution could predict the permissive organs and tissues when infection occurs. RESULTS: Based on the transcriptomic data, the expression profiles of receptor or associated entry factors for viruses capable of causing respiratory, blood, and brain diseases were described in detail. Conserved cellular connectomes and regulomes were also identified, revealing fundamental cell-cell and gene-gene cross-talks from reptiles to humans. CONCLUSIONS: Overall, our study provides a resource of the single-cell atlas of the animal kingdom which could help to identify the potential host range and tissue tropism of viruses and reveal the host-virus co-evolution.

Mutation Profiles, Glycosylation Site Distribution and Codon Usage Bias of Human Papillomavirus Type 16.

Human papillomavirus type 16 (HPV16) is the most prevalent HPV type causing cervical cancers. Herein, using 1597 full genomes, we systemically investigated the mutation profiles, surface protein glycosylation sites and the codon usage bias (CUB) of HPV16 from different lineages and sublineages. Multiple lineage- or sublineage-conserved mutation sites were identified. Glycosylation analysis showed that HPV16 lineage D contained the highest number of different glycosylation sites from lineage A in both L1 and L2 capsid proteins, which might lead to their antigenic distances between the two lineages. CUB analysis showed that the HPV16 open reading frames (ORFs) preferred codons ending with A/T. The CUB of HPV16 ORFs was mainly affected by natural selection except for E1, E5 and L2. HPV16 only shared some of the preferred codons with humans, which might help reduce competition in translational resources. These findings increase our understanding of the heterogeneity between HPV16 lineages and sublineages, and the adaptation mechanism of HPV in human cells. In summary, this study might facilitate HPV classification and improve vaccine development and application.

Genetic signatures for lineage/sublineage classification of HPV16, 18, 52 and 58 variants.

Increasing evidences indicate that high-risk HPV variants are heterogeneous in carcinogenicity and ethnic dispersion. In this work, we identified genetic signatures for convenient determination of lineage/sublineage of HPV16, 18, 52 and 58 variants. Using publicly available genomes, we found that E2 of HPV16, L2 of HPV18, L1 and LCR of HPV52, and L2, LCR and E1 of HPV58 contain the proper genetic signature for lineage/sublineage classification. Sets of hierarchical signature nucleotide positions were further confirmed for high accuracy (>95%) by classifying HPV genomes obtained from Chinese females, which included 117 HPV16 variants, 48 HPV18 variants, 117 HPV52 variants and 89 HPV58 variants. The circulation of HPV variants posing higher cancer risk in Eastern China, such as HPV16 A4 and HPV58 A3, calls for continuous surveillance in this region. The marker genes and signature nucleotide positions may facilitate cost-effective diagnostic detections of HPV variants in clinical settings.

Population Bottlenecks and Intra-host Evolution During Human-to-Human Transmission of SARS-CoV-2.

The emergence of the novel human coronavirus, SARS-CoV-2, causes a global COVID-19 (coronavirus disease 2019) pandemic. Here, we have characterized and compared viral populations of SARS-CoV-2 among COVID-19 patients within and across households. Our work showed an active viral replication activity in the human respiratory tract and the co-existence of genetically distinct viruses within the same host. The inter-host comparison among viral populations further revealed a narrow transmission bottleneck between patients from the same households, suggesting a dominated role of stochastic dynamics in both inter-host and intra-host evolutions.

Intra-host variation and evolutionary dynamics of SARS-CoV-2 populations in COVID-19 patients.

BACKGROUND: Since early February 2021, the causative agent of COVID-19, SARS-CoV-2, has infected over 104 million people with more than 2 million deaths according to official reports. The key to understanding the biology and virus-host interactions of SARS-CoV-2 requires the knowledge of mutation and evolution of this virus at both inter- and intra-host levels. However, despite quite a few polymorphic sites identified among SARS-CoV-2 populations, intra-host variant spectra and their evolutionary dynamics remain mostly unknown. METHODS: Using high-throughput sequencing of metatranscriptomic and hybrid captured libraries, we characterized consensus genomes and intra-host single nucleotide variations (iSNVs) of serial samples collected from eight patients with COVID-19. The distribution of iSNVs along the SARS-CoV-2 genome was analyzed and co-occurring iSNVs among COVID-19 patients were identified. We also compared the evolutionary dynamics of SARS-CoV-2 population in the respiratory tract (RT) and gastrointestinal tract (GIT). RESULTS: The 32 consensus genomes revealed the co-existence of different genotypes within the same patient. We further identified 40 intra-host single nucleotide variants (iSNVs). Most (30/40) iSNVs presented in a single patient, while ten iSNVs were found in at least two patients or identical to consensus variants. Comparing allele frequencies of the iSNVs revealed a clear genetic differentiation between intra-host populations from the respiratory tract (RT) and gastrointestinal tract (GIT), mostly driven by bottleneck events during intra-host migrations. Compared to RT populations, the GIT populations showed a better maintenance and rapid development of viral genetic diversity following the suspected intra-host bottlenecks. CONCLUSIONS: Our findings here illustrate the intra-host bottlenecks and evolutionary dynamics of SARS-CoV-2 in different anatomic sites and may provide new insights to understand the virus-host interactions of coronaviruses and other RNA viruses.

Single cell atlas for 11 non-model mammals, reptiles and birds.

The availability of viral entry factors is a prerequisite for the cross-species transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Large-scale single-cell screening of animal cells could reveal the expression patterns of viral entry genes in different hosts. However, such exploration for SARS-CoV-2 remains limited. Here, we perform single-nucleus RNA sequencing for 11 non-model species, including pets (cat, dog, hamster, and lizard), livestock (goat and rabbit), poultry (duck and pigeon), and wildlife (pangolin, tiger, and deer), and investigated the co-expression of ACE2 and TMPRSS2. Furthermore, cross-species analysis of the lung cell atlas of the studied mammals, reptiles, and birds reveals core developmental programs, critical connectomes, and conserved regulatory circuits among these evolutionarily distant species. Overall, our work provides a compendium of gene expression profiles for non-model animals, which could be employed to identify potential SARS-CoV-2 target cells and putative zoonotic reservoirs.

A Path toward SARS-CoV-2 Attenuation: Metabolic Pressure on CTP Synthesis Rules the Virus Evolution.

In the context of the COVID-19 pandemic, we describe here the singular metabolic background that constrains enveloped RNA viruses to evolve toward likely attenuation in the long term, possibly after a step of increased pathogenicity. Cytidine triphosphate (CTP) is at the crossroad of the processes allowing SARS-CoV-2 to multiply, because CTP is in demand for four essential metabolic steps. It is a building block of the virus genome, it is required for synthesis of the cytosine-based liponucleotide precursors of the viral envelope, it is a critical building block of the host transfer RNAs synthesis and it is required for synthesis of dolichol-phosphate, a precursor of viral protein glycosylation. The CCA 3'-end of all the transfer RNAs required to translate the RNA genome and further transcripts into the proteins used to build active virus copies is not coded in the human genome. It must be synthesized de novo from CTP and ATP. Furthermore, intermediary metabolism is built on compulsory steps of synthesis and salvage of cytosine-based metabolites via uridine triphosphate that keep limiting CTP availability. As a consequence, accidental replication errors tend to replace cytosine by uracil in the genome, unless recombination events allow the sequence to return to its ancestral sequences. We document some of the consequences of this situation in the function of viral proteins. This unique metabolic setup allowed us to highlight and provide a raison d'être to viperin, an enzyme of innate antiviral immunity, which synthesizes 3'-deoxy-3',4'-didehydro-CTP as an extremely efficient antiviral nucleotide.

A compromised specific humoral immune response against the SARS-CoV-2 receptor-binding domain is related to viral persistence and periodic shedding in the gastrointestinal tract.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been redetected after discharge in some coronavirus disease 2019 (COVID-19) patients. The reason for the recurrent positivity of the test and the potential public health concern due to this occurrence are still unknown. Here, we analyzed the viral data and clinical manifestations of 289 domestic Chinese COVID-19 patients and found that 21 individuals (7.3%) were readmitted for hospitalization after detection of SARS-CoV-2 after discharge. First, we experimentally confirmed that the virus was involved in the initial infection and was not a secondary infection. In positive retests, the virus was usually found in anal samples (15 of 21, 71.4%). Through analysis of the intracellular viral subgenomic messenger RNA (sgmRNA), we verified that positive retest patients had active viral replication in their gastrointestinal tracts (3 of 16 patients, 18.7%) but not in their respiratory tracts. Then, we found that viral persistence was not associated with high viral titers, delayed viral clearance, old age, or more severe clinical symptoms during the first hospitalization. In contrast, viral rebound was associated with significantly lower levels of and slower generation of viral receptor-binding domain (RBD)-specific IgA and IgG antibodies. Our study demonstrated that the positive retest patients failed to create a robust protective humoral immune response, which might result in SARS-CoV-2 persistence in the gastrointestinal tract and possibly in active viral shedding. Further exploration of the mechanism underlying the rebound in SARS-CoV-2 in this population will be crucial for preventing virus spread and developing effective vaccines.

Multiple approaches for massively parallel sequencing of SARS-CoV-2 genomes directly from clinical samples.

BACKGROUND: COVID-19 (coronavirus disease 2019) has caused a major epidemic worldwide; however, much is yet to be known about the epidemiology and evolution of the virus partly due to the scarcity of full-length SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) genomes reported. One reason is that the challenges underneath sequencing SARS-CoV-2 directly from clinical samples have not been completely tackled, i.e., sequencing samples with low viral load often results in insufficient viral reads for analyses. METHODS: We applied a novel multiplex PCR amplicon (amplicon)-based and hybrid capture (capture)-based sequencing, as well as ultra-high-throughput metatranscriptomic (meta) sequencing in retrieving complete genomes, inter-individual and intra-individual variations of SARS-CoV-2 from serials dilutions of a cultured isolate, and eight clinical samples covering a range of sample types and viral loads. We also examined and compared the sensitivity, accuracy, and other characteristics of these approaches in a comprehensive manner. RESULTS: We demonstrated that both amplicon and capture methods efficiently enriched SARS-CoV-2 content from clinical samples, while the enrichment efficiency of amplicon outran that of capture in more challenging samples. We found that capture was not as accurate as meta and amplicon in identifying between-sample variations, whereas amplicon method was not as accurate as the other two in investigating within-sample variations, suggesting amplicon sequencing was not suitable for studying virus-host interactions and viral transmission that heavily rely on intra-host dynamics. We illustrated that meta uncovered rich genetic information in the clinical samples besides SARS-CoV-2, providing references for clinical diagnostics and therapeutics. Taken all factors above and cost-effectiveness into consideration, we proposed guidance for how to choose sequencing strategy for SARS-CoV-2 under different situations. CONCLUSIONS: This is, to the best of our knowledge, the first work systematically investigating inter- and intra-individual variations of SARS-CoV-2 using amplicon- and capture-based whole-genome sequencing, as well as the first comparative study among multiple approaches. Our work offers practical solutions for genome sequencing and analyses of SARS-CoV-2 and other emerging viruses.