Genomic and pathogenicity analysis of two novel highly pathogenic recombinant NADC30-like PRRSV strains in China, in 2023.

Porcine reproductive and respiratory syndrome viruses (PRRSVs) exhibit high mutability and recombination, posing challenges to their immunization and control. This study isolated two new PRRSV strains, GD-7 and GX-3, from samples collected in Guangdong and Guangxi in 2023. Whole-genome sequencing, along with phylogenetic and recombination analyses, confirmed that GD-7 and GX-3 are natural novel recombinant strains of NADC30 PRRSV. Moreover, we established a pathogenicity model for piglets and sows based on the two isolates. The results of piglet pathogenicity revealed that both GD-7 and GX-3 caused clinical symptoms such as fever, loss of appetite, depression, and slow weight gain. Moreover, we observed that the mortality rate of GD-7-inoculated group piglets was 33.3%, which was similar to that of piglets infected with other highly pathogenic PRRSV strains and exceeded the mortality rate of most NADC30-like PRRSV. In pregnant sow models, the survival rate of sows in the GD-7 group was 75%, in contrast to the GX-3 group, where no sow mortality was observed, and both strains resulted in abortion, mummified fetuses, and stillbirths. These results highlight the elevated pathogenicity of these recombinant strains in sows, with GD-7 mainly causing sows to abort, and GX-3 mainly causing sows to give birth to mummified fetuses. This study introduces two distinct clinical recombinant PRRSV strains that differ from the prevalent strains in China. This research furthers our understanding of the epidemiology of PRRSV and underscores the significance of ongoing monitoring and research in the face of evolving virus strains. Moreover, these discoveries act as early warnings, underscoring the necessity for active control and immunization against PRRSV.IMPORTANCESince the discovery of NADC30-like PRRSV in China in 2013, it has gradually become the dominant strain of PRRSV in China. NADC30-like PRRSV exhibits high recombination characteristics, constantly recombining with different strains, leading to the emergence of numerous novel strains. Of particular importance is the observation that NADC30-like PRRSV with different recombination patterns exhibits varying pathogenicity, which has a significant impact on the pig farming industry. This emphasizes the necessity of monitoring and responding to evolving PRRSV strains to develop effective immunization and control strategies. In this paper, we conducted pathogenicity studies on the isolated NADC30-like PRRSV and analyzed the differences in the genomes and pathogenicity of the different strains by recording clinical symptoms, temperature changes, detoxification tests, and changes in viremia and histopathology in infected pigs. This was done to provide a theoretical basis for the epidemiological situation and epidemic prevention and control of PRRSV.

Genetic heterogeneity and potential recombination across hosts of Gyrovirus galga1 in central and eastern China during 2021 to 2024.

Gyrovirus galga1 (GyVg1), formerly known as AGV2, was initially identified in chickens in southern Brazil. The prevalence of GyVg1 from 2021 to 2024 in 28 out of the 63 poultry farms located in Jiangsu, Anhui, Henan, Hunan, Shandong, and Hubei provinces in eastern and central China was detected via PCR. The complete genomes of the 28 strains were sequenced and exhibited a full length of 2,376 bp. Similarity analysis of these strains did not suggest definite correlation with evolutionary branching and geographical distribution. Compared with the reference GyVg1 strains, HN2202 shared the highest similarity of 99.71% with HLJ1511 (chicken-originated) from northeastern China in 2015 to 2016. Recombination analysis revealed that AH2102 was a potential recombinant of peafowl-originated HN2019-PF1 and chicken-originated HLJ1506-2, whereas HN2304 was a recombinant of peafowl-originated HN2019-PF1 and the Hungarian ferret strain G13. Mutation site analysis of the capsid protein revealed that highly mutated regions occurred between sites 288 to 316 and 383 to 419. These results indicate that GyVg1 may have undergone an interspecies transmission, which involved complex mutations and recombination. This study may provide a reference for subsequent investigations targeting the molecular epidemiology and viral evolution of GyVg1.

Phylogenetic and Genetic Variation Analysis of Porcine Epidemic Diarrhea Virus in East Central China during 2020-2023.

Porcine epidemic diarrhea virus (PEDV) is a major causative pathogen of a highly contagious, acute enteric viral disease. This study evaluated the emergence of nine variants in Jiangsu and Anhui provinces of China from 2020 to 2023. S gene-based phylogenetic analysis indicated that three variants belong to the G1c subgroup, while the other six strains are clustered within the G2c subgroup. Recombination analyses supported that three variants of the G1c subgroup were likely derived from recombination of parental variants FR0012014 and a donor variant AJ1102. In addition, there are novel mutations on amino acid 141-148 and these likely resulted in changes in antigenicity in the three variants. These results illustrated that the study provides novel insights into the epidemiology, evolution, and transmission of PEDV in China.

First report of goose circovirus identified in ducks from China.

Goose circovirus (GoCV) is a common pathogen that causes immunosuppression and promotes secondary infections with other infectious agents in geese worldwide. In the present study, we identified GoCV in 2 out of 93 duck flocks from China and successfully sequenced the complete genomes of 2 strains (AH22du and HN20du). The whole genome of the two strains shared a high identity of 90.5 to 98.63% with China GoCV reference, and low identity of 58.98% with DuCV reference, respectively. Phylogenetic tree constructed on the two and other genome sequences of GoCV revealed three main branches. Both strains sequenced in this study were distributed on different sub-branches with most other Chinese GoCV strains, and AH22du clustered into an independent sub-branch within the cluster. Recombination analysis predicted that HN20du might potentially recombine from the major parent of yk4 (Zhejiang Province, China, 2007) and minor parent of GD/YJ/g2 (Guangdong Province, China, 2020). To the best of our knowledge, this is the first report of GoCV in ducks from China. This broadened host spectrum of GoCVs requires attention from the waterfowl industry and researchers.

Molecular epidemiological investigation and recombination analysis of Cachavirus prevalent in China.

Chaphamaparvovirus carnivoran1 (canine Chaphamaparvovirus, also known as Cachavirus [CachaV]) is a novel parvovirus first reported in dog feces collected from the United States in 2017 and China in 2019. To continuously track its infection and evolution status, 276 canine anal swabs were obtained from pet hospitals in central, northern, and eastern China between 2021 and 2023 and screened via polymerase chain reaction; subsequently, a systematic study was conducted. Of these samples, nine (3.3%) were positive for CachaV. Using polymerase chain reaction, whole genome sequences of the nine CachaV-positive strains were amplified. The NS1 amino acid sequence identity between CachaV strains from China and other countries was 96.23-99.85%, whereas the VP1 protein sequence identity was 95.83-100%. CHN230521 demonstrated the highest identity for NS1 amino acids (99.85%) and VP1 amino acids (100%) with NWT-W88 and CP-T015. According to the model prediction of CHN220916-VP1 protein, Met64Thr, Thr107Ala, and Phe131Ser mutations may cause tertiary structural changes in VP1 protein. Interestingly, each of the nine CachaV strains harbored the same site mutations in NS1 (Ser252Cys, Gly253Leu, and Gly254Thr). Although no explicit recombination events were predicted, the clustering and branching of the phylogenetic tree were complicated. Based on the evolution trees for VP1 and NS1, the nine CachaV strains identified from 2021 to 2023 were closely related to those identified in gray wolves and cats. This study may be beneficial for evaluating the prevalence of CachaVs in China, thereby understanding the evolution trend of CachaVs.

Investigating pigeon circovirus infection in a pigeon farm: molecular detection, phylogenetic analysis and complete genome analysis.

BACKGROUND: Pigeon circovirus infections in pigeons (Columba livia domestica) have been reported worldwide. Pigeons should be PiCV-free when utilized as qualified experimental animals. However, pigeons can be freely purchased as experimental animals without any clear guidelines to follow. Herein, we investigated the status quo of PiCV infections on a pigeon farm in Beijing, China, which provides pigeons for experimental use. RESULTS: PiCV infection was verified in at least three types of tissues in all forty pigeons tested. A total of 29 full-length genomes were obtained and deposited in GenBank. The whole genome sequence comparison among the 29 identified PiCV strains revealed nucleotide homologies of 85.8-100%, and these sequences exhibited nucleotide homologies of 82.7-98.9% as compared with those of the reference sequences. The cap gene displayed genetic diversity, with a wide range of amino acid homologies ranging from 64.5% to 100%. Phylogenetic analysis of the 29 full-genome sequences revealed that the PiCV strains in this study could be further divided into four clades: A (17.2%), B (10.4%), C (37.9%) and D (34.5%). Thirteen recombination events were also detected in 18 out of the 29 PiCV genomes obtained in this study. Phylogenetic research using the rep and cap genes verified the recombination events, which occurred between clades A/F, A/B, C/D, and B/D among the 18 PiCV strains studied. CONCLUSIONS: In conclusion, PiCV infection, which is highly genetically varied, is extremely widespread on pigeon farms in Beijing. These findings indicate that if pigeons are to be used as experimental animals, it is necessary to evaluate the impact of PiCV infection on the results.

Molecular analysis of Gyrovirus galga1 variants identified from the sera of dogs and cats in China.

Gyrovirus galga1 (GyVg1), a member of the Anelloviridae family and Gyrovirus genus, has been detected in chicken and human tissue samples. In this study, the DNA of GyVg1-related gyroviruses in the sera of six dogs and three cats from Central and Eastern China was identified using PCR. Alignment analysis between the nine obtained and reference GyVg1 strains revealed that the genome identity ranged from 99.20% (DOG03 and DOG04 strains) to 96.17% (DOG01 and DOG06 strains). Six recombination events were predicted in multiple strains, including DOG01, DOG05, DOG06, CAT01, CAT02, and CAT03. The predicted major and minor parents of DOG05 came from Brazil. The DOG06 strain is potentially recombined from strains originating from humans and cats, whereas DOG01 is potentially recombined from G17 (ferret-originated) and Ave3 (chicken-originated), indicating that transmissions across species and regions may occur. Sixteen representative amino acid mutation sites were identified: nine in VP1 (12 R/H, 114S/N, 123I/M, 167 L/P, 231 P/S, 237 P/L, 243 R/W, 335 T/A, and 444S/N), four in VP2 (81 A/P, 103 R/H, 223 R/G, and 228 A/T), and three in VP3 (38 M/I, 61 A/T, and 65 V/A). These mutations were only harbored in strains identified in dogs and cats in this study. Whether this is related to host tropism needs further investigation. In this study, GyVg1 was identified in the sera of dogs and cats, and the molecular characteristics prompted the attention of public health.

Molecular Characteristics of Bean Common Mosaic Virus Occurring in Inner Mongolia, China.

Bean common mosaic virus (BCMV) was detected on common bean (Phaseolus vulgaris) plants showing wrinkled and/or narrow leaves, curling, shrinking and chlorosis of leaves, dwarfing of plants, and mottled pods in Inner Mongolia and named BCMV-22Huhe. Its genome has a size of 10,062 bp and was deposited in GenBank under the accession number OR778613. It is closely related to BCMV-Az (GenBank accession no. KP903372, in China) in the lineage of AzBMV. A recombination event was detected for BCMV-22Huhe among the 99 BCMV isolates published in the NCBI GenBank database, showing that BCMV-CJ25 (MK069986, found in Mexico) was a potential major parent, and the minor parent is unknown. This work is the first description of the occurrence of BCMV in Inner Mongolia, China.

Genome-wide analysis of a Coxsackievirus A8 isolate in Yunnan Province / 中国生物制品学杂志

@#Objective To analyze the whole genome sequence and genetic characteristics of the Yunnan isolate of Coxsackievirus A8(CVA8),A8-1/YN/CHN/2019,in order to understand the prevalence and variation of CVA8.Methods The primers for CVA8 were designed,the viral RNA was extracted,the VP1 sequence was amplified by RT-PCR and sequenced,and the whole genome was spliced using BioEdit 7.2.3. MEGA 7.0 software was used for the phylogenetic analysis,and Simplot 3.5.1 and RDP4 software were used for the recombination analysis.Results Strain A8-1/YN/CHN/2019 was7 396 nt long and encoded 2 188 amino acids. The homologies of nucleotide and amino acid sequences between VP1 and CVA8 prototype strain AF081299/Donovan/USA/1998 were 81. 04% and 95. 24%,respectively. Strain A8-1/YN/CHN/2019 belonged to genotype E,while the other Chinese isolates were located in genotypes D and E. Recombination analysis revealed that strain A8-1/YN/CHN/2019 recombined on segments P2 and P3 of the non-structural coding region.Conclusion Strain A8-1/YN/CHN/2019 is of genotype E and has recombination events in the non-structural regions of the P2 and P3 segments.

Molecular identification of carnivore chaphamaparvovirus 2 (feline chaphamaparvovirus) in cats with diarrhea from China.

Chaphamaparvovirus carnivoran2 (feline chaphamaparvovirus, FeChPV) is a novel feline parvovirus originally detected in Canadian cats in 2019, and it has also been identified in domestic cats in other nations. To evaluate the prevalence and genetic diversity of FeChPV in China, rectal swabs of pet cats from Henan, Guangdong, Anhui, Zhejiang, and Inner Mongolia provinces were collected. Of the 230 samples subjected to nested polymerase chain reaction, 6 (2.6%) tested positive for FeChPV. Although all positive samples were from cats with diarrhea, statistical analyses revealed no correlation between the presence of the virus and clinical symptoms (p > 0.05). Phylogenetic trees of nonstructural protein 1 (NS1) and capsid protein (VP1) demonstrated that these six new strains formed a major branch with other reference FeChPV strains and considerably differed from Chaphamaparvoviru carnivoran1. Moreover, recombination analysis revealed that the FeChPV strain CHN20201025, previously detected in a dog, was a recombinant and strains CHN200228 and CHN180917, identified in this study, were the closest relatives to the parental strains. The findings of this study and a previous study wherein FeChPV was detected in dogs suggest that FeChPV can propagate between species. Additionally, these findings indicate that the genetic diversity of FeChPV can provide an insight into the epidemiological status of FeChPV in China.

Complete Genomic Sequence Analysis of a Sugarcane Streak Mosaic Virus Isolate from Yunnan Province of China.

In recent years, the sugarcane streak mosaic virus (SCSMV) has been the primary pathogen of sugarcane mosaic disease in southern China. In this study, the complete genome of a sugarcane mosaic sample (named YN-21) from Kaiyuan City, Yunnan Province, was amplified and sequenced. By comparing the amino acid sequences of YN-21 and 15 other SCSMV isolates from the NCBI database, the protease recognition site of SCSMV was determined. YN-21 had the highest nucleotide and amino acid identities of 97.66% and 99.30%, respectively, in comparison with the SCSMV isolate (JF488066). The P1 had the highest variability of 83.38-99.72% in the amino acid sequence, and 6K2 was the most conserved, with 97.92-100% amino acid sequence identity. A phylogenetic analysis of nucleotide and amino acid sequences clustered the 16 SCSMV isolates into two groups. All the Chinese isolates were clustered into the same group, and YN-21 was closely related to the Yunnan and Hainan isolates in China. Recombination analysis showed no major recombination sites in YN-21. Selective pressure analysis showed that the dN/dS values of 11 proteins of SCSMV were less than 1, all of which were undergoing negative selection. These results can provide practical guidance for monitoring SCSMV epidemics and genetics.

First complete sequences and genetic variation of plum pox virus T strain in Prunus dulcis and Prunus cerasus.

The complete genome of plum pox virus strain T isolates from five different Prunus spp., including almond (P. dulcis) and sour cherry (P. ceracus) isolates, was fully sequenced using the primer pairs designed in this study. The five isolates were aligned with other 50 PPV-T isolates whose complete genome sequences were available in GenBank and then subjected to phylogenetic and diversity analyses. Recombination analysis showed no significant signal detected in the five newly sequenced isolates while confirming four recombinant isolates reported in a previous study. Nucleotide and amino acid phylogenetic trees clustered the tested isolates into three major groups: Balkan 1, 2, and 3. Strain T isolates shared high nucleotide and amino acid identities among them. Diversity analysis applied different parameters to found that the sequences of P3 and 6K1 genes were more conserved over other genes. In accordance, the highly variable P1 and CP genes were found to experience weaker purifying pressures, ω = 0.127 and 0.219, respectively, than other genes. The three neutrality tests gave negative values to all genes, suggesting that strain T populations have expanding or bottleneck selections. Genetic make-up of the only known sour cherry isolate is highly identical to isolates from other Prunus spp. Therefore, this study has updated our knowledge of T strain diversity in new hosts and provided a clear picture of genetic variation and host relationships. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-023-03746-1.

Global Population Structure of Apple Mosaic Virus (ApMV, Genus Ilarvirus).

The gene sequence data for apple mosaic virus (ApMV) in NCBI GenBank were analyzed to determine the phylogeny and population structure of the virus at a global level. The phylogenies of the movement protein (MP) and coat protein (CP) genes, encoded by RNA3, were shown to be identical and consisted of three lineages but did not closely correlate with those of P1 and P2, suggesting the presence of recombinant isolates. Recombination Detection Program (RDP v.4.56) detected significant recombination signal in the P1 region of K75R1 (KY883318) and Apple (HE574162) and the P2 region of Apple (HE574163) and CITH GD (MN822138). Observation on several diversity parameters suggested that the isolates in group 3 had higher divergence among them, compared to isolates in groups 1 and 2. The neutrality tests assigned positive values to P1, indicating that only this region experiencing balanced or contracting selection. Comparisons of the three phylogroups demonstrated high Fixation index (FST) values and confirmed genetic separation and the lack of gene flow among them. Additionally, ±500 bp of partial MP + 'intergenic region' + partial CP coding regions of two Turkish isolates from apple and seven from hazelnut were sequenced and determined that their phylogenetic positions fell within group 1 and 3, respectively.

Whole Genomic Sequence Analysis of Human Adenovirus Species C Shows Frequent Recombination in Tianjin, China.

Human adenovirus species C (HAdV-C) is frequently detected in China and worldwide. For the first time, 16 HAdV-C strains were isolated from sewage water (14 strains) and hospitalised children with diarrhoea (2 strains,) in Tianjin, China. Nearly complete genome data were successfully obtained for these viruses. Subsequently, genomic and bioinformatics analyses of the 16 HAdV-C strains were performed. A phylogenetic tree of the complete HAdV-C genome divided these strains into three types: HAdV-C1, HAdV-C2, HAdV-C5. Phylogenetic analysis based on the fiber gene showed similar outcomes to analyses of the hexon gene and complete HAdV-C genomes, whereas the penton gene sequences showed more variation than previously reported. Furthermore, analysis of the whole-genome sequencing revealed seven recombination patterns transmitted in Tianjin, of which at least four patterns have not been previously reported. However, the penton base gene sequences of the HAdV-C species had significantly lower heterogeneity than those of the hexon and fiber gene sequences of recombinant isolates; that is, many strains were distinct in origin, but shared hexon and fiber genes. These data illustrate the importance of frequent recombination in the complexity of the HAdV-C epidemic in Tianjin, thus emphasising the necessity for HAdV-C sewage and virological monitoring in China.

Genetic characterization of human adenoviruses in patients using metagenomic next-generation sequencing in Hubei, China, from 2018 to 2019.

Objectives: This study aimed to characterize the genomic epidemiology of human adenoviruses (HAdVs) in Hubei, China, using metagenomic next-generation sequencing (mNGS). Methods: In total, 25 HAdV-positive samples collected from 21 pediatric patients were sequenced and subjected to mNGS using the NextSeq 550 and GenoLab M sequencing platforms. The metagenomic data were assembled de novo for molecular typing, phylogenetic and recombination analyzes. Results: We assembled 50 HAdV genomes, 88% (22/25) genomes from GenoLab M, and 84% (21/25) genomes from NextSeq 550 have perfect alignments to reference genomes with greater than 90%. The most fully assembled 25 genomes were categorized into 7 HAdV genotypes, the most abundant of which were HAdV-B3 (9/25) and HAdV-C2 (6/25). Phylogenetic analyzes revealed that the newly isolated HAdV-B3 strains diverged into separate clusters according to their genotypes. Vigilance is needed that HAdV-B3 isolates have begun to form new distinct clusters. High nucleotide identity was observed in the whole genome level within the same HAdV genotypes, while marked differences of three capsid genes across HAdV genotypes were noted. The high nucleotide diversity regions were concordant with the reported hypervariable regions. Further, three recombinant strains were identified: S64 and S71 originated from the parental strains HAdV-B14 and HAdV-B11, and S28 originated from HAdV-C1, HAdV-C5, and HAdV-CBJ113. GenoLab M and NextSeq 550 showed comparable performance with respect to data yield, duplication rate, human ratio, and assembly completeness. Conclusion: The sequencing quality and assembly accuracy showed that mNGS assembled genomes can be used for subsequently HAdV genotyping and genomic characterization. The high nucleotide diversity of capsid genes and high frequency of recombination events has highlighted the necessity for HAdV epidemiological surveillance in China.

Genome analysis and recombination characterization of duck hepatitis B virus isolated from ducks and geese in central China, 2017 to 2019.

Owing to its high similarity to human hepatitis B virus (HBV), duck HBV (DHBV) is often used as an essential model for HBV research. Although intergenotypic recombination of HBV is common, it remains unclear whether the intergenotypic recombination of human HBV is exactly the same as that of DHBV. In this study, 119 serum samples of duck and goose were collected from 51 farms (29 duck and 22 goose farms) in the central and eastern regions of China. A total of 22 strains isolated from the 22 DHBV positive flock were sequenced. Genome sequence alignment revealed that the duck- and goose-origin strains shared the highest and lowest similarities (99.7 and 90.52%, respectively). The complete genomes of these DHBV and 31 reference strains were analyzed using phylogenetic methods and classified into 3 clusters, which corresponded to the previously identified DHBV-I, DHBV-II, and DHBV-III branches. Recombination analyses of the 53 DHBV genomes indicated 2 major intergenotypic recombination events with high confidence values. These recombination events occurred between the genotypes of the Chinese isolates Y180813HB (Chinese branch [DHBV-Ⅰ]) and E170101AH (Chinese branch [DHBV-Ⅱ]) and the Western isolate DHBV-XY (Western branch [DHBV-Ⅲ]), resulting in the emergence of 2 Chinese recombinant isolates Y190303HN and Y170101HB. In addition, 40% (2/5) goose-origin and 58.8% (10/17) duck-origin DHBV in this study harbored the mutation site of G133E in preS, which promote the pathogenicity of DHBV. This is the first study to report on the genome analysis and recombination characterization of DHBV isolated from Chinese geese. Further, continuous investigation and molecular identification of DHBV should be conducted to attract researchers' attention.

Complete genome comparative analysis of Coxsackievirus A6 causing hand,foot and mouth disease and herpangina in Yunnan Province in 2018 / 中国生物制品学杂志

@#ObjectiveTo analyze the genome-wide characteristics of 17 strains of Coxsackievirus A6(CVA6)that cause hand,foot and mouth disease(HFMD)and herpangina(HA)in Yunnan Province in 2018,and understand the genetic differences between different pathogenic CVA6.MethodsA total of 1 909 stool samples clinically diagnosed as HFMD and HA in Kunming Children′s Hospital in 2018 were randomly selected for detection using enterovirus group A universal primers and screening of CVA6 positive samples. The CVA6 whole genome sequence was amplified with CVA6 whole genome primers,spliced by BioEdit splicting software,and analyzed for the whole genome characteristics by BioEdit,MEGA 7.0,Simplot,Heml 1.0 and Phyre2softwares.ResultsA total of 929 CVA6 positive samples were screened,and 17CVA6 complete gene sequences were obtained(9 of which were clinically diagnosed as HFMD and 8 were clinically diagnosed as HA). All 17 CVA6 strains were in type IV clade on the whole phylogenetic tree. No significant recombination occurred in HA and HFMD representative strains,while mutations occurred in non-structural protein 3D region. HFMD and HA representative strains showed differences in VP1 loci S597T,Q705L and Q663L. Online predictive analysis showed that the secondary structure of VP1 was consistent with that of CVA6 with no change.ConclusionThe 17 CVA6 strains causing HFMD and HA had high genomic homology,as well as nucleotide and amino acid differences,which may affect the replication and adaptability of CVA6.

Molecular Epidemiology Reveals the Co-Circulation of Two Genotypes of Coxsackievirus B5 in China.

Coxsackievirus B5 (CVB5) is an important enterovirus B species (EV-Bs) type. We used the full-length genomic sequences of 53 viral sequences from the national hand, foot, and mouth disease surveillance network in the Chinese mainland (2001-2021). Among them, 69 entire VP1 coding region nucleotide sequences were used for CVB5 genotyping and genetic evolution analysis. Phylogenetic analysis based on a data set of 448 complete VP1 sequences showed that CVB5 could be divided into four genotypes (A-D) worldwide. Sequences from this study belonged to genotypes B and D, which dominated transmission in the Chinese mainland. Two transmission lineages of CVB5 have been discovered in the Chinese mainland, lineage 2 was predominant. Markov chain Monte Carlo analysis indicated that the tMRCA of CVB5 in the Chinese mainland could be traced to 1955, while the global trend could be traced to 1862, 93 years earlier than China. The evolution rate of CVB5 was higher in the Chinese mainland than worldwide. The spatiotemporal dynamics analysis of CVB5 assessed that virus transportation events were relatively active in Central, Northeast, North and Northwest China. Recombination analysis revealed frequent intertypic recombination in the non-structural region of CVB5 genotypes B and D with the other EV-Bs, revealing eight recombination lineages. Our study showed the molecular evolution and phylogeography of CVB5 that could provide valuable information for disease prevention.

Molecular evolution of human coronavirus-NL63, -229E, -HKU1 and -OC43 in hospitalized children in China.

Human coronaviruses (HCoVs) HCoV-NL63, HCoV-229E, HCoV-HKU1 and HCoV-OC43 have been circulated in the human population worldwide, and they are associated with a broad range of respiratory diseases with varying severity. However, there are neither effective therapeutic drugs nor licensed vaccines available for the treatment and prevention of infections by the four HCoVs. In this study, we collected nasopharyngeal aspirates of children hospitalized for respiratory tract infection in China during 2014-2018 and conducted next-generation sequencing. Sequences of four HCoVs were then selected for an in-depth analysis. Genome sequences of 2 HCoV-NL63, 8 HCoV-229E, 2 HCoV-HKU1, and 6 HCoV-OC43 were obtained. Based on the full-length S gene, a strong temporal signal was found in HCoV-229E and the molecular evolutionary rate was 6 × 10-4 substitutions/site/year. Based on the maximum-likelihood (ML) phylogenetic tree of complete S gene, we designated H78 as a new sub-genotype C2 of HCoV-HKU1, and the obtained P43 sequence was grouped into the reported novel genotype K of HCoV-OC43 circulating in Guangzhou, China. Based on the complete genome, potential recombination events were found to occur as two phenomena, namely intraspecies and interspecies. Moreover, we observed two amino acid substitutions in the S1 subunit of obtained HCoV-NL63 (G534V) and HCoV-HKU1 (H512R), while residues 534 and 512 are important for the binding of angiotensin-converting enzyme 2 and neutralizing antibodies, respectively. Our findings might provide a clue for the molecular evolution of the four HCoVs and help in the early diagnosis, treatment and prevention of broad-spectrum HCoV infection.

Genetic and Cross Neutralization Analyses of Coxsackievirus A16 Circulating in Taiwan from 1998 to 2021 Suggest Dominant Genotype B1 can Serve as Vaccine Candidate.

Coxsackievirus A16 (CVA16) is well known for causing hand-foot-and-mouth disease (HFMD) and outbreaks were frequently reported in Taiwan in the past twenty years. The epidemiology and genetic variations of CVA16 in Taiwan from 1998 to 2021 were analyzed in this study. CVA16 infections usually occurred in early summer and early winter, and showed increased incidence in 1998, 2000-2003, 2005, 2007-2008, and 2010 in Taiwan. Little or no CVA16 was detected from 2017 to 2021. CVA16 infection was prevalent in patients between 1 to 3 years old. A total of 69 isolates were sequenced. Phylogenetic analysis based on the VP1 region showed that CVA16 subgenotype B1 was dominantly isolated in Taiwan from 1998 to 2019, and B2 was identified only from isolates collected in 1999 and 2000. There was a high frequency of synonymous mutations in the amino acid sequences of the VP1 region among CVA16 isolates, with the exception of position 145 which showed positive selection. The recombination analysis of the whole genome of CVA16 isolates indicated that the 5'-untranslated region and the non-structural protein region of CVA16 subgenotype B1 were recombined with Coxsackievirus A4 (CVA4) and enterovirus A71 (EVA71) genotype A, respectively. The recombination pattern of subgenotype B2 was similar to B1, however, the 3D region was similar to EVA71 genotype B. Cross-neutralization among CVA16 showed that mouse antisera from various subgenotypes viruses can cross-neutralize different genotype with high neutralizing antibody titers. These results suggest that the dominant CVA16 genotype B1 can serve as a vaccine candidate for CVA16.