An exact nonparametric method for inferring mosaic structure in sequence triplets

Epidemiology, genetic diversity and evolution of pigeon circovirus

With the global prevalence of pigeon racing, pigeon infections with Pigeon circovirus (PiCV) have been reported across the world. PiCV-infected pigeons are susceptible to subsequent infections and exhibit a range of symptoms including weight loss, lethargy, anorexia, respiratory distress and diarrhea, which exerts a huge impact on the pigeon industry. However, there is currently a lack of epidemiological data on PiCV infection in pigeons. Therefore, the aim of this study was to gain a comprehensive understanding of the prevalence, genetic variations, and evolutionary dynamics of PiCV within pigeon populations. In this study, we collected 28 samples from four cities in China to assess the prevalence of PiCV. The results showed that the positive rate was 92.86%, which indicated that PiCV was widespread in Chinese pigeon flocks. Meanwhile, whole genome sequences of PiCV were obtained for partial positive samples. The results showed that there were 12 different clades of PiCV, and the samples collected in this study have been classified as types 1, 4, 6, and 11. The results of the PiCV genome analysis indicate a high incidence of recombination events. Further analysis of the encoded Cap protein at the gene and protein level revealed a high degree of variation in the Cap protein. The amino acids at positions 30-120 of the Cap protein exhibited the greatest frequency of variation. Meanwhile, we found that amino acid sequences 140-180 were relatively conserved, which was related to the immunogenicity of the circovirus Cap protein. The results of the antigenicity analysis demonstrated that amino acid sequences 140-180 exhibited strong antigenicity, indicating the potential of cap to serve as a candidate protein in the production of a PiCV vaccine. The present study provided epidemiological information on PiCV, which promises to facilitate the development of effective control measures to prevent its transmission.

Evolutionary analysis and population dynamics in the global transmission of Kaposi's sarcoma-associated herpesvirus

With the widespread prevalence of Kaposi's sarcoma-associated herpesvirus (KSHV) globally, particularly in sub-Saharan Africa and the Mediterranean region, the open reading frame (ORF) K1 gene, a key gene for distinguishing different subtypes of KSHV, has been extensively studied for its diversity and sequence variations. In this study, we collected K1 gene sequences representing subtypes of KSHV worldwide in order to assess the global distribution of KSHV subtypes and to investigate the recombination and selection history of KSHV. Recombination and gene flow analysis indicated a minimum average recombination rate of 0.41 per site for the K1 gene. Recombination analysis indicated that 11 major recombination events had occurred, predominantly in subtypes A and C, while subtype B showed minimal involvement in recombination processes, consistent with the gene flow analysis. Using tip-dating methods, we estimated that the most recent common ancestor of KSHV emerged in the 12th century, while the currently globally prevalent subtypes appeared within the past three centuries. Its recent origin and rapid evolution indicate that KSHV is now undergoing strong selection and is in the process of adaptation.

Genomic variants and molecular epidemiological characteristics of dengue virus in China revealed by genome-wide analysis

Since its first academic record in 1978, dengue epidemics have occurred in all provinces of China, except Xizang. The epidemiological and molecular features of the whole genome of dengue virus (DENV) have not yet been completely elucidated, interfering with prevention and control strategies for dengue fever in China. Here, we obtained 553 complete genomes of the four serotypes of DENV (DENV1-4) isolated in China from the GenBank database to analyze the phylogeny, recombination, genomic variants, and selection pressure and to estimate the substitution rates of DENV genomes. Phylogenetic analyses indicated that DENV sequences from China did not cluster together and were genetically closer to those from Southeast Asian countries in the maximum likelihood trees, indicating that DENV was not endemic in China. Thirty intra-serotype recombinant sequences were identified for DENV1-4, with the highest frequency in DENV4. Selection pressure analyses revealed that 13 codons under positive selection were located in the C, NS1, NS2A, NS3, and NS5 proteins. For DENV1 to DENV3, the substitution rates evaluated in this study were 9.23 × 10-4, 7.59 × 10-4, and 7.06 × 10-4 substitutions per site per year, respectively. These findings improve our understanding of the evolution of DENV in China.

Avian Migration-Mediated Transmission and Recombination Driving the Diversity of Gammacoronaviruses and Deltacoronaviruses

In the wake of pandemics like COVID-19, which have zoonotic origins, the role of wildlife as reservoirs for emerging infectious diseases has garnered heightened attention. Migratory birds, traversing continents, represent a potent but under-researched vector for the spread of infectious diseases, including novel coronaviruses. This study delves into the genetic diversity and transmission dynamics of coronaviruses in migratory birds, presenting pivotal findings. From April 2019 to April 2023, we screened 5,263 migratory bird samples collected from Shanghai, China, identifying 372 coronavirus-positive samples belonging to five avian-related coronavirus subgenera and subsequently obtaining 120 complete genome sequences. To facilitate further research with a global perspective, the study curated all available 19,000 avian-associated coronaviruses and expanded the original 12 species to 16, including three novel coronavirus species identified in our study and one re-classified species from the public domain. The study illuminates the intricate genetic evolution and transmission dynamics of birds-related coronaviruses on a global scale. A notable aspect of our research is the identification of complex recombination patterns within the spike protein across different virus species and subgenera, highlighting migratory birds as a reservoir of coronavirus. Notably, the coronaviruses found in migratory birds, predominantly from the orders Anseriformes, Charadriiformes, and Pelecaniformes, with domestic ducks from Anseriformes playing a key role in bridging the transmission of coronaviruses between migratory and non-migratory birds. These findings reveal the genetic and recombination characteristics of coronaviruses in migratory birds, emphasizing the critical role of ecologically pivotal bird species in coronavirus transmission and genetic diversity shaping.

Genomic surveillance and evolution of Getah virus

Getah virus (GETV), a member of the Alphaviruses, has spread widely and is expanding its host range worldwide, posing a serious threat to public health safety and the farming industry. However, genetic monitoring of GETV is inadequate, and its evolution and transmission remain unclear. This study employed reverse transcription-polymerase chain reaction to screen pig tissue samples for the presence of GETV. Subsequent steps included DNA sequencing, phylogenetic analysis, and selection pressure assessments to elucidate the evolutionary history and transmission patterns of the virus. A total of 1382 samples were examined, with a positivity rate of 4.12% (95% confidence interval: 3.07%-5.17%) from 2022 to 2023. Subsequently, seven GETV strains were isolated and identified. A phylogenetic tree was constructed, which showed that all seven strains belonged to Group III. Phylodynamic analysis revealed that GETV evolved rapidly. Additionally, eight amino acid sites within the GETV E2 protein were identified as being under positive selection. These data provide insight into the epidemiology and evolution of GETV.

Mamastrovirus species are shaped by recombination and can be reliably distinguished in ORF1b genome region

Astroviruses are a diverse group of small non-enveloped positive sense single-stranded RNA viruses that infect animals and birds. More than half of all known genome sequences of mammalian astroviruses are not assigned to provisional species, and the biological mechanisms that could support segregation of astroviruses into species are not well understood. The systematic analysis of recombination in Mamastrovirus genomes available in GenBank was done to identify mechanisms providing genetic distinction between astroviruses. Recombination breakpoints were present in all Mamastrovirus genome regions, but occurred most commonly at the ORF1b/ORF2 junction. Recombination was ubiquitous within, but never between established and putative new species, and may be suggested as an additional species criterion. The current species criterion for the genus Mamastrovirus based on ORF2 amino acid sequence p-distances did not reliably distinguish several established species and was of limited use to identify distinct groups among unclassified astroviruses that were isolated recently, predominantly from cattle and pigs. A 17% nucleotide sequence distance cut-off in ORF1b fairly distinguished the established species and several groups among the unclassified viruses, providing better correspondence between phylogenetic grouping, reproductive isolation and the virus hosts. Sequence distance criteria (17% in nucleotide sequence of ORF1b and 25% in amino acid sequence of ORF2) and the recombination pattern corresponded fairly well as species criteria, but all had minor exclusions among mammalian astroviruses. A combination of these taxonomic criteria supported the established Mamastrovirus species and suggested redefining a few provisional species that were proposed earlier and introducing at least six novel species among recently submitted rat and bovine astroviruses.

Expanding the genomic diversity of human anelloviruses

Anelloviruses are a group of small, circular, single-stranded DNA viruses that are found ubiquitously across mammalian hosts. Here, we explored a large number of publicly available human microbiome datasets and retrieved a total of 829 anellovirus genomes, substantially expanding the known diversity of these viruses. The majority of new genomes fall within the three major human anellovirus genera: Alphatorquevirus, Betatorquevirus, and Gammatorquevirus, while we also present new genomes of the under-sampled Hetorquevirus, Memtorquevirus, and Samektorquevirus genera. We performed recombination analysis and show evidence of extensive recombination across all human anelloviruses. Interestingly, more than 95% of the detected events are between members of the same genus and only 15 inter-genus recombination events were detected. The breakpoints of recombination cluster in hotspots at the ends and outside of the ORF1 gene, while a recombination coldspot was detected within the gene. Our analysis suggests that anellovirus evolution is governed by homologous recombination; however, events between distant viruses or ones producing chimaeric ORF1s likely lead to nonviable recombinants. The large number of genomes further allowed us to examine how essential genomic features vary across anelloviruses. These include functional domains in the ORF1 protein and the nucleotide motif of the replication loop region, required for the viruses' rolling-circle replication. A subset of the genomes assembled in both this and previous studies are completely lacking these essential elements, opening up the possibility that anellovirus intracellular populations contain nonstandard viral genomes. However, low-read depth of the metagenomically assembled contigs may partly explain the lack of some features. Overall, our study highlights key features of anellovirus genomics and evolution, a largely understudied group of viruses whose potential in virus-based therapeutics is recently being explored.

Comprehensive phylogenomic analysis of Zika virus: Insights into its origin, past evolutionary dynamics, and global spread

BACKGROUND

Zika virus (ZIKV), a Flaviviridae family member, has been linked to severe neurological disorders. Despite detailed studies on recent outbreaks, the early evolutionary history of ZIKV remains partially unclear. This study elucidates ZIKV origin and evolutionary dynamics, focusing on recombination events, early lineage diversification, and virus spread across continents.

METHODS

We assessed recombination using multiple methods. We conducted Bayesian phylogenetic analyses to understand the evolutionary relationships and timing of key diversification events. Model selection was carried out to determine the most appropriate evolutionary model for our dataset.

RESULTS

Our phylogenies revealed recent recombination between Singaporean and African lineages, indicating the co-circulation of diverse lineages during outbreaks. Thailand was identified as a crucial hub in the spread across Asia. The phylogenetic analysis suggests that the ZIKV lineage dates back to the eleventh century, with the first significant diversification occurring in the nineteenth century. The timing of the re-introduction of the Asian lineage into Africa and the delay between probable introduction and outbreak onset were also determined.

CONCLUSIONS

This study provides novel insights into ZIKV's origin and early evolutionary dynamics, highlighting Thailand's role in the spread of the virus in Asia and recent recombination events between distant lineages. These findings emphasize the need for continuous surveillance and a better understanding of ZIKV biology to forecast and mitigate future outbreaks.

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.

Untangling lineage introductions, persistence and transmission drivers of HP-PRRSV sublineage 8.7

Despite a rapid expansion of Porcine reproductive and respiratory syndrome virus (PRRSV) sublineage 8.7 over recent years, very little is known about the patterns of virus evolution, dispersal, and the factors influencing this dispersal. Relying on a national PRRSV surveillance project established over 20 years ago, we expand the available genomic data of sublineage 8.7 from China. We perform independent interlineage and intralineage recombination analyses for the entire study period, which showed a heterogeneous recombination pattern. A series of Bayesian phylogeographic analyses uncover the role of Guangdong as an important infection hub within Asia. The spatial spread of PRRSV is highly linked with a composite of human activities and the heterogeneous provincial distribution of the swine industry, largely propelled by the smaller-scale Chinese rural farming systems in the past years. We sequence all four available modified live vaccines (MLVs) and perform genomic analyses with publicly available data, of which our results suggest a key "leaky" period spanning 2011-2017 with two concurrent amino acid mutations in ORF1a 957 and ORF2 250. Overall, our study provides an in-depth overview of the evolution, transmission dynamics, and potential leaky status of HP-PRRS MLVs, providing critical insights into new MLV development.

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.

The virulence plasmid associated with AHPND in shrimp appears to have originated from Vibrio owensii through a process of homologous recombination of parental plasmids and the transposable insertion of two large fragments

Acute hepatopancreatic necrosis disease (AHPND) is a highly contagious and lethal disease of shrimp caused by Vibrio strains carrying the virulence plasmid (pAHPND) containing the pirAB virulence genes. Through analysis of plasmid sequence similarity, clustering, and phylogeny, a horizontal transfer element similar to IS91 was discovered within the pAHPND plasmid. Additionally, two distinct clades of plasmids related to pAHPND (designated as pAHPND-r1 and pAHPND-r2) were identified, which may serve as potential parental plasmids for pAHPND. The available evidence, including the difference in G+C content between the plasmid and its host, codon usage preference, and plasmid recombination event prediction, suggests that the formation of the pAHPND plasmid in the Vibrio owensii strain was likely due to the synergistic effect of the recombinase RecA and the associated proteins RecBCD on the pAHPND-r1 and pAHPND-r2, resulting in the recombination and formation of the precursor plasmid for pAHPND (pre-pAHPND). The emergence of pAHPND was found to be a result of successive insertions of the horizontal transfer elements of pirAB-Tn903 and IS91-like segment, which led to the deletion of one third of the pre-pAHPND. This plasmid was then able to spread horizontally to other Vibrio strains, contributing to the epidemics of AHPND. These findings shed light on previously unknown mechanisms involved in the emergence of pAHPND and improve our understanding of the disease's spread.

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.

Recombination Events Among SARS-CoV-2 Omicron Subvariants: Impact on Spike Interaction With ACE2 Receptor and Neutralizing Antibodies

The recombination plays a key role in promoting evolution of RNA viruses and emergence of potentially epidemic variants. Some studies investigated the recombination occurrence among SARS-CoV-2, without exploring its impact on virus-host interaction. In the aim to investigate the burden of recombination in terms of frequency and distribution, the occurrence of recombination was first explored in 44 230 Omicron sequences among BQ subvariants and the under investigation "ML" (Multiple Lineages) denoted sequences, using 3seq software. Second, the recombination impact on interaction between the Spike protein and ACE2 receptor as well as neutralizing antibodies (nAbs), was analyzed using docking tools. Recombination was detected in 56.91% and 82.20% of BQ and ML strains, respectively. It took place mainly in spike and ORF1a genes. For BQ recombinant strains, the docking analysis showed that the spike interacted strongly with ACE2 and weakly with nAbs. The mutations S373P, S375F and T376A constitute a residue network that enhances the RBD interaction with ACE2. Thirteen mutations in RBD (S373P, S375F, T376A, D405N, R408S, K417N, N440K, S477N, P494S, Q498R, N501Y, and Y505H) and NTD (Y240H) seem to be implicated in immune evasion of recombinants by altering spike interaction with nAbs. In conclusion, this "in silico" study demonstrated that the recombination mechanism is frequent among Omicron BQ and ML variants. It highlights new key mutations, that potentially implicated in enhancement of spike binding to ACE2 (F376A) and escape from nAbs (RBD: F376A, D405N, R408S, N440K, S477N, P494S, and Y505H; NTD: Y240H). Our findings present considerable insights for the elaboration of effective prophylaxis and therapeutic strategies against future SARS-CoV-2 waves.

Intra- and inter-host origin, evolution dynamics and spatial-temporal transmission characteristics of circoviruses

Introduction

Since their identification in 1974, circoviruses have caused clinicopathological diseases in various animal species, including humans. However, their origin, transmission, and genetic evolution remain poorly understood.

Methods

In this study, the genome sequences of circovirus were obtained from GenBank, and the Bayesian stochastic search variable selection algorithm was employed to analyzed the evolution and origin of circovirus.

Results

Here, the evolutionary origin, mode of transmission, and genetic recombination of the circovirus were determined based on the available circovirus genome sequences. The origin of circoviruses can be traced back to fish circovirus, which might derive from fish genome, and human contributes to transmission of fish circovirus to other species. Furthermore, mosquitos, ticks, bats, and/or rodents might play a role as intermediate hosts in circovirus intra- and inter-species transmission. Two major lineages (A and B) of circoviruses are identified, and frequent recombination events accelerate their variation and spread. The time to the most recent common ancestor of circoviruses can be traced back to around A.D. 600 and has been evolving at a rate of 10-4 substitutions site-1 year-1 for a long time.

Discussion

These comprehensive findings shed light on the evolutionary origin, population dynamics, transmission model, and genetic recombination of the circovirus providing valuable insights for the development of prevention and control strategies against circovirus infections.

Unveiling mungbean yellow mosaic virus: molecular insights and infectivity validation in mung bean (Vigna radiata) via infectious clones

Yellow mosaic disease (YMD) with typical symptoms of alternating bright yellow to green patches associated with stunting, downward cupping, and wrinkling has been observed in mung bean on agricultural farms in Coimbatore, Tamil Nadu, India. PCR using gene-specific primers indicated the presence of the yellow mosaic virus in symptomatic plants. Rolling circle amplification (RCA) followed by restriction digestion detected ~2.7 kb of DNA-A and DNA-B, allowing the identification of a bipartite genome. The full-length genome sequences were deposited in NCBI GenBank with the accession numbers MK317961 (DNA-A) and MK317962 (DNA-B). Sequence analysis of DNA-A showed the highest sequence identity of 98.39% to the DNA-A of mungbean yellow mosaic virus (MYMV)-Vigna radiata (MW736047), while DNA-B exhibited the highest level of identity (98.21%) to the MYMV-Vigna aconitifolia isolate (DQ865203) reported from Tamil Nadu. Recombinant analysis revealed distinct evidence of recombinant breakpoints of DNA-A within the region encoding the open reading frame (ORF) AC2 (transcription activation protein), with the major parent identified as MYMV-PA1 (KC9111717) and the potential minor parent as MYMV-Namakkal (DQ86520.1). Interestingly, a recombination event in the common region (CR) of DNA-B, which encodes the nuclear shuttle protein and the movement protein, was detected. MYMIV-M120 (FM202447) and MYMV-Vigna (AJ132574) were identified as the event's major and minor parents, respectively. This large variation in DNA-B led us to suspect a recombination in DNA-B. Dimeric MYMV infectious clones were constructed, and the infectivity was confirmed through agroinoculation. In future prospects, unless relying on screening using whiteflies, breeders and plant pathologists can readily use this agroinoculation procedure to identify resistant and susceptible cultivars to YMD.

Evolutionary variation in gene conversion at the avian MHC is explained by fluctuating selection, gene copy numbers and life history

The major histocompatibility complex (MHC) multigene family encodes key pathogen-recognition molecules of the vertebrate adaptive immune system. Hyper-polymorphism of MHC genes is de novo generated by point mutations, but new haplotypes may also arise by re-shuffling of existing variation through intra- and inter-locus gene conversion. Although the occurrence of gene conversion at the MHC has been known for decades, we still have limited understanding of its functional importance. Here, I took advantage of extensive genetic resources (~9000 sequences) to investigate broad scale macroevolutionary patterns in gene conversion processes at the MHC across nearly 200 avian species. Gene conversion was found to constitute a universal mechanism in birds, as 83% of species showed footprints of gene conversion at either MHC class and 25% of all allelic variants were attributed to gene conversion. Gene conversion processes were stronger at MHC-II than MHC-I, but inter-specific variation at both MHC classes was explained by similar evolutionary scenarios, reflecting fluctuating selection towards different optima and drift. Gene conversion showed uneven phylogenetic distribution across birds and was driven by gene copy number variation, supporting significant role of inter-locus gene conversion processes in the evolution of the avian MHC. Finally, MHC gene conversion was stronger in species with fast life histories (high fecundity) and in long-distance migrants, likely reflecting variation in population sizes and host-pathogen coevolutionary dynamics. The results provide a robust comparative framework for understanding macroevolutionary variation in gene conversion at the avian MHC and reinforce important contribution of this mechanism to functional MHC diversity.

Characterization of Six Ampeloviruses Infecting Pineapple in Reunion Island Using a Combination of High-Throughput Sequencing Approaches

The cultivation of pineapple (Ananas comosus) is threatened worldwide by mealybug wilt disease of pineapple (MWP), whose etiology is not yet fully elucidated. In this study, we characterized pineapple mealybug wilt-associated ampeloviruses (PMWaVs, family Closteroviridae) from a diseased pineapple plant collected from Reunion Island, using a high-throughput sequencing approach combining Illumina short reads and Nanopore long reads. Reads co-assembly resulted in complete or near-complete genomes for six distinct ampeloviruses, including the first complete genome of pineapple mealybug wilt-associated virus 5 (PMWaV5) and that of a new species tentatively named pineapple mealybug wilt-associated virus 7 (PMWaV7). Short reads data provided high genome coverage and sequencing depths for all six viral genomes, contrary to long reads data. The 5' and 3' ends of the genome for most of the six ampeloviruses could be recovered from long reads, providing an alternative to RACE-PCRs. Phylogenetic analyses did not unveil any geographic structuring of the diversity of PMWaV1, PMWaV2 and PMWaV3 isolates, supporting the current hypothesis that PMWaVs were mainly spread by human activity and vegetative propagation.

Exploring the Accuracy and Limits of Algorithms for Localizing Recombination Breakpoints

Phylogenetic methods are widely used to reconstruct the evolutionary relationships among species and individuals. However, recombination can obscure ancestral relationships as individuals may inherit different regions of their genome from different ancestors. It is, therefore, often necessary to detect recombination events, locate recombination breakpoints, and select recombination-free alignments prior to reconstructing phylogenetic trees. While many earlier studies have examined the power of different methods to detect recombination, very few have examined the ability of these methods to accurately locate recombination breakpoints. In this study, we simulated genome sequences based on ancestral recombination graphs and explored the accuracy of three popular recombination detection methods: MaxChi, 3SEQ, and Genetic Algorithm Recombination Detection. The accuracy of inferred breakpoint locations was evaluated along with the key factors contributing to variation in accuracy across datasets. While many different genomic features contribute to the variation in performance across methods, the number of informative sites consistent with the pattern of inheritance between parent and recombinant child sequences always has the greatest contribution to accuracy. While partitioning sequence alignments based on identified recombination breakpoints can greatly decrease phylogenetic error, the quality of phylogenetic reconstructions depends very little on how breakpoints are chosen to partition the alignment. Our work sheds light on how different features of recombinant genomes affect the performance of recombination detection methods and suggests best practices for reconstructing phylogenies based on recombination-free alignments.

Recombination and selection trajectory of the monkeypox virus during its adaptation in the human population

Monkeypox, caused by the monkeypox virus (MPXV), was historically confined to West and Central Africa but has now spread globally. Recombination and selection play crucial roles in the evolutionary adaptation of MPXV; however, the evolution of MPXV and its relationship with the recent, ground-breaking monkeypox epidemic remains poorly understood. To gain insights into the evolutionary dynamics of MPXV, comprehensive in silico recombination and selection analyses were conducted based on MPXV whole genome sequence data. Three types of recombination were identified: five ancestor-sharing interspecies recombination events, six specific interspecies recombination events and four intraspecies recombination events. The results highlight the prevalent occurrence of recombination in MPXV, with 73.3% occurring in variable regions of the genome. Selection analysis was performed from three dimensions: proteins around recombination regions, proteins from recombinant ancestors and MPXV branches, and whole-genome gene analysis. Results revealed 2 and 7 proteins under positive selection in the first two dimensions, respectively. These proteins are mainly involved in infection immunity, apoptosis regulation and viral virulence. Whole-genome analysis detected 25 genes under positive selection, mainly associated with immune response and viral regulation. Understanding their evolutionary patterns will help predict and prevent cross-species transmission, zoonotic outbreaks and potential human epidemics.

The pigeon circovirus evolution, epidemiology and interaction with the host immune system under One Loft Race rearing conditions

This study was aimed to investigate the frequency of PiCV recombination, the kinetics of PiCV viremia and shedding and the correlation between viral replication and host immune response in young pigeons subclinically infected with various PiCV variants and kept under conditions mimicking the OLR system. Fifteen racing pigeons originating from five breeding facilities were housed together for six weeks. Blood and cloacal swab samples were collected from birds every seven days to recover complete PiCV genomes and determine PiCV genetic diversity and recombination dynamics, as well as to assess virus shedding rate, level of viremia, expression of selected genes and level of anti-PiCV antibodies. Three hundred and eighty-eight complete PiCV genomes were obtained and thirteen genotypes were distinguished. Twenty-five recombination events were detected. Recombinants emerged during the first three weeks of the experiment which was consistent with the peak level of viremia and viral shedding. A further decrease in viremia and shedding partially corresponded with IFN-γ and MX1 gene expression and antibody dynamics. Considering the role of OLR pigeon rearing system in spreading infectious agents and allowing their recombination, it would be reasonable to reflect on the relevance of pigeon racing from both an animal welfare and epidemiological perspective.

Recombinant Viruses from the Picornaviridae Family Occurring in Racing Pigeons

Viruses from Picornaviridae family are known pathogens of poultry, although the information on their occurrence and pathogenicity in pigeons is scarce. In this research, efforts are made to broaden the knowledge on Megrivirus B and Pigeon picornavirus B prevalence, phylogenetic relationship with other avian picornaviruses and their possible connection with enteric disease in racing pigeons. As a result of Oxford Nanopore Sequencing, five Megrivirus and two pigeon picornavirus B-like genome sequences were recovered, among which three recombinant strains were detected. The recombinant fragments represented an average of 10.9% and 25.5% of the genome length of the Pigeon picornavirus B and Megrivirus B reference strains, respectively. The phylogenetic analysis revealed that pigeons are carriers of species-specific picornaviruses. TaqMan qPCR assays revealed 7.8% and 19.0% prevalence of Megrivirus B and 32.2% and 39.7% prevalence of Pigeon picornavirus B in the group of pigeons exhibiting signs of enteropathy and in the group of asymptomatic pigeons, respectively. In turn, digital droplet PCR showed a considerably higher number of genome copies of both viruses in sick than in asymptomatic pigeons. The results of quantitative analysis leave the role of picornaviruses in enteropathies of pigeons unclear.

Genetic recombination-mediated evolutionary interactions between phages of potential industrial importance and prophages of their hosts within or across the domains of Escherichia, Listeria, Salmonella, Campylobacter, and Staphylococcus

BACKGROUND

The in-depth understanding of the role of lateral genetic transfer (LGT) in phage-prophage interactions is essential to rationalizing phage applications for human and animal therapy, as well as for food and environmental safety. This in silico study aimed to detect LGT between phages of potential industrial importance and their hosts.

METHODS

A large array of genetic recombination detection algorithms, implemented in SplitsTree and RDP4, was applied to detect LGT between various Escherichia, Listeria, Salmonella, Campylobacter, Staphylococcus, Pseudomonas, and Vibrio phages and their hosts. PHASTER and RAST were employed respectively to identify prophages across the host genome and to annotate LGT-affected genes with unknown functions. PhageAI was used to gain deeper insights into the life cycle history of recombined phages.

RESULTS

The split decomposition inferences (bootstrap values: 91.3-100; fit: 91.433-100), coupled with the Phi (0.0-2.836E-12) and RDP4 (P being well below 0.05) statistics, provided strong evidence for LGT between certain Escherichia, Listeria, Salmonella, and Campylobacter virulent phages and prophages of their hosts. The LGT events entailed mainly the phage genes encoding for hypothetical proteins, while some of these genetic loci appeared to have been affected even by intergeneric recombination in specific E. coli and S. enterica virulent phages when interacting with their host prophages. Moreover, it is shown that certain L. monocytogenes virulent phages could serve at least as the donors of the gene loci, involved in encoding for the basal promoter specificity factor, for L. monocytogenes. In contrast, the large genetic clusters were determined to have been simultaneously exchanged by many S. aureus prophages and some Staphylococcus temperate phages proposed earlier as potential therapeutic candidates (in their native or modified state). The above genetic clusters were found to encompass multiple genes encoding for various proteins, such as e.g., phage tail proteins, the capsid and scaffold proteins, holins, and transcriptional terminator proteins.

CONCLUSIONS

It is suggested that phage-prophage interactions, mediated by LGT (including intergeneric recombination), can have a far-reaching impact on the co-evolutionary trajectories of industrial phages and their hosts especially when excessively present across microbially rich environments.

A unified classification system for HIV-1 5' long terminal repeats

The HIV-1 provirus mainly consists of internal coding region flanked by 1 long terminal repeats (LTRs) at each terminus. The LTRs play important roles in HIV-1 reverse transcription, integration, and transcription. However, despite of the significant study advances of the internal coding regions of HIV-1 by using definite reference classification, there are no systematic and phylogenetic classifications for HIV-1 5' LTRs, which hinders our elaboration on 5' LTR and a better understanding of the viral origin, spread and therapy. Here, by analyzing all available resources of 5' LTR sequences in public databases following 4 recognized principles for the reference classification, 83 representatives and 14 consensus sequences were identified as representatives of 2 groups, 6 subtypes, 6 sub-subtypes, and 9 CRFs. To test the reliability of the supplemented classification system, the constructed references were applied to identify the 5' LTR assignment of the 22 clinical isolates in China. The results revealed that 16 out of 22 tested strains showed a consistent subtype classification with the previous LTR-independent classification system. However, 6 strains, for which recombination events within 5' LTR were demonstrated, unexpectedly showed a different subtype classification, leading a significant change of binding sites for important transcription factors including SP1, p53, and NF-κB. The binding change of these transcriptional factors would probably affect the transcriptional activity of 5' LTR. This study supplemented a unified classification system for HIV-1 5' LTRs, which will facilitate HIV-1 characterization and be helpful for both basic and clinical research fields.

Dynamics of Potato Virus Y Infection Pressure and Strain Composition in the San Luis Valley, Colorado

The San Luis Valley (SLV), Colorado, is the second-largest fresh-potato-growing region in the United States, which accounts for about 95% of the total production in Colorado. Potato virus Y (PVY) is the leading cause of seed potato rejection in the SLV, which has caused a constant decline in seed potato production over the past two decades. To help potato growers control PVY, we monitored the dynamics of PVY infection pressure over the growing seasons of 2022 and 2023 (May through August) using tobacco bait plants exposed to field infection weekly. PVY infection dynamics were slightly different between the two seasons, but July and August had the highest infection in both years. The first PVY infection was detected in the second half of June, which coincides with the emergence of potato crops in the valley. PVY infection increased toward the beginning of August and declined toward the end of the season. Three PVY strains were identified in tobacco bait plants and potato fields, namely PVYO, PVYN-Wi, and PVYNTN. Unlike other producing areas of the United States, PVYO is still the major strain infecting potato crops in Colorado, comprising ∼40% of total PVY strain composition. This could be explained by the prevalence of the potato cultivar Russet Norkotah that lacks any identified N genes, including the Nytbr that controls PVYO, which imposes no negative selection against this strain. The current study demonstrated the usefulness of bait plants to understand PVY epidemiology and develop more targeted control practices of PVY.

Data-driven recombination detection in viral genomes

Recombination is a key molecular mechanism for the evolution and adaptation of viruses. The first recombinant SARS-CoV-2 genomes were recognized in 2021; as of today, more than ninety SARS-CoV-2 lineages are designated as recombinant. In the wake of the COVID-19 pandemic, several methods for detecting recombination in SARS-CoV-2 have been proposed; however, none could faithfully confirm manual analyses by experts in the field. We hereby present RecombinHunt, an original data-driven method for the identification of recombinant genomes, capable of recognizing recombinant SARS-CoV-2 genomes (or lineages) with one or two breakpoints with high accuracy and within reduced turn-around times. ReconbinHunt shows high specificity and sensitivity, compares favorably with other state-of-the-art methods, and faithfully confirms manual analyses by experts. RecombinHunt identifies recombinant viral genomes from the recent monkeypox epidemic in high concordance with manually curated analyses by experts, suggesting that our approach is robust and can be applied to any epidemic/pandemic virus.

Whole genome sequencing of recombinant viruses obtained from co-infection and superinfection of Vero cells with modified vaccinia virus ankara vectored influenza vaccine and a naturally occurring cowpox virus

Modified vaccinia virus Ankara (MVA) has been widely tested in clinical trials as recombinant vector vaccine against infectious diseases and cancers in humans and animals. However, one biosafety concern about the use of MVA vectored vaccine is the potential for MVA to recombine with naturally occurring orthopoxviruses in cells and hosts in which it multiplies poorly and, therefore, producing viruses with mosaic genomes with altered genetic and phenotypic properties. We previously conducted co-infection and superinfection experiments with MVA vectored influenza vaccine (MVA-HANP) and a feline Cowpox virus (CPXV-No-F1) in Vero cells (that were semi-permissive to MVA infection) and showed that recombination occurred in both co-infected and superinfected cells. In this study, we selected the putative recombinant viruses and performed genomic characterization of these viruses. Some putative recombinant viruses displayed plaque morphology distinct of that of the parental viruses. Our analysis demonstrated that they had mosaic genomes of different lengths. The recombinant viruses, with a genome more similar to MVA-HANP (>50%), rescued deleted and/or fragmented genes in MVA and gained new host ranges genes. Our analysis also revealed that some MVA-HANP contained a partially deleted transgene expression cassette and one recombinant virus contained part of the transgene expression cassette similar to that incomplete MVA-HANP. The recombination in co-infected and superinfected Vero cells resulted in recombinant viruses with unpredictable biological and genetic properties as well as recovery of delete/fragmented genes in MVA and transfer of the transgene into replication competent CPXV. These results are relevant to hazard characterization and risk assessment of MVA vectored biologicals.

Comprehensive Genomics Investigation of Neboviruses Reveals Distinct Codon Usage Patterns and Host Specificity

Neboviruses (NeVs) from the Caliciviridae family have been linked to enteric diseases in bovines and have been detected worldwide. As viruses rely entirely on the cellular machinery of the host for replication, their ability to thrive in a specific host is greatly impacted by the specific codon usage preferences. Here, we systematically analyzed the codon usage bias in NeVs to explore the genetic and evolutionary patterns. Relative Synonymous Codon Usage and Effective Number of Codon analyses indicated a marginally lower codon usage bias in NeVs, predominantly influenced by the nucleotide compositional constraints. Nonetheless, NeVs showed a higher codon usage bias for codons containing G/C at the third codon position. The neutrality plot analysis revealed natural selection as the primary factor that shaped the codon usage bias in both the VP1 (82%) and VP2 (57%) genes of NeVs. Furthermore, the NeVs showed a highly comparable codon usage pattern to bovines, as reflected through Codon Adaptation Index and Relative Codon Deoptimization Index analyses. Notably, yak NeVs showed considerably different nucleotide compositional constraints and mutational pressure compared to bovine NeVs, which appear to be predominantly host-driven. This study sheds light on the genetic mechanism driving NeVs' adaptability, evolution, and fitness to their host species.

Genome biology and evolution of mating-type loci in four cereal rust fungi

Permanent heterozygous loci, such as sex- or mating-compatibility regions, often display suppression of recombination and signals of genomic degeneration. In Basidiomycota, two distinct loci confer mating compatibility. These loci encode homeodomain (HD) transcription factors and pheromone receptor (Pra)-ligand allele pairs. To date, an analysis of genome level mating-type (MAT) loci is lacking for obligate biotrophic basidiomycetes in the Pucciniales, an order containing serious agricultural plant pathogens. Here, we focus on four species of Puccinia that infect oat and wheat, including P. coronata f. sp. avenae, P. graminis f. sp. tritici, P. triticina and P. striiformis f. sp. tritici. MAT loci are located on two separate chromosomes supporting previous hypotheses of a tetrapolar mating compatibility system in the Pucciniales. The HD genes are multiallelic in all four species while the PR locus appears biallelic, except for P. graminis f. sp. tritici, which potentially has multiple alleles. HD loci are largely conserved in their macrosynteny, both within and between species, without strong signals of recombination suppression. Regions proximal to the PR locus, however, displayed signs of recombination suppression and genomic degeneration in the three species with a biallelic PR locus. Our observations support a link between recombination suppression, genomic degeneration, and allele diversity of MAT loci that is consistent with recent mathematical modelling and simulations. Finally, we confirm that MAT genes are expressed during the asexual infection cycle, and we propose that this may support regulating nuclear maintenance and pairing during infection and spore formation. Our study provides insights into the evolution of MAT loci of key pathogenic Puccinia species. Understanding mating compatibility can help predict possible combinations of nuclear pairs, generated by sexual reproduction or somatic recombination, and the potential evolution of new virulent isolates of these important plant pathogens.

The effect of mitochondrial recombination on fertilization success in blue mussels

The presence of doubly uniparental inheritance (DUI) in bivalves represents a unique mode of mitochondrial transmission, whereby paternal (male-transmitted M-type) and maternal (female-transmitted F-type) haplotypes are transmitted to offspring separately. Male embryos retain both haplotypes, but the M-type is selectively removed from females. Due to the presence of heteroplasmy in males, mtDNA can recombine resulting in a 'masculinized' haplotype referred to as Mf-type. While mtDNA recombination is usually rare, it has been recorded in multiple mussel species across the Northern Hemisphere. Given that mitochondria are the powerhouse of the cell, different mtDNA haplotypes may have different selective advantages under diverse environmental conditions. This may be particularly important for sperm fitness and fertilization success. In this study we aimed to i) determine the presence, prevalence of the Mf-type in Australian blue mussels (Mytilus sp.) and ii) investigate the effect of Mf-mtDNA on sperm performance (a fitness correlate). We found a high prevalence of recombined mtDNA (≈35 %) located within the control region of the mitochondrial genome, which occurred only in specimens that contained Southern Hemisphere mtDNA. The presence of two female mitotypes were identified in the studied mussels, one likely originating from the Northern Hemisphere, and the other either representing the endemic M. planulatus species or introduced genotypes from the Southern Hemisphere. Despite having recombination events present in a third of the studied population, analysis of sperm performance indicated no difference in fertilization success related to mitotype.

Genome sequences of six recombinant variants of potato virus Y identified in North American potato cultivars grown in China

Six genome sequences for potato virus Y (PVY) recombinants are reported from two North American potato cultivars grown in China. The coding complete sequences encode a single open reading frame characteristic of potyviruses. The six sequenced PVY isolates represent three distinct recombinants of PVY, namely N-Wi, SYR-I, and SYR-II.

Recent selection created distinctive variability patterns on MHC class II loci in three dolphin species from the Mediterranean Sea

The major histocompatibility complex (MHC) includes highly polymorphic genes involved in antigen presentation, which is crucial for adaptive immune response. They represent fitness related genetic markers particularly informative for populations exposed to environmental challenges. Here we analyse the diversity and evolutionary traits of MHC class II DQA and DQB genes in the dolphins Stenella coeruleoalba and Grampus griseus from the Mediterranean Sea. We found substantial nucleotide and functional diversity, as well as strong evidence of balancing selection indicated by allele and supertype frequencies, Tajima's D statistics and dN/dS tests. The Risso's dolphin, considered the least abundant in the region, showed the effect of divergent allele advantage at the nucleotide and functional-peptide levels. An outstanding polymorphism was found in the striped dolphin, particularly intriguing in the DQA gene where the Ewens-Watterson test detected a selection sweep that occurred in recent history. We hypothesize that morbillivirus, which has recurrently invaded Mediterranean populations over the last decades, exerted the detected selective pressure.

Evaluation of recombination detection methods for viral sequencing

Recombination is a key evolutionary driver in shaping novel viral populations and lineages. When unaccounted for, recombination can impact evolutionary estimations or complicate their interpretation. Therefore, identifying signals for recombination in sequencing data is a key prerequisite to further analyses. A repertoire of recombination detection methods (RDMs) have been developed over the past two decades; however, the prevalence of pandemic-scale viral sequencing data poses a computational challenge for existing methods. Here, we assessed eight RDMs: PhiPack (Profile), 3SEQ, GENECONV, recombination detection program (RDP) (OpenRDP), MaxChi (OpenRDP), Chimaera (OpenRDP), UCHIME (VSEARCH), and gmos; to determine if any are suitable for the analysis of bulk sequencing data. To test the performance and scalability of these methods, we analysed simulated viral sequencing data across a range of sequence diversities, recombination frequencies, and sample sizes. Furthermore, we provide a practical example for the analysis and validation of empirical data. We find that RDMs need to be scalable, use an analytical approach and resolution that is suitable for the intended research application, and are accurate for the properties of a given dataset (e.g. sequence diversity and estimated recombination frequency). Analysis of simulated and empirical data revealed that the assessed methods exhibited considerable trade-offs between these criteria. Overall, we provide general guidelines for the validation of recombination detection results, the benefits and shortcomings of each assessed method, and future considerations for recombination detection methods for the assessment of large-scale viral sequencing data.

Evolution, genetic recombination, and phylogeography of goose parvovirus

Goose parvovirus (GPV) has garnered global attention due to its association with severe symptoms in waterfowl. However, the process underlying the global emergence and spread of GPV remains largely elusive. In this study, we illustrated the evolutionary characteristics of GPVs from a global perspective using phylogenetic analysis, recombination analysis, selection pressure analysis, and phylogeographic analysis. Our findings indicate that GPV and muscovy duck parvovirus (MDPV) diverge into two distinct branches. Within GPV, there are two classifications: classical GPV (C-GPV) and novel GPV (N-GPV), each containing three subgroups, underscoring the significant genetic diversity of GPV. Recombination analysis revealed 11 recombination events, suggesting C-GPV, N-GPV, and MDPV co-infections. Further, phylogeographic analysis revealed that China is an important exporter of GPV and that trade might serve as a potential transmission conduit. Nonetheless, a detailed understanding of its geographic transmission dynamics warrants further investigation due to the limited scope of current genomic data in our study. This study offers novel insights into the evolutionary state and spread of GPV, holding promise for informing preventive and containment strategies against GPV infection.

Espalier: Efficient Tree Reconciliation and Ancestral Recombination Graphs Reconstruction Using Maximum Agreement Forests

In the presence of recombination individuals may inherit different regions of their genome from different ancestors, resulting in a mosaic of phylogenetic histories across their genome. Ancestral recombination graphs (ARGs) can capture how phylogenetic relationships vary across the genome due to recombination, but reconstructing ARGs from genomic sequence data is notoriously difficult. Here, we present a method for reconciling discordant phylogenetic trees and reconstructing ARGs using maximum agreement forests (MAFs). Given two discordant trees, a MAF identifies the smallest possible set of topologically concordant subtrees present in both trees. We show how discordant trees can be reconciled through their MAF in a way that retains discordances strongly supported by sequence data while eliminating conflicts likely attributable to phylogenetic noise. We further show how MAFs and our reconciliation approach can be combined to select a path of local trees across the genome that maximizes the likelihood of the genomic sequence data, minimizes discordance between neighboring local trees, and identifies the recombination events necessary to explain remaining discordances to obtain a fully connected ARG. While heuristic, our ARG reconstruction approach is often as accurate as more exact methods while being much more computationally efficient. Moreover, important demographic parameters such as recombination rates can be accurately estimated from reconstructed ARGs. Finally, we apply our approach to plant infecting RNA viruses in the genus Potyvirus to demonstrate how true recombination events can be disentangled from phylogenetic noise using our ARG reconstruction methods.

Spatial and temporal dynamics of West Nile virus between Africa and Europe

It is unclear whether West Nile virus (WNV) circulates between Africa and Europe, despite numerous studies supporting an African origin and high transmission in Europe. We integrated genomic data with geographic observations and phylogenetic and phylogeographic inferences to uncover the spatial and temporal viral dynamics of WNV between these two continents. We focused our analysis towards WNV lineages 1 (L1) and 2 (L2), the most spatially widespread and pathogenic WNV lineages. Our study shows a Northern-Western African origin of L1, with back-and-forth exchanges between West Africa and Southern-Western Europe; and a Southern African origin of L2, with one main introduction from South Africa to Europe, and no back introductions observed. We also noticed a potential overlap between L1 and L2 Eastern and Western phylogeography and two Afro-Palearctic bird migratory flyways. Future studies linking avian and mosquito species susceptibility, migratory connectivity patterns, and phylogeographic inference are suggested to elucidate the dynamics of emerging viruses.

Identification and characterization of novel alphacoronaviruses in Tadarida brasiliensis (Chiroptera, Molossidae) from Argentina: insights into recombination as a mechanism favoring bat coronavirus cross-species transmission

Bats are reservoirs of various coronaviruses that can jump between bat species or other mammalian hosts, including humans. This article explores coronavirus infection in three bat species (Tadarida brasiliensis, Eumops bonariensis, and Molossus molossus) of the family Molossidae from Argentina using whole viral metagenome analysis. Fecal samples of 47 bats from three semiurban or highly urbanized areas of the province of Santa Fe were investigated. After viral particle enrichment, total RNA was sequenced using the Illumina NextSeq 550 instrument; the reads were assembled into contigs and taxonomically and phylogenetically analyzed. Three novel complete Alphacoronavirus (AlphaCoV) genomes (Tb1-3) and two partial sequences were identified in T. brasiliensis (Tb4-5), and an additional four partial sequences were identified in M. molossus (Mm1-4). Phylogenomic analysis showed that the novel AlphaCoV clustered in two different lineages distinct from the 15 officially recognized AlphaCoV subgenera. Tb2 and Tb3 isolates appeared to be variants of the same virus, probably involved in a persistent infectious cycle within the T. brasiliensis colony. Using recombination analysis, we detected a statistically significant event in Spike gene, which was reinforced by phylogenetic tree incongruence analysis, involving novel Tb1 and AlphaCoVs identified in Eptesicus fuscus (family Vespertilionidae) from the U.S. The putative recombinant region is in the S1 subdomain of the Spike gene, encompassing the potential receptor-binding domain of AlphaCoVs. This study reports the first AlphaCoV genomes in molossids from the Americas and provides new insights into recombination as an important mode of evolution of coronaviruses involved in cross-species transmission. IMPORTANCE This study generated three novel complete AlphaCoV genomes (Tb1, Tb2, and Tb3 isolates) identified in individuals of Tadarida brasiliensis from Argentina, which showed two different evolutionary patterns and are the first to be reported in the family Molossidae in the Americas. The novel Tb1 isolate was found to be involved in a putative recombination event with alphacoronaviruses identified in bats of the genus Eptesicus from the U.S., whereas isolates Tb2 and Tb3 were found in different collection seasons and might be involved in persistent viral infections in the bat colony. These findings contribute to our knowledge of the global diversity of bat coronaviruses in poorly studied species and highlight the different evolutionary aspects of AlphaCoVs circulating in bat populations in Argentina.

Genotype and Phenotype Characterization of Rhinolophus sp. Sarbecoviruses from Vietnam: Implications for Coronavirus Emergence

Bats are a major reservoir of zoonotic viruses, including coronaviruses. Since the emergence of SARS-CoV in 2002/2003 in Asia, important efforts have been made to describe the diversity of Coronaviridae circulating in bats worldwide, leading to the discovery of the precursors of epidemic and pandemic sarbecoviruses in horseshoe bats. We investigated the viral communities infecting horseshoe bats living in Northern Vietnam, and report here the first identification of sarbecoviruses in Rhinolophus thomasi and Rhinolophus siamensis bats. Phylogenetic characterization of seven strains of Vietnamese sarbecoviruses identified at least three clusters of viruses. Recombination and cross-species transmission between bats seemed to constitute major drivers of virus evolution. Vietnamese sarbecoviruses were mainly enteric, therefore constituting a risk of spillover for guano collectors or people visiting caves. To evaluate the zoonotic potential of these viruses, we analyzed in silico and in vitro the ability of their RBDs to bind to mammalian ACE2s and concluded that these viruses are likely restricted to their bat hosts. The workflow applied here to characterize the spillover potential of novel sarbecoviruses is of major interest for each time a new virus is discovered, in order to concentrate surveillance efforts on high-risk interfaces.

Panoramic analysis of coronaviruses carried by representative bat species in Southern China to better understand the coronavirus sphere

Bats, recognized as considerable reservoirs for coronaviruses (CoVs), serve as natural hosts for several highly pathogenic CoVs, including SARS-CoV and SARS-CoV-2. Investigating the bat CoV community provides insights into the origin for highly pathogenic CoVs and highlights bat CoVs with potential spillover risks. This study probes the evolution, recombination, host range, geographical distribution, and cross-species transmission characteristics of bat CoVs across China and its associated CoVs in other regions. Through detailed research on 13,064 bat samples from 14 provinces of China, 1141 CoV strains are found across 10 subgenera and one unclassified Alpha-CoV, generating 399 complete genome sequences. Within bat CoVs, 11 new CoV species are identified and 425 recombination events are detected. Bats in southern China, particularly in Yunnan province, exhibit a pronounced diversity of CoVs. Limited sampling and low detection rates exist for CoVs in Myotacovirus, Nyctacovirus, Hibecovirus, Nobecovirus in China. The genus Myotis is highlighted as a potential ancestral host for Alpha-CoV, with the genus Hipposideros suggested as a likely progenitor host for bat-associated Beta-CoV, indicating the complexity of cross-species transmission dynamics. Through the comprehensive analysis, this study enriches the understanding of bat CoVs and offers a valuable resource for future research.

An Evolutionary Perspective of Codon Usage Pattern, Dinucleotide Composition and Codon Pair Bias in Prunus Necrotic Ringspot Virus

Prunus necrotic ringspot virus (PNRSV) is a significant virus of ornamental plants and fruit trees. It is essential to study this virus due to its impact on the horticultural industry. Several studies on PNRSV diversity and phytosanitary detection technology were reported, but the content on the codon usage bias (CUB), dinucleotide preference and codon pair bias (CPB) of PNRSV is still uncertain. We performed comprehensive analyses on a dataset consisting of 359 coat protein (CP) gene sequences in PNRSV to examine the characteristics of CUB, dinucleotide composition, and CPB. The CUB analysis of PNRSV CP sequences showed that it was not only affected by natural selection, but also affected by mutations, and natural selection played a more significant role compared to mutations as the driving force. The dinucleotide composition analysis showed an over-expression of the CpC/GpA dinucleotides and an under-expression of the UpA/GpC dinucleotides. The dinucleotide composition of the PNRSV CP gene showed a weak association with the viral lineages and hosts, but a strong association with viral codon positions. Furthermore, the CPB of PNRSV CP gene is low and is related to dinucleotide preference and codon usage patterns. This research provides reference for future research on PNRSV genetic diversity and gene evolution mechanism.

Viral Diversity and Epidemiology in Critically Endangered Yangtze Finless Porpoises (Neophocaena asiaeorientalis asiaeorientalis)

The Yangtze finless porpoise (YFP) (Neophocaena asiaeorientalis asiaeorientalis) is a critically endangered freshwater cetacean, with about 1,249 individuals thought to be left in the wild. However, viral entities and viral diseases of YFPs remain obscure. In this study, anal swabs for virome analysis were collected during the physical examination of YFPs in the Tian-E-Zhou Oxbow (TEO) ex situ reserve. A total of 19 eukaryotic viral species belonging to 9 families, including Papillomaviridae, Herpesviridae, Picornaviridae, Picobirnaviridae, Caliciviridae, Retroviridae, Parvoviridae, Virgaviridae, and Narnaviridae, and other unclassified viruses were identified based on metasequencing. Among these detected viruses, a novel herpesvirus (NaHV), two different kobuviruses (NaKV1-2), and six different papillomaviruses (NaPV1 to -6) were considered potential risks to YFPs and confirmed by PCR or reverse transcription-PCR (RT-PCR). Most YFPs sampled were found to harbor one or more kinds of detected viral genomes (52/58 [89.7%]). Surveillance results demonstrated that kobuvirus and herpesvirus displayed obvious age distribution and PVs showed significant gender difference in YFPs. According to species demarcation criteria in individual genera in Papillomaviridae, two novel species (referred to as Omikronpapillomavirus 2 and 3) and four novel isolates of PV were identified in YFPs. Further evolutionary analysis suggested that NaPVs would occupy the mucosal niche and that virus-host codivergence mixed with duplications and host-switching events drives the evolution of cetacean PVs. Divergence times of PVs in YFP and other cetacean reflect the incipient speciation of YFPs. In summary, our findings revealed the potential viral entities, their prevalence, and their evolutionary history in YFPs, which raises an important issue regarding effects of viral infection on the fitness of YFPs. IMPORTANCE The Yangtze finless porpoise (YFP) is the only cetacean species in freshwater following the functional extinction of the baiji (Lipotes vexillifer). Health management, disease treatment, and other special measures are important for maintaining the existing YFP populations, especially in in situ and ex situ reserves. The discovery of potential viral entities and their prevalence in YFPs raises an important issue regarding the effects of viral infection on the fitness of YFPs and may contribute to the conservation of YFPs. The evolutionary history of papillomaviruses in YFP and other cetaceans reflects the phylogeny of their hosts and supports the status of incipient species, opening a window to investigate the evolutionary adaptation of cetaceans to freshwater as well as their phylogeny to remedy the deficiency of fossil evidence.

The 18S rRNA genes of Haemoproteus (Haemosporida, Apicomplexa) parasites from European songbirds with remarks on improved parasite diagnostics

BACKGROUND

The nuclear ribosomal RNA genes of Plasmodium parasites are assumed to evolve according to a birth-and-death model with new variants originating by duplication and others becoming deleted. For some Plasmodium species, it has been shown that distinct variants of the 18S rRNA genes are expressed differentially in vertebrate hosts and mosquito vectors. The central aim was to evaluate whether avian haemosporidian parasites of the genus Haemoproteus also have substantially distinct 18S variants, focusing on lineages belonging to the Haemoproteus majoris and Haemoproteus belopolskyi species groups.

METHODS

The almost complete 18S rRNA genes of 19 Haemoproteus lineages of the subgenus Parahaemoproteus, which are common in passeriform birds from the Palaearctic, were sequenced. The PCR products of 20 blood and tissue samples containing 19 parasite lineages were subjected to molecular cloning, and ten clones in mean were sequenced each. The sequence features were analysed and phylogenetic trees were calculated, including sequence data published previously from eight additional Parahaemoproteus lineages. The geographic and host distribution of all 27 lineages was visualised as CytB haplotype networks and pie charts. Based on the 18S sequence data, species-specific oligonucleotide probes were designed to target the parasites in host tissue by in situ hybridization assays.

RESULTS

Most Haemoproteus lineages had two or more variants of the 18S gene like many Plasmodium species, but the maximum distances between variants were generally lower. Moreover, unlike in most mammalian and avian Plasmodium species, the 18S sequences of all but one parasite lineage clustered into reciprocally monophyletic clades. Considerably distinct 18S clusters were only found in Haemoproteus tartakovskyi hSISKIN1 and Haemoproteus sp. hROFI1. The presence of chimeric 18S variants in some Haemoproteus lineages indicates that their ribosomal units rather evolve in a semi-concerted fashion than according to a strict model of birth-and-death evolution.

CONCLUSIONS

Parasites of the subgenus Parahaemoproteus contain distinct 18S variants, but the intraspecific variability is lower than in most mammalian and avian Plasmodium species. The new 18S data provides a basis for more thorough investigations on the development of Haemoproteus parasites in host tissue using in situ hybridization techniques targeting specific parasite lineages.

Genetic diversity, recombination and cross-species transmission of a waterbird gammacoronavirus in the wild

Viruses emerging from wildlife can cause outbreaks in humans and domesticated animals. Predicting the emergence of future pathogens and mitigating their impacts requires an understanding of what shapes virus diversity and dynamics in wildlife reservoirs. In order to better understand coronavirus ecology in wild species, we sampled birds within a coastal freshwater lagoon habitat across 5 years, focussing on a large population of mute swans (Cygnus olor) and the diverse species that they interact with. We discovered and characterised the full genome of a divergent gammacoronavirus belonging to the Goose coronavirus CB17 species. We investigated the genetic diversity and dynamics of this gammacoronavirus using untargeted metagenomic sequencing of 223 faecal samples from swans of known age and sex, and RT-PCR screening of 1632 additional bird samples. The virus circulated persistently within the bird community; virus prevalence in mute swans exhibited seasonal variations, but did not change with swan age-class or epidemiological year. One whole genome was fully characterised, and revealed that the virus originated from a recombination event involving an undescribed gammacoronavirus species. Multiple lineages of this gammacoronavirus co-circulated within our study population. Viruses from this species have recently been detected in aquatic birds from both the Anatidae and Rallidae families, implying that host species habitat sharing may be important in shaping virus host range. As the host range of the Goose coronavirus CB17 species is not limited to geese, we propose that this species name should be updated to 'Waterbird gammacoronavirus 1'. Non-invasive sampling of bird coronaviruses may provide a tractable model system for understanding the evolutionary and cross-species dynamics of coronaviruses.

A novel approach to finding conserved features in low-variability gene alignments characterises RNA motifs in SARS-CoV and SARS-CoV-2

Collections of genetic sequences belonging to related organisms contain information on the evolutionary constraints to which the organisms have been subjected. Heavily constrained regions can be investigated to understand their roles in an organism's life cycle, and drugs can be sought to disrupt these roles. In organisms with low genetic diversity, such as newly-emerged pathogens, it is key to obtain this information early to develop new treatments. Here, we present methods that ensure we can leverage all the information available in a low-signal, low-noise set of sequences, to find contiguous regions of relatively conserved nucleic acid. We demonstrate the application of these methods by analysing over 5 million genome sequences of the recently-emerged RNA virus SARS-CoV-2 and correlating these results with an analysis of 119 genome sequences of SARS-CoV. We propose the precise location of a previously described packaging signal, and discuss explanations for other regions of high conservation.

Metagenomics reveals the structure of Mastrevirus-host interaction network within an agro-ecosystem

As highly pervasive parasites that sometimes cause disease, viruses are likely major components of all natural ecosystems. An important step towards both understanding the precise ecological roles of viruses and determining how natural communities of viral species are assembled and evolve is obtaining full descriptions of viral diversity and distributions at ecosystem scales. Here, we focused on obtaining such 'community-scale' data for viruses in a single genus. We chose the genus Mastrevirus (family Geminiviridae), members of which have predominantly been found infecting uncultivated grasses (family Poaceae) throughout the tropical and sub-tropical regions of the world. We sampled over 3 years, 2,884 individual Poaceae plants belonging to thirty different species within a 2-ha plot which included cultivated and uncultivated areas on the island of Reunion. Mastreviruses were found in ∼8 per cent of the samples, of which 96 per cent did not have any discernible disease symptoms. The multitude of host-virus associations that we uncovered reveals both the plant species that most commonly host mastreviruses and the mastrevirus species (such as maize streak virus and maize streak Reunion virus) that have especially large host ranges. Our findings are consistent with the hypothesis that perennial plant species capable of hosting years-long mixed mastrevirus infections likely play a disproportionately important role in the generation of inter-species and inter-strain mastrevirus recombinants.

High functional allelic diversity and copy number in both MHC classes in the common buzzard

BACKGROUND

The major histocompatibility complex (MHC), which encodes molecules that recognize various pathogens and parasites and initiates the adaptive immune response in vertebrates, is renowned for its exceptional polymorphism and is a model of adaptive gene evolution. In birds, the number of MHC genes and sequence diversity varies greatly among taxa, believed due to evolutionary history and differential selection pressures. Earlier characterization studies and recent comparative studies suggest that non-passerine species have relatively few MHC gene copies compared to passerines. Additionally, comparative studies that have looked at partial MHC sequences have speculated that non-passerines have opposite patterns of selection on MHC class I (MHC-I) and class II (MHC-II) loci than passerines: namely, greater sequence diversity and signals of selection on MHC-II than MHC-I. However, new sequencing technology is revealing much greater MHC variation than previously expected while also facilitating full sequence variant detection directly from genomic data. Our study aims to take advantage of high-throughput sequencing methods to fully characterize both classes and domains of MHC of a non-passerine bird of prey, the common buzzard (Buteo buteo), to test predictions of MHC variation and differential selection on MHC classes.

RESULTS

Using genetic, genomic, and transcriptomic high-throughput sequencing data, we established common buzzards have at least three loci that produce functional alleles at both MHC classes. In total, we characterize 91 alleles from 113 common buzzard chicks for MHC-I exon 3 and 41 alleles from 125 chicks for MHC-IIB exon 2. Among these alleles, we found greater sequence polymorphism and stronger diversifying selection at MHC-IIB exon 2 than MHC-I exon 3, suggesting differential selection pressures on MHC classes. However, upon further investigation of the entire peptide-binding groove by including genomic data from MHC-I exon 2 and MHC-IIA exon 2, this turned out to be false. MHC-I exon 2 was as polymorphic as MHC-IIB exon 2 and MHC-IIA exon 2 was essentially invariant. Thus, comparisons between MHC-I and MHC-II that included both domains of the peptide-binding groove showed no differences in polymorphism nor diversifying selection between the classes. Nevertheless, selection analysis indicates balancing selection has been acting on common buzzard MHC and phylogenetic inference revealed that trans-species polymorphism is present between common buzzards and species separated for over 33 million years for class I and class II.

CONCLUSIONS

We characterize and confirm the functionality of unexpectedly high copy number and allelic diversity in both MHC classes of a bird of prey. While balancing selection is acting on both classes, there is no evidence of differential selection pressure on MHC classes in common buzzards and this result may hold more generally once more data for understudied MHC exons becomes available.

Comparative phylogenetic analysis and transcriptomic profiling of Dengue (DENV-3 genotype I) outbreak in 2021 in Bangladesh

Background

The next-generation sequencing (NGS) technology facilitates in-depth study of host-pathogen metatranscriptome. We, therefore, implicated phylodynamic and transcriptomic approaches through NGS technology to know/understand the dengue virus (DENV) origin and host response with dengue fever.

Methods

In this study, blood serum RNA was extracted from 21 dengue patients and 3 healthy individuals. Total transcriptomic data were analyzed for phylogenetic, phylodynamic, differential express gene (DEG), and gene ontology (GO) using respective bioinformatics tools.

Results

The viral genome sequence revealed dengue viral genome size ranges 10647 to 10707 nucleotide. Phylogenetic and phylodynamic analysis showed that the 2021 epidemic isolates were DENV-3 genotype-I and maintained as a new clade in compared to 2019 epidemic. Transcriptome analysis showed a total of 2686 genes were DEG in dengue patients compared to control with a q-value < 0.05. DESeq2 plot counts function of the top 24 genes with the smallest q-values of differential gene expression of RNA-seq data showed that 11 genes were upregulated, whereas 13 genes were downregulated. GO analysis showed a significant upregulation (p = < 0.001) in a process of multicellular organismal, nervous system, sensory perception of chemical stimulus, and G protein-coupled receptor signaling pathways in the dengue patients. However, there were a significant downregulation (p = < 0.001) of intracellular component, cellular anatomical entity, and protein-containing complex in dengue patients. Most importantly, there was a significant increase of a class of immunoregulatory proteins in dengue patients in compared to the controls, with increased GO of immune system process. In addition, upregulation of toll receptor (TLR) signaling pathways were found in dengue patients. These TLR pathways were particularly involved for the activation of innate system coupled with adaptive immune system that probably involved the rapid elimination of dengue virus infected cells. These differentially expressed genes could be further investigated for target based prophylactic interventions for dengue.

Conclusion

This is a first report describing DENV complete genomic features and differentially expressed genes in patients in Bangladesh. These genes may have diagnostic and therapeutic values for dengue infection. Continual genomic surveillance is required to further investigate the shift in dominant genotypes in relation to viral pathogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-023-02030-1.

Prevalence, variation, and transmission patterns of human respiratory syncytial virus from pediatric patients in Hubei, China during 2020-2021

Human respiratory syncytial virus (RSV) is a severe threat to children and a main cause of acute lower respiratory tract infections. Nevertheless, the intra-host evolution and inter-regional diffusion of RSV are little known. In this study, we performed a systematic surveillance in hospitalized children in Hubei during 2020-2021, in which 106 RSV-positive samples were detected both clinically and by metagenomic next generation sequencing (mNGS). RSV-A and RSV-B groups co-circulated during surveillance with RSV-B being predominant. About 46 high-quality genomes were used for further analyses. A total of 163 intra-host nucleotide variation (iSNV) sites distributed in 34 samples were detected, and glycoprotein (G) gene was the most enriched gene for iSNVs, with non-synonymous substitutions more than synonymous substitutions. Evolutionary dynamic analysis showed that the evolutionary rates of G and NS2 genes were higher, and the population size of RSV groups changed over time. We also found evidences of inter-regional diffusion from Europe and Oceania to Hubei for RSV-A and RSV-B, respectively. This study highlighted the intra-host and inter-host evolution of RSV, and provided some evidences for understanding the evolution of RSV.

Virological characteristics of the SARS-CoV-2 XBB variant derived from recombination of two Omicron subvariants

In late 2022, SARS-CoV-2 Omicron subvariants have become highly diversified, and XBB is spreading rapidly around the world. Our phylogenetic analyses suggested that XBB emerged through the recombination of two cocirculating BA.2 lineages, BJ.1 and BM.1.1.1 (a progeny of BA.2.75), during the summer of 2022. XBB.1 is the variant most profoundly resistant to BA.2/5 breakthrough infection sera to date and is more fusogenic than BA.2.75. The recombination breakpoint is located in the receptor-binding domain of spike, and each region of the recombinant spike confers immune evasion and increases fusogenicity. We further provide the structural basis for the interaction between XBB.1 spike and human ACE2. Finally, the intrinsic pathogenicity of XBB.1 in male hamsters is comparable to or even lower than that of BA.2.75. Our multiscale investigation provides evidence suggesting that XBB is the first observed SARS-CoV-2 variant to increase its fitness through recombination rather than substitutions.

Emergence of a novel genotype of class II New Castle Disease virus in North Eastern States of India

Outbreaks of New Castle Disease from three north eastern states of India were confirmed by clinico-pathological examination followed by reverse transcription-PCR detection of F gene of ND Virus (NDV). Irrespective of vaccination, the outbreaks resulted 90-100% mortality in the affected flocks. The analysis of fusion protein sequences from ten field isolates revealed them as the velogenic or highly virulent strain. Phylogenetic analyses based on the complete F gene nucleotide sequences of the isolates have characterized only one of the isolate (OK149201) in the genotype XIII.2.2. The rest of the nine isolates are depicted in a distinct monophyletic group with average nucleotide distances from the other 20 genotypes ranged from 10.90 - 20.70. The nine isolates were further divided into two sub branches with the bootstrap support value of 100% at the nodes that define the two subgroups with an average evolutionary nucleotide distance of 6.00between the isolates in the two subgroups. As per the recommendation put forth in recently updated unified phylogenetic classification system for NDV, our findings clearly indicates emergence of a novel genotype of class II NDV in the biodiversity hot spot region of NER, India. The isolates in the newly identified genotype is designated with next available Roman numerals XXII. Further, the two subgroups within the genotype are designated as XXII.1 and XXII.2.