Ecological and Evolutionary Processes Shaping Viral Genetic Diversity

Evolutionary Analysis of a Parrot Bornavirus 2 Detected in a Sulphur-Crested Cockatoo (Cacatua galerita) Suggests a South American Ancestor

Parrot bornavirus (PaBV) is an RNA virus that causes Proventricular Dilatation Disease (PDD), neurological disorders, and death in Psittaciformes. Its diversity in South America is poorly known. We examined a Cacatua galerita presenting neuropathies, PDD, and oculopathies as the main signs. We detected PaBV through reverse transcription polymerase chain reaction (RT-PCR) and partial sequencing of the nucleoprotein (N) and matrix (M) genes. Maximum likelihood and Bayesian phylogenetic inferences classified it as PaBV-2. The nucleotide identity of the sequenced strain ranged from 88.3% to 90.3% against genotype PaBV-2 and from 80.2% to 84.4% against other genotypes. Selective pressure analysis detected signs of episodic diversifying selection in both the N and M genes. No recombination events were detected. Phylodynamic analysis estimated the time to the most recent common ancestor (TMRCA) as the year 1758 for genotype PaBV-2 and the year 1049 for the Orthobornavirus alphapsittaciforme species. Substitution rates were estimated at 2.73 × 10-4 and 4.08 × 10-4 substitutions per year per site for N and M, respectively. The analysis of population dynamics showed a progressive decline in the effective population size during the last century. Timescale phylogeographic analysis revealed a potential South American ancestor as the origin of genotypes 1, 2, and 8. These results contribute to our knowledge of the evolutionary origin, diversity, and dynamics of PaBVs in South America and the world. Additionally, it highlights the importance of further studies in captive Psittaciformes and the potential impact on endangered wild birds.

Signatures of adaptive decreased virulence of deformed wing virus in an isolated population of wild honeybees (Apis mellifera)

Understanding the ecological and evolutionary processes that drive host-pathogen interactions is critical for combating epidemics and conserving species. The Varroa destructor mite and deformed wing virus (DWV) are two synergistic threats to Western honeybee (Apis mellifera) populations across the globe. Distinct honeybee populations have been found to self-sustain despite Varroa infestations, including colonies within the Arnot Forest outside Ithaca, NY, USA. We hypothesized that in these bee populations, DWV has been selected to produce an avirulent infection phenotype, allowing for the persistence of both host and disease-causing agents. To investigate this, we assessed the titre of viruses in bees from the Arnot Forest and managed apiaries, and assessed genomic variation and virulence differences between DWV isolates. Across groups, we found viral abundance was similar, but DWV genotypes were distinct. We also found that infections with isolates from the Arnot Forest resulted in higher survival and lower rates of symptomatic deformed wings, compared to analogous isolates from managed colonies, providing preliminary evidence to support the hypothesis of adaptive decreased viral virulence. Overall, this multi-level investigation of virus genotype and phenotype indicates that host ecological context can be a significant driver of viral evolution and host-pathogen interactions in honeybees.

Modeling the Climatic Suitability of COVID-19 Cases in Brazil

Studies have shown that climate may affect the distribution of coronavirus disease (COVID-19) and its incidence and fatality rates. Here, we applied an ensemble niche modeling approach to project the climatic suitability of COVID-19 cases in Brazil. We estimated the cumulative incidence, mortality rate, and fatality rate of COVID-19 between 2020 and 2021. Seven statistical algorithms (MAXENT, MARS, RF, FDA, CTA, GAM, and GLM) were selected to model the climate suitability for COVID-19 cases from diverse climate data, including temperature, precipitation, and humidity. The annual temperature range and precipitation seasonality showed a relatively high contribution to the models, partially explaining the distribution of COVID-19 cases in Brazil based on the climatic suitability of the territory. We observed a high probability of climatic suitability for high incidence in the North and South regions and a high probability of mortality and fatality rates in the Midwest and Southeast regions. Despite the social, viral, and human aspects regulating COVID-19 cases and death distribution, we suggest that climate may play an important role as a co-factor in the spread of cases. In Brazil, there are regions with a high probability that climatic suitability will contribute to the high incidence and fatality rates of COVID-19 in 2020 and 2021.

Viral suppressor of RNA silencing in vascular plants also interferes with the development of the bryophyte Physcomitrella patens

Plant viruses are important pathogens able to overcome plant defense mechanisms using their viral suppressors of RNA silencing (VSR). Small RNA pathways of bryophytes and vascular plants have significant similarities, but little is known about how viruses interact with mosses. This study elucidated the responses of Physcomitrella patens to two different VSRs. We transformed P. patens plants to express VSR P19 from tomato bushy stunt virus and VSR 2b from cucumber mosaic virus, respectively. RNA sequencing and quantitative PCR were used to detect the effects of VSRs on gene expression. Small RNA (sRNA) sequencing was used to estimate the influences of VSRs on the sRNA pool of P. patens. Expression of either VSR-encoding gene caused developmental disorders in P. patens. The transcripts of four different transcription factors (AP2/erf, EREB-11 and two MYBs) accumulated in the P19 lines. sRNA sequencing revealed that VSR P19 significantly changed the microRNA pool in P. patens. Our results suggest that VSR P19 is functional in P. patens and affects the abundance of specific microRNAs interfering with gene expression. The results open new opportunities for using Physcomitrella as an alternative system to study plant-virus interactions.

Simulated vector transmission differentially influences dynamics of two viral variants of deformed wing virus in honey bees (Apis mellifera)

Understanding how vectors alter the interactions between viruses and their hosts is a fundamental question in virology and disease ecology. In honey bees, transmission of deformed wing virus (DWV) by parasitic Varroa mites has been associated with elevated disease and host mortality, and Varroa transmission has been hypothesized to lead to increased viral titres or select for more virulent variants. Here, we mimicked Varroa transmission by serially passaging a mixed population of two DWV variants, A and B, by injection through in vitro reared honey bee pupae and tracking these viral populations through five passages. The DWV-A and DWV-B variant proportions shifted dynamically through passaging, with DWV-B outcompeting DWV-A after one passage, but levels of both variants becoming equivalent by Passage 5. Sequencing analysis revealed a dominant, recombinant DWV-B strain (DWV-A derived 5′ IRES region with the rest of the genome DWV-B), with low nucleotide diversity that decreased through passaging. DWV-A populations had higher nucleotide diversity compared to DWV-B, but this also decreased through passaging. Selection signatures were found across functional regions of the DWV-A and DWV-B genomes, including amino acid mutations in the putative capsid protein region. Simulated vector transmission differentially impacted two closely related viral variants which could influence viral interactions with the host, demonstrating surprising plasticity in vector-host-viral dynamics.

Viruses infecting a warm water picoeukaryote shed light on spatial co-occurrence dynamics of marine viruses and their hosts

The marine picoeukaryote Bathycoccus prasinos has been considered a cosmopolitan alga, although recent studies indicate two ecotypes exist, Clade BI (B. prasinos) and Clade BII. Viruses that infect Bathycoccus Clade BI are known (BpVs), but not that infect BII. We isolated three dsDNA prasinoviruses from the Sargasso Sea against Clade BII isolate RCC716. The BII-Vs do not infect BI, and two (BII-V2 and BII-V3) have larger genomes (~210 kb) than BI-Viruses and BII-V1. BII-Vs share ~90% of their proteins, and between 65% to 83% of their proteins with sequenced BpVs. Phylogenomic reconstructions and PolB analyses establish close-relatedness of BII-V2 and BII-V3, yet BII-V2 has 10-fold higher infectivity and induces greater mortality on host isolate RCC716. BII-V1 is more distant, has a shorter latent period, and infects both available BII isolates, RCC716 and RCC715, while BII-V2 and BII-V3 do not exhibit productive infection of the latter in our experiments. Global metagenome analyses show Clade BI and BII algal relative abundances correlate positively with their respective viruses. The distributions delineate BI/BpVs as occupying lower temperature mesotrophic and coastal systems, whereas BII/BII-Vs occupy warmer temperature, higher salinity ecosystems. Accordingly, with molecular diagnostic support, we name Clade BII Bathycoccus calidus sp. nov. and propose that molecular diversity within this new species likely connects to the differentiated host-virus dynamics observed in our time course experiments. Overall, the tightly linked biogeography of Bathycoccus host and virus clades observed herein supports species-level host specificity, with strain-level variations in infection parameters.

Evolution of pathogen tolerance and emerging infections: A missing experimental paradigm

Researchers worldwide are repeatedly warning us against future zoonotic diseases resulting from humankind's insurgence into natural ecosystems. The same zoonotic pathogens that cause severe infections in a human host frequently fail to produce any disease outcome in their natural hosts. What precise features of the immune system enable natural reservoirs to carry these pathogens so efficiently? To understand these effects, we highlight the importance of tracing the evolutionary basis of pathogen tolerance in reservoir hosts, while drawing implications from their diverse physiological and life-history traits, and ecological contexts of host-pathogen interactions. Long-term co-evolution might allow reservoir hosts to modulate immunity and evolve tolerance to zoonotic pathogens, increasing their circulation and infectious period. Such processes can also create a genetically diverse pathogen pool by allowing more mutations and genetic exchanges between circulating strains, thereby harboring rare alive-on-arrival variants with extended infectivity to new hosts (i.e., spillover). Finally, we end by underscoring the indispensability of a large multidisciplinary empirical framework to explore the proposed link between evolved tolerance, pathogen prevalence, and spillover in the wild.

Genomic Diversity of the Ostreid Herpesvirus Type 1 Across Time and Location and Among Host Species

The mechanisms underlying virus emergence are rarely well understood, making the appearance of outbreaks largely unpredictable. This is particularly true for pathogens with low per-site mutation rates, such as DNA viruses, that do not exhibit a large amount of evolutionary change among genetic sequences sampled at different time points. However, whole-genome sequencing can reveal the accumulation of novel genetic variation between samples, promising to render most, if not all, microbial pathogens measurably evolving and suitable for analytical techniques derived from population genetic theory. Here, we aim to assess the measurability of evolution on epidemiological time scales of the Ostreid herpesvirus 1 (OsHV-1), a double stranded DNA virus of which a new variant, OsHV-1 μVar, emerged in France in 2008, spreading across Europe and causing dramatic economic and ecological damage. We performed phylogenetic analyses of heterochronous (n = 21) OsHV-1 genomes sampled worldwide. Results show sufficient temporal signal in the viral sequences to proceed with phylogenetic molecular clock analyses and they indicate that the genetic diversity seen in these OsHV-1 isolates has arisen within the past three decades. OsHV-1 samples from France and New Zealand did not cluster together suggesting a spatial structuration of the viral populations. The genome-wide study of simple and complex polymorphisms shows that specific genomic regions are deleted in several isolates or accumulate a high number of substitutions. These contrasting and non-random patterns of polymorphism suggest that some genomic regions are affected by strong selective pressures. Interestingly, we also found variant genotypes within all infected individuals. Altogether, these results provide baseline evidence that whole genome sequencing could be used to study population dynamic processes of OsHV-1, and more broadly herpesviruses.

Novel genomic approaches to study antagonistic coevolution between hosts and parasites

Host-parasite coevolution is ubiquitous, shaping genetic and phenotypic diversity and the evolutionary trajectory of interacting species. With the advances of high throughput sequencing technologies applicable to model and non-model organisms alike, it is now feasible to study in greater detail (a) the genetic underpinnings of coevolution, (b) the speed and type of dynamics at coevolving loci, and (c) the genomic consequences of coevolution. This review focuses on three recently developed approaches that leverage information from host and parasite full genome data simultaneously to pinpoint coevolving loci and draw inference on the coevolutionary history. First, co-genome-wide association study (co-GWAS) methods allow pinpointing the loci underlying host-parasite interactions. These methods focus on detecting associations between genetic variants and the outcome of experimental infection tests or on correlations between genomes of naturally infected hosts and their infecting parasites. Second, extensions to population genomics methods can detect genes under coevolution and infer the coevolutionary history, such as fitness costs. Third, correlations between host and parasite population size in time are indicative of coevolution, and polymorphism levels across independent spatially distributed populations of hosts and parasites can reveal coevolutionary loci and infer coevolutionary history. We describe the principles of these three approaches and discuss their advantages and limitations based on coevolutionary theory. We present recommendations for their application to various host (prokaryotes, fungi, plants, and animals) and parasite (viruses, bacteria, fungi, and macroparasites) species. We conclude by pointing out methodological and theoretical gaps to be filled to extract maximum information from full genome data and thereby to shed light on the molecular underpinnings of coevolution.

Evolution, Ecology, and Zoonotic Transmission of Betacoronaviruses: A Review

Coronavirus infections have been a part of the animal kingdom for millennia. The difference emerging in the twenty-first century is that a greater number of novel coronaviruses are being discovered primarily due to more advanced technology and that a greater number can be transmitted to humans, either directly or via an intermediate host. This has a range of effects from annual infections that are mild to full-blown pandemics. This review compares the zoonotic potential and relationship between MERS, SARS-CoV, and SARS-CoV-2. The role of bats as possible host species and possible intermediate hosts including pangolins, civets, mink, birds, and other mammals are discussed with reference to mutations of the viral genome affecting zoonosis. Ecological, social, cultural, and environmental factors that may play a role in zoonotic transmission are considered with reference to SARS-CoV, MERS, and SARS-CoV-2 and possible future zoonotic events.

Shifts in global bat diversity suggest a possible role of climate change in the emergence of SARS-CoV-1 and SARS-CoV-2

Bats are the likely zoonotic origin of several coronaviruses (CoVs) that infect humans, including SARS-CoV-1 and SARS-CoV-2, both of which have caused large-scale epidemics. The number of CoVs present in an area is strongly correlated with local bat species richness, which in turn is affected by climatic conditions that drive the geographical distributions of species. Here we show that the southern Chinese Yunnan province and neighbouring regions in Myanmar and Laos form a global hotspot of climate change-driven increase in bat richness. This region coincides with the likely spatial origin of bat-borne ancestors of SARS-CoV-1 and SARS-CoV-2. Accounting for an estimated increase in the order of 100 bat-borne CoVs across the region, climate change may have played a key role in the evolution or transmission of the two SARS CoVs.

Genomic Study of COVID-19 Corona Virus Excludes Its Origin from Recombination or Characterized Biological Sources and Suggests a Role for HERVS in Its Wide Range Symptoms

The COVID-19 corona virus has become a world pandemic which started in December 2019 in Wuhan, China with no confirmed biological source. Various countries reported the genomic sequence of different isolates obtained from infected patients. This allowed us to obtain a number of 38 isolates of full genomic sequences. Alignment of nucleotide (nt) sequence was carried out using Clustal Omega multiple alignment service at the EBI website. Alignment of nt sequence and phylogenetic relationship revealed that the COVID-19 is a new viral strain and its biological source has not been yet detected. The expected orf pattern was different among isolates obtained from the same country or different countries as well as from SARS-CoV isolates or bats CoV suggesting different virus human interaction possibilities during infection and severity. All isolates had the main five orfs (1ab, S, M, N, E), whereas they differed in the expected accessory orfs. Being with the biological source of COVID-19 undetected, the role of human endogenous retrovirus (HERVs) in the regulation of the host cell gene expression or the encoding for products that could modulate COVID-19 infection and the spectrum of its symptoms is discussed.

Genetic Diversity and Evolution of Viral Populations

Population genetic diversity plays a prominent role in viral evolution, pathogenesis, immune escape, and drug resistance. Different mechanisms are responsible for creating and maintaining genetic diversity in viruses, including error-prone replication, repair avoidance, and genome editing, among others. This diversity is subsequently modulated by natural selection and random genetic drift, whose action in turn depends on multiple factors including viral genetic architecture, viral demography, and ecology. Understanding these processes should contribute to the development of more efficient control and treatment strategies against viral pathogens.

SARS-CoV 2; Possible alternative virus receptors and pathophysiological determinants

Understanding how severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) highjacks epithelial cells and infiltrates the lung, as well as other organs and tissues, is essential for developing treatment strategies and vaccines against this highly contagious virus. Another major goal is to fully elucidate the mechanisms by which SARS-CoV- 2 bypasses the innate immune system and induces a cytokine storm, and its effects on mortality. Currently, SARS- CoV-2 is thought to evade innate antiviral immunity, undergo endocytosis, and fuse with the host cell membrane by exploiting ACE2 receptors and the protease TMMPRSS2, with cathepsin B/L as alternative protease, for entry into the epithelial cells of tissues vulnerable to developing coronavirus disease 2019 (COVID-19) symptoms. However, the incorporation of new and unique binding sites, i.e., O-linked glycans, and the preservation and augmentation of effective binding sites (N-linked glycans) on the outer membrane of SARS-CoV-2 may represent other strategies of infecting the human host. Here, I will rationalize the possibility that other host molecules-i.e., sugar molecules and the sialic acidsN-glycolylneuraminic acid, N-acetylneuraminic acid, and their derivates could be viable candidates for the use as virus receptors by SARS-CoV-2 and/or serve as determinants for the adherence on ACE2 of SARS-CoV-2.

Evolutionary Genomics of High Fecundity

Natural highly fecund populations abound. These range from viruses to gadids. Many highly fecund populations are economically important. Highly fecund populations provide an important contrast to the low-fecundity organisms that have traditionally been applied in evolutionary studies. A key question regarding high fecundity is whether large numbers of offspring are produced on a regular basis, by few individuals each time, in a sweepstakes mode of reproduction. Such reproduction characteristics are not incorporated into the classical Wright-Fisher model, the standard reference model of population genetics, or similar types of models, in which each individual can produce only small numbers of offspring relative to the population size. The expected genomic footprints of population genetic models of sweepstakes reproduction are very different from those of the Wright-Fisher model. A key, immediate issue involves identifying the footprints of sweepstakes reproduction in genomic data. Whole-genome sequencing data can be used to distinguish the patterns made by sweepstakes reproduction from the patterns made by population growth in a population evolving according to the Wright-Fisher model (or similar models). If the hypothesis of sweepstakes reproduction cannot be rejected, then models of sweepstakes reproduction and associated multiple-merger coalescents will become at least as relevant as the Wright-Fisher model (or similar models) and the Kingman coalescent, the cornerstones of mathematical population genetics, in further discussions of evolutionary genomics of highly fecund populations.