Revisiting Recombination Signal in the Tick-Borne Encephalitis Virus: A Simulation Approach

Complete genome constellation of a dominant Bovine rotavirus genotype circulating in Bangladesh reveals NSP4 intragenic recombination with human strains

Rotavirus A is a leading cause of non-bacterial gastroenteritis in humans and domesticated animals. Despite the vast diversity of bovine Rotavirus A strains documented in South Asian countries, there are very few whole genomes available for phylogenetic study. A cross-sectional study identified a high prevalence of the G6P[11] genotype of bovine Rotavirus A circulating in the commercial cattle population in Bangladesh. Next-generation sequencing and downstream phylogenetic analysis unveiled all 11 complete gene segments of this strain (BD_ROTA_CVASU), classifying it under the genomic constellation G6P[11]-I2-R2-C2-M2-A13-N2-T6-E2-H3, which belongs to a classical DS-1-like genomic backbone. We found strong evidence of intragenic recombination between human and bovine strains in the Non-structural protein 4 (NSP4) gene, which encodes a multifunctional enterotoxin. Our analyses highlight frequent zoonotic transmissions of rotaviruses in diverse human-animal interfaces, which might have contributed to the evolution and pathogenesis of this dominant genotype circulating in the commercial cattle population in Bangladesh.

More than mcr: canonical plasmid- and transposon-encoded mobilized colistin resistance genes represent a subset of phosphoethanolamine transferases

Mobilized colistin resistance genes (mcr) may confer resistance to the last-resort antimicrobial colistin and can often be transmitted horizontally. mcr encode phosphoethanolamine transferases (PET), which are closely related to chromosomally encoded, intrinsic lipid modification PET (i-PET; e.g., EptA, EptB, CptA). To gain insight into the evolution of mcr within the context of i-PET, we identified 69,814 MCR-like proteins present across 256 bacterial genera (obtained by querying known MCR family representatives against the National Center for Biotechnology Information [NCBI] non-redundant protein database via protein BLAST). We subsequently identified 125 putative novel mcr-like genes, which were located on the same contig as (i) ≥1 plasmid replicon and (ii) ≥1 additional antimicrobial resistance gene (obtained by querying the PlasmidFinder database and NCBI's National Database of Antibiotic Resistant Organisms, respectively, via nucleotide BLAST). At 80% amino acid identity, these putative novel MCR-like proteins formed 13 clusters, five of which represented putative novel MCR families. Sequence similarity and a maximum likelihood phylogeny of mcr, putative novel mcr-like, and ipet genes indicated that sequence similarity was insufficient to discriminate mcr from ipet genes. A mixed-effect model of evolution (MEME) indicated that site- and branch-specific positive selection played a role in the evolution of alleles within the mcr-2 and mcr-9 families. MEME suggested that positive selection played a role in the diversification of several residues in structurally important regions, including (i) a bridging region that connects the membrane-bound and catalytic periplasmic domains, and (ii) a periplasmic loop juxtaposing the substrate entry tunnel. Moreover, eptA and mcr were localized within different genomic contexts. Canonical eptA genes were typically chromosomally encoded in an operon with a two-component regulatory system or adjacent to a TetR-type regulator. Conversely, mcr were represented by single-gene operons or adjacent to pap2 and dgkA, which encode a PAP2 family lipid A phosphatase and diacylglycerol kinase, respectively. Our data suggest that eptA can give rise to "colistin resistance genes" through various mechanisms, including mobilization, selection, and diversification of genomic context and regulatory pathways. These mechanisms likely altered gene expression levels and enzyme activity, allowing bona fide eptA to evolve to function in colistin resistance.

The Emergence and Dynamics of Tick-Borne Encephalitis Virus in a New Endemic Region in Southern Germany

Tick-borne encephalitis (TBE) is the most important viral tick-borne infection in Europe and Asia. It is emerging in new areas. The mechanisms of emergence are fairly unknown or speculative. In the Ravensburg district in southern Germany, TBE emerged, mainly over the last five years. Here, we analyzed the underlying epidemiology in humans. The resulting identified natural foci of the causal TBE virus (TBEV) were genetically characterized. We sampled 13 potential infection sites at these foci and detected TBEV in ticks (Ixodes ricinus) at eight sites. Phylogenetic analysis spurred the introduction of at least four distinct TBEV lineages of the European subtype into the Ravensburg district over the last few years. In two instances, a continuous spread of these virus strains over up to 10 km was observed.

Current Methods for Recombination Detection in Bacteria

The role of genetic exchanges, i.e., homologous recombination (HR) and horizontal gene transfer (HGT), in bacteria cannot be overestimated for it is a pivotal mechanism leading to their evolution and adaptation, thus, tracking the signs of recombination and HGT events is importance both for fundamental and applied science. To date, dozens of bioinformatics tools for revealing recombination signals are available, however, their pros and cons as well as the spectra of solvable tasks have not yet been systematically reviewed. Moreover, there are two major groups of software. One aims to infer evidence of HR, while the other only deals with horizontal gene transfer (HGT). However, despite seemingly different goals, all the methods use similar algorithmic approaches, and the processes are interconnected in terms of genomic evolution influencing each other. In this review, we propose a classification of novel instruments for both HR and HGT detection based on the genomic consequences of recombination. In this context, we summarize available methodologies paying particular attention to the type of traceable events for which a certain program has been designed.

LOUPING-ILL VIRUS SEROSURVEY OF WILLOW PTARMIGAN (LAGOPUS LAGOPUS LAGOPUS) IN NORWAY

In Norway, the Willow Ptarmigan (Lagopus lagopus lagopus) is experiencing population declines and is nationally Red Listed as Near Threatened. Although disease has not generally been regarded as an important factor behind population fluctuations for Willow Ptarmigan in Norway, disease occurrence has been poorly investigated. Both louping-ill virus (LIV) and the closely related tick-borne encephalitis virus are found along the southern part of the Norwegian coast. We assessed whether and where Norwegian Willow Ptarmigan populations have been infected with LIV. We expected to find infected individuals in populations in the southernmost part of the country. We did not expect to find infected individuals in populations further north and at higher altitudes because of the absence of the main vector, the sheep tick (Ixodes ricinus). We collected serum samples on Nobuto filter paper and used a hemagglutination inhibition assay for antibodies against LIV. We collected data at both local and country-wide levels. For local sampling, we collected and analyzed 87 hunter-collected samples from one of the southernmost Willow Ptarmigan populations in Norway. Of these birds, only three positives (3.4%) were found. For the country-wide sampling, we collected serum samples from 163 Willow Ptarmigan carcasses submitted from selected locations all over the country. Of these birds, 32% (53) were seropositive for LIV or a cross-reacting virus. Surprisingly, we found seropositive individuals from locations across the whole country, including outside the known distribution of the sheep tick. These results suggest that either LIV or a cross-reacting virus infects ptarmigan in large parts of Norway, including at high altitudes and latitudes.

Mitochondrial DNA in Ixodus ricinus (Acari: Ixodidae) on birds reflects ticks' transportation routes to Lista, Norway

Ticks are important pathogen vectors, and large mammals and birds have the greatest potential for dispersing them. To study tick dispersal by migrating birds, we have analysed genetic variations in mitochondrial DNA control region from Ixodes ricinus from northward migrating blackbird, Turdus merula, and (European) robin, Erithacus rubecula, at the Lista Bird Observatory in southwestern Norway. We compared their genetic structure with that of resident tick populations from areas covering their expected last stop (i.e. Great Britain and Jutland, Denmark) before taking off for southern Norway, and the resident tick population at Lista. The statistical analysis showed that the I. ricinus found on blackbirds differed significantly from those found on robins, which is consistent with the birds' differential migration routes. I. ricinus from robins did not differ genetically from those flagged at Jutland, suggesting that the former mainly originate in continental Europe. Bayesian analysis indicated that most of the blackbirds caught early in the spring (i.e. before or on the 1st of April) carried ticks of a mixed origin from both Great Britain and continental Europe, while blackbirds caught later in the season carried an increasing amount of ticks acquired locally.

Reporting of New tick-borne encephalitis virus strains isolated in Eastern Siberia (Russia) in 1960-2011 and explaining them in an evolutionary context using Bayesian phylogenetic inference

Tick-borne encephalitis virus (TBEV) is one of the main tick-borne viral pathogens of humans. Infection may induce signs of meningitis, encephalitis, paralysis and high fever. TBEV is well studied by molecular phylogenetic methods. The present-day implementation of Bayesian phylogenetic models allows population dynamics to be tracked, providing changes in population size that were not directly observed. However, the description of the past population dynamics of TBEV is rare in the literature. In our investigation, we provide data on the dynamics of viral genetic diversity of TBEV in Zabaikalsky Krai (Eastern Siberia, Russia) revealed by the Bayesian coalescent inference in a BEAST program. As a data set, we used the envelope (E) protein partial gene sequences (1308 nt) of 38 TBEV strains (including six "886-84-like" or Baikalian subtype strains (TBEV-B)), isolated in Zabaikalsky Krai (Eastern Siberia, Russia) in 1960-1963 and 1995-2011. To increase estimations reliability, we compared 9 model combinations by Path sampling and Stepping-stone sampling methods. It has been shown that the genetic diversity decline in the population history of TBEV in the 1950s coincides with the date of the beginning of wide dichlorodiphenyltrichloroethane forest dusting in Siberia. We assumed that the TBEV population on the territory of Siberia went through a genetic bottleneck. Also, we provide data estimating the divergence time of TBEV-B strains and indicate the specific evolution rate of an ancestor lineage of the Baikalian subtype, illustrated on a phylogenetic tree, and reconstructed under a relaxed clock model.

Genomic Analyses of Human Sapoviruses Detected over a 40-Year Period Reveal Disparate Patterns of Evolution among Genotypes and Genome Regions

Human sapovirus is a causative agent of acute gastroenteritis in all age groups. The use of full-length viral genomes has proven beneficial to investigate evolutionary dynamics and transmission chains. In this study, we developed a full-length genome sequencing platform for human sapovirus and sequenced the oldest available strains (collected in the 1970s) to analyse diversification of sapoviruses. Sequence analyses from five major genotypes (GI.1, GI.2, GII.1, GII.3, and GIV.1) showed limited intra-genotypic diversification for over 20–40 years. The accumulation of amino acid mutations in VP1 was detected for GI.2 and GIV.1 viruses, while having a similar rate of nucleotide evolution to the other genotypes. Differences in the phylogenetic clustering were detected between RdRp and VP1 sequences of our archival strains as well as other reported putative recombinants. However, the lack of the parental strains and differences in diversification among genomic regions suggest that discrepancies in the phylogenetic clustering of sapoviruses could be explained, not only by recombination, but also by disparate nucleotide substitution patterns between RdRp and VP1 sequences. Together, this study shows that, contrary to noroviruses, sapoviruses present limited diversification by means of intra-genotype variation and recombination.

RDP5: a computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets

For the past 20 years, the recombination detection program (RDP) project has focused on the development of a fast, flexible, and easy to use Windows-based recombination analysis tool. Whereas previous versions of this tool have relied on considerable user-mediated verification of detected recombination events, the latest iteration, RDP5, is automated enough that it can be integrated within analysis pipelines and run without any user input. The main innovation enabling this degree of automation is the implementation of statistical tests to identify recombination signals that could be attributable to evolutionary processes other than recombination. The additional analysis time required for these tests has been offset by algorithmic improvements throughout the program such that, relative to RDP4, RDP5 will still run up to five times faster and be capable of analyzing alignments containing twice as many sequences (up to 5000) that are five times longer (up to 50 million sites). For users wanting to remove signals of recombination from their datasets before using them for downstream phylogenetics-based molecular evolution analyses, RDP5 can disassemble detected recombinant sequences into their constituent parts and output a variety of different recombination-free datasets in an array of different alignment formats. For users that are interested in exploring the recombination history of their datasets, all the manual verification, data management and data visualization components of RDP5 have been extensively updated to minimize the amount of time needed by users to individually verify and refine the program's interpretation of each of the individual recombination events that it detects.

Intragenic recombination influences rotavirus diversity and evolution

Because of their replication mode and segmented dsRNA genome, homologous recombination is assumed to be rare in the rotaviruses. We analyzed 23,627 complete rotavirus genome sequences available in the NCBI Virus Variation database, and found 109 instances of homologous recombination, at least eleven of which prevailed across multiple sequenced isolates. In one case, recombination may have generated a novel rotavirus VP1 lineage. We also found strong evidence for intergenotypic recombination in which more than one sequence strongly supported the same event, particularly between different genotypes of segment 9, which encodes the glycoprotein, VP7. The recombined regions of many putative recombinants showed amino acid substitutions differentiating them from their major and minor parents. This finding suggests that these recombination events were not overly deleterious, since presumably these recombinants proliferated long enough to acquire adaptive mutations in their recombined regions. Protein structural predictions indicated that, despite the sometimes substantial amino acid replacements resulting from recombination, the overall protein structures remained relatively unaffected. Notably, recombination junctions appear to occur nonrandomly with hot spots corresponding to secondary RNA structures, a pattern seen consistently across segments. In total, we found strong evidence for recombination in nine of eleven rotavirus A segments. Only segments 7 (NSP3) and 11 (NSP5) did not show strong evidence of recombination. Collectively, the results of our computational analyses suggest that, contrary to the prevailing sentiment, recombination may be a significant driver of rotavirus evolution and may influence circulating strain diversity.

Genetic diversity and geographical distribution of the Siberian subtype of the tick-borne encephalitis virus

The tick-borne encephalitis virus (TBEV), a member of the Flaviviridae family, is currently subdivided into three main subtypes-the European (TBEV-Eu), the Far-Eastern (TBEV-FE), and the Siberian (TBEV-Sib). The TBEV-Sib is the most common subtype and found in all regions where TBEV was detected, except for Central and Western Europe. Currently, four genetic lineages have been described within TBEV-Sib. In this study, detailed analysis of TBEV-Sib genetic diversity, geographic distribution, phylogeography and divergence time of different TBEV-Sib genetic lineages based on E gene fragments, complete genome sequences, and all currently available data in the GenBank database was performed. As a result, a novel Bosnia lineage within the TBEV-Sib was identified. It was demonstrated that the Zausaev lineage is the most widely distributed among the TBEV-Sib lineages, and was detected in all studied regions except the Far East. The Vasilchenko lineage was found from Western Siberia to the Far East. The Baltic lineage is presented from Europe to Western Siberia. The Obskaya lineage was found only in Western Siberia. TBEV strains from a newly described Bosnia lineage were detected in Bosnia, the Crimean peninsula, Kyrgyzstan and Kazakhstan. The greatest divergence of the TBEV-Sib genetic variants was observed in Western Siberia. Within the TBEV-Sib, the Obskaya lineage diverged from the common ancestor the earliest, after that the Bosnia lineage was separated, then the Baltic lineage, and the Zausaev and Vasilchenko lineages diverged most recently.

Variability in the 3' untranslated regions of the genomes of the different tick-borne encephalitis virus subtypes

Tick-borne encephalitis viruses (TBEVs) are usually divided into three major subtypes: European (TBEV-Eu), Siberian (TBEV-Sib) and Far Eastern (TBEV-FE). The TBEV-Eu strains have the longest genomes, and TBEV-FE strains have the smallest genomes. Changes in the variable region of the untranslated region (V3' UTR) play a major role in determining the viral genome length. Analyses of the 3' UTRs of the different subtypes of TBEV have revealed significant changes in the secondary structures of the V3' UTR of TBEV. More complex secondary structures of the V3' UTR regions are typical for TBEV-Eu. The Siberian strain Tomsk-PT122 was isolated from birds and has an unusual 3' UTR. Several short fragment (24-26 nucleotides) insertions derived from the viral E (2) and NS4a (1) genes have been found in the V3' UTR of Tomsk-PT122. Additionally, the length of the V3' UTR increases from 21 to 37 nucleotides during passages of the C11-13 strain of TBEV-Sib into PEK, 293 and Neuro-2a cells. The elongation of the V3' UTRs of Tomsk-PT122 and C11-13 is the first direct evidence of an intragenomic 3' UTR modification (insertion) for TBEV. Thus, the obtained results suggest that changing the length of the V3' UTR in the genome is typical for different TBEV subtypes and can play an essential role in effective TBEV replication in different host cells.

A case series of fatal meningoencephalitis in Mongolia: epidemiological and molecular characteristics of tick-borne encephalitis virus

In Mongolia, the incidence and fatality rates of tick-borne encephalitis (TBE) have been increasing. We aimed to identify the epidemiological and molecular characteristics of tick-borne encephalitis virus (TBEV) associated with fatal meningoencephalitis in Mongolia. We conducted a descriptive study of 14 fatal cases of TBE that occurred between 2008 and 2017 in Mongolia. Reverse transcription polymerase chain reaction (RT-PCR) was used to detect viral RNA in brain tissue. RT-PCR products from six patients who died from TBE between 2013 and 2017 were directly sequenced and analysed phylogenetically. Ticks collected from Selenge and Bulgan provinces were also tested for TBEV by RT-PCR. Between 2008 and 2017, there were 14 fatal TBE cases in hospitals in Mongolia. The 14 patients who died reported receiving tick bites in Bulgan or Selenge province; 71.4% of deaths resulted from tick bites in Bulgan province. The TBE case fatality rate was 28.6% for patients in Bulgan province and 2.7% for those in Selenge province. All of the fatalities were men; the median age was 45 ± 12.6 years. Tick bites occurred between April and June in forested areas. In 2013, a 388 base pair fragment of the envelope (E) gene was obtained from a hospitalized patient. The closest relatives of this virus are Far-Eastern TBEV isolates. The case fatality rate differed between two provinces where tick bites occurred. A higher number of TBE cases and the virulent Far-Eastern subtype occurred in patients in Bulgan province. This province should increase vaccination coverage, training, education and investigations.

Mechanisms and consequences of positive-strand RNA virus recombination

Genetic recombination in positive-strand RNA viruses is a significant evolutionary mechanism that drives the creation of viral diversity by the formation of novel chimaeric genomes. The process and its consequences, for example the generation of viruses with novel phenotypes, has historically been studied by analysis of the end products. More recently, with an appreciation that there are both replicative and non-replicative mechanisms at work, and with new approaches and techniques to analyse intermediate products, the viral and cellular factors that influence the process are becoming understood. The major influence on replicative recombination is the fidelity of viral polymerase, although RNA structures and sequences may also have an impact. In replicative recombination the viral polymerase is necessary and sufficient, although roles for other viral or cellular proteins may exist. In contrast, non-replicative recombination appears to be mediated solely by cellular components. Despite these insights, the relative importance of replicative and non-replicative mechanisms is not clear. Using single-stranded positive-sense RNA viruses as exemplars, we review the current state of understanding of the processes and consequences of recombination.