Search

An abstract of the conference paper "Identification of a clade A HIV envelope immunogen from Protocol G that elicits neutralizing antibodies to tier 2 viruses," by S. Hoffenberg and colleagues is presented.

The N1 neuraminidases (NAs) of avian and pandemic human influenza viruses contain tyrosine and asparagine, respectively, at position 347 on the rim of the catalytic site; the biological significance of this difference is not clear. Here, we used molecular dynamics simulation to model the effects of amino acid 347 on N1 NA interactions with sialyllacto-N-tetraoses 6'SLN-LC and 3'SLN-LC, which represent NA substrates in humans and birds, respectively. Our analysis predicted that Y347 plays an important role in the NA preference for the avian-type substrates. The Y347N substitution facilitates hydrolysis of human-type substrates by resolving steric conflicts of the Neu5Ac2-6Gal moiety with the bulky side chain of Y347, decreasing the free energy of substrate binding, and increasing the solvation of the Neu5Ac2-6Gal bond. Y347 was conserved in all N1 NA sequences of avian influenza viruses in the GISAID EpiFlu database with two exceptions. First, the Y347F substitution was present in the NA of a specific H6N1 poultry virus lineage and was associated with the substitutions G228S and/or E190V/L in the receptor-binding site (RBS) of the hemagglutinin (HA). Second, the highly pathogenic avian H5N1 viruses of the Gs/Gd lineage contained sporadic variants with the NA substitutions Y347H/D, which were frequently associated with substitutions in the HA RBS. The Y347N substitution occurred following the introductions of avian precursors into humans and pigs with N/D347 conserved during virus circulation in these hosts. Comparative evolutionary analysis of site 347 revealed episodic positive selection across the entire tree and negative selection within most host-specific groups of viruses, suggesting that substitutions at NA position 347 occurred during host switches and remained under pervasive purifying selection thereafter. Our results elucidate the role of amino acid 347 in NA recognition of sialoglycan substrates and emphasize the significance of substitutions at position 347 as a marker of host range and adaptive evolution of influenza viruses., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2023 Elli, Raffaini, Guerrini, Kosakovsky Pond and Matrosovich.)

The programmed frameshift element (PFE) rerouting translation from ORF1a to ORF1b is essential for the propagation of coronaviruses. The combination of genomic features that make up PFE-the overlap between the two reading frames, a slippery sequence, as well as an ensemble of complex secondary structure elements-places severe constraints on this region as most possible nucleotide substitution may disrupt one or more of these elements. The vast amount of SARS-CoV-2 sequencing data generated within the past year provides an opportunity to assess the evolutionary dynamics of PFE in great detail. Here, we performed a comparative analysis of all available coronaviral genomic data available to date. We show that the overlap between ORF1a and ORF1b evolved as a set of discrete 7, 16, 22, 25, and 31 nucleotide stretches with a well-defined phylogenetic specificity. We further examined sequencing data from over 1,500,000 complete genomes and 55,000 raw read data sets to demonstrate exceptional conservation and detect signatures of selection within the PFE region., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The hemagglutinin (HA) of A/H3N2 pandemic influenza viruses (IAVs) of 1968 differed from its inferred avian precursor by eight amino acid substitutions. To determine their phenotypic effects, we studied recombinant variants of A/Hong Kong/1/1968 virus containing either human-type or avian-type amino acids in the corresponding positions of HA. The precursor HA displayed receptor binding profile and high conformational stability typical for duck IAVs. Substitutions Q226L and G228S, in addition to their known effects on receptor specificity and replication, marginally decreased HA stability. Substitutions R62I, D63N, D81N and N193S reduced HA binding avidity. Substitutions R62I, D81N and A144G promoted viral replication in human airway epithelial cultures. Analysis of HA sequences revealed that substitutions D63N and D81N accompanied by the addition of N-glycans represent common markers of avian H3 HA adaptation to mammals. Our results advance understanding of genotypic and phenotypic changes in IAV HA required for avian-to-human adaptation and pandemic emergence., Competing Interests: The authors have declared that no competing interests exist.

The current state of much of the Wuhan pneumonia virus (severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2]) research shows a regrettable lack of data sharing and considerable analytical obfuscation. This impedes global research cooperation, which is essential for tackling public health emergencies and requires unimpeded access to data, analysis tools, and computational infrastructure. Here, we show that community efforts in developing open analytical software tools over the past 10 years, combined with national investments into scientific computational infrastructure, can overcome these deficiencies and provide an accessible platform for tackling global health emergencies in an open and transparent manner. Specifically, we use all SARS-CoV-2 genomic data available in the public domain so far to (1) underscore the importance of access to raw data and (2) demonstrate that existing community efforts in curation and deployment of biomedical software can reliably support rapid, reproducible research during global health crises. All our analyses are fully documented at https://github.com/galaxyproject/SARS-CoV-2., Competing Interests: Dannon Baker, Dan Blankenberg, Nate Coroar, John Chilton, James Taylor, and Anton Nekrutenko are founders of and hold equity in GalaxyWorks, LLC. The results of the study discussed in this publication could affect the value of GalaxyWorks, LLC.

Transmission of multiple founder variants has been associated with faster HIV disease progression. Many studies have attempted to determine the number of founder variants, mainly by analysis of sequence diversity and/or tree topology from acutely HIV-infected individuals. We hypothesized that adding sequence data collected from source partners might improve resolution and characterization of transmission events. Blood plasma samples were collected from both the source and recipient in thirty epidemiologically- and phylogenetically linked transmission pairs. All were men who have sex with men, sampled on average 70 days (range 11-170) after the recipient's estimated date of infection. Next generation sequencing (454 FLX, Roche) of HIV-1 env (C2-V3) was performed for all samples. Inspection of sequence alignments, highlighter plots, phylogenetic tree topologies and sequence diversity were used to determine the multiplicity of founder viruses with and without the inclusion of source data. Using only recipient sequence data, we were able to resolve multiplicity in twenty-six of the thirty transmission pairs (87 percent). Among them, five presented with a high viral diversity at baseline (>0.10 subst/site), consistent with multiple founders. By incorporating sequence data collected from the source partner, we were able to characterize all thirty transmission pairs. Overall, sixteen transmission events (53.3 percent) involved multiple founders. Results obtained by combining sequence data from recipient and source were congruent for nineteen of the twenty-six (73 percent) cases where conclusions were made using only recipient sequence data. The multiplicity of founders was associated with significantly higher HIV RNA levels ( P  = 0.04). To further evaluate the transmission bottleneck, we focused on single founder transmissions (fourteen of the thirty), and identified four recipients (28.6 percent) that had founder variants that were inferred to arise from minority viral populations in the source. These source clades ranged from 1.0 to 5.4 percent of the sampled population. Incorporating sequence data from the source increased of the ability to determine the multiplicity of founder variants, reduced misclassification, and allowed us to infer the transmission of minority variants.

Broadly neutralizing antibodies (bNAbs) PG9 and PG16 were isolated from an International AIDS Vaccine Initiative (IAVI) Protocol G subject infected with human immunodeficiency virus type 1 (HIV-1) clade A. Both antibodies are highly potent and neutralize greater than 70% of viruses tested. We sought to begin immunogen design based on viral sequences from this patient; however, pseudoviruses prepared with 19 envelope sequences from this subject were resistant to neutralization by PG9 and PG16. Therefore, we used a bioinformatics approach to identify closely related viruses that were potentially sensitive to PG9 and PG16. A most-recent common ancestor (MRCA) sequence for the viral envelope (Env) was determined and aligned with 99 subtype A gp160 sequences from the Los Alamos HIV database. Virus BG505.W6M.ENV.C2 (BG505) was found to have the highest degree of homology (73%) to the MRCA sequence. Pseudoviruses prepared with this Env were sensitive to neutralization with a broad panel of bNAbs, including PG9 and PG16. When expressed by 293T cells as soluble gp120, the BG505 monomer bound well to both PG9 and PG16. We further showed that a point mutation (L111A) enabled more efficient production of a stable gp120 monomer that preserves the major neutralization epitopes. Finally, we showed that an adjuvanted formulation of this gp120 protein elicited neutralizing antibodies in rabbits (following a gp120 DNA vaccine prime) and that the antisera competed with bNAbs from 3 classes of nonoverlapping epitopes. Thus, the BG505 Env protein warrants further investigation as an HIV vaccine candidate, as a stand-alone protein, or as a component of a vaccine vector.

Markov models of codon substitution are powerful inferential tools for studying biological processes such as natural selection and preferences in amino acid substitution. The equilibrium character distributions of these models are almost always estimated using nucleotide frequencies observed in a sequence alignment, primarily as a matter of historical convention. In this note, we demonstrate that a popular class of such estimators are biased, and that this bias has an adverse effect on goodness of fit and estimates of substitution rates. We propose a "corrected" empirical estimator that begins with observed nucleotide counts, but accounts for the nucleotide composition of stop codons. We show via simulation that the corrected estimates outperform the de facto standard estimates not just by providing better estimates of the frequencies themselves, but also by leading to improved estimation of other parameters in the evolutionary models. On a curated collection of sequence alignments, our estimators show a significant improvement in goodness of fit compared to the approach. Maximum likelihood estimation of the frequency parameters appears to be warranted in many cases, albeit at a greater computational cost. Our results demonstrate that there is little justification, either statistical or computational, for continued use of the -style estimators.

The single rate codon model of non-synonymous substitution is ubiquitous in phylogenetic modeling. Indeed, the use of a non-synonymous to synonymous substitution rate ratio parameter has facilitated the interpretation of selection pressure on genomes. Although the single rate model has achieved wide acceptance, we argue that the assumption of a single rate of non-synonymous substitution is biologically unreasonable, given observed differences in substitution rates evident from empirical amino acid models. Some have attempted to incorporate amino acid substitution biases into models of codon evolution and have shown improved model performance versus the single rate model. Here, we show that the single rate model of non-synonymous substitution is easily outperformed by a model with multiple non-synonymous rate classes, yet in which amino acid substitution pairs are assigned randomly to these classes. We argue that, since the single rate model is so easy to improve upon, new codon models should not be validated entirely on the basis of improved model fit over this model. Rather, we should strive to both improve on the single rate model and to approximate the general time-reversible model of codon substitution, with as few parameters as possible, so as to reduce model over-fitting. We hint at how this can be achieved with a Genetic Algorithm approach in which rate classes are assigned on the basis of sequence information content.

How many species inhabit our immediate surroundings? A straightforward collection technique suitable for answering this question is known to anyone who has ever driven a car at highway speeds. The windshield of a moving vehicle is subjected to numerous insect strikes and can be used as a collection device for representative sampling. Unfortunately the analysis of biological material collected in that manner, as with most metagenomic studies, proves to be rather demanding due to the large number of required tools and considerable computational infrastructure. In this study, we use organic matter collected by a moving vehicle to design and test a comprehensive pipeline for phylogenetic profiling of metagenomic samples that includes all steps from processing and quality control of data generated by next-generation sequencing technologies to statistical analyses and data visualization. To the best of our knowledge, this is also the first publication that features a live online supplement providing access to exact analyses and workflows used in the article.