Genomic surveillance of SARS-CoV-2 in Thailand reveals mixed imported populations, a local lineage expansion and a virus with truncated ORF7a

Near complete genome sequences from Southern Vietnam revealed local features of genetic diversity and intergenerational changes in SARS- CoV-2 variants in 2020-2021

BACKGROUND

Since its beginnings in 2019, the COVID-19 pandemic is still a problem of global medical concern. Southern Vietnam is one of the country's vast regions, including 20 provinces and the densely populated metropolis Ho Chi Minh City. A randomized retrospective study was performed to investigate the epidemiology and genetic diversity of COVID-19. Whole-genome sequencing of 126 SARS-CoV-2 samples collected from Southern Vietnam between January 2020 and December 2021 revealed the main circulating variants and their distribution.

METHODS

Epidemiological data were obtained from the Department of Preventive Medicine of the Vietnamese Ministry of Health. To identify circulating variants, RNA, extracted from 126 nasopharyngeal swabs of patients with suspected COVID-19 were sequenced on Illunina MiSeq to obtain near complete genomes SARS-CoV-2.

RESULTS

Due to the effectiveness of restrictive measures in Vietnam, it was possible to keep incidence at a low level. The partial relaxation of restrictive measures, and the spread of Delta lineages, contributed to the beginning of a logarithmic increase in incidence. Lineages 20A-H circulated in Southern Vietnam during 2020. Spread of the Delta lineage in Southern Vietnam began in March 2021, causing a logarithmic rise in the number of COVID-19 cases.

CONCLUSIONS

Pandemic dynamics in Southern Vietnam feature specific variations in incidence, and these reflect the success of the restrictive measures put in place during the early stages of the pandemic.

Evaluation of SARS-CoV-2 ORF7a Deletions from COVID-19-Positive Individuals and Its Impact on Virus Spread in Cell Culture

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing the COVID-19 outbreak, posed a primary concern of public health worldwide. The most common changes in SARS-CoV-2 are single nucleotide substitutions, also reported insertions and deletions. This work investigates the presence of SARS-CoV-2 ORF7a deletions identified in COVID-19-positive individuals. Sequencing of SARS-CoV-2 complete genomes showed three different ORF7a size deletions (190-nt, 339-nt and 365-nt). Deletions were confirmed through Sanger sequencing. The ORF7a∆190 was detected in a group of five relatives with mild symptoms of COVID-19, and the ORF7a∆339 and ORF7a∆365 in a couple of co-workers. These deletions did not affect subgenomic RNAs (sgRNA) production downstream of ORF7a. Still, fragments associated with sgRNA of genes upstream of ORF7a showed a decrease in size when corresponding to samples with deletions. In silico analysis suggests that the deletions impair protein proper function; however, isolated viruses with partial deletion of ORF7a can replicate in culture cells similarly to wild-type viruses at 24 hpi, but with less infectious particles after 48 hpi. These findings on deleted ORF7a accessory protein gene, contribute to understanding SARS-CoV-2 phenotypes such as replication, immune evasion and evolutionary fitness as well insights into the role of SARS-CoV-2_ORF7a in the mechanism of virus-host interactions.

Persistence of a Frameshifting Deletion in SARS-CoV-2 ORF7a for the Duration of a Major Outbreak

Australia experienced widespread COVID-19 outbreaks from infection with the SARS-CoV-2 Delta variant between June 2021 and February 2022. A 17-nucleotide frameshift-inducing deletion in ORF7a rapidly became represented at the consensus level (Delta-ORF7a^Δ17del) in most Australian outbreak cases. Studies from early in the COVID-19 pandemic suggest that frameshift-inducing deletions in ORF7a do not persist for long in the population; therefore, Delta-ORF7a^Δ17del genomes should have disappeared early in the Australian outbreak. In this study, we conducted a retrospective analysis of global Delta genomes to characterise the dynamics of Delta-ORF7a^Δ17del over time, determined the frequency of all ORF7a deletions worldwide, and compared global trends with those of the Australian Delta outbreak. We downloaded all GISAID clade GK Delta genomes and scanned them for deletions in ORF7a. For each deletion we identified, we characterised its frequency, the number of countries it was found in, and how long it persisted. Of the 4,018,216 Delta genomes identified globally, 134,751 (~3.35%) possessed an ORF7a deletion, and ORF7a^Δ17del was the most common. ORF7a^Δ17del was the sole deletion in 28,014 genomes, of which 27,912 (~99.6%) originated from the Australian outbreak. During the outbreak, ~87% of genomes were Delta-ORF7a^Δ17del, and genomes with this deletion were sampled until the outbreak's end. These data demonstrate that, contrary to suggestions early in the COVID-19 pandemic, genomes with frameshifting deletions in ORF7a can persist over long time periods. We suggest that the proliferation of Delta-ORF7a^Δ17del genomes was likely a chance founder effect. Nonetheless, the frequency of ORF7a deletions in SARS-CoV-2 genomes worldwide suggests they might have some benefit for virus transmission.

Emergence and spreading of the largest SARS-CoV-2 deletion in the Delta AY.20 lineage from Uruguay

The genetic variability of SARS-CoV-2 (genus Betacoronavirus, family Coronaviridae) has been scrutinized since its first detection in December 2019. Although the role of structural variants, particularly deletions, in virus evolution is little explored, these genome changes are extremely frequent. They are associated with relevant processes, including immune escape and attenuation. Deletions commonly occur in accessory ORFs and might even lead to the complete loss of one or more ORFs. This scenario poses an interesting question about the origin and spreading of extreme structural rearrangements that persist without compromising virus viability. Here, we analyze the genome of SARS-CoV-2 in late 2021 in Uruguay and identify a Delta lineage (AY.20) that experienced a large deletion (872 nucleotides according to the reference Wuhan strain) that removes the 7a, 7b, and 8 ORFs. Deleted viruses coexist with wild-type (without deletion) AY.20 and AY.43 strains. The Uruguayan deletion is like those identified in Delta strains from Poland and Japan but occurs in a different Delta clade. Besides providing proof of the circulation of this large deletion in America, we infer that the 872-deletion arises by the consecutive occurrence of a 6-nucleotide deletion, characteristic of delta strains, and an 866-nucleotide deletion that arose independently in the AY.20 Uruguayan lineage. The largest deletion occurs adjacent to transcription regulatory sequences needed to synthesize the nested set of subgenomic mRNAs that serve as templates for transcription. Our findings support the role of transcription sequences as a hotspot for copy-choice recombination and highlight the remarkable dynamic of SARS-CoV-2 genomes.

Impairment of antiviral immune response and disruption of cellular functions by SARS-CoV-2 ORF7a and ORF7b

SARS-CoV-2, the causative agent of the present COVID-19 pandemic, possesses eleven accessory proteins encoded in its genome, and some have been implicated in facilitating infection and pathogenesis through their interaction with cellular components. Among these proteins, accessory protein ORF7a and ORF7b functions are poorly understood. In this study, A549 cells were transduced to express ORF7a and ORF7b, respectively, to explore more in depth the role of each accessory protein in the pathological manifestation leading to COVID-19. Bioinformatic analysis and integration of transcriptome results identified defined canonical pathways and functional groupings revealing that after expression of ORF7a or ORF7b, the lung cells are potentially altered to create conditions more favorable for SARS-CoV-2, by inhibiting the IFN-I response, increasing proinflammatory cytokines release, and altering cell metabolic activity and adhesion. Based on these results, it is plausible to suggest that ORF7a or ORF7b could be used as biomarkers of progression in this pandemic.

Mutational induction in SARS-CoV-2 major lineages by experimental exposure to neutralising sera

The ongoing evolution of SARS-CoV-2 and the emergence of new viral variants bearing specific escape mutations responsible for immune evasion from antibody neutralisation has required a more accurate characterisation of the immune response as one of the evolutive forces behind viral adaptation to a largely immunised human population. In this work, culturing in the presence of neutralising sera vigorously promoted mutagenesis leading to the acquisition of known escape mutations on the spike as well as new presumptive escape mutations on structural proteins whose role as target of the neutralizing antibody response might have been thus far widely neglected. From this perspective, this study, in addition to tracing the past evolution of the species back to interactions with neutralising antibody immune response, also offers a glimpse into future evolutive scenarios.

Genomic epidemiology of Delta SARS-CoV-2 during transition from elimination to suppression in Aotearoa New Zealand

New Zealand’s COVID-19 elimination strategy heavily relied on the use of genomics to inform contact tracing, linking cases to the border and to clusters during community outbreaks. In August 2021, New Zealand entered its second nationwide lockdown after the detection of a single community case with no immediately apparent epidemiological link to the border. This incursion resulted in the largest outbreak seen in New Zealand caused by the Delta Variant of Concern. Here we generated 3806 high quality SARS-CoV-2 genomes from cases reported in New Zealand between 17 August and 1 December 2021, representing 43% of reported cases. We detected wide geographical spread coupled with undetected community transmission, characterised by the apparent extinction and reappearance of genomically linked clusters. We also identified the emergence, and near replacement, of genomes possessing a 10-nucleotide frameshift deletion that caused the likely truncation of accessory protein ORF7a. By early October, New Zealand moved from an elimination strategy to a suppression strategy and the role of genomics changed markedly from being used to track and trace, towards population-level surveillance.

Aotearoa New Zealand pursued a COVID-19 elimination strategy until October 2021 when it moved to a suppression strategy. In this genomic surveillance study, the authors describe spread of the virus during the transition between these strategies, with evidence of substantial undetected community transmission.

Role of genomics in combating COVID-19 pandemic

The coronavirus disease 2019 (COVID-19) quickly swept over the world, becoming one of the most devastating outbreaks in human history. Being the first pandemic in the post-genomic era, advancements in genomics contributed significantly to scientific understanding and public health response to COVID-19. Genomic technologies have been employed by researchers all over the world to better understand the biology of SARS-CoV-2 and its origin, genomic diversity, and evolution. Worldwide genomic resources have greatly aided in the investigation of the COVID-19 pandemic. The pandemic has ushered in a new era of genomic surveillance, wherein scientists are tracking the changes of the SARS-CoV-2 genome in real-time at the international and national levels. Availability of genomic and proteomic information enables the rapid development of molecular diagnostics and therapeutics. The advent of high-throughput sequencing and genome editing technologies led to the development of modern vaccines. We briefly discuss the impact of genomics in the ongoing COVID-19 pandemic in this review.

Consecutive deletions in a unique Uruguayan SARS-CoV-2 lineage evidence the genetic variability potential of accessory genes

Deletions frequently occur in the six accessory genes of SARS-CoV-2, but most genomes with deletions are sporadic and have limited spreading capability. Here, we analyze deletions in the ORF7a of the N.7 lineage, a unique Uruguayan clade from the Brazilian B.1.1.33 lineage. Thirteen samples collected during the early SARS-CoV-2 wave in Uruguay had deletions in the ORF7a. Complete genomes were obtained by Illumina next-generation sequencing, and deletions were confirmed by Sanger sequencing and capillary electrophoresis. The N.7 lineage includes several individuals with a 12-nucleotide deletion that removes four amino acids of the ORF7a. Notably, four individuals underwent an additional 68-nucleotide novel deletion that locates 44 nucleotides downstream in the terminal region of the same ORF7a. The simultaneous occurrence of the 12 and 68-nucleotide deletions fuses the ORF7a and ORF7b, two contiguous accessory genes that encode transmembrane proteins with immune-modulation activity. The fused ORF retains the signal peptide and the complete Ig-like fold of the 7a protein and the transmembrane domain of the 7b protein, suggesting that the fused protein plays similar functions to original proteins in a single format. Our findings evidence the remarkable dynamics of SARS-CoV-2 and the possibility that single and consecutive deletions occur in accessory genes and promote changes in the genomic organization that help the virus explore genetic variations and select for new, higher fit changes.

Characterization of SARS-CoV-2 ORF6 deletion variants detected in a nosocomial cluster during routine genomic surveillance, Lyon, France

During routine molecular surveillance of SARS-CoV-2 performed at the National Reference Center of Respiratory Viruses (Lyon, France) (n = 229 sequences collected February-April 2020), two frameshifting deletions were detected in the open reading frame 6, at the same position (27267). While a 26-nucleotide deletion variant (D26) was only found in one nasopharyngeal sample in March 2020, the 34-nucleotide deletion (D34) was found within a single geriatric hospital unit in 5/9 patients and one health care worker in April 2020. Phylogeny analysis strongly suggested a nosocomial transmission of D34, with potential fecal transmission, as also identified in a stool sample. No difference in disease severity was observed between patients hospitalized in the geriatric unit infected with WT or D34. In vitro D26 and D34 characterization revealed comparable replication kinetics with the wild-type (WT), but differential host immune responses. While interferon-stimulated genes were similarly upregulated after infection with WT and ORF6 variants, the latter specifically induced overexpression of 9 genes coding for inflammatory cytokines in the NF-kB pathway, including CCL2/MCP1, PTX3, and TNFα, for which high plasma levels have been associated with severe COVID-19. Our findings emphasize the need to monitor the occurrence of ORF6 deletions and assess their impact on the host immune response.

The Pandemic Experience in Southeast Asia: Interface Between SARS-CoV-2, Malaria, and Dengue

Southeast Asia (SEA) emerged relatively unscathed from the first year of the global SARS-CoV-2 pandemic, but as of July 2021 the region is experiencing a surge in case numbers primarily driven by Alpha (B.1.1.7) and subsequently the more transmissible Delta (B.1.617.2) variants. While initial disease burden was mitigated by swift government responses, favorable cultural and societal factors, the more recent rise in cases suggests an under-appreciation of prior prevalence and over-appreciation of possible cross-protective immunity from exposure to endemic viruses, and highlights the effects of vaccine rollout at varying tempos and of variable efficacy. This burgeoning crisis is further complicated by co-existence of malaria and dengue in the region, with implications of serological cross-reactivity on interpretation of SARS-CoV-2 assays and competing resource demands impacting efforts to contain both endemic and pandemic disease.

Expansion of a SARS-CoV-2 Delta variant with an 872 nt deletion encompassing ORF7a, ORF7b and ORF8, Poland, July to August 2021

Routine genomic surveillance on samples from COVID-19 patients collected in Poland during summer 2021 revealed the emergence of a SARS-CoV-2 Delta variant with a large 872 nt deletion. This change, confirmed by Sanger and deep sequencing, causes complete loss of ORF7a, ORF7b, and ORF8 genes. The index case carrying the deletion is unknown. The standard pipeline for sequencing may mask this deletion with a long stretch of N’s. Effects of this deletion on phenotype or immune evasion needs further study.

Insight into the LFA-1/SARS-CoV-2 Orf7a Complex by Protein-Protein Docking, Molecular Dynamics, and MM-GBSA Calculations

In the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) genome, open reading frames (ORFs) encode for viral accessory proteins. Among these, Orf7a structurally resembles the members of the immunoglobulin (Ig) superfamily and intracellular adhesion molecules (ICAMs), in particular. ICAMs are involved in integrin binding through lymphocyte function-associated antigen 1 (LFA-1). Based on such considerations and on previous findings on SARS-CoV, it has been postulated that the formation of the LFA-1/Orf7a complex could contribute to SARS-CoV-2 infectivity and pathogenicity. With the current work, we aim at providing insight into this macromolecular assembly, taking advantage of the recently reported SARS-CoV-2 Orf7a structure. Protein-protein docking, molecular dynamics (MD) simulations, and a Molecular Mechanical-Generalized Born Surface Area (MM-GBSA)-based stage were enrolled to provide refined models.

SARS-CoV-2 one year on: evidence for ongoing viral adaptation

SARS-CoV-2 is thought to have originated in the human population from a zoonotic spillover event. Infection in humans results in a variety of outcomes ranging from asymptomatic cases to the disease COVID-19, which can have significant morbidity and mortality, with over two million confirmed deaths worldwide as of January 2021. Over a year into the pandemic, sequencing analysis has shown that variants of SARS-CoV-2 are being selected as the virus continues to circulate widely within the human population. The predominant drivers of genetic variation within SARS-CoV-2 are single nucleotide polymorphisms (SNPs) caused by polymerase error, potential host factor driven RNA modification, and insertion/deletions (indels) resulting from the discontinuous nature of viral RNA synthesis. While many mutations represent neutral 'genetic drift' or have quickly died out, a subset may be affecting viral traits such as transmissibility, pathogenicity, host range, and antigenicity of the virus. In this review, we summarise the current extent of genetic change in SARS-CoV-2, particularly recently emerging variants of concern, and consider the phenotypic consequences of this viral evolution that may impact the future trajectory of the pandemic.

Pitfalls of exceptions for COVID-19 travel quarantine: lessons from a dignitary visit to Thailand

A dignitary from a European country made an official visit to Thailand in November 2020. Due to the nature of this visit, a standard 14-day quarantine was not implemented. After a series of meetings with diplomats and hotel staffers, the dignitary was diagnosed with COVID-19. Here, we detailed the transmission and the genomic investigation from the visit, presenting the factors that need to be considered when a standard quarantine was waivered.

Atypical Divergence of SARS-CoV-2 Orf8 from Orf7a within the Coronavirus Lineage Suggests Potential Stealthy Viral Strategies in Immune Evasion

Orf8, one of the most puzzling genes in the SARS lineage of coronaviruses, marks a unique and striking difference in genome organization between SARS-CoV-2 and SARS-CoV-1. Here, using sequence comparisons, we unequivocally reveal the distant sequence similarities between SARS-CoV-2 Orf8 with its SARS-CoV-1 counterparts and the X4-like genes of coronaviruses, including its highly divergent "paralog" gene Orf7a, whose product is a potential immune antagonist of known structure. Supervised sequence space walks unravel identity levels that drop below 10% and yet exhibit subtle conservation patterns in this novel superfamily, characterized by an immunoglobulin-like beta sandwich topology. We document the high accuracy of the sequence space walk process in detail and characterize the subgroups of the superfamily in sequence space by systematic annotation of gene and taxon groups. While SARS-CoV-1 Orf7a and Orf8 genes are most similar to bat virus sequences, their SARS-CoV-2 counterparts are closer to pangolin virus homologs, reflecting the fine structure of conservation patterns within the SARS-CoV-2 genomes. The divergence between Orf7a and Orf8 is exceptionally idiosyncratic, since Orf7a is more constrained, whereas Orf8 is subject to rampant change, a peculiar feature that may be related to hitherto-unknown viral infection strategies. Despite their common origin, the Orf7a and Orf8 protein families exhibit different modes of evolutionary trajectories within the coronavirus lineage, which might be partly attributable to their complex interactions with the mammalian host cell, reflected by a multitude of functional associations of Orf8 in SARS-CoV-2 compared to a very small number of interactions discovered for Orf7a.IMPORTANCE Orf8 is one of the most puzzling genes in the SARS lineage of coronaviruses, including SARS-CoV-2. Using sophisticated sequence comparisons, we confirm its origins from Orf7a, another gene in the lineage that appears as more conserved, compared to Orf8. Orf7a is a potential immune antagonist of known structure, while a deletion of Orf8 was shown to decrease the severity of the infection in a cohort study. The subtle sequence similarities imply that Orf8 has the same immunoglobulin-like fold as Orf7a, confirmed by structure determination. We characterize the subgroups of this superfamily and demonstrate the highly idiosyncratic divergence patterns during the evolution of the virus.

SARS-CoV-2 spread across the Colombian-Venezuelan border

INTRODUCTION

Venezuela and Colombia both adopted measures of containment early in response to the COVID-19 pandemic. However, Venezuela's ongoing humanitarian crisis has decimated its health care system, and forced millions of Venezuelans to flee through its porous border with Colombia. The extensive shared border, and illegal cross-border transit through improvised trails between the two countries are major challenges for public health authorities. We report the first SARS-CoV-2 genomes from Venezuela, and present a snapshot of the SARS-CoV-2 epidemiologic landscape in the Colombian-Venezuelan border region.

METHODS

We sequenced and assembled viral genomes from total RNA extracted from nasopharyngeal (NP) clinical specimens using a custom reference-based analysis pipeline. Three assemblies obtained were subjected to typing using the Phylogenetic Assignment of Named Global Outbreak LINeages 'Pangolin' tool. A total of 376 publicly available SARS-CoV-2 genomes from South America were obtained from the GISAID database to perform comparative genomic analyses. Additionally, the Wuhan-1 strain was used as reference.

RESULTS

We found that two of the SARS-CoV-2 genomes from Venezuela belonged to the B1 lineage, and the third to the B.1.13 lineage. We observed a point mutation in the Spike protein gene (D614G substitution), previously reported to be associated with increased infectivity, in all three Venezuelan genomes. Additionally, three mutations (R203K/G204R substitution) were present in the nucleocapsid (N) gene of one Venezuelan genome.

CONCLUSIONS

Genomic sequencing demonstrates similarity between SARS-CoV-2 lineages from Venezuela and viruses collected from patients in bordering areas in Colombia and from Brazil, consistent with cross-border transit despite administrative measures including lockdowns. The presence of mutations associated with increased infectivity in the 3 Venezuelan genomes we report and Colombian SARS-CoV-2 genomes from neighboring borders areas may pose additional challenges for control of SARS-CoV-2 spread in the complex epidemiological landscape in Latin American countries. Public health authorities should carefully follow the progress of the pandemic and its impact on displaced populations within the region.

Evaluation of NGS-based approaches for SARS-CoV-2 whole genome characterisation

Since the beginning of the COVID-19 outbreak, SARS-CoV-2 whole-genome sequencing (WGS) has been performed at unprecedented rate worldwide with the use of very diverse Next-Generation Sequencing (NGS) methods. Herein, we compare the performance of four NGS-based approaches for SARS-CoV-2 WGS. Twenty-four clinical respiratory samples with a large scale of Ct values (from 10.7 to 33.9) were sequenced with four methods. Three used Illumina sequencing: an in-house metagenomic NGS (mNGS) protocol and two newly commercialised kits including a hybridisation capture method developed by Illumina (DNA Prep with Enrichment kit and Respiratory Virus Oligo Panel, RVOP), and an amplicon sequencing method developed by Paragon Genomics (CleanPlex SARS-CoV-2 kit). We also evaluated the widely used amplicon sequencing protocol developed by ARTIC Network and combined with Oxford Nanopore Technologies (ONT) sequencing. All four methods yielded near-complete genomes (>99%) for high viral loads samples (n = 8), with mNGS and RVOP producing the most complete genomes. For mid viral loads (Ct 20–25), amplicon-based enrichment methods led to genome coverage >99 per cent for all samples while 1/8 sample sequenced with RVOP and 2/8 samples sequenced with mNGS had a genome coverage below 99 per cent. For low viral loads (Ct ≥25), amplicon-based enrichment methods were the most sensitive techniques. All methods were highly concordant in terms of identity in complete consensus sequence. Just one mismatch in three samples was observed in CleanPlex vs the other methods, due to the dedicated bioinformatics pipeline setting a high threshold to call SNP compared to reference sequence. Importantly, all methods correctly identified a newly observed 34nt-deletion in ORF6 but required specific bioinformatic validation for RVOP. Finally, as a major warning for targeted techniques, a loss of coverage in any given region of the genome should alert to a potential rearrangement or a SNP in primer-annealing or probe-hybridizing regions and would require further validation using unbiased metagenomic sequencing.