Identification of shared populations of human immunodeficiency virus type 1 infecting microglia and tissue macrophages outside the central nervous system

A simple phylogenetic approach to analyze hypermutated HIV proviruses reveals insights into their dynamics and persistence during antiretroviral therapy

Hypermutated proviruses, which arise in a single Human Immunodeficiency Virus (HIV) replication cycle when host antiviral APOBEC3 proteins introduce extensive guanine to adenine mutations throughout the viral genome, persist in all people living with HIV receiving antiretroviral therapy (ART). However, hypermutated sequences are routinely excluded from phylogenetic trees because their extensive mutations complicate phylogenetic inference, and as a result, we know relatively little about their within-host evolutionary origins and dynamics. Using >1400 longitudinal single-genome-amplified HIV env-gp120 sequences isolated from six women over a median of 18 years of follow-up-including plasma HIV RNA sequences collected over a median of 9 years between seroconversion and ART initiation, and >500 proviruses isolated over a median of 9 years on ART-we evaluated three approaches for masking hypermutation in nucleotide alignments. Our goals were to (i) reconstruct phylogenies that can be used for molecular dating and (ii) phylogenetically infer the integration dates of hypermutated proviruses persisting during ART. Two of the approaches (stripping all positions containing putative APOBEC3 mutations from the alignment or replacing individual putative APOBEC3 mutations in hypermutated sequences with the ambiguous base R) consistently normalized tree topologies, eliminated erroneous clustering of hypermutated proviruses, and brought env-intact and hypermutated proviruses into comparable ranges with respect to multiple tree-based metrics. Importantly, these corrected trees produced integration date estimates for env-intact proviruses that were highly concordant with those from benchmark trees that excluded hypermutated sequences, supporting the use of these corrected trees for molecular dating. Subsequent molecular dating of hypermutated proviruses revealed that these sequences spanned a wide within-host age range, with the oldest ones dating to shortly after infection. This indicates that hypermutated proviruses, like other provirus types, begin to be seeded into the proviral pool immediately following infection and can persist for decades. In two of the six participants, hypermutated proviruses differed from env-intact ones in terms of their age distributions, suggesting that different provirus types decay at heterogeneous rates in some hosts. These simple approaches to reconstruct hypermutated provirus' evolutionary histories reveal insights into their in vivo origins and longevity toward a more comprehensive understanding of HIV persistence during ART.

Phylogenetic evidence of extensive spatial mixing of diverse HIV-1 group M lineages within Cameroon but not between its neighbours

From the perspective of developing relevant interventions for treating HIV and controlling its spread, it is particularly important to comprehensively understand the underlying diversity of the virus, especially in countries where the virus has been present and evolving since the cross-species transmission event that triggered the global pandemic. Here, we generate and phylogenetically analyse sequences derived from the gag-protease (2010 bp; n = 115), partial integrase (345 bp; n = 36), and nef (719 bp; n = 321) genes of HIV-1 group M (HIV-1M) isolates sampled between 2000 and 2022 from two cosmopolitan cities and 40 remote villages of Cameroon. While 52.4% of all sequenced viruses belonged to circulating recombinant form (CRF) 02_AG (CRF02_AG), the remainder were highly diverse, collectively representing seven subtypes and sub-subtypes, eight CRFs, and 36 highly divergent lineages that fall outside the established HIV-1M classification. Additionally, in 77 samples for which at least two genes were typed, 31% of the studied viruses apparently had fragments from viruses belonging to different clades. Furthermore, we found that the distribution of HIV-1M populations is similar between different regions of Cameroon. In contrast, HIV-1M demographics in Cameroon differ significantly from those in its neighbouring countries in the Congo Basin (CB). In phylogenetic trees, viral sequences cluster according to the countries where they were sampled, suggesting that while there are minimal geographical or social barriers to viral dissemination throughout Cameroon, there is strongly impeded dispersal of HIV-1M lineages between Cameroon and other locations of the CB. This suggests that the apparent stability of highly diverse Cameroonian HIV-1M populations may be attributable to the extensive mixing of human populations within the country and the concomitant trans-national movements of major lineages with very similar degrees of fitness; coupled with the relatively infrequent inter-national transmission of these lineages from neighbouring countries in the CB.

The origin and spread of HIV-1 CRF59_01B epidemic in China: A molecular network and phylogeographic analysis

Human immunodeficiency virus type 1 CRF59_01B, identified in China in 2013, has been detected nationwide, exhibiting notably high prevalence in Guangzhou and its vicinity. This study aimed to unravel its origin and migration. A data set was established, incorporating all available CRF59_01B pol gene sequences and their metadata from Guangzhou and the public database. Bayesian phylogeographic analysis demonstrated that CRF59_01B originated in Shenzhen, the neighboring city of Guangzhou, around 1998 with posterior probability of 0.937. Molecular network analysis detected 1131 transmission links and showed a remarkably high clustering rate (78.9%). Substantial inter-city transmissions (26.5%, 300/1131) were observed between Shenzhen and Guangzhou while inter-region transmissions linked Guangzhou with South (46) and Southwest (64) China. The centre of Guangzhou was the hub of CRF59_01B transmission, including the inflow from Shenzhen (3.57 events/year) and outflow to the outskirts of Guangzhou (>2 events/year). The large-scale analysis revealed significant migration from Shenzhen to Guangzhou (5.08 events/year) and North China (0.59 events/year), and spread from Guangzhou to Central (0.47 events/year), East (0.42 events/year), South (0.76 events/year), Southwest China (0.76 events/year) and Shenzhen (1.89 events/year). Shenzhen and Guangzhou served as the origin and the hub of CRF59_01B circulation, emphasizing inter-city cooperation and data sharing to confine its nationwide diffusion.

The replication-competent HIV reservoir is a genetically restricted, younger subset of the overall pool of HIV proviruses persisting during therapy, which is highly genetically stable over time

Within-host HIV populations continually diversify during untreated infection, and this diversity persists within infected cell reservoirs during antiretroviral therapy (ART). Achieving a better understanding of on-ART proviral evolutionary dynamics, and a better appreciation of how the overall persisting pool of (largely genetically defective) proviruses differs from the much smaller replication-competent HIV reservoir, is critical to HIV cure efforts. We reconstructed within-host HIV evolutionary histories in blood from seven participants of the Women's Interagency HIV Study who experienced HIV seroconversion, and used these data to characterize the diversity, lineage origins, and ages of proviral env-gp120 sequences sampled longitudinally up to 12 years on ART. We also studied HIV sequences emerging from the reservoir in two participants. We observed that proviral clonality generally increased over time on ART, with clones frequently persisting long term. While on-ART proviral integration dates generally spanned the duration of untreated infection, HIV emerging in plasma was exclusively younger (i.e., dated to the years immediately pre-ART). The genetic and age distributions of distinct proviral sequences remained stable during ART in all but one participant, in whom there was evidence that younger proviruses had been preferentially eliminated after 12 years on ART. Analysis of the gag region in three participants corroborated our env-gp120-based observations, indicating that our observations are not influenced by the HIV region studied. Our results underscore the remarkable genetic stability of the distinct proviral sequences that persist in blood during ART. Our results also suggest that the replication-competent HIV reservoir is a genetically restricted, younger subset of this overall proviral pool.IMPORTANCECharacterizing the genetically diverse HIV sequences that persist in the reservoir despite antiretroviral therapy (ART) is critical to cure efforts. Our observations confirm that proviruses persisting in blood on ART, which are largely genetically defective, broadly reflect the extent of within-host HIV evolution pre-ART. Moreover, on-ART clonal expansion is not appreciably accompanied by the loss of distinct proviral lineages. In fact, on-ART proviral genetic composition remained stable in all but one participant, in whom, after 12 years on ART, proviruses dating to around near ART initiation had been preferentially eliminated. We also identified recombinant proviruses between parental sequence fragments of different ages. Though rare, such sequences suggest that reservoir cells can be superinfected with HIV from another infection era. Overall, our finding that the replication-competent reservoir in blood is a genetically restricted, younger subset of all persisting proviruses suggests that HIV cure strategies will need to eliminate a reservoir that differs in key respects from the overall proviral pool.

Ecological Strategies for Resource Use by Three Bromoviruses in Anthropic and Wild Plant Communities

Ecological strategies for resource utilisation are important features of pathogens, yet have been overshadowed by stronger interest in genetic mechanisms underlying disease emergence. The purpose of this study is to ask whether host range and transmission traits translate into ecological strategies for host-species utilisation in a heterogeneous ecosystem, and whether host utilisation corresponds to genetic differentiation among three bromoviruses. We combine high-throughput sequencing and population genomics with analyses of species co-occurrence to unravel the ecological strategies of the viruses across four habitat types. The results show that the bromoviruses that were more closely related genetically did not share similar ecological strategies, but that the more distantly related pair did. Shared strategies included a broad host range and more frequent co-occurrences, which both were habitat-dependent. Each habitat thus presents as a barrier to gene flow, and each virus has an ecological strategy to navigate limitations to colonising non-natal habitats. Variation in ecological strategies could therefore hold the key to unlocking events that lead to emergence.

The replication-competent HIV reservoir is a genetically restricted, younger subset of the overall pool of HIV proviruses persisting during therapy, which is highly genetically stable over time

Within-host HIV populations continually diversify during untreated infection, and members of these diverse forms persist within infected cell reservoirs, even during antiretroviral therapy (ART). Characterizing the diverse viral sequences that persist during ART is critical to HIV cure efforts, but our knowledge of on-ART proviral evolutionary dynamics remains incomplete, as does our understanding of the differences between the overall pool of persisting proviral DNA (which is largely genetically defective) and the subset of intact HIV sequences capable of reactivating. Here, we reconstructed within-host HIV evolutionary histories in blood from seven participants of the Women's Interagency HIV Study (WIHS) who experienced HIV seroconversion. We measured diversity, lineage origins and ages of proviral sequences (env-gp120) sampled up to four times, up to 12 years on ART. We used the same techniques to study HIV sequences emerging from the reservoir in two participants. Proviral clonality generally increased over time on ART, with clones frequently persisting across multiple time points. The integration dates of proviruses persisting on ART generally spanned the duration of untreated infection (though were often skewed towards years immediately pre-ART), while in contrast, reservoir-origin viremia emerging in plasma was exclusively "younger" (i.e., dated to the years immediately pre-ART). The genetic and age distributions of distinct proviral sequences remained highly stable during ART in all but one participant in whom, after 12 years, there was evidence that "younger" proviruses had been preferentially eliminated. Analysis of within-host recombinant proviral sequences also suggested that HIV reservoirs can be superinfected with virus reactivated from an older era, yielding infectious viral progeny with mosaic genomes of sequences with different ages. Overall, results underscore the remarkable genetic stability of distinct proviral sequences that persist on ART, yet suggest that replication-competent HIV reservoir represents a genetically-restricted and overall "younger" subset of the overall persisting proviral pool in blood.

Population genetics and phylogeography of alfalfa mosaic virus in China and a comparison with other regional epidemics based on the cp gene

Alfalfa mosaic virus (AMV) is the most pervasive epidemic virus affecting alfalfa production. However, detailed investigations on the molecular population genetics and evolutionary dynamics of AMV are scarce. This study aimed to report on a large-scale long-term survey of genetic variability in AMV populations from China and perform a comparative analysis of AMV population genetics in the three most thoroughly studied countries to date: China, Iran, and Spain. The study was based on the analysis of the coat protein gene (cp) using two analytical approaches: an analysis of molecular variance (AMOVA) and a Bayesian Markov Chain Monte Carlo approach that investigates the association between geographical origin and phylogeny. Both analytical approaches found significant genetic differentiation within localities, but not among localities nor among provinces. This observation might result from inappropriate agronomical practices involving extensive exchange of plant materials followed by rapid viral diversification within localities. In the Chinese population, both methods found that genetic diversification in AMV was strongly associated with different bioclimatic zones. Rates of molecular evolution were similar in the three countries. The estimated epidemic exponential population size and growth rate suggest that the epidemics grew faster and with higher incidence in Iran, followed by Spain and China. Estimates of the time to the most recent common ancestors suggest that AMV was first seen in Spain by the beginning of the twentieth century and later on in eastern and central Eurasia. After ruling out the existence of recombination breakpoints within the cp gene, a codon-based selection analysis per population was performed and identified many codons under significant negative selection and a few under significant positive selection; the latter varied among countries, suggesting regional differences in selective pressures.

Persistence of envelopes in different CD4+ T-cell subsets in antiretroviral therapy-suppressed people with HIV

OBJECTIVES

Despite suppressive antiretroviral therapy (ART), HIV can persist in a diverse range of CD4+ T-cell subsets. Through longitudinal env sampling from people with HIV (PWH) on ART, we characterized the persistence and phenotypic properties of HIV envs over two time-points (T1 and T2).

METHODS

Longitudinal blood and lymphoid tissue samples were obtained from eight PWH on suppressive ART. Single genome amplification (SGA) was performed on env to understand the genetic diversity and degree of clonal expansions over time. A subset of envs were used to generate pseudovirus particles to assess sensitivity to autologous plasma IgG and broadly neutralizing antibodies (bNAbs).

RESULTS

Identical env sequences indicating clonal expansion persisted between T1 and T2 and within multiple T-cell subsets. At both time-points, CXCR4-tropic (X4) Envs were more prevalent in naive and central memory cells; the proportion of X4 Envs did not significantly change in each subset between T1 and T2. Autologous purified plasma IgG showed variable neutralization of Envs, with no significant difference in neutralization between R5 and X4 Envs. X4 Envs were more sensitive to neutralization with clinical bNAbs, with CD4-binding site bNAbs demonstrating high breadth and potency against Envs.

CONCLUSION

Our data suggest the viral reservoir in PWH on ART was predominantly maintained over time through proliferation and potentially differentiation of infected cells. We found the humoral immune response to Envs within the latent reservoir was variable between PWH. Finally, we identified coreceptor usage can influence bNAb sensitivity and may need to be considered for future bNAb immunotherapy approaches.

The diversity and evolutionary relationships of ticks and tick-borne bacteria collected in China

Background

Ticks (order Ixodida) are ectoparasites, vectors and reservoirs of many infectious agents affecting humans and domestic animals. However, the lack of information on tick genomic diversity leaves significant gaps in the understanding of the evolution of ticks and associated bacteria.

Results

We collected > 20,000 contemporary and historical (up to 60 years of preservation) tick samples representing a wide range of tick biodiversity across diverse geographic regions in China. Metagenomic sequencing was performed on individual ticks to obtain the complete or near-complete mitochondrial (mt) genome sequences from 46 tick species, among which mitochondrial genomes of 23 species were recovered for the first time. These new mt genomes data greatly expanded the diversity of many tick groups and revealed five cryptic species. Utilizing the same metagenomic sequence data we identified divergent and abundant bacteria in Haemaphysalis, Ixodes, Dermacentor and Carios ticks, including nine species of pathogenetic bacteria and potentially new species within the genus Borrelia. We also used these data to explore the evolutionary relationship between ticks and their associated bacteria, revealing a pattern of long-term co-divergence relationship between ticks and Rickettsia and Coxiella bacteria.

Conclusions

In sum, our study provides important new information on the genetic diversity of ticks based on an analysis of mitochondrial DNA as well as on the prevalence of tick-borne pathogens in China. It also sheds new light on the long-term evolutionary and ecological relationships between ticks and their associated bacteria.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05485-3.

Discordance between HIV-1 Population in Plasma at Rebound after Structured Treatment Interruption and Archived Provirus Population in Peripheral Blood Mononuclear Cells

Antiretroviral therapy (ART) can sustain the suppression of plasma viremia to below detection levels. Infected individuals undergoing a treatment interruption exhibit rapid viral rebound in plasma viremia which is fueled by cellular reservoirs such as CD4+ T cells, myeloid cells, and potentially uncharacterized cellular sources. Interrogating the populations of viruses found during analytical treatment interruption (ATI) can give insights into the biologically competent reservoirs that persist under effective ART as well as the nature of the cellular reservoirs that enable viral persistence under ART. We interrogated plasma viremia from four rare cases of individuals undergoing sequential ATIs. We performed next-generation sequencing (NGS) on cell-associated viral DNA and cell-free virus to understand the interrelationship between sequential ATIs as well as the relationship between viral genomes in circulating peripheral blood mononuclear cells (PBMCs) and RNA from rebound plasma. We observed population differences between viral populations recrudescing at sequential ATIs as well as divergence between viral sequences in plasma and those in PBMCs. This indicated that viruses in PBMCs were not a major source of post-ATI viremia and highlights the role of anatomic reservoirs in post-ATI viremia and viral persistence. IMPORTANCE Even with effective ART, HIV-1 persists at undetectable levels and rebounds in individuals who stop treatment. Cellular and anatomical reservoirs ignite viral rebound upon treatment interruption, remaining one of the key obstacles for HIV-1 cure. To further examine HIV-1 persistence, a better understanding of the distinct populations that fuel viral rebound is necessary to identify and target reservoirs and the eradication of HIV-1. This study investigates the populations of viruses found from proviral genomes from PBMCs and plasma at rebound from a unique cohort of individuals who underwent multiple rounds of treatment interruption. Using NGS, we characterized the subtypes of viral sequences and found divergence in viral populations between plasma and PBMCs at each rebound, suggesting that distinct viral populations appear at each treatment interruption.

Statistical Challenges in Tracking the Evolution of SARS-CoV-2

Genomic surveillance of SARS-CoV-2 has been instrumental in tracking the spread and evolution of the virus during the pandemic. The availability of SARS-CoV-2 molecular sequences isolated from infected individuals, coupled with phylodynamic methods, have provided insights into the origin of the virus, its evolutionary rate, the timing of introductions, the patterns of transmission, and the rise of novel variants that have spread through populations. Despite enormous global efforts of governments, laboratories, and researchers to collect and sequence molecular data, many challenges remain in analyzing and interpreting the data collected. Here, we describe the models and methods currently used to monitor the spread of SARS-CoV-2, discuss long-standing and new statistical challenges, and propose a method for tracking the rise of novel variants during the epidemic.

A Daily-Updated Database and Tools for Comprehensive SARS-CoV-2 Mutation-Annotated Trees

The vast scale of SARS-CoV-2 sequencing data has made it increasingly challenging to comprehensively analyze all available data using existing tools and file formats. To address this, we present a database of SARS-CoV-2 phylogenetic trees inferred with unrestricted public sequences, which we update daily to incorporate new sequences. Our database uses the recently proposed mutation-annotated tree (MAT) format to efficiently encode the tree with branches labeled with parsimony-inferred mutations, as well as Nextstrain clade and Pango lineage labels at clade roots. As of June 9, 2021, our SARS-CoV-2 MAT consists of 834,521 sequences and provides a comprehensive view of the virus' evolutionary history using public data. We also present matUtils-a command-line utility for rapidly querying, interpreting, and manipulating the MATs. Our daily-updated SARS-CoV-2 MAT database and matUtils software are available at http://hgdownload.soe.ucsc.edu/goldenPath/wuhCor1/UShER_SARS-CoV-2/ and https://github.com/yatisht/usher, respectively.

Rotavirus A Genome Segments Show Distinct Segregation and Codon Usage Patterns

Reassortment of the Rotavirus A (RVA) 11-segment dsRNA genome may generate new genome constellations that allow RVA to expand its host range or evade immune responses. Reassortment may also produce phylogenetic incongruities and weakly linked evolutionary histories across the 11 segments, obscuring reassortment-specific epistasis and changes in substitution rates. To determine the co-segregation patterns of RVA segments, we generated time-scaled phylogenetic trees for each of the 11 segments of 789 complete RVA genomes isolated from mammalian hosts and compared the segments’ geodesic distances. We found that segments 4 (VP4) and 9 (VP7) occupied significantly different tree spaces from each other and from the rest of the genome. By contrast, segments 10 and 11 (NSP4 and NSP5/6) occupied nearly indistinguishable tree spaces, suggesting strong co-segregation. Host-species barriers appeared to vary by segment, with segment 9 (VP7) presenting the weakest association with host species. Bayesian Skyride plots were generated for each segment to compare relative genetic diversity among segments over time. All segments showed a dramatic decrease in diversity around 2007 coinciding with the introduction of RVA vaccines. To assess selection pressures, codon adaptation indices and relative codon deoptimization indices were calculated with respect to different host genomes. Codon usage varied by segment with segment 11 (NSP5) exhibiting significantly higher adaptation to host genomes. Furthermore, RVA codon usage patterns appeared optimized for expression in humans and birds relative to the other hosts examined, suggesting that translational efficiency is not a barrier in RVA zoonosis.

A daily-updated database and tools for comprehensive SARS-CoV-2 mutation-annotated trees

The vast scale of SARS-CoV-2 sequencing data has made it increasingly challenging to comprehensively analyze all available data using existing tools and file formats. To address this, we present a database of SARS-CoV-2 phylogenetic trees inferred with unrestricted public sequences, which we update daily to incorporate new sequences. Our database uses the recently-proposed mutation-annotated tree (MAT) format to efficiently encode the tree with branches labeled with parsimony-inferred mutations as well as Nextstrain clade and Pango lineage labels at clade roots. As of June 9, 2021, our SARS-CoV-2 MAT consists of 834,521 sequences and provides a comprehensive view of the virus’ evolutionary history using public data. We also present matUtils – a command-line utility for rapidly querying, interpreting and manipulating the MATs. Our daily-updated SARS-CoV-2 MAT database and matUtils software are available at http://hgdownload.soe.ucsc.edu/goldenPath/wuhCor1/UShER_SARS-CoV-2/ and https://github.com/yatisht/usher, respectively.

The Molecular Epidemiology and Clinical Phylogenetics of Rhinoviruses Among Paediatric Cases in Sydney, Australia

OBJECTIVES

Rhinoviruses (RV) represent the most common aetiological agent of all acute respiratory tract infections across all age groups and a significant burden of disease among children. Recent studies have shown that RV-A and RV-C species are associated with increased disease severity. In order to better understand the potential associations between RV species and clinical features among paediatric cases, this study aimed to integrate genetic and epidemiological data using Bayesian phylogenetic methods.

METHODS

Potential associations between RV species and subtypes, and clinical disease severity using a matched dataset of 52 RV isolates sampled from children (< 18 years) in Sydney, Australia, between 2006 and 2009 were uncovered using epidemiological and phylogenetic methods.

RESULTS

It was found that RV-C was significantly more likely to be isolated from paediatric cases aged < 2 years compared with RV-A, although no significant differences in recorded symptoms were observed. Significant phylogenetic-trait associations between age and the VP4/VP2 capsid protein phylogeny suggest that age-specific variations in infectivity among subtypes may may be possible.

CONCLUSION

This study adds to the growing body of epidemiological evidence concerning RV. Improving surveillance and testing for RV, including routine whole genome sequencing, may improve understanding of the varied disease outcomes of RV species and subtypes. Future studies could aim to identify specific genetic markers associated with age-specific infectivity of RV, which could inform treatment practices and public health surveillance of RV.

St. Louis Encephalitis Virus in the Southwestern United States: A Phylogeographic Case for a Multi-Variant Introduction Event

Since the reemergence of St. Louis Encephalitis (SLE) Virus (SLEV) in the Southwest United States, identified during the 2015 outbreak in Arizona, SLEV has been seasonally detected within Culex spp. populations throughout the Southwest United States. Previous work revealed the 2015 outbreak was caused by an importation of SLEV genotype III, which had only been detected previously in Argentina. However, little is known about when the importation occurred or the transmission and genetic dynamics since its arrival into the Southwest. In this study, we sought to determine whether the annual detection of SLEV in the Southwest is due to enzootic cycling or new importations. To address this question, we analyzed 174 SLEV genomes (142 sequenced as part of this study) using Bayesian phylogenetic analyses to estimate the date of arrival into the American Southwest and characterize the underlying population structure of SLEV. Phylogenetic clustering showed that SLEV variants circulating in Maricopa and Riverside counties form two distinct populations with little evidence of inter-county transmission since the onset of the outbreak. Alternatively, it appears that in 2019, Yuma and Clark counties experienced annual importations of SLEV that originated in Riverside and Maricopa counties. Finally, the earliest representatives of SLEV genotype III in the Southwest form a polytomy that includes both California and Arizona samples. We propose that the initial outbreak most likely resulted from the importation of a population of SLEV genotype III variants, perhaps in multiple birds, possibly multiple species, migrating north in 2013, rather than a single variant introduced by one bird.

Extensive C->U transition biases in the genomes of a wide range of mammalian RNA viruses; potential associations with transcriptional mutations, damage- or host-mediated editing of viral RNA

The rapid evolution of RNA viruses has been long considered to result from a combination of high copying error frequencies during RNA replication, short generation times and the consequent extensive fixation of neutral or adaptive changes over short periods. While both the identities and sites of mutations are typically modelled as being random, recent investigations of sequence diversity of SARS coronavirus 2 (SARS-CoV-2) have identified a preponderance of C->U transitions, proposed to be driven by an APOBEC-like RNA editing process. The current study investigated whether this phenomenon could be observed in datasets of other RNA viruses. Using a 5% divergence filter to infer directionality, 18 from 36 datasets of aligned coding region sequences from a diverse range of mammalian RNA viruses (including Picornaviridae, Flaviviridae, Matonaviridae, Caliciviridae and Coronaviridae) showed a >2-fold base composition normalised excess of C->U transitions compared to U->C (range 2.1x-7.5x), with a consistently observed favoured 5' U upstream context. The presence of genome scale RNA secondary structure (GORS) was the only other genomic or structural parameter significantly associated with C->U/U->C transition asymmetries by multivariable analysis (ANOVA), potentially reflecting RNA structure dependence of sites targeted for C->U mutations. Using the association index metric, C->U changes were specifically over-represented at phylogenetically uninformative sites, potentially paralleling extensive homoplasy of this transition reported in SARS-CoV-2. Although mechanisms remain to be functionally characterised, excess C->U substitutions accounted for 11-14% of standing sequence variability of structured viruses and may therefore represent a potent driver of their sequence diversification and longer-term evolution.

SARS-CoV-2 within-host diversity and transmission

Extensive global sampling and sequencing of the pandemic virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have enabled researchers to monitor its spread and to identify concerning new variants. Two important determinants of variant spread are how frequently they arise within individuals and how likely they are to be transmitted. To characterize within-host diversity and transmission, we deep-sequenced 1313 clinical samples from the United Kingdom. SARS-CoV-2 infections are characterized by low levels of within-host diversity when viral loads are high and by a narrow bottleneck at transmission. Most variants are either lost or occasionally fixed at the point of transmission, with minimal persistence of shared diversity, patterns that are readily observable on the phylogenetic tree. Our results suggest that transmission-enhancing and/or immune-escape SARS-CoV-2 variants are likely to arise infrequently but could spread rapidly if successfully transmitted.

HIV RNA Rebound in Seminal Plasma after Antiretroviral Treatment Interruption

If strategies currently in development succeed in eradicating HIV reservoirs in peripheral blood and lymphoid tissues, residual sources of virus may remain in anatomic compartments. Paired blood and semen samples were collected from 12 individuals enrolled in a randomized, double-blind, placebo-controlled therapeutic vaccine clinical trial in people with HIV (PWH) who began antiretroviral therapy (ART) during acute or early infection (ClinicalTrials registration no. NCT01859325). After the week 56 visit (postintervention), all participants interrupted ART. At the first available time points after viral rebound, we sequenced HIV-1 env (C2-V3), gag (p24), and pol (reverse transcriptase) regions amplified from cell-free HIV RNA in blood and seminal plasma using the MiSeq Illumina platform. Comprehensive sequence and phylogenetic analyses were performed to evaluate viral population structure, compartmentalization, and viral diversity in blood and seminal plasma. Compared to that in blood, HIV RNA rebound in semen occurred significantly later (median of 66 versus 42 days post-ART interruption, P < 0.01) and reached lower levels (median 164 versus 16,090 copies/ml, P < 0.01). Three of five participants with available sequencing data presented compartmentalized viral rebound between blood and semen in one HIV coding region. Despite early ART initiation, HIV RNA molecular diversity was higher in semen than in blood in all three coding regions for most participants. Higher HIV RNA molecular diversity in the genital tract (compared to that in blood plasma) and evidence of compartmentalization illustrate the distinct evolutionary dynamics between these two compartments after ART interruption. Future research should evaluate whether the genital compartment might contribute to viral rebound in some PWH interrupting ART.IMPORTANCE To cure HIV, we likely need to target the reservoirs in all anatomic compartments. Here, we used sophisticated statistical and phylogenetic methods to analyze blood and semen samples collected from 12 persons with HIV who began antiretroviral therapy (ART) during very early HIV infection and who interrupted their ART as part of a clinical trial. First, we found that HIV RNA rebound in semen occurred significantly later and reached lower levels than in blood. Second, we found that the virus in semen was genetically different in some participants compared to that in blood. Finally, we found increased HIV RNA molecular diversity in semen compared to that in blood in almost all study participants. These data suggest that the HIV RNA populations emerging from the genital compartment after ART interruption might not be the same as those emerging from blood plasma. Future research should evaluate whether the genital compartment might contribute to viral rebound in some people with HIV (PWH) interrupting ART.

The Glutamate System as a Crucial Regulator of CNS Toxicity and Survival of HIV Reservoirs

Glutamate (Glu) is the most abundant excitatory neurotransmitter in the central nervous system (CNS). HIV-1 and viral proteins compromise glutamate synaptic transmission, resulting in poor cell-to-cell signaling and bystander toxicity. In this study, we identified that myeloid HIV-1-brain reservoirs survive in Glu and glutamine (Gln) as a major source of energy. Thus, we found a link between synaptic compromise, metabolomics of viral reservoirs, and viral persistence. In the current manuscript we will discuss all these interactions and the potential to achieve eradication and cure using this unique metabolic profile.

HIV persists throughout deep tissues with repopulation from multiple anatomical sources

BACKGROUNDUnderstanding HIV dynamics across the human body is important for cure efforts. This goal has been hampered by technical difficulties and the challenge of obtaining fresh tissues.METHODSThis observational study evaluated 6 individuals with HIV (n = 4 with viral suppression using antiretroviral [ART] therapy; n = 2 with rebound viremia after stopping ART), who provided serial blood samples before death and their bodies for rapid autopsy. HIV reservoirs were characterized by digital droplet PCR, single-genome amplification, and sequencing of full-length (FL) envelope HIV. Phylogeographic methods were used to reconstruct HIV spread, and generalized linear models were tested for viral factors associated with dispersal.RESULTSAcross participants, HIV DNA levels varied from approximately 0 to 659 copies/106 cells (IQR: 22.9-126.5). A total of 605 intact FL env sequences were recovered in antemortem blood cells and across 28 tissues (IQR: 5-9). Sequence analysis showed (a) the emergence of large, identical, intact HIV RNA populations in blood after cessation of therapy, which repopulated tissues throughout the body; (b) that multiple sites acted as hubs for HIV dissemination but that blood and lymphoid tissues were the main source; (c) that viral exchanges occurred within brain areas and across the blood-brain barrier; and (d) that migration was associated with low HIV divergence between sites and greater diversity at the recipient site.CONCLUSIONHIV reservoirs persisted in all deep tissues, and blood was the main source of dispersal. This may explain why eliminating HIV susceptibility in circulating T cells via bone marrow transplants allowed some individuals with HIV to experience therapy-free remission, even though deeper tissue reservoirs were not targeted.TRIAL REGISTRATIONNot applicable.FUNDINGNIH grants P01 AI31385, P30 AI036214, AI131971-01, AI120009AI036214, HD094646, AI027763, AI134295, and AI68636.

Genetic Strain Diversity of Multi-Host RNA Viruses that Infect a Wide Range of Pollinators and Associates is Shaped by Geographic Origins

Emerging viruses have caused concerns about pollinator population declines, as multi-host RNA viruses may pose a health threat to pollinators and associated arthropods. In order to understand the ecology and impact these viruses have, we studied their host range and determined to what extent host and spatial variation affect strain diversity. Firstly, we used RT-PCR to screen pollinators and associates, including honey bees (Apis mellifera) and invasive Argentine ants (Linepithema humile), for virus presence and replication. We tested for the black queen cell virus (BQCV), deformed wing virus (DWV), and Kashmir bee virus (KBV) that were initially detected in bees, and the two recently discovered Linepithema humile bunya-like virus 1 (LhuBLV1) and Moku virus (MKV). DWV, KBV, and MKV were detected and replicated in a wide range of hosts and commonly co-infected hymenopterans. Secondly, we placed KBV and DWV in a global phylogeny with sequences from various countries and hosts to determine the association of geographic origin and host with shared ancestry. Both phylogenies showed strong geographic rather than host-specific clustering, suggesting frequent inter-species virus transmission. Transmission routes between hosts are largely unknown. Nonetheless, avoiding the introduction of non-native species and diseased pollinators appears important to limit spill overs and disease emergence.

Tracing the Impact of Public Health Interventions on HIV-1 Transmission in Portugal Using Molecular Epidemiology

BACKGROUND

Estimation of temporal changes in human immunodeficiency virus (HIV) transmission patterns can help to elucidate the impact of preventive strategies and public health policies.

METHODS

Portuguese HIV-1 subtype B and G pol genetic sequences were appended to global reference data sets to identify country-specific transmission clades. Bayesian birth-death models were used to estimate subtype-specific effective reproductive numbers (Re). Discrete trait analysis (DTA) was used to quantify mixing among transmission groups.

RESULTS

We identified 5 subtype B Portuguese clades (26-79 sequences) and a large monophyletic subtype G Portuguese clade (236 sequences). We estimated that major shifts in HIV-1 transmission occurred around 1999 (95% Bayesian credible interval [BCI], 1998-2000) and 2000 (95% BCI, 1998-2001) for subtypes B and G, respectively. For subtype B, Re dropped from 1.91 (95% BCI, 1.73-2.09) to 0.62 (95% BCI,.52-.72). For subtype G, Re decreased from 1.49 (95% BCI, 1.39-1.59) to 0.72 (95% BCI, .63-.8). The DTA suggests that people who inject drugs (PWID) and heterosexuals were the source of most (>80%) virus lineage transitions for subtypes G and B, respectively.

CONCLUSIONS

The estimated declines in Re coincide with the introduction of highly active antiretroviral therapy and the scale-up of harm reduction for PWID. Inferred transmission events across transmission groups emphasize the importance of prevention efforts for bridging populations.

Swift evolutionary response of microbes to a rise in anthropogenic mercury in the Northern Hemisphere

Anthropogenic mercury remobilization has considerably increased since the Industrial Revolution in the late 1700s. The Minamata Convention on Mercury is a United Nations treaty (2017) aiming at curbing mercury emissions. Unfortunately, evaluating the effectiveness of such a global treaty is hampered by our inability to determine the lag in aquatic ecosystem responses to a change in atmospheric mercury deposition. Whereas past metal concentrations are obtained from core samples, there are currently no means of tracking historical metal bioavailability or toxicity. Here, we recovered DNA from nine dated sediment cores collected in Canada and Finland, and reconstructed the past demographics of microbes carrying genes coding for the mercuric reductase (MerA)—an enzyme involved in Hg detoxification—using Bayesian relaxed molecular clocks. We found that the evolutionary dynamics of merA exhibited a dramatic increase in effective population size starting from 1783.8 ± 3.9 CE, which coincides with both the Industrial Revolution, and with independent measurements of atmospheric Hg concentrations. We show that even low levels of anthropogenic mercury affected the evolutionary trajectory of microbes in the Northern Hemisphere, and that microbial DNA encoding for detoxification determinants stored in environmental archives can be used to track historical pollutant toxicity.

Distinct rates and patterns of spread of the major HIV-1 subtypes in Central and East Africa

Since the ignition of the HIV-1 group M pandemic in the beginning of the 20th century, group M lineages have spread heterogeneously throughout the world. Subtype C spread rapidly through sub-Saharan Africa and is currently the dominant HIV lineage worldwide. Yet the epidemiological and evolutionary circumstances that contributed to its epidemiological expansion remain poorly understood. Here, we analyse 346 novel pol sequences from the DRC to compare the evolutionary dynamics of the main HIV-1 lineages, subtypes A1, C and D. Our results place the origins of subtype C in the 1950s in Mbuji-Mayi, the mining city of southern DRC, while subtypes A1 and D emerged in the capital city of Kinshasa, and subtypes H and J in the less accessible port city of Matadi. Following a 15-year period of local transmission in southern DRC, we find that subtype C spread at least three-fold faster than other subtypes circulating in Central and East Africa. In conclusion, our results shed light on the origins of HIV-1 main lineages and suggest that socio-historical rather than evolutionary factors may have determined the epidemiological fate of subtype C in sub-Saharan Africa.

HIV Rebound Is Predominantly Fueled by Genetically Identical Viral Expansions from Diverse Reservoirs

Viral rebound upon stopping combined antiretroviral therapy poses a major barrier toward an HIV cure. Cellular and anatomical sources responsible for reinitiating viral replication remain a subject of ardent debate, despite extensive research efforts. To unravel the source of rebounding viruses, we conducted a large-scale HIV-STAR (HIV-1 sequencing before analytical treatment interruption to identify the anatomically relevant HIV reservoir) clinical trial. We collected samples from 11 participants and compared the genetic composition of (pro)viruses collected under treatment from different cellular and anatomical compartments with that of plasma viruses sampled during analytical treatment interruption. We found a remarkably heterogeneous source of viral rebound. In addition, irrespective of the compartment or cell subset, genetically identical viral expansions played a significant role in viral rebound. Our study suggests that although there does not seem to be a primary source for rebound HIV, cellular proliferation is an important driver of HIV persistence and should therefore be considered in future curative strategies.

Feline coronavirus isolates from a part of Brazil: insights into molecular epidemiology and phylogeny inferred from the 7b gene

The Feline coronavirus (FCoV) can lead to Feline infectious peritonitis (FIP), which the precise cause is still unknown. The theory of internal mutation suggests that a less virulent biotype of FCoV (FECV) would lead to another more pathogenic biotype (FIPV) capable of causing FIP. In this work, the 7b gene was amplified from 51 domestic cat plasma samples by semi-nested PCR and tested through phylogenetic and phylogeographical approaches. The 7b gene of Brazilian isolates displayed high conservation, a strong correlation between the geographic origin of the viral isolates and their genealogy, and its evolution was possibly shaped by a combination of high rates of nucleotide substitution and purifying selection.

Phylogeographic Analyses Reveal the Early Expansion and Frequent Bidirectional Cross-Border Transmissions of Non-pandemic HIV-1 Subtype B Strains in Hispaniola

The human immunodeficiency virus-type 1 (HIV-1) subtype B has probably been circulating on the island of Hispaniola since the 1960s, but information about the early viral history on this Caribbean island is scarce. In this study, we reconstruct the dissemination dynamics of early divergent non-pandemic subtype B lineages (designated B_CAR) on Hispaniola by analyzing a country-balanced dataset of HIV-1 B_CARpol sequences from Haiti (n = 103) and the Dominican Republic (n = 123). Phylogenetic analyses supported that B_CAR strains from Haiti and the Dominican Republic were highly intermixed between each other, although the null hypothesis of completely random mixing was rejected. Bayesian phylogeographic analyses placed the ancestral B_CAR virus in Haiti and the Dominican Republic with the same posterior probability support. These analyses estimate frequent viral transmissions between Haiti and the Dominican Republic since the early 1970s onwards, and the presence of local B_CAR transmission networks in both countries before first AIDS cases was officially recognized. Demographic reconstructions point that the B_CAR epidemic in Hispaniola grew exponentially until the 1990s. These findings support that the HIV-1 epidemics in Haiti and the Dominican Republic have been connected by a recurrent bidirectional viral flux since the initial phase, which poses a great challenge in tracing the geographic origin of the B_CAR epidemic within Hispaniola using only genetic data. These data also reinforce the notion that prevention programs have successfully reduced the rate of new HIV-1 transmissions in Hispaniola since the end of the 1990s.

Virus dynamics and phyloanatomy: Merging population dynamic and phylogenetic approaches

In evolutionary biology and epidemiology, phylodynamic methods are widely used to infer population biological characteristics, such as the rates of replication, death, migration, or, in the epidemiological context, pathogen spread. More recently, these methods have been used to elucidate the dynamics of viruses within their hosts. Especially the application of phylogeographic approaches has the potential to shed light on anatomical colonization pathways and the exchange of viruses between distinct anatomical compartments. We and others have termed this phyloanatomy. Here, we review the promise and challenges of phyloanatomy, and compare them to more classical virus dynamics and population genetic approaches. We argue that the extremely strong selection pressures that exist within the host may represent the main obstacle to reliable phyloanatomic analysis.

Phylogenetic inference for the study of within-host HIV-1 dynamics and persistence on antiretroviral therapy

Although antiretroviral therapy (ART) is highly effective at inhibiting HIV-1 replication and preventing AIDS, it cannot eradicate the infection. Many studies have used viral genetic information from single-genome and deep sequencing of blood and tissue samples to investigate the mechanisms that sustain the HIV-1 reservoir. Sequence data are analysed by use of measurements of population diversity and divergence and by exploration of phylogenetic associations. The study of intrahost HIV-1 populations on ART requires specific considerations as their dynamics can be shaped by host factors such as cell death and proliferation. Hence, understanding both the biology of HIV-1 persistence and the phylogenetic methods that can be applied to this field is crucial. We conclude that the most suitable phylogenetic methods and evolutionary models for characterising HIV-1 populations on ART include using neighbour-joining trees to identify identical proviral sequences that might result from T-cell proliferation, and using maximum-likelihood analysis to investigate the possibility of ongoing viral replication on ART. Characterising the reservoir for HIV-1 on ART is a high priority for the design of curative interventions.

Intrahost Dynamics of Human Cytomegalovirus Variants Acquired by Seronegative Glycoprotein B Vaccinees

Human cytomegalovirus (HCMV) is the most common congenital infection worldwide and a frequent cause of hearing loss and debilitating neurologic disease in newborn infants. Thus, a vaccine to prevent HCMV-associated congenital disease is a public health priority. One potential strategy is vaccination of women of child bearing age to prevent maternal HCMV acquisition during pregnancy. The glycoprotein B (gB) plus MF59 adjuvant subunit vaccine is the most efficacious tested clinically to date, demonstrating 50% protection against primary HCMV infection in a phase 2 clinical trial. Yet, the impact of gB/MF59-elicited immune responses on the population of viruses acquired by trial participants has not been assessed. In this analysis, we employed quantitative PCR as well as multiple sequencing methodologies to interrogate the magnitude and genetic composition of HCMV populations infecting gB/MF59 vaccinees and placebo recipients. We identified several differences between the viral dynamics in acutely infected vaccinees and placebo recipients. First, viral load was reduced in the saliva of gB vaccinees, though not in whole blood, vaginal fluid, or urine. Additionally, we observed possible anatomic compartmentalization of gB variants in the majority of vaccinees compared to only a single placebo recipient. Finally, we observed reduced acquisition of genetically related gB1, gB2, and gB4 genotype "supergroup" HCMV variants among vaccine recipients, suggesting that the gB1 genotype vaccine construct may have elicited partial protection against HCMV viruses with antigenically similar gB sequences. These findings suggest that gB immunization had a measurable impact on viral intrahost population dynamics and support future analysis of a larger cohort.IMPORTANCE Though not a household name like Zika virus, human cytomegalovirus (HCMV) causes permanent neurologic disability in one newborn child every hour in the United States, which is more than that for Down syndrome, fetal alcohol syndrome, and neural tube defects combined. There are currently no established effective measures to prevent viral transmission to the infant following HCMV infection of a pregnant mother. However, the glycoprotein B (gB)/MF59 vaccine, which aims to prevent pregnant women from acquiring HCMV, is the most successful HCMV vaccine tested clinically to date. Here, we used viral DNA isolated from patients enrolled in a gB vaccine trial who acquired HCMV and identified several impacts that this vaccine had on the size, distribution, and composition of the in vivo viral population. These results have increased our understanding of why the gB/MF59 vaccine was partially efficacious, and such investigations will inform future rational design of a vaccine to prevent congenital HCMV.

Phylogenetic and morphological analyses of the mycoparasitic genus Piptocephalis

The Piptocephalidaceae (Zoopagales, Zoopagomycota) contains three genera of mycoparasitic, haustoria-forming fungi: Kuzuhaea, Piptocephalis, and Syncephalis. Although the species in this family are diverse and ubiquitous in soil and dung, they are among the least studied fungi. Co-cultures of Piptocephalis and their hosts are relatively easy to isolate from soil and dung samples across the globe, making them a good model taxon for the order Zoopagales. This study focuses on the systematics of the genus Piptocephalis. Despite the fact that there are approximately 40 described Piptocephalis species, there are no modern taxonomic or molecular phylogenetic treatments of this group. Minimal sequence data are available, and relatively little is known about the true diversity or biogeography of the genus. Our study addresses two aspects: Piptocephalis systematics and analyses of the length and inter- and infraspecific variation of the nuc rDNA internal transcribed spacer (ITS1-5.8S-ITS2 = ITS) region. First, we generated a large subunit (28S) nuc rDNA phylogeny and evaluated several morphological characters by testing their correlation with the phylogeny using Bayesian Tip-association Significance testing (BaTS). We found monophyly of Piptocephalis species identified based on morphological traits, but morphological character states were not conserved across clades, suggesting that there have been multiple gains and losses of morphological characters. We also found that Kuzhuaea is nested within Piptocephalis. Second, we amplified the ITS from many Piptocephalis isolates, created a sequence alignment, and measured the lengths using the software ITSx. Piptocephalis species had ITS regions that were longer than the average for most Dikarya but were similar in length to those of the related genus Syncephalis.

The Grass was Greener: Repeated Evolution of Specialized Morphologies and Habitat Shifts in Ghost Spiders Following Grassland Expansion in South America

While grasslands, one of Earth's major biomes, are known for their close evolutionary ties with ungulate grazers, these habitats are also paramount to the origins and diversification of other animals. Within the primarily South American spider subfamily Amaurobioidinae (Anyphaenidae), several species are found living in the continent's grasslands, with some displaying putative morphological adaptations to dwelling unnoticed in the grass blades. Herein, a dated molecular phylogeny provides the backbone for analyses revealing the ecological and morphological processes behind these spiders' grassland adaptations. The multiple switches from Patagonian forests to open habitats coincide with the expansion of South America's grasslands during the Miocene, while the specialized morphology of several grass-dwelling spiders originated at least three independent times and is best described as the result of different selective regimes operating on macroevolutionary timescales. Although grass-adapted lineages evolved towards different peaks in adaptive landscape, they all share one characteristic: an anterior narrowing of the prosoma allowing spiders to extend the first two pairs of legs, thus maintaining a slender resting posture in the grass blade. By combining phylogenetic, morphological, and biogeographic perspectives we disentangle multiple factors determining the evolution of a clade of terrestrial invertebrate predators alongside their biomes.

Geography but not alternative host species explain the spread of raccoon rabies virus in Vermont

In North America, the raccoon-associated variant of rabies virus (RRV) is of special concern, given its relatively rapid spread throughout the eastern USA and its potential public health impact due to high raccoon host densities in urban areas. Northward expansion of this epizootic included an outbreak in the Canadian province of Quebec in 2006-2009 due to trans-border spread from the State of Vermont. To inform a more proactive approach to future control efforts, this study uses phylogenetic analyses to explore the role of geography and alternative carnivore hosts in the dynamics of RRV spread within Vermont. Specifically, we sought to examine whether striped skunks, a species frequently infected by RRV, could be part of the maintenance host community. Whole genome sequencing of 160 RRV samples from Vermont and neighbouring US states were used for fine-scale phylogeographic analyses. Results, together with the complete surveillance record of raccoon rabies since its entry into Vermont in 1994, document incursions by two distinct viral lineages and identify topographical features of the landscape which have significantly influenced viral spread, resulting in a complex distribution pattern of viral variants throughout the state. Results of phylogenetic cluster analysis and discrete state reconstruction contained some evidence of skunk-to-skunk and skunk-to-raccoon transmission but overall failed to support a role for skunks as alternative maintenance hosts.

The evolutionary history of vertebrate RNA viruses

Our understanding of the diversity and evolution of vertebrate RNA viruses is largely limited to those found in mammalian and avian hosts and associated with overt disease. Here, using a large-scale meta-transcriptomic approach, we discover 214 vertebrate-associated viruses in reptiles, amphibians, lungfish, ray-finned fish, cartilaginous fish and jawless fish. The newly discovered viruses appear in every family or genus of RNA virus associated with vertebrate infection, including those containing human pathogens such as influenza virus, the Arenaviridae and Filoviridae families, and have branching orders that broadly reflected the phylogenetic history of their hosts. We establish a long evolutionary history for most groups of vertebrate RNA virus, and support this by evaluating evolutionary timescales using dated orthologous endogenous virus elements. We also identify new vertebrate-specific RNA viruses and genome architectures, and re-evaluate the evolution of vector-borne RNA viruses. In summary, this study reveals diverse virus-host associations across the entire evolutionary history of the vertebrates.

Landscape attributes governing local transmission of an endemic zoonosis: Rabies virus in domestic dogs

Landscape heterogeneity plays an important role in disease spread and persistence, but quantifying landscape influences and their scale dependence is challenging. Studies have focused on how environmental features or global transport networks influence pathogen invasion and spread, but their influence on local transmission dynamics that underpin the persistence of endemic diseases remains unexplored. Bayesian phylogeographic frameworks that incorporate spatial heterogeneities are promising tools for analysing linked epidemiological, environmental and genetic data. Here, we extend these methodological approaches to decipher the relative contribution and scale-dependent effects of landscape influences on the transmission of endemic rabies virus in Serengeti district, Tanzania (area ~4,900 km2 ). Utilizing detailed epidemiological data and 152 complete viral genomes collected between 2004 and 2013, we show that the localized presence of dogs but not their density is the most important determinant of diffusion, implying that culling will be ineffective for rabies control. Rivers and roads acted as barriers and facilitators to viral spread, respectively, and vaccination impeded diffusion despite variable annual coverage. Notably, we found that landscape effects were scale-dependent: rivers were barriers and roads facilitators on larger scales, whereas the distribution of dogs was important for rabies dispersal across multiple scales. This nuanced understanding of the spatial processes that underpin rabies transmission can be exploited for targeted control at the scale where it will have the greatest impact. Moreover, this research demonstrates how current phylogeographic frameworks can be adapted to improve our understanding of endemic disease dynamics at different spatial scales.

Barriers to HIV Cure

Since the beginning of the epidemic, more than 70 million people have been infected with human immunodeficiency virus (HIV) and about 38 million have died from acquired immune deficiency syndrome (AIDS)-related illnesses. While the discovery of highly active antiretroviral therapy (HAART) in the mid 90's has saved millions of lives, a complete eradication of HIV is still not possible as HIV can persist for decades in a small reservoir of latently infected cells. Once reactivated, these latently infected cells can actively produce viral particles. Recent studies suggest that several sanctuaries exist within infected individuals where HIV can remain undetected by the immune system. These cellular, anatomical and microanatomical viral reservoirs represent a major obstacle for the eradication of HIV. Here we review recent findings on potential sanctuaries of HIV and address potential avenues to overcome these immunological barriers.

Insights into the Impact of CD8+ Immune Modulation on Human Immunodeficiency Virus Evolutionary Dynamics in Distinct Anatomical Compartments by Using Simian Immunodeficiency Virus-Infected Macaque Models of AIDS Progression

A thorough understanding of the role of human immunodeficiency virus (HIV) intrahost evolution in AIDS pathogenesis has been limited by the need for longitudinally sampled viral sequences from the vast target space within the host, which are often difficult to obtain from human subjects. CD8⁺ lymphocyte-depleted macaques infected with simian immunodeficiency virus (SIV) provide an increasingly utilized model of pathogenesis due to clinical manifestations similar to those for HIV-1 infection and AIDS progression, as well as a characteristic rapid disease onset. Comparison of this model with SIV-infected non-CD8⁺ lymphocyte-depleted macaques also provides a unique opportunity to investigate the role of CD8⁺ cells in viral evolution and population dynamics throughout the duration of infection. Using several different phylogenetic methods, we analyzed viral gp120 sequences obtained from extensive longitudinal sampling of multiple tissues and enriched leukocyte populations from SIVmac251-infected macaques with or without CD8⁺ lymphocyte depletion. SIV evolutionary and selection patterns in non-CD8⁺ lymphocyte-depleted animals were characterized by sequential population turnover and continual viral adaptation, a scenario readily comparable to intrahost evolutionary patterns during human HIV infection in the absence of antiretroviral therapy. Alternatively, animals that were depleted of CD8⁺ lymphocytes exhibited greater variation in population dynamics among tissues and cell populations over the course of infection. Our findings highlight the major role for CD8⁺ lymphocytes in prolonging disease progression through continual control of SIV subpopulations from various anatomical compartments and the potential for greater independent viral evolutionary behavior among these compartments in response to immune modulation.IMPORTANCE Although developments in combined antiretroviral therapy (cART) strategies have successfully prolonged the time to AIDS onset in HIV-1-infected individuals, a functional cure has yet to be found. Improvement of drug interventions for a virus that is able to infect a wide range of tissues and cell types requires a thorough understanding of viral adaptation and infection dynamics within this target milieu. Although it is difficult to accomplish in the human host, longitudinal sampling of multiple anatomical locations is readily accessible in the SIV-infected macaque models of neuro-AIDS. The significance of our research is in identifying the impact of immune modulation, through differing immune selective pressures, on viral evolutionary behavior in a multitude of anatomical compartments. The results provide evidence encouraging the development of a more sophisticated model that considers a network of individual viral subpopulations within the host, with differing infection and transmission dynamics, which is necessary for more effective treatment strategies.

Molecular epidemiology of hepatitis B virus infection in Tanzania

Despite the significant public health problems associated with hepatitis B virus (HBV) in sub-Saharan Africa, many countries in this region do not have systematic HBV surveillance or genetic information on HBV circulating locally. Here, we report on the genetic characterization of 772 HBV strains from Tanzania. Phylogenetic analysis of the S-gene sequences showed prevalence of HBV genotype A (HBV/A, n=671, 86.9 %), followed by genotypes D (HBV/D, n=95, 12.3 %) and E (HBV/E, n=6, 0.8 %). All HBV/A sequences were further classified into subtype A1, while the HBV/D sequences were assigned to a new cluster. Among the Tanzanian sequences, 84 % of HBV/A1 and 94 % of HBV/D were unique. The Tanzanian and global HBV/A1 sequences were compared and were completely intermixed in the phylogenetic tree, with the Tanzanian sequences frequently generating long terminal branches, indicating a long history of HBV/A1 infections in the country. The time to the most recent common ancestor was estimated to be 188 years ago [95 % highest posterior density (HPD): 132 to 265 years] for HBV/A1 and 127 years ago (95 % HPD: 79 to 192 years) for HBV/D. The Bayesian skyline plot showed that the number of transmissions 'exploded' exponentially between 1960-1970 for HBV/A1 and 1970-1990 for HBV/D, with the effective population of HBV/A1 having expanded twice as much as that of HBV/D. The data suggest that Tanzania is at least a part of the geographic origin of the HBV/A1 subtype. A recent increase in the transmission rate and significant HBV genetic diversity should be taken into consideration when devising public health interventions to control HBV infections in Tanzania.

Assessment of contemporary genetic diversity and inter-taxa/inter-region exchange of avian paramyxovirus serotype 1 in wild birds sampled in North America

BACKGROUND

Avian paramyxovirus serotype 1 (APMV-1) viruses are globally distributed, infect wild, peridomestic, and domestic birds, and sometimes lead to outbreaks of disease. Thus, the maintenance, evolution, and spread of APMV-1 viruses are relevant to avian health.

METHODS

In this study we sequenced the fusion gene from 58 APMV-1 isolates recovered from thirteen species of wild birds sampled throughout the USA during 2007-2014. We analyzed sequence information with previously reported data in order to assess contemporary genetic diversity and inter-taxa/inter-region exchange of APMV-1 in wild birds sampled in North America.

RESULTS

Our results suggest that wild birds maintain previously undescribed genetic diversity of APMV-1; however, such diversity is unlikely to be pathogenic to domestic poultry. Phylogenetic analyses revealed that APMV-1 diversity detected in wild birds of North America has been found in birds belonging to numerous taxonomic host orders and within hosts inhabiting multiple geographic regions suggesting some level of viral exchange. However, our results also provide statistical support for associations between phylogenetic tree topology and host taxonomic order/region of sample origin which supports restricted exchange among taxa and geographical regions of North America for some APMV-1 sub-genotypes.

CONCLUSIONS

We identify previously unrecognized genetic diversity of APMV-1 in wild birds in North America which is likely a function of continued viral evolution in reservoir hosts. We did not, however, find support for the emergence or maintenance of APMV-1 strains predicted to be pathogenic to poultry in wild birds of North America outside of the order Suliformes (i.e., cormorants). Furthermore, genetic evidence suggests that ecological drivers or other mechanisms may restrict viral exchange among taxa and regions of North America. Additional and more systematic sampling for APMV-1 in North America would likely provide further inference on viral dynamics for this infectious agent in wild bird populations.

Resolving relationships at the animal-fungal divergence: A molecular phylogenetic study of the protist trichomycetes (Ichthyosporea, Eccrinida)

Trichomycetes is a group of microorganisms that was considered a class of fungi comprising four orders of commensal, gut-dwelling endosymbionts obligately associated with arthropods. Since molecular phylogenies revealed two of those orders (Amoebidiales and Eccrinales="protist trichos") to be closely related to members of the protist class Ichthyosporea (=Mesomycetozoea), trichomycetes have been considered an ecological association of both early-diverging fungi and protists. Understanding of the taxonomy, evolution, and diversity of the protist trichos is lacking largely due to the difficulties inherent in species collection that have contributed to undersampling and understudy. The most recent classification divides the protist trichos between two families, Amoebidiidae and Eccrinidae (suborder Trichomycina, order Eccrinida). However, there is no comprehensive molecular phylogeny available for this group and major questions about the systematics of protist trichos remain unanswered. Therefore, we generated 18S and 28S rDNA sequences for 106 protist tricho samples and combined them with publicly available Eccrinida sequences for phylogenetic analyses. We also sequenced a conserved protein-coding gene (heat-shock 70 protein) to obtain a multigene data set. We conducted ancestral state reconstruction (ASR) and Bayesian tip-association significance test (BaTS) analyses by mapping six morphological and ecological characters onto the resulting phylogenetic trees. Our results demonstrate: (1) several ecological and morphological character states (habitat, host type, host stage at time of infestation, location within host, spore production, and growth form) are significantly correlated with the phylogeny, and (2) two additional protist tricho families should be incorporated into the taxonomy to reflect phylogenetic relationships. Our data suggest that an integrated strategy that combines morphological, ecological, and molecular characters is needed to further resolve and clarify the systematics of the Eccrinida.

A molecular portrait of maternal sepsis from Byzantine Troy

Pregnancy complications are poorly represented in the archeological record, despite their importance in contemporary and ancient societies. While excavating a Byzantine cemetery in Troy, we discovered calcified abscesses among a woman’s remains. Scanning electron microscopy of the tissue revealed ‘ghost cells’, resulting from dystrophic calcification, which preserved ancient maternal, fetal and bacterial DNA of a severe infection, likely chorioamnionitis. Gardnerella vaginalis and Staphylococcus saprophyticus dominated the abscesses. Phylogenomic analyses of ancient, historical, and contemporary data showed that G. vaginalis Troy fell within contemporary genetic diversity, whereas S. saprophyticus Troy belongs to a lineage that does not appear to be commonly associated with human disease today. We speculate that the ecology of S. saprophyticus infection may have differed in the ancient world as a result of close contacts between humans and domesticated animals. These results highlight the complex and dynamic interactions with our microbial milieu that underlie severe maternal infections.

DOI:http://dx.doi.org/10.7554/eLife.20983.001

Why and how have some bacteria evolved to cause illness in humans? One way to study bacterial evolution is to search for ancient samples of bacteria and use DNA sequencing technology to investigate how modern bacteria have changed from their ancestors. Understanding the evolution process may help researchers to understand how some bacteria become resistant to the antibiotics designed to kill them.

Complications that occur during pregnancy, including bacterial infections, have long been a major cause of death for women. Now, Devault, Mortimer et al. have been able to sequence the DNA of bacteria found in tissue collected from a woman buried 800 years ago in a cemetery in Troy. Some of the woman’s tissues had been well preserved because they had calcified (probably as the result of infection), which preserved their structure in a mineralized layer. Two mineralized “nodules” in the body appear to be the remains of abscesses. Some of the human DNA in the nodules came from a male, suggesting that the woman was pregnant with a boy and that the abscesses formed in placental tissue.

Sequencing the DNA of the bacteria in the abscess allowed Devault, Mortimer et al. to diagnose the woman’s infection, which was caused by two types of bacteria. One species, called Gardnerella vaginalis, is found in modern pregnancy-related infections. The DNA of the ancient samples was similar to that of modern bacteria. The other bacteria species was an ancient form of Staphylococcus saprophyticus, a type of bacteria that causes urinary tract infections. However, the DNA of the ancient S. saprophyticus bacteria is quite different to that of the bacteria found in modern humans. Instead, their DNA sequence appears more similar to forms of the bacteria that infect currently livestock. As humans lived closely with their livestock at the time the woman lived, her infection may be due to a type of bacteria that passed easily between humans and animals.

Overall, the results suggest that the disease-causing properties of bacteria can arise from a wide range of sources. In addition, Devault, Mortimer et al. have demonstrated that certain types of tissue found in archeological remains are a potential gold mine of information about the evolution of bacteria and other microbes found in the human body.

DOI:http://dx.doi.org/10.7554/eLife.20983.002

Molecular Evidence of Transmission of Influenza A/H1N1 2009 on a University Campus

Background

In the recent years, the data on the molecular epidemiology of influenza viruses have expanded enormously because of the availability of cutting-edge sequencing technologies. However, much of the information is from the temperate regions with few studies from tropical regions such as South-east Asia. Despite the fact that influenza has been known to transmit rapidly within semi-closed communities, such as military camps and educational institutions, data are limited from these communities.

Objectives

To determine the phylogeography of influenza viruses on a university campus, we examined the spatial distribution of influenza virus on the National University of Singapore (NUS) campus.

Methods

Consenting students from the NUS who sought medical attention at the UHC provided two nasopharyngeal swabs and demographic data. PCR was used for detection of influenza viruses. 34 full-genomes of pH1N1/09 viruses were successfully sequenced by Sanger method and concatenated using Geneious R7. Phylogenetic analysis was conducted using these 34 sequences and 1518 global sequences. Phylogeographic analysis was done using BaTS software and Association index and Fitch parsimony scores were determined.

Results

Integrating whole genome sequencing data with epidemiological data, we found strong evidence of influenza transmission on campus as isolates from students residing on-campus were highly similar to each other (AI, P value = 0.009; PS, P value = 0.04). There was also evidence of multiple introductions from the community.

Conclusions

Such data are useful in formulating pandemic preparedness plans which can use these communities as sentinel sites for detection and monitoring of emerging respiratory viral infections.

Sequential CCR5-Tropic HIV-1 Reactivation from Distinct Cellular Reservoirs following Perturbation of Elite Control

Background

HIV Elite Controllers may reveal insights into virus persistence given they harbour small reservoir sizes, akin to HIV non-controllers treated early with combination antiretroviral therapy. Both groups of patients represent the most promising candidates for interventions aimed at sustained remission or ‘cure’. Analytic treatment interruption (ATI) in the latter group leads to stochastic rebound of virus, though it is unclear whether loss of elite control is also associated with similar rebound characteristics.

Methods

We studied three discrete periods of virus rebound during myeloma related immune disruption over 2.5 years in an elite controller who previously underwent autologous stem cell transplantation (ASCT) in the absence of any antiretroviral therapy. Single genome sequencing of the V1-V4 region of env in PBMC and plasma was performed and phylogenies reconstructed. Average pairwise distance (APD) was calculated and non-parametric methods used to assess compartmentalisation. Coreceptor usage was predicted based on genotypic algorithms.

Results

122 single genome sequences were obtained (median 26 sequences per rebound). The initial rebounding plasma env sequences following ASCT represented two distinct lineages, and clustered with proviral DNA sequences isolated prior to ASCT. One of the lineages was monophyletic, possibly indicating reactivation from clonally expanded cells. The second rebound occurred 470 days after spontaneous control of the first rebound and was phylogenetically distinct from the first, confirmed by compartmentalisation analysis, with a different cellular origin rather than ongoing replication. By contrast, third rebound viruses clustered with second rebound viruses, with evidence for ongoing evolution that was associated with lymphopenia and myeloma progression. Following ASCT a shift in tropism from CXCR4-tropic viruses to a CCR5-tropic population was observed to persist through to the third rebound.

Conclusions

Our data highlight similarities in the viral reservoir between elite and non-controllers undergoing ATI following allogeneic transplantation. The lack of propagation of CXCR4 using viruses following transplantation warrants further study.

Evolution of Neuroadaptation in the Periphery and Purifying Selection in the Brain Contribute to Compartmentalization of Simian Immunodeficiency Virus (SIV) in the Brains of Rhesus Macaques with SIV-Associated Encephalitis

The emergence of a distinct subpopulation of human or simian immunodeficiency virus (HIV/SIV) sequences within the brain (compartmentalization) during infection is hypothesized to be linked to AIDS-related central nervous system (CNS) neuropathology. However, the exact evolutionary mechanism responsible for HIV/SIV brain compartmentalization has not been thoroughly investigated. Using extensive viral sampling from several different peripheral tissues and cell types and from three distinct regions within the brain from two well-characterized rhesus macaque models of the neurological complications of HIV infection (neuroAIDS), we have been able to perform in-depth evolutionary analyses that have been unattainable in HIV-infected subjects. The results indicate that, despite multiple introductions of virus into the brain over the course of infection, brain sequence compartmentalization in macaques with SIV-associated CNS neuropathology likely results from late viral entry of virus that has acquired through evolution in the periphery sufficient adaptation for the distinct microenvironment of the CNS.

IMPORTANCE

HIV-associated neurocognitive disorders remain prevalent among HIV type 1-infected individuals, whereas our understanding of the critical components of disease pathogenesis, such as virus evolution and adaptation, remains limited. Building upon earlier findings of specific viral subpopulations in the brain, we present novel yet fundamental results concerning the evolutionary patterns driving this phenomenon in two well-characterized animal models of neuroAIDS and provide insight into the timing of entry of virus into the brain and selective pressure associated with viral adaptation to this particular microenvironment. Such knowledge is invaluable for therapeutic strategies designed to slow or even prevent neurocognitive impairment associated with AIDS.

Phylodynamics of influenza A(H3N2) in South America, 1999-2012

The limited influenza A(H3N2) genetic data available from the Southern Hemisphere (particularly from Africa and Latin America), constrains the accurate reconstruction of viral dissemination dynamics within those regions. Our objective was to describe the spatial dissemination dynamics of influenza A(H3N2) within South America. A total of 469 sequences of the HA1 portion of the hemagglutinin gene (HA) from influenza A(H3N2) viruses sampled in temperate and tropical South American countries between 1999 and 2012 were combined with available contemporary sequences from Australia, Hong Kong, United Kingdom and the United States. Phylogenetic analyses revealed that influenza A(H3N2) sequences from South America were highly intermixed with sequences from other geographical regions, although a clear geographic virus population structure was detected globally. We identified 14 clades mostly (≥80%) composed of influenza sequences from South American countries. Bayesian phylogeographic analyses of those clades support a significant role of both temperate and tropical regions in the introduction and dissemination of new influenza A(H3N2) strains within South America and identify an intensive bidirectional viral exchange between different geographical areas. These findings indicate that seasonal influenza A(H3N2) epidemics in South America are seeded by both the continuous importation of viral variants from other geographic regions and the short-term persistence of local lineages. This study also supports a complex metapopulation model of influenza A(H3N2) dissemination in South America, with no preferential direction in viral movement between temperate and tropical regions.

Phylogenetics and Phyloanatomy of HIV/SIV Intra-Host Compartments and Reservoirs: The Key Role of the Central Nervous System

BACKGROUND

The ability of the human immunodeficiency virus type 1 (HIV-1) to persist in anatomic compartments and cellular reservoirs is a major obstacle for eradication of replicationcompetent virus in the infected host.

APPROACH

We extensively review recent advancements in phylogenetic and phylogeographic techniques that provide a unique opportunity for studies of intra-host HIV-1 compartmentalization and the detection of potential reservoirs.

CONCLUSION

We show that infected macrophages in the central nervous system (CNS) harbor viral subpopulations that play a key role in the emergence of escape variants and viral rebound following discontinuation of antiretroviral therapy. An HIV cure, therefore, cannot be achieved without the effective targeting of the virus in the CNS, for which in depth knowledge of viral population dynamics contributing to the development and maintenance of latent reservoirs is critical.

Host ecology determines the dispersal patterns of a plant virus

Since its isolation in 1966 in Kenya, rice yellow mottle virus (RYMV) has been reported throughout Africa resulting in one of the economically most important tropical plant emerging diseases. A thorough understanding of RYMV evolution and dispersal is critical to manage viral spread in tropical areas that heavily rely on agriculture for subsistence. Phylogenetic analyses have suggested a relatively recent expansion, perhaps driven by the intensification of agricultural practices, but this has not yet been examined in a coherent statistical framework. To gain insight into the historical spread of RYMV within Africa rice cultivations, we analyse a dataset of 300 coat protein gene sequences, sampled from East to West Africa over a 46-year period, using Bayesian evolutionary inference. Spatiotemporal reconstructions date the origin of RMYV back to 1852 (1791-1903) and confirm Tanzania as the most likely geographic origin. Following a single long-distance transmission event from East to West Africa, separate viral populations have been maintained for about a century. To identify the factors that shaped the RYMV distribution, we apply a generalised linear model (GLM) extension of discrete phylogenetic diffusion and provide strong support for distances measured on a rice connectivity landscape as the major determinant of RYMV spread. Phylogeographic estimates in continuous space further complement this by demonstrating more pronounced expansion dynamics in West Africa that are consistent with agricultural intensification and extensification. Taken together, our principled phylogeographic inference approach shows for the first time that host ecology dynamics have shaped the historical spread of a plant virus.