Cost-Effective Method to Perform SARS-CoV-2 Variant Surveillance: Detection of Alpha, Gamma, Lambda, Delta, Epsilon, and Zeta in Argentina

SARS-CoV-2 variants with concerning characteristics have emerged since the end of 2020. Surveillance of SARS-CoV-2 variants was performed on a total of 4,851 samples from the capital city and 10 provinces of Argentina, during 51 epidemiological weeks (EWs) that covered the end of the first wave and the ongoing second wave of the COVID-19 pandemic in the country (EW 44/2020 to EW 41/2021). The surveillance strategy was mainly based on Sanger sequencing of a Spike coding region that allows the identification of signature mutations associated with variants. In addition, whole-genome sequences were obtained from 637 samples. The main variants found were Gamma and Lambda, and to a lesser extent, Alpha, Zeta, and Epsilon, and more recently, Delta. Whereas, Gamma dominated in different regions of the country, both Gamma and Lambda prevailed in the most populated area, the metropolitan region of Buenos Aires. The lineages that circulated on the first wave were replaced by emergent variants in a term of a few weeks. At the end of the ongoing second wave, Delta began to be detected, replacing Gamma and Lambda. This scenario is consistent with the Latin American variant landscape, so far characterized by a concurrent increase in Delta circulation and a stabilization in the number of cases. The cost-effective surveillance protocol presented here allowed for a rapid response in a resource-limited setting, added information on the expansion of Lambda in South America, and contributed to the implementation of public health measures to control the disease spread in Argentina.

Near-Complete Genome of SARS-CoV-2 Delta (AY.3) Variant Identified in a Dog in Kansas, USA

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) descriptions of infection and transmission have been increasing in companion animals in the past year. Although canine susceptibility is generally considered low, their role in the COVID-19 disease cycle remains unknown. In this study, we detected and sequenced a delta variant (AY.3) from a 12-year-old Collie living with owners that previously tested positive for SARS-CoV-2. It is unclear if the dogs’ symptoms were related to SARS-CoV-2 infection or underlying conditions. The whole genome sequence obtained from the dog sample had several unique consensus level changes not previously identified in a SARS-CoV-2 genome that may play a role in the rapid adaptation from humans to dogs. Within the spike coding region, 5/7 of the subconsensus variants identified in the dog sequence were also identified in the closest in-house human reference case. Taken together, the whole genome sequence, and phylogenetic and subconsensus variant analyses indicate the virus infecting the animal originated from a local outbreak cluster. The results of these analyses emphasize the importance of rapid detection and characterization of SARS-CoV-2 variants of concern in companion animals.

B.1.617.3 SARS CoV-2 spike E156G/Δ157-158 mutations contribute to reduced neutralization sensitivity and increased infectivity

SARS CoV-2 variants raise significant concerns due to their ability to cause vaccine breakthrough infections. Here, we sequence-characterized the spike gene, isolated from a breakthrough infection, that corresponded to B.1.617.3 lineage. Delineating the functional impact of spike mutations using reporter pseudoviruses (PV) revealed that N-terminal domain (NTD)-specific E156G/Δ157-158 contributed to increased infectivity and reduced sensitivity to ChAdOx1 nCoV-19 vaccine (CovishieldTM)-elicited neutralizing antibodies. A six-nucleotide deletion (467-472) in the spike coding region introduced this change in the NTD. We confirmed the presence of E156G/Δ157-158 in the RT-PCR-positive cases concurrently screened, in addition to other circulating spike (S1) mutations like T19R, T95I, L452R, E484Q, and D614G. Notably, E156G/Δ157-158 was present in more than 85% of the sequences reported from the USA, UK, and India in August 2021. The spike PV bearing combination of E156G/Δ157-158 and L452R further promoted infectivity and conferred immune evasion. Additionally, increased cell-to-cell fusion was observed when spike harbored E156G/Δ157-158, L452R, and E484Q, suggesting a combinatorial effect of these mutations. Notwithstanding, the plasma from a recovered individual robustly inhibited mutant spike PV, indicating the increased breadth of neutralization post-recovery. Our data highlights the importance of spike NTD-specific changes in determining infectivity and immune escape of variants.

Surveillance of SARS-CoV-2 variants in Argentina: detection of Alpha, Gamma, Lambda, Epsilon and Zeta in locally transmitted and imported cases

Molecular surveillance of SARS-CoV-2 variants was performed on a total of 2,406 samples from the capital city and nine provinces of Argentina, during 30 epidemiological weeks (EW) that covered the end of the first wave and the beginning of the ongoing second wave of the COVID-19 pandemic in the country (EW 44/2020 to EW 20/2021). The surveillance strategy was mainly based on Sanger sequencing of a Spike coding region that allows the simultaneous identification of signature mutations associated with worldwide circulating variants. In addition, whole SARS-CoV-2 genome sequences were obtained from 456 samples. The main variants found were Gamma, Lambda and Alpha, and to a lesser extent, Zeta and Epsilon. Whereas Gamma dominated in different regions of the country, both Gamma and Lambda prevailed in the most populated area, the metropolitan region of Buenos Aires (MABA), although showing a heterogeneous distribution along this region. This cost-effective surveillance protocol allowed for a rapid response in a limited access to resources scenario, added information on the expansion of the Lambda variant in South America and contributed to the implementation of public health measures to control the disease spread in Argentina.

SARS-CoV-2 and other human coronavirus show genome patterns previously associated to reduced viral recognition and altered immune response

A new pandemic caused by the betacoronavirus SARS-CoV-2 originated in China in late 2019. Although often asymptomatic, a relevant percentage of affected people can develop severe pneumonia. Initial evidence suggests that dysregulation of the immune response could contribute to the pathogenesis, as previously demonstrated for SARS-CoV. The presence of genome composition features involved in delaying viral recognition is herein investigated for human coronaviruses (HCoVs), with a special emphasis on SARS-CoV-2. A broad collection of HCoVs polyprotein, envelope, matrix, nucleocapsid and spike coding sequences was downloaded and several statistics representative of genome composition and codon bias were investigated. A model able to evaluate and test the presence of a significant under- or over-representation of dinucleotide pairs while accounting for the underlying codon bias and protein sequence was also implemented. The study revealed the significant under-representation of CpG dinucleotide pair in all HcoV, but especially in SARS-CoV and even more in SARS-CoV-2. The presence of forces acting to minimize CpG content was confirmed by relative synonymous codon usage pattern. Codons containing the CpG pair were severely under-represented, primarily in the polyprotein and spike coding sequences of SARS-CoV-2. Additionally, a significant under-representation of the TpA pair was observed in the N and S region of SARS-CoV and SARS-CoV-2. Increasing experimental evidence has proven that CpG and TpA are targeted by innate antiviral host defences, contributing both to RNA degradation and RIG-1 mediated interferon production. The low content of these dinucleotides could contribute to a delayed interferon production, dysregulated immune response, higher viral replication and poor outcome. Significantly, the RIG-1 signalling pathway was proven to be defective in elderlies, suggesting a likely interaction between limited viral recognition and lower responsiveness in interferon production that could justify the higher disease severity and mortality in older patients.

SARS-CoV-2 And Other Human Coronavirus Show Genome Patterns Previously Associated to Reduced Viral Recognition and Altered Immune Response

A new pandemic caused by the betacoronavirus SARS-CoV-2 originated in China in late 2019. Although often asymptomatic, a relevant percentage of affected people can develop severe pneumonia. Initial evidence suggests that dysregulation of the immune response could contribute to the pathogenesis, as previously demonstrated for SARS-CoV. The presence of genome composition features involved in delaying viral recognition is herein investigated for human coronaviruses (HCoVs), with a special emphasis on SARS-CoV-2. A broad collection of HCoVs polyprotein, envelope, matrix, nucleocapsid and spike coding sequences was downloaded and several statistics representative of genome composition and codon bias were investigated. A model able to evaluate and test the presence of a significant under- or over-representation of dinucleotide pairs while accounting for the underlying codon bias and protein sequence was also implemented. The study revealed the significant under-representation of CpG dinucleotide pair in all HcoV, but especially in SARS-CoV and even more in SARS-CoV-2. The presence of forces acting to minimize CpG content was confirmed by relative synonymous codon usage pattern. Codons containing the CpG pair were severely under-represented, primarily in the polyprotein and spike coding sequences of SARS-CoV-2. Additionally, a significant under-representation of the TpA pair was observed in the N and S region of SARS-CoV and SARS-CoV-2. Increasing experimental evidence has proven that CpG and TpA are targeted by innate antiviral host defences, contributing both to RNA degradation and RIG-1 mediated interferon production. The low content of these dinucleotides could contribute to a delayed interferon production, dysregulated immune response, higher viral replication and poor outcome. Significantly, the RIG-1 signalling pathway was proven to be defective in elderlies, suggesting a likely interaction between limited viral recognition and lower responsiveness in interferon production that could justify the higher disease severity and mortality in older patients.

The heterogeneous landscape and early evolution of pathogen-associated CpG dinucleotides in SARS-CoV-2

COVID-19 can lead to acute respiratory syndrome, which can be due to dysregulated immune signaling. We analyze the distribution of CpG dinucleotides, a pathogen-associated molecular pattern, in the SARS-CoV-2 genome. We characterize CpG content by a CpG force that accounts for statistical constraints acting on the genome at the nucleotidic and amino-acid levels. The CpG force, as the CpG content, is overall low compared to other pathogenic betacoronaviruses; however it widely fluctuates along the genome, with a particularly low value, comparable to the circulating seasonal HKU1, in the spike coding region and a greater value, comparable to SARS and MERS, in the highly expressed nucleocapside coding region (N ORF), whose transcripts are relatively abundant in the cytoplasm of infected cells and present in the 3’UTRs of all subgenomic RNA. This dual nature of CpG content could confer to SARS-CoV-2 the ability to avoid triggering pattern recognition receptors upon entry, while eliciting a stronger response during replication. We then investigate the evolution of synonymous mutations since the outbreak of the COVID-19 pandemic, finding a signature of CpG loss in regions with a greater CpG force. Sequence motifs preceding the CpG-loss-associated loci in the N ORF match recently identified binding patterns of the Zinc finger Anti-viral Protein. Using a model of the viral gene evolution under human host pressure, we find that synonymous mutations seem driven in the SARS-CoV-2 genome, and particularly in the N ORF, by the viral codon bias, the transition-transversion bias and the pressure to lower CpG content.