The first SARS-CoV-2 genetic variants of concern (VOC) in Poland: The concept of a comprehensive approach to monitoring and surveillance of emerging variants

PurposeWe analyzed the SARS-CoV-2 genome using our integrated genome analysis system and present the concept of a comprehensive approach to monitoring and surveillance of emerging variants.Material/methodsA total of 69 SARS-CoV-2 positive samples (with Ct value ​≤ ​28) were tested. Samples included in this study were selected from 7 areas of eastern Poland. All samples were sequenced on an Illumina MiSeq platform using a 300-cycle MiSeq Reagent Kit v2. BWA was used for reads mapping on the reference SARS-CoV-2 sequence. SAMTools were used for post-processing of reads to genome assembly. Pango lineage and Nexstrain were used to identify variants and amino acid mutations. Statistical analysis was performed with R 4.0.2.ResultsThis study shows the first confirmed case of SARS-CoV-2 in Poland with the lineage B.1.351 (known as 501Y.V2 South African variant), as well as another 18 cases with epidemiologically relevant lineage B.1.1.7, known as British variant. Supplementary analysis of SARS-CoV-2 sequences deposited in GISAID shows that the share of a new variant can change rapidly within one month. In addition, we show a complete, integrated concept of a networked system for analyzing the variability of the SARS-CoV-2 genome, which, used in the present study, generated data and a variant report within 6 days.ConclusionThe analyzed viral genomes showed considerable variability with simultaneous clear distinction of local clusters of genomes showing high similarity. Implementing real-time monitoring of new SARS-CoV-2 variants in Poland is urgently needed, and our developed system is available to be implemented on a large scale.     

Comparative analysis of SARS-CoV-2 Omicron BA.2.12.1 and BA.5.2 variants

The initial severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron subvariants, BA.1 and BA.2, are being progressively displaced by BA.5 in many countries. To provide insight on the replacement of BA.2 by BA.5 as the dominant SARS-CoV-2 variant, we performed a comparative analysis of Omicron BA.2.12.1 and BA.5.2 variants in cell culture and hamster models. We found that BA.5.2 exhibited enhanced replicative kinetics over BA.2.12.1 in vitro and in vivo, which is evidenced by the dominant BA.5.2 viral genome detected at different time points, regardless of immune selection pressure with vaccine-induced serum antibodies. Utilizing reverse genetics, we constructed a mutant SARS-CoV-2 carrying spike F486V substitution, which is an uncharacterized mutation that concurrently discriminates Omicron BA.5.2 from BA.2.12.1 variant. We noticed that the 486th residue does not confer viral replication advantage to the virus. We also found that 486V displayed generally reduced immune evasion capacity when compared with its predecessor, 486F. However, the surge of fitness in BA.5.2 over BA.2.12.1 was not due to stand-alone F486V substitution but as a result of the combination of multiple mutations. Our study upholds the urgency for continuous monitoring of SARS-CoV-2 Omicron variants with enhanced replication fitness.     

Zoonotic spill-over of SARS-CoV-2: mink-adapted virus in humans

ObjectivesThis work aimed to analyse possible zoonotic spill-over of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We report the spill-over of mink-adapted SARS-CoV-2 from farmed mink to humans after adaptation that lasted at least 3 months.MethodsNext-generation sequencing and a bioinformatic approach were applied to analyse the data.ResultsIn an isolate obtained from an asymptomatic patient testing positive for SARS-CoV-2, we found four distinguishing mutations in the S gene that gave rise to the mink-adapted variant (G75V, M177T, Y453F, and C1247F) and others.ConclusionsZoonotic spill-over of SARS-CoV-2 can occur from mink to human.     

Evolutionary history,potential intermediate animal host,and cross-species analyses of SARS-CoV-2

To investigate the evolutionary history of the recent outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in China, a total of 70 genomes of virus strains from China and elsewhere with sampling dates between 24 December 2019 and 3 February 2020 were analyzed. To explore the potential intermediate animal host of the SARS-CoV-2 virus, we reanalyzed virome data sets from pangolins and representative SARS-related coronaviruses isolates from bats, with particular attention paid to the spike glycoprotein gene. We performed phylogenetic, split network, transmission network, likelihood-mapping, and comparative analyses of the genomes. Based on Bayesian time-scaled phylogenetic analysis using the tip-dating method, we estimated the time to the most recent common ancestor and evolutionary rate of SARS-CoV-2, which ranged from 22 to 24 November 2019 and 1.19 to 1.31 × 10⁻³ substitutions per site per year, respectively. Our results also revealed that the BetaCoV/bat/Yunnan/RaTG13/2013 virus was more similar to the SARS-CoV-2 virus than the coronavirus obtained from the two pangolin samples (SRR10168377 and SRR10168378). We also identified a unique peptide (PRRA) insertion in the human SARS-CoV-2 virus, which may be involved in the proteolytic cleavage of the spike protein by cellular proteases, and thus could impact host range and transmissibility. Interestingly, the coronavirus carried by pangolins did not have the RRAR motif. Therefore, we concluded that the human SARS-CoV-2 virus, which is responsible for the recent outbreak of COVID-19, did not come directly from pangolins.     

Early phylogenetic estimate of the effective reproduction number of SARS-CoV-2

To reconstruct the evolutionary dynamics of the 2019 novel-coronavirus recently causing an outbreak in Wuhan, China, 52 SARS-CoV-2 genomes available on 4 February 2020 at Global Initiative on Sharing All Influenza Data were analyzed. The two models used to estimate the reproduction number (coalescent-based exponential growth and a birth-death skyline method) indicated an estimated mean evolutionary rate of 7.8 × 10⁻⁴ subs/site/year (range, 1.1 × 10⁻⁴-15 × 10⁻⁴) and a mean tMRCA of the tree root of 73 days. The estimated R value was 2.6 (range, 2.1-5.1), and increased from 0.8 to 2.4 in December 2019. The estimated mean doubling time of the epidemic was between 3.6 and 4.1 days. This study proves the usefulness of phylogeny in supporting the surveillance of emerging new infections even as the epidemic is growing.     

The first Italian outbreak of SARS-CoV-2 B.1.1.7 lineage in Corzano,Lombardy

In December 2020, Italy experienced the first case of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) B.1.1.7 lineage. In January 2021, we identified 21 cases of this variant in Corzano, defining the first outbreak of SARS-CoV-2 B.1.1.7 lineage in Italy. The high transmissibility of the B.1.1.7 variant represented an important benefit for the virus, which became rapidly dominant on the territory. Containment measures induced the epidemic curve onto a decreasing trajectory underlining the importance of appropriate control and surveillance for restraint of virus spread. Highlights The first Italian outbreak of SARS-CoV-2 B.1.1.7 lineage occurred in Lombardy in January 2021. The outbreak originated by a single introduction of the B.1.1.7 lineage. The genomic sequencing revealed, for the first time, the presence of the V551F mutation in the B.1.1.7 lineage in Italy. Surveillance, prompt sequencing and tracing efforts were fundamental to identify and to quickly contain the outbreak.     

Genomic differences between the diabetogenic and nondiabetogenic variants of encephalomyocarditis virus

Plaque purification of the M variant of encephalomyocarditis (EMC-M) virus resulted in the isolation of two stable variants. One is a highly diabetogenic D variant (EMC-D) and the other is a nondiabetogenic B variant (EMC-B). The cDNA of EMC-D and EMC-B genomes were cloned and seven overlapping cDNA clones were selected to cover the entire genome except the 5'-end 310 bases which were determined by RNA-dependent DNA sequencing and enzymatic RNA sequencing. Each clone was restriction-mapped, subcloned, and sequenced. The genomes of EMC-D and EMC-B are composed of 7829 and 7825 bases, respectively. Both genomes contain a long open reading frame of 6876 nucleotides starting at position 830 on the consensus sequence, which encodes a polyprotein of 2292 amino acids. The sequences of EMC-D and EMC-B differ by two deletions, one insertion, and eight point mutations. The first deletion of 3 nucleotides is located in the 5' poly(C) tract where EMC-B has 127 nucleotides compared with 130 nucleotides in EMC-D. The second deletion in EMC-B involves 2 nucleotides at the 3'-end polyadenylation site. A single base insertion of U occurs at the 5' noncoding region of EMC-B. The eight point mutations are located in the polyprotein coding region. Two are silent and are each located in the structural gene 1B and in the nonstructural gene 2B. The remaining six mutations, one on the L gene and the other five on the 1D gene, introduce respective amino acid changes. It is concluded that the diabetogenic EMC-D viral genome (7829 bases) differs from the nondiabetogenic EMC-B viral genome (7825 bases) by 14 nucleotides out of 7829.     

Comparative structural analyses of selected spike protein-RBD mutations in SARS-CoV-2 lineages

The severity of COVID-19 has been observed throughout the world as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) globally claimed more than 2 million lives and left a devastating impact worldwide. Recently several virulent mutant strains of this virus, such as the B.1.1.7, B.1.351, and P1 lineages, have emerged with initial predominance in UK, South Africa, and Brazil. Another extremely pathogenic B.1.617 lineage and its sub-lineages, first detected in India, are now affecting some countries at notably stronger spread-rates. The present paper computationally examines the time-based structures of B.1.1.7, B.1.351, and P1 lineages with selected spike protein mutations. The mutations in the more recently found B.1.617 lineage and its sub-lineages are explored, and the implications for multiple point mutations of the spike protein’s receptor-binding domain (RBD) are described. The selected S1 mutations within the highly contagious B.1.617.2 sub-lineage, also known as the delta variant, are examined as well.

     

HHV-8 subtypes in South Africa: identification of a case suggesting a novel B variant.

Human herpesvirus 8 (HHV-8) has been identified as the causative agent for all forms of Kaposi's sarcoma and is also associated with the development of body cavity-based B-cell lymphomas and multicentric Castleman's disease. HHV-8 genomes are now classified into five major subtypes (A-E) that reflect sequence heterogeneity in the highly variable open reading frame (ORF) K1. To identify HHV-8 subtypes associated with different forms of Kaposi's sarcoma, we compared the ORF 26 and ORF-K1 gene sequences from South African patients with the prototype strains of the major subtypes, as well as published sequences from other African strains. DNA prepared from Kaposi's sarcoma biopsies and/or peripheral blood lymphocytes were available from 14 patients with postrenal transplant (iatrogenic) Kaposi's sarcoma, six patients with the African endemic form, and one patient with AIDS-related body cavity-based B-cell lymphoma. We identified a B2 subtype in six patients, four of whom also had a novel B5 type ORF 26 polymorphism. Two patients had B2 type patterns for both the ORF 26 and ORF-K1 genes. The ORF-K1 subtype A5 was identified in samples from three patients with a B3/C2 type polymorphism in the ORF 26 gene. A novel ORF-K1 B variant strain was identified in a patient with African endemic Kaposi's sarcoma, who also had a B3/C2 class ORF 26 pattern. In 58.3% of iatrogenic Kaposi's sarcoma patients, a B5-type ORF 26 gene sequence pattern was identified. No association was found among particular subtypes, geographical origin of patients, or clinical presentation.     

Clinical outcomes in COVID-19 patients infected with different SARS-CoV-2 variants in Marseille,France

ObjectivesTo compare the clinical and epidemiological aspects associated with different predominant lineages circulating in Marseille from March 2020 to January 2021.MethodsIn this single-centre retrospective cohort study, characteristics of patients infected with four different severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants were documented from medical files. The outcome was the occurrence of clinical failure, defined as hospitalization (for outpatients), transfer to the intensive care unit (inpatients) and death (all).ResultsA total of 254 patients were infected with clade 20A (20AS), 85 with Marseille-1 (M1V), 190 with Marseille-4 (M4V) and 211 with N501Y (N501YV) variants. 20AS presented a bell-shaped epidemiological curve and nearly disappeared around May 2020. M1V reached a very weak peak, then disappeared after six weeks. M4V appeared in July presented an atypical wave form for 7 months. N501YV has only recently appeared. Compared with 20AS, patients infected with M1V were less likely to report dyspnoea (adjusted odds ratio (OR) 0.50, p 0.04), rhinitis (aOR 0.57, p 0.04) and to be hospitalized (aOR 0.22, p 0.002). Patients infected with M4V were more likely to report fever than those with 20AS and M1V (aOR 2.49, p < 0.0001 and aOR 2.30, p 0.007, respectively) and to be hospitalized than those with M1V (aOR 4.81, p 0.003). Patients infected with N501YV reported lower rate of rhinitis (aOR 0.50, p 0.001) and anosmia (aOR 0.57, p 0.02), compared with those infected with 20AS. A lower rate of hospitalization was associated with N501YV infection compared with 20AS and M4V (aOR 0.33, p < 0.0001 and aOR 0.27, p < 0.0001, respectively).ConclusionsThe four lineages have presentations that differ from one another, epidemiologically and clinically. This supports SARS-CoV-2 genomic surveillance through next-generation sequencing.     

Clinical severity of Omicron subvariants BA.1, BA.2, and BA.5 in a population-based cohort study in British Columbia,Canada

The SARS-CoV-2 variant Omicron emerged in late 2021. In British Columbia (BC), Canada, and globally, three genetically distinct subvariants of Omicron, BA.1, BA.2, and BA.5, emerged and became dominant successively within an 8-month period. SARS-CoV-2 subvariants continue to circulate in the population, acquiring new mutations that have the potential to alter infectivity, immunity, and disease severity. Here, we report a propensity-matched severity analysis from residents of BC over the course of the Omicron wave, including 39,237 individuals infected with BA.1, BA.2, or BA.5 based on paired high-quality sequence data and linked to comprehensive clinical outcomes data between December 23, 2021 and August 31, 2022. Relative to BA.1, BA.2 cases were associated with a 15% and 28% lower risk of hospitalization and intensive care unit (ICU) admission (aHR_hospital = 1.17; 95% confidence interval [CI] = 1.096–1.252; aHR_ICU = 1.368; 95% CI = 1.152–1.624), whereas BA.5 infections were associated with an 18% higher risk of hospitalization (aHR_hospital = 1.18; 95% CI = 1.133–1.224) after accounting for age, sex, comorbidities, vaccination status, geography, and social determinants of health. Phylogenetic analysis revealed no specific subclades associated with more severe clinical outcomes for any Omicron subvariant. In summary, BA.1, BA.2, and BA.5 subvariants were associated with differences in clinical severity, emphasizing how variant-specific monitoring programs remain critical components of patient and population-level public health responses as the pandemic continues.     

A global analysis of replacement of genetic variants of SARS-CoV-2 in association with containment capacity and changes in disease severity

ObjectivesTo examine severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant replacement in association with containment capacity and changes in case fatality at country level.MethodsAltogether, 69 571 full SARS-CoV-2 genomes collected globally within the first 6 months of the pandemic were examined. The correlation between variant replacement and containment capacity was examined by logistic regression models using the WHO International Health Regulation (IHR) score, the Oxford COVID-19 Government Response Tracker (OxCGRT) and the vulnerability index INFORM as proxies, while correlation with changes in monthly crude case fatality ratios was examined by a mixed effect model.ResultsAt the global level, variant lineage G1, characterized by the S-D614G mutation, replaced the older lineages L and S in March 2020. European countries—including Finland, France and Italy—were the first to reach a 50% increment of G1, whereas only Singapore and South Korea had non-G1 persisting throughout the first 6 months. Countries with higher IHR scores (β-coefficient –0.001, 95%CI –0.016, –0.001; p 0.034) and higher stringency indexes (OxCGRT) (β-coefficient –0.011, 95%CI –0.020, –0.001; p 0.035) were associated with lower levels of G1 replacement, whereas higher vulnerability indexes (INFORM) (β-coefficient 0.049, 95%CI 0.001, 0.097; p 0.044) were associated with higher replacement levels. Crude case fatality ratio showed a positive correlation with G1 replacement (β-coefficient: 0.034, 95%CI 0.011, 0.058; p 0.004), even after adjusting for testing capacity and other country-specific characteristics.ConclusionsSARS-CoV-2 variant lineage G1 (S-D614G) replaced older lineages more efficiently in countries with lower containment capacity, and its possible association with increased disease severity deserves further investigation.     

SARS-CoV-2: Structural diversity,phylogeny, and potential animal host identification of spike glycoprotein

To investigate the evolutionary history of the current pandemic outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a total of 137 genomes of coronavirus strains with release dates between January 2019 and 25 March 2020, were analyzed. To investigate the potential intermediate host of the SARS-CoV-2, we analyzed spike glycoprotein sequences from different animals, with particular emphasis on bats. We performed phylogenetic analysis and structural reconstruction of the spike glycoproteins with subsequent alignment and comparison. Our phylogenetic results revealed that SARS-CoV-2 was more similar to the bats' betacoronavirus isolates: HKU5-related from Pipistrellus abramus and HKU4-related from Tylonycteris pachypus. We also identified a yak betacoronavirus strain, YAK/HY24/CH/2017, as the closest match in the comparison of the structural models of spike glycoproteins. Interestingly, a set of unique features has been described for this particular strain of the yak betacoronavirus. Therefore, our results suggest that the human SARS-CoV-2, responsible for the current outbreak of COVID-19, could also come from yak as an intermediate host.     

Cross neutralization of SARS-CoV-2 omicron subvariants after repeated doses of COVID-19 mRNA vaccines

We have measured the humoral response to messenger RNA (mRNA) vaccines in COVID-19 naïve and convalescent individuals. Third doses of mRNA COVID-19 vaccines induced a significant increase in potency and breadth of neutralization against SARS-CoV-2 variants of concern (VoC) including Omicron subvariants BA.1, BA.2, and BA.2.12.1, that were cross-neutralized at comparable levels and less for BA.4/5. This booster effect was especially important in naïve individuals that only after the third dose achieved a level that was comparable with that of vaccinated COVID-19 convalescents except for BA.4/5. Avidity of RBD-binding antibodies was also significantly increased in naïve individuals after the third dose, indicating an association between affinity maturation and cross neutralization of VoC. These results suggest that at least three antigenic stimuli by infection or vaccination with ancestral SARS-CoV-2 sequences are required to induce high avidity cross-neutralizing antibodies. Nevertheless, the circulation of new subvariants such as BA.4/5 with partial resistance to neutralization will have to be closely monitored and eventually consider for future vaccine developments.     

Genetic variants and source of introduction of SARS-CoV-2 in South America

After more than 4 months of the COVID-19 pandemics with genomic information of SARS-CoV-2 around the globe, there are more than 1000 complete genomes of this virus. We used 691 genomes from the GISAID database. Several studies have been reporting mutations and hotspots according to viral evolution. Our work intends to show and compare positions that have variants in 30 complete viral genomes from South American countries. We classified strains according to point alterations and portray the source where strains came into this region. Most viruses entered South America from Europe, followed by Oceania. Only Chilean isolates demonstrated a relationship with Asian isolates. Some changes in South American genomes are near to specific domains related to viral replication or the S protein. Our work contributes to the global understanding of which sort of strains are spreading throughout South America, and the differences among them according to the first isolates introduced to this region.     

SARS-CoV-2 Omicron variant (B.1.1.529): A concern with immune escape

Omicron, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant that is now spreading across the world, is the most altered version to emerge so far, with mutations comparable to changes reported in earlier variants of concern linked with increased transmissibility and partial resistance to vaccine-induced immunity. This article provides an overview of the SARS-CoV-2 variant Omicron (B.1.1.529) by reviewing the literature from major scientific databases. Although clear immunological and clinical data are not yet available, we extrapolated from what is known about mutations present in the Omicron variant of SARS-CoV-2 and offer preliminary indications on transmissibility, severity, and immune escape through existing research and databases.     

Analysis of genome determinants of virulence and setting-up of a library of full-size genome copies of the attenuated virus strain of the porcine reproductive and respiratory syndrome virus North American type

Primary genome structures of 3 variants of the NADC-8 North American virulent strain of porcine reproductive and respiratory syndrome virus (PRRSV) were compared for the purpose of detecting any potential genetic virulence determinants of genus Arterivirus. Apart from the virulent variant, we also investigated the attenuated variant, obtained after 251 passages in cell culture, and the intermediate variant isolated from a pig after a partial reversion of the attenuated virus. The attenuated variant genome acquired a 3-nucleotide deletion and 50 mutations versus its virulent precursor. A comparison of the attenuated and intermediary virus variants denoted 8 nucleotide mutations entailing substitutions of 6 amino acids in 3 open reading frames (ORF1a, ORF1b and ORF6). A 32-clone library was constructed in the pACYC177 plasmid vector, which comprised full-size copies of the genome of the NADC-8 attenuated variant strain (251), virus PPCC, for the purpose of experimentally verifying the functional role of the obtained mutations. Full-size analogues ((+)-chain of RNA) of the viral genome, comprising the CAP-structures and polyadenylated ones were obtained in vitro on the basis of the cloned DNA. Seven of the 8 analyzed clones of the viral genome were infected and their insertion into the MARC-145 cell resulted in obtaining of infectious PRRSVs. Four of the constructed recombinant viruses had delayed growth parameters, and 3 of them were similar to the parental strain. The described technology (inverse genetics) would make it possible to introduce changes into the viral genome in applied and fundamental research of Arteriviruses.     

Genomic characterization and phylogenetic analysis of SARS-COV-2 in Italy

This report describes the isolation, molecular characterization, and phylogenetic analysis of the first three complete genomes of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) isolated from three patients involved in the first outbreak of COVID-19 in Lombardy, Italy. Early molecular epidemiological tracing suggests that SARS-CoV-2 was present in Italy weeks before the first reported cases of infection.     

An Outbreak of Breakthrough Infections by the SARS-CoV-2 Delta Variant in a Psychiatric Closed Ward

BackgroundA rapid decline in immunity and low neutralizing activity against the delta variant in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccinees has been observed. This study describes an outbreak of coronavirus disease 2019 (COVID-19) breakthrough infections caused by the SARS-CoV-2 delta variant in a psychiatric closed ward.MethodsData from epidemic intelligence service officers were utilized to obtain information regarding demographic, vaccination history, and clinical data along with SARS-CoV-2 PCR test results for a COVID-19 outbreak that occurred in a closed psychiatric ward.ResultsAmong the 164 residents, 144 (87.8%) received two doses of vaccines and 137 (95.1%) of them received ChAdOx1 nCoV-19 vaccine. The mean interval between the second vaccination and COVID-19 diagnosis was 132.77 ± 40.68 days. At the time of detection of the index case, SARS-CoV-2 had spread throughout the ward, infecting 162 of 164 residents. The case-fatality ratio was lower than that in the previously reported outbreak before the vaccination (1.2%, 2/162 vs. 6.9%, P = 0.030). Prolonged hospitalization occurred in 17 patients (11.1%) and was less prevalent in the vaccinated group than in the unvaccinated group (8.5% vs. 25.0%, P = 0.040).ConclusionThe findings of this study highlight that while vaccination can reduce mortality and the duration of hospitalization, it is not sufficient to prevent an outbreak of the SARS-CoV-2 delta variant in the present psychiatric hospital setting.