Next Article in Journal
Genotypic and Phenotypic Characteristics of Moraxella catarrhalis from Patients and Healthy Asymptomatic Participants among Preschool Children
Previous Article in Journal
Increased Levels of C5a in Gingival Crevicular Fluid and Saliva of Patients with Periodontal Disease
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pathogens
Volume 11
Issue 9
10.3390/pathogens11090985
pathogens-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Strain Variation Based on Spike Glycoprotein Gene of SARS-CoV-2 in Kuwait from 2020 to 2021

by
Nada Madi
1,*,
Mohammad Sadeq
2,
Sahar Essa
1,
Hussain A. Safar
3,
Anfal Al-Adwani
1 and
Marwa Al-Khabbaz
1
1
Virology Unit, Department of Microbiology, Faculty of Medicine, Kuwait University, Safat 13110, Kuwait
2
Jaber Al-Ahmad Armed Forces Hospital, Sabhan 91710, Kuwait
3
Research Core Facility and OMICS Research Unit, Faculty of Medicine, Kuwait University, Safat 13110, Kuwait
*
Author to whom correspondence should be addressed.
Pathogens 2022, 11(9), 985; https://doi.org/10.3390/pathogens11090985
Submission received: 31 July 2022 / Revised: 22 August 2022 / Accepted: 26 August 2022 / Published: 29 August 2022
(This article belongs to the Section Viral Pathogens)
Browse Figures
Versions Notes

Abstract

:
Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is the causative agent of coronavirus disease 2019 (COVID-19), which was first identified in Wuhan, China, in December 2019. With the global transmission of the virus, many SARS-CoV-2 variants have emerged due to the alterations of the spike glycoprotein. Therefore, the S glycoprotein encoding gene has widely been used for the molecular analysis of SARS-Co-2 due to its features affecting antigenicity and immunogenicity. We analyzed the S gene sequences of 35 SARS-CoV-2 isolates in Kuwait from March 2020 to February 2021 using the Sanger method and MinION nanopore technology to confirm novel nucleotide alterations. Our results show that the Kuwaiti strains from clade 19A and B were the dominant variants early in the pandemic, while clade 20I (Alpha, V1) was the dominant variant from February 2021 onward. Besides the known mutations, 21 nucleotide deletions in the S glycoprotein in one Kuwaiti strain were detected, which might reveal a recombinant SARS-CoV-2 with the defective viral genome (DVG). This study emphasizes the importance of closely perceiving the emerging clades with these mutations during this continuous pandemic as some may influence the specificity of diagnostic tests, such as RT-PCR and even vaccine design directing these positions.

1. Introduction

The novel severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), the causative agent of coronavirus disease (COVID-19), has become a global threat since its first identification in Wuhan City, China, in December 2020 [1]. At the time of writing, the death toll has surpassed 6 million among 596,020,485 confirmed cases around the Globe (until 21 August 2022) [2].
According to the phylogenetic analysis, SARS-CoV-2 is a member of the genus Betacoronavirus, including the SARS coronavirus (SARS-CoV) [3,4]. Studies have shown that SARS-CoV-2 has a close relationship with bat SARS-like CoVs [5]; however, the intermediate host for the transmission to humans remains uncertain. Many progenies of SARS-CoV-2 have already been produced from the original Wuhan strain due to a high number of mutations in the genome, thus escaping the host immune responses [6,7,8,9,10]. SARS-CoV-2 is an enveloped, single-stranded positive-sense RNA with a genome of approximately 29.9 kb encoding 27 proteins from 14 ORFs; 15 non-structural, eight accessory, and four major structural proteins [4]. The four major structural proteins—spike (S) protein, nucleoprotein (N), membrane (M) protein, and envelope (E) protein—are required for the completion of successful infection events of viral entry, assembly, packaging, and release of new viral particles inside human cells [11,12,13].
There is evidence that the spike (S) protein, a 1273 amino acid long glycoprotein with multiple domains, plays a significant role in SARS-CoV-2 pathogenesis. Like other coronaviruses, the S protein of SARS-CoV-2 facilitates receptor recognition, cell attachment, and fusion during viral infection [14]. It binds to the angiotensin-converting enzyme 2 (ACE2) receptor [15]. The S protein consists of two functional subunits, S1 and S2, with a critical rule in viral entry. The S1 subunit recognizes and binds to the host receptor through the receptor-binding domain (RBD), while S2 is responsible for fusion with the host cell membrane [16,17]. Moreover, recent studies have shown that neutralizing antibodies are generated in response to the entry and fusion of surface-exposed S protein (mainly RBD), making it an important target for vaccine development [18,19]. Several variants of SARS-CoV-2 have emerged depending on amino acid mutations in S gene sequences, and molecular analysis based on the S gene can deliver visions into antigenicity, immunogenicity, or evolutionary trends [20]. Therefore, in this study, we analyzed the S gene sequences of SARS-CoV-2 to understand this virus better. Furthermore, detecting and understanding the sequence variations of S glycoprotein among the circulating viruses is essential for better understanding the biology of SARS-CoV-2 infection, developing antiviral treatments, and designing potent vaccines.
In this study, we analyzed spike glycoprotein gene sequences of SARS-CoV-2 circulating in Kuwait to identify mutations causing amino acid variations in the glycoprotein S. The identification of the diversification in SARS-CoV-2 during the evolution of the virus is essential for the early diagnosis and control of the disease.

2. Results

2.1. Characteristics of the Patients

Between March 2020 and February 2021, 35 samples from SARS-CoV-2 positive patients were collected and analyzed. These samples consisted of residual nasopharyngeal swabs that tested positive for SARS-CoV-2 by RT-PCR. The average age of the enrolled patients was 46 years (range (2–91)), and the gender distribution was 49% (17/35) female and 51% (18/35) male, 54% (19/35) were Kuwaiti among the patients, and 46% (16/35) were non-Kuwaiti. Symptomatic individuals represented 60%, more than half of the study population. According to SARS-CoV-2 [21] guidelines, symptomatic patients enrolled in this study had mild to critical SARS-CoV-2 infection; yet, 14 (40%) of the patients were admitted to ICU, while two (5.7%) died. Details on the clinical conditions and symptoms of the symptomatic patients are presented in Table 1.

2.2. Identification of Mutations and Deletions in the SARS-CoV-2 Spike Glycoprotein Gene

Based on mutation analysis of the SARS-CoV-2 isolates, 32 of 35 isolates from Kuwait contained a signature 23063A>G (D614G) mutation in the spike glycoprotein, which represents the most frequent mutation (91%) (Table 2). The second frequent mutations found at a similar rate were 23271 C>A (A570D), 23604 C>A (P681H), and 24914 G>C (D1118H), found in 40% (14/35) of the isolates. On the other hand, deletions at nt position 21766–21771 (ACATGT); (H69-V70del) were the most frequent deletions, found in 40% (14/35) of the isolates in Kuwait, while the second most frequent deletion, which was found in 37% (13/35) of the isolates, was at nt position 21992–21994 (TAT); (Y144-del). Mutation at nt position 23709 C>T (T716I) and mutation at nt position 24506 T>G (S982A) were found less frequently (37%, 13/35) in the Kuwaiti isolates. In addition, single nucleotide–polymorphism (SNP) at position 23063 A>T (N501Y) was found in 34% (12/35) of the isolates, while SNP at position 22622 A>G (N354D) was the least frequent mutation detected in the isolates, found in only 6% (2/35) of the isolates.
Interestingly, two novel synonymous mutations (22498 C>T [G311G], 22051 G>T [D1163D]) were each found in only one isolate (3%). Additionally, novel and unique deletions at nt position 22331–22351delGGTTGGACAGCTGGTGCTGCA (G257-A263del) represented a defective SARS-CoV-2 stain and were detected in only one Kuwaiti isolate (Table 2). This unique deletion was confirmed with whole-genome sequencing of SARS-CoV-2 using MinION Nanopore technology. Raw data of the sequences, including the whole genome, are provided as Supplementary Materials.

2.3. Genomic Variations of SARS-CoV-2 Based on the S Glycoprotein Gene

To define the circulating clades of SARS-CoV-2 in Kuwait, we applied the Nextstrain nomenclature that includes the original clade 19A and its derivation 19B and the emerging clade 20B and its derivation 20I. The phylogenetic analytical tree of 35 SARS-CoV-2 isolates from Kuwait and the global strains based on the S glycoprotein gene was performed by Nextclade beta of Nextstrain platform, which is based on banded Smith–Waterman alignment with an affine gap-penalty. The evolutional distances showed that all the SARS-CoV-2 isolates circulating in Kuwait between March 2020 and March 2021 belong to three main clades: 19A, 19B, and 20I (Alpha/V1) (Figure 1A, Table 1). Clade 19B constituted 51% (18/35) of the Kuwaiti strains, while 40% (14/35) of the strains belonged to clade 20I (Alpha, V1), which was first identified in the UK in the fall of 2020. The least abundant clade was 19A (9%; 3/35), the original clade of SARS-CoV-2 identified in Wuhan in December 2019. As for the time evolution, the data showed that clades 19A and 19B were dominant from March to December 2020, while from February 2021 onwards, clade 20I was the dominant (Table 1). Figure 1B shows that there is a divergence of 28–30% at the nucleotides level from the original clade 19A (NC_045512) and the Kuwaiti strains of clade 20I/Alpha/501Y.V1, as well as 7–9% divergence in the nucleotides levels from clade 19A (NC_045512) and the Kuwaiti strains of clade 19B (Figure 1C).
More detailed phylogenetic analysis of the Kuwaiti strains of SARS-CoV-2 based on spike glycoprotein encoded gene sequence grounded on the evolutional tree constructed by MEGA 7 software demonstrated that the Kuwaiti strains are clustered into two main lineages with nucleotide identity 98.8% between all Kuwaiti SARS-CoV-2 strains (Figure 2). Lineage A contained 21/35 (60%) of the strains that are subdivided into three sub-lineage; the first contained 14 Kuwaiti strains (KW-35M, KW-14A, KW-18M, KW-55A, KW-11M, KW-35A, KW-59A, KW-26A, KW-94M, KW-25A, KW-76M, KW-31A, KW-61A, KW-9M) that clustered with the reference strains; MW521759.1/USA (clade 19B), OB999872.1/UK (clade 20E), MW093524.1/USA (clade 20C), and MW155890/Australia (clade 20F) with 0–0.1% divergence among the strains, whereas the second sub-lineage contained four identical Kuwaiti strains (KW-12M, KW-32M, KW-33M, and KW-39 M) with 0–0.1% divergence at the nucleotide level with MN997409.1/USA (clade 19B) and NC_045512.2/Wuhan (clade 19A). The third sub-lineage of lineage A compromised three Kuwaiti strains (KW-2M, KW-4M, KW-51M) that were 100% identical at the nucleotide level to the reference strains MN997409.1/USA (clade 19B) and NC_045512.2/Wuhan (clade 19A) (Figure 2). Lineage B contained 14 Kuwaiti strains subdivided into four sub-lineages that showed 0–0.2% divergence at the nucleotide level with clade 20I/Alpha/501Y.V1.

2.4. Detection of 21 nt Deletion in S Glycoprotein Gene of SARS-CoV-2

One Kuwaiti strain (KW-35A) showed a substantial deletion of 21 nucleotides at positions 22,331–22,351 (GGTTGGACAGCTGGTGCTGCA), which resulted in the deletion of amino acids G257–A263. This deletion is in the S1 subunit/N-terminal domain (14–305) of the spike glycoprotein (Figure 3). To eliminate the possibility that these findings were caused by errors in PCR amplification or Sanger sequencing, whole-genome sequencing of this strain was performed using the Oxford Nanopore Technologies MinION platform. The sequencing results confirmed these deletions, which may suggest the presence of a recombinant SARS-CoV-2 strain with a defective viral genome (DVG).

3. Discussion

This study, to our knowledge, is the first in Kuwait on the SARS-CoV-2 strain variation circulating in Kuwait from March 2020 to February 2021 based on the analysis of spike protein sequences. The first cases of SARS-CoV-2 were documented in February 2020 from 10 travelers arriving from Iran. After that, more cases were imported from other parts of the world, mainly from Europe (Italy), followed by additional importations from other European countries and worldwide until March 2020, when most air traffic was suspended. The analysis of the study population demographics indicates that males averaging 46 years showed higher positivity for SARS-CoV-2. Generally, SARS-CoV-2 infection causes flu-like symptoms; however, the symptoms may become severe and lead to ICU admission, whereas others die from complications [22]. Our study showed that fever, cough, pneumonia, SOB, reduced oxygen saturation, and bilateral infiltration of the lungs through X-ray imaging were frequent events among the admitted symptomatic patients. Patients admitted to the ICU had more severe symptoms that required intubation, and some died from complications of COVID-19 (Table 1). Those tested positive patients at the drive-through and lodged in the quarantine stations were asymptomatic and not followed up; therefore, it was unknown if they developed symptoms afterwards.
To explore SARS-CoV-2 variations, we identified positions in the spike glycoprotein-coding gene that were frequently altered and were known clade-defining positions in addition to novel positions detected for the first time in the Kuwaiti strains. According to the sequence analysis of the spike glycoprotein-encoding gene, the most frequent missense mutation (91%) identified in Kuwaiti strains was at nucleotide position 23,403, resulting in A>G alteration and causing an amino acid change from aspartate to glycine residue at position 614 (D614G). This mutation was frequently found globally (frequency in GISAID; 97.7%) and was either a solo mutation or combined with other mutations in 31 Kuwaiti strains from clades 19B and 20I (Alpha/V1). D614G is located in B-cell epitopes of the S1 subunit that contain numerous major and minor neutralizing antibody epitopes [16,23]. These amino acid changes are believed to alter immunogenicity, which results in the elimination of the B-cell epitopes. B-cells, however, play a significant role in recognizing pathogens and stimulating adaptive immunity in the immune response against virus infection. Therefore, this elimination will likely reduce the immunogenicity by obstructing the immune cell recognition of the virus [24]. Korber et al. showed that this mutation is associated with greater infectivity and higher viral loads but does not appear to cause severe disease [25]. Our results have also shown other missense mutations and deletions that were detected at relatively similar rates and were found as the second most frequent alterations; 23271 C>A (A570D), 23604 C>A (P681H), 24914 G>C (D1118H), 23709 C>T (T716I), 24506 T>G (S982A), 23063 A>T (N501Y), (H69-V70del), and (Y144-del) (Table 2). The least frequent mutations detected were: the missense mutation 22622 A>G (N354D), two novel synonymous mutations 22498 C>T (G311G) and 22051 G>T (D1163D), and the novel 21 nucleotide deletions (G257-A263).
It is well known that any mutation in the S gene coding for one of the structural proteins will affect the attachment of the virus and the transmission of the disease [26]. The previous alterations were distributed along with the spike protein, including the N-terminal domain (NTD) (14–305 residues), receptor-binding domain (RBD, 319–541 residues), heptapeptide repeat sequence 1 (HR1) (912–984 residues), and HR2 (11163–1213 residues) [27,28]. The S protein binds to ACE-2 (angiotensin-converting enzyme 2) via the NTD and RBD region of the S1 subunit, facilitating viral attachment to host cells in the form of a trimer; therefore, alterations in this position might alter the interaction of viral RBD with the host receptor [4,29]. HR1, on the other hand, is found in the S2 domain and is responsible for viral fusion and entry by forming helical bundles with HR2. Mutations in the HR1 region are predominantly responsible for conferring resistance to mouse hepatitis coronaviruses against HR2-derived peptide entry inhibitors [4,30].
Sequencing and analysis of the SARS-CoV-2 spike glycoprotein gene of SARS-CoV-2 strains revealed that different clades were circulated in Kuwait between March 2020 and February 2021, including 19A, 19B, and 20I (Alpha, V1). Of note, clades 19A and 19B emerged in Wuhan and have dominated worldwide, including Kuwait in the early outbreak, while 20I emerged over the summer of 2020 and included two “variants of concern” (VOC) with signature mutations S: N501Y. However, clade 20I quickly became dominant in Kuwait in February 2021, similar to observations across Europe and other parts of the world (Figure 1A, Table 1). Clade 19B, which was dominant in Kuwait from March to December 2020, characterized by the spike D614G amino acid changes, was the only site identified in early March 2020. At that time, the D614G mutation was rare globally and in Kuwait but gained distinction and became dominant afterwards. Before March 1, 2020, it was found in 10% of 997 global sequences; between March 1 and March 31, 2020, it represented 67% of 14,951 sequences, and between April 1 and May 18, 2020, it represented 78% of 12,194 sequences [25], however, until June 2021, it represented 99.7% of the global sequences (according to GISAID data).
According to Nextclade sequence analysis, the 19B variant in Kuwait showed 7–9% divergence at the nucleotide level from the original clade 19A (Figure 1C). Clade 20I (Alpha, V1) was detected as a new variant from clade 20B in Wales, the United Kingdom (UK), in September 2020, and it was detected for the first time in Kuwait in January 2021 in 11 travelers from the UK, which then became the dominant clade in Kuwait until May 2021 when the new variant 20A (Delta) started to emerge (data are not shown). Of note, clade 20I (Alpha, V1) is defined by many genetic changes, with at least 24 mutations, including 14 non-synonymous mutations, four deletions, and six synonymous mutations in ORF1ab, ORF8, nucleocapsid, and spike proteins [31]. In addition, sequence analysis based on the Nextclade platform demonstrated that 20I (Alpha, V1) variant in Kuwait showed 28–30% divergence at the nucleotide level from the original clade 19A (Figure 1B). Much evidence indicates that clade 20I is more efficiently transmitted than other SARS-CoV-2 variants circulating in the UK, and index patients with 20I infections are infected with a higher proportion of contacts than index patients infected with other variants [32,33]. Currently, there are no known differences in clinical outcomes associated with the SARS-CoV-2 variants; however, a higher transmission rate will lead to more detected cases, increasing the number of persons who need clinical care schemes, resulting in more deaths.
Genetic variability among Kuwaiti SARS-CoV-2 strains was also studied. The analysis demonstrated that 16 out of 3822 nucleotide positions were variable, with five (31%) being singleton sites and 11 (69%) positions being parsimony-informative sites. Our data have also shown 98.8% nucleotide identity between all Kuwait strains. The Kuwaiti strains were divided into two main lineages; the first lineage showed 0–0.1% divergence at the nucleotide level from clade 19A and 19B, while the second lineage showed 0–0.2% divergence at the nucleotide level with clade 20I (Alpha).
Remarkably, substantial deletions in one SARS-CoV-2 strain were detected using whole-genome sequencing with the Oxford Nanopore Technologies MinION platform. The strain contained multiple deletions at position 22,331–22,351 in the spike glycoprotein gene while retaining intact 5′ and 3′ genomic untranslated regions (5′ and 3′ UTRs). These deletions may suggest the presence of a recombinant strain with DVG. However, to support this assumption, in vitro assays must be performed. DVGs of SARS-CoV-2 are believed to be amplified by coronavirus replication-transcription complex (RTC) machinery supplied by co-infecting full-length helper coronaviruses [34,35]. Additionally, DVGs in respiratory viruses can act as pathogen-associated molecular patterns (PAMPs) and stimulate the innate immune system [36]. So far, the biological role of DVGs is mainly unknown, though some DVGs interfere with viral replication [37,38]. Therefore, coronaviruses infrequently perform recombination during their replication, producing complex populations of recombined RNA molecules.

4. Materials and Methods

4.1. Study Population

Thirty-five nasopharyngeal swab samples from confirmed cases of SARS-CoV-2 were collected between March 2020 and February 2021 by two hospitals in Kuwait, Al-Adan and Mubarak Al-Kabeer hospitals. Diagnostic purposes for the presence of SARS-CoV-2, based on real-time reverse transcriptase PCR assays (RT-PCR), were followed as per the Ministry of Health guidelines in Kuwait. Swab samples from Al-Adan hospital were collected from patients admitted to the ICU for COVID-19 complications. Swab samples collected from Mubarak Al-Kabeer hospital were from three sources; patients admitted to the hospital’s wards and ICU, individuals screened for SARS-CoV-2 infection in the Sabhan drive-through station for COVID-19 testing, and individuals screened for SARS-CoV-2 infection lodged in quarantine stations. The swab samples were sent to Mubarak Al-Kabeer hospital for SARS-CoV-2 RT-PCR testing.

4.2. S Gene Sequencing and Sequence Analysis

The complete length spike glycoprotein gene from the confirmed cases of SARS-CoV-2 was sequenced and analyzed. For this analysis, the original specimens from the patients were sequenced using the Sanger method. The PCR amplification and DNA sequencing of the entire length of the S gene was achieved using eight pairs of in-house designed primers (Table S1) [39]. The amplification profile comprised an initial denaturation at 95 °C for 5 min, followed by 35 amplification cycles (94 °C for 1 min, 65 °C for 1 min, and 72 °C for 2 min), with a final extension at 72 °C for 10 min. The amplicons were analyzed by electrophoresis in a 1.5% (W/V) agarose gel. The PCR products were purified with a QIAquick gel extraction kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Subsequently, the purified products were sequenced in the forward and the reverse directions by the primers. Sequencing was performed using the ABI 3500/3500× L genetic analyzer (PE Applied Biosystems, Inc., Foster City, CA, USA) with the ABI Prism BigDye Terminator Cycle Sequencing Ready Reaction kit (PE Applied Biosystems, Inc., Foster City, CA, USA). Mutations specific to the Kuwaitis SARS-CoV-2 isolates were identified using multiple analysis tools using the Wuhan SARS-CoV-2 virus as a reference sequence (NC_045512). Sequence variations analysis of Kuwaiti isolates’ S glycoprotein gene sequences and the reference genome was performed using the Clustal W method of MEGA software version 7. Clade assignment, mutation calling, and sequence quality checks were performed using the Nextstrain platform (https://nextstrain.org/, accessed on 21 August 2022) and CoVsurver mutations App-enabled by GISAID (https://www.gisaid.org/, accessed on 21 August 2022). Genome Detective Coronavirus Typing Tool Version 1.13 (https://www.genomedetective.com/app/typingtool/cov/, accessed on 21 August 2022) [40] was also used for variant analysis of SARS-CoV-2 isolates based on S gene sequences. The phylogenic tree was constructed using the Neighbor-Joining method of MEGA 7 based on the Kimura 2-parameter method as a substitution model with uniform rates among sites and complete gap deletion with 1000 replicates. Bootstrap pre-sampling (1000 replications) was used to measure the reliability of individual nodes in the phylogenetic tree. As for comparative phylogenetic analysis, different SARS-CoV-2 variants sequences published in the National Centre for Biotechnology Information (NCBI) database were used. Spike protein gene sequences of the Kuwaiti SARS-CoV-2 are publicly available in GISAID with accession IDs in Table 1.

4.3. SARS-CoV-2 Whole-Genome Sequencing Using MinION Nanopore Technology

The whole genome of SARS-CoV-2 from the sample “KW-35A” was sequenced using the Oxford Nanopore sequencing technology following the ARTIC network protocol [41]. Briefly, viral RNA was reverse-transcripted with LunaScript™ RT SuperMix Kit (New England Biolabs, Ipswich, MA, USA), and the transcribed cDNA was used as the primary input for overlapping tiled PCR reactions using ARTIC Network V3 primer pools with Q5 Hot Start High-Fidelity DNA polymerase (New England Biolabs, Ipswich, MA, USA). The amplicons were cleaned up with AMPure XP beads (Beckman Coulter Diagnostics, CA, USA), and the library was prepared using the ligation sequencing kit (SQK-LSK109) from Oxford Nanopore Technologies (Oxford, UK). Then, the library was quantified using QUBIT 1X dsDNA HS Assay Kit (Invitrogen, Waltham, Massachusetts, United States), and 15 ng of the prepared library was loaded into Oxford Nanopore MinION SpotON Flow Cells FLO-MIN106D, R9.4.1 (Oxford Nanopore Technologies, Oxford, UK). The FastQ files generated by the Mk1C device were used for analysis following the ARTIC Network analysis workflow and EP2ME-lab. The generated consensus file was used to align the sequence of KW-35A with the reference strain (NC_045512) using MEGA software version 7, while sequence variations analysis and mutation detection were performed using the Nextstrain platform and CoVsurver mutations App-enabled by GISAID.

5. Conclusions

The SARS-CoV-2 pandemic has caused a considerable impact on the health and economy worldwide. Thus, understanding the genetic diversity of the virus has become a priority in the fight against the disease. Frequently mutated positions were identified in Kuwaiti strains, most correlated with known clade-defining mutations and observed in sequences worldwide. Most of these mutations were alterations that caused a change in the AA attribute group, which has a bigger chance of affecting the protein. Therefore, it is crucial to closely perceive the emerging clades with these mutations during this continuous pandemic as some may influence the specificity of diagnostic tests such as RT-PCR and even vaccine design directing these positions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens11090985/s1, File S1: S glycoprotein gene sequences of the 35 SARS-CoV-2 Kuwaiti isolates; File S2: Whole genome sequence of KW-35A SARS-CoV-2 isolate; Table S1: S gene primers.

Author Contributions

Conceptualization, N.M.; Methodology, M.S., S.E., H.A.S., A.A.-A. and M.A.-K.; Validation, N.M.; Formal Analysis, N.M.; Data Curation analysis, N.M.; Writing—Original Draft Preparation, N.M.; Writing—Review and Editing, N.M.; Supervision, N.M.; Project Administration, N.M. All authors have read and agreed to the published version of the manuscript.

Funding

This project was unfunded.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki. And the Ministry of Health, Kuwait Ethics Committee approved this study (project no. 1461/2020).

Informed Consent Statement

Since this is a retrospective study of archived samples, the sample names were anonymized, and informed consent from the participants was not obtained and considered unnecessary, according to the Ethics Committee of the Ministry of Health of Kuwait.

Data Availability Statement

Upon request to the corresponding author.

Acknowledgments

We thank all the medical staff at Adan and Mubarak Al-Kabeer hospitals for their outstanding efforts in providing the clinical samples. Great appreciation goes to all laboratory technicians at the two hospitals for devoting time and effort during the SARS-CoV-2 pandemic in Kuwait.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ralph, R.; Lew, J.; Zeng, T.; Francis, M.; Xue, B.; Roux, M.; Ostadgavahi, A.T.; Rubino, S.; Dawe, N.J.; Al-Ahdal, M.N.; et al. 2019-NCoV (Wuhan Virus), a Novel Coronavirus: Human-to-Human Transmission, Travel-Related Cases, and Vaccine Readiness. J. Infect. Dev. Ctries. 2020, 14, 3–17. [Google Scholar] [CrossRef] [PubMed]
  2. COVID-19 Map—Johns Hopkins Coronavirus Resource Center. Available online: https://coronavirus.jhu.edu/map.html (accessed on 26 April 2021).
  3. Gorbalenya, A.E.; Baker, S.C.; Baric, R.S.; de Groot, R.J.; Drosten, C.; Gulyaeva, A.A.; Haagmans, B.L.; Lauber, C.; Leontovich, A.M.; Neuman, B.W.; et al. The Species Severe Acute Respiratory Syndrome-Related Coronavirus: Classifying 2019-NCoV and Naming It SARS-CoV-2. Nat. Microbiol. 2020, 5, 536–544. [Google Scholar]
  4. Mohammad Lokman, S.; Rasheduzzaman, M.; Salauddin, A.; Barua, R.; Yeasmin Tanzina, A.; Hasan Rumi, M.; Imran Hossain, M.; Zonaed Siddiki, A.; Mannan, A.; Mahbub Hasan, M. Exploring the Genomic and Proteomic Variations of SARS-CoV-2 Spike Glycoprotein: A Computational Biology Approach. Infect. Genet. Evol. 2020, 84, 104389. [Google Scholar] [CrossRef]
  5. Zhou, P.; Yang, X.-L.; Wang, X.-G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.-R.; Zhu, Y.; Li, B.; Huang, C.-L.; et al. A Pneumonia Outbreak Associated with a New Coronavirus of Probable Bat Origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef]
  6. DeDiego, M.L.; Pewe, L.; Alvarez, E.; Rejas, M.T.; Perlman, S.; Enjuanes, L. Pathogenicity of Severe Acute Respiratory Coronavirus Deletion Mutants in HACE-2 Transgenic Mice. Virology 2008, 376, 379–389. [Google Scholar] [CrossRef]
  7. Rahman, M.S.; Islam, M.R.; Hoque, M.N.; Alam, A.S.M.R.U.; Akther, M.; Puspo, J.A.; Akter, S.; Anwar, A.; Sultana, M.; Hossain, M.A. Comprehensive Annotations of the Mutational Spectra of SARS-CoV-2 Spike Protein: A Fast and Accurate Pipeline. Transbound. Emerg. Dis. 2020, 68, 1625–1638. [Google Scholar] [CrossRef] [PubMed]
  8. Frampton, D.; Rampling, T.; Cross, A.; Bailey, H.; Heaney, J.; Byott, M.; Scott, R.; Sconza, R.; Price, J.; Margaritis, M.; et al. Genomic Characteristics and Clinical Effect of the Emergent SARS-CoV-2 B.1.1.7 Lineage in London, UK: A Whole-Genome Sequencing and Hospital-Based Cohort Study. Lancet Infect. Dis. 2021, 21, 1246–1256. [Google Scholar] [CrossRef]
  9. Tegally, H.; Wilkinson, E.; Lessells, R.J.; Giandhari, J.; Pillay, S.; Msomi, N.; Mlisana, K.; Bhiman, J.N.; von Gottberg, A.; Walaza, S.; et al. Sixteen Novel Lineages of SARS-CoV-2 in South Africa. Nat. Med. 2021, 27, 440–446. [Google Scholar] [CrossRef]
  10. Voloch, C.M.; da Silva Francisco, R.; de Almeida, L.G.P.; Cardoso, C.C.; Brustolini, O.J.; Gerber, A.L.; de Guimarães, C.A.P.; Mariani, D.; da Costa, R.M.; Ferreira, O.C.; et al. Genomic Characterization of a Novel SARS-CoV-2 Lineage from Rio de Janeiro, Brazil. J. Virol. 2021, 95, e00119–e00121. [Google Scholar] [CrossRef]
  11. Rahman, M.S.; Hoque, M.N.; Islam, M.R.; Akter, S.; Rubayet-Ul-Alam, A.; Siddique, M.A.; Saha, O.; Rahaman, M.M.; Sultana, M.; Crandall, K.A.; et al. Epitope-Based Chimeric Peptide Vaccine Design against S, M and E Proteins of SARS-CoV-2 Etiologic Agent of Global Pandemic COVID-19: An in Silico Approach. PeerJ 2020, 8, e9572. [Google Scholar] [CrossRef]
  12. Rahman, M.S.; Hoque, M.N.; Islam, M.R.; Islam, I.; Mishu, I.D.; Rahaman, M.M.; Sultana, M.; Hossain, M.A. Mutational Insights into the Envelope Protein of SARS-CoV-2. Gene Rep. 2021, 22, 100997. [Google Scholar] [CrossRef] [PubMed]
  13. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.-M.; Wang, W.; Song, Z.-G.; Hu, Y.; Tao, Z.-W.; Tian, J.-H.; Pei, Y.-Y.; et al. A New Coronavirus Associated with Human Respiratory Disease in China. Nature 2020, 579, 265–269. [Google Scholar] [CrossRef] [PubMed]
  14. Huang, Y.; Yang, C.; Xu, X.-F.; Xu, W.; Liu, S.-W. Structural and Functional Properties of SARS-CoV-2 Spike Protein: Potential Antivirus Drug Development for COVID-19. Acta Pharmacol. Sin. 2020, 41, 1141–1149. [Google Scholar] [CrossRef] [PubMed]
  15. Karpiński, T.M.; Ożarowski, M.; Seremak-Mrozikiewicz, A.; Wolski, H.; Wlodkowic, D. The 2020 Race towards SARS-CoV-2 Specific Vaccines. Theranostics 2021, 11, 1690. [Google Scholar] [CrossRef] [PubMed]
  16. Walls, A.C.; Park, Y.-J.; Tortorici, M.A.; Wall, A.; McGuire, A.T.; Veesler, D. Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein. Cell 2020, 181, 281–292.e6. [Google Scholar] [CrossRef]
  17. Wrapp, D.; Wang, N.; Corbett, K.S.; Goldsmith, J.A.; Hsieh, C.-L.; Abiona, O.; Graham, B.S.; McLellan, J.S. Cryo-EM Structure of the 2019-NCoV Spike in the Prefusion Conformation. Science 2020, 367, 1260–1263. [Google Scholar] [CrossRef]
  18. Tian, X.; Li, C.; Huang, A.; Xia, S.; Lu, S.; Shi, Z.; Lu, L.; Jiang, S.; Yang, Z.; Wu, Y.; et al. Potent Binding of 2019 Novel Coronavirus Spike Protein by a SARS Coronavirus-Specific Human Monoclonal Antibody. Emerg. Microbes Infect. 2020, 9, 382–385. [Google Scholar] [CrossRef]
  19. Chen, Y.; Guo, Y.; Pan, Y.; Zhao, Z.J. Structure Analysis of the Receptor Binding of 2019-NCoV. Biochem. Biophys. Res. Commun. 2020, 525, 135–140. [Google Scholar] [CrossRef]
  20. Kim, S.-J.; Nguyen, V.-G.; Park, Y.-H.; Park, B.-K.; Chung, H.-C. A Novel Synonymous Mutation of SARS-CoV-2: Is This Possible to Affect Their Antigenicity and Immunogenicity? Vaccines 2020, 8, 220. [Google Scholar] [CrossRef]
  21. Infection Control: Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2)|CDC. Available online: https://www.cdc.gov/coronavirus/2019-ncov/hcp/infection-control-recommendations.html (accessed on 2 August 2021).
  22. Ortiz-Prado, E.; Simbaña-Rivera, K.; Gómez- Barreno, L.; Rubio-Neira, M.; Guaman, L.P.; Kyriakidis, N.C.; Muslin, C.; Jaramillo, A.M.G.; Barba-Ostria, C.; Cevallos-Robalino, D.; et al. Clinical, Molecular, and Epidemiological Characterization of the SARS-CoV-2 Virus and the Coronavirus Disease 2019 (COVID-19), a Comprehensive Literature Review. Diagn. Microbiol. Infect. Dis. 2020, 98, 115094. [Google Scholar] [CrossRef]
  23. He, Y.; Li, J.; Heck, S.; Lustigman, S.; Jiang, S. Antigenic and Immunogenic Characterization of Recombinant Baculovirus-Expressed Severe Acute Respiratory Syndrome Coronavirus Spike Protein: Implication for Vaccine Design. J. Virol. 2006, 80, 5757–5767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Nagata, S.; Pastan, I. Removal of B Cell Epitopes as a Practical Approach for Reducing the Immunogenicity of Foreign Protein-Based Therapeutics. Adv. Drug Deliv. Rev. 2009, 61, 977–985. [Google Scholar] [CrossRef] [PubMed]
  25. Korber, B.; Fischer, W.M.; Gnanakaran, S.; Yoon, H.; Theiler, J.; Abfalterer, W.; Hengartner, N.; Giorgi, E.E.; Bhattacharya, T.; Foley, B.; et al. Tracking Changes in SARS-CoV-2 Spike: Evidence That D614G Increases Infectivity of the COVID-19 Virus. Cell 2020, 182, 812–827.e19. [Google Scholar] [CrossRef] [PubMed]
  26. Hulswit, R.J.G.; de Haan, C.A.M.; Bosch, B.-J. Coronavirus Spike Protein and Tropism Changes. Adv. Virus Res. 2016, 96, 29–57. [Google Scholar]
  27. Xia, S.; Zhu, Y.; Liu, M.; Lan, Q.; Xu, W.; Wu, Y.; Ying, T.; Liu, S.; Shi, Z.; Jiang, S.; et al. Fusion Mechanism of 2019-NCoV and Fusion Inhibitors Targeting HR1 Domain in Spike Protein. Cell. Mol. Immunol. 2020, 17, 765–767. [Google Scholar] [CrossRef]
  28. Zhu, Y.; Yu, D.; Yan, H.; Chong, H.; He, Y. Design of Potent Membrane Fusion Inhibitors against SARS-CoV-2, an Emerging Coronavirus with High Fusogenic Activity. J. Virol. 2020, 94. [Google Scholar] [CrossRef]
  29. Ortega, J.T.; Serrano, M.L.; Pujol, F.H.; Rangel, H.R. Role of Changes in SARS-CoV-2 Spike Protein in the Interaction with the Human ACE2 Receptor: An in Silico Analysis. EXCLI J. 2020, 19, 410–417. [Google Scholar] [CrossRef]
  30. Cui, J.; Li, F.; Shi, Z.-L. Origin and Evolution of Pathogenic Coronaviruses. Nat. Rev. Microbiol. 2019, 17, 181–192. [Google Scholar] [CrossRef]
  31. Leung, K.; Shum, M.H.H.; Leung, G.M.; Lam, T.T.Y.; Wu, Y.T. Early Transmissibility Assessment of the N501Y Mutant Strains of SARS-CoV-2 in the United Kingdom, October to November 2020. Eurosurveillance 2021, 26, 2002106. [Google Scholar] [CrossRef]
  32. Public Health England. Investigation of Novel SARS-CoV-2 Variant: Variant of Concern 202012/01. Available online: https://www.gov.uk/government/publications/investigation-of-novel-sars-cov-2-variant-variant-of-concern-20201201 (accessed on 21 August 2022).
  33. Volz, E.; Mishra, S.; Chand, M.; Barrett, J.C.; Johnson, R.; Geidelberg, L.; Hinsley, W.R.; Laydon, D.J.; Dabrera, G.; O’Toole, Á.; et al. Assessing transmissibility of SARS-CoV-2 lineage B.1.1.7 in England The COVID-19 Genomics UK (COG-UK) consortium. Nature 2021, 593, 266–269. [Google Scholar] [CrossRef]
  34. Brian, D.A.; Spaan, W.J.M. Recombination and Coronavirus Defective Interfering RNAs. Semin. Virol. 1997, 8, 101–111. [Google Scholar] [CrossRef] [PubMed]
  35. Makino, S.; Fujioka, N.; Fujiwara, K. Structure of the Intracellular Defective Viral RNAs of Defective Interfering Particles of Mouse Hepatitis Virus. J. Virol. 1985, 54, 329–336. [Google Scholar] [CrossRef] [Green Version]
  36. Sun, Y.; Jain, D.; Koziol-White, C.J.; Genoyer, E.; Gilbert, M.; Tapia, K.; Panettieri, R.A.; Hodinka, R.L.; López, C.B. Immunostimulatory Defective Viral Genomes from Respiratory Syncytial Virus Promote a Strong Innate Antiviral Response during Infection in Mice and Humans. PLoS Pathog. 2015, 11, e1005122. [Google Scholar] [CrossRef] [PubMed]
  37. Méndez, A.; Smerdou, C.; Izeta, A.; Gebauer, F.; Enjuanes, L. Molecular Characterization of Transmissible Gastroenteritis Coronavirus Defective Interfering Genomes: Packaging and Heterogeneity. Virology 1996, 217, 495–507. [Google Scholar] [CrossRef]
  38. Pénzes, Z.; Wroe, C.; Brown, T.D.; Britton, P.; Cavanagh, D. Replication and Packaging of Coronavirus Infectious Bronchitis Virus Defective RNAs Lacking a Long Open Reading Frame. J. Virol. 1996, 70, 8660–8668. [Google Scholar] [CrossRef]
  39. Mak, G.C.K.; Lau, A.W.L.; Chan, A.M.Y.; Chan, D.Y.W.; Tsang, D.N.C. The D614G Substitution in the S Gene and Clinical Information for Patients with COVID-19 Detected in Hong Kong. J. Clin. Virol. 2020, 130, 104550. [Google Scholar] [CrossRef]
  40. Cleemput, S.; Dumon, W.; Fonseca, V.; Abdool Karim, W.; Giovanetti, M.; Alcantara, L.C.; Deforche, K.; de Oliveira, T. Genome Detective Coronavirus Typing Tool for Rapid Identification and Characterization of Novel Coronavirus Genomes. Bioinformatics 2020, 36, 3552–3555. [Google Scholar] [CrossRef]
  41. Pater, A.A.; Bosmeny, M.S.; White, A.A.; Sylvain, R.J.; Eddington, S.B.; Parasrampuria, M.; Ovington, K.N.; Metz, P.E.; Yinusa, A.O.; Barkau, C.L.; et al. High Throughput Nanopore Sequencing of SARS-CoV-2 Viral Genomes from Patient Samples. J. Biol. Methods 2021, 8, e155. [Google Scholar] [CrossRef] [PubMed]
Pathogens 11 00985 g001a 550Pathogens 11 00985 g001b 550
Figure 1. Phylogenetic tree of Kuwaiti sequences based on s gene sequencing compared to the global strains. (A) Maximum clade credibility tree constructed by Nextclade (https://clades.nextstrain.org/ accessed on 21 August 2022) of S gene sequences from Kuwait and all sequences available in GISAID. Global clades are colored differently, and clades corresponding to Kuwaiti strains are circled in red. (B) Divergence of Kuwait strains clade 20I (Alpha, V1) from 19A clade. (C) Divergence of Kuwait strains clade 19B from 19A clade.
Figure 1. Phylogenetic tree of Kuwaiti sequences based on s gene sequencing compared to the global strains. (A) Maximum clade credibility tree constructed by Nextclade (https://clades.nextstrain.org/ accessed on 21 August 2022) of S gene sequences from Kuwait and all sequences available in GISAID. Global clades are colored differently, and clades corresponding to Kuwaiti strains are circled in red. (B) Divergence of Kuwait strains clade 20I (Alpha, V1) from 19A clade. (C) Divergence of Kuwait strains clade 19B from 19A clade.
Pathogens 11 00985 g001aPathogens 11 00985 g001b
Pathogens 11 00985 g002 550
Figure 2. Sequence phylogeny of SARS-CoV-2 spike gene strains in Kuwait. Phylogenic tree representing SARS-CoV-2 circulated in Kuwait and other SARS-CoV-2 variants, retrieved from the NCBI database using the Neighbor-Joining method of MEGA 7 based on the Kimura 2-parameter method substitution model, and tree structure was confirmed by running analysis on 1000 bootstraps. The branch length is indicated in the scale bar. Reference strains are indicated by ●.
Figure 2. Sequence phylogeny of SARS-CoV-2 spike gene strains in Kuwait. Phylogenic tree representing SARS-CoV-2 circulated in Kuwait and other SARS-CoV-2 variants, retrieved from the NCBI database using the Neighbor-Joining method of MEGA 7 based on the Kimura 2-parameter method substitution model, and tree structure was confirmed by running analysis on 1000 bootstraps. The branch length is indicated in the scale bar. Reference strains are indicated by ●.
Pathogens 11 00985 g002
Pathogens 11 00985 g003 550
Figure 3. Deletions identified in the SARS-CoV-2 strain (KW-35A). High-throughput sequencing of the KW-35A strain revealed 21 nt deletions (22331-22351/GGTTGGACAGCTGGTGCTGCA) in the S gene corresponding to aa deletions (G257-A263) in spike glycoprotein. NC-045512.2 was used as a reference strain. Diagram was generated by the Genome Detective Coronavirus Typing Tool Version 1.13 (https://www.genomedetective.com/app/typingtool/cov/ accessed on 21 August 2022).
Figure 3. Deletions identified in the SARS-CoV-2 strain (KW-35A). High-throughput sequencing of the KW-35A strain revealed 21 nt deletions (22331-22351/GGTTGGACAGCTGGTGCTGCA) in the S gene corresponding to aa deletions (G257-A263) in spike glycoprotein. NC-045512.2 was used as a reference strain. Diagram was generated by the Genome Detective Coronavirus Typing Tool Version 1.13 (https://www.genomedetective.com/app/typingtool/cov/ accessed on 21 August 2022).
Pathogens 11 00985 g003
Table
Table 1. Details on the Kuwaiti SARS-CoV-2 strains identified in the study population.
Table 1. Details on the Kuwaiti SARS-CoV-2 strains identified in the study population.
Patient IdentifierStrain NumberCollection DateNationalityAgeSexClinical DataLocationGISAID.
Accession No.		Virus Clade
1KW-61A8 March 2020NK55MFever, cough, pneumonia, CXR: bilateral infiltration.ICU-AH ▲EPI_ISL_327679719B
2KW-76M23 March 2020NK57MFever, cough, loss of taste and smell sensations, SOB, SPO2 80% in right atrium, CXR: B/L infiltrates.ICU-MKH ●EPI_ISL_327686519B
3KW-31A7 April 2020NK72MFever, cough, pneumonia, CXR: bilateral infiltration.ICU-AHEPI_ISL_329829819B
4KW-2M9 May 2020K77MSevere productive cough, SOB, wheeze, SPO2 90% in right atrium, generalized fatigue and tiredness, sore throat, CXR: B/L infiltrates.ICU-MKHEPI_ISL_329829919A
5KW-4M10 May 2020K62FDM type 2, progressive SOB, SPO2 88% in right atrium, high BP, CXR: B/L lung infiltrates and haziness >50% more on the right side.ICU-MKHEPI_ISL_329830119A
6KW-94M18 May 2020K91FReduced level of consciousness, generalized fatigue, vomiting, intubated. Ketoacidosis. CXR: fair entry B/L. Patient deceased.ICU-MKHEPI_ISL_329830219B
7KW-9M18 May 2020NK62FCOPD, smoker, dry cough, fever, SPO2 72% on arrival, CXR: B/L infiltrates and pulmonary congestion.ICU-MKHEPI_ISL_329831119B
8KW-18M26 May 2020NK69MPara-umbilical hernia, fever, cough, SOB, SPO2 85% in right atraium, blood acidosis. CXR: diffuse B/L infiltrates more on the left side. Put on mechanical ventilation—patient deceased.ICU-MKHEPI_ISL_329837419B
9KW-26A23 June 2020NK.49MFever, pneumonia, CXR: bilateral infiltration/Retroperotonial hematoma.ICU-AHEPI_ISL_329837519B
10KW-25A24 June 2020NK.59MFever, cough, pneumonia, CXR: bilateral infiltration.ICU-AHEPI_ISL_329837619B
11KW-11M1 July 2020K43FFever, SOB, noisy chest pain, cough (progressive, green sputum), SPO2 >95% right atrium, chest: B/L wheeze. CXR: B/L haziness.Ward 28-MKH ♣EPI_ISL_329837719B
12KW-55A27 July 2020NK55FFever, cough, pneumonia, CXR: bilateral infiltration.ICU-AHEPI_ISL_329837919B
13KW-12M28 July 2020K34FExacerbation of B.A., headache, fever, low BP. Patient on bronchodilators, Chest: Normal. CXR: Normal.Ward 28-MKHEPI_ISL_329838119B
14KW-59A31 July 2020K70FFever, pneumonia, CXR: bilateral infiltration.ICU-AHEPI_ISL_329868619B
15KW-33M8 August 2020K84MCOPD, HTN, ex-smoker. 9 days H/O SOB, dry cough, 2 days fever, SPO2 66% on arrival —> 98%, CXR: B/L infiltrates and pulmonary congestion.Ward 28-MKHEPI_ISL_329868819B
16KW-32M11 August 2020NK42MDM, HTN on treatment, headache, blurred vision, abdominal pain. Chest: Normal. CXR: Normal.Ward 28-MKHEPI_ISL_329869119B
17KW-39M11 August 2020NK37FFever, 3 days cough, 1 day SOB, pleuritic chest pain, SPO2 87% in right atrium, CXR: B/L infiltrates.Ward 28-MKHEPI_ISL_329876319B
18KW-35M19 August 2020K36FDry cough, SOB, chills, runny nose, no fever. CXR: B/L scattered infiltrates. Chest exam: B/L scattered expiratory wheeze.Ward 28-MKHEPI_ISL_329876519B
19KW-51M5 October 2020K62FUncontrolled DM type 2, HTN, fever, cough productive of sputum, pleuritic chest pain, SOB, SPO 88% in right atrium, loss of taste. Chest: reduced AE B/L. CXR: B/L haziness all over. CT of chest: lung fibrotic-like findings.ICU-MKHEPI_ISL_329876719A
20KW-35A7 October 2020NK31MCough, pneumonia, CXR: bill reticular/miliary opacities.ICU-AHEPI_ISL_329876919B
21KW-14A7 December 2020NK65MFever, cough, pneumonia, SOB, desaturating.ICU-AHEPI_ISL_330195019B
22KW-2V2 February 2021K24FNoneSabhan-Drive through ■EPI_ISL_330224520I (Alpha, V1)
23KW-4V3 February 2021K25FNoneQuarantine lodgeEPI_ISL_330248520I (Alpha, V1)
24KW-6V3 February 2021K26FNoneQuarantine lodgeEPI_ISL_330270020I (Alpha, V1)
25KW-7V4 February 2021NK.58MNoneQuarantine lodgeEPI_ISL_330293520I (Alpha, V1)
26KW-8V4 February 2021K35FNoneSabhan-Drive throughEPI_ISL_330315720I (Alpha, V1)
27KW-9V4 February 2021K34MNoneSabhan-Drive throughEPI_ISL_330338020I (Alpha, V1)
28KW-11V4 February 2021NK38FNoneSabhan-Drive throughEPI_ISL_330363120I (Alpha, V1)
29KW-12V4 February 2021K8MNoneSabhan-Drive throughEPI_ISL_330372620I (Alpha, V1)
30KW-13V4 February 2021K33FNoneSabhan-Drive throughEPI_ISL_330372720I (Alpha, V1)
31KW-15V4 February 2021K16MNoneSabhan-Drive throughEPI_ISL_330372920I (Alpha, V1)
32KW-10V4 February 2021K2MNoneSabhan-Drive throughEPI_ISL_330393220I (Alpha, V1)
33KW-17V5 February 2021NK30MNoneQuarantine lodgeEPI_ISL_330393320I (Alpha, V1)
34KW-19V5 February 2021K38MNoneSabhan-Drive throughEPI_ISL_330393420I (Alpha, V1)
35KW-20V5 February 2021NK30FNoneSabhan-Drive throughEPI_ISL_330393520I (Alpha, V1)
▲ Intensive Care Unit-Al-Adan hospital; ● Intensive Care Unit-Mubarak Al-Kabeer hospital; ♣ Ward 28-Mubarak Al-Kabeer Hospital; ■ Sabhan Drive-through station for COVID-19 testing; NK: Non-Kuwaiti; K: Kuwaiti; SOB: Shortness of breath; DM: Diabetes Mellitus; SPO2: Oxygen saturation; BP: blood pressure; CXR: Chest x-ray; BA: Bronchial Asthma; HTN: Hypertension; COPD: Chronic obstructive pulmonary disease.
Table
Table 2. Variants frequencies of 35 Kuwaiti SARS-CoV-2 isolates.
Table 2. Variants frequencies of 35 Kuwaiti SARS-CoV-2 isolates.
Genome ChangeAmino Acid ChangeType of MutationNumber of IsolatesFrequency in GISAID (%)
23403 A>GD614GMissense3297.7
23271 C>AA570DMissense1445
23604 C>AP681HMissense1447.1
24914 G>CD1118HMissense1444.1
21766-21771delACATGTH69-V70delIn frame deletion1445.5, 44.8
23709 C>TT716IMissense1344.8
24506 T>GS982AMissense1344
21992-21994delTATY144-delIn frame deletion1344.1
23063 A>TN501YMissense1247.5
22622 A>GN354DMissense20.01
22498 C>TG311GSynonymous10
22051 G>TD1163DSynonymous10
22331-22351delGGTTGGACAGCTGGTGCTGCAG257-A263delIn frame deletion10
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Madi, N.; Sadeq, M.; Essa, S.; Safar, H.A.; Al-Adwani, A.; Al-Khabbaz, M. Strain Variation Based on Spike Glycoprotein Gene of SARS-CoV-2 in Kuwait from 2020 to 2021. Pathogens 2022, 11, 985. https://doi.org/10.3390/pathogens11090985

AMA Style

Madi N, Sadeq M, Essa S, Safar HA, Al-Adwani A, Al-Khabbaz M. Strain Variation Based on Spike Glycoprotein Gene of SARS-CoV-2 in Kuwait from 2020 to 2021. Pathogens. 2022; 11(9):985. https://doi.org/10.3390/pathogens11090985

Chicago/Turabian Style

Madi, Nada, Mohammad Sadeq, Sahar Essa, Hussain A. Safar, Anfal Al-Adwani, and Marwa Al-Khabbaz. 2022. "Strain Variation Based on Spike Glycoprotein Gene of SARS-CoV-2 in Kuwait from 2020 to 2021" Pathogens 11, no. 9: 985. https://doi.org/10.3390/pathogens11090985

APA Style

Madi, N., Sadeq, M., Essa, S., Safar, H. A., Al-Adwani, A., & Al-Khabbaz, M. (2022). Strain Variation Based on Spike Glycoprotein Gene of SARS-CoV-2 in Kuwait from 2020 to 2021. Pathogens, 11(9), 985. https://doi.org/10.3390/pathogens11090985

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop