RDP4: Detection and analysis of recombination patterns in virus genomes

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.

Isolation and characterization of a novel S1-gene insertion porcine epidemic diarrhea virus with low pathogenicity in newborn piglets

Porcine epidemic diarrhea virus (PEDV) causes diarrhea and vomiting in piglets, leading to a mortality rate of 100%. Due to the high frequency of mutation, it is important to monitor the evolution of PEDV and develop potential vaccine candidates. In this study, two PEDV strains (ZJ2022 and ZQ2022) were identified by PCR. These strains were subsequently isolated, and their genome sequences, growth characteristics, and pathogenicity were compared. Phylogenetic and recombination analyses revealed that both strains belonged to GIIa-subgroup, and ZQ2022 was identified as a recombinant strain derived from ZJ2022. Further sequence analysis showed that the ZJ2022 strain had a modified top region of the S1 protein due to a three amino acid insertion (T380_Y380insGGE) in the S1 gene. According to the virus growth curve, ZJ2022 exhibited better cellular adaptation than ZQ2022, with higher viral titers from 8 hpi to 24 hpi. Additionally, ZQ2022 exhibited a high level of pathogenicity, causing severe diarrhea in piglets at 36 hpi and a 100% mortality rate by 96 hpi. In contrast, ZJ2022 showed lower pathogenicity, inducing severe diarrhea in piglets at 60 hpi, with a mortality rate of 60% at 96 hpi and 100% at 120 hpi. In summary, our findings provided evidence of the undergoing mutations in Chinese PEDV strains. Furthermore, the S gene insertion strain ZJ2022 exhibited strong cellular adaptability and low pathogenicity, making it a potential candidate strain for vaccine development.

Canine parvovirus in North-East India: a phylogenetic and evolutionary analysis

Canine parvovirus type 2 (CPV-2) infection in dogs is considered as one of the most common cause of morbidity and mortality in young dogs and continues to occur with high incidence worldwide. Despite a single-stranded DNA virus, CPV-2 possesses a high mutation rate which has led to the development of new variants from time to time. These variants are classically classified based on the amino acid markers present in the VP2 gene. In this study, we examined 20 different cases of CPV-2 infection from seven different states of the North East region (NER) of India. The near-complete genome sequences of all these isolates were subjected to phylodynamic and phylogeographic analysis to evaluate the genetic diversity and geographical spread of CPV-2 variants. Analysis of the deduced amino acid sequences revealed residues characteristic of the 'Asian CPV-2c lineage' in all the 20 sequences confirming it as the dominant strain circulating in NER, India. The phylogenetic analysis based on the whole genome showed that all 20 sequences formed a monophyletic clade together with other Asian CPV-2c sequences. Furthermore, phylogeographic analysis based on the VP2 gene showed the likely introduction of Asian CPV-2c strain to India from China. This study marks the first comprehensive report elucidating the molecular epidemiology of CPV-2 in India.

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.

Evidence of circulating recombinants between deformed wing virus and Varroa destructor virus-1 in honey bee colonies in Türkiye

Deformed wing virus (DWV), which is an important honey bee virus transmitted by Varroa destructor (V. destructor), causes colony losses in honey bee colonies. This study aimed to investigate the prevalence and genetic diversity of DWV in honey bees in Türkiye and to determine the role of V. destructor in the transmission of the genetic variants of DWV. Honey bee samples were collected from 62 apiaries, by simple random sampling, during March 2022 and April 2023. The presence of V. destructor in collected bee samples was examined using a stereo microscope. Real-time RT-PCR was used for the detection of DWV-A and DWV-B (Varroa destructor virus-1 (VDV-1)) viruses. Genetic characterisation of the positive samples was conducted by sequencing polyprotein genomic region. Considering the V. destructor infestation rate of 3% as relevant, out of the 62 apiaries examined, 17 (27.4%) were positive. However, DWV-A and VDV-1 specific RNA was not detected in V. destructor samples. VDV-1 specific RNA was detected in 6.5% (4/62) of the apiaries, whereas DWV-A was not detected in the sampled apiaries. Phylogenetic analysis showed that isolates detected in this study were located in a separate cluster from previously characterised DWV-A and VDV-1 isolates. According to RDP4 and GARD analyses, DWV-VDV-1 recombination breakpoints were detected in field isolates. To the best our knowledge, this is the first report of the presence of VDV-1-DWV recombinants in Türkiye. Further studies are needed to determine the impact of VDV-1-DWV recombinants and their virological and antigenic properties.

Genetic characterization of immune adaptor molecule MyD88 in Culex pipiens complex (Diptera: Culicidae) mosquitoes from China

Mosquitoes of the Culex (Cx.) pipiens complex are vectors of severe diseases including West Nile fever by West Nile virus, Japanese encephalitis by Japanese encephalitis virus, and Lymphatic filariasis by filarial nematode Wuchereria bancrofti. As a major portion of mosquito immune system, the Toll pathway implicates in response against infections of mosquito-borne pathogens and biocontrol agents. The genetic diversity of immune-related molecules is expected to be a feasible and effective introduction to expand our knowledge of the mosquito-microbe interplay. However, a comprehensive description is currently lacking regarding the genetic characteristic of the Toll pathway molecules in Cx. pipiens complex mosquitoes. In the present study, genetic changes in Cx. pipiens complex MyD88 (Myeloid differentiation primary response protein 88) were analyzed as a precedent for the Toll pathway molecules in this taxon. MyD88 is a critical adaptor of the pathway transducing signals from TIR-containing receptors to downstream death domain-containing molecules. Our results revealed that adaptive selection has influenced the genetic changes of the molecule, giving rise to acceleration of diversity at a number of amino acid sites. The adaptively selected sites lie in the death domain, intermediate domain, and C-terminal extension. The characteristics of the genetic changes shed insights into the prominent molecular-level structural basis and the involvement strategy of the adaptor in the arms race against exogenous challenges. This finding would be beneficial for further exploration and deeper understanding of the mosquitoes' vectorial capacity and facilitating the effectiveness and sustainability of the biocontrol agents.

Novel imported clades accelerated the RSV surge in Beijing, China, 2023-2024

OBJECTIVES

Despite the optimization of the zero-COVID policy in late 2022, there was a subsequent increase noted in the number of respiratory syncytial virus (RSV) cases in Northern China. In this study, we investigated and characterized the dynamics of this surge at the genomic level in Beijing, China.

METHODS

Patients with acute respiratory tract infections (ARTIs) were enrolled from 35 sentinel hospitals in Beijing, China. Epidemiological investigations, G gene amplification, and whole-genome sequencing were performed, followed by epidemiological analysis, imported clade detection, and mutation identification. We also combined global data to illustrate the biological and epidemiological characteristics of the emerging clades.

RESULTS

A total of 60,423 patients with ARTIs were recruited between January 2015 and January 2024. The RSV peak observed in the winter of 2023 was the highest in the past 9 years. Two novel imported clades, A.D.5.2 and B.D.E.1, were detected for the first time in China. This surge was mainly driven by B.D.E.1, which exhibited a significantly higher proportion of older individuals both in Beijing and globally. Seven non-synonymous mutations in B.D.E.1 were found in Beijing. B.D.E.1 had more sites suffering from positive selection than its parent.

CONCLUSIONS

The novel imported clade B.D.E.1 accelerated an unprecedented RSV surge in Beijing, presenting noteworthy epidemiological and biological characteristics. Continuous RSV genome monitoring has important implications for RSV outbreak identification, genetic diversity tracking, vaccine development, and strategy implementation.

Molecular detection and genetic characteristics of porcine reproductive and respiratory syndrome virus in central China

Porcine reproductive and respiratory syndrome caused by porcine reproductive and respiratory syndrome virus (PRRSV) is one of the most economically devastating viral diseases in the global pork industry. To further clarify the epidemic characteristics of the virus, 365 clinical samples were collected from diseased pigs suffering from abortion and respiratory disease from 2018 to 2023 on 63 pig farms in Henan and Shanxi provinces, and screened for the presence of PRRSV using reverse transcription-polymerase chain reaction (RT-PCR). A total of 62 clinical samples (62/365, 16.99 %) were positive for PRRSV, and subsequently, full-length ORF5 gene sequences of 29 PRRSV strains and the complete genome sequence of one PRRSV HeN-HC isolate were obtained and analyzed. Phylogenetic analysis based on the ORF5 gene showed that 22 of the 29 PRRSV2 strains belonged to sublineage 1.8 (NADC30-like), 5 belonged to sublineage 8.5 (HP-PRRSV), and 2 belonged to sublineage 5.1 (VR-2332-like), indicating that both HP-PRRSV and NADC30-like strains were mainly circulating in Henan and Shanxi provinces. Compared to VR-2332 strain, different types of amino acid mutations were found in the GP5 protein of these 29 strains, and the amino acid deletions were displayed in the Nsp2 protein of the HeN-HC isolate, leading to the variation of protein structures. It is noteworthy that recombination events were identified in the HeN-Ping and HeN-B strains. In addition, a total of 60, 094 pig serum samples from Henan province were collected, and the positive rate of specific antibodies against PRRSV was 86.37 % from 2019 to 2022, and 86.66 %, 84.85 %, 87.54 % and 86.30 % in 2019, 2020, 2021 and 2022, respectively. Overall, this study provides valuable insights into the molecular epidemiology and evolution of PRRSV circulating in central China.

Genotyping and phylogeographic dynamics of coxsackievirus A16

Coxsackievirus A16 (CV-A16) is one of the major pathogens of Hand, Foot and Mouth disease. Here, we analyzed 287 full-length genome sequences of CV-A16 found worldwide from 1994 to 2019 to see the genomic evolution characteristics. Full-length genome-based phylogenetic tree divided the viruses into five different genotypes, G-a to G-e. The CV-A16 strains circulating in China dominate G-a and G-c, but can also be found in other genotypes including G-b and G-e. Phylogeographic analysis showed a high diversity of CV-A16 distribution. In addition, recombination was shown to drive the genomic evolution of CV-A16 during past decades. However, the structural proteins still remain relative conserved while there is extensive genomic recombination. This study updates the phylogenetic and phylogeographic information of CV-A16 and provides insights into the genetic characteristics of CV-A16.

Domestic cat hepadnavirus genotype B is present in Southern Brazil

Domestic cat hepadnavirus (DCH) (Orthohepadnavirus felisdomestici) is an emerging virus related to the hepatitis B virus (HBV) already reported in many countries. The molecular prevalence of DCH varies widely in the regions investigated so far. In the present work, we reported the presence of DCH in Brazil. Sixty cat serum samples tested by DCH presence using PCR and 1.67% (1/60) were positive, similar to the low positive molecular rates reported in United States and Japan. The DCH full-length genome was classified in genotype B, which is uncommon since this genotype was only reported once in Japan. The DCH-positive sample was obtained in a stray cat female apparently healthy, presenting ALT, AST, and ALKP normal values, and negative for FIV and FeLV. Due the low positivity rate detected, some factors as alteration in hepatic enzymes and FIV/FeLV infection could not be evaluated. Other works are necessary to statistically validate these observations in Brazil.

Genomic characteristics of human respiratory syncytial virus from children in China during 2017-2020

Acute lower respiratory tract infections (ALRTIs) are a leading cause of mortality in young children worldwide due to human respiratory syncytial virus (RSV). The aim of this study was to monitor genetic variations in RSV and provide genomic data support for RSV prevention and control. A total of 105 complete RSV genome sequences were determined during 2017-2020. Phylogenetic analysis showed that all of the RSVA sequences were of genotype ON1, and all of the RSVB sequences were of genotype BA9. Notably, a phylogenetic tree based on the whole genome had more branches than a tree based on the G gene. In comparison to the RSV prototype sequences, 71.43% (50/70) of the ON1 sequences had five amino acid substitutions (T113I, V131N, N178G, H258Q, and H266L) that occurred simultaneously, and 68.57% (24/35) of the BA9 genotype sequences had 12 amino acid substitutions, four of which (A131T, T137I, T288I, and T310I) occurred simultaneously. In the F gene, there were 19 amino acid substitutions, which were mainly located in the antigenic sites Ø, II, V, and VII. Other amino acid substitutions were found in the NS1, NS2, P, SH, and L proteins. No significant evidence of recombination was found in any of the sequences. These findings provide important data that will be useful for prevention, control, and vaccine development against RSV.

Genomic surveillance reveals a dengue 2 virus epidemic lineage with a marked decrease in sensitivity to Mosnodenvir

Dengue fever is the most important arbovirosis for public health, with more than 5 million cases worldwide in 2023. Mosnodenvir is the first anti-dengue compound with very high preclinical pan-serotype activity, currently undergoing phase 2 clinical evaluation. Here, by analyzing dengue virus (DENV) genomes from the 2023-2024 epidemic in the French Caribbean Islands, we show that they all exhibit mutation NS4B:V91A, previously associated with a marked decrease in sensitivity to mosnodenvir in vitro. Using antiviral activity tests on four clinical and reverse-genetic strains, we confirm a marked decrease in mosnodenvir sensitivity for DENV-2 ( > 1000 fold). Finally, combining phylogenetic analysis and experimental testing for resistance, we find that virus lineages with low sensitivity to mosnodenvir due to the V91A mutation likely emerged multiple times over the last 30 years in DENV-2 and DENV-3. These results call for increased genomic surveillance, in particular to track lineages with resistance mutations. These efforts should allow to better assess the activity profile of DENV treatments in development against circulating strains.

Full Genome Characterization of the First Oropouche Virus Isolate Imported in Europe from Cuba

On 27 May 2024, the Cuban Ministry of Health reported the first outbreak of Oropouche fever on the island. The etiologic agent, Oropouche virus (OROV), is a poorly understood arbovirus that has been known since the 1960s and represents a public health burden in Latin America. We report the whole-genome characterization of the first European OROV isolate from a returning traveler from Cuba with Oropouche fever-like symptoms. The isolate was obtained from the patient's serum; whole-genome sequencing was performed by next-generation sequencing, followed by phylogenetic analysis and genetic variability studies. The analysis showed that the most closely related sequence was from the French Guiana 2020 outbreak. Interestingly, our isolate is a reassortant virus, included in a highly supported monophyletic clade containing recent OROV cases (Brazil 2015-Colombia 2021), separated from the other four previously known genotypes. More deeply, it was found to be included in a distinct branch containing the sequences of the Brazil 2022-2024 outbreak. The reassortment event involved the S and L segments, which have high similarity with sequences belonging to a new cluster (here defined as OROV_SCDC_2024), while the M segment shows high similarity with older sequences. These results likely describe the viral strain responsible for the current outbreak in Cuba, which may also reflect the ongoing outbreak in Latin America. Further studies are needed to understand how OROV evolves towards traits that facilitate its spread and adaptation outside its original basin, and to track its spread and evolution in the European continent.

Landscape Heterogeneity Drives Genetic Diversity in the Highly Dispersive Moss Funaria hygrometrica Hedw

Funaria hygrometrica, a cosmopolitan moss species known for its remarkable dispersal capacity, was selected as the focal organism to investigate the relationship between landscape features and genetic diversity. Our study encompassed samples collected from two distinct regions: the Spanish Sierra Nevada Mountains (SN), characterized by a diverse landscape with an altitudinal difference of nearly 3500 m within a short distance, and the Murcia Region (MU) in Southeast Spain, characterized by a uniform landscape akin to the lowlands of Sierra Nevada. Genotyping analysis targeted three genetic regions: the nuclear ribosomal internal transcribed spacer (nrITS), the chloroplast rps3-rpl16 region, and the mitochondrial rpl5-rpl16 spacer. Through this analysis, we aimed to assess genetic variability and population structure across these environmentally contrasting regions. The Sierra Nevada populations exhibited significantly higher haplotype diversity (Hd = 0.78 in the highlands and 0.67 overall) and nucleotide diversity (π% = 0.51 for ITS1) compared to the Murcia populations (Hd = 0.35, π% = 0.14). Further investigation unveiled that samples from the lowlands of Sierra Nevada showed a closer genetic affinity to Murcia than to the highlands of Sierra Nevada. Furthermore, the genetic differentiation between highland and lowland populations was significant (ΦST = 0.55), with partial Mantel tests and ResistanceGA analysis revealing a strong correlation between ITS1-based genetic diversity and landscape features, including altitude and bioclimatic variables. Our study elucidated potential explanations for the observed genetic structuring within F. hygrometrica samples' populations. These included factors such as a high selfing rate within restricted habitats, a limited average dispersal distance of spores, hybrid depression affecting partially incompatible genetic lineages, and recent migration facilitated via human activities into formerly unoccupied areas of the dry zones of Southeast Spain.

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 characteristics of human parainfluenza viruses 1-4 associated with acute lower respiratory tract infection in Chinese children, during 2015-2021

Human parainfluenza viruses (HPIVs) are a significant cause of acute lower respiratory tract infections (ALRTIs) among young children and elderly individuals worldwide. The four types of HPIVs (HPIV1-4) can cause recurrent infections and pose a significant economic burden on health care systems globally. However, owing to the limited availability of complete genome sequences, the genetic evolution of these viruses and the development of vaccines and antiviral treatments are hampered. To address this issue, this study utilized next-generation sequencing to obtain 156 complete genome sequences of HPIV1-4, which were isolated from hospitalized children with ALRTIs in six regions of China between 2015 and 2021. This study revealed multiple clades, lineages, or sublineages of HPIVs circulating in mainland China, with a novel clade D of HPIV1 identified as geographically restricted to China. Moreover, this study identified the endemic dominant genotype of HPIV3, lineage C3, which has widely spread and continuously circulated in China. Bioinformatic analysis of the genome sequences revealed that the proteins of HPIV3 possessed the most variable sites, with the P protein showing more diversity than the other proteins among all types of HPIVs. The HN proteins of HPIV1-3 are all under negative/purifying selection, and two amino acid substitutions in the HN proteins correspond to known mAb neutralizing sites in the two HPIV3 strains. These findings provide crucial insights into the genetic diversity and evolutionary dynamics of HPIVs circulating among children in China and may facilitate research on the molecular diagnosis, vaccine development, and surveillance of HPIVs.IMPORTANCEPhylogenetic analysis revealed the prevalence of multiple clades, lineages, or sublineages of human parainfluenza viruses (HPIVs) circulating in mainland China. Notably, a unique evolutionary branch of HPIV1 containing only Chinese strains was identified and designated clade D. Furthermore, in 2023, HPIV3 strains from Pakistan and Russia formed a new lineage within clade C, named C6. The first HPIV4b sequence obtained in this study from China belongs to lineage C2. Evolutionary rate assessments revealed that both the HN and whole-genome sequences of HPIV3 presented the lowest evolutionary rates compared with those of the other HPIV types, with rates of 6.98E-04 substitutions/site/year (95% HPD: 5.87E-04 to 8.25E-03) and 5.85E-04 substitutions/site/year (95% HPD: 5.12E-04 to 6.62E-04), respectively. Recombination analysis revealed a potential recombination event in the F gene of an HPIV1 strain in this study. Additionally, all the newly obtained HPIV1-3 strains exhibited negative selection pressure, and two mutations were identified in the HN protein of two HPIV3 strains at monoclonal antibody-binding sites.

The complete genome of equid herpesvirus-1 (EHV-1) field isolates from Argentina reveals an interspecific recombinant strain

The Equid alphaherpesvirus type 1 (EHV-1) infection can have devastating economic consequences in the horse industry due to large-scale outbreaks of abortions, perinatal foal mortality, and myeloencephalopathy. The present study analyzed the genome of two isolates obtained from aborted fetuses in Argentina, E/745/99 and E/1297/07. The E745/99 genome shares 98.2% sequence identity with Ab4, a reference EHV-1 strain. The E/1297/07 genome shares 99.8% identity with NY03, a recombinant strain containing part of ORF64 and part of the intergenic region from Equid alphaherpesvirus-4 (EHV-4). The E/1297/07 genome has the same breakpoints as other United States and Japanese recombinants, including NY03. The recombinant regions have varying numbers of tandem repeat sequences and different minor parental sequences (EHV-4), suggesting distinct origins of the recombinant events. These are the first complete genomes of EHV-1 from Argentina and South America available in the Databases.

Farmed fur animals harbour viruses with zoonotic spillover potential

Animals such as raccoon dogs, mink and muskrats are farmed for fur and are sometimes used as food or medicinal products1,2, yet they are also potential reservoirs of emerging pathogens3. Here we performed single-sample metatranscriptomic sequencing of internal tissues from 461 individual fur animals that were found dead due to disease. We characterized 125 virus species, including 36 that were novel and 39 at potentially high risk of cross-species transmission, including zoonotic spillover. Notably, we identified seven species of coronaviruses, expanding their known host range, and documented the cross-species transmission of a novel canine respiratory coronavirus to raccoon dogs and of bat HKU5-like coronaviruses to mink, present at a high abundance in lung tissues. Three subtypes of influenza A virus-H1N2, H5N6 and H6N2-were detected in the lungs of guinea pig, mink and muskrat, respectively. Multiple known zoonotic viruses, such as Japanese encephalitis virus and mammalian orthoreovirus4,5, were detected in guinea pigs. Raccoon dogs and mink carried the highest number of potentially high-risk viruses, while viruses from the Coronaviridae, Paramyxoviridae and Sedoreoviridae families commonly infected multiple hosts. These data also reveal potential virus transmission between farmed animals and wild animals, and from humans to farmed animals, indicating that fur farming represents an important transmission hub for viral zoonoses.

Coinfection of human adenovirus and recombinant human astrovirus in a case of acute gastroenteritis: A report from China

Diarrhea is one of the major public health issues worldwide. Although the infections of individual enteric virus have been extensively studied, elucidation of the coinfection involving multiple viruses is still limited. In this study, we identified the coinfection of human adenovirus (HAdV) and human astrovirus (HAstV) in a child with acute gastroenteritis, analyzed their genotypes and molecular evolution characteristics. The sample was collected and identified using RT-PCR and subjected to whole-genome sequencing on the NovaSeq (Illumina) platform. Obtained sequences were assembled into the complete genome of HAdV and the ORF1 of HAstV. We conducted phylogenetic analysis using IQ-TREE software and conducted recombination analysis with the Recombination Detection Program. The sequenced HAdV was confirmed to be genotype 41, and was genetically close to some European strains. Phylogenetic analysis revealed that the HAstV was genetically close to both HAstV-2 and HAstV-4 and was different from the genotype prevalent in Shenzhen before. The recombination analysis confirmed that the sequenced HAstV strain is a recombinant of HAstV-2 and HAstV-4. Our analysis has shown that the strains in this coinfection are both uncommon variants in this geographical region, instead of dominant subtypes that have prevailed for years. This study presents a coinfection of HAdV and HAstV and conducts an evolutionary analysis on involved viruses, which reveals the genetic diversity of epidemic strains in Southern China and offers valuable insights into vaccine and medical research.

Comparative analysis of two novel complete genomes of myxoma virus vaccine strains

Myxoma virus (MYXV) is a double-stranded DNA-containing virus of the family Poxviridae, genus Leporipoxvirus. MYXV is an important model virus for evolutionary and immunological research and a promising oncolytic. In this study, we sequenced and analyzed two complete genomes of MYXV virus vaccine strains B-82 and Rabbivac-B, which are widely used for vaccine production in Russia. Here, we first show that MYXV vaccine strains B-82 and Rabbivac-B share a common origin with the American recombinant MYXV MAV vaccine strain. In addition, our data suggest that the MYXV B-82 and Rabbivac-B strains contain a number of genes at the 5' and 3' ends that are identical to the virulent MYXV Lausanne strain. Several unique genetic signatures were identified in the M013L, M017L, M023, and M121R genes, helping to achieve high genetic resolution between vaccine strains. Overall, these findings highlight the evolutionary flexibility of certain genes in the MYXV genome and provide insights into the molecular epidemiology of the virus and subsequent vaccine development.

Identification of a novel monopartite begomovirus associated with leaf curl disease of Citharexylum spinosum in India

The present study reports the complete genome of a novel monopartite begomovirus, named tentatively as "Citharexylum leaf curl virus" (CitLCuV), associated with leaf curl disease of Citharexylum spinosum in India. CitLCuV genome (2767 nucleotide) contained the typical genome organization of Old World begomoviruses and shared the maximum nucleotide sequence identity of 89.7% with a papaya leaf crumple virus (PaLCrV) isolate. In addition, two small non-canonical open reading frames (C5 and C6) were determined in the complementary strand of CitLCuV genome. Phylogenetic analysis revealed the relatedness of CitLCuV to PaLCrV and rose leaf curl virus. Recombination analysis detected a possible recombination event in CitLCuV genome. Based on begomovirus species demarcation criteria, CitLCuV can be regarded as a novel begomoviral species.

Surveillance of Wildlife Viruses: Insights from South Australia's Monitoring of Rabbit Haemorrhagic Disease Virus (RHDV GI.1 and GI.2)

Surveillance of wildlife virus impacts can be passive or active. Both approaches have their strengths and weaknesses, especially regarding cost and knowledge that can be gained. Monitoring of rabbit haemorrhagic disease virus (GI.1 and GI.2) in South Australia has utilised both strategies and their methods and gained insights are discussed. Active strategies to monitor the continuing impact of rabbit haemorrhagic disease virus 2 (GI.2) on susceptible lagomorphs in countries such as the USA, Mexico, South Africa, Spain, France and Portugal are encouraged to gain critical insights into the evolution, spread and impact of this virus. Furthermore, there are lessons here for the international monitoring of diseases in wildlife, particularly where there is a risk of them becoming zoonotic.

Inferring the origin of new D-loop haplotypes of loggerhead sea turtles (Testudinata: Cheloniidae) from the Southwest Atlantic lineage

The populations of the loggerhead turtles, Caretta caretta, present four main D-loop mitochondrial haplogroups that are distributed across the Indo-Pacific, Mediterranean, and Atlantic oceans. The Southwestern Atlantic (SWA) is one of the Regional Management Units (RMUs) of loggerheads, characterized by unique haplotypes, high nest density, and distinct life history traits. Detecting new D-loop haplogroups is important, particularly endemic ones, as they can enhance our understanding of their life history within the RMUs and contribute to the resolution of mixed stock analysis. In this study, we conducted a series of phylogenetic delimitation and network analyses to identify, validate, and infer the origin of four new D-loop haplotypes detected in the loggerhead populations from the SWA. Our findings demonstrate that these new D-loop haplotypes are valid and unique to the SWA lineage, potentially aiding in the delimitation of individuals' origins and the inference of their lineage.

Feline bocaviruses found in Thailand have undergone genetic recombination for their evolutions

Feline bocaviruses (FBoVs) have been discovered for a decade and are often detected in feces, possibly associated with diarrhea in cats. Studies on FBoV evolution remain limited and have mainly focused on prevalence and genetic characterization. Although genetic recombination serves as a potential mechanism in bocavirus evolution, research on this process for FBoVs has been scarce. In this study, we characterized 19 complete coding sequences of FBoVs obtained from Thai cats, revealing that FBoV-1, -2, and -3 were endemic in Thailand. Genetic characterizations showed that most Thai FBoVs were closely related to previously detected strains in Thailand and China. Recombination analyses indicated intragenic, intraspecies recombination in all FBoV species, with recombination breakpoints commonly found in the NP1 and VP1/2 genes, highlighting these genes may be hotspots for FBoV recombination. However, no interspecies recombination was detected. Selective pressure analysis of various FBoV genes revealed that these viruses underwent purifying selection. Although the VP1/2 gene of all FBoV species was under the strongest negative selection pressure, positive selection sites were only found in FBoV-1 and FBoV-3. This study is the first to identify natural recombination in FBoV-2 and FBoV-3 and provides evidence that genetic recombination is a potential driver of FBoV evolutions. Additionally, this study offers up-to-date information on the genetic characteristics, evolutionary dynamics, and selective pressure status of FBoVs, which should be continuously monitored.

Next Generation Sequencing and Genetic Analyses Reveal Factors Driving Evolution of Sweetpotato Viruses in Uganda

Sweetpotato (Ipomoea batatas L.) is an essential food crop globally, especially for farmers facing resource limitations. Like other crops, sweetpotato cultivation faces significant production challenges due to viral infections. This study aimed to identify and characterize viruses affecting sweetpotato crops in Uganda, mostly those associated with sweetpotato virus disease (SPVD). Infected leaf samples were collected from farmers' fields in multiple districts spanning three regions in Uganda. MiSeq, a next-generation sequencing platform, was used to generate reads from the viral nucleic acid. The results revealed nine viruses infecting sweetpotato crops in Uganda, with most plants infected by multiple viral species. Sweet potato pakakuy and sweet potato symptomless virus_1 are reported in Uganda for the first time. Phylogenetic analyses demonstrated that some viruses have evolved to form new phylogroups, likely due to high mutations and recombination, particularly in the coat protein, P1 protein, cylindrical inclusion, and helper component proteinase regions of the potyvirus. The sweet potato virus C carried more codons under positive diversifying selection than the closely related sweet potato feathery mottle virus, particularly in the P1 gene. This study provides valuable insights into the viral species infecting sweetpotato crops, infection severity, and the evolution of sweet potato viruses in Uganda.

Whole genome characteristics of hedgehog coronaviruses from Poland and analysis of the evolution of the Spike protein for its interspecies transmission potential

BACKGROUND

The hedgehogs have been recently identified as possible reservoir of Middle East respiratory syndrome coronavirus like (MERS-CoV-like). These viruses were classified as a distinct Betacoronavirus erinacei (BCoV-Eri) species within the MerBCoV-Eriirus subgenus. As coronaviruses are known for their ability to jump between different hosts, including humans, this can pose a particular threat to people in direct contact with hedgehogs, such as those working at animal asylums. Our previous studies have shown the presence of BCoV-Eri strains in animals collected in the wildlife rehabilitation centre. This study aimed to investigate the presence of CoV in subsequent hedgehogs collected from the urban area of Poland and their molecular characteristics.

RESULTS

Monitoring for the presence of coronavirus infection in hedgehogs revealed five positive individuals. The presence of BCoV-Eri was found in a total of 20% of animals tested. Our analyses revealed no correlation between CoVs positivity and animal health conditions but a higher probability of such infection in juveniles and females. The whole genome of two Polish Hedgehog coronavirus 1 strains were sequenced and compared with available counterparts from European and Asian countries. Phylogenetic analysis showed that both CoV strains formed common cluster with other similar MerBCoV-Eriirus, but they were also found to be genetically variable and most changes in the S protein were identified. Our analysis revealed that some S protein sites of the Hedgehog coronavirus 1 strains evolved under positive selection pressure and of five such sites, three are in the S1 region while the other two in the S2 region of the Spike.

CONCLUSIONS

BCoV-Eri is to some extent prevalent in wildlife asylums in Poland. Given that the S protein of BCoVs-Eri is highly variable and that some sites of this protein evolve under positive selection pressure, these strains could potentially acquire a favourable feature for cross-species transmission. Consequently, the threat to humans working in such asylums is particularly high. Adequate biosecurity safeguards, but also human awareness of such risks, are therefore essential.

New isolate of sweet potato virus 2 from Ipomoea nil: molecular characterization, codon usage bias, and phylogenetic analysis based on complete genome

BACKGROUND

Viral diseases of sweet potatoes are causing severe crop losses worldwide. More than 30 viruses have been identified to infect sweet potatoes among which the sweet potato latent virus (SPLV), sweet potato mild speckling virus (SPMSV), sweet potato virus G (SPVG) and sweet potato virus 2 (SPV2) have been recognized as distinct species of the genus Potyvirus in the family Potyviridae. The sweet potato virus 2 (SPV2) is a primary pathogen affecting sweet potato crops.

METHODS

In this study, we detected an SPV2 isolate (named SPV2-LN) in Ipomoea nil in China. The complete genomic sequence of SPV2-LN was obtained using sequencing of small RNAs, RT-PCR, and RACE amplification. The codon usage, phylogeny, recombination analysis and selective pressure analysis were assessed on the SPV2-LN genome.

RESULTS

The complete genome of SPV2-LN consisted of 10,606 nt (GenBank No. OR842902), encoding 3425 amino acids. There were 28 codons in the SPV2-LN genome with a relative synonymous codon usage (RSCU) value greater than 1, of which 21 end in A/U. Among the 12 proteins of SPV2, P3 and P3N-PIPO exhibited the highest variability in their amino acid sequences, while P1 was the most conserved, with an amino acid sequence identity of 87-95.3%. The phylogenetic analysis showed that 21 SPV2 isolates were clustered into four groups, and SPV2-LN was clustered together with isolate yu-17-47 (MK778808) in group IV. Recombination analysis indicated no major recombination sites in SPV2-LN. Selective pressure analysis showed dN/dS of the 12 proteins of SPV2 were less than 1, indicating that all were undergoing negative selection, except for P1N-PISPO.

CONCLUSION

This study identified a sweet potato virus, SPV2-LN, in Ipomoea nil. Sequence identities and genome analysis showed high similarity between our isolate and a Chinese isolate, yu-17-47, isolated from sweet potato. These results will provide a theoretical basis for understanding the genetic evolution and viral spread of SPV2.

The Discovery of a Citrus Yellow Vein Clearing Virus Hacienda Heights Isolate Diversifies the Geological Origins of the Virus in California, United States

The citrus yellow vein clearing virus (CYVCV) is an emerging threat to the U.S. citrus industry. Reports from China shows it cause significant reductions in fruit yield and growth, particularly in lemon trees. In 2022, CYVCV was detected in a wide range of citrus cultivars in localized urban properties in Tulare, California. In 2024, a CYVCV-infected lemon tree was detected in Hacienda Heights in Los Angeles County, California, geographically separated from the Tulare foci. Through long-read sequencing technology, the whole-genome sequence of a CYVCV isolate from Hacienda Heights (designated as CYVCV-CA-HH1, Accession number PP840891.1) was obtained. Sequence alignments and neighbornet analysis strongly suggested that the CYVCV-CA-HH1 isolate has a different origin than the Tulare CYVCV (CYVCV CA-TL) isolates. The CYVCV CA-TL isolates were grouped with those from South Asia (India and Pakistan) and the Middle East (Türkiye), while the CYVCV-CA-HH1 isolate was grouped with isolates from East Asia (China and South Korea). Maximum likelihood phylogenetic analysis further supports this finding, showing that the CYVCV-CA-HH1 isolate shares the most recent common ancestor with a South Korean lineage, which derives from Chinese isolates. Together, our data suggest a diverse geological origin of CYVCV isolates in California.

Independent repeated mutations within the alphaviruses Ross River virus and Barmah Forest virus indicates convergent evolution and past positive selection in ancestral populations despite ongoing purifying selection

Ross River virus (RRV) and Barmah Forest virus (BFV) are arthritogenic arthropod-borne viruses (arboviruses) that exhibit generalist host associations and share distributions in Australia and Papua New Guinea (PNG). Using stochastic mapping and discrete-trait phylogenetic analyses, we profiled the independent evolution of RRV and BFV signature mutations. Analysis of 186 RRV and 88 BFV genomes demonstrated their viral evolution trajectories have involved repeated selection of mutations, particularly in the nonstructural protein 1 (nsP1) and envelope 3 (E3) genes suggesting convergent evolution. Convergent mutations in the nsP1 genes of RRV (residues 248 and 441) and BFV (residues 297 and 447) may be involved with catalytic enzyme mechanisms and host membrane interactions during viral RNA replication and capping. Convergent E3 mutations (RRV site 59 and BFV site 57) may be associated with enzymatic furin activity and cleavage of E3 from protein precursors assisting viral maturation and infectivity. Given their requirement to replicate in disparate insect and vertebrate hosts, convergent evolution in RRV and BFV may represent a dynamic link between their requirement to selectively 'fine-tune' intracellular host interactions and viral replicative enzymatic processes. Despite evidence of evolutionary convergence, selection pressure analyses did not reveal any RRV or BFV amino acid sites under strong positive selection and only weak positive selection for nonstructural protein sites. These findings may indicate that their alphavirus ancestors were subject to positive selection events which predisposed ongoing pervasive convergent evolution, and this largely supports continued purifying selection in RRV and BFV populations during their replication in mosquito and vertebrate hosts.

Wolbachia of phylogenetic supergroup K identified in oribatid mite Nothrus anauniensis (Acari: Oribatida: Nothridae)

Heritable endosymbionts widely occur in arthropod and nematode hosts. Among these endosymbionts, Wolbachia has been extensively detected in many arthropods, such as insects and crustaceans. Maternal inheritance is the most basic and dominant mode of transmission of Wolbachia, and it might regulate the reproductive system of the host in four ways: feminization, parthenogenesis, male killing, and cytoplasmic incompatibility. There is a relatively high percentage (10%) of thelytokous species in Oribatida, a suborder under the subclass Acari of arthropods, but the study of the endosymbionts in oribatid mites is almost negligible. In this paper, we detected endosymbiotic bacteria in two parthenogenetic oribatid species, Nothrus anauniensis Canestrini and Fanzago, 1877, which has never been tested for endosymbionts, and Oppiella nova, in which Wolbachia and Cardinium have been reported before. The results showed that Wolbachia was first found in N. anauniensis with an infection rate of 100% across three populations. Phylogenetic analysis showed that Wolbachia in N. anauniensis belonged to the supergroup K, marking the second supergroup of Wolbachia found in oribatid mites. Unlike previous studies, our study did not detect Wolbachia in O. nova, leading to the exclusion of Wolbachia's role in mediating thelytoky in this species.

Bayesian phylogeographic inference of wheat yellow mosaic virus in China and Japan suggests that the virus migration history coincided with historical events

Wheat yellow mosaic virus (WYMV) is one of the most serious viral pathogens causing reductions in wheat yield in East Asia. We investigated the phylodynamics of WYMV by analysing the CP, VPg and P1 genes to understand the origin and dispersal of the virus. A Bayesian phylogenetic analysis revealed that the most recent common WYMV ancestor occurred in approximately 1742 (95% credibility interval, 1439-1916) CE (Common Era), and the evolutionary rates of the VPg, CP and P1 genes were 6.669 × 10-4 (95% credibility interval: 4.575 × 10-4-8.927 × 10-4), 2.468 × 10-4 (95% credibility interval: 1.667 × 10-4-3.338 × 10-4) and 5.765 × 10-5 (95% credibility interval: 3.285 × 10-6-1.252 × 10-4), respectively. Our phylogeographic analysis indicated that the WYMV population may have originated in Henan Province, China, first spreading to Japan in the mid-19th century and stopping after the Japanese surrender in World War II. The second wave spread to Japan from Shandong Province, China, in approximately 1977, a few years after the establishment of diplomatic relations between China and Japan. Before the founding of the People's Republic of China, Henan Province was the emigration centre of WYMV in East Asia, and after the late 20th century, Jiangsu and Shandong Provinces were also the virus emigration centres in East Asia. In addition, there were two migration pathways from Japan to Jiangsu and Shandong Provinces, China, in approximately 1918 and approximately 1999 respectively. Our results suggest that the wide spread of WYMV in East Asia is strongly related to human factors.

The evolutionary and transmission dynamics of HIV-1 CRF08_BC

HIV-1 CRF08_BC is a significant subtype in China, though its origin and spread remain incompletely understood. Previous studies using partial genomic data have provided insights but lack comprehensive analysis. Here, we investigate the early evolutionary and spatiotemporal dynamics of HIV-1 CRF08_BC in China and Myanmar using near-complete genome sequences. We analyzed 28 near-complete HIV-1 CRF08_BC genomes from China and Myanmar (1997-2013). Phylogenetic, molecular clock, and Bayesian discrete trait analyses were performed to infer the virus's origin, spread, and associated risk groups. Based on Bayesian time-scaled inference with the best-fitting combination of models determined by marginal likelihood estimation (MLE), we inferred the time to the most recent common ancestor (TMRCA) and evolutionary rate of HIV-1 CRF08_BC to be at 3 October 1991 (95% HPD: 22 February1989-27 November 1993) and 2.30 × 10-3 substitutions per site per year (95% HPD: 1.96 × 10-3-2.63 × 10-3), respectively. Our analysis suggests that HIV-1 CRF08_BC originated in Yunnan Province, China, among injecting drug users, and subsequently spread to other regions. This study provides valuable insights into the early dynamics of HIV-1 CRF08_BC through combined genomic and epidemiological data, which may inform effective prevention and mitigation efforts. However, the limited genomic data influenced the extent of our findings, and challenges in collecting accurate risk group information during surveillance were evident.

The classification, origin, and evolutionary dynamics of severe fever with thrombocytopenia syndrome virus circulating in East Asia

The classification of severe fever with thrombocytopenia syndrome virus (SFTSV) lacked consistency due to limited virus sequences used across previous studies, and the origin and transmission dynamics of the SFTSV remains not fully understood. In this study, we analyzed the diversity and phylodynamics of SFTSV using the most comprehensive and largest dataset publicly available for a better understanding of SFTSV classification and transmission. A total of 1267 L segments, 1289 M segments, and 1438 S segments collected from China, South Korea, and Japan were included in this study. Maximum likelihood trees were reconstructed to classify the lineages. Discrete phylogeographic analysis was conducted to infer the phylodynamics of SFTSV. We found that the L, M, and S segments were highly conserved, with mean pairwise nucleotide distances of 2.80, 3.36, and 3.35% and could be separated into 16, 13, and 15 lineages, respectively. The evolutionary rate for L, M, and the S segment was 0.61 × 10-4 (95% HPD: 0.48-0.73 × 10-4), 1.31 × 10-4 (95% HPD: 0.77-1.77 × 10-4) and 1.27 × 10-4 (95% HPD: 0.65-1.85 × 10-4) subs/site/year. The SFTSV most likely originated from South Korea around the year of 1617.6 (95% HPD: 1513.1-1724.3), 1700.4 (95% HPD: 1493.7-1814.0), and 1790.1 (95% HPD: 1605.4-1887.2) for L, M, and S segments, respectively. Hubei Province in China played a critical role in the geographical expansion of the SFTSV. The effective population size of SFTSV peaked around 2010 to 2013. We also identified several codons under positive selection in the RdRp, Gn-Gc, and NS genes. By leveraging the largest dataset of SFTSV, our analysis could provide new insights into the evolution and dispersal of SFTSV, which may be beneficial for the control and prevention of severe fever with thrombocytopenia syndrome.

Movement of the A-strain maize streak virus in and out of Madagascar

The maize streak virus belongs in the genus Mastrevirus, in the family Geminiviridae. The A-strain of the virus (MSV-A) is recognised as the principal causative agent of the most severe manifestation of maize streak disease (MSD). This disease continues to be a persistent limitation on maize output across sub-Saharan Africa and the nearby Indian Ocean islands. Irrespective of the causes behind the spread of MSV-A, we can determine the paths and speeds with which MSV-A spreads by analysing MSV genome sequence data along with information on when and where samples were taken. This information is valuable for identifying the geographical origins of viral strains that cause sporadic MSD epidemics in specific places and the geographical regions where viruses remain in reservoirs and contribute to prolonged epidemics during outbreaks. Our aim is to utilise these analyses to estimate the timing and origin of the MSV-A that arrived on the island of Madagascar in the Indian Ocean. Specifically, we employ model-based phylogeographic analyses on 524 complete MSV-A genome sequences, which consist of 56 newly obtained genomes from infected maize plants collected in Madagascar. These studies allow us to reconstruct the most likely paths of MSV-A to Madagascar. We found strong evidence for the existence of at least four separate movements of MSV-A variants from East and southern Africa to Madagascar. These movements took place between roughly 1979 (with a 95% highest probability density interval [HPD] ranging from 1976 to 1982) and 2003 (with a 95% HPD ranging from 2002 to 2003). While we inferred that MSV-A variants are spreading at an average rate of 38.9 km/year (with a 95% highest posterior density interval of 34.0-44.4) across their geographical range. Since their arrival in Madagascar, MSV-A variants have been migrating at an average rate of 47.6 km/year (with a 95% highest posterior density interval of 36.05-61.70). Human influences are likely significant contributors to both sporadic long-range movements of MSV-A between mainland Africa and Madagascar, as well as shorter to medium range movements within the island.

Molecular characterization to study the genetic diversity of begomoviruses occurring in the major chilli growing areas of Tamil Nadu state of India

Chilli leaf curl disease (ChiLCD), which is a significant problem in chilli cultivation, is caused by begomoviruses that are transmitted by the whitefly Bemisia tabaci. This disease leads to severe impacts on crop yields. To determine the incidence of begomovirus in the chilli crop, infected chilli leaf samples exhibiting symptoms such as curling, yellowing, reduced leaf size, and overall stunted growth were collected from various districts of Tamil Nadu, namely, Coimbatore, Dharmapuri, Kancheepuram, Karur, Salem, Krishnagiri, Thoothukudi, Thiruvallur, Tiruchirappalli, Virudhunagar, Tiruvannamalai, Tenkasi, and Vellore, during the years 2018-2022. To determine the complete genome sequence of the begomoviruses, the rolling circle amplification (RCA) method was used to clone and sequence the begomovirus genomes from the chilli samples collected from various districts of Tamil Nadu. Here we characterized 17 DNA A genome sequences and 12 betasatellite sequences. BLAST results of the DNA A genome sequences revealed nucleotide identities ranging from 94.2 to 99.7% with five distinct begomovirus species of chilli, namely, chilli leaf curl Salem virus (HM007119), chilli leaf curl virus Bhavanisagar (NC_055130), chilli leaf curl Ahmedabad virus (MW795666), chilli leaf curl virus (NC_055131), and chilli leaf curl Sri Lanka virus (JN555600). BLAST results of the betasatellite sequences showed nucleotide identities of 96 to 98.8% with the tomato leaf curl Bangladesh betasatellite (MZ151286). In the present study, five distinct begomovirus species and one associated betasatellite were found to infect chilli crops in Tamil Nadu. This finding indicates a changing pattern of begomovirus occurrence in the different districts of Tamil Nadu. This study highlights the prevalence of chilli-infecting begomoviruses in the major chilli growing districts of Tamil Nadu, the identification of begomovirus species, and the significance of understanding and managing these viruses to safeguard chilli cultivation in the region.

Genetic Diversity and Pathogenicity Characterization of Tomato-Infecting Begomoviruses in Taiwan

The tomato yellow leaf curl disease (TYLCD) caused by whitefly (Bemisia tabaci)-transmitted begomoviruses (Geminiviridae) has constrained tomato production in Taiwan since 1981. Lisianthus enation leaf curl virus (LELCV), tomato leaf curl Taiwan virus (ToLCTV), and tomato yellow leaf curl Thailand virus (TYLCTHV) were the major viruses associated with TYLCD. In 2019 to 2020, we investigated TYLCD throughout Taiwan, with a 10 to 100% incidence on tomato fields. Begomovirus sequences were detected in 321 out of 506 collected samples by PCR with primers PAL1v1978B and PAR1c715H. In 2015 to 2016, 59 out of 99 samples collected in Hualien-Taitung areas were also found to have begomovirus sequences. Based on the analysis of 68 viral genomic sequences, six begomoviruses were identified, including LELCV, ToLCTV, TYLCTHV, tomato leaf curl Hsinchu virus, and two new begomoviruses, tentatively named tomato leaf curl Chiayi virus (ToLCCYV) and tomato leaf curl Nantou virus (ToLCNTV). Various isolates of LELCV and TYLCTHV were grouped into four and two strains, respectively. Recombinants were detected in LELCV-A, -C, and -D, ToLCCYV, ToLCNTV, and TYLCTHV-F. Based on virus-specific detection, the majority of TYLCD-associated viruses were mixed-infected by TYLCTHV-B with TYLCTHV-F, LELCV-A, -B, or -D, and/or ToLCTV. Meanwhile, viral DNA-B was mostly associated with TYLCTHV, and all identified DNA-Bs were highly homologous with previous TYLCTHV DNA-B. The pathogenicity of selected begomoviruses was confirmed through agroinfection and whitefly transmission. All tomato plants carrying Ty-1/3 and Ty-2 resistant genes were infected by all LELCV strains and ToLCCYV, although they appeared symptomless, suggesting these viruses could be managed through the use of the resistance pyramid.

Identification of a Novel Parvovirus in the Arctic Wolf (Canis lupus arctos)

A novel virus, temporarily named "Arctic wolf parvovirus" (AWPV), was discovered in a pharyngeal metagenomic library derived from an Arctic wolf (Canis lupus arctos) in China. The genome sequence was assigned GenBase accession number C_AA071902.1. AWPV has a genome comprised of 4,920 base pairs with a nucleotide composition of 36.4% A, 23.4% T, 18.2% G, and 22.0% C, with a GC content of 40.2%. Its structure resembles parvoviruses, containing two open reading frames: the nonstructural (NS) region encoding replication enzymes and the structural (VP) region encoding capsid protein. Pairwise sequence comparison and phylogenetic analysis suggest AWPV may represent a novel species within the genus Protoparvovirus. This discovery enhances our understanding of mammalian virus ecology and potential future infectious diseases.

Bacillus cereus

Bacillus cereus is an emerging food-borne pathogen responsible for two types of food poisoning: emetic and diarrhoeal type. Here we report an emetic type of food-borne illness attributable to Bacillus cereus. On 2nd February, 2021, 202 people suffered from pain in abdomen and vomiting after consuming the rice provided during a public gathering in Diphu, Assam. Culture of leftover fried rice showed growth of Bacillus cereus group of organisms. Molecular detection of enterotoxin and emetic toxin genes was done in the isolated strains by polymerase chain reaction. Multi locus sequence typing (MLST) and phylogenetic analysis was done to characterise the isolated strains. A total of five strains of Bacillus cereus were isolated. Ces gene was found in isolates GMC22 & GMC24 and other enterotoxins producing genes were found in isolates GMC23 and GMC24. MLST identified four sequence types (STs) (ST1051, ST1616, ST998 and ST1000). Phylogenetic analysis clustered ST-1051 assigned to the GMC22 strain into the previously defined clade I and was in close relation with ST-144, representing a new cereulide-producing emetic cluster. As Bacillus cereus is a common contaminant of foods, it is essential to evaluate the pathogenic potential of the bacteria for a definite link between causative agents and the illness. MLST can be used to characterize the Bacillus cereus strains isolated from outbreak samples in order to link the probable pathogens with the illness. In this outbreak, we suggest that ST-1051 is the strain responsible for the food-borne illness, which was predominantly of emetic type.

Genetic diversity, variation and recombination among the human papillomaviruses (HPVs) genomes isolated in China: a comparative genomic and phylogenetic analysis

Human papillomaviruses (HPVs) are widespread, sexually transmitted group of viruses that infect most individuals at some stage, causing genital warts and cancers. They are members of the Papillomaviridae family, which contains about 400 HPV types. China is among the high HPV burden countries with reported infections of multiple HPV types, accounting for 17.3% of global deaths and 18.2% of global new cases. Thus, understanding the genetic variation and geographic diversity characteristics of HPVs isolated in China is critical for global HPV prevention strategies. Thus, we analyzed the available HPV genome sequences isolated in China that grouped into two categories (alpha- and gamma-papillomaviruses) based on full-length genomes. The most common were HPV-16, -6, -58, and -52 respectively. In addition, four of the novel strains isolated in China, e.g. TG550, JDFY01, CH2, and L55 clustered with the HPV-mSK 159, 244, 201, and 200 respectively. Our phylogeographic network analysis indicated that the L55, TG550, and CH2 are genetically identical to the mSK 200, 046, and 201 respectively, while JDFY01 appeared separately, connected to the mSK-040 following five mutational steps. Also, we found ten recombination events among HPV-6/11 types within their E1, E2, E7, L1/L2 proteins, and Long Control Region ORFs. We achieved the consensus amino acid sequences of HPV proteins and found a conserved stretch of amino acids within E5A of all HPVs circulating in China. These findings offer valued insights into the genetic relationships, distribution, and evolution of the HPVs in China that may assist in adapting effective HPV preventive measures.

Molecular Epidemiology of Western Equine Encephalitis Virus, South America, 2023-2024

Western equine encephalitis virus (WEEV) is a mosquitoborne virus that reemerged in December 2023 in Argentina and Uruguay, causing a major outbreak. We investigated the outbreak using epidemiologic, entomological, and genomic analyses, focusing on WEEV circulation near the Argentina‒Uruguay border in Rio Grande do Sul state, Brazil. During November 2023‒April 2024, the outbreak in Argentina and Uruguay resulted in 217 human cases, 12 of which were fatal, and 2,548 equine cases. We determined cases on the basis of laboratory and clinical epidemiologic criteria. We characterized 3 fatal equine cases caused by a novel WEEV lineage identified through a nearly complete coding sequence analysis, which we propose as lineage C. Our findings highlight the importance of continued surveillance and equine vaccination to control future WEEV outbreaks in South America.

Comparative genomics of canine parvovirus in South America: Diversification patterns in local populations

Canine parvovirus (CPV) is a significant pathogen in domestic dogs worldwide, causing a severe and often fatal disease. CPV comprises three antigenic variants (2a, 2b, and 2c) distributed unevenly among several phylogenetic groups. The present study compared genetic variability and evolutionary patterns in South American CPV populations. We collected samples from puppies suspected of CPV infection in the neighboring Argentina and Uruguay. Antigenic variants were preliminarily characterized using PCR-RFLP and partial vp2 sequencing. Samples collected in Argentina during 2008-2018 were mainly of the 2c variant. In the Uruguayan strains (2012-2019), the 2a variant wholly replaced the 2c from 2014. Full-length coding genome and vp2 sequences were compared with global strains. The 2c and 2a strains fell by phylogenetic analysis into two phylogroups (Europe I and Asia I). The 2c strains from Argentina and Uruguay clustered in the Europe I group, with strains from America, Europe, Asia, and Oceania. Europe I is widely distributed in South America in the dog population and is also being detected in the wildlife population. The 2a strains from Uruguay formed the distinct Asia I group with strains from Asia, Africa, America, and Oceania. This Asia I group is increasing its distribution in South America and worldwide. Our research reveals high genetic variability in adjacent synchronic samples and different evolutionary patterns in South American CPV. We also highlight the importance of ancestral migrations and local diversification in the evolution of global CPV strains.

Detection and Molecular Characterization of Porcine Teschoviruses in India: Identification of New Genotypes

Porcine Teschoviruses (PTVs) are ubiquitous enteric viral pathogens that infect pigs and wild boars worldwide. PTVs have been responsible for causing the severe clinical disease (Teschen disease) to asymptomatic infections. However, to date, limited information is available on large-scale epidemiological data and molecular characterization of PTVs in several countries. In this study, we report epidemiological data on PTVs based on screening of 534 porcine fecal samples from different states of India and a RT-PCR based detection of PTVs shows a percent positivity of 8.24% (44/534). The PTV prevalence varied among different regions of the country with the highest detection rates observed in the state of Karnataka (38.1%). Phylogenetic analysis based on VP1 gene reveals the presence of PTV genotype 6 and 13 along with some unassigned novel genotypes which did not cluster with any of the established PTV genotypes (PTV 1-PTV 13). Indian PTV 6 strains are genetically closest to the Spanish strains (85.7-94.4%) whereas PTV 13 and novel genotype strains were found to be more similar to the Chinese strains (88.1-99.1%). Using recombination detection software, no Indian PTVs found to be recombinant on VP1 gene and selection pressure analysis revealed the purifying selection in the several sites of the VP1 gene of PTVs. The Bayesian analysis of Indian PTVs shows 1.16 × 10-4 substitution/site/year as the mean evolutionary rate. Further, isolation of the novel PTV strains from India and more detailed investigation much needed to know the evolutionary history of PTV strains circulating in porcine populations in India.

Genetic substructure and host-specific natural selection trend across vaccine-candidate ORF-2 capsid protein of hepatitis-E virus

Hepatitis E virus is a primary cause of acute hepatitis worldwide. The present study attempts to assess the genetic variability and evolutionary divergence among HEV genotypes. A vaccine promising capsid-protein coding ORF-2 gene sequences of HEV was evaluated using phylogenetics, model-based population genetic methods and principal component analysis. The analyses unveiled nine distinct clusters as subpopulations for six HEV genotypes. HEV-3 genotype samples stratified into four different subgroups, while HEV-4 stratified into three additional subclusters. Rabbit-infectious HEV-3ra samples constitute a distinct cluster. Pairwise analysis identified marked genetic distinction of HEV-4c and HEV-4i subgenotypes compared to other genotypes. Numerous admixed, inter and intragenotype recombinant strains were detected. The MEME method identified several ORF-2 codon sites under positive selection. Some selection signatures lead to amino acid substitutions within ORF-2, resulting in altered physicochemical features. Moreover, a pattern of host-specific adaptive signatures was identified among HEV genotypes. The analyses conclusively depict that recombination and episodic positive selection events have shaped the observed genetic diversity among different HEV genotypes. The significant genetic diversity and stratification of HEV-3 and HEV-4 genotypes into subgroups, as identified in the current study, are noteworthy and may have implications for the efficacy of anti-HEV vaccines.

Epidemics and diversity of norovirus variants with acute gastroenteritis outbreak in Hongshan District, Wuhan City, China, 2021-2023

BACKGROUND

Norovirus is the predominant pathogen causing foodborne illnesses and acute gastroenteritis (AGE) outbreaks worldwide, imposing a significant disease burden. This study aimed to investigate the epidemiological characteristics and genotypic diversity of norovirus outbreaks in Hongshan District, Wuhan City.

METHODS

A total of 463 AGE cases from 39 AGE-related outbreaks in Hongshan District between January 1, 2021, and June 30, 2023, were included in the study. Reverse transcription-polymerase chain reaction (RT-PCR) was used to identify norovirus types GI and GII in anal swab samples from all cases. Norovirus-positive samples were sequenced and analyzed for the open reading frame (ORF) 1/ORF2 hinge region.

RESULTS

26 norovirus infectious outbreaks were reported among 39 acute diarrheal outbreaks, including 14 outbreaks in kindergartens, 8 in elementary schools, and 4 in universities. Based on clinical symptoms and epidemiological investigations, a total of 1295 individuals were identified as having been exposed to norovirus, yielding an attack rate of 35.75 %. A higher proportion of outbreaks was observed during the winter and spring seasons (38.46 %). Additionally, norovirus-positive samples were subjected to sequencing and analysis of the open reading frame (ORF) 1/ORF2 hinge region. Genotypic data for norovirus was successfully obtained from 18 (69.23 %) of the infectious outbreaks, revealing 10 distinct recombinant genotypes. GII.4 Sydney 2012 [P31] and GII.17[P17] were the predominant strains in 2021 and 2022, GII.3 [P12] emerged as the dominant strain in 2023.

CONCLUSION

Norovirus outbreaks in Hongshan District predominantly occurred in crowded educational institutions, with peaks in the cold season and a high attack rate in universities. GII.3 [P12] has become the locally predominant strain.

Codon Usage Analysis Reveals Distinct Evolutionary Patterns and Host Adaptation Strategies in Duck Hepatitis Virus 1 (DHV-1) Phylogroups

Duck hepatitis virus 1 (DHV-1) is a major threat to the global poultry industry, causing significant economic losses due to high mortality rates in young ducklings. To better understand the evolution and host adaptation strategies of DHV-1, we conducted a comprehensive codon usage analysis of DHV-1 genomes. Our phylogenetic analysis revealed three well-supported DHV-1 phylogroups (Ia, Ib, and II) with distinct genetic diversity patterns. Comparative analyses of the codon usage bias and dinucleotide abundance uncovered a strong preference for A/U-ended codons and a biased pattern of dinucleotide usage in the DHV-1 genome, with CG dinucleotides being extremely underrepresented. Effective number of codons (ENC) analysis indicated a low codon usage bias in the DHV-1 ORF sequences, suggesting adaptation to host codon usage preferences. PR2 bias, ENC plot, and neutrality analyses revealed that both mutation pressure and natural selection influence the codon usage patterns of DHV-1. Notably, the three DHV-1 phylogroups exhibited distinct evolutionary trends, with phylogroups Ia and Ib showing evidence of neutral evolution accompanied by selective pressure, while the phylogroup II evolution was primarily driven by random genetic drift. Comparative analysis of the codon usage indices (CAI, RCDI, and SiD) among the phylogroups highlighted significant differences between subgroups Ia and Ib, suggesting distinct evolutionary pressures or adaptations influencing their codon usage. These findings contribute to our understanding of DHV-1 evolution and host adaptation, with potential implications for the development of effective control measures and vaccines.

Genome sequencing of canine distemper virus isolates from unvaccinated dogs in Mongolia

Canine distemper virus (CDV) triggers a severe, often fatal disease in dogs and wildlife known as canine distemper (CD). Prior research has noted significant genetic diversity and recombination among CDV isolates from different geographical regions, potentially contributing to vaccine failures. Despite this, no genetic characterization of Mongolian CDVs has been conducted. This study, isolated CDVs from three unvaccinated dogs: two 10-month-old mixed-breeds and an 18-month-old Samoyed. All exhibited CD symptoms and subsequently died. Virus isolation was conducted using Vero/dog SLAM cells, with genome sequencing performed via nanopore technology. The mixed-breed dogs were infected with non-recombinant CDV isolates, forming a sister clade to the Asia-1 lineage prevalent in Asia. The Samoyed was infected with a non-recombinant CDV isolate, classifying as Asia-4 lineage sporadically reported in some Asian countries. This sequencing data offers foundational information on genetic diversity, aiding CD control measure development and benefiting future Eurasia and Asian studies.

Molecular characteristic, evolution, and pathogenicity analysis of avian infectious bronchitis virus isolates associated with QX type in China

Infectious bronchitis virus (IBV) is one of the major avian pathogens plaguing the global poultry industry. Although vaccination is the primary preventive measure for IBV infection, the emergence of virus variants with mutations and recombination has resulted in IBV circulating globally, presenting a challenge for IB control. Here, we isolated 3 IBV strains (CZ200515, CZ210840, and CZ211063) from suspected sick chickens vaccinated with IBV live attenuated vaccines (H120, 4/91, or QXL87). Phylogenetic analysis of the S1 gene sequence of the spike (S) revealed that the 3 isolates belonged to the QX-type (GI-19 lineage). Whole genome sequencing and recombination analysis indicated that CZ200515 and CZ210840 contained genetic material from 4/91 and Scyz3 (QX-type), possibly due to recombination between the circulating strain and the 4/91 vaccine strain, while no evidence of recombination was found in CZ211063. Pathogenicity analysis in 1-day-old specific pathogen-free (SPF) chickens demonstrated that all 3 isolates caused severe tissue damage and varying degrees of mortality. Virus cross-neutralization assay revealed decreased antigen relatedness between the isolates and the QX-type vaccine strain (QXL87). Amino acid sequence homology analysis of S1 revealed 5%-6.5% variances between the isolates and QXL87. Analysis of the S1 subunit structure revealed that mutations of amino acid residues in the hypervariable region (HVR) and the neutralizing epitope region resulted in antigenic variation in isolates by changing the antigen conformation. Our data indicate antigenicity variances between QX isolates and QXL87 vaccine strains, potentially resulting in immune evasion occurrences. Overall, these results offer crucial insights into the epidemiology and pathogenicity of QX-type IBV, facilitating improved selection and formulation of vaccines for disease management.

The tissue virome of black-spotted frogs reveals a diversity of uncharacterized viruses

Amphibians are an essential class in the maintenance of global ecosystem equilibrium, but they face serious extinction risks driven by climate change and infectious diseases. Unfortunately, the virus diversity harbored by these creatures has been rarely investigated. By profiling the virus flora residing in different tissues of 100 farmed black-spotted frogs (Rana nigromaculata) using a combination of DNA and RNA viromic methods, we captured 28 high-quality viral sequences covering at least 11 viral families. Most of these sequences were remarkably divergent, adding at least 10 new species and 4 new genera within the families Orthomyxoviridae, Adenoviridae, Nodaviridae, Phenuiviridae, and Picornaviridae. We recovered five orthomyxovirus segments, with three distantly neighboring two Chinese fish-related viruses. The recombination event of frog virus 3 occurred among the frog and turtle strains. The relative abundance and molecular detection revealed different tissue tropisms of these viruses, with the orthomyxovirus and adenoviruses being enteric and probably also neurotropic, but the new astrovirus and picornavirus being hepatophilic. These results expand the spectrum of viruses harbored by anurans, highlighting the necessity to continuously monitor these viruses and to investigate the virus diversity in a broader area with more diverse amphibian species.

High-throughput sequencing detected a virus-viroid complex in a single pokeweed plant

In this study, total RNA high-throughput sequencing (HTS) of a single symptomatic Phytolacca americana plant enabled the obtention of a nearly complete genome of two new isolates of turnip yellows virus (TuYV), named TuYV-ITA1 and TuYV-ITA2, and revealed a mixed infection with a new variant of citrus exocortis viroid (CEVd), named CEVd-ITA1. The TuYV-ITA2 isolate diverged from the known virus isolates of TuYV and showed variability in the P0 and P5 readthrough domain. Recombination analysis revealed its recombinant nature between TuYV and an unidentified polerovirus. The putative recombination event was identified in the P5 readthrough domain of the TuYMV-ITA2 isolate. Our results thus represent the first report of TuYV in Italy and some molecular evidence for the possible natural co-existence of TuYV and CEVd in a new natural host for both infectious entities. This study is adding further knowledge about the role of weed plants as virus reservoirs, and thus additional biological and impact studies would be desirable to determine in particular the role of P. americana in the spread of TuYV and if this virus should be considered a new threat for the susceptible Italian crops.

First reports of several viruses and a viroid including a novel vitivirus in Japan, found through virome analysis of bulk grape genetic resources

Virome analysis was performed on 174 grape genetic resources from the National Agriculture and Food Research Organization, Japan. A total of 20 bulk samples was prepared by grouping the vines into batches of 6-10 plants. Each of the bulk samples was analyzed using high-throughput sequencing, which detected 27 viruses and 5 viroids, including six viruses and one viroid reported in Japan for the first time (grapevine viruses F, L, and T, grapevine Kizil Sapak virus, grapevine Syrah virus 1, grapevine satellite virus, and grapevine yellow speckle viroid 2). In addition, a novel vitivirus was detected with a maximum nucleotide sequence identity of only 58% to its closest relative, grapevine virus A (GVA). The genome of this novel virus was 7,461 nucleotides in length and encoded five open reading frames showing the typical genomic structure of vitiviruses. Phylogenetic trees of vitiviruses placed it in a distinct position nearest to GVA or grapevine virus F (GVF) in genomes and amino acids of deduced replication-associated protein (RAP) and coat protein (CP). The amino acid sequence identities of RAP and CP with GVA, GVF, and other vitiviruses were a maximum of 53% and 73%, respectively, which were significantly below the species demarcation threshold of 80% in the genus. The low identity and phylogenetic analyses indicate the discovery of a novel vitivirus species provisionally named grapevine virus P.