An exact nonparametric method for inferring mosaic structure in sequence triplets

Phylogeography reveals association between swine trade and the spread of porcine epidemic diarrhea virus in China and across the world

The ongoing SARS (severe acute respiratory syndrome)-CoV (coronavirus)-2 pandemic has exposed major gaps in our knowledge on the origin, ecology, evolution, and spread of animal coronaviruses. Porcine epidemic diarrhea virus (PEDV) is a member of the genus Alphacoronavirus in the family Coronaviridae that may have originated from bats and leads to significant hazards and widespread epidemics in the swine population. The role of local and global trade of live swine and swine-related products in disseminating PEDV remains unclear, especially in developing countries with complex swine production systems. Here, we undertake an in-depth phylogeographic analysis of PEDV sequence data (including 247 newly sequenced samples) and employ an extension of this inference framework that enables formally testing the contribution of a range of predictor variables to the geographic spread of PEDV. Within China, the provinces of Guangdong and Henan were identified as primary hubs for the spread of PEDV, for which we estimate live swine trade to play a very important role. On a global scale, the United States and China maintain the highest number of PEDV lineages. We estimate that, after an initial introduction out of China, the United States acted as an important source of PEDV introductions into Japan, Korea, China, and Mexico. Live swine trade also explains the dispersal of PEDV on a global scale. Given the increasingly global trade of live swine, our findings have important implications for designing prevention and containment measures to combat a wide range of livestock coronaviruses.

The phylodynamics of emerging porcine deltacoronavirus (PDCoV) in Southeast Asia

Porcine deltacoronavirus (PDCoV), a recently emerging pathogen, causes diarrhea in pigs. A previous phylogenetic analysis based on spike genes suggested that PDCoV was divided into 3 different groups, including China, US, and Southeast Asia (SEA). SEA PDCoV, however, is genetically separated from China and US but shares a common ancestor. Its origin and evolution have yet been identified. Herein, phylodynamic analyses based on the full-length genome were performed to investigate the origin and evolution of SEA PDCoV. In the study, 18 full-length genome sequences of SEA PDCoV identified in 2013-2016 together with PDCoV from other regions were used in analyses. The results demonstrated that PDCoV was classified into 2 genogroups including G1 and G2. G1 is further evolved into G1a (China), G1b (US). G2 (SEA) group is further evolved into 3 clades, including SEA-1 (Thailand), SEA-2 (Vietnam), and SEA-2r (Vietnam recombinant) clades. The time to the most recent common ancestor (MRCA) of global PDCoV was estimated to be approximately 1989-1990 and possibly have been circulated in SEA more than a decade. SEA PDCoV is genetically diverse compared to China and US PDCoV. The substitution rate of SEA PDCoV was lower than those of China and US, but the recombination rate of SEA was higher. Recombination analyses revealed 4 potential recombinant events in SEA PDCoV, suggesting that they were derived from the same ancestor of China PDCoV. The SEA-2r subgroup was potentially recombinant between SEA-2 and US strains. In conclusion, major mechanisms driving the complex evolution and genetic diversity of SEA PDCoV were multiple introductions of exotic PDCoV strains followed by recombination. This article is protected by copyright. All rights reserved.

Molecular Characterization of Small Ruminant Lentiviruses in Polish Mixed Flocks Supports Evidence of Cross Species Transmission, Dual Infection, a Recombination Event, and Reveals the Existence of New Subtypes within Group A

Small ruminant lentiviruses (SRLVs) are a group of highly divergent viruses responsible for global infection in sheep and goats. In a previous study we showed that SRLV strains found in mixed flocks in Poland belonged to subtype A13 and A18, but this study was restricted only to the few flocks from Małopolska region. The present work aimed at extending earlier findings with the analysis of SRLVs in mixed flocks including larger numbers of animals and flocks from different part of Poland. On the basis of gag and env sequences, Polish SRLVs were assigned to the subtypes B2, A5, A12, and A17. Furthermore, the existence of a new subtypes, tentatively designed as A23 and A24, were described for the first time. Subtypes A5 and A17 were only found in goats, subtype A24 has been detected only in sheep while subtypes A12, A23, and B2 have been found in both sheep and goats. Co-infection with strains belonging to different subtypes was evidenced in three sheep and two goats originating from two flocks. Furthermore, three putative recombination events were identified within gag and env SRLVs sequences derived from three sheep. Amino acid (aa) sequences of immunodominant epitopes in CA protein were well conserved while Major Homology Region (MHR) had more alteration showing unique mutations in sequences of subtypes A5 and A17. In contrast, aa sequences of surface glycoprotein exhibited higher variability confirming type-specific variation in the SU5 epitope. The number of potential N-linked glycosylation sites (PNGS) ranged from 3 to 6 in respective sequences and were located in different positions. The analysis of LTR sequences revealed that sequences corresponding to the TATA box, AP-4, AML-vis, and polyadenylation signal (poly A) were quite conserved, while considerable alteration was observed in AP-1 sites. Interestingly, our results revealed that all sequences belonging to subtype A17 had unique substitution T to A in the fifth position of TATA box and did not have a 11 nt deletion in the R region which was noted in other sequences from Poland. These data revealed a complex picture of SRLVs population with ovine and caprine strains belonging to group A and B. We present strong and multiple evidence of dually infected sheep and goats in mixed flocks and present evidence that these viruses can recombine in vivo.

Metagenomic and Recombination Analyses of Antimicrobial Resistance Genes from Recreational Waters of Black Sea Coastal Areas and Other Marine Environments Unveil Extensive Evidence for Their both Intrageneric and Intergeneric Transmission across Genetically Very Diverse Microbial Communities

Microbial communities of marine coastal recreation waters have become large reservoirs of AMR genes (ARGs), contributing to the emergence and transmission of various zoonotic, foodborne and other infections that exhibit resistance to various antibiotics. Thus, it is highly imperative to determine ARGs assemblages as well as mechanisms and trajectories of their transmission across these microbial communities for our better understanding of the evolutionary trends of AMR (AMR). In this study, using metagenomics approaches, we screened for ARGs in recreation waters of the Black Sea coastal areas of the Batumi City (Georgia). Also, a large array of the recombination detection algorithms of the SplitsTree, RDP4, and GARD was applied to elucidate genetic recombination of ARGs and trajectories of their transmission across various marine microbial communities. The metagenomics analyses of sea water samples, obtained from across the above marine sites, could identify putative ARGs encoding for multidrug resistance efflux transporters mainly from the Major Facilitator and Resistance Nodulation Division superfamilies. The data, generated by SplitsTree (fit ≥95.619; bootstrap values ≥ 95; Phi p ≤ 0.0494), RDP4 (p ≤ 0.0490), and GARD, provided strong statistical evidence not only for intrageneric recombination of these ARGs, but also for their intergeneric recombination across fairly large and diverse microbial communities of marine environment. These bacteria included both human pathogenic and nonpathogenic species, exhibiting collectively the genera of Vibrio, Aeromonas, Synechococcus, Citromicrobium, Rhodobacteraceae, Pseudoalteromonas, Altererythrobacter, Erythrobacter, Altererythrobacter, Marivivens, Xuhuaishuia, and Loktanella. The above nonpathogenic bacteria are strongly suggested to contribute to ARGs transmission in marine ecosystems.

Applicability of several rooted phylogenetic network algorithms for representing the evolutionary history of SARS-CoV-2

Background

Rooted phylogenetic networks are used to display complex evolutionary history involving so-called reticulation events, such as genetic recombination. Various methods have been developed to construct such networks, using for example a multiple sequence alignment or multiple phylogenetic trees as input data. Coronaviruses are known to recombine frequently, but rooted phylogenetic networks have not yet been used extensively to describe their evolutionary history. Here, we created a workflow to compare the evolutionary history of SARS-CoV-2 with other SARS-like viruses using several rooted phylogenetic network inference algorithms. This workflow includes filtering noise from sets of phylogenetic trees by contracting edges based on branch length and bootstrap support, followed by resolution of multifurcations. We explored the running times of the network inference algorithms, the impact of filtering on the properties of the produced networks, and attempted to derive biological insights regarding the evolution of SARS-CoV-2 from them.

Results

The network inference algorithms are capable of constructing rooted phylogenetic networks for coronavirus data, although running-time limitations require restricting such datasets to a relatively small number of taxa. Filtering generally reduces the number of reticulations in the produced networks and increases their temporal consistency. Taxon bat-SL-CoVZC45 emerges as a major and structural source of discordance in the dataset. The tested algorithms often indicate that SARS-CoV-2/RaTG13 is a tree-like clade, with possibly some reticulate activity further back in their history. A smaller number of constructed networks posit SARS-CoV-2 as a possible recombinant, although this might be a methodological artefact arising from the interaction of bat-SL-CoVZC45 discordance and the optimization criteria used.

Conclusion

Our results demonstrate that as part of a wider workflow and with careful attention paid to running time, rooted phylogenetic network algorithms are capable of producing plausible networks from coronavirus data. These networks partly corroborate existing theories about SARS-CoV-2, and partly produce new avenues for exploration regarding the location and significance of reticulate activity within the wider group of SARS-like viruses. Our workflow may serve as a model for pipelines in which phylogenetic network algorithms can be used to analyse different datasets and test different hypotheses.

Are ecological communities the seat of endosymbiont horizontal transfer and diversification? A case study with soil arthropod community

Maternally inherited endosymbionts of arthropods are one of the most abundant and diverse group of bacteria. These bacterial endosymbionts also show extensive horizontal transfer to taxonomically unrelated hosts and widespread recombination in their genomes. Such horizontal transfers can be enhanced when different arthropod hosts come in contact like in an ecological community. Higher rates of horizontal transfer can also increase the probability of recombination between endosymbionts, as they now share the same host cytoplasm. However, reports of community-wide endosymbiont data are rare as most studies choose few host taxa and specific ecological interactions among the hosts. To better understand endosymbiont spread within host populations, we investigated the incidence, diversity, extent of horizontal transfer, and recombination of three endosymbionts (Wolbachia, Cardinium, and Arsenophonus) in a specific soil arthropod community. Wolbachia strains were characterized with MLST genes whereas 16S rRNA gene was used for Cardinium and Arsenophonus. Among 3,509 individual host arthropods, belonging to 390 morphospecies, 12.05% were infected with Wolbachia, 2.82% with Cardinium and 2.05% with Arsenophonus. Phylogenetic incongruence between host and endosymbiont indicated extensive horizontal transfer of endosymbionts within this community. Three cases of recombination between Wolbachia supergroups and eight incidences of within-supergroup recombination were also found. Statistical tests of similarity indicated supergroup A Wolbachia and Cardinium show a pattern consistent with extensive horizontal transfer within the community but not for supergroup B Wolbachia and Arsenophonus. We highlight the importance of extensive community-wide studies for a better understanding of the spread of endosymbionts across global arthropod communities.

Molecular characterization of porcine epidemic diarrhoea virus (PEDV) in Poland reveals the presence of swine enteric coronavirus (SeCoV) sequence in S gene

Porcine epidemic diarrhoea (PED) is a highly contagious enteric viral disease of pigs with a high morbidity and mortality rate, which ultimately results in huge economic losses in the pig production sector. The etiological agent of this disease is the porcine epidemic diarrhoea virus (PEDV) which is an enveloped, positive single-stranded RNA virus. The aim of this study was to perform molecular characterization of PEDV to identify the strains circulating in Poland. In this study, 662 faecal samples from 2015 to 2021 were tested with reverse transcription quantitative real-time PCR (RT-qPCR) and the results showed that 3.8% of the tested samples revealed a positive result for PEDV. A phylogenetic analysis of the complete genome and complete S gene sequences showed that Polish PEDV strains belonged to the G1b (S-INDEL) subgroup and were closely related to the European PEDV strains isolated from 2014 to 2019. Furthermore, RDP4 analysis revealed that the Polish PEDV strains harboured a recombinant fragment of ~400 nt in the 5' end of S gene with PEDV and swine enteric coronavirus (SeCoV) being the major and minor parents, respectively. Antigenic analysis showed that the aa sequences of neutralizing epitopes were conserved among the Polish PEDV strains. Only one strain, #0100/5P, had a unique substitution in the COE epitope. However, Polish PEDV strains showed several substitutions, especially in the COE antigen, as compared to the classical strain CV777. To the best of our knowledge, this is the first report concerning the molecular characterization of porcine epidemic diarrhoea virus strains, as well as the first phylogenetic analysis for PEDV in Poland.

Genetic variability and evolutionary dynamics of atypical Papaya ringspot virus infecting Papaya

Papaya ringspot virus biotype-P is a detrimental pathogen of economically important papaya and cucurbits worldwide. The mutation prone feature of this virus perhaps accounts for its geographical dissemination. In this study, investigations of the atypical PRSV-P strain was conducted based on phylogenetic, recombination and genetic differentiation analyses considering of it's likely spread across India and Bangladesh. Full length genomic sequences of 38 PRSV isolates and 35 CP gene sequences were subjected to recombination analysis. A total of 61 recombination events were detected in aligned complete PRSV genome sequences. 3 events were detected in complete genome of PRSV strain PK whereas one was in its CP gene sequence. The PRSV-PK appeared to be recombinant of a major parent from Bangladesh. However, the genetic differentiation based on full length genomic sequences revealed less frequent gene flow between virus PRSV-PK and the population from America, India, Colombia, other Asian Countries and Australia. Whereas, frequent gene flow exists between Pakistan and Bangladesh virus populations. These results provided evidence correlating geographical position and genetic distances. We speculate that the genetic variations and evolutionary dynamics of this virus may challenge the resistance developed in papaya against PRSV and give rise to virus lineage because of its atypical emergence where geographic spread is already occurring.

Generation and transmission of interlineage recombinants in the SARS-CoV-2 pandemic

We present evidence for multiple independent origins of recombinant SARS-CoV-2 viruses sampled from late 2020 and early 2021 in the United Kingdom. Their genomes carry single-nucleotide polymorphisms and deletions that are characteristic of the B.1.1.7 variant of concern but lack the full complement of lineage-defining mutations. Instead, the remainder of their genomes share contiguous genetic variation with non-B.1.1.7 viruses circulating in the same geographic area at the same time as the recombinants. In four instances, there was evidence for onward transmission of a recombinant-origin virus, including one transmission cluster of 45 sequenced cases over the course of 2 months. The inferred genomic locations of recombination breakpoints suggest that every community-transmitted recombinant virus inherited its spike region from a B.1.1.7 parental virus, consistent with a transmission advantage for B.1.1.7's set of mutations.

Diversity and transmission of Aleutian mink disease virus in feral and farmed American mink and native mustelids

Aleutian mink disease virus (AMDV), which causes Aleutian disease, is widely spread both in farmed mink and wild mustelids. However, only limited data are available on the role of wild animals in AMDV transmission and spread. Our aim was to shed light on AMDV transmission among wild mustelids and estimate the effect of intense farming practices on the virus circulation by studying AMDV prevalence and genetic diversity among wild mustelids in Poland. We compared AMDV seroprevalence and proportion of PCR-positive individuals in American mink, polecats, otters, stone martens, and pine martens and used the phylogenetic analysis of the NS1 region to study transmission. In addition, we used a metagenomic approach to sequence complete AMDV genomes from tissue samples. In eastern Poland, AMDV seroprevalence in wild mustelids varied from 22 per cent in otters to 62 per cent and 64 per cent in stone martens and feral mink, respectively. All studied antibody-positive mink were also PCR positive, whereas only 10, 15, and 18 per cent of antibody-positive polecats, pine martens, and stone martens, respectively, were PCR positive, suggesting lower virus persistence among these animal species as compared to feral mink. In phylogenetic analysis, most sequences from feral mink formed region-specific clusters that have most likely emerged through multiple introductions of AMDV to feral mink population over decades. However, virus spread between regions was also observed. Virus sequences derived from farmed and wild animals formed separate subclusters in the phylogenetic tree, and no signs of recent virus transmission between farmed and wild animals were observed despite the frequent inflow of farmed mink escapees to wild populations. These results provide new information about the role of different mustelid species in AMDV transmission and about virus circulation among the wild mustelids. In addition, we pinpoint gaps of knowledge, where more studies are needed to achieve a comprehensive picture of AMDV transmission.

Variation and Selection in the Putative Sperm-Binding Region of ZP3 in Muroid Rodents: A Comparison between Cricetids and Murines

In mammals, the zona pellucida glycoprotein 3 (ZP3) is considered a primary sperm receptor of the oocyte and is hypothesized to be involved in reproductive isolation. We investigated patterns of diversity and selection in the putative sperm-binding region (pSBR) of mouse ZP3 across Cricetidae and Murinae, two hyperdiverse taxonomic groups within muroid rodents. In murines, the pSBR is fairly conserved, in particular the serine-rich stretch containing the glycosylation sites proposed as essential for sperm binding. In contrast, cricetid amino acid sequences of the pSBR were much more variable and the serine-rich motif, typical of murines, was generally substantially modified. Overall, our results suggest a general lack of species specificity of the pSBR across the two muroid families. We document statistical evidence of positive selection acting on exons 6 and 7 of ZP3 and identified several amino acid sites that are likely targets of selection, with most positively selected sites falling within or adjacent to the pSBR.

Genome-wide population structure inferences of human coxsackievirus-A; insights the genotypes diversity and evolution

Coxsackievirus-A (CV-A) is a causative agent of Hand Foot Mouth Disease (HFMD) worldwide. It belongs to the Human Enterovirus genus of the family Picornaviridae. The genomics data availability of CV-A samples, isolated from human host across different continental regions, provide an excellent opportunity to study its genetic composition, diversity, and evolutionary events. The complete genome sequences of 424 CV-A isolates were analyzed through a model-based population genetic approach implemented in the STRUCTURE program. Twelve genetically distinct sub-populations were identified for CV-A isolates with a marked Fst distinction of 0.76991 (P-value = 0.00000). Besides, genetically admixed strains were characterized in the G-Id, G-IIIb clusters constituted by the CV-A12 and CV-A6 enterovirus serotypes. The serotypes depicted inter/intra-genotype recombination and episodic positive selection signatures in the structural and non-structural protein-coding regions. The observed genetic composition of CV-A samples was also deduced by the phylogenetic tree analyses, where a uniform genetic structure was inferred for most of the CV-A genotypes. However, the CV-A6 serotype samples genetically stratified into three sub-populations that may lead to the emergence of new lineages in future. These informations may implicate in planning the effective strategies to combat the coxsackievirus-A-mediated infection.

Pirahy virus: Identification of a new and potential emerging arbovirus in South Brazil

Genomic and epidemiological surveillance are paramount for the discovery of new viruses with the potential to cross species barriers. Here, we present a new member of the genus Alphavirus found in Trichoprosopon and Wyeomia mosquitoes, tentatively named Pirahy virus (PIRAV). PIRAV was isolated from mosquito pools collected in a rural area of Piraí do Sul, South Brazil. In vitro assays revealed that PIRAV replicates and causes cytopathic effects in vertebrate cell lines such as Vero E6, SH-SY5Y, BHK-21 and UMNSAH/DF-1. Genomic signature analysis supports these results showing a dinucleotide and codon usage balance compatible with several hosts. Phylogenetic analyses placed PIRAV basal to the Venezuelan equine encephalitis complex. Genome analyses, electron microscopy, and biological characterization show findings that may alert for the emergence of a new arbovirus in South America.

Feline Leukemia Virus (FeLV) Endogenous and Exogenous Recombination Events Result in Multiple FeLV-B Subtypes during Natural Infection

Feline leukemia virus (FeLV) is associated with a range of clinical signs in felid species. Differences in disease processes are closely related to genetic variation in the envelope (env) region of the genome of six defined subgroups. The primary hosts of FeLV are domestic cats of the Felis genus that also harbor endogenous FeLV (enFeLV) elements stably integrated in their genomes. EnFeLV elements display 86% nucleotide identity to exogenous, horizontally transmitted FeLV (FeLV-A). Variation between enFeLV and FeLV-A is primarily in the long terminal repeat (LTR) and env regions, which potentiates generation of the FeLV-B recombinant subgroup during natural infection. The aim of this study was to examine recombination behavior of exogenous FeLV (exFeLV) and enFeLV in a natural FeLV epizootic. We previously described that of 65 individuals in a closed colony, 32 had productive FeLV-A infection, and 22 of these individuals had detectable circulating FeLV-B. We cloned and sequenced the env gene of FeLV-B, FeLV-A, and enFeLV spanning known recombination breakpoints and examined between 1 and 13 clones in 22 animals with FeLV-B to assess sequence diversity and recombination breakpoints. Our analysis revealed that FeLV-A sequences circulating in the population, as well as enFeLV env sequences, are highly conserved. We documented many recombination breakpoints resulting in the production of unique FeLV-B genotypes. More than half of the cats harbored more than one FeLV-B variant, suggesting multiple recombination events between enFeLV and FeLV-A. We concluded that FeLV-B was predominantly generated de novo within each host, although we could not definitively rule out horizontal transmission, as nearly all cats harbored FeLV-B sequences that were genetically highly similar to those identified in other individuals. This work represents a comprehensive analysis of endogenous-exogenous retroviral interactions with important insights into host-virus interactions that underlie disease pathogenesis in a natural setting. IMPORTANCE Feline leukemia virus (FeLV) is a felid retrovirus with a variety of disease outcomes. Exogenous FeLV-A is the virus subgroup almost exclusively transmitted between cats. Recombination between FeLV-A and endogenous FeLV analogues in the cat genome may result in emergence of largely replication-defective but highly virulent subgroups. FeLV-B is formed when the 3' envelope (env) region of endogenous FeLV (enFeLV) recombines with that of the exogenous FeLV (exFeLV) during viral reverse transcription and integration. Both domestic cats and wild relatives of the Felis genus harbor enFeLV, which has been shown to limit FeLV-A disease outcome. However, enFeLV also contributes genetic material to the recombinant FeLV-B subgroup. This study evaluates endogenous-exogenous recombination outcomes in a naturally infected closed colony of cats to determine mechanisms and risk of endogenous retroviral recombination during exogenous virus exposure that leads to enhanced virulence. While FeLV-A and enFeLV env regions were highly conserved from cat to cat, nearly all individuals with emergent FeLV-B had unique combinations of genotypes, representative of a wide range of recombination sites within env. The findings provide insight into unique recombination patterns for emergence of new pathogens and can be related to similar viruses across species.

Interspecies Recombination Has Driven the Macroevolution of Cassava Mosaic Begomoviruses

Begomoviruses (family Geminiviridae, genus Begomovirus) significantly hamper crop production and threaten food security around the world. The frequent emergence of new begomovirus genotypes is facilitated by high mutation frequencies and the propensity to recombine and reassort. Homologous recombination has been especially implicated in the emergence of novel cassava mosaic begomovirus (CMB) genotypes, which cause cassava mosaic disease (CMD). Cassava (Manihot esculenta) is a staple food crop throughout Africa and an important industrial crop in Asia, two continents where production is severely constrained by CMD. The CMD species complex is comprised of 11 bipartite begomovirus species with ample distribution throughout Africa and the Indian subcontinent. While recombination is regarded as a frequent occurrence for CMBs, a revised, systematic assessment of recombination and its impact on CMB phylogeny is currently lacking. We assembled data sets of all publicly available, full-length DNA-A (n = 880) and DNA-B (n = 369) nucleotide sequences from the 11 recognized CMB species. Phylogenetic networks and complementary recombination detection methods revealed extensive recombination among the CMB sequences. Six out of the 11 species descended from unique interspecies recombination events. Estimates of recombination and mutation rates revealed that all species experience mutation more frequently than recombination, but measures of population divergence indicate that recombination is largely responsible for the genetic differences between species. Our results support that recombination has significantly impacted the CMB phylogeny and has driven speciation in the CMD species complex. IMPORTANCE Cassava mosaic disease (CMD) is a significant threat to cassava production throughout Africa and Asia. CMD is caused by a complex comprised of 11 recognized virus species exhibiting accelerated rates of evolution, driven by high frequencies of mutation and genetic exchange. Here, we present a systematic analysis of the contribution of genetic exchange to cassava mosaic virus species-level diversity. Most of these species emerged as a result of genetic exchange. This is the first study to report the significant impact of genetic exchange on speciation in a group of viruses.

Homologous recombination as a mechanism of genetic changes in bovine parainfluenza-3 virus

Bovine parainfluenza-3 virus (BPIV-3) is one of the main viruses associated with bovine respiratory disease complex (BRDC) worldwide. BPIV-3 infect the bovine respiratory tract causing from subclinical infections to severe pneumonia with significant economic losses in the cattle industry. BPIV-3 is a RNA virus with high genetic variability, nevertheless, the contribution of recombination events to its variability has not been assessed so far. In this study the 25 complete genome sequences (CGS) reported so far and 215 partial sequences of different viral genes of BPIV-3 were analyzed to determine their genotypes and subgenotypes, distribution, and the existence of potential recombination events. Based on the analysis of the HN, M, N, and P genes one hypothetical subgenotype was found (subgenotype A4). Four recombination events between sequences of swine and cattle were detected by RDP4 analysis in conjunction with phylogenetic incongruences in the L gene. In addition, 9 sequences reported from Argentina were found to be miss-classified. These results reveal that homologous recombination events have a relevant role in the evolution of BPIV-3 and highlight the importance of implement advanced molecular characterization to better understand the variability and evolution of BPIV-3 as a component of BRDC.

New members to Arctic-like lineage of canine distemper virus from Turkey

Canine distemper virus (CDV) causes a multisystemic fatal disease, briefly named as distemper, in domestic and wild animals. Molecular characterization studies serve to identify local strains, accordingly, helps to determine the scope of vaccination in prevention of distemper. We aimed with this study to update the molecular status of CDV in domestic dogs in Turkey. Sequence analysis of the H gene revealed that novel Turkish sequences formed a separated clade in Arctic-like lineage. Italian clade which mainly included strains originated from wild canid or non-canid localized nearly to novel Turkish clade. Codons 530th and 549th determining the affinity of domestic or wild animals to distemper were Asparagine and Tyrosine, respectively. This report presented the presence of CDV strains belonging to Arctic-like lineage for the first time in domestic dogs in Turkey. The findings pave the way for the reassessment of the circulation and geographical shifting of Arctic-like lineages of CDV.

Distinct evolutionary trajectories of MHC class I and class II genes in Old World finches and buntings

Major histocompatibility complex (MHC) genes code for key proteins of the adaptive immune system, which present antigens from intra-cellular (MHC class I) and extra-cellular (MHC class II) pathogens. Because of their unprecedented diversity, MHC genes have long been an object of scientific interest, but due to methodological difficulties in genotyping of duplicated loci, our knowledge on the evolution of the MHC across different vertebrate lineages is still limited. Here, we compared the evolution of MHC class I and class II genes in three sister clades of common passerine birds, finches (Fringillinae and Carduelinae) and buntings (Emberizidae) using a uniform methodological (genotyping and data processing) approach and uniform sample sizes. Our analyses revealed contrasting evolutionary trajectories of the two MHC classes. We found a stronger signature of pervasive positive selection and higher allele diversity (allele numbers) at the MHC class I than class II. In contrast, MHC class II genes showed greater allele divergence (in terms of nucleotide diversity) and a much stronger recombination (gene conversion) signal. Gene copy numbers at both MHC class I and class II evolved via fluctuating selection and drift (Brownian Motion evolution), but the evolutionary rate was higher at class I. Our study constitutes one of few existing examples, where evolution of MHC class I and class II genes was directly compared using a multi-species approach. We recommend that re-focusing MHC research from single-species and single-class approaches towards multi-species analyses of both MHC classes can substantially increase our understanding MHC evolution in a broad phylogenetic context.

Complete genome characterization and evolutionary analysis of dengue viruses isolated during 2016-2017 in Pune, India

Dengue is the most common mosquito-borne viral infection in tropical and sub-tropical countries. In the recent years, frequent dengue outbreaks are being reported in many parts of India. DENV circulates as four independent serotypes posing a major public health threat around the globe. Phylogenetic and full genome sequence analyses of 19 complete DENV genome sequences presenting all the four serotypes in Pune, India (2016-2017) revealed no change in the circulating genotypes i.e., genotype V clade C (D1), genotype IVB (D2), genotype III lineage III (D3) and genotype I clade D (D4). Additionally, unique amino acid substitutions that may potentially influence viral fitness and virulence in host cells were identified. Mapping of the unique amino acid substitutions onto the T cell epitopes of the reference strains revealed that 8/10 (D1), 14/15 (D2), 3/4 (D3) and 21/74 (D4), amino acids were involved in T-cell epitope presentation for a maximum number of HLA alleles associated with disease outcome. Selection pressure analysis documented a positive selection pressure to be acting on few amino acid sites indicating continuous evolutionary changes in the viral RNA. Overall, the evolutionary and selection pressure data generated during this study may help in better understanding of DENV evolution and epidemiology.

Bi- and Multi-directional Gene Transfer in the Natural Populations of Polyvalent Bacteriophages, and Their Host Species Spectrum Representing Foodborne Versus Other Human and/or Animal Pathogens

Unraveling the trends of phage-host versus phage-phage coevolution is critical for avoiding possible undesirable outcomes from the use of phage preparations intended for therapeutic, food safety or environmental safety purposes. We aimed to investigate a phenomenon of intergeneric recombination and its trajectories across the natural populations of phages predominantly linked to foodborne pathogens. The results from the recombination analyses, using a large array of the recombination detection algorithms imbedded in SplitsTree, RDP4, and Simplot software packages, provided strong evidence (fit: 100; P  ≤ 0.014) for both bi- and multi-directional intergeneric recombination of the genetic loci involved collectively in phage morphogenesis, host specificity, virulence, replication, and persistence. Intergeneric recombination was determined to occur not only among conspecifics of the virulent versus temperate phages but also between the phages with these different lifestyles. The recombining polyvalent phages were suggested to interact with fairly large host species networks, including sometimes genetically very distinct species, such as e.g., Salmonella enterica and/or Escherichia coli versus Staphylococcus aureus or Yersinia pestis. Further studies are needed to understand whether phage-driven intergeneric recombination can lead to undesirable changes of intestinal and other microbiota in humans and animals.

Genetic Patterns Found in the Nuclear Localization Signals (NLSs) Associated with EBV-1 and EBV-2 Provide New Insights into Their Contribution to Different Cell-Type Specificities

Simple Summary

The Epstein–Barr virus (EBV) has been implicated in several human neoplastic diseases. The EBV-1 can transform B cells into LCL more efficiently than EBV-2, and EBV-2 preferentially infects T-cell lymphocytes. The EBNA3A oncoprotein has an essential role in B-cell transformation. The six peptide motifs called nuclear localization signals (NLSs) from EBNA3A ensure nucleocytoplasmic protein trafficking. Multiple NLSs have been suggested to enhance EBNA3 function or different specificities to different cell types; however, a comprehensive assessment of their genetic variability has not been addressed. Our objective was to study the NLSs’ variability and their relationship with EBV types. Based on a comprehensive analysis of over a thousand EBNA3A sequences from different clinical manifestations and geographic locations, we found that EBNA3A from EBV-2 has two of the six NLSs altered, and genetic patterns in the NLSs are associated with EBV-1 and EBV-2.

Abstract

The Epstein–Barr virus (EBV) is a globally dispersed pathogen involved in several human cancers of B-cell and non-B-cell origin. EBV has been classified into EBV-1 and EBV-2, which have differences in their transformative ability. EBV-1 can transform B-cells into LCL more efficiently than EBV-2, and EBV-2 preferentially infects T-cell lymphocytes. The EBNA3A oncoprotein is a transcriptional regulator of virus and host cell genes, and is required in order to transform B-cells. EBNA3A has six peptide motifs called nuclear localization signals (NLSs) that ensure nucleocytoplasmic protein trafficking. The presence of multiple NLSs has been suggested to enhance EBNA3 function or different specificities in different cell types. However, studies about the NLS variability associated with EBV types are scarce. Based on a systematic sequence analysis considering more than a thousand EBNA3A sequences of EBV from different human clinical manifestations and geographic locations, we found differences in NLSs’ nucleotide structures among EBV types. Compared with the EBNA3A EBV-1, EBNA3A EBV-2 has two of the six NLSs altered, and these mutations were possibly acquired by recombination. These genetic patterns in the NLSs associated with EBV-1 and EBV-2 provide new information about the traits of EBNA3A in EBV biology.

Genetic and evolutionary analysis of emerging HoBi-like pestivirus

HoBi-like pestivirus, an emerging species within the Pestivirus genus, is an important pathogen associated with a variety of clinical manifestations of ruminants, especially cattle. HoBi-like pestiviruses were identified in several countries and from different hosts, and raised concerns with regard to their acute and persistent infections, which is implicated in economic losses for cattle farmers. However, the transmission path, codon usage bias, and host adaptation of the virus has not been studied. Hence, we performed the analysis the spatio-temporal transmission based on the available 5'-UTR sequences of HoBi-like pestivirus, and then conducted codon analysis of the complete coding sequence of the virus. The results show the virus appeared in 1952 (95% HPD: 1905-1985) and may have originated in India. In addition, Italy is the hub for the spread of the virus. Moreover, six potential recombination events and two complex recombination events were discovered. Analysis of codon usage patterns revealed that the effective number of codon (ENC) values with an average of 50.85, and the codon usage bias is greatly affected by natural selection, which is different from the previous BVDV-1, 2. Finally, codon adaptation index (CAI) analysis shows that pigs may be the potential origin species of the HoBi-like pestivirus. These findings will contribute to more effective control of the spread of the virus, extend the knowledge about the genetic and evolutionary features of HoBi-like viruses and provide some information for vaccine research.

Tryptophan Operon Diversity Reveals Evolutionary Trends among Geographically Disparate Chlamydia trachomatis Ocular and Urogenital Strains Affecting Tryptophan Repressor and Synthase Function

The obligate intracellular pathogen Chlamydia trachomatis (Ct) is the leading cause of bacterial sexually transmitted infections and blindness globally. To date, Ct urogenital strains are considered tryptophan prototrophs, utilizing indole for tryptophan synthesis within a closed-conformation tetramer comprised of two α (TrpA)- and two β (TrpB)-subunits. In contrast, ocular strains are auxotrophs due to mutations in TrpA, relying on host tryptophan pools for survival. It has been speculated that there is strong selective pressure for urogenital strains to maintain a functional operon. Here, we performed genetic, phylogenetic, and novel functional modeling analyses of 595 geographically diverse Ct ocular, urethral, vaginal, and rectal strains with complete operon sequences. We found that ocular and urogenital, but not lymphogranuloma venereum, TrpA-coding sequences were under positive selection. However, vaginal and urethral strains exhibited greater nucleotide diversity and a higher ratio of nonsynonymous to synonymous substitutions [Pi(a)/Pi(s)] than ocular strains, suggesting a more rapid evolution of beneficial mutations. We also identified nonsynonymous amino acid changes for an ocular isolate with a urogenital backbone in the intergenic region between TrpR and TrpB at the exact binding site for YtgR-the only known iron-dependent transcription factor in Chlamydia-indicating that selective pressure has disabled the response to fluctuating iron levels. In silico effects on protein stability, ligand-binding affinity, and tryptophan repressor (TrpR) affinity for single-stranded DNA (ssDNA) measured by calculating free energy changes (ΔΔG) between Ct reference and mutant tryptophan operon proteins were also analyzed. We found that tryptophan synthase function was likely suboptimal compared to other bacterial tryptophan prototrophs and that a diversity of urogenital strain mutations rendered the synthase nonfunctional or inefficient. The novel mutations identified here affected active sites in an orthosteric manner but also hindered α- and β-subunit allosteric interactions from distant sites, reducing efficiency of the tryptophan synthase. Importantly, strains with mutant proteins were inclined toward energy conservation by exhibiting an altered affinity for their respective ligands compared to reference strains, indicating greater fitness. This is not surprising as l-tryptophan is one of the most energetically costly amino acids to synthesize. Mutations in the tryptophan repressor gene (trpR) among urogenital strains were similarly detrimental to function. Our findings indicate that urogenital strains are evolving more rapidly than previously recognized with mutations that impact tryptophan operon function in a manner that is energetically beneficial, providing a novel host-pathogen evolutionary mechanism for intracellular survival.IMPORTANCE Chlamydia trachomatis (Ct) is a major global public health concern causing sexually transmitted and ocular infections affecting over 130 million and 260 million people, respectively. Sequelae include infertility, preterm birth, ectopic pregnancy, and blindness. Ct relies on available host tryptophan pools and/or substrates to synthesize tryptophan to survive. Urogenital strains synthesize tryptophan from indole using their intact tryptophan synthase (TS). Ocular strains contain a trpA frameshift mutation that encodes a truncated TrpA with loss of TS function. We found that TS function is likely suboptimal compared to other tryptophan prototrophs and that urogenital stains contain diverse mutations that render TS nonfunctional/inefficient, evolve more rapidly than previously recognized, and impact operon function in a manner that is energetically beneficial, providing an alternative host-pathogen evolutionary mechanism for intracellular survival. Our research has broad scientific appeal since our approach can be applied to other bacteria that may explain evolution/survival in host-pathogen interactions.

KwARG: Parsimonious reconstruction of ancestral recombination graphs with recurrent mutation

MOTIVATION

The reconstruction of possible histories given a sample of genetic data in the presence of recombination and recurrent mutation is a challenging problem, but can provide key insights into the evolution of a population. We present KwARG, which implements a parsimony-based greedy heuristic algorithm for finding plausible genealogical histories (ancestral recombination graphs) that are minimal or near-minimal in the number of posited recombination and mutation events.

RESULTS

Given an input dataset of aligned sequences, KwARG outputs a list of possible candidate solutions, each comprising a list of mutation and recombination events that could have generated the dataset; the relative proportion of recombinations and recurrent mutations in a solution can be controlled via specifying a set of 'cost' parameters. We demonstrate that the algorithm performs well when compared against existing methods.

AVAILABILITY

The software is available at https://github.com/a-ignatieva/kwarg.

SUPPLEMENTARY INFORMATION

Supplementary materials are available at Bioinformatics online.

The strength of selection is consistent across both domains of the MHC class I peptide-binding groove in birds

Background

The Major Histocompatibility Complex (MHC) codes for the key vertebrate immune receptors responsible for pathogen recognition. Foreign antigens are recognized via their compatibility to hyper-variable region of the peptide-binding groove (PBR), which consists of two separate protein domains. Specifically, the PBR of the MHC class I receptors, which recognize intra-cellular pathogens, has two α domains encoded by exon 2 (α₁) and exon 3 (α₂) of the same gene. Most research on avian MHC class I polymorphism has traditionally focused exclusively on exon 3 and comparisons of selection between the two domains have been hampered by the scarcity of molecular data for exon 2. Thus, it is not clear whether the two domains vary in their specificity towards different antigens and whether they are subject to different selective pressure.

Results

Here, we took advantage of rapidly accumulating genomic resources to test for the differences in selection patterns between both MHC class I domains of the peptide-binding groove in birds. For this purpose, we compiled a dataset of MHC class I exon 2 and 3 sequences for 120 avian species from 46 families. Our phylogenetically-robust approach provided strong evidence for highly consistent levels of selection on the α₁ and α₂ domains. There were strong correlations in all selection measures (number of positively/negatively selected residues and dN/dS ratios) between both PBR exons. Similar positive associations were found for the level of amino acid polymorphism across the two domains.

Conclusions

We conclude that the strength of selection and the level of polymorphism are highly consistent between both peptide-binding domains (α₁ and α₂) of the avian MHC class I.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01812-x.

The Landscape of Recombination Events That Create Nonribosomal Peptide Diversity

Nonribosomal peptides (NRP) are crucial molecular mediators in microbial ecology and provide indispensable drugs. Nevertheless, the evolution of the flexible biosynthetic machineries that correlates with the stunning structural diversity of NRPs is poorly understood. Here, we show that recombination is a key driver in the evolution of bacterial NRP synthetase (NRPS) genes across distant bacterial phyla, which has guided structural diversification in a plethora of NRP families by extensive mixing and matching of biosynthesis genes. The systematic dissection of a large number of individual recombination events did not only unveil a striking plurality in the nature and origin of the exchange units but allowed the deduction of overarching principles that enable the efficient exchange of adenylation (A) domain substrates while keeping the functionality of the dynamic multienzyme complexes. In the majority of cases, recombination events have targeted variable portions of the Acore domains, yet domain interfaces and the flexible Asub domain remained untapped. Our results strongly contradict the widespread assumption that adenylation and condensation (C) domains coevolve and significantly challenge the attributed role of C domains as stringent selectivity filter during NRP synthesis. Moreover, they teach valuable lessons on the choice of natural exchange units in the evolution of NRPS diversity, which may guide future engineering approaches.

Identification of the Begomoviruses Squash Leaf Curl Virus and Watermelon Chlorotic Stunt Virus in Various Plant Samples in North America

Geminiviruses are a group of plant-infecting viruses with single-stranded DNA genomes. Within this family, viruses in the genus Begomovirus are known to have a worldwide distribution causing a range of severe diseases in a multitude of dicotyledonous plant species. Begomoviruses are transmitted by the whitefly Bemisia tabaci, and their ssDNA genomes can be either monopartite or bipartite. As part of a viral survey, various plants including those in the families Alliaceae, Amaranthaceae, Apiaceae, Asteraceae, Brassicaceae, Cactaceae, Cucurbitaceae, Lamiaceae, Lauraceae, Malvaceae, Oleaceae and Solanaceae were sampled and screened for begomoviruses using both a high-throughput sequencing and a begomovirus-specific primer pair approach. Based on the sequences derived using these approaches, the full-length genome of various begomoviruses were amplified from plants using abutting primers. Squash leaf curl virus (SLCV) and watermelon chlorotic stunt virus (WCSV) were identified in Cactaceae (n = 25), Solanaceae (n = 7), Cucurbitaceae (n = 2) and Lamiaceae (n = 1) samples. WCSV is an Old World bipartite begomovirus that has only recently been discovered infecting watermelons in the Americas. Our discovery of WCSV in the USA is the first indication that it has reached this country and indicates that this virus might be widespread throughout North America. Phylogenetic analysis suggests WCSV was introduced to the New World twice. The detection of begomoviruses in cactus plants suggests possible spillover events from agricultural areas into native vegetation. Since WCSV and SLCV have previously been found in mixed infections, pseudo-recombination infection experiments were conducted. We demonstrate that WCSV DNA-B is successfully trans-replicated by SLCV DNA-A despite very low degree of similarity between the replication-associated iterative sequences present in their common region, an essential feature for binding of the replication associated protein. This study highlights the importance of viral surveys for the detection of spillover events into native vegetation, but also suggests the need for more surveillance of WCSV in the USA, as this virus is a serious threat to watermelon cultivation in the Middle East.

Genome-based local dynamics of canine rabies virus epidemiology, transmission, and evolution in Davao City, Philippines, 2018-2019

Rabies is a fatal zoonotic and neglected tropical disease caused by the rabies virus (RABV) and is associated with neuronal dysfunction and death, with dogs as the predominant carrier. The Philippines plans to eradicate rabies by 2022, but this is challenged with sub-optimal coverage of vaccination programs coupled with sustained transmission chains, making it unable to eradicate the disease. We investigated the dynamics of canine rabies in the highly urbanized Davao City of the Philippines and its neighboring localities by assessing genetic relationships, transmission patterns, selection pressure, and recombination events using the whole genome sequence of 49 RABV cases from June 2018 to May 2019, majority of which (46%) were from the district of Talomo, Davao City. Although phylogeographic clustering was observed, local variants also exhibited genetic sub-lineages. Phylogenetic and spatial transmission analysis provided evidence for intra- and inter-city transmission predominantly through the Talomo district of Davao City. Around 84% of the cases were owned dogs, but the genetic similiarity of RABVs from stray and owned dogs further alluded to the role of the former as transmission vectors. The high rate of improper vaccination among the affected dogs (80%) was also a likely contributor to transmission. The RABV population under Investigation is generally under strong purifying selection with no evidence of vaccine evasion due to the genetic homogeneity of viruses from vaccinated and improperly vaccinated dogs. However, some homologous recombination (HR) events were identified along the G and L genes, also predominantly associated with viruses from Talomo. The complementary findings on epidemiology, transmission, and recombination for Talomo suggest that high incidence areas can be seeds for virus dispersal and evolution. We recommend further Investigations on the possibility of HR in future large-scale genome studies. Finally, districts associated with these phenomena can be targeted for evidence-based local strategies that can help break RABV transmission chains and prevent emergence of novel strains in Davao City.

New World Cactaceae Plants Harbor Diverse Geminiviruses

The family Cactaceae comprises a diverse group of typically succulent plants that are native to the American continent but have been introduced to nearly all other continents, predominantly for ornamental purposes. Despite their economic, cultural, and ecological importance, very little research has been conducted on the viral community that infects them. We previously identified a highly divergent geminivirus that is the first known to infect cacti. Recent research efforts in non-cultivated and asymptomatic plants have shown that the diversity of this viral family has been under-sampled. As a consequence, little is known about the effects and interactions of geminiviruses in many plants, such as cacti. With the objective to expand knowledge on the diversity of geminiviruses infecting cacti, we used previously acquired high-throughput sequencing results to search for viral sequences using BLASTx against a viral RefSeq protein database. We identified two additional sequences with similarity to geminiviruses, for which we designed abutting primers and recovered full-length genomes. From 42 cacti and five scale insects, we derived 42 complete genome sequences of a novel geminivirus species that we have tentatively named Opuntia virus 2 (OpV2) and 32 genomes of an Opuntia-infecting becurtovirus (which is a new strain of the spinach curly top Arizona virus species). Interspecies recombination analysis of the OpV2 group revealed several recombinant regions, in some cases spanning half of the genome. Phylogenetic analysis demonstrated that OpV2 is a novel geminivirus more closely related to viruses of the genus Curtovirus, which was further supported by the detection of three recombination events between curtoviruses and OpV2. Both OpV2 and Opuntia becurtoviruses were identified in mixed infections, which also included the previously characterized Opuntia virus 1. Viral quantification of the co-infected cactus plants compared with single infections did not show any clear trend in viral dynamics that might be associated with the mixed infections. Using experimental Rhizobium-mediated inoculations, we found that the initial accumulation of OpV2 is facilitated by co-infection with OpV1. This study shows that the diversity of geminiviruses that infect cacti is under-sampled and that cacti harbor diverse geminiviruses. The detection of the Opuntia becurtoviruses suggests spill-over events between viruses of cultivated species and native vegetation. The threat this poses to cacti needs to be further investigated.

Molecular Characterization of MHC Class I Genes in Four Species of the Turdidae Family to Assess Genetic Diversity and Selection

In vertebrate animals, the molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immunity. MHC class I deals with intracellular pathogens (virus) in birds. MHC class I diversity depends on the consequence of local and global environment selective pressure and gene flow. Here, we evaluated the MHC class I gene in four species of the Turdidae family from a broad geographical area of northeast China. We isolated 77 MHC class I sequences, including 47 putatively functional sequences and 30 pseudosequences from 80 individuals. Using the method based on analysis of cloned amplicons (n = 25) for each species, we found two and seven MHC I sequences per individual indicating more than one MHC I locus identified in all sampled species. Results revealed an overall elevated genetic diversity at MHC class I, evidence of different selection patterns among the domains of PBR and non-PBR. Alleles are found to be divergent with overall polymorphic sites per species ranging between 58 and 70 (out of 291 sites). Moreover, transspecies alleles were evident due to convergent evolution or recent speciation for the genus. Phylogenetic relationships among MHC I show an intermingling of alleles clustering among the Turdidae family rather than between other passerines. Pronounced MHC I gene diversity is essential for the existence of species. Our study signifies a valuable tool for the characterization of evolutionary relevant difference across a population of birds with high conservational concerns.

The Phylogeography of Potato Virus X Shows the Fingerprints of Its Human Vector

Potato virus X (PVX) occurs worldwide and causes an important potato disease. Complete PVX genomes were obtained from 326 new isolates from Peru, which is within the potato crop's main domestication center, 10 from historical PVX isolates from the Andes (Bolivia, Peru) or Europe (UK), and three from Africa (Burundi). Concatenated open reading frames (ORFs) from these genomes plus 49 published genomic sequences were analyzed. Only 18 of them were recombinants, 17 of them Peruvian. A phylogeny of the non-recombinant sequences found two major (I, II) and five minor (I-1, I-2, II-1, II-2, II-3) phylogroups, which included 12 statistically supported clusters. Analysis of 488 coat protein (CP) gene sequences, including 128 published previously, gave a completely congruent phylogeny. Among the minor phylogroups, I-2 and II-3 only contained Andean isolates, I-1 and II-2 were of both Andean and other isolates, but all of the three II-1 isolates were European. I-1, I-2, II-1 and II-2 all contained biologically typed isolates. Population genetic and dating analyses indicated that PVX emerged after potato's domestication 9000 years ago and was transported to Europe after the 15th century. Major clusters A-D probably resulted from expansions that occurred soon after the potato late-blight pandemic of the mid-19th century. Genetic comparisons of the PVX populations of different Peruvian Departments found similarities between those linked by local transport of seed potato tubers for summer rain-watered highland crops, and those linked to winter-irrigated crops in nearby coastal Departments. Comparisons also showed that, although the Andean PVX population was diverse and evolving neutrally, its spread to Europe and then elsewhere involved population expansion. PVX forms a basal Potexvirus genus lineage but its immediate progenitor is unknown. Establishing whether PVX's entirely Andean phylogroups I-2 and II-3 and its Andean recombinants threaten potato production elsewhere requires future biological studies.

Circular DNA enrichment sequencing reveals the viral/satellites genetic diversity associated with the third epidemic of cotton leaf curl disease

Cotton leaf curl disease (CLCuD) is the most important limiting factor for cotton production in Pakistan. The CLCuD passed through two major epidemics in this region with distinct begomoviruses/satellites complexes. Since 2015 the disease has again started to appear in epidemic form, causing heavy losses to cotton crop, which we termed as the “third epidemic”. We applied CIDER-seq (Circular DNA Enrichment Sequencing), a recently developed sequencing method for PCR-free virus enrichment to produce a full length read of a single circular viral genome coupled with Sanger sequencing to explore the genetic diversity of the disease complex. We identified a highly recombinant strain of Cotton leaf curl Multan virus and a recently evolved strain of Cotton leaf curl Multan betasatellite that are dominant in all major cotton growing regions in the country. Moreover, we also identified multiple species of alphasatellites with one distinct species, Mesta yellow vein mosaic alphasatellite (MeYVMA) for the first time in cotton. Relative abundance of virus and associated satellites was also determined by real-time quantitative PCR. To the best of our knowledge, this is the first study that determined the CLCuD complex associated with its third epidemic.

Recombination Events Shape the Genomic Evolution of Infectious Bronchitis Virus in Europe

Infectious bronchitis of chicken is a high morbidity and mortality viral disease affecting the poultry industry worldwide; therefore, a better understanding of this pathogen is of utmost importance. The primary aim of this study was to obtain a deeper insight into the genomic diversity of field infectious bronchitis virus (IBV) strains using phylogenetic and recombination analysis. We sequenced the genome of 20 randomly selected strains from seven European countries. After sequencing, we created a genome sequence data set that contained 36 European origin field isolates and 33 vaccine strains. When analyzing these 69 IBV genome sequences, we identified 215 recombination events highlighting that some strains had multiple recombination breaking points. Recombination hot spots were identified mostly in the regions coding for non-structural proteins, and multiple recombination hot spots were identified in the nsp2, nsp3, nsp8, and nsp12 coding regions. Recombination occurred among different IBV genotypes and involved both field and vaccine IBV strains. Ninety percent of field strains and nearly half of vaccine strains showed evidence of recombination. Despite the low number and the scattered geographical and temporal origin of whole-genome sequence data collected from European Gammacoronaviruses, this study underlines the importance of recombination as a major evolutionary mechanism of IBVs.

Recent Out-of-Africa Migration of Human Herpes Simplex Viruses

Herpes simplex virus types 1 and 2 (HSV-1 and HSV-2) are ubiquitous human pathogens. Both viruses evolved from simplex viruses infecting African primates and they are thus thought to have left Africa during early human migrations. We analyzed the population structure of HSV-1 and HSV-2 circulating strains. Results indicated that HSV-1 populations have limited geographic structure and the most evident clustering by geography is likely due to recent bottlenecks. For HSV-2, the only level of population structure is accounted for by the so-called "worldwide" and "African" lineages. Analysis of ancestry components and nucleotide diversity, however, did not support the view that the worldwide lineage followed early humans during out-of-Africa dispersal. Although phylogeographic analysis confirmed an African origin for both viruses, molecular dating with a method that corrects for the time-dependent rate phenomenon indicated that HSV-1 and HSV-2 migrated from Africa in relatively recent times. In particular, we estimated that the HSV-2 worldwide lineage left the continent in the 18th century, which corresponds to the height of the transatlantic slave trade, possibly explaining the high prevalence of HSV-2 in the Americas (second highest after Africa). The limited geographic clustering of HSV-1 makes it difficult to date its exit from Africa. The split between the basal clade, containing mostly African sequences, and all other strains was dated at ∼5,000 years ago. Our data do not imply that herpes simplex viruses did not infect early humans but show that the worldwide distribution of circulating strains is the result of relatively recent events.

Characterization of HIV-1 recombinant and subtype B near full-length genome among men who have sex with men in South Korea

In Korea, subtype B is the predominant variant of HIV-1, but full genome sequencing and analysis of its viral variants are lacking. We performed near full-length genome (NFLG) sequencing and phylogenetic and recombination analyses of fifty plasma samples from HIV-positive men who have sex with men (MSM) from a Korea HIV/AIDS cohort study. Viral genomes were amplified and the near-full-length sequences were determined using next-generation sequencing (NGS) and Sanger sequencing. We focused on the HIV-1 subtype classification and identification of HIV recombinants. Twelve HIV-1 NFLGs were determined: ten were subtyped as pure HIV-1 subtype B and two recombinant strains as a common subtype CRF07_BC, and a novel subtype CRF43_02G recombined with CRF02_AG again, or a new CRF02_AG and subtype G recombinant. For the ten NFLGs determined by NGS, “the novel recombinant emerged at approximately 2003 and the other nine subtype B about 2004 or 2005”. This is the first report analyzing HIV-1 NFLG, including recombinants and clinical characteristics, by subtype among MSM in Korea. Our results provide novel insights for understanding the recombinants in the HIV-1 epidemic in Korea.

High Diversity of Human Non-Polio Enterovirus Serotypes Identified in Contaminated Water in Nigeria

Human enteroviruses (EVs) are highly prevalent in sewage and have been associated with human diseases with complications leading to severe neurological syndromes. We have used a recently developed molecular method to investigate the presence of EVs in eight samples collected in 2017–2018 from water streams contaminated by drainage channels in three different locations in Nigeria. A total of 93 human EV strains belonging to 45 different serotypes were identified, far exceeding the number of strains and serotypes found in similar samples in previous studies. Next generation sequencing analysis retrieved whole-capsid genomic nucleotide sequences of EV strains belonging to all four A, B, C, and D species. Our results further demonstrate the value of environmental surveillance for the detection of EV transmission of both serotypes commonly associated with clinical syndromes, such as EV-A71, and those that appear to circulate silently but could eventually cause outbreaks and disease. Several uncommon serotypes, rarely reported elsewhere, were detected such as EV-A119, EV-B87, EV-C116, and EV-D111. Ten EV serotypes were detected in Nigeria for the first time and two of them, CV-A12 and EV-B86, firstly described in Africa. This method can be expanded to generate whole-genome EV sequences as we show here for one EV-D111 strain. Our data revealed phylogenetic relationships of Nigerian sewage strains with EV strains reported elsewhere, mostly from African origin, and provided new insights into the whole-genome structure of emerging serotype EV-D111 and recombination events among EV-D serotypes.

The First Whole Genome Sequence and Characterisation of Avian Nephritis Virus Genotype 3

Avian nephritis virus (ANV) is classified in the Avastroviridae family with disease associations with nephritis, uneven flock growth and runting stunting syndrome (RSS) in chicken and turkey flocks, and other avian species. The whole genome of ANV genotype 3 (ANV-3) of 6959 nucleotides including the untranslated 5′ and 3′ regions and polyadenylated tail was detected in a metagenomic virome investigation of RSS-affected chicken broiler flocks. This report characterises the ANV-3 genome, identifying partially overlapping open reading frames (ORFs), ORF1a and ORF1b, and an opposing secondary pseudoknot prior to a ribosomal frameshift stemloop structure, with a separate ORF2, whilst observing conserved astrovirus motifs. Phylogenetic analysis of the Avastroviridae whole genome and ORF2 capsid polyprotein classified the first complete whole genome of ANV-3 within Avastroviridae genogroup 2.

Phylodynamics and Codon Usage Pattern Analysis of Broad Bean Wilt Virus 2

Broad bean wilt virus 2 (BBWV-2), which belongs to the genus Fabavirus of the family Secoviridae, is an important pathogen that causes damage to broad bean, pepper, yam, spinach and other economically important ornamental and horticultural crops worldwide. Previously, only limited reports have shown the genetic variation of BBWV2. Meanwhile, the detailed evolutionary changes, synonymous codon usage bias and host adaptation of this virus are largely unclear. Here, we performed comprehensive analyses of the phylodynamics, reassortment, composition bias and codon usage pattern of BBWV2 using forty-two complete genome sequences of BBWV-2 isolates together with two other full-length RNA1 sequences and six full-length RNA2 sequences. Both recombination and reassortment had a significant influence on the genomic evolution of BBWV2. Through phylogenetic analysis we detected three and four lineages based on the ORF1 and ORF2 nonrecombinant sequences, respectively. The evolutionary rates of the two BBWV2 ORF coding sequences were 8.895 × 10⁻⁴ and 4.560 × 10⁻⁴ subs/site/year, respectively. We found a relatively conserved and stable genomic composition with a lower codon usage choice in the two BBWV2 protein coding sequences. ENC-plot and neutrality plot analyses showed that natural selection is the key factor shaping the codon usage pattern of BBWV2. Strong correlations between BBWV2 and broad bean and pepper were observed from similarity index (SiD), codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analyses. Our study is the first to evaluate the phylodynamics, codon usage patterns and adaptive evolution of a fabavirus, and our results may be useful for the understanding of the origin of this virus.

Genomic Analysis of an Indian G8P[1] Caprine Rotavirus-A Strain Revealing Artiodactyl and DS-1-Like Human Multispecies Reassortment

The surveillance studies for the presence of caprine rotavirus A (RVA) are limited in India, and the data for the whole-genome analysis of the caprine RVA is not available. This study describes the whole-genome-based analysis of a caprine rotavirus A strain, RVA/Goat-wt/IND/K-98/2015, from a goat kid in India. The genomic analysis revealed that the caprine RVA strain K-98, possess artiodactyl-like and DS-1 human-like genome constellation G8P[1]-I2-R2-C2-M2-A3-N2-T6-E2-H3. The three structural genes (VP2, VP4, and VP7) were close to caprine host having nucleotide-based identity range between 97.5 and 98.9%. Apart from them, other gene segments showed similarity with either bovine or human like genes, ultimately pointing toward a common evolutionary origin having an artiodactyl-type backbone of strain K-98. Phylogenetically, the various genes of the current study isolate also clustered inside clades comprising Human-Bovine-Caprine isolates from worldwide. The current findings add to the knowledge on caprine rotaviruses and might play a substantial role in designing future vaccines or different alternative strategies combating such infections having public health significance. To the best of our knowledge, this is the first report on the whole-genome characterization of a caprine RVA G8P[1] strain from India. Concerning the complex nature of the K-98 genome, whole-genome analyses of more numbers of RVA strains from different parts of the country are needed to comprehend the genomic nature and genetic diversity among caprine RVA.

Synonymous mutations and the molecular evolution of SARS-CoV-2 origins

Human severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is most closely related, by average genetic distance, to two coronaviruses isolated from bats, RaTG13 and RmYN02. However, there is a segment of high amino acid similarity between human SARS-CoV-2 and a pangolin-isolated strain, GD410721, in the receptor-binding domain (RBD) of the spike protein, a pattern that can be caused by either recombination or by convergent amino acid evolution driven by natural selection. We perform a detailed analysis of the synonymous divergence, which is less likely to be affected by selection than amino acid divergence, between human SARS-CoV-2 and related strains. We show that the synonymous divergence between the bat-derived viruses and SARS-CoV-2 is larger than between GD410721 and SARS-CoV-2 in the RBD, providing strong additional support for the recombination hypothesis. However, the synonymous divergence between pangolin strain and SARS-CoV-2 is also relatively high, which is not consistent with a recent recombination between them, instead, it suggests a recombination into RaTG13. We also find a 14-fold increase in the dN /dS ratio from the lineage leading to SARS-CoV-2 to the strains of the current pandemic, suggesting that the vast majority of nonsynonymous mutations currently segregating within the human strains have a negative impact on viral fitness. Finally, we estimate that the time to the most recent common ancestor of SARS-CoV-2 and RaTG13 or RmYN02 based on synonymous divergence is 51.71 years (95% CI, 28.11-75.31) and 37.02 years (95% CI, 18.19-55.85), respectively.

Potato Virus A Isolates from Three Continents: Their Biological Properties, Phylogenetics, and Prehistory

Forty-seven potato virus A (PVA) isolates from Europe, Australia, and South America's Andean region were subjected to high-throughput sequencing, and 46 complete genomes from Europe (n = 9), Australia (n = 2), and the Andes (n = 35) obtained. These and 17 other genomes gave alignments of 63 open reading frames 9,180 nucleotides long; 9 were recombinants. The nonrecombinants formed three tightly clustered, almost equidistant phylogroups; A comprised 14 Peruvian potato isolates; W comprised 37 from potato in Peru, Argentina, and elsewhere in the world; and T contained three from tamarillo in New Zealand. When five isolates were inoculated to a potato cultivar differential, three strain groups (= pathotypes) unrelated to phylogenetic groupings were recognized. No temporal signal was detected among the dated nonrecombinant sequences, but PVA and potato virus Y (PVY) are from related lineages and ecologically similar; therefore, "relative dating" was obtained using a single maximum-likelihood phylogeny of PVA and PVY sequences and PVY's well-supported 157 CE "time to most common recent ancestor". The PVA datings obtained were supported by several independent historical coincidences. The PVA and PVY populations apparently arose in the Andes approximately 18 centuries ago, and were taken to Europe during the Columbian Exchange, radiating there after the mid-19th century potato late blight pandemic. PVA's phylogroup A population diverged more recently in the Andean region, probably after new cultivars were bred locally using newly introduced Solanum tuberosum subsp. tuberosum as a parent. Such cultivars became widely grown, and apparently generated the A × W phylogroup recombinants. Phylogroup A, and its interphylogroup recombinants, might pose a biosecurity risk.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Genomic recombination events may reveal the evolution of coronavirus and the origin of SARS-CoV-2

To trace the evolution of coronaviruses and reveal the possible origin of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which causes the coronavirus disease 2019 (COVID-19), we collected and thoroughly analyzed 29,452 publicly available coronavirus genomes, including 26,312 genomes of SARS-CoV-2 strains. We observed coronavirus recombination events among different hosts including 3 independent recombination events with statistical significance between some isolates from humans, bats and pangolins. Consistent with previous records, we also detected putative recombination between strains similar or related to Bat-CoV-RaTG13 and Pangolin-CoV-2019. The putative recombination region is located inside the receptor-binding domain (RBD) of the spike glycoprotein (S protein), which may represent the origin of SARS-CoV-2. Population genetic analyses provide estimates suggesting that the putative introduced DNA within the RBD is undergoing directional evolution. This may result in the adaptation of the virus to hosts. Unsurprisingly, we found that the putative recombination region in S protein was highly diverse among strains from bats. Bats harbor numerous coronavirus subclades that frequently participate in recombination events with human coronavirus. Therefore, bats may provide a pool of genetic diversity for the origin of SARS-CoV-2.

Analysis of Whole-Genome Sequences of Infectious laryngotracheitis Virus Isolates from Poultry Flocks in Canada: Evidence of Recombination

Infectious laryngotracheitis virus (ILTV) is a herpes virus that causes an acute respiratory disease of poultry known as infectious laryngotracheitis (ILT). Chicken embryo origin (CEO) and tissue culture origin (TCO) live attenuated vaccines are routinely used for the control of ILT. However, vaccine virus is known to revert to virulence, and it has been recently shown that ILT field viral strains can undergo recombination with vaccinal ILTV and such recombinant ILT viruses possess greater transmission and pathogenicity potential. Based on complete or partial genes of the ILTV genome, few studies genotyped ILTV strains circulating in Canada, and so far, information is scarce on whole-genome sequencing or the presence of recombination in Canadian ILTV isolates. The objective of this study was to genetically characterize the 14 ILTV isolates that originated from three provinces in Canada (Alberta, British Columbia and Quebec). To this end, a phylogenetic analysis of 50 ILTV complete genome sequences, including 14 sequences of Canadian origin, was carried out. Additional phylogenetic analysis of the unique long, unique short and inverted repeat regions of the ILTV genome was also performed. We observed that 71%, 21% and 7% of the ILTV isolates were categorized as CEO revertant, wild-type and TCO vaccine-related, respectively. The sequences were also analyzed for potential recombination events, which included evidence in the British Columbia ILTV isolate. This event involved two ILTV vaccine (CEO) strains as parental strains. Recombination analysis also identified that one ILTV isolate from Alberta as a potential parental strain for a United States origin ILTV isolate. The positions of the possible recombination breakpoints were identified. These results indicate that the ILTV wild-type strains can recombine with vaccinal strains complicating vaccine-mediated control of ILT. Further studies on the pathogenicity of these ILTV strains, including the recombinant ILTV isolate are currently ongoing.

Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic

There are outstanding evolutionary questions on the recent emergence of human coronavirus SARS-CoV-2 including the role of reservoir species, the role of recombination and its time of divergence from animal viruses. We find that the sarbecoviruses-the viral subgenus containing SARS-CoV and SARS-CoV-2-undergo frequent recombination and exhibit spatially structured genetic diversity on a regional scale in China. SARS-CoV-2 itself is not a recombinant of any sarbecoviruses detected to date, and its receptor-binding motif, important for specificity to human ACE2 receptors, appears to be an ancestral trait shared with bat viruses and not one acquired recently via recombination. To employ phylogenetic dating methods, recombinant regions of a 68-genome sarbecovirus alignment were removed with three independent methods. Bayesian evolutionary rate and divergence date estimates were shown to be consistent for these three approaches and for two different prior specifications of evolutionary rates based on HCoV-OC43 and MERS-CoV. Divergence dates between SARS-CoV-2 and the bat sarbecovirus reservoir were estimated as 1948 (95% highest posterior density (HPD): 1879-1999), 1969 (95% HPD: 1930-2000) and 1982 (95% HPD: 1948-2009), indicating that the lineage giving rise to SARS-CoV-2 has been circulating unnoticed in bats for decades.

Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic

There are outstanding evolutionary questions on the recent emergence of human coronavirus SARS-CoV-2 including the role of reservoir species, the role of recombination and its time of divergence from animal viruses. We find that the sarbecoviruses-the viral subgenus containing SARS-CoV and SARS-CoV-2-undergo frequent recombination and exhibit spatially structured genetic diversity on a regional scale in China. SARS-CoV-2 itself is not a recombinant of any sarbecoviruses detected to date, and its receptor-binding motif, important for specificity to human ACE2 receptors, appears to be an ancestral trait shared with bat viruses and not one acquired recently via recombination. To employ phylogenetic dating methods, recombinant regions of a 68-genome sarbecovirus alignment were removed with three independent methods. Bayesian evolutionary rate and divergence date estimates were shown to be consistent for these three approaches and for two different prior specifications of evolutionary rates based on HCoV-OC43 and MERS-CoV. Divergence dates between SARS-CoV-2 and the bat sarbecovirus reservoir were estimated as 1948 (95% highest posterior density (HPD): 1879-1999), 1969 (95% HPD: 1930-2000) and 1982 (95% HPD: 1948-2009), indicating that the lineage giving rise to SARS-CoV-2 has been circulating unnoticed in bats for decades.

Comprehensive codon usage analysis of rice black-streaked dwarf virus based on P8 and P10 protein coding sequences

Rice black-streaked dwarf virus (RBSDV) belongs to the genus Fijivirus of the family Reoviridae and is an important pathogen that damages rice, maize and wheat worldwide. Previously, several reports have described the genetic variation and population structure of RBSDV. However, the details of the evolutionary changes, synonymous codon usage patterns and host adaptation of the virus are largely unclear. Here, we performed a detailed analysis of the codon usage and host adaptability of RBSDV based on 130 full-length P8 and 234 full-length P10 sequences. Infrequent recombination and frequent segment reassortment influence the genomic evolution of RBSDV. Our phylogenetic analysis found three and four lineages based on the P8 and P10 non-recombinant sequences respectively. We found relatively stable and conserved genomic composition with lower codon usage choice in the RBSDV P8 and P10 protein coding sequences. Both ENC-plot and neutrality-plot analyses showed that natural selection is the key factor that shapes the codon usage pattern of RBSDV. Codon adaptation index (CAI), relative codon deoptimization index (RCDI) and similarity index (SiD) analyses indicated strong correlation between RBSDV and rice rather than maize, wheat or Laodelphax striatellus. Our study provides deep insight into the evaluation of the codon usage pattern and adaptive evolution of RBSDV based on P8 and P10 sequences and should be taken into consideration for the prevention and control of this virus.

Recombination shapes African swine fever virus serotype-specific locus evolution

The recombination is one of the most frequently identified drivers of double-stranded DNA viruses evolution. However, the recombination events in African swine fever virus (ASFV) genomes have been poorly annotated. We hypothesize that the genetic determinants of ASFV variability are potential hot-spots for recombination. Here, we analyzed ASFV serotype-specific locus (C-type lectin (EP153R) and CD2v (EP402R)) in order to allocate the recombination breakpoints in these immunologically important proteins and reveal driving forces of virus evolution. The recombinations were found in both proteins, mostly among ASFV strains from East Africa, where multiple virus transmission cycles are notified. The recombination events were essentially associated with the domain organization of proteins. The phylogenetic analysis demonstrated the lack of clonal evolution for African strains which conclusively support the significance of recombinations in the serotype-specific locus. In addition, the signature of adaptive evolution of these two genes, pN/pS > 1, was demonstrated. These results have implications for the interpretation of cross-protection potential between evolutionary distant ASFV strains and strongly suggest that C-type lectin and CD2v may experience substantial selective pressure than previously thought.

Ongoing genetic evolution of H9N2 avian influenza viruses in Iranian industrial poultry farms

Despite the use of wide-scale vaccination programmes against the H9N2 virus, enzootic outbreaks of H9N2 avian influenza (AI) have often occurred and caused serious nationwide economic losses, particularly in broiler chickens. In this study, the haemagglutinin (HA) and neuraminidase (NA) genes of nine recent H9N2s and a common vaccine strain were fully sequenced and compared with other representative Iranian viruses. Phylogenetic analysis revealed that all Iranian viruses were grouped into the G1 sub-lineage with different clusters in which recent isolates (2014-2017) formed a distinct cluster compared to the vaccine group (1998-2004). All Iranian H9N2s exhibited low pathogenicity AI connecting peptide feature with an R/KSSR motif. Amino acid 226, located in the 220 loop of the receptor binding site, was leucine among the recent Iranian viruses, a characteristic of human influenza viruses. With an overall gradual increase in the genetic diversity of H9N2s, Bayesian skyline plots of Iranian HA and NA genes depicted a fluctuation and a relative stable situation, respectively. It is recommended to apply constant surveillance to assess any increase in viral human adaptation and evolutionary changes in circulating field H9N2s. Moreover, antigenic characterisation of the prevailing H9N2 viruses seems to be necessary for evaluating the possible antigenic drift from the vaccine strain.