Evaluation of methods for detecting recombination from DNA sequences: computer simulations

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.

Comparative genomics revealed adaptive admixture in Cryptosporidium hominis in Africa

Cryptosporidiosis is a major cause of diarrhoeal illness among African children, and is associated with childhood mortality, malnutrition, cognitive development and growth retardation. Cryptosporidium hominis is the dominant pathogen in Africa, and genotyping at the glycoprotein 60 (gp60) gene has revealed a complex distribution of different subtypes across this continent. However, a comprehensive exploration of the metapopulation structure and evolution based on whole-genome data has yet to be performed. Here, we sequenced and analysed the genomes of 26 C. hominis isolates, representing different gp60 subtypes, collected at rural sites in Gabon, Ghana, Madagascar and Tanzania. Phylogenetic and cluster analyses based on single-nucleotide polymorphisms showed that isolates predominantly clustered by their country of origin, irrespective of their gp60 subtype. We found a significant isolation-by-distance signature that shows the importance of local transmission, but we also detected evidence of hybridization between isolates of different geographical regions. We identified 37 outlier genes with exceptionally high nucleotide diversity, and this group is significantly enriched for genes encoding extracellular proteins and signal peptides. Furthermore, these genes are found more often than expected in recombinant regions, and they show a distinct signature of positive or balancing selection. We conclude that: (1) the metapopulation structure of C. hominis can only be accurately captured by whole-genome analyses; (2) local anthroponotic transmission underpins the spread of this pathogen in Africa; (3) hybridization occurs between distinct geographical lineages; and (4) genetic introgression provides novel substrate for positive or balancing selection in genes involved in host-parasite coevolution.

Analysis of the beak and feather disease viral genome indicates evidence of multiple introduction events into Saudi Arabia

Psittacine beak and feather disease (PBFD), caused by beak and feather disease virus (BFDV) is a highly contagious disease in wild and captive psittacine populations and has an almost global presence. However, the BFDV infection in Saudi Arabia remains largely unknown. In the present study, we report the full genome sequence of BFDV strains from Saudi Arabia and its genetic diversity. The complete genome sequences were analyzed for 14 BFDV-infected birds representing 6 psittacine species. The complete genome sequence of BFDV strains was compared with 201 previously reported sequences to evaluate their diversity and possible recombination events, if any. Our analysis revealed that newly sequenced BFDV genomes from Saudi Arabia belonged to six different strains. Phylogenetic analysis suggested that the isolated BFDV genomes were highly recombinant with a high degree of diversity. It is evident from the study that psittacine species in Saudi Arabia are at risk from the spread of BFDV. As per the CITES trade database, about 190,000 parrots have been imported to Saudi Arabia since 1975 over a thousand instances. Presumably, during any of these trade events or unregulated trade of birds has predisposed the introduction of BFDV to Saudi Arabia. Understanding the epidemiology of BFDV is necessitated to address the threat posed by the virus to the psittacine population of Saudi Arabia.

Genetic Diversity of Potyviruses Associated with Tulip Breaking Syndrome

Tulip breaking is economically the most important viral disease of modern-day tulip growing. It is characterized by irregular flame and feather-like patterns in the flowers and mosaic on the foliage. Thirty-two leaf samples were collected from cultivated tulip plants showing tulip breaking syndrome from Hungary in 2017 and 2018. Virus identification was performed by serological (ELISA) and molecular (RT-PCR) methods. All samples proved to be infected with a potyvirus and evidence was provided that three potyvirus species could be identified in the samples: Lily mottle virus (LMoV), Tulip breaking virus (TBV) and Rembrandt tulip-breaking virus (ReTBV). Recombination prediction accomplished with Recombination Detection Program (RDP) v4.98 revealed potential intraspecies recombination in the case of TBV and LMoV. Phylogenetic analyses of the coat protein (CP) regions proved the monophyletic origin of these viruses and verified them as three different species according to current International Committee on Taxonomy of Viruses (ICTV) species demarcation criteria. Based on these results, we analyzed taxonomic relations concerning potyviruses associated with tulip breaking syndrome. We propose the elevation of ReTBV to species level, and emergence of two new subgroups in ReTBV.

Genomic signatures of recombination in a natural population of the bdelloid rotifer Adineta vaga

Sexual reproduction is almost ubiquitous among extant eukaryotes. As most asexual lineages are short-lived, abandoning sex is commonly regarded as an evolutionary dead end. Still, putative anciently asexual lineages challenge this view. One of the most striking examples are bdelloid rotifers, microscopic freshwater invertebrates believed to have completely abandoned sexual reproduction tens of Myr ago. Here, we compare whole genomes of 11 wild-caught individuals of the bdelloid rotifer Adineta vaga and present evidence that some patterns in its genetic variation are incompatible with strict clonality and lack of genetic exchange. These patterns include genotype proportions close to Hardy-Weinberg expectations within loci, lack of linkage disequilibrium between distant loci, incongruent haplotype phylogenies across the genome, and evidence for hybridization between divergent lineages. Analysis of triallelic sites independently corroborates these findings. Our results provide evidence for interindividual genetic exchange and recombination in A. vaga, a species previously thought to be anciently asexual.

Several lineages of porcine epidemic diarrhea virus in Colombia during the 2014 and 2016 epidemic

Porcine epidemic diarrhea virus (PEDV) is a significant global, enteric coronavirus in pigs and was first reported in Colombia in 2014. However, the epidemiology, genetic and antigenic characteristics of the virus have yet to be investigated. In this study, we investigated the dissemination of PEDV by testing 536 samples by RT-PCR over a 33-month period. The 35.8% of positive samples (n = 192) was significantly different (p < .01) between months over time, with a higher number of positives samples occurring at the beginning of the epidemic and during the second epidemic wave within the main pork producing region. The complete PEDV genomes were generated for 21 strains, which shared a high nucleotide and amino acid sequence identity, except for the spike (S) gene. Recombinant regions were identified within the Colombian strains and between Colombian and Asian PEDV strains. Phylogenetic analysis of the 21 Colombian strains demonstrated the presence of 7 lineages that shared common ancestors with PEDV strains from the United States. Moreover, the antigenic analysis demonstrated residue differences in the neutralizing epitopes in the spike and nucleocapsid proteins. Our results illustrated the simultaneous introduction of the two PEDV genotypes (GIIa American pandemic and S-INDEL) into the Colombian swine industry during the 2014 PEDV epidemic and enhanced our understanding of the epidemiology and molecular diversity of PEDV in Colombia.

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.

Bovine Coronavirus: Variability, Evolution, and Dispersal Patterns of a No Longer Neglected Betacoronavirus

Bovine coronavirus (BoCV) is an important pathogen of cattle, causing severe enteric disease and playing a role in the bovine respiratory disease complex. Similar to other coronaviruses, a remarkable variability characterizes both its genome and biology. Despite their potential relevance, different aspects of the evolution of BoCV remain elusive. The present study reconstructs the history and evolution of BoCV using a phylodynamic approach based on complete genome and spike protein sequences. The results demonstrate high mutation and recombination rates affecting different parts of the viral genome. In the spike gene, this variability undergoes significant selective pressures—particularly episodic pressure—located mainly on the protein surface, suggesting an immune-induced selective pressure. The occurrence of compensatory mutations was also identified. On the contrary, no strong evidence in favor of host and/or tissue tropism affecting viral evolution has been proven. The well-known plasticity is thus ascribable to the innate broad viral tropism rather than mid- or long-term adaptation. The evaluation of the geographic spreading pattern clearly evidenced two clusters: a European cluster and an American–Asian cluster. While a relatively dense and quick migration network was identified in the former, the latter was dominated by the primary role of the United States (US) as a viral exportation source. Since the viral spreading pattern strongly mirrored the cattle trade, the need for more intense monitoring and preventive measures cannot be underestimated as well as the need to enforce the vaccination of young animals before international trade, to reduce not only the clinical impact but also the transferal and mixing of BoCV strains.

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.

Identifying Candidate Genetic Markers of CDV Cross-Species Pathogenicity in African Lions

Canine distemper virus (CDV) is a multi-host pathogen with variable clinical outcomes of infection across and within species. We used whole-genome sequencing (WGS) to search for viral markers correlated with clinical distemper in African lions. To identify candidate markers, we first documented single-nucleotide polymorphisms (SNPs) differentiating CDV strains associated with different clinical outcomes in lions in East Africa. We then conducted evolutionary analyses on WGS from all global CDV lineages to identify loci subject to selection. SNPs that both differentiated East African strains and were under selection were mapped to a phylogenetic tree representing global CDV diversity to assess if candidate markers correlated with documented outbreaks of clinical distemper in lions (n = 3). Of 54 SNPs differentiating East African strains, ten were under positive or episodic diversifying selection and 20 occurred in the clinical strain despite strong purifying selection at those loci. Candidate markers were in functional domains of the RNP complex (n = 19), the matrix protein (n = 4), on CDV glycoproteins (n = 5), and on the V protein (n = 1). We found mutations at two loci in common between sequences from three CDV outbreaks of clinical distemper in African lions; one in the signaling lymphocytic activation molecule receptor (SLAM)-binding region of the hemagglutinin protein and another in the catalytic center of phosphodiester bond formation on the large polymerase protein. These results suggest convergent evolution at these sites may have a functional role in clinical distemper outbreaks in African lions and uncover potential novel barriers to pathogenicity in this species.

You Will Never Walk Alone: Codispersal of JC Polyomavirus with Human Populations

JC polyomavirus (JCPyV) is one of the most prevalent human viruses. Findings based on the geographic distribution of viral subtypes suggested that JCPyV codiverged with human populations. This view was however challenged by data reporting a much more recent origin and expansion of JCPyV. We collected information on ∼1,100 worldwide strains and we show that their geographic distribution roughly corresponds to major human migratory routes. Bayesian phylogeographic analysis inferred a Subsaharan origin for JCPyV, although with low posterior probability. High confidence inference at internal nodes provided strong support for a long-standing association between the virus and human populations. In line with these data, pairwise FST values for JCPyV and human mtDNA sampled from the same areas showed a positive and significant correlation. Likewise, very strong relationships were found when node ages in the JCPyV phylogeny were correlated with human population genetic distances (nuclear-marker based FST). Reconciliation analysis detected a significant cophylogenetic signal for the human population and JCPyV trees. Notably, JCPyV also traced some relatively recent migration events such as the expansion of people from the Philippines/Taiwan area into Remote Oceania, the gene flow between North-Eastern Siberian and Ainus, and the Koryak contribution to Circum-Arctic Americans. Finally, different molecular dating approaches dated the origin of JCPyV in a time frame that precedes human out-of-Africa migration. Thus, JCPyV infected early human populations and accompanied our species during worldwide dispersal. JCPyV typing can provide reliable geographic information and the virus most likely adapted to the genetic background of human populations.

Emergence and Spread of Piscine orthoreovirus Genotype 3

Piscine orthoreovirus (PRV) is a relevant pathogen for salmonid aquaculture worldwide. In 2015, a new genotype of PRV (genotype 3, PRV-3) was discovered in Norway, and in 2017 PRV-3 was detected for first time in Denmark in association with complex disease cases in rainbow trout in recirculating aquaculture systems (RAS). To explore the epidemiology of PRV-3 in Denmark, a surveillance study was conducted in 2017 to 2019. Fifty-three farms, including both flow through and RAS, were screened for PRV-3. Of the farms examined, PRV-3 was detected in thirty-eight (71.7%), with the highest prevalence in grow-out farms. Notably, in Denmark disease outbreaks were only observed in RAS. Additionally, wild Atlantic salmon and brown trout populations were included in the screening, and PRV-3 was not detected in the three years where samples were obtained (2016, 2018, and 2019). Historical samples in the form of archived material at the Danish National Reference Laboratory for Fish Diseases were also tested for the presence of PRV-3, allowing us to establish that the virus has been present in Denmark at least since 1995. Sequence analyses of segment S1 and M2, as well as full genome analyses of selected isolates, did not reveal clear association between genetic makeup in these two segments and virulence in the form of disease outbreaks in the field.

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.

Chicken Astrovirus (CAstV) Molecular Studies Reveal Evidence of Multiple Past Recombination Events in Sequences Originated from Clinical Samples of White Chick Syndrome (WCS) in Western Canada

In this study, we aimed to molecularly characterize 14 whole genome sequences of chicken astrovirus (CAstV) isolated from samples obtained from white chick syndrome (WCS) outbreaks in Western Canada during the period of 2014–2019. Genome sequence comparisons showed all these sequences correspond to the novel Biv group from which no confirmed representatives were published in GenBank. Molecular recombination analyses using recombination detection software (i.e., RDP5 and SimPlot) and phylogenetic analyses suggest multiple past recombination events in open reading frame (ORF)1a, ORF1b, and ORF2. Our findings suggest that recombination events and the accumulation of point mutations may have contributed to the substantial genetic variation observed in CAstV and evidenced by the current seven antigenic sub-clusters hitherto described. This is the first paper that describes recombination events in CAstV following analysis of complete CAstV sequences originated in Canada.

Genotypic characterization and molecular evolution of avian reovirus in poultry flocks from Brazil

Avian reovirus (ARV) is one of the main causes of infectious arthritis/tenosynovitis and malabsorption syndrome (MAS) in poultry. ARVs have been disseminated in Brazilian poultry flocks in the last years. This study aimed to genotype ARVs and to evaluate the molecular evolution of the more frequent ARV lineages detected in Brazilian poultry-producing farms. A total of 100 poultry flocks with clinical signs of tenosynovitis/MAS, from all Brazilian poultry-producing regions were positive for ARV by PCR. Seventeen bird tissues were submitted to cell culture and ARV RNA detection/genotyping by two PCRs. The phylogenetic classification was based on σC gene alignment using a dataset with other Brazilian and worldwide ARVs sequences. ARVs were specifically detected by both PCRs from the 17 cell cultures, and σC gene partial fragments were sequenced. All these sequences were aligned with a total of 451 ARV σC gene data available in GenBank. Phylogenetic analysis demonstrated five well-defined clusters that were classified into lineages I, II, III, IV, and V. Three lineages could be further divided into sub-lineages: I (I vaccine, Ia, Ib), II (IIa, IIb, IIc) and IV (IVa and IVb). Brazilian ARVs were from four lineages/sub-lineages: Ib (48.2%), IIb (22.2%), III (3.7%) and V (25.9%). The Bayesian analysis demonstrated that the most frequent sub-lineage Ib emerged in the world around 1968 and it was introduced into Brazil in 2010, with increasing spread soon after. In conclusion, four different ARV lineages are circulating in Brazilian poultry flocks, all associated with clinical diseases. RESEARCH HIGHLIGHTS One-hundred ARV-positive flocks were detected in all main poultry-producing regions from Brazil. A large dataset of 468 S1 sequences was constructed and divided ARVs into five lineages. Four lineages/sub-lineages (Ib, IIb, III and V) were detected in commercial poultry flocks from Brazil. Brazilian lineages shared a low identity with the commercial vaccine lineage (I vaccine). Sub-lineage Ib emerged around 1968 and was introduced into Brazil in 2010.

DNA sequence reconstruction based on innovated hybridization technique of probabilistic cellular automata and particle swarm optimization

DNA sequence reconstruction is a challenging research problem in the computational biology field. The evolution of the DNA is too complex to be characterized by a few parameters. Therefore, there is a need for a modeling approach for analyzing DNA patterns. In this paper, we proposed a novel framework for DNA pattern analysis. The proposed framework consists of two main stages. The first stage is for analyzing the DNA sequences evolution, whereas the other stage is for the reconstruction process. We utilized cellular automata (CA) rules for analyzing and predicting the DNA sequence. Then, a modified procedure for the reconstruction process is introduced, which is based on the Probabilistic Cellular Automata (PCA) integrated with Particle Swarm Optimization (PSO) algorithm. This integration makes the proposed framework more efficient and achieves optimum transition rules. Our innovated model leans on the hypothesis that mutations are probabilistic events. As a result, their evolution can be simulated as a PCA model. The main objective of this paper is to analyze various DNA sequences to predict the changes that occur in DNA during evolution (mutations). We used a similarity score as a fitness measure to detect symmetry relations, which is appropriate for numerous extremely long sequences. Results are given for the CpG-methylation-deamination processes, which are regions of DNA where a guanine nucleotide follows a cytosine nucleotide in the linear sequence of bases. The DNA evolution is handled as the evolved colored paradigms. Therefore, incorporating probabilistic components help to produce a tool capable of foretelling the likelihood of specific mutations. Besides, it shows their capabilities in dealing with complex relations.

Occurrence, genetic diversity, and recombination of maize lethal necrosis disease-causing viruses in Kenya

Maize is the most important food crop in Kenya accounting for more than 51 % of all staples grown in the country. Out of Kenya's 5.3 million ha total crops area, more than 2.1 million ha is occupied by maize which translates to 40 % of all crops area. However, with the emergence of maize lethal necrosis (MLN) disease in 2011, the average yields plummeted to all-time lows with severely affected counties recording 90-100% yield loss in 2013 and 2014. The disease is mainly caused by Maize chlorotic mottle virus (MCMV) in combination with Sugarcane mosaic virus (SCMV) or other potyviruses. In this study, a country-wide survey was carried out to assess the MLN causing viruses in Kenya, their distribution, genetic diversity, and recombination. The causative viruses of MLN were determined by RT-PCR using virus-specific primers and DAS-ELISA. Next-generation sequencing (NGS) data was generated, viral sequences identified, genetic diversity of MLN viruses was determined, and recombination was evaluated. MCMV and SCMV were detected in all the maize growing regions at varying levels of incidence, and severity while MaYMV, a polerovirus was detected in some samples through NGS. However, there were some samples in this study where only MCMV was detected with severe MLN symptoms. SCMV Sequences were highly diverse while MCMV sequences exhibited low variability. Potential recombination events were detected only in SCMV explaining the elevated level of diversity and associated risk of this virus in Kenya and the eastern Africa region.

Does the virus cross the road? Viral phylogeographic patterns among bobcat populations reflect a history of urban development

Urban development has major impacts on connectivity among wildlife populations and is thus likely an important factor shaping pathogen transmission in wildlife. However, most investigations of wildlife diseases in urban areas focus on prevalence and infection risk rather than potential effects of urbanization on transmission itself. Feline immunodeficiency virus (FIV) is a directly transmitted retrovirus that infects many felid species and can be used as a model for studying pathogen transmission at landscape scales. We investigated phylogenetic relationships among FIV isolates sampled from five bobcat (Lynx rufus) populations in coastal southern California that appear isolated due to major highways and dense urban development. Divergence dates among FIV phylogenetic lineages in several cases reflected historical urban growth and construction of major highways. We found strong FIV phylogeographic structure among three host populations north-west of Los Angeles, largely coincident with host genetic structure. In contrast, relatively little FIV phylogeographic structure existed among two genetically distinct host populations south-east of Los Angeles. Rates of FIV transfer among host populations did not vary significantly, with the lack of phylogenetic structure south-east of Los Angeles unlikely to reflect frequent contemporary transmission among populations. Our results indicate that major barriers to host gene flow can also act as barriers to pathogen spread, suggesting potentially reduced susceptibility of fragmented populations to novel directly transmitted pathogens. Infrequent exchange of FIV among host populations suggests that populations would best be managed as distinct units in the event of a severe disease outbreak. Phylogeographic inference of pathogen transmission is useful for estimating the ability of geographic barriers to constrain disease spread and can provide insights into contemporary and historical drivers of host population connectivity.

Potato Virus Y (PVY) Isolates from Solanum betaceum Represent Three Novel Recombinants Within the PVYN Strain Group and Are Unable to Systemically Spread in Potato

Tamarillo, or tree tomato (Solanum betaceum), is a perennial small tree or shrub species cultivated in subtropical areas for fresh fruit and juice production. In Ecuador, tamarillo orchards are affected by several viruses, with one previously identified as potato virus Y (PVY); however, the specific strain composition of PVY in tamarillo was not determined. In 2015 and 2016, eight tamarillo plants exhibiting symptoms of leaf drop, mosaic, and mottled fruit were sampled near Tumbaco and Quito, Ecuador. These tamarillo PVY isolates were able to systemically infect tobacco, Nicotiana benthamiana, naranjilla, and tamarillo. Seven of the eight PVY isolates from tamarillo exhibited N-serotype, while one of the PVY isolates studied, Tam15, had no identifiable serotype. One isolate, Tam17, had N-serotype but produced asymptomatic systemic infection in tobacco. In tamarillo, four tamarillo isolates induced mosaic and slight growth retardation and were unable to systemically infect pepper or potato. Tamarillo, on the other hand, was unable to support systemic infection of PVY isolates belonging to the PVY^O and PVY^Eu-N strains. The whole genomes of eight PVY isolates were sequenced from a series of overlapping RT-PCR fragments. Phylogenetically, tamarillo PVY isolates were found to belong to the large PVY^N lineage, in a new tamarillo clade. Recombination analysis revealed that these tamarillo PVY isolates represent at least three novel recombinant types not reported before. The combination of the biological and molecular properties found in these eight PVY isolates suggested the existence of a new tamarillo strain of PVY that may have coevolved with S. betaceum.

Genetic Diversity of Nine Non-Recombinant Potato virus Y Isolates From Three Biological Strain Groups: Historical and Geographical Insights

Potato virus Y (PVY) isolates from potato currently exist as a complex of six biologically defined strain groups all containing nonrecombinant isolates and at least 14 recombinant minor phylogroups. Recent studies on eight historical UK potato PVY isolates preserved since 1984 found only nonrecombinants. Here, four of five PVY isolates from cultivated potato or wild Solanum spp. collected recently in Australia, Mexico, and the U.S.A. were typed by inoculation to tobacco plants and/or serological testing using monoclonal antibodies. Next, these five modern isolates and four additional historical UK isolates belonging to biological strain groups PVY^C, PVY^Z, or PVY^N obtained from cultivated potato in 1943 to 1984 were sequenced. None of the nine complete PVY genomes obtained were recombinants. Phylogenetic analysis revealed that the four historical UK isolates were in minor phylogroups PVY^C1 (Y^C-R), PVY^O-O (Y^Z-CM1), PVY^NA-N (Y^N-M), or PVY^Eu-N (Y^N-RM), Australian isolate Y^O-BL2 was in minor phylogroup PVY^O-O5, and both Mexican isolate Y^N-Mex43 and U.S.A. isolates Y^N-MT12_Oth288, Y^N-MT12_Oth295, and Y^N-WWAA150131G42 were in minor phylogroup PVY^Eu-N. When combined, these new findings and those from the eight historical UK isolates sequenced earlier provide important historical insights concerning the diversity of early PVY populations in Europe and the appearance of recombinants in that part of the world. They and four recent Australian isolates sequenced earlier also provide geographical insights about the geographical distribution and diversity of PVY populations in Australia and North America.

Serpentoviruses: More than Respiratory Pathogens

In recent years, nidoviruses have emerged as important respiratory pathogens of reptiles, affecting captive python populations. In pythons, nidovirus (recently reclassified as serpentovirus) infection induces an inflammation of the upper respiratory and alimentary tract which can develop into a severe, often fatal proliferative pneumonia. We observed pyogranulomatous and fibrinonecrotic lesions in organ systems other than the respiratory tract during full postmortem examinations on 30 serpentovirus reverse transcription-PCR (RT-PCR)-positive pythons of varying species originating from Switzerland and Spain. The observations prompted us to study whether this not yet reported wider distribution of lesions is associated with previously unknown serpentoviruses or changes in the serpentovirus genome. RT-PCR and inoculation of Morelia viridis cell cultures served to recruit the cases and obtain virus isolates. Immunohistochemistry and immunofluorescence staining against serpentovirus nucleoprotein demonstrated that the virus infects not only a broad spectrum of epithelia (respiratory and alimentary epithelium, hepatocytes, renal tubules, pancreatic ducts, etc.), but also intravascular monocytes, intralesional macrophages, and endothelial cells. With next-generation sequencing we obtained a full-length genome for a novel serpentovirus species circulating in Switzerland. Analysis of viral genomes recovered from pythons showing serpentovirus infection-associated respiratory or systemic disease did not reveal sequence association to phenotypes; however, functional studies with different strains are needed to confirm this observation. The results indicate that serpentoviruses have a broad cell and tissue tropism, further suggesting that the course of infection could vary and involve lesions in a broad spectrum of tissues and organ systems as a consequence of monocyte-mediated viral systemic spread.IMPORTANCE During the last years, python nidoviruses (now reclassified as serpentoviruses) have become a primary cause of fatal disease in pythons. Serpentoviruses represent a threat to captive snake collections, as they spread rapidly and can be associated with high morbidity and mortality. Our study indicates that, different from previous evidence, the viruses do not only affect the respiratory tract, but can spread in the entire body with blood monocytes, have a broad spectrum of target cells, and can induce a variety of lesions. Nidovirales is an order of animal and human viruses that comprises important zoonotic pathogens such as Middle East respiratory syndrome coronavirus (MERS-CoV), severe acute respiratory syndrome coronavirus (SARS-CoV), and SARS-CoV-2. Serpentoviruses belong to the same order as the above-mentioned human viruses and show similar characteristics (rapid spread, respiratory and gastrointestinal tropism, etc.). The present study confirms the relevance of natural animal diseases to better understand the complexity of viruses of the order Nidovirales.

2000-year-old pathogen genomes reconstructed from metagenomic analysis of Egyptian mummified individuals

BACKGROUND

Recent advances in sequencing have facilitated large-scale analyses of the metagenomic composition of different samples, including the environmental microbiome of air, water, and soil, as well as the microbiome of living humans and other animals. Analyses of the microbiome of ancient human samples may provide insights into human health and disease, as well as pathogen evolution, but the field is still in its very early stages and considered highly challenging.

RESULTS

The metagenomic and pathogen content of Egyptian mummified individuals from different time periods was investigated via genetic analysis of the microbial composition of various tissues. The analysis of the dental calculus' microbiome identified Red Complex bacteria, which are correlated with periodontal diseases. From bone and soft tissue, genomes of two ancient pathogens, a 2200-year-old Mycobacterium leprae strain and a 2000-year-old human hepatitis B virus, were successfully reconstructed.

CONCLUSIONS

The results show the reliability of metagenomic studies on Egyptian mummified individuals and the potential to use them as a source for the extraction of ancient pathogen DNA.

Global allele polymorphism indicates a high rate of allele genesis at a locus under balancing selection

When selection favours rare alleles over common ones (balancing selection in the form of negative frequency-dependent selection), a locus may maintain a large number of alleles, each at similar frequency. To better understand how allelic richness is generated and maintained at such loci, we assessed 201 sequences of the complementary sex determiner (csd) of the Asian honeybee (Apis cerana), sampled from across its range. Honeybees are haplodiploid; hemizygotes at csd develop as males and heterozygotes as females, while homozygosity is lethal. Thus, csd is under strong negative frequency-dependent selection because rare alleles are less likely to end up in the lethal homozygous form. We find that in A. cerana, as in other Apis, just a few amino acid differences between csd alleles in the hypervariable region are sufficient to trigger female development. We then show that while allelic lineages are spread across geographical regions, allelic differentiation is high between populations, with most csd alleles (86.3%) detected in only one sample location. Furthermore, nucleotide diversity in the hypervariable region indicates an excess of recently arisen alleles, possibly associated with population expansion across Asia since the last glacial maximum. Only the newly invasive populations of the Austral-Pacific share most of their csd alleles. In all, the geographic patterns of csd diversity in A. cerana indicate that high mutation rates and balancing selection act together to produce high rates of allele genesis and turnover at the honeybee sex locus, which in turn leads to its exceptionally high local and global polymorphism.

Diverse variola virus (smallpox) strains were widespread in northern Europe in the Viking Age

Smallpox, one of the most devastating human diseases, killed between 300 million and 500 million people in the 20th century alone. We recovered viral sequences from 13 northern European individuals, including 11 dated to ~600-1050 CE, overlapping the Viking Age, and reconstructed near-complete variola virus genomes for four of them. The samples predate the earliest confirmed smallpox cases by ~1000 years, and the sequences reveal a now-extinct sister clade of the modern variola viruses that were in circulation before the eradication of smallpox. We date the most recent common ancestor of variola virus to ~1700 years ago. Distinct patterns of gene inactivation in the four near-complete sequences show that different evolutionary paths of genotypic host adaptation resulted in variola viruses that circulated widely among humans.

Representative Genotyping, Recombination and Evolutionary Dynamics Analysis of TSA56 Gene Segment of Orientia tsutsugamushi

Scrub typhus is a zoonotic disease caused by Orientia tsutsugamushi (O. tsutsugamushi). Orientia tsutsugamushi has various genotypes and more new strains with difference in sequences increasingly appeared. Whether the accurateness of one special nested PCR method which amplifies segment instead of entire open reading frame (ORF) sequence meets the current work of identifying new strains and classifying genotypes remains to be confirmed. And the origins and evolution of this organism have not been thoroughly elucidated. Accordingly, in this study, segments and the entire ORF of the 56-kDa type-specific antigen (TSA56) gene of O. tsutsugamushi were collected, including 209 clinically isolated strains in Guangzhou, China from 2012 to 2016 and 139 reference strains worldwide. By performing phylogenetic analysis, we proved that the accurateness of the particular PCR method which almost met detection need. This re-grouping result showed that segments perfectly represented and identified strains of Karp, Boryong, Gilliam, TA763, Kawasaki and part of Kato genotype, and this accuracy is not restricted by region and time. Sequence diversification of Shimokoshi and some Kato strains made their genotyping need to consider entire ORF sequences, but their weak recognition might not be due to recombination. The frequent genetic recombination and high point mutations contributed to genetic diversification of the TSA56 gene. Major overlapping regions of most recombination events occurred between strains of the same genotype, especially Karp and Kato genotype. And cross-genotype overlapping events occurred between Karp and Boryong/Gilliam/TA763/Kato, Kato and Kawasaki/Gilliam/TA686, Boryong and TA686, and Gilliam and Kawasaki. But Segment has quite low recombination frequency and stable mutation trend from 1943 to 2016. So segment is a relatively conserved part of the TSA56 ORF as for its stable trend of genetic diversity, and it may anchor and represent the entire TSA56 ORF gene. And genetic diversity is rejected as one potential reason for the increased incidence of scrub typhus. But an occasional recombination event created an unrecognized genotype which might be due to the breakage of VD II and AD II. Additionally, strains in Guangzhou were homologous and Karp genotype was detected as a dominant.

Structural diversity, natural selection and intragenic recombination in the Plasmodium vivax merozoite surface protein 9 locus in Thailand

The merozoite surface protein 9 (MSP9) of malarial parasite forms co-ligand complex with the 19 kDa fragment of merozoite surface protein 1 (MSP1) prior to erythrocyte invasion. Interruption of this process could hamper subsequent asexual erythrocytic development of malaria parasites; therefore, these proteins are considered potential vaccine candidates. In Plasmodium vivax, MSP9 (PvMSP9) contains both conserved and polymorphic repetitive domains that were immunogenic upon natural malaria exposure and conferred protection in vaccination studies in animal models. To investigate the extent of sequence diversity at this locus, 104 P. vivax isolates from 4 major malaria endemic areas of Thailand were analyzed. Results revealed that pvmsp9 contained 3 repeat domains (R1-R3) flanked by conserved domains. Repeat domains exhibit extensive sequence and length variation, in which 14, 39 and 16 haplotypes for domains R1-R3, respectively, circulated in this country. Sequence diversity in pvmsp9 among P. vivax isolates from each endemic area displayed population structure. The extent of sequence diversity in pvmsp9 isolates from the provinces of Tak, Chanthaburi, Ubon Ratchathani and Prachuap Khiri Khan in northwestern, eastern, northeastern and southwestern areas, respectively, was almost comparable and was remarkably higher than that from Yala/Narathiwat population in southern Thailand. Evidence for intragenic recombination in this locus was observed within each P. vivax population except among isolates from Yala and Narathiwat. Synonymous nucleotide diversity significantly exceeded nonsynonymous nucleotide diversity in domains R2 and R3, indicating purifying selection. However, micro-scale signatures of positive and negative selections occurred in both conserved and repeat domains, implying two opposing forces, probably from functional or structural constraint and host immune pressure, could have influenced diversity at this locus. The immunodominant T and B cell epitopes so far identified were invariant or highly conserved across isolates. Further analysis of global isolates is warranted for vaccine design based on this protein.

The Ingenuity of Bacterial Genomes

The genomes of bacteria contain fewer genes and substantially less noncoding DNA than those of eukaryotes, and as a result, they have much less raw material to invent new traits. Yet, bacteria are vastly more taxonomically diverse, numerically abundant, and globally successful in colonizing new habitats compared to eukaryotes. Although bacterial genomes are generally considered to be optimized for efficient growth and rapid adaptation, nonadaptive processes have played a major role in shaping the size, contents, and compact organization of bacterial genomes and have allowed the establishment of deleterious traits that serve as the raw materials for genetic innovation.

Phylogenetic and Timescale Analysis of Barmah Forest Virus as Inferred from Genome Sequence Analysis

Barmah Forest virus (BFV) is a medically important mosquito-borne alphavirus endemic to Australia. Symptomatic disease can be a major cause of morbidity, associated with fever, rash, and debilitating arthralgia. BFV disease is similar to that caused by Ross River virus (RRV), the other major Australian alphavirus. Currently, just four BFV whole-genome sequences are available with no genome-scale phylogeny in existence to robustly characterise genetic diversity. Thirty novel genome sequences were derived for this study, for a final 34-taxon dataset sampled over a 44 year period. Three distinct BFV genotypes were characterised (G1-3) that have circulated in Australia and Papua New Guinea (PNG). Evidence of spatio-temporal co-circulation of G2 and G3 within regions of Australia was noted, including in the South West region of Western Australia (WA) during the first reported disease outbreaks in the state's history. Compared with RRV, the BFV population appeared more stable with less frequent emergence of novel lineages. Preliminary in vitro assessment of RRV and BFV replication kinetics found that RRV replicates at a significantly faster rate and to a higher, more persistent titre compared with BFV, perhaps indicating mosquitoes may be infectious with RRV for longer than with BFV. This investigation resolved a greater diversity of BFV, and a greater understanding of the evolutionary dynamics and history was attained.

Analysis of codon usage bias in potato virus Y non-recombinant strains

Potato virus Y (PVY) is a member of the genus Potyvirus, family Potyviridae, is considered one of the most devastating pest affecting economically important crops, such as potato, tobacco, tomato and pepper, representing a serious threat due to high incidence and worldwide distribution. Its economic significance as well as it biological and molecular complexities have aroused great attention, thus several studies have explore it genetic characteristics. However, little is known about PVY codon usage. To shed light on the relation of codon usage among viruses and their hosts is extremely important to understand virus survival, fitness and evolution. In this study, we performed a comprehensive analysis of codon usage and composition of PVY non-recombinant strains (PVY^N-NA, PVY^Eu-N, PVY^O, PVY^O5, PVY^C) based on 130 complete open reading frame sequences extracted from public databases. Furthermore, similarities between the synonymous codon usage of PVY and its main hosts were investigated. The results obtained in the current study suggest that the overall codon usage among PVY genotypes is similar and slightly biased. PVY codon usage is strongly influenced by mutational bias, but also by G + C compositional constraint and dinucleotide composition. Furthermore, similarities among codon usage preferences between PVY strains and analyzed hosts were observed.

Complete genome characterization of porcine circovirus 3 recovered from wild boars in Southern Brazil

In the present study, the complete nucleotide sequence of porcine circovirus 3 (PCV3) recovered from wild boars lymph nodes is described. The full genome was named PCV3-wb/Br/RS and comprises 2,000 nucleotides with two open reading frames (ORFs) with a stem-loop motif in intergenic region. The ORFs are oriented in opposite directions and encode the putative capsid (Cap) and replicase (Rep) proteins. Based on amino acid motif analysis, PCV3-wb/Br/RS as well as most of the sequences from wild boars are classified as PCV3b. Phylogenetic analysis including 97 PCV3 sequences available in databases showed that the PCV3-wb/Br/RS genome is more closely related to genomes recovered in Spain, China, Germany and Denmark. Phylogenetic inferences among PCV3-wb/Br/RS and other circoviruses confirmed that these seem to have a most recent common ancestor with bat-associated circoviruses. In addition, PCV3 infection was investigated by real-time PCR in a cohort of 80 wild boars in Southern Brazil. A total of 29 animals (36.3%) were PCV3-positive leading the conclusion that PCV3 is circulating in the wild boar population in Southern Brazil. The role played by PCV3-like infections in wild boars and the risk these could pose to commercial swine production within that region remains to be further investigated.

Viral Ecogenomics of Arctic Cryopeg Brine and Sea Ice

Arctic regions, which are changing rapidly as they warm 2 to 3 times faster than the global average, still retain microbial habitats that serve as natural laboratories for understanding mechanisms of microbial adaptation to extreme conditions. Seawater-derived brines within both sea ice (sea-ice brine) and ancient layers of permafrost (cryopeg brine) support diverse microbes adapted to subzero temperatures and high salinities, yet little is known about viruses in these extreme environments, which, if analogous to other systems, could play important evolutionary and ecosystem roles. Here, we characterized viral communities and their functions in samples of cryopeg brine, sea-ice brine, and melted sea ice. Viral abundance was high in cryopeg brine (1.2 × 108 ml-1) and much lower in sea-ice brine (1.3 × 105 to 2.1 × 105 ml-1), which roughly paralleled the differences in cell concentrations in these samples. Five low-input, quantitative viral metagenomes were sequenced to yield 476 viral populations (i.e., species level; ≥10 kb), only 12% of which could be assigned taxonomy by traditional database approaches, indicating a high degree of novelty. Additional analyses revealed that these viruses: (i) formed communities that differed between sample type and vertically with sea-ice depth; (ii) infected hosts that dominated these extreme ecosystems, including Marinobacter, Glaciecola, and Colwellia; and (iii) encoded fatty acid desaturase (FAD) genes that likely helped their hosts overcome cold and salt stress during infection, as well as mediated horizontal gene transfer of FAD genes between microbes. Together, these findings contribute to understanding viral abundances and communities and how viruses impact their microbial hosts in subzero brines and sea ice.IMPORTANCE This study explores viral community structure and function in remote and extreme Arctic environments, including subzero brines within marine layers of permafrost and sea ice, using a modern viral ecogenomics toolkit for the first time. In addition to providing foundational data sets for these climate-threatened habitats, we found evidence that the viruses had habitat specificity, infected dominant microbial hosts, encoded host-derived metabolic genes, and mediated horizontal gene transfer among hosts. These results advance our understanding of the virosphere and how viruses influence extreme ecosystems. More broadly, the evidence that virally mediated gene transfers may be limited by host range in these extreme habitats contributes to a mechanistic understanding of genetic exchange among microbes under stressful conditions in other systems.

Microevolution in the major outer membrane protein OmpA of Acinetobacter baumannii

Acinetobacter baumannii is nowadays a relevant nosocomial pathogen characterized by multidrug resistance (MDR) and concomitant difficulties to treat infections. OmpA is the most abundant A. baumannii outer membrane (OM) protein, and is involved in virulence, host-cell recognition, biofilm formation, regulation of OM stability, permeability and antibiotic resistance. OmpA members are two-domain proteins with an N-terminal eight-stranded β-barrel domain with four external loops (ELs) interacting with the environment, and a C-terminal periplasmic domain binding non-covalently to the peptidoglycan. Here, we combined data from genome sequencing, phylogenetic and multilocus sequence analyses from 975 strains/isolates of the Acinetobacter calcoaceticus/Acinetobacter baumannii complex (ACB), 946 from A. baumannii, to explore ompA microevolutionary divergence. Five major ompA variant groups were identified (V1 to V5) in A. baumannii, encompassing 52 different alleles coding for 23 different proteins. Polymorphisms were concentrated in five regions corresponding to the four ELs and the C-terminal end, and provided evidence for intra-genic recombination. ompA variants were not randomly distributed across the A. baumannii phylogeny, with the most frequent V1(lct)a1 allele found in most clonal complex 2 (CC2) strains and the second most frequent V2(lct)a1 allele in the majority of CC1 strains. Evidence was found for assortative exchanges of ompA alleles not only between separate A. baumannii lineages, but also different ACB species. The overall results have implications for A. baumannii evolution, epidemiology, virulence and vaccine design.

Phylogenetic and recombinant analyses of complete coding sequences of DENV-1 from field-caught mosquitoes in Thailand

Dengue diseases are a group of infectious diseases that have been widespread throughout the world for several decades. Dengue outbreaks have occurred in many parts of the world, including Southeast Asia. The outbreak and the severity of the diseases depend on many factors including dengue genotype. Data on the genetic variation of dengue virus is highly informative for dengue protection plans and vaccine development. In this study, we focused on the analyses of genetic variation and amino acid changes of the whole coding sequences in two dengue strains isolated from Aedes mosquitoes in Bangkok, Thailand, an endemic area. The strains were identified as belonging to dengue virus serotype 1 (DENV-1) genotype I, and have unique nucleotide sequences. In the recombinant analysis, these strains were identified as recombinants derived from Chinese counterparts for both the major (DENV-1 genotype I) and the minor (DENV-1 genotype IV) parental strains. The recombination event occurred within the prM and E genes. This corresponded with the result of multiple alignments where several amino acid residues between the recombination breakpoints in the strains were identical to those in DENV-1 genotype IV. Several of the amino acid substitutions also have changed the amino acid properties, which might affect viral infection and antigenicity. These results provide insight into the genetic variation of DENV in this endemic area, which might have been involved in the dengue outbreak and high numbers of dengue haemorrhagic fever and dengue shock syndrome cases in Thailand in 2015.

Identification of Reptarenaviruses, Hartmaniviruses, and a Novel Chuvirus in Captive Native Brazilian Boa Constrictors with Boid Inclusion Body Disease

Boid inclusion body disease (BIBD) is a transmissible viral disease of captive snakes that causes severe losses in snake collections worldwide. It is caused by reptarenavirus infection, which can persist over several years without overt signs but is generally associated with the eventual death of the affected snakes. Thus far, reports have confirmed the existence of reptarenaviruses in captive snakes in North America, Europe, Asia, and Australia, but there is no evidence that it also occurs in wild snakes. BIBD affects boa species within the subfamily Boinae and pythons in the family Pythonidae, the habitats of which do not naturally overlap. Here, we studied Brazilian captive snakes with BIBD using a metatranscriptomic approach, and we report the identification of novel reptarenaviruses, hartmaniviruses, and a new species in the family Chuviridae The reptarenavirus L segments identified are divergent enough to represent six novel species, while we found only a single novel reptarenavirus S segment. Until now, hartmaniviruses had been identified only in European captive boas with BIBD, and the present results increase the number of known hartmaniviruses from four to six. The newly identified chuvirus showed 38.4%, 40.9%, and 48.1% amino acid identity to the nucleoprotein, glycoprotein, and RNA-dependent RNA polymerase, respectively, of its closest relative, Guangdong red-banded snake chuvirus-like virus. Although we cannot rule out the possibility that the found viruses originated from imported snakes, the results suggest that the viruses could circulate in indigenous snake populations.IMPORTANCE Boid inclusion body disease (BIBD), caused by reptarenavirus infection, affects captive snake populations worldwide, but the reservoir hosts of reptarenaviruses remain unknown. Here, we report the identification of novel reptarenaviruses, hartmaniviruses, and a chuvirus in captive Brazilian boas with BIBD. Three of the four snakes studied showed coinfection with all three viruses, and one of the snakes harbored three novel reptarenavirus L segments and one novel S segment. The samples originated from collections with Brazilian indigenous snakes only, which could indicate that these viruses circulate in wild snakes. The findings could further indicate that boid snakes are the natural reservoir of reptarena- and hartmaniviruses commonly found in captive snakes. The snakes infected with the novel chuvirus all suffered from BIBD; it is therefore not possible to comment on its potential pathogenicity and contribution to the observed changes in the present case material.

Recombination Between High-Risk Human Papillomaviruses and Non-Human Primate Papillomaviruses: Evidence of Ancient Host Switching Among Alphapapillomaviruses

We use all the currently known 405 Papillomavirus (PV) sequences, 343 curated PV sequences from both humans and animals from the PAVE data base, to analyse the recombination dynamics of these viruses at the whole genome levels. After showing some evidence of human and non-human primate PV recombination, we report a comprehensive recombination analysis of all currently known 82 Alphapapillomaviruses (Alpha-PVs). We carried out an exploratory study and found novel recombination events between High-Risk HPV Types and Macaca fascicularis PV1 (MfPV1), Macaca Fuscata PV2 (MfuPV2) and Pan Paniscus PV1 (PpPV1) Papillomaviruses. This is the first evidence of interactions between PVs from different hosts and hence postulates the likelihood of ancient host switching among Alpha-PVs. Notwithstanding these results should be interpreted with caution because the major and minor parents indicated by RDP4 program are simply the sequences in the alignment that most closely resemble the actual parents. We found statistically significant differences between the phylogenies of the PV sequences with recombination regions and PV sequences without recombination regions using the Shimodaira–Hasegawa phylogenetic incongruence testing. We show that not more than 76MYA Alpha-PVs were in the same biological niche, a pre-requisite for recombination, and as the hosts evolved and diversified, the viruses adapted to specific host niches which eventually led to coevolution with specific hosts before speciation of primate species. Thus providing evidence that in ancient times no earlier than the Cretaceous period of the Mesozoic age, Alpha-PVs recombined and switched hosts, but whether this host switching is occurring currently is unknown. However, a clearer picture of the PVs evolutionary landscape can only be achieved with the incremental discovery of PV sequences, especially from the animal kingdom.

Coding potential and sequence conservation of SARS-CoV-2 and related animal viruses

In December 2019, a novel human-infecting coronavirus (SARS-CoV-2) was recognized in China. In a few months, SARS-CoV-2 has caused thousands of disease cases and deaths in several countries. Phylogenetic analyses indicated that SARS-CoV-2 clusters with SARS-CoV in the Sarbecovirus subgenus and viruses related to SARS-CoV-2 were identified from bats and pangolins. Coronaviruses have long and complex genomes with high plasticity in terms of gene content. To date, the coding potential of SARS-CoV-2 remains partially unknown. We thus used available sequences of bat and pangolin viruses to determine the selective events that shaped the genome structure of SARS-CoV-2 and to assess its coding potential. By searching for signals of significantly reduced variability at synonymous sites (dS), we identified six genomic regions, one of these corresponding to the programmed −1 ribosomal frameshift. The most prominent signal of dS reduction was observed within the E gene. A genome-wide analysis of conserved RNA structures indicated that this region harbors a putative functional RNA element that is shared with the SARS-CoV lineage. Additional signals of reduced dS indicated the presence of internal ORFs. Whereas the presence ORF9a (internal to N) was previously proposed by homology with a well characterized protein of SARS-CoV, ORF3h (for hypothetical, within ORF3a) was not previously described. The predicted product of ORF3h has 90% identity with the corresponding predicted product of SARS-CoV and displays features suggestive of a viroporin. Finally, analysis of the putative ORF10 revealed high dN/dS (3.82) in SARS-CoV-2 and related coronaviruses. In the SARS-CoV lineage, the ORF is predicted to encode a truncated protein and is neutrally evolving. These data suggest that ORF10 encodes a functional protein in SARS-CoV-2 and that positive selection is driving its evolution. Experimental analyses will be necessary to validate and characterize the coding and non-coding functional elements we identified.

•

We analyzed the coding region of SARS-CoV-2 and related bat/pangolin viruses.

•

We identified six regions of significantly low variability at sysnonymous sites.

•

One of these corresponds to a conserved RNA structure shared with the SARS-CoV lineage.

•

The dS reduction within ORF3a corresponds to a potential ORF encoding a viroporin.

•

In SARS-CoV-2 and related viruses, the putative 3′ terminal ORF10 has high dN/dS.

Genetic changes and host adaptability in sugarcane mosaic virus based on complete genome sequences

Sugarcane mosaic virus (SCMV), a member of the genus Potyvirus in the family Potyviridae, is an important pathogen that causes mosaic diseases in maize, sugarcane, canna and other graminaceous species worldwide. Previously, several reports have showed the genetic variation and population structure of SCMV. However, the evolutionary dynamics, synonymous codon usage pattern and adaptive evolution of the virus is unclear. In this study, we performed comprehensive analyses of phylodynamics, composition bias and codon usage of SCMV using 108 complete genomic sequences. Our phylogenetic analysis found six host- and geographically confined phylogenetic lineages within the SCMV non-recombinant isolates. We found a relatively stable and conserved genomic composition with a lower codon usage choice in the SCMV protein coding sequences. Mutation pressure and natural selection have shaped the codon usage patterns of the SCMV protein coding sequences with natural selection being the dominant factor. The codon adaptation index (CAI), relative codon deoptimization index (RCDI) and similarity index (SiD) analyses revealed a stronger correlation between SCMV and maize than between SCMV and sugarcane or canna. Our study is the first to evaluate the codon usage pattern of SCMV based on complete sequences and may provide a better understanding of the origin of SCMV and its evolutionary patterns for future research.

The Genetic Diversity of HIV-1 Quasispecies Within Primary Infected Individuals

HIV has remarkable genetic diversity among populations. The diversity has critical impacts on transmission, immune escape, pathogenesis, and clinical management. HIV-1 diversity originates from frequent mutation and recombination during reverse transcription. This work focuses on the quasispecies genetic dynamics within individuals with primary infections. Eleven men who have sex with men from the Beijing PRIMO Clinical Cohort were identified as primary infection and had three or four series of their anticoagulant blood samples collected. Viral RNA was extracted and amplified using single-genome amplification. Products of the gp120 gene that met single-genome amplification requirements were sequenced. Subtype assortment of all collected sequences was performed using both the jumping profile hidden Markov model (jpHMM) and REGA. Quasispecies diversity at each time was estimated using Mega 6. Intrapatient recombination was analyzed using RDP4. According to the Fiebig classification system, YA-81 belongs to stage III and YA-113 belongs to stage IV. The other samples are all associated with the infection stage of V/VI. YA113 had a dual infection with subtype B and a new unique recombinant form involving CRF01_AE and C. The other eight were infected with CRF01_AE, one was infected with B/C recombinant, and the last one with B. Of the 10 single infections, 8 were caused by 1 founder virus. They all displayed a sharp increase of quasispecies diversity during the sampling times. Two were caused by at least two founder viruses. The diversity of these strains starts at a significantly high level and is followed by a relatively steady trend. Critically, the separate subtypes YA113-B and YA113-CRF01_AE/C both showed a similar trend to those infected by a single founder virus. Recombination analysis revealed that 5 of 11 cases underwent detectable intrapatient recombination. These findings indicate that tracing the dynamics of HIV-1 quasispecies during early infection may be relevant and valuable for understanding pathways of viral diversification and immune escape.

FLAVI: An Amino Acid Substitution Model for Flaviviruses

Amino acid substitution models represent substitution rates among amino acids during the evolution. The models play an important role in analyzing protein sequences, especially inferring phylogenies. The rapid evolution of flaviviruses is expanding the threat in public health. A number of models have been estimated for some viruses, however, they are unable to properly represent amino acid substitution patterns of flaviviruses. In this study, we collected protein sequences from the flavivirus genus to specifically estimate an amino acid substitution model, called FLAVI, for flaviviruses. Experiments showed that the collected dataset was sufficient to estimate a stable model. More importantly, the FLAVI model was remarkably better than other existing models in analyzing flavivirus protein sequences. We recommend researchers to use the FLAVI model when studying protein sequences of flaviviruses or closely related viruses.

Evolutionary analysis of exogenous and integrated HHV-6A/HHV-6B populations

Human betaherpesviruses 6A and 6B (HHV-6A and HHV-6B) are highly prevalent in human populations. The genomes of these viruses can be stably integrated at the telomeres of human chromosomes and be vertically transmitted (inherited chromosomally integrated HHV-6A/HHV-6B, iciHHV-6A/iciHHV-6B). We reconstructed the population structures of HHV-6A and HHV-6B, showing that HHV-6A diverged less than HHV-6B genomes from the projected common ancestral population. Thus, HHV-6B genomes experienced stronger drift, as also supported by calculation of nucleotide diversity and Tajima's D. Analysis of ancestry proportions indicated that HHV-6A exogenous viruses and iciHHV-6A derived most of their genomes from distinct ancestral sources. Conversely, ancestry proportions were similar in exogenous HHV-6B viruses and iciHHV-6B. In line with previous indications, this suggests the distinct exogenous viral populations that originated iciHHV-6B in subjects with European and Asian ancestry are still causing infections in the corresponding geographic areas. Notably, for both iciHHV-6A and iciHHV-6B, we found that European and American sequences tend to have high proportions of ancestry from viral populations that experienced considerable drift, suggesting that they underwent one or more bottlenecks followed by population expansion. Finally, analysis of HHV-6B exogenous viruses sampled in Japan indicated that proportions of ancestry components of most of these viruses are different from the majority of those sampled in the USA. More generally, we show that, in both viral species, both integrated and exogenous viral genomes have different ancestry components, partially depending on geographic location. It would be extremely important to determine whether such differences account for the diversity of HHV-6A/HHV-6B-associated clinical symptoms and epidemiology. Also, the sequencing of additional exogenous and integrated viral genomes will be instrumental to confirm and expand our conclusions, which are based on a relatively small number of genomes, sequenced with variable quality, and with unequal sampling in terms of geographic origin.

Population structure of Lassa Mammarenavirus in West Africa

Lassa mammarenavirus (LASV) is the etiologic agent of Lassa fever. In endemic regions in West Africa, LASV genetic diversity tends to cluster by geographic area. Seven LASV lineages are recognized, but the role of viral genetic determinants on disease presentation in humans is uncertain. We investigated the geographic structure and distribution of LASV in West Africa. We found strong spatial clustering of LASV populations, with two major east–west and north–south diversity gradients. Analysis of ancestry components indicated that known LASV lineages diverged from an ancestral population that most likely circulated in Nigeria, although alternative locations, such as Togo, cannot be excluded. Extant sequences carrying the largest contribution of this ancestral population include the prototype Pinneo strain, the Togo isolates, and a few viruses isolated in Nigeria. The LASV populations that experienced the strongest drift circulate in Mali and the Ivory Coast. By focusing on sequences form a single LASV sublineage (IIg), we identified an ancestry component possibly associated with protection from a fatal disease outcome. Although the same ancestry component tends to associate with lower viral loads in plasma, the small sample size requires that these results are treated with extreme caution.

In silico Analysis of Genetic Diversity of Human Hepatitis B Virus in Southeast Asia, Australia and New Zealand

The extent of whole genome diversity amongst hepatitis B virus (HBV) genotypes is not well described. This study aimed to update the current distribution of HBV types and to investigate mutation rates and nucleotide diversity between genotypes in Southeast Asia, Australia and New Zealand. We retrieved 930 human HBV complete genomes from these regions from the NCBI nucleotide database for genotyping, detection of potential recombination, serotype prediction, mutation identification and comparative genome analyses. Overall, HBV genotypes B (44.1%) and C (46.2%) together with predicted serotypes adr (36%), adw2 (29%) and ayw1 (19.9%) were the most commonly circulating HBV types in the studied region. The three HBV variants identified most frequently were p.V5L, c.1896G>A and double mutation c.1762A>T/c.1764G>A, while genotypes B and C had the widest range of mutation types. The study also highlighted the distinct nucleotide diversity of HBV genotypes for whole genome and along the genome length. Therefore, this study provided a robust update to HBV currently circulating in Southeast Asia, Australia and New Zealand as well as an insight into the association of HBV genetic hypervariability and prevalence of well reported mutations.

A Novel Divergent Geminivirus Identified in Asymptomatic New World Cactaceae Plants

Cactaceae comprise a diverse and iconic group of flowering plants which are almost exclusively indigenous to the New World. The wide variety of growth forms found amongst the cacti have led to the trafficking of many species throughout the world as ornamentals. Despite the evolution and physiological properties of these plants having been extensively studied, little research has focused on cactus-associated viral communities. While only single-stranded RNA viruses had ever been reported in cacti, here we report the discovery of cactus-infecting single-stranded DNA viruses. These viruses all apparently belong to a single divergent species of the family Geminiviridae and have been tentatively named Opuntia virus 1 (OpV1). A total of 79 apparently complete OpV1 genomes were recovered from 31 different cactus plants (belonging to 20 different cactus species from both the Cactoideae and Opuntioideae clades) and from nine cactus-feeding cochineal insects (Dactylopius sp.) sampled in the USA and Mexico. These 79 OpV1 genomes all share > 78.4% nucleotide identity with one another and < 64.9% identity with previously characterized geminiviruses. Collectively, the OpV1 genomes display evidence of frequent recombination, with some genomes displaying up to five recombinant regions. In one case, recombinant regions span ~40% of the genome. We demonstrate that an infectious clone of an OpV1 genome can replicate in Nicotiana benthamiana and Opuntia microdasys. In addition to expanding the inventory of viruses that are known to infect cacti, the OpV1 group is so distantly related to other known geminiviruses that it likely represents a new geminivirus genus. It remains to be determined whether, like its cactus hosts, its geographical distribution spans the globe.

Evolution of major histocompatibility complex class I genes in the sable Martes zibellina (Carnivora, Mustelidae)

The molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immune response among vertebrates. We investigated the molecular evolution of MHC class I genes in the sable Martes zibellina. We isolated 26 MHC class I sequences, including 12 putatively functional sequences and 14 pseudogene sequences, from 24 individuals from two geographic areas of northeast China. The number of putatively functional sequences found in a single individual ranged from one to five, which might be at least 1-3 loci. We found that both balancing selection and recombination contribute to evolution of MHC class I genes in M. zibellina. In addition, we identified a candidate nonclassical MHC class I lineage in Carnivora, which may have preceded the divergence (about 52-57 Mya) of Caniformia and Feliformia. This may contribute to further understanding of the origin and evolution of nonclassical MHC class I genes. Our study provides important immune information of MHC for M. zibellina, as well as other carnivores.

Intra-strain biological and epidemiological characterization of plum pox virus

Plum pox virus (PPV) is one of the most important plant viruses causing serious economic losses. Thus far, strain typing based on the definition of 10 monophyletic strains with partially differentiable biological properties has been the sole approach used for epidemiological characterization of PPV. However, elucidating the genetic determinants underlying intra-strain biological variation among populations or isolates remains a relevant but unexamined aspect of the epidemiology of the virus. In this study, based on complete nucleotide sequence information of 210 Japanese and 47 non-Japanese isolates of the PPV-Dideron (D) strain, we identified five positively selected sites in the PPV-D genome. Among them, molecular studies showed that amino acid substitutions at position 2,635 in viral replicase correlate with viral titre and competitiveness at the systemic level, suggesting that amino acid position 2,635 is involved in aphid transmission efficiency and symptom severity. Estimation of ancestral genome sequences indicated that substitutions at amino acid position 2,635 were reversible and peculiar to one of two genetically distinct PPV-D populations in Japan. The reversible amino acid evolution probably contributes to the dissemination of the virus population. This study provides the first genomic insight into the evolutionary epidemiology of PPV based on intra-strain biological variation ascribed to positive selection.

Interspecies Transmission, Genetic Diversity, and Evolutionary Dynamics of Pseudorabies Virus

BACKGROUND

Pseudorabies virus (PRV) causes Aujeszky's disease in pigs and can be transmitted to other mammals, including humans. In the current study, we systematically studied the interspecies transmission and evolutionary history of PRV.

METHODS

We performed comprehensive analysis on the phylodynamics, selection, and structural biology to summarize the phylogenetic and adaptive evolution of PRV based on all available full-length and major glycoprotein sequences.

RESULTS

PRV can be divided into 2 main clades with frequent interclade and intraclade recombination. Clade 2.2 (variant PRV) is currently the most prevalent genotype worldwide, and most commonly involved in cross-species transmission events (including humans). We also found that the population size of clade 2.2 has increased since 2011, and the effective reproduction number was >1 from 2011 to 2016, indicating that PRV may be still circulating in swine herds and is still a risk in relation with cross-species transmission in China. Of note, we identified amino acid sites in some important glycoproteins gB, gC, gD, and gE that may be associated with PRV adaptation to new hosts and immune escape to vaccines.

CONCLUSIONS

Our study provides important genetic insight into the interspecies transmission and evolution of PRV within and between different hosts that warrant additional surveillance.

Genomic epidemiological characteristics of dengue fever in Guangdong province, China from 2013 to 2017

Dengue fever, a mosquito-borne viral disease in humans, has been endemic in many Southeast Asian countries. Since its first outbreak in 1978 in Foshan, Guangdong province, China, dengue has been continually epidemic in recent years in Guangdong, which raised the concern whether dengue infection is endemic in Guangdong. In this study, we performed phylogenetic, recombinant, and nucleotide variation analyses of 114 complete genome sequences of dengue virus serotypes 1–4 (DENV1-4) collected from 2013 to 2017 in 18 of 21 cities of Guangdong. Phylogenetic analyses revealed that DENV sequences did not form a single cluster, indicating that dengue fever was not endemic in Guangdong, although DENV1-4 co-circulated in Guangdong. Twenty intra-serotype recombinant isolates involving DENV1-4 were detected, but no inter-serotype recombinant events were identified in this study. Additionally, the most recombinant events were detected simultaneously in the gene NS3 of DENV1-4. Nucleotide variation analyses showed that no significant intra-serotype differences were observed, whereas more significant inter-subtype differences were discovered in non-structural genes than in structural genes. Our investigation will facilitate the understanding of the current prevalent status of dengue fever in Guangdong and contribute to designing more effective preventive and control strategies for dengue infection.

In 1978, dengue fever was first reported in Guangdong province, China, and this has been continuously prevalent in Guangdong in recent years. This is responsible for the heavy burden on the control of dengue, and raises the concern about whether dengue outbreaks have become endemic in Guangdong. Previous studies based on single E gene or few full-length genome sequences were inconclusive. In this study, we sequenced 114 DENV complete genomes of DENV1-4 obtained from 2013 to 2017 in Guangdong and further analyzed the epidemiological and molecular characteristics. Phylogenetic analyses revealed that dengue fever was not endemic in Guangdong, which was indirectly supported by results of our recombination analyses. Nucleotide variation analyses indicated that purification selection shaped dengue virus population. Our investigation will facilitate the development of more effective epidemiological surveillance strategies for dengue infection.