A novel astrovirus identified in wild rhesus monkey feces in China.

The discovery and analysis of pathogens carried by non-human primates are important for understanding zoonotic infections in humans. We identified a highly divergent astrovirus (AstV) from fecal matter from a rhesus monkey in China, which has been tentatively named "monkey-feces-associated AstV" (MkAstV). The full-length genome of MkAstV was determined to be 7377 nt in length. It exhibits the standard genomic AstV organization of three open reading frames (ORFs) and is most closely related to duck AstV (28%, 49%, and 35% amino acid sequence identity in ORF1a, ORF1b, and ORF2, respectively). Coincidentally, while this report was being prepared, an astrovirus sequence from Hainan black-spectacled toad became available in the GenBank database, showing 95%, 94% and 92% aa sequence identity in ORF1a, ORF1b and ORF2, respectively, to the corresponding ORFs of MkAstV. Phylogenetic analysis of ORF1a, ORF1b, and ORF2 indicated that MkAstV and the amphibian-related astroviruses formed an independent cluster in the genus Avastrovirus. The host of MkAstV remains unknown. Epidemiological and serological studies of this novel virus should be undertaken in primates, including humans.

Experimental infection of Marmota monax with a novel hepatitis A virus.

To establish an animal model for the newly identified Marmota Himalayana hepatovirus, MHHAV, so as to develop a better understanding of the infection of hepatitis A viruses. Five experimental woodchucks (Marmota monax) were inoculated intravenously with the purified MHHAV from wild woodchuck feces. One animal injected with PBS was defined as a control. Feces and blood were routinely collected. After the animals were subjected to necropsy, different tissues were collected. The presence of viral RNA and negative sense viral RNA was analyzed in all the samples and histopathological and in situ hybridization analysis was performed for the tissues. MHHAV infection caused fever but no severe symptoms or death. Virus was shed in feces beginning at 2 dpi, and MHHAV RNA persisted in feces for ~2 months, with a biphasic increase, and in blood for ~30 days. Viral RNA was detected in all the tissues, with high levels in the liver and spleen. Negative-strand viral RNA was detected only in the liver. Furthermore, the animals showed histological signs of hepatitis at 45 dpi. MHHAV can infect M. monax and is associated with hepatic disease. Therefore, this animal can be used as a model of HAV pathogenesis and to evaluate antiviral and anticancer therapeutics.

Diverse novel astroviruses identified in wild Himalayan marmots.

With advances in viral surveillance and next-generation sequencing, highly diverse novel astroviruses (AstVs) and different animal hosts had been discovered in recent years. However, the existence of AstVs in marmots had yet to be shown. Here, we identified two highly divergent strains of AstVs (tentatively named Qinghai Himalayanmarmot AstVs, HHMAstV1 and HHMAstV2), by viral metagenomic analysis in liver tissues isolated from wild Marmota himalayana in China. Overall, 12 of 99 (12.1 %) M. himalayana faecal samples were positive for the presence of genetically diverse AstVs, while only HHMAstV1 and HHMAstV2 were identified in 300 liver samples. The complete genomic sequences of HHMAstV1 and HHMAstV2 were 6681 and 6610 nt in length, respectively, with the typical genomic organization of AstVs. Analysis of the complete ORF 2 sequence showed that these novel AstVs are most closely related to the rabbit AstV, mamastrovirus 23 (with 31.0 and 48.0 % shared amino acid identity, respectively). Phylogenetic analysis of the amino acid sequences of ORF1a, ORF1b and ORF2 indicated that HHMAstV1 and HHMAstV2 form two distinct clusters among the mamastroviruses, and may share a common ancestor with the rabbit-specific mamastrovirus 23. These results suggest that HHMAstV1 and HHMAstV2 are two novel species of the genus Mamastrovirus in the Astroviridae. The remarkable diversity of these novel AstVs will contribute to a greater understanding of the evolution and ecology of AstVs, although additional studies will be needed to understand the clinical significance of these novel AstVs in marmots, as well as in humans.

Complete genome sequence of a novel human papillomavirus identified by metagenomic analysis from a child with diarrhea in China.

A highly divergent human papillomavirus, named HPV-CH2, was identified in fecal samples from children with diarrhea in China by 454 high-throughput sequencing. Here, we report the complete genome sequence and genetic organization of the virus. The full-length nucleotide sequence of HPV-CH2 shares the highest sequence similarity with HPV-156, with 72 % nucleotide sequence identity. The L1 gene of the HPV-CH2 shared <70 % nucleotide identity with previously reported HPVs, suggesting HPV-CH2 as a new type of papillomavirus. Phylogenetic analysis revealed that the HPV-CH2 belongs to the genus Gammapapillomavirus. No HPV-CH2 was detected by PCR in samples from children with both gastroenteritis and respiratory infection.

Detection of human bocavirus in children with acute respiratory tract infections in Lanzhou and Nanjing, China.

OBJECTIVE: The aim of this study was to explore the prevalent characteristics of HBoV1 and its co-infection. METHODS: PCR was used to detect HBoV1-DNA (HBoV1) and other viruses. A multivariate logistic regression model was used to explore possibility of co-detected for related viruses. RESULTS: The positivity rates in Nanjing and Lanzhou were 9.38% (74/789) and 11.62% (161/1386), respectively (P>0.05). The HBoV1 positive group was younger than negative group (P<0.05). Seasonal differences were noted, with a higher frequency of infection in December and July. HBoV1-positive children [72.34% (169/235)] were co-infected with other respiratory viruses. Multifactorial analysis showed no correlations between HBoV1 and the clinical classification, region, gender, age, or treatment as an outpatient or in a hospital. Correlations were identified between HBoV1 infections with ADV (OR=1.53, 95% CI 1.03-2.28), RSV (OR=0.71, 95% CI 0.52-0.98), and IFVA (OR=1.77, 95% CI 1.00-3.13). CONCLUSION: Presence of HBoV1 in nasopharyngeal aspirates did not correlate with region or gender, although the prevalence of HBoV1 was higher in younger children. There were no correlations between HBoV1 and other variables, except for the season and ADV, RSV, or IFVA infections.

Coexistence of adeno-associated virus 2 with adenovirus 18 or herpesvirus may be associated with severe lingual papillomatosis in an immunocompromised individual.

Next-generation sequencing (NGS) has significantly improved the accuracy and efficiency of pathogen diagnosis for a wide range of diseases. In this study, viral metagenomics analysis was conducted on fecal and tissue samples from a 13-year-old recipient of hematopoietic stem cell transplantation (HSCT) afflicted with severe lingual papillomatosis. The analysis revealed a high abundance of adeno-associated virus 2 (AAV2), alongside potential helper viruses, herpesvirus type 1 (HSV-1), and the uncommon adenovirus serotype 18 (AdV18). Although a direct causal relationship was not definitively established, the concurrence of these viruses indicated a plausible link to the development of severe lingual papillomatosis in immunocompromised individuals. Notably, the study generated a complete genome sequence of AdV18, offering insights into adenovirus genetic variability, origin, and pathogenicity. Noteworthy findings include three amino acid substitutions in the polymerase and one in the hexon, distinguishing them from previously published strains of AdV18. Phylogenetic analysis unveiled a close relationship between both the polymerase and hexon regions of AdV18 in our study and previously reported AdV18 sequences. This study underscores the pivotal role of comprehensive viral scrutiny in elucidating infections among HSCT patients with lingual papillomatosis.

Detection of novel viruses in porcine fecal samples from China.

BACKGROUND: Pigs are well known source of human infectious disease. To better understand the spectrum of viruses present in pigs, we utilized the 454 Life Sciences GS-FLX high-throughput sequencing platform to sequence stool samples from healthy pigs. FINDINGS: Total nucleic acid was extracted from stool samples of healthy piglets and randomly amplified. The amplified materials were pooled and processed using a high-throughput pyrosequencing technique. The raw sequences were deconvoluted on the basis of the barcode and then processed through a standardized bioinformatics pipeline. The unique reads (348, 70 and 13) had limited similarity to known astroviruses, bocaviruses and parechoviruses. Specific primers were synthesized to assess the prevalence of the viruses in healthy piglets. Our results indicate extremely high rates of positivity. CONCLUSIONS: Several novel astroviruses, bocaviruses and Ljungan-like viruses were identified in stool samples from healthy pigs. The rates of isolation for the new viruses were high. The high detection rate, diverse sequences and categories indicate that pigs are well-established reservoirs for and likely sources of different enteric viruses.

Tracking of Mutational Signature of SARS-CoV-2 Omicron on Distinct Continents and Little Difference was Found.

The Omicron variant is currently ravaging the world, raising serious concern globally. Monitoring genomic variations and determining their influence on biological features are critical for tracing its ongoing transmission and facilitating effective measures. Based on large-scale sequences from different continents, this study found that: (i) The genetic diversity of Omicron is much lower than that of the Delta variant. Still, eight deletions (Del 1-8) and 1 insertion, as well as 130 SNPs, were detected on the Omicron genomes, with two deletions (Del 3 and 4) and 38 SNPs commonly detected on all continents and exhibiting high-occurring frequencies. (ii) Four groups of tightly linked SNPs (linkage I-IV) were detected, among which linkage I, containing 38 SNPs, with 6 located in the RBD, increased its occurring frequency remarkably over time. (iii) The third codons of the Omicron shouldered the most mutation pressures, while the second codons presented the least flexibility. (iv) Four major mutants with amino acid substitutions in the RBD were detected, and further structural analysis suggested that the substitutions did not alter the viral receptor binding ability greatly. It was inferred that though the Omicron genome harbored great changes in antigenicity and remarkable ability to evade immunity, it was immune-pressure selected. This study tracked mutational signatures of Omicron variant and the potential biological significance of the SNPs, and the linkages await further functional verification.

Comparative analysis of human respiratory syncytial virus evolutionary patterns during the COVID-19 pandemic and pre-pandemic periods.

The COVID-19 pandemic has resulted in the implementation of strict mitigation measures that have impacted the transmission dynamics of human respiratory syncytial virus (HRSV). The measures also have the potential to influence the evolutionary patterns of the virus. In this study, we conducted a comprehensive analysis comparing genomic variations and evolving characteristics of its neutralizing antigens, specifically F and G proteins, before and during the COVID-19 pandemic. Our findings showed that both HRSV A and B exhibited an overall chronological evolutionary pattern. For the sequences obtained during the pandemic period (2019-2022), we observed that the HRSV A distributed in A23 genotype, but formed into three subclusters; whereas the HRSV B sequences were relatively concentrated within genotype B6. Additionally, multiple positively selected sites were detected on F and G proteins but none were located at neutralizing antigenic sites of the F protein. Notably, amino acids within antigenic site III, IV, and V of F protein remained strictly conserved, while some substitutions occurred over time on antigenic site Ø, I, II and VIII; substitution S389P on antigenic site I of HRSV B occurred during the pandemic period with nearly 50% frequency. However, further analysis revealed no substitutions have altered the structural conformations of the antigenic sites, the vial antigenicity has not been changed. We inferred that the intensive public health interventions during the COVID-19 pandemic did not affect the evolutionary mode of HRSV.

Simultaneous detection of seven enteric viruses associated with acute gastroenteritis by a multiplexed Luminex-based assay.

Rapid and broad diagnostic methods are needed for the identification of viral agents of gastroenteritis. In this study, we used Luminex xMAP technology to develop a multiplexed assay for the simultaneous identification of major enteric viral pathogens, including rotavirus A (RVA), noroviruses (NoVs) (including genogroups GI and GII), sapoviruses (SaV), human astrovirus (HAstV), enteric adenoviruses (EAds), and human bocavirus 2 (HBoV2). The analytical sensitivity allowed detection of 10(3) (EAds, HBoV2, and RVA) and 10(4) (NoV GI and GII, SaV, and HAstV) copies per reaction mixture. Compared to conventional PCR, the Luminex-based assay yielded greater than 75% sensitivity and 97% specificity for each virus, and the kappa correlation for detection of all viruses ranged from 0.75 to 1.00. In conclusion, this multiplexed Luminex-based assay provides a potentially rapid, high-throughput, and maneuverable diagnostic tool for major viral pathogens associated with gastroenteritis.

Genome characterization of a novel porcine bocavirus.

Using a high-throughput DNA sequencing method, one DNA sequence (contig01006), suspected to belong to a novel porcine bocavirus (PBoV), was found with a high rate of detection (19.6 %) in fecal samples from healthy piglets. Moreover, a novel PBoV (tentatively named PBoV3C) with a nearly complete genome sequence (5235 bp) was identified. PBoV3C exhibits typical genome characteristics of bocaviruses and shows the highest genomic sequence identity (78 % to 81 %) to PBoV3A/B (PBoV3/4-UK) and PBoV3D/E (PBoV3/4-HK), respectively. Phylogenetic and recombination analysis indicated high diversity, prevalence and complexity among the PBoVs. The phospholipase A2 (PLA2) site of VP1 and the secondary structure of VP2 of PBoV3C were also analyzed. Additionally, we propose a uniform method of PBoV nomenclature based on the VP1 gene.

What is the Potential Cause for the Predominance of GII.2[P16] Norovirus in Acute Gastroenteritis Outbreaks in China?

GII.2[P16] noroviruses (NoV) reemerged and rapidly became the main epidemic strain in acute gastroenteritis (AGE) outbreaks in Asian countries since 2016. The current GII.2 [P16] NoV showed the same antigenicity to the ones before 2016, but several unique amino acid substitutions existed in the RNA dependent RNA polymerase (RdRp) and other non-structural proteins, and the viral load of the current GII.2[P16] NoV was higher than those of other genotypes, it was estimated that the viral replication ability may have improved. However, other genotypes, such as GII.1 and GII.3, also had recombination with the novel RdRp, were not prevalent in AGE-outbreaks; thus, it was inferred that the capsid proteins also played an important role in the enhanced replication process. The viral infection could also be affected by other factors, such as the population genetic background, the climate and environment, and people's lifestyles. Continued surveillance on genetic diversity and evolutionary pattern for the GII.2[P16] NoV is necessary.

Evolutionary analyses of emerging GII.2[P16] and GII.4 Sydney [P16] noroviruses.

GII.2[P16] and GII.4 Sydney [P16] are currently the two predominant norovirus genotypes. This study sought to clarify their evolutionary patterns by analyzing the major capsid VP1 and RNA-dependent RNA polymerase (RdRp) genes. Sequence diversities were analyzed at both nucleotide and amino acid levels. Selective pressures were evaluated with the Hyphy package in different models. Phylogenetic trees were constructed by the maximum likelihood method from full VP1 sequences, and evolutionary rates were estimated by the Bayesian Markov Chain Monte Carlo approach. The results showed that (1) several groups of tightly linked mutations between the RdRp and VP1 genes were detected in the GII.2[P16] and GII.4[P16] noroviruses, and most of these mutations were synonymous, which may lead to a better viral fitness to the host; (2) although the pattern of having new GII.4 variants every 2-4 years has been broken, both the pre- and the post-2015 Sydney VP1 had comparable evolutionary rates to previously epidemic GII.4 variants, and half of the major antigenic sites on GII.4 Sydney had residue substitutions and several caused obvious changes in the carbohydrate-binding surface that may potentially alter the property of the virus; and (3) GII.4 Sydney variants during 2018-21 showed geographical specificity in East Asia, South Asia, and North America; the antigenic sites of GII.2 are strictly conserved, but the GII.2 VP1 chronologically evolved into nine different sublineages over time, with sublineage IX being the most prevalent one since 2018. This study suggested that both VP1 and RdRp of the GII.2[P16] and GII.4 Sydney [P16] noroviruses exhibited different evolutionary directions. GII.4[P16] is likely to generate potential novel epidemic variants by accumulating mutations in the P2 domain, similar to previously epidemic GII.4 variants, while GII.2[P16] has conserved predicted antigenicity and may evolve by changing the properties of nonstructural proteins, such as polymerase replicational fidelity and efficiency. This study expands the understanding of the evolutionary dynamics of GII.2[P16] and GII.4[P16] noroviruses and may predict the emergence of new variants.

The protective immunity induced by intranasally inoculated serotype 63 chimpanzee adenovirus vector expressing human respiratory syncytial virus prefusion fusion glycoprotein in BALB/c mice.

Human respiratory syncytial virus (RSV) is a ubiquitous pediatric pathogen causing serious lower respiratory tract disease worldwide. No licensed vaccine is currently available. In this work, the coding gene for mDS-Dav1, the full-length and prefusion conformation RSV fusion glycoprotein (F), was designed by introducing the stabilized prefusion F (preF) mutations from DS-Cav1 into the encoding gene of wild-type RSV (wtRSV) F protein. The recombinant adenovirus encoding mDS-Cav1, rChAd63-mDS-Cav1, was constructed based on serotype 63 chimpanzee adenovirus vector and characterized in vitro. After immunizing mice via intranasal route, the rChAd63-mDS-Cav1 induced enhanced neutralizing antibody and F-specific CD8+ T cell responses as well as good immune protection against RSV challenge with the absence of enhanced RSV disease (ERD) in BALB/c mice. The results indicate that rChAd63-mDS-Cav1 is a promising mucosal vaccine candidate against RSV infection and warrants further development.

Analysis and characterization of the complete genome of a member of a new species of kobuvirus associated with swine.

A virus belonging to a new species in the genus Kobuvirus, family Picornaviridae, was first isolated in 2008 from apparently healthy pigs in Hungary and China. We report the complete genome sequence and the genetic organization of the novel porcine kobuvirus strain Y-1-CHI, which was identified in China. The RNA genome of strain Y-1-CHI contains 8210 nucleotides (nt) and has an organization similar to that of other picornaviruses. The full-length nucleotide sequence of Y-1-CHI was 88.62%, 58.66%, and 48.86% identical to those of S-1-HUN, U-1, and Aichi virus, respectively. No positive results were found in 454 stool samples from children with acute gastroenteritis. Dendrograms indicated that Y-1-CHI and S-1-HUN are most closely related to each other and belong to the same species. Our results suggest that members of this novel species have the typical genome characteristics of members of the genus Kobuvirus and may be distributed globally in swine.

Genetic diversity and molecular evolution of human respiratory syncytial virus A and B.

Human respiratory syncytial viruses (RSVs) are classified into two major groups (A and B) based on antigenic differences in the G glycoprotein. To investigate circulating characteristics and phylodynamic history of RSV, we analyzed the genetic variability and evolutionary pattern of RSVs from 1977 to 2019 in this study. The results revealed that there was no recombination event of intergroup. Single nucleotide polymorphisms (SNPs) were observed through the genome with the highest occurrence rate in the G gene. Five and six sites in G protein of RSV-A and RSV-B, respectively, were further identified with a strong positive selection. The mean evolutionary rates for RSV-A and -B were estimated to be 1.48 × 10-3 and 1.92 × 10-3 nucleotide substitutions/site/year, respectively. The Bayesian skyline plot showed a constant population size of RSV-A and a sharp expansion of population size of RSV-B since 2005, and an obvious decrease 5 years later, then became stable again. The total population size of RSVs showed a similar tendency to that of RSV-B. Time-scaled phylogeny suggested a temporal specificity of the RSV-genotypes. Monitoring nucleotide changes and analyzing evolution pattern for RSVs could give valuable insights for vaccine and therapy strategies against RSV infection.

Biological Significance of the Genomic Variation and Structural Dynamics of SARS-CoV-2 B.1.617.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants have been emerging and circulating globally since the start of the COVID-19 pandemic, of which B.1.617 lineage that was first reported in India at the end of 2020, soon became predominant. Tracing genomic variations and understanding their impact on the viral properties are the foundations for the vaccine and drug development and for the mitigation measures to be taken or lifted. In this study, 1,051 near-complete genomes and 1,559 spike (S) sequences belonging to the B.1.617 were analyzed. A genome-wide spread of single nucleotide polymorphisms (SNPs) was identified. Of the high frequency mutations identified, 61% (11/18) involved structural proteins, despite two third of the viral genome encoding nonstructural proteins. There were 22 positive selection sites, mostly distributed across the S protein, of which 16 were led by non-C to U transition and should be of a special attention. Haplotype network revealed that a large number of daughter haplotypes were continually derived throughout the pandemic, of which H177, H181 H219 and H286 from the ancestor haplotype H176 of B.1.617.2 were widely prevalent. Besides the well known substitutions of L452R, P681R and deletions of E156 and F157, as well as the potential biological significance, structural analysis in this study still indicated that new amino acid changes in B.1.617, such as E484Q and N501Y, had reshaped the viral bonding network, and increasingly sequenced N501Y mutant with a potential enhanced binding ability was detected in many other countries in the follow-up monitoring. Although we can't conclude the properties of all the mutants including N501Y thoroughly, it merits focusing on their spread epidemically and biologically.

Evaluation of the Safety and Immune Efficacy of Recombinant Human Respiratory Syncytial Virus Strain Long Live Attenuated Vaccine Candidates.

Human respiratory syncytial virus (RSV) infection is the leading cause of lower respiratory tract illness (LRTI), and no vaccine against LRTI has proven to be safe and effective in infants. Our study assessed attenuated recombinant RSVs as vaccine candidates to prevent RSV infection in mice. The constructed recombinant plasmids harbored (5' to 3') a T7 promoter, hammerhead ribozyme, RSV Long strain antigenomic cDNA with cold-passaged (cp) mutations or cp combined with temperature-sensitive attenuated mutations from the A2 strain (A2cpts) or further combined with SH gene deletion (A2cptsΔSH), HDV ribozyme (δ), and a T7 terminator. These vectors were subsequently co-transfected with four helper plasmids encoding N, P, L, and M2-1 viral proteins into BHK/T7-9 cells, and the recovered viruses were then passaged in Vero cells. The rescued recombinant RSVs (rRSVs) were named rRSV-Long/A2cp, rRSV-Long/A2cpts, and rRSV-Long/A2cptsΔSH, respectively, and stably passaged in vitro, without reversion to wild type (wt) at sites containing introduced mutations or deletion. Although rRSV-Long/A2cpts and rRSV-Long/A2cptsΔSH displayed  temperature-sensitive (ts) phenotype in vitro and in vivo, all rRSVs were significantly attenuated in vivo. Furthermore, BALB/c mice immunized with rRSVs produced Th1-biased immune response, resisted wtRSV infection, and were free from enhanced respiratory disease. We showed that the combination of ΔSH with attenuation (att) mutations of cpts contributed to improving att phenotype, efficacy, and gene stability of rRSV. By successfully introducing att mutations and SH gene deletion into the RSV Long parent and producing three rRSV strains, we have laid an important foundation for the development of RSV live attenuated vaccines.

Complete genome of a novel recombinant human astrovirus and its quasispecies in patients following hematopoietic stem cell transplantation.

Human astroviruses (HAstVs) were first identified in 1975 and can be classified into three clades: classic HAstVs (HAstV 1-8), MLB (MLB1-3) and VA (VA1-5), with MLB and VA were newly identified. Recombination and a high mutation rate make HAstV as one of the rapidly evolving infectious agents. This study reported a novel identified recombinant human astrovirus (Y/1-CHN) and its long existence in two immunocompromised patients with diarrhea following allogeneic hematopoietic stem cell transplantation (allo-HSCT). The identified Yu/1-CHN genome contains 6801 base pairs encoding three open reading frames, with ORF1a best hit to the HAstV1 (Pune strain, 97 % nucleotide identity), while ORF1b and ORF2 best hit to HAstV-5 (DL30 strain, 99 % nucleotide identity). Possible recombination breakpoint was predicted to be located in the boundary of ORF1a and ORF1b. Different quasispecies were found in the host, and the dN/dS ratios of the S and P domains were determined to be 1.189 and 1.444, respectively, suggesting a positive selection existed. Fecal samples collected in different clinical phases from the two patients were all positive for Yu/1-CHN, suggesting a long existence of the virus in the host. It was indicated that immunocompromised patients may a reservoir for astrovirus, their excreta should be monitored even after discharge from hospital.

Intra-Host Evolution of Norovirus GII.4 in a Chronic Infected Patient With Hematopoietic Stem Cell Transplantation.

Human noroviruses (NVs) are the leading cause of acute gastroenteritis outbreaks worldwide. The majority of outbreaks are caused by genogroup II.4 (GII.4), with new variants emerging every 2 to 4 years. Immunocompromised patients are hypothesized to be important reservoirs where new NV variants emerge. Here, we examined intra-host NV variants and assessed immune-driven NV evolution in chronically infected immunocompromised hosts. Three NV GII.4-positive samples were collected from the same patient in different clinical phases following allogeneic hematopoietic stem cell transplantation, and had viral RNA concentrations of 2.46 × 106, 1.47 × 106, and 2.26 × 106 genome copies/mL. The non-synonymous (dN) and synonymous (dS) substitution ratio of the sequences in the partial P domain were >1, indicating strong positive selection in the patient. Both the number and the frequency of the single nucleotide variants increased over time in the patient. Also, the majority of capsid amino acid changes were located at blocking epitopes and histo-blood group antigen (HBGA)-binding sites, and 11 positive selection sites were found in the capsid region, of which 8 sites were presented in blocking epitopes or HBGA-binding sites. Homodimeric P-domain capsid models also suggested a structural change in the epitopes and HBGA-binding sites. The results suggested that novel variants of NV GII.4 with HBGA and antigenic site changes were produced in the immunocompromised patient. Further functional and epidemiological studies are needed to determine whether the new variants are a risk to public health.