Tracing HIV-1 strains that imprint broadly neutralizing antibody responses.

Understanding the determinants of broadly neutralizing antibody (bNAb) evolution is crucial for the development of bNAb-based HIV vaccines1. Despite emerging information on cofactors that promote bNAb evolution in natural HIV-1 infections, in which the induction of bNAbs is genuinely rare2, information on the impact of the infecting virus strain on determining the breadth and specificity of the antibody responses to HIV-1 is lacking. Here we analyse the influence of viral antigens in shaping antibody responses in humans. We call the ability of a virus strain to induce similar antibody responses across different hosts its antibody-imprinting capacity, which from an evolutionary biology perspective corresponds to the viral heritability of the antibody responses. Analysis of 53 measured parameters of HIV-1-binding and neutralizing antibody responses in a cohort of 303 HIV-1 transmission pairs (individuals who harboured highly related HIV-1 strains and were putative direct transmission partners or members of an HIV-1 transmission chain) revealed that the effect of the infecting virus on the outcome of the bNAb response is moderate in magnitude but highly significant. We introduce the concept of bNAb-imprinting viruses and provide evidence for the existence of such viruses in a systematic screening of our cohort. The bNAb-imprinting capacity can be substantial, as indicated by a transmission pair with highly similar HIV-1 antibody responses and strong bNAb activity. Identification of viruses that have bNAb-imprinting capacities and their characterization may thus provide the potential to develop lead immunogens.

Ethnic homophily affects vaccine prioritization strategies.

People are more likely to interact with other people of their ethnicity-a phenomenon known as ethnic homophily. In the United States, people of color are known to hold proportionately more high-contact jobs and are thus more at risk of virus infection. At the same time, these ethnic groups are on average younger than the rest of the population. This gives rise to interesting disease dynamics and non-trivial trade-offs that should be taken into consideration when developing prioritization strategies for future mass vaccine roll-outs. Here, we study the spread of COVID-19 through the US population, stratified by age, ethnicity, and occupation, using a detailed, previously-developed compartmental disease model. Based on historic data from the US mass COVID-19 vaccine roll-out that began in December 2020, we show, (i) how ethnic homophily affects the choice of optimal vaccine allocation strategy, (ii) that, notwithstanding potential ethical concerns, differentiating by ethnicity in these strategies can improve outcomes (e.g., fewer deaths), and (iii) that the most likely social context in the United States is very different from the standard assumptions made by models which do not account for ethnicity and this difference affects which allocation strategy is optimal. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

Phylogenetic Cluster Analysis Identifies Virological and Behavioral Drivers of Human Immunodeficiency Virus Transmission in Men Who Have Sex With Men.

BACKGROUND: Identifying local outbreaks and their drivers is a key step toward curbing human immunodeficiency virus (HIV) transmission and potentially achieving HIV elimination. Such outbreaks can be identified as transmission clusters extracted from phylogenetic trees constructed of densely sampled viral sequences. In this study, we combined phylogenetic transmission clusters with extensive data on virological suppression and behavioral risk of cluster members to quantify the drivers of ongoing transmission over 10 years. METHODS: Using the comprehensive Swiss HIV Cohort Study and its drug-resistance database, we reconstructed phylogenetic trees for each year between 2007 and 2017. We identified HIV transmission clusters dominated by men who have sex with men (MSM) and determined their annual growth. We used Poisson regression to assess if cluster growth was associated with a per-cluster infectivity and behavioral risk score. RESULTS: Both infectivity and behavioral risk scores were significantly higher in growing MSM transmission clusters compared to nongrowing clusters (P ≤ .01). The fraction of transmission clusters without infectious members acquiring new infections increased significantly over the study period. The infectivity score was significantly associated with per-capita incidence of MSM transmission clusters in 8 years, while the behavioral risk score was significantly associated with per-capita incidence of MSM transmission clusters in 3 years. CONCLUSIONS: We present a phylogenetic method to identify hotspots of ongoing transmission among MSM. Our results demonstrate the effectiveness of treatment as prevention at the population level. However, the significantly increasing number of new infections among transmission clusters without infectious members highlights a relative shift from diagnosed to undiagnosed individuals as drivers of HIV transmission in Swiss MSM.

A model-based evaluation of the efficacy of COVID-19 social distancing, testing and hospital triage policies.

A stochastic compartmental network model of SARS-CoV-2 spread explores the simultaneous effects of policy choices in three domains: social distancing, hospital triaging, and testing. Considering policy domains together provides insight into how different policy decisions interact. The model incorporates important characteristics of COVID-19, the disease caused by SARS-CoV-2, such as heterogeneous risk factors and asymptomatic transmission, and enables a reliable qualitative comparison of policy choices despite the current uncertainty in key virus and disease parameters. Results suggest possible refinements to current policies, including emphasizing the need to reduce random encounters more than personal contacts, and testing low-risk symptomatic individuals before high-risk symptomatic individuals. The strength of social distancing of symptomatic individuals affects the degree to which asymptomatic cases drive the epidemic as well as the level of population-wide contact reduction needed to keep hospitals below capacity. The relative importance of testing and triaging also depends on the overall level of social distancing.

Phenotypic deficits in the HIV-1 envelope are associated with the maturation of a V2-directed broadly neutralizing antibody lineage.

Broadly neutralizing antibodies (bnAbs) to HIV-1 can evolve after years of an iterative process of virus escape and antibody adaptation that HIV-1 vaccine design seeks to mimic. To enable this, properties that render HIV-1 envelopes (Env) capable of eliciting bnAb responses need to be defined. Here, we followed the evolution of the V2 apex directed bnAb lineage VRC26 in the HIV-1 subtype C superinfected donor CAP256 to investigate the phenotypic changes of the virus populations circulating before and during the early phases of bnAb induction. Longitudinal viruses that evolved from the VRC26-resistant primary infecting (PI) virus, the VRC26-sensitive superinfecting (SU) virus and ensuing PI-SU recombinants revealed substantial phenotypic changes in Env, with a switch in Env properties coinciding with early resistance to VRC26. Decreased sensitivity of SU-like viruses to VRC26 was linked with reduced infectivity, altered entry kinetics and lower sensitivity to neutralization after CD4 attachment. VRC26 maintained neutralization activity against cell-associated CAP256 virus, indicating that escape through the cell-cell transmission route is not a dominant escape pathway. Reduced fitness of the early escape variants and sustained sensitivity in cell-cell transmission are both features that limit virus replication, thereby impeding rapid escape. This supports a scenario where VRC26 allowed only partial viral escape for a prolonged period, possibly increasing the time window for bnAb maturation. Collectively, our data highlight the phenotypic plasticity of the HIV-1 Env in evading bnAb pressure and the need to consider phenotypic traits when selecting and designing Env immunogens. Combinations of Env variants with differential phenotypic patterns and bnAb sensitivity, as we describe here for CAP256, may maximize the potential for inducing bnAb responses by vaccination.

Importance of routine viral load monitoring: higher levels of resistance at ART failure in Uganda and Lesotho compared with Switzerland.

Objectives: Emerging resistance to antiretroviral drugs may jeopardize the achievements of improved access to ART. We compared the prevalence of different resistance mutations in HIV-infected adults with virological failure in a cohort with regular routine viral load (VL) monitoring (Switzerland) and cohorts with limited access to VL testing (Uganda and Lesotho). Methods: We considered individuals who had genotypic resistance testing (GRT) upon virological failure (≥1000 copies/mL) and were on ART consisting of at least one NNRTI and two NRTIs. From Lesotho, individuals with two subsequent VLs ≥1000 copies/mL despite enhanced adherence counselling (n = 58) were included in the analysis. From Uganda, individuals with a single VL ≥1000 copies/mL (n = 120) were included in the analysis. From the Swiss HIV Cohort Study (SHCS), a population without history of monotherapy or dual therapy with the first GRT upon virological failure (n = 61) was selected. Results: We found that 50.8% of individuals in the SHCS, 72.5% in Uganda and 81.0% in Lesotho harboured HIV with high-level resistance to at least two drugs from their current regimen. Stanford resistance scores were higher in Uganda compared with Switzerland for all drugs used in first-line treatment except zidovudine and tenofovir (P < 0.01) and higher in Lesotho compared with Uganda for all drugs used in first-line treatment except zidovudine (P < 0.01). Conclusions: Frequent VL monitoring and possibly pretreatment GRT as done in the SHCS pays off by low levels of resistance even when treatment failure occurs. The high-level resistance patterns in Lesotho compared with Uganda could reflect a selection of strains with multiple resistance during enhanced adherence counselling.

Delineating CD4 dependency of HIV-1: Adaptation to infect low level CD4 expressing target cells widens cellular tropism but severely impacts on envelope functionality.

A hallmark of HIV-1 infection is the continuously declining number of the virus' predominant target cells, activated CD4+ T cells. With diminishing CD4+ T cell levels, the capacity to utilize alternate cell types and receptors, including cells that express low CD4 receptor levels such as macrophages, thus becomes crucial. To explore evolutionary paths that allow HIV-1 to acquire a wider host cell range by infecting cells with lower CD4 levels, we dissected the evolution of the envelope-CD4 interaction under in vitro culture conditions that mimicked the decline of CD4high target cells, using a prototypic subtype B, R5-tropic strain. Adaptation to CD4low targets proved to severely alter envelope functions including trimer opening as indicated by a higher affinity to CD4 and loss in shielding against neutralizing antibodies. We observed a strikingly decreased infectivity on CD4high target cells, but sustained infectivity on CD4low targets, including macrophages. Intriguingly, the adaptation to CD4low targets altered the kinetic of the entry process, leading to rapid CD4 engagement and an extended transition time between CD4 and CCR5 binding during entry. This phenotype was also observed for certain central nervous system (CNS) derived macrophage-tropic viruses, highlighting that the functional perturbation we defined upon in vitro adaptation to CD4low targets occurs in vivo. Collectively, our findings suggest that CD4low adapted envelopes may exhibit severe deficiencies in entry fitness and shielding early in their evolution. Considering this, adaptation to CD4low targets may preferentially occur in a sheltered and immune-privileged environment such as the CNS to allow fitness restoring compensatory mutations to occur.

Parent-offspring regression to estimate the heritability of an HIV-1 trait in a realistic setup.

BACKGROUND: Parent-offspring (PO) regression is a central tool to determine the heritability of phenotypic traits; i.e., the relative extent to which those traits are controlled by genetic factors. The applicability of PO regression to viral traits is unclear because the direction of viral transmission-who is the donor (parent) and who is the recipient (offspring)-is typically unknown and viral phylogenies are sparsely sampled. METHODS: We assessed the applicability of PO regression in a realistic setting using Ornstein-Uhlenbeck simulated data on phylogenies built from 11,442 Swiss HIV Cohort Study (SHCS) partial pol sequences and set-point viral load (SPVL) data from 3293 patients. RESULTS: We found that the misidentification of donor and recipient plays a minor role in estimating heritability and showed that sparse sampling does not influence the mean heritability estimated by PO regression. A mixed-effect model approach yielded the same heritability as PO regression but could be extended to clusters of size greater than 2 and allowed for the correction of confounding effects. Finally, we used both methods to estimate SPVL heritability in the SHCS. We employed a wide range of transmission pair criteria to measure heritability and found a strong dependence of the heritability estimates to these criteria. For the most conservative genetic distance criteria, for which heritability estimates are conceptually expected to be closest to true heritability, we found estimates ranging from 32 to 46% across different bootstrap criteria. For less conservative distance criteria, we found estimates ranging down to 8%. All estimates did not change substantially after adjusting for host-demographic factors in the mixed-effect model (±2%). CONCLUSIONS: For conservative transmission pair criteria, both PO regression and mixed-effect models are flexible and robust tools to estimate the contribution of viral genetic effects to viral traits under real-world settings. Overall, we find a strong effect of viral genetics on SPVL that is not confounded by host demographics.

Canalization reduces the nonlinearity of regulation in biological networks.

Biological networks, such as gene regulatory networks, possess desirable properties. They are more robust and controllable than random networks. This motivates the search for structural and dynamical features that evolution has incorporated into biological networks. A recent meta-analysis of published, expert-curated Boolean biological network models has revealed several such features, often referred to as design principles. Among others, the biological networks are enriched for certain recurring network motifs, the dynamic update rules are more redundant, more biased, and more canalizing than expected, and the dynamics of biological networks are better approximable by linear and lower-order approximations than those of comparable random networks. Since most of these features are interrelated, it is paramount to disentangle cause and effect, that is, to understand which features evolution actively selects for, and thus truly constitute evolutionary design principles. Here, we compare published Boolean biological network models with different ensembles of null models and show that the abundance of canalization in biological networks can almost completely explain their recently postulated high approximability. Moreover, an analysis of random N-K Kauffman models reveals a strong dependence of approximability on the dynamical robustness of a network.

Diffusion and Chemical Degradation of Vitamin B6 in Chickpeas (Cicer arietinum L.) during Hydrothermal Treatments: A Kinetic Approach.

Chickpeas are more sustainable than other food systems and have high a nutritional value, especially regarding their vitamin composition. One of the main vitamins in chickpeas is vitamin B6, which is very important for several human metabolic functions. Since chickpeas are consumed after cooking, our goal was to better understand the role of leaching (diffusion) and thermal degradation of vitamin B6 in chickpeas during hydrothermal processing. Kinetics were conducted at four temperatures, ranging from 25 to 85 °C, carried out for 4 h in an excess of water for the diffusion kinetics, or in hermetic bags for the thermal degradation kinetics. Thermal degradation was modeled according to a first-order reaction, and diffusion was modeled according to a modified version of Fick's second law. Diffusivity constants varied from 4.76 × 10-14 m2/s at 25 °C to 2.07 × 10-10 m2/s at 85 °C; the temperature had an impact on both the diffusivity constant and the residual vitamin B6. The kinetic constant ranged from 9.35 × 10-6 at 25 °C to 54.9 × 10-6 s-1 at 85 °C, with a lower impact of the temperature. In conclusion, vitamin B6 is relatively stable to heat degradation; loss is mainly due to diffusion, especially during shorter treatment times.

Learning from the COVID-19 pandemic: A systematic review of mathematical vaccine prioritization models.

As the world becomes ever more connected, the chance of pandemics increases as well. The recent COVID-19 pandemic and the concurrent global mass vaccine roll-out provides an ideal setting to learn from and refine our understanding of infectious disease models for better future preparedness. In this review, we systematically analyze and categorize mathematical models that have been developed to design optimal vaccine prioritization strategies of an initially limited vaccine. As older individuals are disproportionately affected by COVID-19, the focus is on models that take age explicitly into account. The lower mobility and activity level of older individuals gives rise to non-trivial trade-offs. Secondary research questions concern the optimal time interval between vaccine doses and spatial vaccine distribution. This review showcases the effect of various modeling assumptions on model outcomes. A solid understanding of these relationships yields better infectious disease models and thus public health decisions during the next pandemic.

Learning from the COVID-19 pandemic: a systematic review of mathematical vaccine prioritization models.

As the world becomes ever more connected, the chance of pandemics increases as well. The recent COVID-19 pandemic and the concurrent global mass vaccine roll-out provides an ideal setting to learn from and refine our understanding of infectious disease models for better future preparedness. In this review, we systematically analyze and categorize mathematical models that have been developed to design optimal vaccine prioritization strategies of an initially limited vaccine. As older individuals are disproportionately affected by COVID-19, the focus is on models that take age explicitly into account. The lower mobility and activity level of older individuals gives rise to non-trivial trade-offs. Secondary research questions concern the optimal time interval between vaccine doses and spatial vaccine distribution. This review showcases the effect of various modeling assumptions on model outcomes. A solid understanding of these relationships yields better infectious disease models and thus public health decisions during the next pandemic.

A meta-analysis of Boolean network models reveals design principles of gene regulatory networks.

Gene regulatory networks (GRNs) play a central role in cellular decision-making. Understanding their structure and how it impacts their dynamics constitutes thus a fundamental biological question. GRNs are frequently modeled as Boolean networks, which are intuitive, simple to describe, and can yield qualitative results even when data are sparse. We assembled the largest repository of expert-curated Boolean GRN models. A meta-analysis of this diverse set of models reveals several design principles. GRNs exhibit more canalization, redundancy, and stable dynamics than expected. Moreover, they are enriched for certain recurring network motifs. This raises the important question why evolution favors these design mechanisms.

Modular Control of Biological Networks.

The concept of control is central to understanding and applications of biological network models. Some of their key structural features relate to control functions, through gene regulation, signaling, or metabolic mechanisms, and computational models need to encode these. Applications of models often focus on model-based control, such as in biomedicine or metabolic engineering. This paper presents an approach to model-based control that exploits two common features of biological networks, namely their modular structure and canalizing features of their regulatory mechanisms. The paper focuses on intracellular regulatory networks, represented by Boolean network models. A main result of this paper is that control strategies can be identified by focusing on one module at a time. This paper also presents a criterion based on canalizing features of the regulatory rules to identify modules that do not contribute to network control and can be excluded. For even moderately sized networks, finding global control inputs is computationally very challenging. The modular approach presented here leads to a highly efficient approach to solving this problem. This approach is applied to a published Boolean network model of blood cancer large granular lymphocyte (T-LGL) leukemia to identify a minimal control set that achieves a desired control objective.

Projecting social contact matrices to populations stratified by binary attributes with known homophily.

Contact networks are heterogeneous. People with similar characteristics are more likely to interact, a phenomenon called assortative mixing or homophily. Empirical age-stratified social contact matrices have been derived by extensive survey work. We lack however similar empirical studies that provide social contact matrices for a population stratified by attributes beyond age, such as gender, sexual orientation, or ethnicity. Accounting for heterogeneities with respect to these attributes can have a profound effect on model dynamics. Here, we introduce a new method, which uses linear algebra and non-linear optimization, to expand a given contact matrix to populations stratified by binary attributes with a known level of homophily. Using a standard epidemiological model, we highlight the effect homophily can have on model dynamics, and conclude by briefly describing more complicated extensions. The available Python source code enables any modeler to account for the presence of homophily with respect to binary attributes in contact patterns, ultimately yielding more accurate predictive models.

Modularity of biological systems: a link between structure and function.

This paper addresses two topics in systems biology, the hypothesis that biological systems are modular and the problem of relating structure and function of biological systems. The focus here is on gene regulatory networks, represented by Boolean network models, a commonly used tool. Most of the research on gene regulatory network modularity has focused on network structure, typically represented through either directed or undirected graphs. But since gene regulation is a highly dynamic process as it determines the function of cells over time, it is natural to consider functional modularity as well. One of the main results is that the structural decomposition of a network into modules induces an analogous decomposition of the dynamic structure, exhibiting a strong relationship between network structure and function. An extensive simulation study provides evidence for the hypothesis that modularity might have evolved to increase phenotypic complexity while maintaining maximal dynamic robustness to external perturbations.

Modularity of biological systems: a link between structure and function.

This paper addresses two topics in systems biology, the hypothesis that biological systems are modular and the problem of relating structure and function of biological systems. The focus here is on gene regulatory networks, represented by Boolean network models, a commonly used tool. Most of the research on gene regulatory network modularity has focused on network structure, typically represented through either directed or undirected graphs. But since gene regulation is a highly dynamic process as it determines the function of cells over time, it is natural to consider functional modularity as well. One of the main results is that the structural decomposition of a network into modules induces an analogous decomposition of the dynamic structure, exhibiting a strong relationship between network structure and function. An extensive simulation study provides evidence for the hypothesis that modularity might have evolved to increase phenotypic complexity while maintaining maximal dynamic robustness to external perturbations.

Effect of homophily and correlation of beliefs on COVID-19 and general infectious disease outbreaks.

Contact between people with similar opinions and characteristics occurs at a higher rate than among other people, a phenomenon known as homophily. The presence of clusters of unvaccinated people has been associated with increased incidence of infectious disease outbreaks despite high population-wide vaccination rates. The epidemiological consequences of homophily regarding other beliefs as well as correlations among beliefs or circumstances are poorly understood, however. Here, we use a simple compartmental disease model as well as a more complex COVID-19 model to study how homophily and correlation of beliefs and circumstances in a social interaction network affect the probability of disease outbreak and COVID-19-related mortality. We find that the current social context, characterized by the presence of homophily and correlations between who vaccinates, who engages in risk reduction, and individual risk status, corresponds to a situation with substantially worse disease burden than in the absence of heterogeneities. In the presence of an effective vaccine, the effects of homophily and correlation of beliefs and circumstances become stronger. Further, the optimal vaccination strategy depends on the degree of homophily regarding vaccination status as well as the relative level of risk mitigation high- and low-risk individuals practice. The developed methods are broadly applicable to any investigation in which node attributes in a graph might reasonably be expected to cluster or exhibit correlations.

Evaluation of the United States COVID-19 vaccine allocation strategy.

BACKGROUND: Anticipating an initial shortage of vaccines for COVID-19, the Centers for Disease Control (CDC) in the United States developed priority vaccine allocations for specific demographic groups in the population. This study evaluates the performance of the CDC vaccine allocation strategy with respect to multiple potentially competing vaccination goals (minimizing mortality, cases, infections, and years of life lost (YLL)), under the same framework as the CDC allocation: four priority vaccination groups and population demographics stratified by age, comorbidities, occupation and living condition (congested or non-congested). METHODS AND FINDINGS: We developed a compartmental disease model that incorporates key elements of the current pandemic including age-varying susceptibility to infection, age-varying clinical fraction, an active case-count dependent social distancing level, and time-varying infectivity (accounting for the emergence of more infectious virus strains). The CDC allocation strategy is compared to all other possibly optimal allocations that stagger vaccine roll-out in up to four phases (17.5 million strategies). The CDC allocation strategy performed well in all vaccination goals but never optimally. Under the developed model, the CDC allocation deviated from the optimal allocations by small amounts, with 0.19% more deaths, 4.0% more cases, 4.07% more infections, and 0.97% higher YLL, than the respective optimal strategies. The CDC decision to not prioritize the vaccination of individuals under the age of 16 was optimal, as was the prioritization of health-care workers and other essential workers over non-essential workers. Finally, a higher prioritization of individuals with comorbidities in all age groups improved outcomes compared to the CDC allocation. CONCLUSION: The developed approach can be used to inform the design of future vaccine allocation strategies in the United States, or adapted for use by other countries seeking to optimize the effectiveness of their vaccine allocation strategies.

Distinct conformations of the HIV-1 V3 loop crown are targetable for broad neutralization.

The V3 loop of the HIV-1 envelope (Env) protein elicits a vigorous, but largely non-neutralizing antibody response directed to the V3-crown, whereas rare broadly neutralizing antibodies (bnAbs) target the V3-base. Challenging this view, we present V3-crown directed broadly neutralizing Designed Ankyrin Repeat Proteins (bnDs) matching the breadth of V3-base bnAbs. While most bnAbs target prefusion Env, V3-crown bnDs bind open Env conformations triggered by CD4 engagement. BnDs achieve breadth by focusing on highly conserved residues that are accessible in two distinct V3 conformations, one of which resembles CCR5-bound V3. We further show that these V3-crown conformations can, in principle, be attacked by antibodies. Supporting this conclusion, analysis of antibody binding activity in the Swiss 4.5 K HIV-1 cohort (n = 4,281) revealed a co-evolution of V3-crown reactivities and neutralization breadth. Our results indicate a role of V3-crown responses and its conformational preferences in bnAb development to be considered in preventive and therapeutic approaches.