Contemporary HIV-1 consensus Env with AI-assisted redesigned hypervariable loops promote antibody binding.

An effective HIV-1 vaccine must elicit broadly neutralizing antibodies (bnAbs) against highly diverse Envelope glycoproteins (Env). Since Env with the longest hypervariable (HV) loops is more resistant to the cognate bnAbs than Env with shorter HV loops, we redesigned hypervariable loops for updated Env consensus sequences of subtypes B and C and CRF01_AE. Using modeling with AlphaFold2, we reduced the length of V1, V2, and V5 HV loops while maintaining the integrity of the Env structure and glycan shield, and modified the V4 HV loop. Spacers are designed to limit strain-specific targeting. All updated Env are infectious as pseudoviruses. Preliminary structural characterization suggests that the modified HV loops have a limited impact on Env's conformation. Binding assays show improved binding to modified subtype B and CRF01_AE Env but not to subtype C Env. Neutralization assays show increases in sensitivity to bnAbs, although not always consistently across clades. Strikingly, the HV loop modification renders the resistant CRF01_AE Env sensitive to 10-1074 despite the absence of a glycan at N332.

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.

Evolution of HIV-1 envelope towards reduced neutralization sensitivity, as demonstrated by contemporary HIV-1 subtype B from the United States.

Subtype B HIV-1 has been the primary driver of the HIV-1 epidemic in the United States (U.S.) for over forty years and is also a prominent subtype in the Americas, Europe, Australia, the Middle East and North Africa. In this study, the neutralization profiles of contemporary subtype B Envs from the U.S. were assessed to characterize changes in neutralization sensitivities over time. We generated a panel of 30 contemporary pseudoviruses (PSVs) and demonstrated continued diversification of subtype B Env from the 1980s up to 2018. Neutralization sensitivities of the contemporary subtype B PSVs were characterized using 31 neutralizing antibodies (NAbs) and were compared with strains from earlier in the HIV-1 pandemic. A significant reduction in Env neutralization sensitivity was observed for 27 out of 31 NAbs for the contemporary as compared to earlier-decade subtype B PSVs. A decline in neutralization sensitivity was observed across all Env domains; the NAbs that were most potent early in the pandemic suffered the greatest decline in potency over time. A meta-analysis demonstrated this trend across multiple subtypes. As HIV-1 Env diversification continues, changes in Env antigenicity and neutralization sensitivity should continue to be evaluated to inform the development of improved vaccine and antibody products to prevent and treat HIV-1.

Shark nanobodies with potent SARS-CoV-2 neutralizing activity and broad sarbecovirus reactivity.

Despite rapid and ongoing vaccine and therapeutic development, SARS-CoV-2 continues to evolve and evade, presenting a need for next-generation diverse therapeutic modalities. Here we show that nurse sharks immunized with SARS-CoV-2 recombinant receptor binding domain (RBD), RBD-ferritin (RFN), or spike protein ferritin nanoparticle (SpFN) immunogens elicit a set of new antigen receptor antibody (IgNAR) molecules that target two non-overlapping conserved epitopes on the spike RBD. Representative shark antibody variable NAR-Fc chimeras (ShAbs) targeting either of the two epitopes mediate cell-effector functions, with high affinity to all SARS-CoV-2 viral variants of concern, including the divergent Omicron strains. The ShAbs potently cross-neutralize SARS-CoV-2 WA-1, Alpha, Beta, Delta, Omicron BA.1 and BA.5, and SARS-CoV-1 pseudoviruses, and confer protection against SARS-CoV-2 challenge in the K18-hACE2 transgenic mouse model. Structural definition of the RBD-ShAb01-ShAb02 complex enabled design and production of multi-specific nanobodies with enhanced neutralization capacity, and picomolar affinity to divergent sarbecovirus clade 1a, 1b and 2 RBD molecules. These shark nanobodies represent potent immunotherapeutics both for current use, and future sarbecovirus pandemic preparation.

Design of a pan-betacoronavirus vaccine candidate through a phylogenetically informed approach.

Coronaviruses are a diverse family of viruses that crossed over into humans at least seven times, precipitating mild to catastrophic outcomes. The severe acute respiratory syndrome coronavirus 2 pandemic renewed efforts to identify strains with zoonotic potential and to develop pan-coronavirus vaccines. The analysis of 2181 coronavirus genomes (from 102 host species) confirmed the limited sequence conservation across genera (alpha-, beta-, delta-, and gammacoronavirus) and proteins. A phylogenetically informed pan-coronavirus vaccine was not feasible because of high genetic heterogeneity across genera. We focused on betacoronaviruses and identified nonhuman-infecting receptor binding domain (RBD) sequences that were more genetically similar to human coronaviruses than expected given their phylogenetic divergence. These human-like RBDs defined three phylogenetic clusters. A vaccine candidate based on a representative sequence for each cluster covers the diversity estimated to protect against existing and future human-infecting betacoronaviruses. Our findings emphasize the potential value of conceptualizing prophylaxis against zoonoses in terms of genetic, rather than species, diversity.

Optimal sequence-based design for multi-antigen HIV-1 vaccines using minimally distant antigens.

The immense global diversity of HIV-1 is a significant obstacle to developing a safe and effective vaccine. We recently showed that infections established with multiple founder variants are associated with the development of neutralization breadth years later. We propose a novel vaccine design strategy that integrates the variability observed in acute HIV-1 infections with multiple founder variants. We developed a probabilistic model to simulate this variability, yielding a set of sequences that present the minimal diversity seen in an infection with multiple founders. We applied this model to a subtype C consensus sequence for the Envelope (Env) (used as input) and showed that the simulated Env sequences mimic the mutational landscape of an infection with multiple founder variants, including diversity at antibody epitopes. The derived set of multi-founder-variant-like, minimally distant antigens is designed to be used as a vaccine cocktail specific to a HIV-1 subtype or circulating recombinant form and is expected to promote the development of broadly neutralizing antibodies.

Does HIV-1 virulence matter in the ART era?

Are highly virulent HIV-1 strains circulating? Wymant et al.1 described the first concrete example of a virulent HIV-1 strain associated with high viremia and fast CD4 T cell decline. This is certainly not an isolated case, emphasizing the need for thorough HIV-1 surveillance and universal access to treatment.

HIV-1 infections with multiple founders associate with the development of neutralization breadth.

Eliciting broadly neutralizing antibodies (bnAbs) is a cornerstone of HIV-1 vaccine strategies. Comparing HIV-1 envelope (env) sequences from the first weeks of infection to the breadth of antibody responses observed several years after infection can help define viral features critical to vaccine design. We investigated the relationship between HIV-1 env genetics and the development of neutralization breadth in 70 individuals enrolled in a prospective acute HIV-1 cohort. Half of the individuals who developed bnAbs were infected with multiple HIV-1 founder variants, whereas all individuals with limited neutralization breadth had been infected with single HIV-1 founders. Accordingly, at HIV-1 diagnosis, env diversity was significantly higher in participants who later developed bnAbs compared to those with limited breadth (p = 0.012). This association between founder multiplicity and the subsequent development of neutralization breadth was also observed in 56 placebo recipients in the RV144 vaccine efficacy trial. In addition, we found no evidence that neutralization breath was heritable when analyzing env sequences from the 126 participants. These results demonstrate that the presence of slightly different HIV-1 variants in acute infection could promote the induction of bnAbs, suggesting a novel vaccine strategy, whereby an initial immunization with a cocktail of minimally distant antigens would be able to initiate bnAb development towards breadth.

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.

Monocyte-derived transcriptome signature indicates antibody-dependent cellular phagocytosis as a potential mechanism of vaccine-induced protection against HIV-1.

A gene signature was previously found to be correlated with mosaic adenovirus 26 vaccine protection in simian immunodeficiency virus and simian-human immunodeficiency virus challenge models in non-human primates. In this report, we investigated the presence of this signature as a correlate of reduced risk in human clinical trials and potential mechanisms of protection. The absence of this gene signature in the DNA/rAd5 human vaccine trial, which did not show efficacy, strengthens our hypothesis that this signature is only enriched in studies that demonstrated protection. This gene signature was enriched in the partially effective RV144 human trial that administered the ALVAC/protein vaccine, and we find that the signature associates with both decreased risk of HIV-1 acquisition and increased vaccine efficacy (VE). Total RNA-seq in a clinical trial that used the same vaccine regimen as the RV144 HIV vaccine implicated antibody-dependent cellular phagocytosis (ADCP) as a potential mechanism of vaccine protection. CITE-seq profiling of 53 surface markers and transcriptomes of 53,777 single cells from the same trial showed that genes in this signature were primarily expressed in cells belonging to the myeloid lineage, including monocytes, which are major effector cells for ADCP. The consistent association of this transcriptome signature with VE represents a tool both to identify potential mechanisms, as with ADCP here, and to screen novel approaches to accelerate the development of new vaccine candidates.

Factors influencing estimates of HIV-1 infection timing using BEAST.

While large datasets of HIV-1 sequences are increasingly being generated, many studies rely on a single gene or fragment of the genome and few comparative studies across genes have been done. We performed genome-based and gene-specific Bayesian phylogenetic analyses to investigate how certain factors impact estimates of the infection dates in an acute HIV-1 infection cohort, RV217. In this cohort, HIV-1 diagnosis corresponded to the first RNA positive test and occurred a median of four days after the last negative test, allowing us to compare timing estimates using BEAST to a narrow window of infection. We analyzed HIV-1 sequences sampled one week, one month and six months after HIV-1 diagnosis in 39 individuals. We found that shared diversity and temporal signal was limited in acute infection, and insufficient to allow timing inferences in the shortest HIV-1 genes, thus dated phylogenies were primarily analyzed for env, gag, pol and near full-length genomes. There was no one best-fitting model across participants and genes, though relaxed molecular clocks (73% of best-fitting models) and the Bayesian skyline (49%) tended to be favored. For infections with single founders, the infection date was estimated to be around one week pre-diagnosis for env (IQR: 3-9 days) and gag (IQR: 5-9 days), whilst the genome placed it at a median of 10 days (IQR: 4-19). Multiply-founded infections proved problematic to date. Our ability to compare timing inferences to precise estimates of HIV-1 infection (within a week) highlights that molecular dating methods can be applied to within-host datasets from early infection. Nonetheless, our results also suggest caution when using uniform clock and population models or short genes with limited information content.

A SARS-CoV-2 vaccine candidate would likely match all currently circulating variants.

The magnitude of the COVID-19 pandemic underscores the urgency for a safe and effective vaccine. Many vaccine candidates focus on the Spike protein, as it is targeted by neutralizing antibodies and plays a key role in viral entry. Here we investigate the diversity seen in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequences and compare it to the sequence on which most vaccine candidates are based. Using 18,514 sequences, we perform phylogenetic, population genetics, and structural bioinformatics analyses. We find limited diversity across SARS-CoV-2 genomes: Only 11 sites show polymorphisms in >5% of sequences; yet two mutations, including the D614G mutation in Spike, have already become consensus. Because SARS-CoV-2 is being transmitted more rapidly than it evolves, the viral population is becoming more homogeneous, with a median of seven nucleotide substitutions between genomes. There is evidence of purifying selection but little evidence of diversifying selection, with substitution rates comparable across structural versus nonstructural genes. Finally, the Wuhan-Hu-1 reference sequence for the Spike protein, which is the basis for different vaccine candidates, matches optimized vaccine inserts, being identical to an ancestral sequence and one mutation away from the consensus. While the rapid spread of the D614G mutation warrants further study, our results indicate that drift and bottleneck events can explain the minimal diversity found among SARS-CoV-2 sequences. These findings suggest that a single vaccine candidate should be efficacious against currently circulating lineages.

Molecular dating and viral load growth rates suggested that the eclipse phase lasted about a week in HIV-1 infected adults in East Africa and Thailand.

Most HIV-1 infected individuals do not know their infection dates. Precise infection timing is crucial information for studies that document transmission networks or drug levels at infection. To improve infection timing, we used the prospective RV217 cohort where the window when plasma viremia becomes detectable is narrow: the last negative visit occurred a median of four days before the first detectable HIV-1 viremia with an RNA test, referred below as diagnosis. We sequenced 1,280 HIV-1 genomes from 39 participants at a median of 4, 32 and 170 days post-diagnosis. HIV-1 infections were dated by using sequence-based methods and a viral load regression method. Bayesian coalescent and viral load regression estimated that infections occurred a median of 6 days prior to diagnosis (IQR: 9-3 and 11-4 days prior, respectively). Poisson-Fitter, which analyzes the distribution of hamming distances among sequences, estimated a median of 7 days prior to diagnosis (IQR: 15-4 days) based on sequences sampled 4 days post-diagnosis, but it did not yield plausible results using sequences sampled at 32 days. Fourteen participants reported a high-risk exposure event at a median of 8 days prior to diagnosis (IQR: 12 to 6 days prior). These different methods concurred that HIV-1 infection occurred about a week before detectable viremia, corresponding to 20 days (IQR: 34-15 days) before peak viral load. Together, our methods comparison helps define a framework for future dating studies in early HIV-1 infection.

Characterizing and Comparing Phylogenetic Trait Data from Their Normalized Laplacian Spectrum.

The dissection of the mode and tempo of phenotypic evolution is integral to our understanding of global biodiversity. Our ability to infer patterns of phenotypes across phylogenetic clades is essential to how we infer the macroevolutionary processes governing those patterns. Many methods are already available for fitting models of phenotypic evolution to data. However, there is currently no comprehensive nonparametric framework for characterizing and comparing patterns of phenotypic evolution. Here, we build on a recently introduced approach for using the phylogenetic spectral density profile (SDP) to compare and characterize patterns of phylogenetic diversification, in order to provide a framework for nonparametric analysis of phylogenetic trait data. We show how to construct the SDP of trait data on a phylogenetic tree from the normalized graph Laplacian. We demonstrate on simulated data the utility of the SDP to successfully cluster phylogenetic trait data into meaningful groups and to characterize the phenotypic patterning within those groups. We furthermore demonstrate how the SDP is a powerful tool for visualizing phenotypic space across traits and for assessing whether distinct trait evolution models are distinguishable on a given empirical phylogeny. We illustrate the approach in two empirical data sets: a comprehensive data set of traits involved in song, plumage, and resource-use in tanagers, and a high-dimensional data set of endocranial landmarks in New World monkeys. Considering the proliferation of morphometric and molecular data collected across the tree of life, we expect this approach will benefit big data analyses requiring a comprehensive and intuitive framework.

A non-parametric analytic framework for within-host viral phylogenies and a test for HIV-1 founder multiplicity.

Phylogenetics is a powerful tool for understanding the diversification dynamics of viral pathogens. Here we present an extension of the spectral density profile of the modified graph Laplacian, which facilitates the characterization of within-host molecular evolution of viruses and the direct comparison of diversification dynamics between hosts. This approach is non-parametric and therefore fast and model-free. We used simulations of within-host evolutionary scenarios to evaluate the efficiency of our approach and to demonstrate the significance of interpreting a viral phylogeny by its spectral density profile in terms of diversification dynamics. The key features that are captured by the profile are positive selection on the viral gene (or genome), temporal changes in substitution rates, mutational fitness, and time between sampling. Using sequences from individuals infected with HIV-1, we showed the utility of this approach for characterizing within-host diversification dynamics, for comparing dynamics between hosts, and for charting disease progression in infected individuals sampled over multiple years. We furthermore propose a heuristic test for assessing founder heterogeneity, which allows us to classify infections with single and multiple HIV-1 founder viruses. This non-parametric approach can be a valuable complement to existing parametric approaches.

Author Correction: Clade-specific diversification dynamics of marine diatoms since the Jurassic.

In the version of this Article originally published, the authors did not give credit to David G. Mann for the four microscopic images used in Fig. 1a. This has now been amended in all versions of the Article.

Clade-specific diversification dynamics of marine diatoms since the Jurassic.

Diatoms are one of the most abundant and diverse groups of phytoplankton and play a major role in marine ecosystems and the Earth's biogeochemical cycles. Here we combine DNA metabarcoding data from the Tara Oceans expedition with palaeoenvironmental data and phylogenetic models of diversification to analyse the diversity dynamics of marine diatoms. We reveal a primary effect of variation in carbon dioxide partial pressure (pCO2) on early diatom diversification, followed by a major burst of diversification in the late Eocene epoch, after which diversification is chiefly affected by sea level, an influx of silica availability and competition with other planktonic groups. Our results demonstrate a remarkable heterogeneity of diversification dynamics across diatoms and suggest that a changing climate will favour some clades at the expense of others.

Inferring Evolutionary Process From Neuroanatomical Data.

Brain evolution has interested neuroanatomists for over a century. These interests often fall on how free the brain is to evolve independently of the body, how free brain regions are to evolve independently of each other, and how different environmental and ecological factors affect the brain over evolutionary time. But despite major advances in phylogenetic methods, comparative neuroanatomists have tended to limit their macroevolutionary toolbox to regression-based analyses and ignored the scope of evolutionary process-based models at their disposal. This Review summarizes the history of comparative neuroanatomy and highlights the pitfalls of the methodologies traditionally used. It provides an overview of evolutionary process-based modeling approaches for investigating univariate and multivariate data, as well as more sophisticated methods that incorporate hypotheses about biotic and abiotic pressures that may drive brain evolution. The benefits of evolutionary process-based models, and shortcomings of regression-based ones, are illustrated with widely used neuroanatomical data. Ultimately, the intent of this Review is to be a guide for subsuming macroevolutionary methods not typically used in comparative neuroanatomy, in order to improve our understanding of how the brain evolves.

Detecting Environment-Dependent Diversification From Phylogenies: A Simulation Study and Some Empirical Illustrations.

Understanding the relative influence of various abiotic and biotic variables on diversification dynamics is a major goal of macroevolutionary studies. Recently, phylogenetic approaches have been developed that make it possible to estimate the role of various environmental variables on diversification using time-calibrated species trees, paleoenvironmental data, and maximum-likelihood techniques. These approaches have been effectively employed to estimate how speciation and extinction rates vary with key abiotic variables, such as temperature and sea level, and we can anticipate that they will be increasingly used in the future. Here we compile a series of biotic and abiotic paleodatasets that can be used as explanatory variables in these models and use simulations to assess the statistical properties of the approach when applied to these paleodatasets. We demonstrate that environment-dependent models perform well in recovering environment-dependent speciation and extinction parameters, as well as in correctly identifying the simulated environmental model when speciation is environment-dependent. We explore how the strength of the environment-dependency, tree size, missing taxa, and characteristics of the paleoenvironmental curves influence the performance of the models. Finally, using these models, we infer environment-dependent diversification in two empirical phylogenies: temperature-dependence in Cetacea and $\delta^{13}C$-dependence in Ruminantia. We illustrate how to evaluate the relative importance of abiotic and biotic variables in these two clades and interpret these results in light of macroevolutionary hypotheses. Given the important role paleoenvironments are presumed to have played in species evolution, our statistical assessment of how environment-dependent models behave is crucial for their utility in macroevolutionary analysis.