Virulence and pathogenesis of HIV-1 infection: an evolutionary perspective

Identification of antibody targets associated with lower HIV viral load and viremic control

BACKGROUND

High HIV viral loads (VL) are associated with increased morbidity, mortality, and on-going transmission. HIV controllers maintain low VLs in the absence of antiretroviral therapy (ART). We previously used a massively multiplexed antibody profiling assay (VirScan) to compare antibody profiles in HIV controllers and persons living with HIV (PWH) who were virally suppressed on ART. In this report, we used VirScan to evaluate whether antibody reactivity to specific HIV targets and broad reactivity across the HIV genome was associated with VL and controller status 1-2 years after infection.

METHODS

Samples were obtained from participants who acquired HIV infection in a community-randomized trial in Africa that evaluated an integrated strategy for HIV prevention (HPTN 071 PopART). Controller status was determined using VL and antiretroviral (ARV) drug data obtained at the seroconversion visit and 1 year later. Viremic controllers had VLs <2,000 copies/mL at both visits; non-controllers had VLs >2,000 copies/mL at both visits. Both groups had no ARV drugs detected at either visit. VirScan testing was performed at the second HIV-positive visit (1-2 years after HIV infection).

RESULTS

The study cohort included 13 viremic controllers and 64 non-controllers. We identified ten clusters of homologous peptides that had high levels of antibody reactivity (three in gag, three in env, two in integrase, one in protease, and one in vpu). Reactivity to 43 peptides (eight unique epitopes) in six of these clusters was associated with lower VL; reactivity to six of the eight epitopes was associated with HIV controller status. Higher aggregate antibody reactivity across the eight epitopes (more epitopes targeted, higher mean reactivity across all epitopes) and across the HIV genome was also associated with lower VL and controller status.

CONCLUSIONS

We identified HIV antibody targets associated with lower VL and HIV controller status 1-2 years after infection. Robust aggregate responses to these targets and broad antibody reactivity across the HIV genome were also associated with lower VL and controller status. These findings provide novel insights into the relationship between humoral immunity and viral containment that could help inform the design of antibody-based approaches for reducing HIV VL.

Low prevalence of HIV in the northern Cameroon: contribution of some AIDS restriction genes and potential implications for gene therapy

Background

HIV infection and its progression to AIDS depend on several factors including host genetic factors. The immunological mechanisms of host resistance to HIV infection greatly influence the prevalence of HIV in a given region. Worldwide, Cameroon not exempted, the frequency of AIDS-associated genes varies and may influence this prevalence. The North and Far North Regions of Cameroon have had the lowest HIV prevalence in the country for many years despite risky behaviors associated with their customs and habits. In this work, we seek to explore the contribution of host genes to the HIV low prevalence in these regions.

Methodology

Five genes variants previously described as HIV AIDS related were studied. These genes are: CCR5Δ32, CCR5promoter59029G, CCR2-64I, SDF1-3'A and Trim5α(R136Q). A total of 384 consented participants were included in this study. The HIV serological status was confirmed using national algorithm. Genomic DNA was extracted from the buffy coats and used for genotyping. The results obtained were compiled in Excel 2016, Epi Info 7.1 and snpStats software and Chi two tests allowed us to compare the frequencies of the AIDS related alleles in the North with those in other Regions of Cameroon and to measure the impact of these ARGs on protection against HIV.

Results

The frequency of protective alleles CCR5Δ32, CCR5promoter59029G, CCR2-64I, SDF1-3'A and Trim5α(R136Q) was the allelic frequencies should be expressed as percentages i.e. 0.52%; 37.56%; 36.46%; 25.19% and 69.33%. These allelic frequencies exhibited a significant difference when compared to those obtained in other regions of Cameroon (p < 0.01). Protective alleles were predominant in the Northern region compared to others and were associated with resistance to HIV [(p < 0.0001); OR = 2.02 CI, 95%].

Conclusion

The higher frequency of HIV-protective alleles in the northern regions may be a contributing factor to the lower prevalence of HIV. Nevertheless, this should be reinforced by other preventive and surveillance methods to guarantee the sustained low prevalence. HIV can develop resistance through the process of mutation, but the host targets themselves are genetically stable. The study of these host genetic restriction factors is of great value in the design of a practical cure for HIV infection or an effective vaccine.

Antiviral capacity of the early CD8 T-cell response is predictive of natural control of SIV infection: Learning in vivo dynamics using ex vivo data

While most individuals suffer progressive disease following HIV infection, a small fraction spontaneously controls the infection. Although CD8 T-cells have been implicated in this natural control, their mechanistic roles are yet to be established. Here, we combined mathematical modeling and analysis of previously published data from 16 SIV-infected macaques, of which 12 were natural controllers, to elucidate the role of CD8 T-cells in natural control. For each macaque, we considered, in addition to the canonical in vivo plasma viral load and SIV DNA data, longitudinal ex vivo measurements of the virus suppressive capacity of CD8 T-cells. Available mathematical models do not allow analysis of such combined in vivo-ex vivo datasets. We explicitly modeled the ex vivo assay, derived analytical approximations that link the ex vivo measurements with the in vivo effector function of CD8-T cells, and integrated them with an in vivo model of virus dynamics, thus developing a new learning framework that enabled the analysis. Our model fit the data well and estimated the recruitment rate and/or maximal killing rate of CD8 T-cells to be up to 2-fold higher in controllers than non-controllers (p = 0.013). Importantly, the cumulative suppressive capacity of CD8 T-cells over the first 4-6 weeks of infection was associated with virus control (Spearman's ρ = -0.51; p = 0.05). Thus, our analysis identified the early cumulative suppressive capacity of CD8 T-cells as a predictor of natural control. Furthermore, simulating a large virtual population, our model quantified the minimum capacity of this early CD8 T-cell response necessary for long-term control. Our study presents new, quantitative insights into the role of CD8 T-cells in the natural control of HIV infection and has implications for remission strategies.

Phylogenetic signatures reveal multilevel selection and fitness costs in SARS-CoV-2

Background

Large-scale sequencing of SARS-CoV-2 has enabled the study of viral evolution during the COVID-19 pandemic. Some viral mutations may be advantageous to viral replication within hosts but detrimental to transmission, thus carrying a transient fitness advantage. By affecting the number of descendants, persistence times and growth rates of associated clades, these mutations generate localised imbalance in phylogenies. Quantifying these features in closely-related clades with and without recurring mutations can elucidate the tradeoffs between within-host replication and between-host transmission.

Methods

We implemented a novel phylogenetic clustering algorithm ( mlscluster, https://github.com/mrc-ide/mlscluster) to systematically explore time-scaled phylogenies for mutations under transient/multilevel selection. We applied this method to a SARS-CoV-2 time-calibrated phylogeny with >1.2 million sequences from England, and characterised these recurrent mutations that may influence transmission fitness across PANGO-lineages and genomic regions using Poisson regressions and summary statistics.

Results

We found no major differences across two epidemic stages (before and after Omicron), PANGO-lineages, and genomic regions. However, spike, nucleocapsid, and ORF3a were proportionally more enriched for transmission fitness polymorphisms (TFP)-homoplasies than other proteins. We provide a catalog of SARS-CoV-2 sites under multilevel selection, which can guide experimental investigations within and beyond the spike protein.

Conclusions

This study provides empirical evidence for the existence of important tradeoffs between within-host replication and between-host transmission shaping the fitness landscape of SARS-CoV-2. This method may be used as a fast and scalable means to shortlist large sequence databases for sites under putative multilevel selection which may warrant subsequent confirmatory analyses and experimental confirmation.

Proliferation in malaria parasites: How resource limitation can prevent evolution of greater virulence

For parasites, robust proliferation within hosts is crucial for establishing the infection and creating opportunities for onward transmission. While faster proliferation enhances transmission rates, it is often assumed to curtail transmission duration by killing the host (virulence), a trade-off constraining parasite evolution. Yet in many diseases, including malaria, the preponderance of infections with mild or absent symptoms suggests that host mortality is not a sufficient constraint, raising the question of what restrains evolution toward faster proliferation. In malaria infections, the maximum rate of proliferation is determined by the burst size, the number of daughter parasites produced per infected red blood cell. Larger burst sizes should expand the pool of infected red blood cells that can be used to produce the specialized transmission forms needed to infect mosquitoes. We use a within-host model parameterized for rodent malaria parasites (Plasmodium chabaudi) to project the transmission consequences of burst size, focusing on initial acute infection where resource limitation and risk of host mortality are greatest. We find that resource limitation restricts evolution toward higher burst sizes below the level predicted by host mortality alone. Our results suggest resource limitation could represent a more general constraint than virulence-transmission trade-offs, preventing evolution towards faster proliferation.

The association between single-nucleotide polymorphisms within type 1 interferon pathway genes and human immunodeficiency virus type 1 viral load in antiretroviral-naïve participants

BACKGROUND

Human genetic contribution to HIV progression remains inadequately explained. The type 1 interferon (IFN) pathway is important for host control of HIV and variation in type 1 IFN genes may contribute to disease progression. This study assessed the impact of variations at the gene and pathway level of type 1 IFN on HIV-1 viral load (VL).

METHODS

Two cohorts of antiretroviral (ART) naïve participants living with HIV (PLWH) with either early (START) or advanced infection (FIRST) were analysed separately. Type 1 IFN genes (n = 17) and receptor subunits (IFNAR1, IFNAR2) were examined for both cumulated type 1 IFN pathway analysis and individual gene analysis. SKAT-O was applied to detect associations between the genotype and HIV-1 study entry viral load (log10 transformed) as a proxy for set point VL; P-values were corrected using Bonferroni (P < 0.0025).

RESULTS

The analyses among those with early infection included 2429 individuals from five continents. The median study entry HIV VL was 14,623 (IQR 3460-45100) copies/mL. Across 673 SNPs within 19 type 1 IFN genes, no significant association with study entry VL was detected. Conversely, examining individual genes in START showed a borderline significant association between IFNW1, and study entry VL (P = 0.0025). This significance remained after separate adjustments for age, CD4+ T-cell count, CD4+/CD8+ T-cell ratio and recent infection. When controlling for population structure using linear mixed effects models (LME), in addition to principal components used in the main model, this was no longer significant (p = 0.0244). In subgroup analyses stratified by geographical region, the association between IFNW1 and study entry VL was only observed among African participants, although, the association was not significant when controlling for population structure using LME. Of the 17 SNPs within the IFNW1 region, only rs79876898 (A > G) was associated with study entry VL (p = 0.0020, beta = 0.32; G associated with higher study entry VL than A) in single SNP association analyses. The findings were not reproduced in FIRST participants.

CONCLUSION

Across 19 type 1 IFN genes, only IFNW1 was associated with HIV-1 study entry VL in a cohort of ART-naïve individuals in early stages of their infection, however, this was no longer significant in sensitivity analyses that controlled for population structures using LME.

The germ theory revisited: A noncentric view on infection outcome

The germ theory states that pathogenic microorganisms are responsible for causing infectious diseases. The theory is inherently microbe-centric and does not account for variability in disease severity among individuals and asymptomatic carriership-two phenomena indicating an important role for host variability in infection outcome. The basic tenet of the germ theory was recently challenged, and a radically host-centric paradigm referred to as the "full-blown host theory" was proposed. According to this view, the pathogen is reduced to a passive environmental trigger, and the development of disease is instead due to pre-existing immunodeficiencies of the host. Here, we consider the factors that determine disease severity using established knowledge concerning evolutionary biology, microbial pathogenesis, and host-pathogen interactions. We note that the available data support a noncentric view that recognizes key roles for both the causative microbe and the host in dictating infection outcome.

How robust are estimates of key parameters in standard viral dynamic models?

Mathematical models of viral infection have been developed, fitted to data, and provide insight into disease pathogenesis for multiple agents that cause chronic infection, including HIV, hepatitis C, and B virus. However, for agents that cause acute infections or during the acute stage of agents that cause chronic infections, viral load data are often collected after symptoms develop, usually around or after the peak viral load. Consequently, we frequently lack data in the initial phase of viral growth, i.e., when pre-symptomatic transmission events occur. Missing data may make estimating the time of infection, the infectious period, and parameters in viral dynamic models, such as the cell infection rate, difficult. However, having extra information, such as the average time to peak viral load, may improve the robustness of the estimation. Here, we evaluated the robustness of estimates of key model parameters when viral load data prior to the viral load peak is missing, when we know the values of some parameters and/or the time from infection to peak viral load. Although estimates of the time of infection are sensitive to the quality and amount of available data, particularly pre-peak, other parameters important in understanding disease pathogenesis, such as the loss rate of infected cells, are less sensitive. Viral infectivity and the viral production rate are key parameters affecting the robustness of data fits. Fixing their values to literature values can help estimate the remaining model parameters when pre-peak data is missing or limited. We find a lack of data in the pre-peak growth phase underestimates the time to peak viral load by several days, leading to a shorter predicted growth phase. On the other hand, knowing the time of infection (e.g., from epidemiological data) and fixing it results in good estimates of dynamical parameters even in the absence of early data. While we provide ways to approximate model parameters in the absence of early viral load data, our results also suggest that these data, when available, are needed to estimate model parameters more precisely.

HIV-1 Gag, Pol, and Env diversified with limited adaptation since the 1980s

Knowledge of HIV-1 global sequence diversity is critical for developing an effective prophylactic against HIV-1 infection. We developed the Hervé platform to analyze and visualize trends in HIV-1 diversification. Using Hervé, we analyzed 4,830 Env, 4,407 Gag, and 3,002 Pol publicly available independent sequences corresponding to subtypes A1, A6, B, C, D, F1, and G and circulating recombinant forms (CRFs) 01_AE, 02_AG, and 07_BC; sequences were sampled between 1980 and 2020 from 82 countries. HIV-1 diversified with a median of 1.82 amino acid substitutions per year in Env, 0.297 in Gag, and 0.779 in Pol. Yet, Env subtype B diversification plateaued post-2000. Pairwise diversity within subtypes and CRFs increased by 41.82% (range = 24.85%-54.41%) in Env, 56.93% (15.38%-89.16%) in Gag, and 46.12% (11.70%-70.57%) in Pol. Consensus sequences based on sequences sampled in each decade remained relatively stable over time. Similarly, at antibody epitope sites, only 0-8 residues that were minority variants became consensus over time in any subtype/CRF and only one known drug resistance mutation site differed from the reference (subtype G). The apparent contradiction between the fast diversification of HIV-1 and its limited adaptation illustrates that HIV-1 evolution is not directional and its consensus is at the intersection of millions of within-host selective processes occurring in a star-like manner. While a consensus sequence is a better representation of HIV-1 diversity than any individual sequence, consensus sequences have progressively become more distant from the circulating sequences they represent.

IMPORTANCE

Global surveillance of HIV-1 sequences is critical for designing relevant prophylactic and therapeutic interventions to infection. We designed an open-source platform, Hervé, for analyzing and visualizing the diversification dynamics of HIV-1 protein sequences. We characterized the evolution of over 12,000 HIV-1 Env, Gag, and Pol protein sequences from 1980-2020 and found that, despite a steady increase in intra-subtype and circulating recombinant form diversity, the most frequent residue at each site, i.e., the consensus, has varied only moderately.

Discovery, Crystallographic Studies, and Mechanistic Investigations of Novel Phenylalanine Derivatives Bearing a Quinazolin-4-one Scaffold as Potent HIV Capsid Modulators

Optimization of compound 11L led to the identification of novel HIV capsid modulators, quinazolin-4-one-bearing phenylalanine derivatives, displaying potent antiviral activities against both HIV-1 and HIV-2. Notably, derivatives 12a2 and 21a2 showed significant improvements, with 2.5-fold over 11L and 7.3-fold over PF74 for HIV-1, and approximately 40-fold over PF74 for HIV-2. The X-ray co-crystal structures confirmed the multiple pocket occupation of 12a2 and 21a2 in the binding site. Mechanistic studies revealed a dual-stage inhibition profile, where the compounds disrupted capsid-host factor interactions at the early stage and promoted capsid misassembly at the late stage. Remarkably, 12a2 and 21a2 significantly promoted capsid misassembly, outperforming 11L, PF74, and LEN. The substitution of easily metabolized amide bond with quinolin-4-one marginally enhanced the stability of 12a2 in human liver microsomes compared to controls. Overall, 12a2 and 21a2 highlight their potential as potent HIV capsid modulators, paving the way for future advancements in anti-HIV drug design.

Genetic- and Fiber-Diet-Mediated Changes in Antibiotic Resistance Genes in Pig Colon Contents and Feces and Their Driving Factors

Comprehensive studies on the effects of genetics and fiber diets on antibiotic resistance genes (ARGs) remain scarce. In this study, we analyzed the profiles of ARGs in colonic contents and fecal samples of Taoyuan, Duroc, and Xiangcun pigs (n = 10) fed at different fiber levels. Through macrogenomic analysis, we identified a total of 850 unique types of ARGs and classified them into 111 drug resistance classes. The abundance of partially drug-resistant ARGs was higher in the colonic contents of local pig breeds under a large-scale farming model. ARGs were found to be widely distributed among a variety of bacteria, predominantly in the phyla Firmicutes, Proteobacteria, and Bacteroidetes. Fiber diets reduce the abundance of ARGs in colonic contents and feces, and mobile genetic elements (MGEs) and short-chain fatty acids (SCFAs) are important drivers in mediating the effect of fiber diets on the abundance of ARGs. In vitro fermentation experiments confirmed that butyric acid significantly reduced the abundance of ARGs. In summary, the results of this study enhanced our understanding of the distribution and composition of ARGs in the colon of different breeds of pigs and revealed that a fiber diet can reduce ARGs in feces through its Butyric acid, providing reference data for environmental safety.

Comprehensive profiling of pre-infection antibodies identifies HIV targets associated with viremic control and viral load

Background

High HIV viral load (VL) is associated with increased transmission risk and faster disease progression. HIV controllers achieve viral suppression without antiretroviral (ARV) treatment. We evaluated viremic control in a community-randomized trial with >48,000 participants.

Methods

A massively multiplexed antibody profiling system, VirScan, was used to quantify pre- and post-infection antibody reactivity to HIV peptides in 664 samples from 429 participants (13 controllers, 135 viremic non-controllers, 64 other non-controllers, 217 uninfected persons). Controllers had VLs <2,000 copies/mL with no ARV drugs detected at the first HIV-positive visit and one year later. Viremic non-controllers had VLs 2,000 copies/mL with no ARV drugs detected at the first HIV-positive visit. Other non-controllers had either ARV drugs detected at the first HIV-positive visit (n=47) or VLs <2,000 copies/mL with no ARV drugs detected at only one HIV-positive visit (n=17).

Results

We identified pre-infection HIV antibody reactivities that correlated with post-infection VL. Pre-infection reactivity to an epitope in the HR2 domain of gp41 was associated with controller status and lower VL. Pre-infection reactivity to an epitope in the C2 domain of gp120 was associated with non-controller status and higher VL. Different patterns of antibody reactivity were observed over time for these two epitopes.

Conclusion

These studies suggest that pre-infection HIV antibodies are associated with controller status and modulation of HIV VL. These findings may inform research on antibody-based interventions for HIV treatment.

Associations of functional human leucocyte antigen class I groups with HIV viral load in a heterogeneous cohort

OBJECTIVE

Human leucocyte antigen (HLA) class I alleles are the main host genetic factors involved in controlling HIV-1 viral load (VL). Nevertheless, HLA diversity has proven a significant challenge in association studies. We assessed how accounting for binding affinities of HLA class I alleles to HIV-1 peptides facilitate association testing of HLA with HIV-1 VL in a heterogeneous cohort.

DESIGN

Cohort from the Strategic Timing of AntiRetroviral Treatment (START) study.

METHODS

We imputed HLA class I alleles from host genetic data (2546 HIV+ participants) and sampled immunopeptidomes from 2079 host-paired viral genomes (targeted amplicon sequencing). We predicted HLA class I binding affinities to HIV-1 and unspecific peptides, grouping alleles into functional clusters through consensus clustering. These functional HLA class I clusters were used to test associations with HIV VL.

RESULTS

We identified four clades totaling 30 HLA alleles accounting for 11.4% variability in VL. We highlight HLA-B∗57:01 and B∗57:03 as functionally similar but yet overrepresented in distinct ethnic groups, showing when combined a protective association with HIV+ VL (log, β -0.25; adj. P-value < 0.05). We further demonstrate only a slight power reduction when using unspecific immunopeptidomes, facilitating the use of the inferred functional HLA groups in other studies.

CONCLUSION

The outlined computational approach provides a robust and efficient way to incorporate HLA function and peptide diversity, aiding clinical association studies in heterogeneous cohorts. To facilitate access to the proposed methods and results we provide an interactive application for exploring data.

Extraordinary parasite multiplication rates in human malaria infections

For pathogenic organisms, faster rates of multiplication promote transmission success, the potential to harm hosts, and the evolution of drug resistance. Parasite multiplication rates (PMRs) are often quantified in malaria infections, given the relative ease of sampling. Using modern and historical human infection data, we show that established methods return extraordinarily - and implausibly - large PMRs. We illustrate how inflated PMRs arise from two facets of malaria biology that are far from unique: (i) some developmental ages are easier to sample than others; (ii) the distribution of developmental ages changes over the course of infection. The difficulty of accurately quantifying PMRs demonstrates a need for robust methods and a subsequent re-evaluation of what is known even in the well-studied system of malaria.

Death is overrated: the potential role of detection in driving virulence evolution

A common assumption in the evolution of virulence theory literature is that pathogens transmit better when they exploit their host more heavily, but by doing so, they impose a greater risk of killing their host, thus truncating infectious periods and reducing their own opportunities for transmission. Here, I derive an equation for the magnitude of this cost in terms of the infection fatality rate, and in doing so, I show that there are many cases where mortality costs are too small to plausibly constrain increases in host exploitation by pathogens. I propose that pathogen evolution may often be constrained by detection costs, whereby hosts alter their behaviour when infection is detectable, and thus reduce pathogen opportunities for onward transmission. I then derive an inequality to illustrate when mortality costs or detection costs impose stronger constraints on pathogen evolution, and I use empirical data from the literature to demonstrate that detection costs are frequently large in both human and animal populations. Finally, I give examples of how evolutionary predictions can change depending on whether costs of host exploitation are borne out through mortality or detection.

Host-Parasite Coevolution in Primates

Organisms adapt to their environment through evolutionary processes. Environments consist of abiotic factors, but also of other organisms. In many cases, two or more species interact over generations and adapt in a reciprocal way to evolutionary changes in the respective other species. Such coevolutionary processes are found in mutualistic and antagonistic systems, such as predator-prey and host-parasite (including pathogens) relationships. Coevolution often results in an "arms race" between pathogens and hosts and can significantly affect the virulence of pathogens and thus the severity of infectious diseases, a process that we are currently witnessing with SARS-CoV-2. Furthermore, it can lead to co-speciation, resulting in congruent phylogenies of, e.g., the host and parasite. Monkeys and other primates are no exception. They are hosts to a large number of pathogens that have shaped not only the primate immune system but also various ecological and behavioral adaptions. These pathogens can cause severe diseases and most likely also infect multiple primate species, including humans. Here, we briefly review general aspects of the coevolutionary process in its strict sense and highlight the value of cophylogenetic analyses as an indicator for coevolution.

Why does age at HIV infection correlate with set point viral load? An evolutionary hypothesis

BACKGROUND

Set-point viral load (SPVL) correlates with the age at which people acquire HIV. Although immunosenescence may seem like a parsimonious explanation for this, it does not easily explain the observation that the relationship between age and SPVL attenuates when accounting for source partner SPVL. Here we propose an alternative explanation that encompasses this latter finding: that decreasing risk of acquisition with older age generates a selection bottleneck that selects for more virulent strains with age.

METHODS

We adapted a previously published model of HIV transmission and evolution (EvoNetHIV), parameterized here for men who have sex with men (MSM). We conducted a series of simulation experiments that vary seven behavioral or clinical parameters that affect exposure risk as people age. We conducted regressions to determine the mean increase in SPVL per 10-year increase in seroconversion age, with and without source SPVL in the model.

RESULTS

All runs generated significant relationships between seroconversion age and SPVL when not including source SPVL. All saw attenuated relationships, most to near 0, with source SPVL included. Four of our behavioral measures (relational duration, age-related homophily, coital frequency, and mean age at relationship formation) had clear effects on this relationship, all in the hypothesized direction. Combining multiple forms of behavioral heterogeneity yielded an increase of 0.056 log10 copies/mL SPVL per 10-year increase in seroconversion age, nearly as large as that seen in two empirical studies of age-SPVL correlations in MSM.

CONCLUSION

The higher virulence of HIV among those infected later in life may be partly explained by a combination of selective bottlenecks and behavioral heterogeneity by age. Variation in the strength of this effect across populations may be in part due to different behavioral, epidemiological and clinical conditions, and not require assumptions about differences in patterns of immunosenescence among populations.

CRISPR/Cas9: a tool to eradicate HIV-1

The development of antiretroviral therapy (ART) has been effective in suppressing HIV replication. However, severe drug toxicities due to the therapy and its failure in targeting the integrated proviral genome have led to the introduction of a new paradigm of gene-based therapies. With its effective inhibition and high precision, clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein-9 nuclease (Cas9) or CRISPR/Cas9 has emerged as an effective genome editing tool in the last decade. Mediated by guide RNAs (gRNAs), Cas9 endonuclease acts like genetic scissors that can modify specific target sites. With this concept, CRISPR/Cas9 has been used to target the integrated proviral HIV-1 genome both in in vitro as well as in vivo studies including non-human primates. The CRISPR has also been tested for targeting latent HIV-1 by modulating the proviral transcription with the help of a specialized Cas9 mutant. Overcoming the limitations of the current therapy, CRISPR has the potential to become the primary genome editing tool for eradicating HIV-1 infection. In this review, we summarize the recent advancements of CRISPR to target the proviral HIV-1 genome, the challenges and future prospects.

HTLV-1 infection of donor-derived T cells might promote acute graft-versus-host disease following liver transplantation

Acute graft versus host disease (aGVHD) is a rare, but severe complication of liver transplantation (LT). It is caused by the activation of donor immune cells in the graft against the host shortly after transplantation, but the contributing pathogenic factors remain unclear. Here we show that human T cell lymphotropic virus type I (HTLV-1) infection of donor T cells is highly associated with aGVHD following LT. The presence of HTLV-1 in peripheral blood and tissue samples from a discovery cohort of 7 aGVHD patients and 17 control patients is assessed with hybridization probes (TargetSeq), mass cytometry (CyTOF), and multiplex immunohistology (IMC). All 7 of our aGVHD patients display detectable HTLV-1 Tax signals by IMC. We identify donor-derived cells based on a Y chromosome-specific genetic marker, EIF1AY. Thus, we confirm the presence of CD4⁺Tax⁺EIF1AY⁺ T cells and Tax⁺CD68⁺EIF1AY⁺ antigen-presenting cells, indicating HTLV-1 infection of donor immune cells. In an independent cohort of 400 patients, we verify that HTLV-1 prevalence correlates with aGVHD incidence, while none of the control viruses shows significant associations. Our findings thus provide new insights into the aetio-pathology of liver-transplantation-associated aGVHD and raise the possibility of preventing aGVHD prior to transplantation.

Tuning Rex rules HTLV-1 pathogenesis

HTLV-1 is an oncovirus causing ATL and other inflammatory diseases such as HAM/TSP and HU in about 5% of infected individuals. It is also known that HTLV-1-infected cells maintain a disease-free, immortalized, latent state throughout the lifetimes of about 95% of infected individuals. We believe that the stable maintenance of disease-free infected cells in the carrier is an intrinsic characteristic of HTLV-1 that has been acquired during its evolution in the human life cycle. We speculate that the pathogenesis of the virus is ruled by the orchestrated functions of viral proteins. In particular, the regulation of Rex, the conductor of viral replication rate, is expected to be closely related to the viral program in the early active viral replication followed by the stable latency in HTLV-1 infected T cells. HTLV-1 and HIV-1 belong to the family Retroviridae and share the same tropism, e.g., human CD4+ T cells. These viruses show significant similarities in the viral genomic structure and the molecular mechanism of the replication cycle. However, HTLV-1 and HIV-1 infected T cells show different phenotypes, especially in the level of virion production. We speculate that how the activity of HTLV-1 Rex and its counterpart HIV-1 Rev are regulated may be closely related to the properties of respective infected T cells. In this review, we compare various pathological aspects of HTLV-1 and HIV-1. In particular, we investigated the presence or absence of a virally encoded "regulatory valve" for HTLV-1 Rex or HIV-1 Rev to explore its importance in the regulation of viral particle production in infected T cells. Finally, wereaffirm Rex as the key conductor for viral replication and viral pathogenesis based on our recent study on the novel functional aspects of Rex. Since the activity of Rex is closely related to the viral replication rate, we hypothesize that the "regulatory valve" on the Rex activity may have been selectively evolved to achieve the "scenario" with early viral particle production and the subsequent long, stable deep latency in HTLV-1 infected cells.

Contribution of the HIV-1 Envelope Glycoprotein to AIDS Pathogenesis and Clinical Progression

In the absence of antiviral therapy, HIV-1 infection progresses to a wide spectrum of clinical manifestations that are the result of an entangled contribution of host, immune and viral factors. The contribution of these factors is not completely established. Several investigations have described the involvement of the immune system in the viral control. In addition, distinct HLA-B alleles, HLA-B27, -B57-58, were associated with infection control. The combination of these elements and antiviral host restriction factors results in different clinical outcomes. The role of the viral proteins in HIV-1 infection has been, however, less investigated. We will review contributions dedicated to the pathogenesis of HIV-1 infection focusing on studies identifying the function of the viral envelope glycoprotein (Env) in the clinical progression because of its essential role in the initial events of the virus life-cycle. Some analysis showed that inefficient viral Envs were dominant in non-progressor individuals. These poorly-functional viral proteins resulted in lower cellular activation, viral replication and minor viral loads. This limited viral antigenic production allows a better immune response and a lower immune exhaustion. Thus, the properties of HIV-1 Env are significant in the clinical outcome of the HIV-1 infection and AIDS pathogenesis.

A phylogeny-aware GWAS framework to correct for heritable pathogen effects on infectious disease traits

Infectious diseases are particularly challenging for genome-wide association studies (GWAS) because genetic effects from two organisms (pathogen and host) can influence a trait. Traditional GWAS assume individual samples are independent observations. However, pathogen effects on a trait can be heritable from donor to recipient in transmission chains. Thus, residuals in GWAS association tests for host genetic effects may not be independent due to shared pathogen ancestry. We propose a new method to estimate and remove heritable pathogen effects on a trait based on the pathogen phylogeny prior to host GWAS, thus restoring independence of samples. In simulations, we show this additional step can increase GWAS power to detect truly associated host variants when pathogen effects are highly heritable, with strong phylogenetic correlations. We applied our framework to data from two different host-pathogen systems, HIV in humans and X. arboricola in A. thaliana. In both systems, the heritability and thus phylogenetic correlations turn out to be low enough such that qualitative results of GWAS do not change when accounting for the pathogen shared ancestry through a correction step. This means that previous GWAS results applied to these two systems should not be biased due to shared pathogen ancestry. In summary, our framework provides additional information on the evolutionary dynamics of traits in pathogen populations and may improve GWAS if pathogen effects are highly phylogenetically correlated amongst individuals in a cohort.

Parasite-induced shifts in host movement may explain the transient coexistence of high- and low-pathogenic disease strains

Many parasites induce decreased host movement, known as lethargy, which can impact disease spread and the evolution of virulence. Mathematical models have investigated virulence evolution when parasites cause host death, but disease-induced decreased host movement has received relatively less attention. Here, we consider a model where, due to the within-host parasite replication rate, an infected host can become lethargic and shift from a moving to a resting state, where it can die. We find that when the lethargy and disease-induced mortality costs to the parasites are not high, then evolutionary bistability can arise, and either moderate or high virulence can evolve depending on the initial virulence and the magnitude of mutation. These results suggest, firstly, the coexistence of strains with different virulence, which may explain the transient coexistence of low- and high-pathogenic strains of avian influenza viruses, and secondly, that medical interventions to treat the symptoms of lethargy or prevent disease-induced host deaths can result in a large jump in virulence and the rapid evolution of high virulence. In complement to existing results that show bistability when hosts are heterogeneous at the population level, we show that evolutionary bistability may arise due to transmission heterogeneity at the individual host level.

Functional and Highly Cross-Linkable HIV-1 Envelope Glycoproteins Enriched in a Pretriggered Conformation

Binding to the receptor, CD4, drives the pretriggered, "closed" (state-1) conformation of the human immunodeficiency virus type 1 (HIV-1) envelope glycoprotein (Env) trimer into more "open" conformations (states 2 and 3). Broadly neutralizing antibodies, which are elicited inefficiently, mostly recognize the state-1 Env conformation, whereas the more commonly elicited poorly neutralizing antibodies recognize states 2/3. HIV-1 Env metastability has created challenges for defining the state-1 structure and developing immunogens mimicking this labile conformation. The availability of functional state-1 Envs that can be efficiently cross-linked at lysine and/or acidic amino acid residues might assist these endeavors. To that end, we modified HIV-1_AD8 Env, which exhibits an intermediate level of triggerability by CD4. We introduced lysine/acidic residues at positions that exhibit such polymorphisms in natural HIV-1 strains. Env changes that were tolerated with respect to gp120-gp41 processing, subunit association, and virus entry were further combined. Two common polymorphisms, Q114E and Q567K, as well as a known variant, A582T, additively rendered pseudoviruses resistant to cold, soluble CD4, and a CD4-mimetic compound, phenotypes indicative of stabilization of the pretriggered state-1 Env conformation. Combining these changes resulted in two lysine-rich HIV-1_AD8 Env variants (E.2 and AE.2) with neutralization- and cold-resistant phenotypes comparable to those of natural, less triggerable tier 2/3 HIV-1 isolates. Compared with these and the parental Envs, the E.2 and AE.2 Envs were cleaved more efficiently and exhibited stronger gp120-trimer association in detergent lysates. These highly cross-linkable Envs enriched in a pretriggered conformation should assist characterization of the structure and immunogenicity of this labile state. IMPORTANCE The development of an efficient vaccine is critical for combating HIV-1 infection worldwide. However, the instability of the pretriggered shape (state 1) of the viral envelope glycoprotein (Env) makes it difficult to raise neutralizing antibodies against HIV-1. Here, by introducing multiple changes in Env, we derived two HIV-1 Env variants that are enriched in state 1 and can be efficiently cross-linked to maintain this shape. These Env complexes are more stable in detergent, assisting their purification. Thus, our study provides a path to a better characterization of the native pretriggered Env, which should assist vaccine development.

The Characteristics of the HIV-1 Env Glycoprotein Are Linked With Viral Pathogenesis

The understanding of HIV-1 pathogenesis and clinical progression is incomplete due to the variable contribution of host, immune, and viral factors. The involvement of viral factors has been investigated in extreme clinical phenotypes from rapid progressors to long-term non-progressors (LTNPs). Among HIV-1 proteins, the envelope glycoprotein complex (Env) has been concentrated on in many studies for its important role in the immune response and in the first steps of viral replication. In this study, we analyzed the contribution of 41 Envs from 24 patients with different clinical progression rates and viral loads (VLs), LTNP-Elite Controllers (LTNP-ECs); Viremic LTNPs (vLTNPs), and non-controller individuals contemporary to LTNPs or recent, named Old and Modern progressors. We studied the Env expression, the fusion and cell-to-cell transfer capacities, as well as viral infectivity. The sequence and phylogenetic analysis of Envs were also performed. In every functional characteristic, the Envs from subjects with viral control (LTNP-ECs and vLTNPs) showed significant lower performance compared to those from the progressor individuals (Old and Modern). Regarding sequence analysis, the variable loops of the gp120 subunit of the Env (i.e., V2, V4, and mainly V5) of the progressor individuals showed longer and more glycosylated sequences than controller subjects. Therefore, HIV-1 Envs from virus of patients presenting viremic control and the non-progressor clinical phenotype showed poor viral functions and shorter sequences, whereas functional Envs were associated with virus of patients lacking virological control and with progressor clinical phenotypes. These correlations support the role of Env genotypic and phenotypic characteristics in the in vivo HIV-1 infection and pathogenesis.

Phylogenetic estimation of the viral fitness landscape of HIV-1 set-point viral load

Set-point viral load (SPVL), a common measure of human immunodeficiency virus (HIV)-1 virulence, is partially determined by viral genotype. Epidemiological evidence suggests that this viral property has been under stabilising selection, with a typical optimum for the virus between 104 and 105 copies of viral RNA per ml. Here we aimed to detect transmission fitness differences between viruses from individuals with different SPVLs directly from phylogenetic trees inferred from whole-genome sequences. We used the local branching index (LBI) as a proxy for transmission fitness. We found that LBI is more sensitive to differences in infectiousness than to differences in the duration of the infectious state. By analysing subtype-B samples from the Bridging the Evolution and Epidemiology of HIV in Europe project, we inferred a significant positive relationship between SPVL and LBI up to approximately 105 copies/ml, with some evidence for a peak around this value of SPVL. This is evidence of selection against low values of SPVL in HIV-1 subtype-B strains, likely related to lower infectiousness, and perhaps a peak in the transmission fitness in the expected range of SPVL. The less prominent signatures of selection against higher SPVL could be explained by an inherent limit of the method or the deployment of antiretroviral therapy.

Epstein-Barr Virus in Burkitt Lymphoma in Africa Reveals a Limited Set of Whole Genome and LMP-1 Sequence Patterns: Analysis of Archival Datasets and Field Samples From Uganda, Tanzania, and Kenya

Epstein-Barr virus (EBV) is associated with endemic Burkitt lymphoma (eBL), but the contribution of EBV variants is ill-defined. Studies of EBV whole genome sequences (WGS) have identified phylogroups that appear to be distinct for Asian versus non-Asian EBV, but samples from BL or Africa, where EBV was first discovered, are under-represented. We conducted a phylogenetic analysis of EBV WGS and LMP-1 sequences obtained primarily from BL patients in Africa and representative non-African EBV from other conditions or regions using data from GenBank, Sequence Read Archive, or Genomic Data Commons for the Burkitt Lymphoma Genome Sequencing Project (BLGSP) to generate data to support the use of a simpler biomarker of geographic or phenotypic associations. We also investigated LMP-1 patterns in 414 eBL cases and 414 geographically matched controls in the Epidemiology of Burkitt Lymphoma in East African children and minors (EMBLEM) study using LMP-1 PCR and Sanger sequencing. Phylogenetic analysis revealed distinct genetic patterns of African versus Asian EBV sequences. We identified 281 single nucleotide variations (SNVs) in LMP-1 promoter and coding region, which formed 12 unique patterns (A to L). Nine patterns (A, AB, C, D, F, I, J, K and L) predominated in African EBV, of which four were found in 92% of BL samples (A, AB, D, and H). Predominant patterns were B and G in Asia and H in Europe. EBV positivity in peripheral blood was detected in 95.6% of EMBLEM eBL cases versus 79.2% of the healthy controls (odds ratio [OR] =3.83; 95% confidence interval 2.06-7.14). LMP-1 was successfully sequenced in 66.7% of the EBV DNA positive cases but in 29.6% of the controls (ORs ranging 5-11 for different patterns). Four LMP-1 patterns (A, AB, D, and K) were detected in 63.1% of the cases versus 27.1% controls (ORs ranges: 5.58-11.4). Dual strain EBV infections were identified in WGS and PCR-Sanger data. In conclusion, EBV from Africa is phylogenetically separate from EBV in Asia. Genetic diversity in LMP-1 formed 12 patterns, which showed promising geographic and phenotypic associations. Presence of multiple strain infection should be considered in efforts to refine or improve EBV markers of ancestry or phenotype.

Lay Summary

Epstein-Barr virus (EBV) infection, a ubiquitous infection, contributes to the etiology of both Burkitt Lymphoma (BL) and nasopharyngeal carcinoma, yet their global distributions vary geographically with no overlap. Genomic variation in EBV is suspected to play a role in the geographical patterns of these EBV-associated cancers, but relatively few EBV samples from BL have been comprehensively studied. We sought to compare phylogenetic patterns of EBV genomes obtained from BL samples in Africa and from tumor and non-tumor samples from elsewhere. We concluded that EBV obtained from BL in Africa is genetically separate from EBV in Asia. Through comprehensive analysis of nucleotide variations in EBV's LMP-1 gene, we describe 12 LMP-1 patterns, two of which (B and G) were found mostly in Asia. Four LMP-1 patterns (A, AB, D, and F) accounted for 92% of EBVs sequenced from BL in Africa. Our results identified extensive diversity of EBV, but BL in Africa was associated with a limited number of variants identified, which were different from those identified in Asia. Further research is needed to optimize the use of PCR and sequencing to study LMP-1 diversity for classification of EBV variants and for use in epidemiologic studies to characterize geographic and/or phenotypic associations of EBV variants with EBV-associated malignancies, including eBL.

Transmitted HIV-1 is more virulent in heterosexual individuals than men-who-have-sex-with-men

Transmission bottlenecks introduce selection pressures on HIV-1 that vary with the mode of transmission. Recent studies on small cohorts have suggested that stronger selection pressures lead to fitter transmitted/founder (T/F) strains. Manifestations of this selection bias at the population level have remained elusive. Here, we analysed early CD4 cell count measurements reported from ∼340,000 infected heterosexual individuals (HET) and men-who-have-sex-with-men (MSM), across geographies, ethnicities and calendar years. The reduction in CD4 counts early in infection is reflective of the virulence of T/F strains. MSM and HET use predominant modes of transmission, namely, anal and penile-vaginal, with among the largest differences in the selection pressures at transmission across modes. Further, in most geographies, the groups show little inter-mixing, allowing for the differential selection bias to be sustained and amplified. We found that the early reduction in CD4 counts was consistently greater in HET than MSM (P<0.05). To account for inherent variations in baseline CD4 counts, we constructed a metric to quantify the extent of progression to AIDS as the ratio of the reduction in measured CD4 counts from baseline and the reduction associated with AIDS. We found that this progression corresponding to the early CD4 measurements was ∼68% for MSM and ∼87% for HET on average (P<10⁻⁴; Cohen’s d, d_s = 0.36), reflecting the more severe disease caused by T/F strains in HET than MSM at the population level. Interestingly, the set-point viral load was not different between the groups (d_s<0.12), suggesting that MSM were more tolerant and not more resistant to their T/F strains than HET. This difference remained when we controlled for confounding factors using multivariable regression. We concluded that the different selection pressures at transmission have resulted in more virulent T/F strains in HET than MSM. These findings have implications for our understanding of HIV-1 pathogenesis, evolution, and epidemiology.

HIV-1 encounters a key bottleneck at the time of its transmission from one individual to another. This transmission bottleneck can differ between modes of transmission. The stronger this bottleneck is, the more fit the virus has to be to be successfully transmitted. Accordingly, the transmitted/founder (T/F) strains of HIV-1 may have different fitness in risk groups that use different modes of transmission. While studies on small cohorts do support this notion, observations of the manifestations of this differential selection bias at the population level have been lacking. Here, we examined reported early CD4 count measurements from ∼340,000 HET and MSM, across geographies, ethnicities, and calendar years. Early CD4 counts are a measure of the severity of the infection due to T/F strains. HET and MSM transmit predominantly via penile-vaginal and anal modes, respectively, and do not inter-mix significantly. Remarkably, we found that HET consistently had lower early CD4 counts than MSM. This difference could not be attributed to potential confounding factors, such as set-point viral load. The difference thus provided evidence that T/F strains had evolved to be more virulent in HET than MSM at the population level. Intervention strategies may benefit from accounting for this difference between risk groups.

Evolution during primary HIV infection does not require adaptive immune selection

Modern HIV research depends crucially on both viral sequencing and population measurements. To directly link mechanistic biological processes and evolutionary dynamics during HIV infection, we developed multiple within-host phylodynamic models of HIV primary infection for comparative validation against viral load and evolutionary dynamics data. The optimal model of primary infection required no positive selection, suggesting that the host adaptive immune system reduces viral load but surprisingly does not drive observed viral evolution. Rather, the fitness (infectivity) of mutant variants is drawn from an exponential distribution in which most variants are slightly less infectious than their parents (nearly neutral evolution). This distribution was not largely different from either in vivo fitness distributions recorded beyond primary infection or in vitro distributions that are observed without adaptive immunity, suggesting the intrinsic viral fitness distribution may drive evolution. Simulated phylogenetic trees also agree with independent data and illuminate how phylogenetic inference must consider viral and immune-cell population dynamics to gain accurate mechanistic insights.

Characterization of Near Full-Length Transmitted/Founder HIV-1 Subtype D and A/D Recombinant Genomes in a Heterosexual Ugandan Population (2006–2011)

Detailed characterization of transmitted HIV-1 variants in Uganda is fundamentally important to inform vaccine design, yet studies on the transmitted full-length strains of subtype D viruses are limited. Here, we amplified single genomes and characterized viruses, some of which were previously classified as subtype D by sub-genomic pol sequencing that were transmitted in Uganda between December 2006 to June 2011. Analysis of 5′ and 3′ half genome sequences showed 73% (19/26) of infections involved single virus transmissions, whereas 27% (7/26) of infections involved multiple variant transmissions based on predictions of a model of random virus evolution. Subtype analysis of inferred transmitted/founder viruses showed a high transmission rate of inter-subtype recombinants (69%, 20/29) involving mainly A1/D, while pure subtype D variants accounted for one-third of infections (31%, 9/29). Recombination patterns included a predominance of subtype D in the gag/pol region and a highly recombinogenic envelope gene. The signal peptide-C1 region and gp41 transmembrane domain (Tat2/Rev2 flanking region) were hotspots for A1/D recombination events. Analysis of a panel of 14 transmitted/founder molecular clones showed no difference in replication capacity between subtype D viruses (n = 3) and inter-subtype mosaic recombinants (n = 11). However, individuals infected with high replication capacity viruses had a faster CD4 T cell loss. The high transmission rate of unique inter-subtype recombinants is striking and emphasizes the extraordinary challenge for vaccine design and, in particular, for the highly variable and recombinogenic envelope gene, which is targeted by rational designs aimed to elicit broadly neutralizing antibodies.

When viruses become more virulent

[Figure: see text].

A highly virulent variant of HIV-1 circulating in the Netherlands

We discovered a highly virulent variant of subtype-B HIV-1 in the Netherlands. One hundred nine individuals with this variant had a 0.54 to 0.74 log₁₀ increase (i.e., a ~3.5-fold to 5.5-fold increase) in viral load compared with, and exhibited CD4 cell decline twice as fast as, 6604 individuals with other subtype-B strains. Without treatment, advanced HIV-CD4 cell counts below 350 cells per cubic millimeter, with long-term clinical consequences-is expected to be reached, on average, 9 months after diagnosis for individuals in their thirties with this variant. Age, sex, suspected mode of transmission, and place of birth for the aforementioned 109 individuals were typical for HIV-positive people in the Netherlands, which suggests that the increased virulence is attributable to the viral strain. Genetic sequence analysis suggests that this variant arose in the 1990s from de novo mutation, not recombination, with increased transmissibility and an unfamiliar molecular mechanism of virulence.

Second Messenger 2'3'-cyclic GMP-AMP (2'3'-cGAMP):Synthesis, transmission, and degradation

Cyclic GMP-AMP synthase (cGAS) senses foreign DNA to produce 2'3'-cyclic GMP-AMP (2'3'-cGAMP). 2'3'-cGAMP is a second messenger that binds and activates the adaptor protein STING, which triggers the innate immune response. As a STING agonist, the small molecule 2'3'-cGAMP plays pivotal roles in antiviral defense and has adjuvant applications, and anti-tumor effects. 2'3'-cGAMP and its analogs are thus putative targets for immunotherapy and are currently being testedin clinical trials to treat solid tumors. However, several barriers to further development have emerged from these studies, such as evidence of immune and inflammatory side-effects, poor pharmacokinetics, and undesirable biodistribution. Here, we review the status of 2'3'-cGAMP research and outline the role of 2'3'-cGAMP in immune signaling, adjuvant applications, and cancer immunotherapy, as well as various 2'3'-cGAMP detection methods.

Sensitivity Analysis in an Immuno-Epidemiological Vector-Host Model

Sensitivity Analysis (SA) is a useful tool to measure the impact of changes in model parameters on the infection dynamics, particularly to quantify the expected efficacy of disease control strategies. SA has only been applied to epidemic models at the population level, ignoring the effect of within-host virus-with-immune-system interactions on the disease spread. Connecting the scales from individual to population can help inform drug and vaccine development. Thus the value of understanding the impact of immunological parameters on epidemiological quantities. Here we consider an age-since-infection structured vector-host model, in which epidemiological parameters are formulated as functions of within-host virus and antibody densities, governed by an ODE system. We then use SA for these immuno-epidemiological models to investigate the impact of immunological parameters on population-level disease dynamics such as basic reproduction number, final size of the epidemic or the infectiousness at different phases of an outbreak. As a case study, we consider Rift Valley Fever Disease utilizing parameter estimations from prior studies. SA indicates that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\%$$\end{document}1% increase in within-host pathogen growth rate can lead up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8\%$$\end{document}8% increase in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal R_0,$$\end{document}R0, up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1 \%$$\end{document}1% increase in steady-state infected host abundance, and up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4\%$$\end{document}4% increase in infectiousness of hosts when the reproduction number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal R_0$$\end{document}R0 is larger than one. These significant increases in population-scale disease quantities suggest that control strategies that reduce the within-host pathogen growth can be important in reducing disease prevalence.

SARS-CoV-2 virulence evolution: Avirulence theory, immunity and trade-offs

The COVID-19 pandemic has led to a resurgence of the debate on whether host-parasite interactions should evolve towards avirulence. In this review, we first show that SARS-CoV-2 virulence is evolving, before explaining why some expect the mortality caused by the epidemic to converge towards that of human seasonal alphacoronaviruses. Leaning on existing theory, we then include viral evolution into the picture and discuss hypotheses explaining why the virulence has increased since the beginning of the pandemic. Finally, we mention some potential scenarios for the future.

RV144 vaccine imprinting constrained HIV-1 evolution following breakthrough infection

The scale of the HIV-1 epidemic underscores the need for a vaccine. The multitude of circulating HIV-1 strains together with HIV-1’s high evolvability hints that HIV-1 could adapt to a future vaccine. Here, we wanted to investigate the effect of vaccination on the evolution of the virus post-breakthrough infection. We analyzed 2,635 HIV-1 env sequences sampled up to a year post-diagnosis from 110 vaccine and placebo participants who became infected in the RV144 vaccine efficacy trial. We showed that the Env signature sites that were previously identified to distinguish vaccine and placebo participants were maintained over time. In addition, fewer sites were under diversifying selection in the vaccine group than in the placebo group. These results indicate that HIV-1 would possibly adapt to a vaccine upon its roll-out.

Virologic Outcomes Among People Living With Human Immunodeficiency Virus With High Pretherapy Viral Load Burden Initiating on Common Core Agents

Background

People living with human immunodeficiency virus (PLWH) initiating antiretroviral therapy (ART) with viral loads (VLs) ≥100 000 copies/mL are less likely to achieve virologic success, but few studies have characterized real-world treatment outcomes.

Methods

ART-naive PLWH with VLs ≥100 000 copies/mL initiating dolutegravir (DTG), elvitegravir (EVG), raltegravir (RAL), or darunavir (DRV) between 12 August 2013 and 31 July 2017 were identified from the OPERA database. Virologic failure was defined as (i) 2 consecutive VLs ≥200 copies/mL after 36 weeks of ART; (ii) 1 VL ≥200 copies/mL with core agent discontinuation after 36 weeks; (iii) 2 consecutive VLs ≥200 copies/mL after suppression (≤50 copies/mL) before 36 weeks; or (iv) 1 VL ≥200 copies/mL with discontinuation after suppression before 36 weeks. Cox modeling estimated the association between regimen and virologic failure.

Results

There were 2038 ART-naive patients with high VL who initiated DTG (36%), EVG (46%), DRV (16%), or RAL (2%). Median follow-up was 18.1 (interquartile range, 12.4–28.9) months. EVG and DTG initiators were similar at baseline, but RAL initiators were older and more likely to be female with low CD4 cell counts while DRV initiators differed notably on factors associated with treatment failure. Virologic failure was experienced by 9.2% DTG, 13.2% EVG, 18.4% RAL, and 18.8% DRV initiators. Compared to DTG, the adjusted hazard ratio (95% confidence interval) was 1.46 (1.05–2.03) for EVG, 2.24 (1.50–3.34) for DRV, and 4.13 (1.85–9.24) for RAL.

Conclusions

ART-naive PLWH with high VLs initiating on DTG were significantly less likely to experience virologic failure compared to EVG, RAL, and DRV initiators.

Antiretroviral therapy-naïve people living with HIV (PLWH) initiating therapy with viral loads ≥100,000 copies/mL varied markedly at baseline. In adjusted models, PLWH initiating dolutegravir-based regimens were less likely to experience virologic failure as compared to elvitegravir, raltegravir and darunavir initiators.

Influence of sexual risk behaviour and STI co-infection dynamics on the evolution of HIV set point viral load in MSM

HIV viral load (VL) is an important predictor of HIV progression and transmission. Anti-retroviral therapy (ART) has been reported to reduce HIV transmission by lowering VL. However, apart from this beneficial effect, increased levels of population mean set-point viral load (SPVL), an estimator for HIV virulence, have been observed in men who have sex with men (MSM) in the decade following the introduction of ART in The Netherlands. Several studies have been devoted to explain these counter-intuitive trends in SPVL. However, to our knowledge, none of these studies has investigated an explanation in which it arises as the result of a sexually transmitted infection (STI) co-factor in detail. In this study, we adapted an event-based, individual-based model to investigate how STI co-infection and sexual risk behaviour affect the evolution of HIV SPVL in MSM before and after the introduction of ART. The results suggest that sexual risk behaviour has an effect on SPVL and indicate that more data are needed to test the effect of STI co-factors on SPVL. Furthermore, the observed trends in SPVL cannot be explained by sexual risk behaviour and STI co-factors only. We recommend to develop mathematical models including also factors related to viral evolution as reported earlier in the literature. However, this requires more complex models, and the collection of more data for parameter estimation than what is currently available.

Human immunotypes impose selection on viral genotypes through viral epitope specificity

BACKGROUND

Understanding the genetic interplay between human hosts and infectious pathogens is crucial for how we interpret virulence factors. Here, we tested for associations between HIV and host genetics, and interactive genetic effects on viral load (VL) in HIV+ ART-naive clinical trial participants.

METHODS

HIV genomes were sequenced and the encoded amino acid (AA) variants were associated with VL, human single nucleotide polymorphisms (SNPs) and imputed HLA alleles, using generalized linear models with Bonferroni correction.

RESULTS

Human (388,501 SNPs) and HIV (3,010 variants) genetic data was available for 2,122 persons. Four HIV variants were associated with VL (p-values<1.66×10 -5). Twelve HIV variants were associated with a range of 1-512 human SNPs (p-value<4.28×10 -11). We found 46 associations between HLA alleles and HIV variants (p-values<1.29×10 -7). We found HIV variants and immunotypes when analyzed separately, were associated with lower VL, whereas the opposite was true when analyzed in concert. Epitope binding prediction showed HLA alleles to be weaker binders of associated HIV AA variants relative to alternative variants on the same position.

CONCLUSIONS

Our results show the importance of immunotype specificity on viral antigenic determinants, and the identified genetic interplay puts emphasis that viral and human genetics should be studied in the context of each other.

Virus shedding kinetics and unconventional virulence tradeoffs

Tradeoff theory, which postulates that virulence provides both transmission costs and benefits for pathogens, has become widely adopted by the scientific community. Although theoretical literature exploring virulence-tradeoffs is vast, empirical studies validating various assumptions still remain sparse. In particular, truncation of transmission duration as a cost of virulence has been difficult to quantify with robust controlled in vivo studies. We sought to fill this knowledge gap by investigating how transmission rate and duration were associated with virulence for infectious hematopoietic necrosis virus (IHNV) in rainbow trout (Oncorhynchus mykiss). Using host mortality to quantify virulence and viral shedding to quantify transmission, we found that IHNV did not conform to classical tradeoff theory. More virulent genotypes of the virus were found to have longer transmission durations due to lower recovery rates of infected hosts, but the relationship was not saturating as assumed by tradeoff theory. Furthermore, the impact of host mortality on limiting transmission duration was minimal and greatly outweighed by recovery. Transmission rate differences between high and low virulence genotypes were also small and inconsistent. Ultimately, more virulent genotypes were found to have the overall fitness advantage, and there was no apparent constraint on the evolution of increased virulence for IHNV. However, using a mathematical model parameterized with experimental data, it was found that host culling resurrected the virulence tradeoff and provided low virulence genotypes with the advantage. Human-induced or natural culling, as well as host population fragmentation, may be some of the mechanisms by which virulence diversity is maintained in nature. This work highlights the importance of considering non-classical virulence tradeoffs.

Functional analysis of SARS-CoV-2 proteins in Drosophila identifies Orf6-induced pathogenic effects with Selinexor as an effective treatment

BACKGROUND

SARS-CoV-2 causes COVID-19 with a widely diverse disease profile that affects many different tissues. The mechanisms underlying its pathogenicity in host organisms remain unclear. Animal models for studying the pathogenicity of SARS-CoV-2 proteins are lacking.

METHODS

Using bioinformatic analysis, we found that 90% of the virus-host interactions involve human proteins conserved in Drosophila. Therefore, we generated a series of transgenic fly lines for individual SARS-CoV-2 genes, and used the Gal4-UAS system to express these viral genes in Drosophila to study their pathogenicity.

RESULTS

We found that the ubiquitous expression of Orf6, Nsp6 or Orf7a in Drosophila led to reduced viability and tissue defects, including reduced trachea branching as well as muscle deficits resulting in a "held-up" wing phenotype and poor climbing ability. Furthermore, muscles in these flies showed dramatically reduced mitochondria. Since Orf6 was found to interact with nucleopore proteins XPO1, we tested Selinexor, a drug that inhibits XPO1, and found that it could attenuate the Orf6-induced lethality and tissue-specific phenotypes observed in flies.

CONCLUSIONS

Our study established Drosophila as a model for studying the function of SARS-CoV2 genes, identified Orf6 as a highly pathogenic protein in various tissues, and demonstrated the potential of Selinexor for inhibiting Orf6 toxicity using an in vivo animal model system.

Characterization of SARS-CoV-2 proteins reveals Orf6 pathogenicity, subcellular localization, host interactions and attenuation by Selinexor

BACKGROUND

SARS-CoV-2 causes COVID-19 which has a widely diverse disease profile. The mechanisms underlying its pathogenicity remain unclear. We set out to identify the SARS-CoV-2 pathogenic proteins that through host interactions cause the cellular damages underlying COVID-19 symptomatology.

METHODS

We examined each of the individual SARS-CoV-2 proteins for their cytotoxicity in HEK 293 T cells and their subcellular localization in COS-7 cells. We also used Mass-Spec Affinity purification to identify the host proteins interacting with SARS-CoV-2 Orf6 protein and tested a drug that could inhibit a specific Orf6 and host protein interaction.

RESULTS

We found that Orf6, Nsp6 and Orf7a induced the highest toxicity when over-expressed in human 293 T cells. All three proteins showed membrane localization in COS-7 cells. We focused on Orf6, which was most cytotoxic and localized to the endoplasmic reticulum, autophagosome and lysosomal membranes. Proteomics revealed Orf6 interacts with nucleopore proteins (RAE1, XPO1, RANBP2 and nucleoporins). Treatment with Selinexor, an FDA-approved inhibitor for XPO1, attenuated Orf6-induced cellular toxicity in human 293 T cells.

CONCLUSIONS

Our study revealed Orf6 as a highly pathogenic protein from the SARS-CoV-2 genome, identified its key host interacting proteins, and Selinexor as a drug candidate for directly targeting Orf6 host protein interaction that leads to cytotoxicity.

A Novel High Throughput, Parallel Infection Assay for Determining the Replication Capacities of 346 Primary HIV-1 Isolates of the Zurich Primary HIV-1 Infection Study in Primary Cells

HIV-1 replication capacity is an important characteristic to understand the replication competence of single variants or virus populations. It can further aid in the understanding of HIV-1 pathogenicity, disease progression, and drug resistance mutations. To effectively study RC, many assays have been established. However, there is still demand for a high throughput replication capacity assay using primary cells which is robust and reproducible. In this study, we established such an assay and validated it using 346 primary HIV-1 isolates from patients enrolled in the Zurich Primary HIV Infection study (ZPHI) and two control viruses, HIV-1 JR-CSFWT and HIV-1 JR-CSFK65R_M184V. Replication capacity was determined by measuring the viral growth on PBMCs over 10 days by longitudinally transferring cell culture supernatant to TZM-bl reporter cells. By utilizing the TZM-bl luciferase reporter assay, we determined replication capacity by measuring viral infectivity. The simplicity of the experimental setup allowed for all 346 primary HIV-1 isolates to be replicated at one time. Although the infectious input dose for each virus was normalized, a broad range of replication capacity values over 4 logs was observed. The approach was confirmed by two repeated experiments and we demonstrated that the reproducibility of the replication capacity values is statistically comparable between the two separate experiments. In summary, these results endorse our high throughput replication capacity assay as reproducible and robust and can be utilized for large scale HIV-1 replication capacity experiments in primary cells.

The influence of human genetic variation on Epstein-Barr virus sequence diversity

Epstein-Barr virus (EBV) is one of the most common viruses latently infecting humans. Little is known about the impact of human genetic variation on the large inter-individual differences observed in response to EBV infection. To search for a potential imprint of host genomic variation on the EBV sequence, we jointly analyzed paired viral and human genomic data from 268 HIV-coinfected individuals with CD4 + T cell count < 200/mm3 and elevated EBV viremia. We hypothesized that the reactivated virus circulating in these patients could carry sequence variants acquired during primary EBV infection, thereby providing a snapshot of early adaptation to the pressure exerted on EBV by the individual immune response. We searched for associations between host and pathogen genetic variants, taking into account human and EBV population structure. Our analyses revealed significant associations between human and EBV sequence variation. Three polymorphic regions in the human genome were found to be associated with EBV variation: one at the amino acid level (BRLF1:p.Lys316Glu); and two at the gene level (burden testing of rare variants in BALF5 and BBRF1). Our findings confirm that jointly analyzing host and pathogen genomes can identify sites of genomic interactions, which could help dissect pathogenic mechanisms and suggest new therapeutic avenues.

Physical 'strength' of the multi-protein chain connecting immune cells: Does the weakest link limit antibody affinity maturation?: The weakest link in the multi-protein chain facilitating antigen acquisition by B cells in germinal centres limits antibody affinity maturation

The affinities of antibodies (Abs) for their target antigens (Ags) gradually increase in vivo following an infection or vaccination, but reach saturation at values well below those realisable in vitro. This 'affinity ceiling' could in many cases restrict our ability to fight infections and compromise vaccines. What determines the affinity ceiling has been an unresolved question for decades. Here, we argue that it arises from the strength of the chain of protein complexes that is pulled by B cells during the process of Ag acquisition. The affinity ceiling is determined by the strength of the weakest link in the chain. We identify the weakest link and show that the resulting affinity ceiling can explain the Ab affinities realized in vivo, providing a conceptual understanding of Ab affinity maturation. We explore plausible evolutionary underpinnings of the affinity ceiling, examine supporting evidence and alternative hypotheses and discuss implications for vaccination strategies.

Treating symptomatic infections and the co-evolution of virulence and drug resistance

Antimicrobial therapeutic treatments are by definition applied after the onset of symptoms, which tend to correlate with infection severity. Using mathematical epidemiology models, I explore how this link affects the coevolutionary dynamics between the virulence of an infection, measured via host mortality rate, and its susceptibility to chemotherapy. I show that unless resistance pre-exists in the population, drug-resistant infections are initially more virulent than drug-sensitive ones. As the epidemic unfolds, virulence is more counter-selected in drug-sensitive than in drug-resistant infections. This difference decreases over time and, eventually, the exact shape of genetic trade-offs govern long-term evolutionary dynamics. Using adaptive dynamics, I show that two types of evolutionary stable strategies (ESS) may be reached in the context of this simple model and that, depending on the parameter values, an ESS may only be locally stable. In general, the more the treatment rate increases with virulence, the lower the ESS value. Overall, both on the short-term and long-term, having treatment rate depend on infection virulence tend to favour less virulent strains in drug-sensitive infections. These results highlight the importance of the feedbacks between epidemiology, public health policies and parasite evolution, and have implications for the monitoring of virulence evolution.

Selective Pressure-Free Treatments for COVID-19

The new severe acute respiratory syndrome coronavirus (SARS-CoV-2) has caused more than 40 million human infections since December 2019, when a cluster of unexplained pneumonia cases was first reported in Wuhan, China. Just a few days after the coronavirus was officially recognized, it was identified as the causative agent of this mysterious pneumonia. This paper discusses the pros and cons of antiviral drugs from the selective pressure and possible drug resistance point of view. We also address the key advantages of potential selective pressure-free treatment methods such as the use of sparsely and densely ionizing low-dose radiation (LDR). It is known that LDR has the capacity to modulate excessive inflammatory responses, regulate lymphocyte counts and control bacterial co-infections in patients with COVID-19 and different modalities. Substantial evidence shows that viruses are constantly mutating and evolving. When an antiviral immune response is unable to eliminate a virus, viral evolution is promoted. Therefore, it is of crucial importance to limit the use of antivirals/vaccines against SARS-CoV-2 when their effects on viral fitness are not fully understood. Furthermore, to limit the spread of the virus, it is essential to develop a vaccine that is available for as many people as possible. However, with the advent of vaccines or new therapies, the new situation may force the virus to evolve. Given this consideration, selective pressure-free treatments for COVID-19 are of great importance.

Pre-existing resistance in the latent reservoir can compromise VRC01 therapy during chronic HIV-1 infection

Passive immunization with broadly neutralizing antibodies (bNAbs) of HIV-1 appears a promising strategy for eliciting long-term HIV-1 remission. When administered concomitantly with the cessation of antiretroviral therapy (ART) to patients with established viremic control, bNAb therapy is expected to prolong remission. Surprisingly, in clinical trials on chronic HIV-1 patients, the bNAb VRC01 failed to prolong remission substantially. Identifying the cause of this failure is important for improving VRC01-based therapies and unraveling potential vulnerabilities of other bNAbs. In the trials, viremia resurged rapidly in most patients despite suppressive VRC01 concentrations in circulation, suggesting that VRC01 resistance was the likely cause of failure. ART swiftly halts viral replication, precluding the development of resistance during ART. If resistance were to emerge post ART, virological breakthrough would have taken longer than without VRC01 therapy. We hypothesized therefore that VRC01-resistant strains must have been formed before ART initiation, survived ART in latently infected cells, and been activated during VRC01 therapy, causing treatment failure. Current assays preclude testing this hypothesis experimentally. We developed a mathematical model based on the hypothesis and challenged it with available clinical data. The model integrated within-host HIV-1 evolution, stochastic latency reactivation, and viral dynamics with multiple-dose VRC01 pharmacokinetics. The model predicted that single but not higher VRC01-resistant mutants would pre-exist in the latent reservoir. We constructed a virtual patient population that parsimoniously recapitulated inter-patient variations. Model predictions with this population quantitatively captured data of VRC01 failure from clinical trials, presenting strong evidence supporting the hypothesis. We attributed VRC01 failure to single-mutant VRC01-resistant proviruses in the latent reservoir triggering viral recrudescence, particularly when VRC01 was at trough levels. Pre-existing resistant proviruses in the latent reservoir may similarly compromise other bNAbs. Our study provides a framework for designing bNAb-based therapeutic protocols that would avert such failure and maximize HIV-1 remission.

Antiretroviral therapy (ART) can control but not eradicate HIV-1. Stopping ART leads to rapid viral resurgence and progressive disease. ART is therefore administered lifelong. Tremendous efforts are ongoing to devise strategies that will enable stopping ART and yet prevent viral resurgence. One such strategy involves the administration of broadly neutralizing antibodies (bNAbs) of HIV-1 at the time of stopping ART. This strategy is expected to delay if not prevent viral resurgence. Surprisingly, treatment with VRC01, a potent bNAb, resulted in hardly any improvement in viral remission. In this study, we elucidate the cause of this failure. We hypothesized that VRC01-resistant strains may pre-exist in latently infected cells, which are unaffected by ART. They can thus outlast ART and get reactivated, triggering VRC01 failure. We built a detailed mathematical model based on this hypothesis and showed that it quantitatively captured observations of VRC01 failure in clinical trials on chronic HIV-1 patients. Our study thus identifies a potential vulnerability of bNAbs, namely, bNAb-resistant strains pre-existing in latently infected cells. Our model offers a framework for predicting bNAb-based treatment protocols that would preclude failure due to pre-existing resistance and maximally prolong remission.

Heritability of the HIV-1 reservoir size and decay under long-term suppressive ART

The HIV-1 reservoir is the major hurdle to curing HIV-1. However, the impact of the viral genome on the HIV-1 reservoir, i.e. its heritability, remains unknown. We investigate the heritability of the HIV-1 reservoir size and its long-term decay by analyzing the distribution of those traits on viral phylogenies from both partial-pol and viral near full-length genome sequences. We use a unique nationwide cohort of 610 well-characterized HIV-1 subtype-B infected individuals on suppressive ART for a median of 5.4 years. We find that a moderate but significant fraction of the HIV-1 reservoir size 1.5 years after the initiation of ART is explained by genetic factors. At the same time, we find more tentative evidence for the heritability of the long-term HIV-1 reservoir decay. Our findings indicate that viral genetic factors contribute to the HIV-1 reservoir size and hence the infecting HIV-1 strain may affect individual patients’ hurdle towards a cure.

The HIV reservoir is a major hurdle for a cure of HIV, but the factors determining its size and dynamics remain unclear. Here the authors show in a large cohort of 610 HIV-1 infected individuals, who are on suppressive ART for a median of 5.4 years, that viral genetic factors contribute substantially to the HIV-1 reservoir size.