DISSEQT-DIStribution-based modeling of SEQuence space Time dynamics.

Rapidly evolving microbes are a challenge to model because of the volatile, complex, and dynamic nature of their populations. We developed the DISSEQT pipeline (DIStribution-based SEQuence space Time dynamics) for analyzing, visualizing, and predicting the evolution of heterogeneous biological populations in multidimensional genetic space, suited for population-based modeling of deep sequencing and high-throughput data. The pipeline is openly available on GitHub (https://github.com/rasmushenningsson/DISSEQT.jl, accessed 23 June 2019) and Synapse (https://www.synapse.org/#!Synapse: syn11425758, accessed 23 June 2019), covering the entire workflow from read alignment to visualization of results. Our pipeline is centered around robust dimension and model reduction algorithms for analysis of genotypic data with additional capabilities for including phenotypic features to explore dynamic genotype-phenotype maps. We illustrate its utility and capacity with examples from evolving RNA virus populations, which present one of the highest degrees of genetic heterogeneity within a given population found in nature. Using our pipeline, we empirically reconstruct the evolutionary trajectories of evolving populations in sequence space and genotype-phenotype fitness landscapes. We show that while sequence space is vastly multidimensional, the relevant genetic space of evolving microbial populations is of intrinsically low dimension. In addition, evolutionary trajectories of these populations can be faithfully monitored to identify the key minority genotypes contributing most to evolution. Finally, we show that empirical fitness landscapes, when reconstructed to include minority variants, can predict phenotype from genotype with high accuracy.

Plasma RNA viral load is not associated with intrapatient quasispecies heterogeneity in HIV-1 infection.

The human immunodeficiency virus type 1 (HIV-1) viral set point has been associated with the rate of, disease progression and with the level of HIV-specific immune response. The analysis of the possible association between viral set point and quasispecies heterogeneity has important consequences in the understanding of HIV-1 in vivo evolution. In this study, we analyzed the association between intrapatient viral diversity and RNA viral load in 16 antiretroviral therapy-naive HIV-1-infected patients at a single time point, during the disease free period. Patients were separated into low and high viral load groups according to plasma RNA values. HIV-1 quasispecies complexity was assessed in the C2-V5 env region. The average intrapatient quasispecies heterogeneity in both groups was not significantly different (t-test, P > 0.05). However, while within the low viral load group both synonymous and non-synonymous mutations contribute to the variation observed, in the heterogeneity observed in the high viral load group there was an increase in the contribution of the non-synonymous mutations. Thus, this study show that although intrapatient quasispecies heterogeneity is not associated with viral set point in HIV-1 infection, some differences exist between the two groups in the pattern of mutation accumulation.

Transmission bottlenecks and the evolution of fitness in rapidly evolving RNA viruses.

We explored the evolutionary importance of two factors in the adaptation of RNA viruses to their cellular hosts, size of viral inoculum used to initiate a new infection, and mode of transmission (horizontal versus vertical). Transmission bottlenecks should occur in natural populations of viruses and their profound effects on viral adaptation have been previously documented. However, the role of transmission mode has not received the same attention. Here we used a factorial experimental design to test the combined effects of inoculum (bottleneck) size and mode of transmission in evolution of vesicular stomatitis virus (VSV) in tissue culture, and compared our results to the predictions of a recent theoretical model. Our data were in accord with basic genetic principles concerning the balance between mutation, selection and genetic drift. In particular, attenuation of vertically transmitted viruses was a consequence of the random accumulation of deleterious mutations, whereas horizontally transmitted viruses experiencing similar bottlenecks did not suffer the same fitness losses because effective bottleneck size was actually determined by the number of host individuals. In addition, high levels of viral fitness in horizontally transmitted populations were explained by competition among viral variants.