Emergence and outcomes of the SARS-CoV-2 'Marseille-4' variant

Health Status and COVID-19 Epidemiology in an Inland Region of Portugal: A Retrospective Study

Multiple factors, from socioeconomic development to genetic background, can affect the regional impact of some diseases, and this has also been seen during the COVID-19 pandemic. The objective of this retrospective study was to characterize a population in the interior of Portugal regarding health status and COVID-19 epidemiology. Between October 2021 and January 2023, 1553 subjects residing in Beira Baixa, Portugal, were included. Using a self-report approach, demographic and clinical data were obtained. Blood group, blood pressure, peripheral oxygen saturation and anti-spike protein immunoglobulin concentration were also analyzed. Statistical analysis was performed using IBM SPSS Statistics. The average age of the participants was 48.95 (±14.43) years, with 64% being male and 36% being female. The most prevalent comorbidities were hypertension (19.2%), dyslipidemia (12.6%) and diabetes mellitus (6.6%). Half of the population was overweight, and more than half of the subjects had no history of tobacco consumption. Among the participants, 33% were infected with SARS-CoV-2: 70.1% had mild disease, 14.1% moderate disease and 1.4% severe disease. There was a very significant adherence to vaccination (97%). Previously infected or vaccinated people had higher anti-spike protein immunoglobulin values; this value depended on the vaccine administered (p < 0.001). Patients with autoimmune diseases and smokers had lower levels of anti-S IgG antibodies (p = 0.030 and p = 0.024, respectively). The severity of COVID-19 did not affect the concentration of anti-S IgG (p = 0.430). This study highlights the general health statuses and the impact of COVID-19 on a population in the Portuguese interior. Knowledge of the circulation and impact of the virus in this specific population can alert and assist in better interventions being conducted by health authorities.

Description of nasopharyngeal bacterial pathogens associated with different SARS-CoV-2 variants

The emergence of the coronavirus pandemic facilitated the acquisition of mutations in the SARS-CoV-2 genome, resulting in the appearance of new variants over the past three years. We previously identified several taxa associated with the clinical outcome of COVID-19 disease in a retrospective study involving 120 patients (infected patients and negative subjects). However, little is known about whether the different variants could influence variations in the composition of the nasopharyngeal microbiota. In this study, we used multiplex pathogen-specific PCR to analyse the presence of nasopharyngeal bacterial pathogens from 400 SARS-CoV-2 positive patients (equally distributed in the four SARS-CoV-2 variants studied: B.1.1.7 (Alpha), B.1 0.617.2 (Delta), B.1.160 (Marseille-4), and B.1.1.529 (omicron)). We then compared them to 400 patients who tested negative for all respiratory viruses tested in this study, including SARS-CoV-2. We first observed an enrichment of Staphylococcus aureus (P ≤ .05) and Corynebacterium propinquum (P ≤ .05) in COVID-19-positive patients, regardless of the variant, compared to negative subjects. We specifically highlighted a significantly higher frequency of S. aureus (P ≤ .0001), C. propinquum (P ≤ .0001), and Klebsiella pneumoniae (P ≤ .0001), in patients infected with the omicron variant, whereas that of Haemophilus influenzae was higher in patients infected with Marseille-4 (P ≤ .001) and Alpha (P ≤ .01) variants. Our results suggest that the nasopharyngeal bacterial pathogens have their own specificity according to the SARS-CoV-2 variant and independently of the season. Additional studies are needed to determine the role of these pathogens in the evolution of the clinical outcome of patients.

[Corticoid treatment and SARS-CoV-2 variants: two independent factors associated with COVID-19 mortality in a Spanish regional hospital]

BACKGROUND

Pandemic inter-wave hospital admissions and COVID-19-related mortality rates vary greatly. Some of the factors that may be playing part in this are the profile of the patients, viral variants, pharmacological treatments, or preventive measures. This work aimed to analyze the factors associated with mortality in COVID-19 patients admitted to hospital during 2020-2021.

METHODS

Retrospective cohort study with COVID-19 patients admitted to Hospital de Barbastro (Spain) during 2020-2021. Data were collected from the Spanish Conjunto Mínimo Básico de Datos and microbiology and electronic prescription records.

RESULTS

During the study period, 908 patients were consecutively admitted for COVID-19 (median age 70 years, 57.2% males); 162 (17.8%) patients died. We identified seven successive epidemiological waves. The following variables significantly associated to higher mortality: age, arterial hypertension, chronic renal failure, dementia, chronic obstructive pulmonary disease, heart failure, prior stroke, Charlson index, and wave 2; wave 4 was associated to greater survival. The multivariate analysis showed that age (OR=1.11; 95% CI: 1.09-1.14), chronic obstructive pulmonary disease (OR=2.33; 95% CI: 1.18-4.57), wave 2 (OR=2.57; 95% CI: 1.10-6.00), and wave 3 (OR=2.94; 95% CI: 1.17-7.38) associated with higher mortality. Glucocorticoid treatment was the only protective factor (OR=0.29; 95%CI: 0.14-0.62).

CONCLUSIONS

This study confirms the therapeutic utility of glucocorticoids to reduce in-hospital mortality due to COVID-19. Heterogeneous mortality rates between the different COVID-19 waves suggest a direct role of viral variants as determinants of lethality, regardless of the patient's history.

High Throughput SARS-CoV-2 Genome Sequencing from 384 Respiratory Samples Using the Illumina COVIDSeq Protocol

The emergence of the Coronavirus Disease 2019 (COVID-19) pandemic has fostered the use of high-throughput techniques to sequence the entire severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome and track its evolution. The present study proposes a rapid and relatively less expensive sequencing protocol for 384 samples by adapting the use of an Illumina NovaSeq library to an Illumina MiSeq flow cell instrument. The SARS-CoV-2 genome sequences obtained with Illumina NovaSeq and those obtained using MiSeq instruments were compared with the objective to validate the new, modified protocol. A total of 356 (94.6%) samples yielded interpretable sequences using the modified Illumina COVIDSeq protocol, with an average coverage of 91.6%. By comparison, 357 (94.9%) samples yielded interpretable sequences with the standard COVIDSeq protocol, with an average coverage of 95.6%. Our modified COVIDSeq protocol could save 14,155 euros per run and yield results from 384 samples in 53.5 h, compared to four times 55.5 h with the standard Illumina MiSeq protocol. The modified COVIDSeq protocol thus provides high quality results comparable to those obtained with the standard COVIDSeq protocol, four times faster, while saving money.

mRNA COVID-19 Vaccines-Facts and Hypotheses on Fragmentation and Encapsulation

BACKGROUND

The adventure of the mRNA vaccine began thirty years ago in the context of influenza. This consisted in encapsulating the mRNA coding for a viral protein in a lipid particle. We show how the mRNA encoding S protein has been modified for that purpose in the context of the anti-SARS-CoV-2 vaccination.

RESULTS

by using data coming from genetic and epidemiologic databases, we show the theoretical possibility of fragmentation of this mRNA into small RNA sequences capable of inhibiting important bio-syntheses such as the production of beta-globin.

DISCUSSION

we discuss two aspects related to mRNA vaccine: (i) the plausibility of mRNA fragmentation, and (ii) the role of liposomal nanoparticles (LNPs) used in the vaccine and their impact on mRNA biodistribution.

CONCLUSION

we insist on the need to develop lipid nanoparticles allowing personalized administration of vaccines and avoiding adverse effects due to mRNA fragmentation and inefficient biodistribution. Hence, we recommend (i) adapting the mRNA of vaccines to the least mutated virus proteins and (ii) personalizing its administration to the categories of chronic patients at risk most likely to suffer from adverse effects.

Quasispecies Analysis of SARS-CoV-2 of 15 Different Lineages during the First Year of the Pandemic Prompts Scratching under the Surface of Consensus Genome Sequences

The tremendous majority of SARS-CoV-2 genomic data so far neglected intra-host genetic diversity. Here, we studied SARS-CoV-2 quasispecies based on data generated by next-generation sequencing (NGS) of complete genomes. SARS-CoV-2 raw NGS data had been generated for nasopharyngeal samples collected between March 2020 and February 2021 by the Illumina technology on a MiSeq instrument, without prior PCR amplification. To analyze viral quasispecies, we designed and implemented an in-house Excel file (“QuasiS”) that can characterize intra-sample nucleotide diversity along the genomes using data of the mapping of NGS reads. We compared intra-sample genetic diversity and global genetic diversity available from Nextstrain. Hierarchical clustering of all samples based on the intra-sample genetic diversity was performed and visualized with the Morpheus web application. NGS mapping data from 110 SARS-CoV-2-positive respiratory samples characterized by a mean depth of 169 NGS reads/nucleotide position and for which consensus genomes that had been obtained were classified into 15 viral lineages were analyzed. Mean intra-sample nucleotide diversity was 0.21 ± 0.65%, and 5357 positions (17.9%) exhibited significant (>4%) diversity, in ≥2 genomes for 1730 (5.8%) of them. ORF10, spike, and N genes had the highest number of positions exhibiting diversity (0.56%, 0.34%, and 0.24%, respectively). Nine hot spots of intra-sample diversity were identified in the SARS-CoV-2 NSP6, NSP12, ORF8, and N genes. Hierarchical clustering delineated a set of six genomes of different lineages characterized by 920 positions exhibiting intra-sample diversity. In addition, 118 nucleotide positions (0.4%) exhibited diversity at both intra- and inter-patient levels. Overall, the present study illustrates that the SARS-CoV-2 consensus genome sequences are only an incomplete and imperfect representation of the entire viral population infecting a patient, and that quasispecies analysis may allow deciphering more accurately the viral evolutionary pathways.

Limited permissibility of ENL-R and Mv-1-Lu mink cell lines to SARS-CoV-2

The SARS-CoV-2 pandemic started in the end of 2019 in Wuhan, China, which highlighted the scenario of frequent cross-species transmission events. From the outbreak possibly initiated by viral spill-over into humans from an animal reservoir, now we face the human host moving globally while interacting with domesticated and peridomestic animals. The emergence of a new virus into the ecosystem leads to selecting forces and species-specific adaptations. The adaptation of SARS-CoV-2 to other animals represents a risk to controlling the dissemination of this coronavirus and the emergence of new variants. Since 2020, several mink farms in Europe and the United States have had SARS-CoV-2 outbreaks with human-mink and mink-human transmission, where the mink-selected variants possibly hold evolutionary concerning advantages. Here we investigated the permissibility of mink lung-derived cells using two cell lines, Mv-1-Lu and ENL-R, against several lineages of SARS-CoV-2, including some classified as variants of concern. The viral release rate and the infectious titers indicate that these cells support infections by different SARS-CoV-2 lineages. The viral production occurs in the first few days after infection with the low viral release by these mink cells, which is often absent for the omicron variant for lung cells. The electron microscopy reveals that during the viral replication cycle, the endomembrane system of the mink-host cell undergoes typical changes while the viral particles are produced, especially in the first days of infection. Therefore, even if limited, mink lung cells may represent a selecting source for SARS-CoV-2 variants, impacting their transmissibility and pathogenicity and making it difficult to control this new coronavirus.

The emergence, spread and vanishing of a French SARS-CoV-2 variant exemplifies the fate of RNA virus epidemics and obeys the Mistigri rule

The nature and dynamics of mutations associated with the emergence, spread, and vanishing of SARS-CoV-2 variants causing successive waves are complex. We determined the kinetics of the most common French variant ("Marseille-4") for 10 months since its onset in July 2020. Here, we analyzed and classified into subvariants and lineages 7453 genomes obtained by next-generation sequencing. We identified two subvariants, Marseille-4A, which contains 22 different lineages of at least 50 genomes, and Marseille-4B. Their average lifetime was 4.1 ± 1.4 months, during which 4.1 ± 2.6 mutations accumulated. Growth rate was 0.079 ± 0.045, varying from 0.010 to 0.173. Most of the lineages exhibited a bell-shaped distribution. Several beneficial mutations at unpredicted sites initiated a new outbreak, while the accumulation of other mutations resulted in more viral heterogenicity, increased diversity and vanishing of the lineages. Marseille-4B emerged when the other Marseille-4 lineages vanished. Its ORF8 gene was knocked out by a stop codon, as reported in SARS-CoV-2 of mink and in the Alpha variant. This subvariant was associated with increased hospitalization and death rates, suggesting that ORF8 is a nonvirulence gene. We speculate that the observed heterogenicity of a lineage may predict the end of the outbreak.

A Single Early Introduction Governed Viral Diversity in the Second Wave of SARS-CoV-2 Epidemic in Hungary

Retrospective evaluation of past waves of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemic is key for designing optimal interventions against future waves and novel pandemics. Here, we report on analysing genome sequences of SARS-CoV-2 from the first two waves of the epidemic in 2020 in Hungary, mirroring a suppression and a mitigation strategy, respectively. Our analysis reveals that the two waves markedly differed in viral diversity and transmission patterns. Specifically, unlike in several European areas or in the USA, we have found no evidence for early introduction and cryptic transmission of the virus in the first wave of the pandemic in Hungary. Despite the introduction of multiple viral lineages, extensive community spread was prevented by a timely national lockdown in March 2020. In sharp contrast, the majority of the cases in the much larger second wave can be linked to a single transmission lineage of the pan-European B.1.160 variant. This lineage was introduced unexpectedly early, followed by a 2-month-long cryptic transmission before a soar of detected cases in September 2020. Epidemic analysis has revealed that the dominance of this lineage in the second wave was not associated with an intrinsic transmission advantage. This finding is further supported by the rapid replacement of B.1.160 by the alpha variant (B.1.1.7) that launched the third wave of the epidemic in February 2021. Overall, these results illustrate how the founder effect in combination with the cryptic transmission, instead of repeated international introductions or higher transmissibility, can govern viral diversity.

First evidence of human-to-dog transmission of SARS-CoV-2 B.1.160 variant in France

Since the start of the coronavirus disease of 2019 (COVID-19) pandemic, several episodes of human-to-animal severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission have been described in different countries. The role of pets, especially domestic dogs, in the COVID-19 epidemiology is highly questionable and needs further investigation. In this study, we report a case of COVID-19 in a French dog living in close contact with its owners who were COVID-19 patients. The dog presented rhinitis and was sampled 1 week after its owners (a man and a woman) were tested positive for COVID-19. The nasal swabs for the dog tested remained positive for SARS-CoV-2 by reverse transcription quantitative real-time PCR (RT-qPCR) 1 month following the first diagnosis. Specific anti-SARS-CoV-2 antibodies were detectable 12 days after the first diagnosis and persisted for at least 5 months as tested using enzyme-linked immunoassay (ELISA) and automated western blotting. The whole-genome sequences from the dog and its owners were 99%-100% identical (with the man and the woman's sequences, respectively) and matched the B.1.160 variant of concern (Marseille-4 variant), the most widespread in France at the time the dog was infected. This study documents the first detection of B.1.160 in pets (a dog) in France, and the first canine genome recovery of the B.1.160 variant of global concern. Moreover, given the enhanced infectivity and transmissibility of the Marseille-4 variant for humans, this case also highlights the risk that pets may potentially play a significant role in SARS-CoV-2 outbreaks and may transmit the infection to humans. We have evidence of human-to-dog transmission of the Marseille-4 variant since the owners were first to be infected. Finally, owners and veterinarians must be vigilent for canine COVID-19 when dogs are presented with respiratory clinical signs.

Sequential Appearance and Isolation of a SARS-CoV-2 Recombinant between Two Major SARS-CoV-2 Variants in a Chronically Infected Immunocompromised Patient

Genetic recombination is a major evolutionary mechanism among RNA viruses, and it is common in coronaviruses, including those infecting humans. A few SARS-CoV-2 recombinants have been reported to date whose genome harbored combinations of mutations from different mutants or variants, but only a single patient's sample was analyzed, and the virus was not isolated. Here, we report the gradual emergence of a hybrid genome of B.1.160 and Alpha variants in a lymphoma patient chronically infected for 14 months, and we isolated the recombinant virus. The hybrid genome was obtained by next-generation sequencing, and the recombination sites were confirmed by PCR. This consisted of a parental B.1.160 backbone interspersed with two fragments, including the spike gene, from an Alpha variant. An analysis of seven sequential samples from the patient decoded the recombination steps, including the initial infection with a B.1.160 variant, then a concurrent infection with this variant and an Alpha variant, the generation of hybrid genomes, and eventually the emergence of a predominant recombinant virus isolated at the end of the patient's follow-up. This case exemplifies the recombination process of SARS-CoV-2 in real life, and it calls for intensifying the genomic surveillance in patients coinfected with different SARS-CoV-2 variants, and more generally with several RNA viruses, as this may lead to the appearance of new viruses.

Screening and Whole Genome Sequencing of SARS-CoV-2 Circulating During the First Three Waves of the COVID-19 Pandemic in Libreville and the Haut-Ogooué Province in Gabon

Since the onset of the COVID-19 pandemic, the SARS-CoV-2 viral dynamics in Africa have been less documented than on other continents. In Gabon, a Central African country, a total number of 37,511 cases of COVID-19 and 281 deaths have been reported as of December 8, 2021. After the first COVID-19 case was reported on March 12, 2020, in the capital Libreville, the country experienced two successive waves. The first one, occurred in March 2020 to August 2020, and the second one in January 2021 to May 2021. The third wave began in September 2021 and ended in November 2021. In order to reduce the data gap regarding the dynamics of SARS-CoV-2 in Central Africa, we performed a retrospective genotyping study using 1,006 samples collected from COVID-19 patients in Gabon from 2020 to 2021. Using SARS-CoV-2 variant screening by Real-Time Quantitative Reverse Transcription PCR (qRT-PCR) and whole genome sequencing (WGS), we genotyped 809 SARS-CoV-2 samples through qRT-PCR and identified to generated 291 new genomes. It allowed us to describe specific mutations and changes in the SARS-CoV-2 variants in Gabon. The qRT-PCR screening of 809 positive samples from March 2020 to September 2021 showed that 119 SARS-CoV-2 samples (14.7%) were classified as VOC Alpha (Pangolin lineage B.1.1.7), one (0.1%) was a VOC Beta (B.1.351), and 198 (24.5 %) were VOC Delta (B.1.617.2), while 491 samples (60.7%) remained negative for the variants sought. The B1.1 variant was predominant during the first wave while the VOC Alpha dominated the second wave. The B1.617.2 Delta variant is currently the dominant variant of the third wave. Similarly, the analysis of the 291 genome sequences indicated that the dominant variant during the first wave was lineage B.1.1, while the dominant variants of the second wave were lineages B.1.1.7 (50.6%) and B.1.1.318 (36.4%). The third wave started with the circulation of the Delta variant (B.1.617). Finally, we compared these results to the SARS-CoV-2 sequences reported in other African, European, American and Asian countries. Sequences of Gabonese SARS-CoV-2 strains presented the highest similarities with those of France, Belgium and neighboring countries of Central Africa, as well as West Africa.

Who Were Hospitalized Deceased Patients from COVID-19 During the First Year of Pandemic? Retrospective Analysis of 1104 Deceased Patients in South of France

Introduction

Following the first year of the COVID-19 pandemic, a complete analysis of the characteristics of the deceased hospitalized patients was performed, to identify factors related to premature mortality and to compare patient profiles according to the epidemic periods.

Methods

Retrospective analysis of 1104 deceased patients in two University Hospitals in South-eastern France, between March 1, 2020 and March 12, 2021 from Hospital’s electronic medical records was performed.

Results

Mean age was 80 years (± 11.1) and 10% of the deceased were younger than 65 years with specific comorbidities, e.g., genetic conditions, metastatic cancer, or massive obesity. Among the three clusters identified, two clusters (75% of deceased patients) include very elderly patients with numerous comorbidities, and differ by their proportion of dependent institutionalized patients. The third cluster is made up of younger patients with fewer but severe comorbidities. Deceased patients’ profiles varied according to the epidemic periods: during the first period (March–June 2020), more patients were institutionalized. The second period (September–December2020) coincided with a higher mortality rate.

Conclusions

This study confirmed that most patients hospitalized and dying from COVID-19 were frail, i.e., elderly and/or highly comorbid and that the small proportion of young patients had severe comorbidities.

Supplementary Information

The online version contains supplementary material available at 10.1007/s44197-022-00039-3.

Antiviral Activity of Repurposing Ivermectin against a Panel of 30 Clinical SARS-CoV-2 Strains Belonging to 14 Variants

Over the past two years, several variants of SARS-CoV-2 have emerged and spread all over the world. However, infectivity, clinical severity, re-infection, virulence, transmissibility, vaccine responses and escape, and epidemiological aspects have differed between SARS-CoV-2 variants. Currently, very few treatments are recommended against SARS-CoV-2. Identification of effective drugs among repurposing FDA-approved drugs is a rapid, efficient and low-cost strategy against SARS-CoV-2. One of those drugs is ivermectin. Ivermectin is an antihelminthic agent that previously showed in vitro effects against a SARS-CoV-2 isolate (Australia/VI01/2020 isolate) with an IC50 of around 2 µM. We evaluated the in vitro activity of ivermectin on Vero E6 cells infected with 30 clinically isolated SARS-CoV-2 strains belonging to 14 different variants, and particularly 17 strains belonging to six variants of concern (VOC) (variants related to Wuhan, alpha, beta, gamma, delta and omicron). The in vitro activity of ivermectin was compared to those of chloroquine and remdesivir. Unlike chloroquine (EC50 from 4.3 ± 2.5 to 29.3 ± 5.2 µM) or remdesivir (EC50 from 0.4 ± 0.3 to 25.2 ± 9.4 µM), ivermectin showed a relatively homogeneous in vitro activity against SARS-CoV-2 regardless of the strains or variants (EC50 from 5.1 ± 0.5 to 6.7 ± 0.4 µM), except for one omicron strain (EC50 = 1.3 ± 0.5 µM). Ivermectin (No. EC50 = 219, mean EC50 = 5.7 ± 1.0 µM) was, overall, more potent in vitro than chloroquine (No. EC50 = 214, mean EC50 = 16.1 ± 9.0 µM) (p = 1.3 × 10-34) and remdesivir (No. EC50 = 201, mean EC50 = 11.9 ± 10.0 µM) (p = 1.6 × 10-13). These results should be interpreted with caution regarding the potential use of ivermectin in SARS-CoV-2-infected patients: it is difficult to translate in vitro study results into actual clinical treatment in patients.

The Polarity and Specificity of Antiviral T Lymphocyte Responses Determine Susceptibility to SARS-CoV-2 Infection in Patients with Cancer and Healthy Individuals

Vaccination against coronavirus disease 2019 (COVID-19) relies on the in-depth understanding of protective immune responses to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). We characterized the polarity and specificity of memory T cells directed against SARS-CoV-2 viral lysates and peptides to determine correlates with spontaneous, virus-elicited, or vaccine-induced protection against COVID-19 in disease-free and cancer-bearing individuals. A disbalance between type 1 and 2 cytokine release was associated with high susceptibility to COVID-19. Individuals susceptible to infection exhibited a specific deficit in the T helper 1/T cytotoxic 1 (Th1/Tc1) peptide repertoire affecting the receptor binding domain of the spike protein (S1-RBD), a hotspot of viral mutations. Current vaccines triggered Th1/Tc1 responses in only a fraction of all subject categories, more effectively against the original sequence of S1-RBD than that from viral variants. We speculate that the next generation of vaccines should elicit Th1/Tc1 T-cell responses against the S1-RBD domain of emerging viral variants.

SIGNIFICANCE

This study prospectively analyzed virus-specific T-cell correlates of protection against COVID-19 in healthy and cancer-bearing individuals. A disbalance between Th1/Th2 recall responses conferred susceptibility to COVID-19 in both populations, coinciding with selective defects in Th1 recognition of the receptor binding domain of spike. See related commentary by McGary and Vardhana, p. 892. This article is highlighted in the In This Issue feature, p. 873.

Analysis of SARS-CoV-2 Variants From 24,181 Patients Exemplifies the Role of Globalization and Zoonosis in Pandemics

After the end of the first epidemic episode of SARS-CoV-2 infections, as cases began to rise again during the summer of 2020, we at IHU Méditerranée Infection in Marseille, France, intensified the genomic surveillance of SARS-CoV-2, and described the first viral variants. In this study, we compared the incidence curves of SARS-CoV-2-associated deaths in different countries and reported the classification of SARS-CoV-2 variants detected in our institute, as well as the kinetics and sources of the infections. We used mortality collected from a COVID-19 data repository for 221 countries. Viral variants were defined based on ≥5 hallmark mutations along the whole genome shared by ≥30 genomes. SARS-CoV-2 genotype was determined for 24,181 patients using next-generation genome and gene sequencing (in 47 and 11% of cases, respectively) or variant-specific qPCR (in 42% of cases). Sixteen variants were identified by analyzing viral genomes from 9,788 SARS-CoV-2-diagnosed patients. Our data show that since the first SARS-CoV-2 epidemic episode in Marseille, importation through travel from abroad was documented for seven of the new variants. In addition, for the B.1.160 variant of Pangolin classification (a.k.a. Marseille-4), we suspect transmission from farm minks. In conclusion, we observed that the successive epidemic peaks of SARS-CoV-2 infections are not linked to rebounds of viral genotypes that are already present but to newly introduced variants. We thus suggest that border control is the best mean of combating this type of introduction, and that intensive control of mink farms is also necessary to prevent the emergence of new variants generated in this animal reservoir.

Occurrence of a substitution or deletion of SARS-CoV-2 spike amino acid 677 in various lineages in Marseille, France

Great concerns have been raised about SARS-CoV-2 variants over the past six months. At the end of 2020, an increasing incidence of spike substitutions Q677H/P was described in the USA, which involved six independent lineages. We searched for changes to this amino acid in the sequence database of SARS-CoV-2 genomes obtained at the IHU Méditerranée Infection (Marseille, France) from 3634 patients sampled between February 2020 and April 2021. In seven genomes (0.2%), we found a deletion of five amino acids at spike positions 675-679 (QTQTN) including Q677, and in 76 genomes (2.3%) we found a Q677H substitution. The 83 genomes were classified in ten different Pangolin lineages. Genomes with a spike Q677 deletion were obtained from respiratory samples collected in six cases between 28 March 2020 and 12 October 2020 and in one case on 1 February 2021. The Q677H substitution was found in genomes all obtained from respiratory samples collected from 19 January 2021 and were classified in seven different lineages. Most of these genomes (41 cases) were of UK variant. Two others were classified in the B.1.160 Pangolin lineage (Marseille-4 variant) which was first detected in July 2020 in our institute but was devoid of this substitution until 19 January 2021. Also, eight genomes were classified in the A.27/Marseille-501 lineage which was first detected in our institute in January 2021 and which either harboured or did not harbour the Q677H substitution. Thus, the spike Q677H substitution should be considered as another example of convergent evolution, as it is the case of spike substitutions L18F, E484K, L452R, and N501Y which also independently appeared in various lineages.

SARS-CoV-2 Variants and Clinical Outcomes: A Systematic Review

Background: From the start of the COVID-19 pandemic, new SARS-CoV-2 variants have emerged that potentially affect transmissibility, severity, and immune evasion in infected individuals. In the present systematic review, the impact of different SARS-CoV-2 variants on clinical outcomes is analyzed. Methods: A systematic review was performed according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) 2020. Two databases (PubMed and ScienceDirect) were searched for original articles published from 1 January 2020 to 23 November 2021. The articles that met the selection criteria were appraised according to the Newcastle–Ottawa Quality Assessment Scale. Results: Thirty-three articles were included, involving a total of 253,209 patients and 188,944 partial or complete SARS-CoV-2 sequences. The most reported SARS-CoV-2 variants showed changes in the spike protein, N protein, RdRp and NSP3. In 28 scenarios, SARS-CoV-2 variants were found to be associated with a mild to severe or even fatal clinical outcome, 15 articles reported such association to be statistically significant. Adjustments in eight of them were made for age, sex and other covariates. Conclusions: SARS-CoV-2 variants can potentially have an impact on clinical outcomes; future studies focused on this topic should consider several covariates that influence the clinical course of the disease.

Nosocomial transmission clusters and lineage diversity characterized by SARS-CoV-2 genomes from two large hospitals in Paris, France, in 2020

France went through three deadly epidemic waves due to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing major public health and socioeconomic issues. We proposed to study the course of the pandemic along 2020 from the outlook of two major Parisian hospitals earliest involved in the fight against COVID-19. Genome sequencing and phylogenetic analysis were performed on samples from patients and health care workers (HCWs) from Bichat (BCB) and Pitié-Salpêtrière (PSL) hospitals. A tree-based phylogenetic clustering method and epidemiological data were used to investigate suspected nosocomial transmission clusters. Clades 20A, 20B and 20C were prevalent during the spring wave and, following summer, clades 20A.EU2 and 20E.EU1 emerged and took over. Phylogenetic clustering identified 57 potential transmission clusters. Epidemiological connections between participants were found for 17 of these, with a higher proportion of HCWs. The joint presence of HCWs and patients suggest viral contaminations between these two groups. We provide an enhanced overview of SARS-CoV-2 phylogenetic changes over 2020 in the Paris area, one of the regions with highest incidence in France. Despite the low genetic diversity displayed by the SARS-CoV-2, we showed that phylogenetic analysis, along with comprehensive epidemiological data, helps to identify and investigate healthcare associated clusters.

Emerging SARS-CoV-2 Genotypes Show Different Replication Patterns in Human Pulmonary and Intestinal Epithelial Cells

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) quickly spread worldwide following its emergence in Wuhan, China, and hit pandemic levels. Its tremendous incidence favoured the emergence of viral variants. The current genome diversity of SARS-CoV-2 has a clear impact on epidemiology and clinical practice, especially regarding transmission rates and the effectiveness of vaccines. In this study, we evaluated the replication of different SARS-CoV-2 isolates representing different virus genotypes which have been isolated throughout the pandemic. We used three distinct cell lines, including Vero E6 cells originating from monkeys; Caco-2 cells, an intestinal epithelium cell line originating from humans; and Calu-3 cells, a pulmonary epithelium cell line also originating from humans. We used RT-qPCR to replicate different SARS-CoV-2 genotypes by quantifying the virus released in the culture supernatant of infected cells. We found that the different viral isolates replicate similarly in Caco-2 cells, but show very different replicative capacities in Calu-3 cells. This was especially highlighted for the lineages B.1.1.7, B.1.351 and P.1, which are considered to be variants of concern. These results underscore the importance of the evaluation and characterisation of each SARS-CoV-2 isolate in order to establish the replication patterns before performing tests, and of the consideration of the ideal SARS-CoV-2 genotype-cell type pair for each assay.

Emerging SARS-CoV-2 Variants: Genetic Variability and Clinical Implications

The sudden rise in COVID-19 cases in 2020 and the incessant emergence of fast-spreading variants have created an alarming situation worldwide. Besides the continuous advancements in the design and development of vaccines to combat this deadly pandemic, new variants are frequently reported, possessing mutations that rapidly outcompeted an existing population of circulating variants. As concerns grow about the effects of mutations on the efficacy of vaccines, increased transmissibility, immune escape, and diagnostic failures are few other apprehensions liable for more deadly waves of COVID-19. Although the phenomenon of antigenic drift in new variants of SARS-CoV-2 is still not validated, it is conceived that the virus is acquiring new mutations as a fitness advantage for rapid transmission or to overcome immunological resistance of the host cell. Considerable evolution of SARS-CoV-2 has been observed since its first appearance in 2019, and despite the progress in sequencing efforts to characterize the mutations, their impacts in many variants have not been analyzed. The present article provides a substantial review of literature explaining the emerging variants of SARS-CoV-2 circulating globally, key mutations in viral genome, and the possible impacts of these new mutations on prevention and therapeutic strategies currently administered to combat this pandemic. Rising infections, mortalities, and hospitalizations can possibly be tackled through mass vaccination, social distancing, better management of available healthcare infrastructure, and by prioritizing genome sequencing for better serosurveillance studies and community tracking.

Clinical outcomes in patients infected with different SARS-CoV-2 variants at one hospital during three phases of the COVID-19 epidemic in Marseille, France

OBJECTIVES

To compare the demographics, clinical characteristics and severity of patients infected with nine different SARS-CoV-2 variants, during three phases of the COVID-19 epidemic in Marseille.

METHODS

A single centre retrospective cohort study was conducted in 1760 patients infected with SARS-CoV-2 of Nextstrain clades 20A, 20B, and 20C (first phase, February-May 2020), Pangolin lineages B.1.177 (we named Marseille-2) and B.1.160 (Marseille-4) variants (second phase, June-December 2020), and B.1.1.7 (alpha), B.1.351 (beta), P.1 (gamma) and A.27 (Marseille-501) variants (third phase, January 2021-today). Outcomes were the occurrence of clinical failures, including hospitalisation, transfer to the intensive-care unit, and death.

RESULTS

During each phase, no major differences were observed with regards to age and gender distribution, the prevalence of chronic diseases, and clinical symptoms between variants circulating in a given phase. The B.1.177 and B.1.160 variants were associated with more severe outcomes. Infections occurring during the second phase were associated with a higher rate of death as compared to infections during the first and third phases. Patients in the second phase were more likely to be hospitalised than those in the third phase. Patients infected during the third phase were more frequently obese than others.

CONCLUSION

A large cohort study is recommended to evaluate the transmissibility and to better characterise the clinical severity of emerging variants.

High Individual Heterogeneity of Neutralizing Activities against the Original Strain and Nine Different Variants of SARS-CoV-2

BACKGROUND

Since the beginning of the COVID-19 pandemic, several SARS-CoV-2 variants have sequentially emerged. In France, most cases were due to spike D641G-harbouring viruses that descended initially from the Wuhan strain, then by the variant of B.1.160 lineage we called Marseille-4 since the summer of 2020, which was followed by the Alpha and Beta variants in early 2021, then the Delta variant currently.

METHODS

We determined the neutralising antibody (nAb) titres in sera from convalescent individuals previously infected by these four major local variants and from vaccine recipients to the original Wuhan strain and nine variants, including two recent circulating Delta isolates.

RESULTS

The results show high inter-individual heterogeneity in nAbs, especially according to the variant tested. The major variations among nAbs are based on the genotype responsible for the infection. Patients previously infected with the beta and B.1.160 variants had the lowest nAb titres. We show that this heterogeneity is well explained by spike protein mutants modelling using in silico approaches. The highest titres were observed in individuals vaccinated with the Pfizer/BioNTech COVID-19 vaccine, even against the delta variant.

CONCLUSIONS

Immunity acquired naturally after infection is highly dependent on the infecting variant, and, unexpectedly, mRNA-based vaccine efficacy was shown to be often better than natural immunity in eliciting neutralising antibodies.

Whole Genome Sequencing of SARS-CoV-2 Strains in COVID-19 Patients From Djibouti Shows Novel Mutations and Clades Replacing Over Time

Since the start of COVID-19 pandemic the Republic of Djibouti, in the horn of Africa, has experienced two epidemic waves of the virus between April and August 2020 and between February and May 2021. By May 2021, COVID-19 had affected 1.18% of the Djiboutian population and caused 152 deaths. Djibouti hosts several foreign military bases which makes it a potential hot-spot for the introduction of different SARS-CoV-2 strains. We genotyped fifty three viruses that have spread during the two epidemic waves. Next, using spike sequencing of twenty-eight strains and whole genome sequencing of thirteen strains, we found that Nexstrain clades 20A and 20B with a typically European D614G substitution in the spike and a frequent P2633L substitution in nsp16 were the dominant viruses during the first epidemic wave, while the clade 20H South African variants spread during the second wave characterized by an increase in the number of severe forms of COVID-19.

Spread of Mink SARS-CoV-2 Variants in Humans: A Model of Sarbecovirus Interspecies Evolution

The rapid spread of SARS-CoV-2 variants has quickly spanned doubts and the fear about their ability escape vaccine protection. Some of these variants initially identified in caged were also found in humans. The claim that these variants exhibited lower susceptibility to antibody neutralization led to the slaughter of 17 million minks in Denmark. SARS-CoV-2 prevalence tests led to the discovery of infected farmed minks worldwide. In this study, we revisit the issue of the circulation of SARS-CoV-2 variants in minks as a model of sarbecovirus interspecies evolution by: (1) comparing human and mink angiotensin I converting enzyme 2 (ACE2) and neuropilin 1 (NRP-1) receptors; (2) comparing SARS-CoV-2 sequences from humans and minks; (3) analyzing the impact of mutations on the 3D structure of the spike protein; and (4) predicting linear epitope targets for immune response. Mink-selected SARS-CoV-2 variants carrying the Y453F/D614G mutations display an increased affinity for human ACE2 and can escape neutralization by one monoclonal antibody. However, they are unlikely to lose most of the major epitopes predicted to be targets for neutralizing antibodies. We discuss the consequences of these results for the rational use of SARS-CoV-2 vaccines.

Spreading of a new SARS-CoV-2 N501Y spike variant in a new lineage

OBJECTIVES

Surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomic epidemiology led us to detect several variants since summer 2020. We report the recent spread of a new SARS-CoV-2 spike 501Y variant.

METHODS

SARS-CoV-2 sequences obtained from human nasopharyngeal samples by Illumina next-generation sequencing were analysed using Nextclade and an in-house Python script and were compared using BLASTn to the GISAID database. Phylogeny was investigated using the IQ-TREE software.

RESULTS

We identified that SARS-CoV-2 genomes from four patients diagnosed in our institute harboured a new set of amino acid substitutions including L18F, L452R, N501Y, A653V, H655Y, D796Y, G1219V ± Q677H. These spike N501Y genomes are the first of Nextstrain clade 19B. We obtained partial spike gene sequences of this genotype for an additional 43 patients. All patients infected with this genotype were diagnosed since mid-January 2021. We detected 42 other genomes of this genotype in GISAID, which were obtained from samples collected in December 2020 in four individuals and in 2021 in 38 individuals. The 89 sequences obtained in our institute or other laboratories originated from the Comoros archipelago, western European countries (mostly metropolitan France), Turkey and Nigeria.

CONCLUSION

These findings warrant further studies to investigate the spread, epidemiological and clinical features, and sensitivity to immune responses of this variant.

The next phase of SARS-CoV-2 surveillance: real-time molecular epidemiology

The current coronavirus disease 2019 (COVID-19) pandemic is the first to apply whole-genome sequencing near to real time, with over 2 million severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) whole-genome sequences generated and shared through the GISAID platform. This genomic resource informed public health decision-making throughout the pandemic; it also allowed detection of mutations that might affect virulence, pathogenesis, host range or immune escape as well as the effectiveness of SARS-CoV-2 diagnostics and therapeutics. However, genotype-to-phenotype predictions cannot be performed at the rapid pace of genomic sequencing. To prepare for the next phase of the pandemic, a systematic approach is needed to link global genomic surveillance and timely assessment of the phenotypic characteristics of novel variants, which will support the development and updating of diagnostics, vaccines, therapeutics and nonpharmaceutical interventions. This Review summarizes the current knowledge on key viral mutations and variants and looks to the next phase of surveillance of the evolving pandemic.

A Possible Role of Remdesivir and Plasma Therapy in the Selective Sweep and Emergence of New SARS-CoV-2 Variants

Since summer 2020, SARS-CoV-2 strains at the origin of the COVID-19 pandemic have suddenly been replaced by new SARS-CoV-2 variants, some of which are highly transmissible and spread at a high rate. These variants include the Marseille-4 lineage (Nextclade 20A.EU2) in Europe, the 20I/501Y.V1 variant first detected in the UK, the 20H/501Y.V2 variant first detected in South Africa, and the 20J/501Y.V3 variant first detected in Brazil. These variants are characterized by multiple mutations in the viral spike protein that is targeted by neutralizing antibodies elicited in response to infection or vaccine immunization. The usual coronavirus mutation rate through genetic drift alone cannot account for such rapid changes. Recent reports of the occurrence of such mutations in immunocompromised patients who received remdesivir and/or convalescent plasma or monoclonal antibodies to treat prolonged SARS-CoV-2 infections led us to hypothesize that experimental therapies that fail to cure the patients from COVID-19 could favor the emergence of immune escape SARS-CoV-2 variants. We review here the data that support this hypothesis and urge physicians and clinical trial promoters to systematically monitor viral mutations by whole-genome sequencing for patients who are administered these treatments.

Evaluation of Strategies to Fight COVID-19: The French Paradigm

(1) Background: We collected COVID-19 mortality data and the age distribution of the deceased in France and other European countries, as well as specifically in the cities of Paris and Marseille, and compared them. (2) Methods: Data on mortality related to COVID-19 and the associated age distribution were collected from government institutions in various European countries. In France, data were obtained from INSEE and Santé Publique France. All-cause mortality was also examined in order to study potential excess mortality using EuroMOMO. The Marseille data came from the epidemiological surveillance system. (3) Results: France is one of the European countries most impacted by COVID-19. Its proportion of deaths in people under 60 years of age is higher (6.5%) than that of Italy (4.6%) or Spain (4.7%). Excess mortality (5% more deaths) was also observed. Ile-de-France and the Grand-Est are the two French regions with the highest mortality. The proportion of deaths in the under-60 age group was considerable in Ile-de-France (9.9% vs. 4.5% in the Southern region). There are significantly higher numbers of patients hospitalized, in intensive care and deceased in Paris than in Marseille. (4) Conclusions: No patient management, i.e., from screening to diagnosis, including biological assessment and clinical examination, likely explains the high mortality associated with COVID-19.

Rapid Isothermal Amplification for the Buccal Detection SARS-CoV-2 in the Context of Out-Patient COVID-19 Screening

A commercially available isothermal amplification of SARS-CoV-2 RNA was applied to self-collected saliva samples using dry dental cotton rolls, which were held in the mouth for two minutes. Of 212 tests, isothermal amplification yielded three (0.14%) invalid results, 120 (56.6%) positive results and 89 (42%) negative results. Compared to reference RT-PCR assays routinely performed simultaneously on nasopharyngeal swabs, excluding the three invalid isothermal amplification assays and one RT-PCR invalid assay, these figures indicated that 119/123 (96.7%) samples were positive in both methods and 85/85 samples were negative in both methods. Four positive buccal swabs which were missed by the isothermal amplification, exhibited Ct values of 26-34 in reference RT-PCR assays. Positive isothermal amplification detection was achieved in less than 10 min. Supervision of the self-sampling procedure was key to achieve these performances. These data support the proposal to use the protocol reported in this paper, including supervised buccal self-sampling, to screen people suspected of having COVID-19 at the point of care.

SARS-CoV-2 Infectivity and Severity of COVID-19 According to SARS-CoV-2 Variants: Current Evidence

BACKGROUND

We conducted this review to summarize the relation between viral mutation and infectivity of SARS-CoV-2 and also the severity of COVID-19 in vivo and in vitro.

METHOD

Articles were identified through a literature search until 31 May 2021, in PubMed, Web of Science and Google Scholar.

RESULTS

Sixty-three studies were included. To date, most studies showed that the viral mutations, especially the D614G variant, correlate with a higher infectivity than the wild-type virus. However, the evidence of the association between viral mutation and severity of the disease is scant. A SARS-CoV-2 variant with a 382-nucleotide deletion was associated with less severe infection in patients. The 11,083G > U mutation was significantly associated with asymptomatic patients. By contrast, ORF1ab 4715L and S protein 614G variants were significantly more frequent in patients from countries where high fatality rates were also reported. The current evidence showed that variants of concern have led to increased infectivity and deteriorating epidemiological situations. However, the relation between this variant and severity of COVID-19 infection was contradictory.

CONCLUSION

The COVID-19 pandemic continues to spread worldwide. It is necessary to anticipate large clinical cohorts to evaluate the virulence and transmissibility of SARS-CoV-2 mutants.

The Spike of Concern-The Novel Variants of SARS-CoV-2

The high sequence identity of the first SARS-CoV-2 samples collected in December 2019 at Wuhan did not foretell the emergence of novel variants in the United Kingdom, North and South America, India, or South Africa that drive the current waves of the pandemic. The viral spike receptor possesses two surface areas of high mutagenic plasticity: the supersite in its N-terminal domain (NTD) that is recognised by all anti-NTD antibodies and its receptor binding domain (RBD) where 17 residues make contact with the human Ace2 protein (angiotensin I converting enzyme 2) and many neutralising antibodies bind. While NTD mutations appear at first glance very diverse, they converge on the structure of the supersite. The mutations within the RBD, on the other hand, hone in on only a small number of key sites (K417, L452, E484, N501) that are allosteric control points enabling spike to escape neutralising antibodies while maintaining or even gaining Ace2-binding activity. The D614G mutation is the hallmark of all variants, as it promotes viral spread by increasing the number of open spike protomers in the homo-trimeric receptor complex. This review discusses the recent spike mutations as well as their evolution.

Global Discrepancies between Numbers of Available SARS-CoV-2 Genomes and Human Development Indexes at Country Scales

It has now been over a year since SARS-CoV-2 first emerged in China, in December 2019, and it has spread rapidly around the world. Some variants are currently considered of great concern. We aimed to analyze the numbers of SARS-CoV-2 genome sequences obtained in different countries worldwide until January 2021. On 28 January 2021, we downloaded the deposited genome sequence origin from the GISAID database, and from the "Our world in data" website we downloaded numbers of SARS-CoV-2-diagnosed cases, numbers of SARS-CoV-2-associated deaths, population size, life expectancy, gross domestic product (GDP) per capita, and human development index per country. Files were merged and data were analyzed using Microsoft Excel software. A total of 450,968 SARS-CoV-2 genomes originating from 135 countries on the 5 continents were available. When considering the 19 countries for which the number of genomes per 100 deaths was >100, six were in Europe, while eight were in Asia, three were in Oceania and two were in Africa. Six (30%) of these countries are beyond rank 75, regarding the human development index and four (20%) are beyond rank 80 regarding GDP per capita. Moreover, the comparisons of the number of genomes sequenced per 100 deaths to the human development index by country show that some Western European countries have released similar or lower numbers of genomes than many African or Asian countries with a lower human development index. Previous data highlight great discrepancies between the numbers of available SARS-CoV-2 genomes per 100 cases and deaths and the ranking of countries regarding wealth and development.