Stability of SARS-CoV-2 phylogenies.

Funder: Alfred P. Sloan Foundation; funder-id: http://dx.doi.org/10.13039/100000879Funder: European Molecular Biology Laboratory (EMBL)The SARS-CoV-2 pandemic has led to unprecedented, nearly real-time genetic tracing due to the rapid community sequencing response. Researchers immediately leveraged these data to infer the evolutionary relationships among viral samples and to study key biological questions, including whether host viral genome editing and recombination are features of SARS-CoV-2 evolution. This global sequencing effort is inherently decentralized and must rely on data collected by many labs using a wide variety of molecular and bioinformatic techniques. There is thus a strong possibility that systematic errors associated with lab-or protocol-specific practices affect some sequences in the repositories. We find that some recurrent mutations in reported SARS-CoV-2 genome sequences have been observed predominantly or exclusively by single labs, co-localize with commonly used primer binding sites and are more likely to affect the protein-coding sequences than other similarly recurrent mutations. We show that their inclusion can affect phylogenetic inference on scales relevant to local lineage tracing, and make it appear as though there has been an excess of recurrent mutation or recombination among viral lineages. We suggest how samples can be screened and problematic variants removed, and we plan to regularly inform the scientific community with our updated results as more SARS-CoV-2 genome sequences are shared (https://virological.org/t/issues-with-sars-cov-2-sequencing-data/473 and https://virological.org/t/masking-strategies-for-sars-cov-2-alignments/480). We also develop tools for comparing and visualizing differences among very large phylogenies and we show that consistent clade- and tree-based comparisons can be made between phylogenies produced by different groups. These will facilitate evolutionary inferences and comparisons among phylogenies produced for a wide array of purposes. Building on the SARS-CoV-2 Genome Browser at UCSC, we present a toolkit to compare, analyze and combine SARS-CoV-2 phylogenies, find and remove potential sequencing errors and establish a widely shared, stable clade structure for a more accurate scientific inference and discourse

Stability of SARS-CoV-2 Phylogenies

AbstractThe SARS-CoV-2 pandemic has led to unprecedented, nearly real-time genetic tracing due to the rapid community sequencing response. Researchers immediately leveraged these data to infer the evolutionary relationships among viral samples and to study key biological questions, including whether host viral genome editing and recombination are features of SARS-CoV-2 evolution. This global sequencing effort is inherently decentralized and must rely on data collected by many labs using a wide variety of molecular and bioinformatic techniques. There is thus a strong possibility that systematic errors associated with lab-specific practices affect some sequences in the repositories. We find that some recurrent mutations in reported SARS-CoV-2 genome sequences have been observed predominantly or exclusively by single labs, co-localize with commonly used primer binding sites and are more likely to affect the protein coding sequences than other similarly recurrent mutations. We show that their inclusion can affect phylogenetic inference on scales relevant to local lineage tracing, and make it appear as though there has been an excess of recurrent mutation and/or recombination among viral lineages. We suggest how samples can be screened and problematic mutations removed. We also develop tools for comparing and visualizing differences among phylogenies and we show that consistent clade- and tree-based comparisons can be made between phylogenies produced by different groups. These will facilitate evolutionary inferences and comparisons among phylogenies produced for a wide array of purposes. Building on the SARS-CoV-2 Genome Browser at UCSC, we present a toolkit to compare, analyze and combine SARS-CoV-2 phylogenies, find and remove potential sequencing errors and establish a widely shared, stable clade structure for a more accurate scientific inference and discourse.ForewordWe wish to thank all groups that responded rapidly by producing these invaluable and essential sequence data. Their contributions have enabled an unprecedented, lightning-fast process of scientific discovery---truly an incredible benefit for humanity and for the scientific community. We emphasize that most lab groups with whom we associate specific suspicious alleles are also those who have produced the most sequence data at a time when it was urgently needed. We commend their efforts. We have already contacted each group and many have updated their sequences. Our goal with this work is not to highlight potential errors, but to understand the impacts of these and other kinds of highly recurrent mutations so as to identify commonalities among the suspicious examples that can improve sequence quality and analysis going forward.</jats:sec

Predicting transfer RNA gene activity from sequence and genome context

Transfer RNA (tRNA) genes are among the most highly transcribed genes in the genome owing to their central role in protein synthesis. However, there is evidence for a broad range of gene expression across tRNA loci. This complexity, combined with difficulty in measuring transcript abundance and high sequence identity across transcripts, has severely limited our collective understanding of tRNA gene expression regulation and evolution. We establish sequence-based correlates to tRNA gene expression and develop a tRNA gene classification method that does not require, but benefits from, comparative genomic information and achieves accuracy comparable to molecular assays. We observe that guanine + cytosine (G + C) content and CpG density surrounding tRNA loci is exceptionally well correlated with tRNA gene activity, supporting a prominent regulatory role of the local genomic context in combination with internal sequence features. We use our tRNA gene activity predictions in conjunction with a comprehensive tRNA gene ortholog set spanning 29 placental mammals to estimate the evolutionary rate of functional changes among orthologs. Our method adds a new dimension to large-scale tRNA functional prediction and will help prioritize characterization of functional tRNA variants. Its simplicity and robustness should enable development of similar approaches for other clades, as well as exploration of functional diversification of members of large gene families.</jats:p

Altering the Immunogenicity of Hemagglutinin Immunogens by Hyperglycosylation and Disulfide Stabilization

Influenza virus alters glycosylation patterns on its surface exposed glycoproteins to evade host adaptive immune responses. The viral hemagglutinin (HA), in particular the H3 subtype, has increased its overall surface glycosylation since its introduction in 1968. We previously showed that modulating predicted N-linked glycosylation sites on H3 A/Hong Kong/1/1968 HA identified a conserved epitope at the HA interface. This epitope is occluded on the native HA trimer but is likely exposed during HA “breathing” on the virion surface. Antibodies directed to this site are protective via an ADCC-mediated mechanism. This glycan engineering strategy made an otherwise subdominant epitope dominant in the murine model. Here, we asked whether cysteine stabilization of the hyperglycosylated HA trimer could reverse this immunodominance by preventing access to the interface epitope and focus responses to the HA receptor binding site (RBS). While analysis of serum responses from immunized mice did not show a redirection to the RBS, cysteine stabilization did result in an overall reduction in immunogenicity of the interface epitope. Thus, glycan engineering and cysteine stabilization are two strategies that can be used together to alter immunodominance patterns to HA. These results add to rational immunogen design approaches used to manipulate immune responses for the development of next-generation influenza vaccines.</jats:p

Image_1_Altering the Immunogenicity of Hemagglutinin Immunogens by Hyperglycosylation and Disulfide Stabilization.tif

Influenza virus alters glycosylation patterns on its surface exposed glycoproteins to evade host adaptive immune responses. The viral hemagglutinin (HA), in particular the H3 subtype, has increased its overall surface glycosylation since its introduction in 1968. We previously showed that modulating predicted N-linked glycosylation sites on H3 A/Hong Kong/1/1968 HA identified a conserved epitope at the HA interface. This epitope is occluded on the native HA trimer but is likely exposed during HA “breathing” on the virion surface. Antibodies directed to this site are protective via an ADCC-mediated mechanism. This glycan engineering strategy made an otherwise subdominant epitope dominant in the murine model. Here, we asked whether cysteine stabilization of the hyperglycosylated HA trimer could reverse this immunodominance by preventing access to the interface epitope and focus responses to the HA receptor binding site (RBS). While analysis of serum responses from immunized mice did not show a redirection to the RBS, cysteine stabilization did result in an overall reduction in immunogenicity of the interface epitope. Thus, glycan engineering and cysteine stabilization are two strategies that can be used together to alter immunodominance patterns to HA. These results add to rational immunogen design approaches used to manipulate immune responses for the development of next-generation influenza vaccines.</p

Image_3_Altering the Immunogenicity of Hemagglutinin Immunogens by Hyperglycosylation and Disulfide Stabilization.tif

Influenza virus alters glycosylation patterns on its surface exposed glycoproteins to evade host adaptive immune responses. The viral hemagglutinin (HA), in particular the H3 subtype, has increased its overall surface glycosylation since its introduction in 1968. We previously showed that modulating predicted N-linked glycosylation sites on H3 A/Hong Kong/1/1968 HA identified a conserved epitope at the HA interface. This epitope is occluded on the native HA trimer but is likely exposed during HA “breathing” on the virion surface. Antibodies directed to this site are protective via an ADCC-mediated mechanism. This glycan engineering strategy made an otherwise subdominant epitope dominant in the murine model. Here, we asked whether cysteine stabilization of the hyperglycosylated HA trimer could reverse this immunodominance by preventing access to the interface epitope and focus responses to the HA receptor binding site (RBS). While analysis of serum responses from immunized mice did not show a redirection to the RBS, cysteine stabilization did result in an overall reduction in immunogenicity of the interface epitope. Thus, glycan engineering and cysteine stabilization are two strategies that can be used together to alter immunodominance patterns to HA. These results add to rational immunogen design approaches used to manipulate immune responses for the development of next-generation influenza vaccines.</p

Image_2_Altering the Immunogenicity of Hemagglutinin Immunogens by Hyperglycosylation and Disulfide Stabilization.tif

Influenza virus alters glycosylation patterns on its surface exposed glycoproteins to evade host adaptive immune responses. The viral hemagglutinin (HA), in particular the H3 subtype, has increased its overall surface glycosylation since its introduction in 1968. We previously showed that modulating predicted N-linked glycosylation sites on H3 A/Hong Kong/1/1968 HA identified a conserved epitope at the HA interface. This epitope is occluded on the native HA trimer but is likely exposed during HA “breathing” on the virion surface. Antibodies directed to this site are protective via an ADCC-mediated mechanism. This glycan engineering strategy made an otherwise subdominant epitope dominant in the murine model. Here, we asked whether cysteine stabilization of the hyperglycosylated HA trimer could reverse this immunodominance by preventing access to the interface epitope and focus responses to the HA receptor binding site (RBS). While analysis of serum responses from immunized mice did not show a redirection to the RBS, cysteine stabilization did result in an overall reduction in immunogenicity of the interface epitope. Thus, glycan engineering and cysteine stabilization are two strategies that can be used together to alter immunodominance patterns to HA. These results add to rational immunogen design approaches used to manipulate immune responses for the development of next-generation influenza vaccines.</p

Image_4_Altering the Immunogenicity of Hemagglutinin Immunogens by Hyperglycosylation and Disulfide Stabilization.tif

Influenza virus alters glycosylation patterns on its surface exposed glycoproteins to evade host adaptive immune responses. The viral hemagglutinin (HA), in particular the H3 subtype, has increased its overall surface glycosylation since its introduction in 1968. We previously showed that modulating predicted N-linked glycosylation sites on H3 A/Hong Kong/1/1968 HA identified a conserved epitope at the HA interface. This epitope is occluded on the native HA trimer but is likely exposed during HA “breathing” on the virion surface. Antibodies directed to this site are protective via an ADCC-mediated mechanism. This glycan engineering strategy made an otherwise subdominant epitope dominant in the murine model. Here, we asked whether cysteine stabilization of the hyperglycosylated HA trimer could reverse this immunodominance by preventing access to the interface epitope and focus responses to the HA receptor binding site (RBS). While analysis of serum responses from immunized mice did not show a redirection to the RBS, cysteine stabilization did result in an overall reduction in immunogenicity of the interface epitope. Thus, glycan engineering and cysteine stabilization are two strategies that can be used together to alter immunodominance patterns to HA. These results add to rational immunogen design approaches used to manipulate immune responses for the development of next-generation influenza vaccines.</p

Stability of SARS-CoV-2 phylogenies.

The SARS-CoV-2 pandemic has led to unprecedented, nearly real-time genetic tracing due to the rapid community sequencing response. Researchers immediately leveraged these data to infer the evolutionary relationships among viral samples and to study key biological questions, including whether host viral genome editing and recombination are features of SARS-CoV-2 evolution. This global sequencing effort is inherently decentralized and must rely on data collected by many labs using a wide variety of molecular and bioinformatic techniques. There is thus a strong possibility that systematic errors associated with lab-or protocol-specific practices affect some sequences in the repositories. We find that some recurrent mutations in reported SARS-CoV-2 genome sequences have been observed predominantly or exclusively by single labs, co-localize with commonly used primer binding sites and are more likely to affect the protein-coding sequences than other similarly recurrent mutations. We show that their inclusion can affect phylogenetic inference on scales relevant to local lineage tracing, and make it appear as though there has been an excess of recurrent mutation or recombination among viral lineages. We suggest how samples can be screened and problematic variants removed, and we plan to regularly inform the scientific community with our updated results as more SARS-CoV-2 genome sequences are shared (https://virological.org/t/issues-with-sars-cov-2-sequencing-data/473 and https://virological.org/t/masking-strategies-for-sars-cov-2-alignments/480). We also develop tools for comparing and visualizing differences among very large phylogenies and we show that consistent clade- and tree-based comparisons can be made between phylogenies produced by different groups. These will facilitate evolutionary inferences and comparisons among phylogenies produced for a wide array of purposes. Building on the SARS-CoV-2 Genome Browser at UCSC, we present a toolkit to compare, analyze and combine SARS-CoV-2 phylogenies, find and remove potential sequencing errors and establish a widely shared, stable clade structure for a more accurate scientific inference and discourse