Analysis of Immune Escape Variants from Antibody-Based Therapeutics against COVID-19: A Systematic Review

The accelerated SARS-CoV-2 evolution under selective pressure by massive deployment of neutralizing antibody-based therapeutics is a concern with potentially severe implications for public health. We review here reports of documented immune escape after treatment with monoclonal antibodies and COVID-19-convalescent plasma (CCP). While the former is mainly associated with specific single amino acid mutations at residues within the receptor-binding domain (e.g., E484K/Q, Q493R, and S494P), a few cases of immune evasion after CCP were associated with recurrent deletions within the N-terminal domain of the spike protein (e.g., ΔHV69-70, ∆LGVY141-144 and ΔAL243-244). The continuous genomic monitoring of non-responders is needed to better understand immune escape frequencies and the fitness of emerging variants.

Human TRMT112-Methyltransferase Network Consists of Seven Partners Interacting with a Common Co-Factor

Methylation is an essential epigenetic modification mainly catalysed by S-Adenosyl methionine-dependent methyltransferases (MTases). Several MTases require a cofactor for their metabolic stability and enzymatic activity. TRMT112 is a small evolutionary conserved protein that acts as a co-factor and activator for different MTases involved in rRNA, tRNA and protein methylation. Using a SILAC screen, we pulled down seven methyltransferases—N6AMT1, WBSCR22, METTL5, ALKBH8, THUMPD2, THUMPD3 and TRMT11—as interaction partners of TRMT112. We showed that TRMT112 stabilises all seven MTases in cells. TRMT112 and MTases exhibit a strong mutual feedback loop when expressed together in cells. TRMT112 interacts with its partners in a similar way; however, single amino acid mutations on the surface of TRMT112 reveal several differences as well. In summary, mammalian TRMT112 can be considered as a central “hub” protein that regulates the activity of at least seven methyltransferases.

Omicron variant of SARS-CoV-2 harbors a unique insertion mutation of putative viral or human genomic origin

The emergence of a heavily mutated SARS-CoV-2 variant (B.1.1.529, Omicron) and it’s spread to 6 continents within a week of initial discovery has set off a global public health alarm. Characterizing the mutational profile of Omicron is necessary to interpret its shared or distinctive clinical phenotypes with other SARS-CoV-2 variants. We compared the mutations of Omicron with prior variants of concern (Alpha, Beta, Gamma, Delta), variants of interest (Lambda, Mu, Eta, Iota and Kappa), and all 1523 SARS-CoV-2 lineages constituting 5.4 million SARS-CoV-2 genomes. Omicron’s Spike protein has 26 amino acid mutations (23 substitutions, two deletions and one insertion) that are distinct compared to other variants of concern. Whereas the substitution and deletion mutations have appeared in previous SARS-CoV-2 lineages, the insertion mutation (ins214EPE) has not been previously observed in any SARS-CoV-2 lineage other than Omicron. The nucleotide sequence encoding for ins214EPE could have been acquired by template switching involving the genomes of other viruses that infect the same host cells as SARS-CoV-2 or the human transcriptome of host cells infected with SARS-CoV-2. For instance, given recent clinical reports of co-infections in COVID-19 patients with seasonal coronaviruses (e.g. HCoV-229E), single cell RNA-sequencing data showing co-expression of the SARS-CoV-2 and HCoV-229E entry receptors (ACE2 and ANPEP) in respiratory and gastrointestinal cells, and HCoV genomes harboring sequences homologous to the nucleotide sequence that encodes ins214EPE, it is plausible that the Omicron insertion could have evolved in a co-infected individual. There is a need to understand the function of the Omicron insertion and whether human host cells are being exploited by SARS-CoV-2 as an ‘evolutionary sandbox’ for host-virus and inter-viral genomic interplay.

A global analysis of conservative and non-conservative mutations in SARS-CoV-2 detected in the first year of the COVID-19 world-wide diffusion

The ability of SARS-CoV-2 to rapidly mutate represents a remarkable complicancy. Quantitative evaluations of the effects that these mutations have on the virus structure/function is of great relevance and the availability of a large number of SARS-CoV-2 sequences since the early phases of the pandemic represents a unique opportunity to follow the adaptation of the virus to humans. Here, we evaluated the SARS-CoV-2 amino acid mutations and their progression by analyzing publicly available viral genomes at three stages of the pandemic (2020 March 15th and October 7th, 2021 February 7th). Mutations were classified in conservative and non-conservative based on the probability to be accepted during the evolution according to the Point Accepted Mutation substitution matrices and on the similarity of the encoding codons. We found that the most frequent substitutions are T > I, L > F, and A > V and we observe accumulation of hydrophobic residues. These findings are consistent among the three stages analyzed. We also found that non-conservative mutations are less frequent than conservative ones. This finding may be ascribed to a progressive adaptation of the virus to the host. In conclusion, the present study provides indications of the early evolution of the virus and tools for the global and genome-specific evaluation of the possible impact of mutations on the structure/function of SARS-CoV-2 variants.

The Global Epidemic of the SARS-CoV-2 Delta Variant, Key Spike Mutations and Immune Escape

During the COVID-19 pandemic, SARS-CoV-2 variants have emerged and spread worldwide. The Delta (B.1.617.2) variant was first reported in India in October 2020 and was classified as a “variant of concern (VOC)” by the WHO on 11 May, 2021. Compared to the wild-type strain, several studies have shown that the Delta variant is more transmissible and has higher viral loads in infected samples. COVID-19 patients infected with the Delta variant have a higher risk of hospitalization, intensive care unit (ICU) admission, and mortality. The Delta variant is becoming the dominant strain in many countries around the world. This review summarizes and analyses the biological characteristics of key amino acid mutations, the epidemic characteristics, and the immune escape of the Delta variant. We hope to provide scientific reference for the monitoring and prevention measures of the SARS-CoV-2 Delta variant and the development strategy of a second-generation vaccine.

A Routine Sanger Sequencing Target Specific Mutation Assay for SARS-CoV-2 Variants of Concern and Interest

As SARS-CoV-2 continues to spread among human populations, genetic changes occur and accumulate in the circulating virus. Some of these genetic changes have caused amino acid mutations, including deletions, which may have a potential impact on critical SARS-CoV-2 countermeasures, including vaccines, therapeutics, and diagnostics. Considerable efforts have been made to categorize the amino acid mutations of the angiotensin-converting enzyme 2 (ACE2) receptor binding domain (RBD) of the spike (S) protein, along with certain mutations in other regions within the S protein as specific variants, in an attempt to study the relationship between these mutations and the biological behavior of the virus. However, the currently used whole genome sequencing surveillance technologies can test only a small fraction of the positive specimens with high viral loads and often generate uncertainties in nucleic acid sequencing that needs additional verification for precision determination of mutations. This article introduces a generic protocol to routinely sequence a 437-bp nested RT-PCR cDNA amplicon of the ACE2 RBD and a 490-bp nested RT-PCR cDNA amplicon of the N-terminal domain (NTD) of the S gene for detection of the amino acid mutations needed for accurate determination of all variants of concern and variants of interest according to the definitions published by the U.S. Centers for Disease Control and Prevention. This protocol was able to amplify both nucleic acid targets into cDNA amplicons to be used as templates for Sanger sequencing on all 16 clinical specimens that were positive for SARS-CoV-2.

Simulated vector transmission differentially influences dynamics of two viral variants of deformed wing virus in honey bees (Apis mellifera)

Understanding how vectors alter the interactions between viruses and their hosts is a fundamental question in virology and disease ecology. In honey bees, transmission of deformed wing virus (DWV) by parasitic Varroa mites has been associated with elevated disease and host mortality, and Varroa transmission has been hypothesized to lead to increased viral titres or select for more virulent variants. Here, we mimicked Varroa transmission by serially passaging a mixed population of two DWV variants, A and B, by injection through in vitro reared honey bee pupae and tracking these viral populations through five passages. The DWV-A and DWV-B variant proportions shifted dynamically through passaging, with DWV-B outcompeting DWV-A after one passage, but levels of both variants becoming equivalent by Passage 5. Sequencing analysis revealed a dominant, recombinant DWV-B strain (DWV-A derived 5′ IRES region with the rest of the genome DWV-B), with low nucleotide diversity that decreased through passaging. DWV-A populations had higher nucleotide diversity compared to DWV-B, but this also decreased through passaging. Selection signatures were found across functional regions of the DWV-A and DWV-B genomes, including amino acid mutations in the putative capsid protein region. Simulated vector transmission differentially impacted two closely related viral variants which could influence viral interactions with the host, demonstrating surprising plasticity in vector-host-viral dynamics.

Identification of Key Interactions of Benzimidazole Resistance-Associated Amino Acid Mutations in Ascaris β-Tubulins By Molecular Docking Simulations

Ascaris species are soil-transmitted helminths that infect humans and livestock mainly in low and middle-income countries. Benzimidazole (BZ) class drugs have predominated for many years in the treatment of Ascaris infections, but persistent use of BZs has already led to widespread resistance in other nematodes, and treatment failure is emerging for Ascaris. Benzimidazoles act by binding to β-tubulin proteins and destabilising microtubules. Three mutations in the β-tubulin protein family are associated with BZ resistance. Seven shared β-tubulin isotypes were identified in Ascaris lumbricoides and A. suum genomes. Benzimidazoles were predicted to bind to all β-tubulin isotypes using in silico docking, demonstrating that the selectivity of BZs to interact with one or two β-tubulin isotypes is likely the result of isotype expression levels affecting the frequency of interaction. Ascaris β-tubulin isotype A clusters with helminth β-tubulins previously shown to interact with BZ. Molecular dynamics simulations using β-tubulin isotype A highlighted the key role of amino acid E198 in BZ-β-tubulin interactions. Simulations indicated that mutations at amino acids E198A and F200Y alter binding of BZ, whereas there was no obvious effect of the F167Y mutation. In conclusion, the key interactions vital for BZ binding with β-tubulins have been identified and show how mutations can lead to resistance in nematodes.

De novo assembly and functional annotation of blood transcriptome of loggerhead turtle, and in silico characterization of peroxiredoxins and thioredoxins

The aim of this study was to generate and analyze the atlas of the loggerhead turtle blood transcriptome by RNA-seq, as well as identify and characterize thioredoxin (Tnxs) and peroxiredoxin (Prdxs) antioxidant enzymes of the greatest interest in the control of peroxide levels and other biological functions. The transcriptome of loggerhead turtle was sequenced using the Illumina Hiseq 2000 platform and de novo assembly was performed using the Trinity pipeline. The assembly comprised 515,597 contigs with an N50 of 2,631 bp. Contigs were analyzed with CD-Hit obtaining 374,545 unigenes, of which 165,676 had ORFs encoding putative proteins longer than 100 amino acids. A total of 52,147 (31.5%) of these transcripts had significant homology matches in at least one of the five databases used. From the enrichment of GO terms, 180 proteins with antioxidant activity were identified, among these 28 Prdxs and 50 putative Tnxs. The putative proteins of loggerhead turtles encoded by the genes Prdx1, Prdx3, Prdx5, Prdx6, Txn and Txnip were predicted and characterized in silico. When comparing Prdxs and Txns of loggerhead turtle with homologous human proteins, they showed 18 (9%), 52 (18%) 94 (43%), 36 (16%), 35 (33%) and 74 (19%) amino acid mutations respectively. However, they showed high conservation in active sites and structural motifs (98%), with few specific modifications. Of these, Prdx1, Prdx3, Prdx5, Prdx6, Txn and Txnip presented 0, 25, 18, three, six and two deleterious changes. This study provides a high quality blood transcriptome and functional annotation of loggerhead sea turtles.

ANALYSIS OF IMMUNE ESCAPE VARIANTS FROM ANTIBODY-BASED THERAPEUTICS AGAINST COVID-19.

Accelerated SARS-CoV-2 evolution under selective pressure by massive deployment of neutralizing antibody-based therapeutics is a concern with potentially severe implications for public health. We review here reports of documented immune escape after treatment with monoclonal antibodies and COVID19 convalescent plasma (CCP). While the former is mainly associated with specific single amino acid mutations at residues within the receptor-binding domain (e.g., E484K/Q, Q493R, and S494P), the few cases of immune evasion after CCP were associated with recurrent deletions within the N-terminal domain of Spike protein (e.g, delHV69-70, delLGVY141-144 and delAL243-244). Continuous genomic monitoring of non-responders is needed to better understand immune escape frequencies and fitness of emerging variants.