Next Article in Journal
Biomechanical Behavior Characterization and Constitutive Models of Porcine Trabecular Tibiae
Previous Article in Journal
Exploiting DNA Endonucleases to Advance Mechanisms of DNA Repair
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 10
Issue 6
10.3390/biology10060531
biology-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

SARS-CoV-2 Variants in Lebanon: Evolution and Current Situation

by
Nancy Fayad
1,
Walid Abi Habib
1,
Ahmed Kandeil
2,
Rabeh El-Shesheny
2,3,
Mina Nabil Kamel
2,
Youmna Mourad
4,
Jacques Mokhbat
5,
Ghazi Kayali
6,7,
Jimi Goldstein
6,8,* and
Jad Abdallah
1,*
1
School of Pharmacy, Lebanese American University, P.O. Box 36, Byblos, Lebanon
2
Center of Scientific Excellence for Influenza Viruses, National Research Centre, Giza 12622, Egypt
3
St. Jude Children’s Research Hospital, 262 Danny Thomas Place, Memphis, TN 38105, USA
4
Al Hadi Laboratory and IVF Center, P.O. Box 44, Beirut, Lebanon
5
School of Medicine, Lebanese American University, P.O. Box 36, Byblos, Lebanon
6
Human Link, Dubai 971, United Arab Emirates
7
Department of Epidemiology, Human Genetics, and Environmental Sciences, University of Texas, Houston, TX 77030, USA
8
School of Engineering and Technology, University of Hertfordshire, Hatfield, Hertfordshire AL10 9AB, UK
*
Authors to whom correspondence should be addressed.
Biology 2021, 10(6), 531; https://doi.org/10.3390/biology10060531
Submission received: 15 May 2021 / Revised: 30 May 2021 / Accepted: 8 June 2021 / Published: 14 June 2021
(This article belongs to the Section Infection Biology)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Thanks to whole genome sequencing approaches, the proteins that compose the SARS-CoV-2 virus that causes COVID-19 are compared, and mutations and differences highlighted. This allows strains to be classified into clades or lineages, according to marker mutations in their proteins as well as relatedness to certain epidemiological events. In this study, 58 SARS-CoV-2 Lebanese strains were analyzed. They were classified into four GISAID clades and 11 Pango lineages. Moreover, the mutational survey revealed several mutations in the encoded proteins, particularly the structural ones, in which 22 uncommon mutations were found, 21 of which were in strains sequenced within this study. The latter also showed an interesting combination of mutations in their spike protein, the key element in the virus’ interaction with cellular receptors. In summary, this study highlights the key features of sequenced Lebanese SARS-CoV-2 genomes, including their classification, phylogenetic relationship, and mutations.

Abstract

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has seen a worldwide spread since its emergence in 2019, including to Lebanon, where 534,968 confirmed cases (8% of the population) and 7569 deaths have been reported as of 14 May 2021. With the genome sequencing of strains from various countries, several classification systems were established via genome comparison. For instance, the GISAID clades classification highlights key mutations in the encoded proteins that could potentially affect the virus’ infectivity and transmission rates. In this study, 58 genomes of Lebanese SARS-CoV-2 strains were analyzed, 28 of which were sequenced for this study, and 30 retrieved from the GISAID and GenBank databases. We aimed to classify these strains, establish their phylogenetic relationships, and extract the mutations causing amino acid substitutions within, particularly, the structural proteins. The sequenced Lebanese SARS-COV-2 strains were classified into four GISAID clades and 11 Pango lineages. Moreover, 21 uncommon mutations in the structural proteins were found in the newly sequenced strains, underlining interesting combinations of mutations in the spike proteins. Hence, this study constitutes an observation and description of the current SARS-CoV-2 genetic and clade situation in Lebanon according to the available sequenced strains.

Graphical Abstract

1. Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causal agent of coronavirus disease (COVID-19), emerged in the Wuhan district, China, in December 2019 [1]. Since then, this virus has spread around the globe and, on 11 March 2020, was officially declared a global pandemic by the World Health Organization (WHO) [2]. Globally, several strategies were adopted to fight SARS-CoV-2-induced alveolar damages and reduce the symptoms’ severity. A strategy based on combinational drug therapy involving antimalarial, antifilarial, and antiviral drugs was adopted in several countries [3]. Nonetheless, the recent development of several vaccines and the ongoing inoculation of the population constitutes a hope for an end to this pandemic.
SARS-CoV-2 is an enveloped virus, belonging to the Coronaviridae family, with a nonsegmented, positive-sense, single-stranded RNA. Its genome holds 5′ and 3′ untranslated regions (UTR), in addition to ten open reading frames (ORF) encoding structural, nonstructural, and accessory proteins [4]. ORF1ab is proteolytically cleaved into 16 nonstructural proteins (NSP1-16), whereas ORFs 2, 4, 5, and 9 encode the four structural proteins S (spike), E (envelope), M (membrane), and N (nucleocapsid), respectively [5,6]. In general, RNA viruses are prone to a high mutation rate, up to a million times that of their host cells. This property is key to their survival and evolution, where some mutations can affect the virus’ transmission and/or pathogenicity [7].
Similar to other RNA viruses, more than 10,000 single nucleotide polymorphisms (SNPs) were reported for SARS-CoV-2 in only a few months, following the comparison of its genomes [8,9]. This number is increased by (I) the naturally high mutation rate, (II) community transmission, and (III) global transmission. These SNPs and their effect on the encoded proteins lead to the emergence of new SARS-CoV-2 variants and the subsequent establishment of SARS-CoV-2 classification systems. The Global Initiative on Sharing All Influenza Data (GISAID, [10]) introduced a nomenclature structure, which encompasses eight phylogenetic clades based on marker mutations. Clade “L”, which holds the original SARS-CoV-2 strain, Wuhan-Hu-1 (Refseq accession number: NC_045512.2; [11]), forms alongside clade “S”, the two major clades that emerged at the start of the pandemic in 2019. A subsequent division of clade “L” into “V” and “G”, was brought forth by the marker mutation L37F in NS6 and G251V in NS3 for clade “V”, and mutation D614G in S for clade “G”. Then, “G” was split into “GV”, “GH”, and “GR”, due to the mutations NS3 W57H, S A222V, and N G204R, respectively. Clade “GRY” then emerged from clade “GR” in December 2020, because of a combination of ten marker mutations in S. To these eight clades, an “O” clade was added, which is a generic, less clearly defined group [12,13,14,15].
Amino acid changes, particularly in the structural proteins, can affect the infectivity of the strain as well as the severity of the symptoms. For instance, a recent study in Singapore concluded that, on the one hand, strains belonging to clades L and V showed an increased production of proinflammatory cytokines, hence a more pronounced systemic inflammatory reaction in the lungs and upper airways [16]. On the other hand, strains belonging to the clade G or one of its derivates, i.e., presenting the D614G mutation in S, were shown to have higher infectivity and increased viral loads. This is likely due to a change in the conformation of the S protein, resulting in better binding with the host’s receptor protein, the angiotensin-converting 2 (ACE2) protein [17,18,19].
The GISAID clade classification provides a broad overview, in contrast to a more detailed phylogenetic assignment, the Pango lineage, which identifies a cluster according to the sequence relatedness and the association with an epidemiological event [20]. This dynamic nomenclature system uses a phylogenetic framework to classify and track lineages that contribute to active spread, while simultaneously flagging those that are no longer observed and are probably inactive. Two primary lineages, A and B, were assigned based on the divergence of the few first SARS-CoV-2 strains. Then, further specifications were given (e.g., A.2 or B.1) based on ancestor/descendant relationship between strains. When a maximum of three sublevels is reached, a new letter is assigned to indicate that a new level has emerged. In addition to being indicative of the phylogenetic relationships between the strains, this classification allows a reversion to the country of origin [20].
Lebanon confirmed its first case of COVID-19 on 21 February 2020, after a 45-year-old woman returning from Iran tested positive [21]. In the following months, the number of cases increased but remained limited because of strict national confinement, closed borders, and precautionary measures. Then, the country witnessed a surge in the number of confirmed COVID-19 positive cases due to the deteriorating socioeconomic situation and the Beirut port explosion that occurred on 4 August 2020, among other aspects, to reach 534,968 positive cases and 7569 deaths on 14 May 2021 [21,22,23,24,25].
The genomes of 11 Lebanese strains were previously analyzed and showed a number of mutations in the ORF of structural and nonstructural proteins [26,27]. A total of 40 single nucleotide variations were detected in the analyzed strains and resulted in amino acid changes, a stop codon, a frameshift, and a deletion. These variations were found in the ORF1ab, S, NS3 (ORF3a), E, ORF7a, and N proteins. Mutations S D614G and NS3 Q57H were the most frequent among these strains, found in 7 and 6 of the 11 strains, respectively [26,27].
We analyzed the genome of 58 Lebanese circulating strains collected between February 2020 and January 2021, out of which 28 were sequenced for this study and 30 were previously published sequences. Consequently, we generated the phylogenetic relationships between the Lebanese strains and strains from around the world. Moreover, we studied the mutations present in these 58 strains, particularly the mutations located in the four structural proteins S, N, M, and E. This study thus allowed us to highlight the SARS-CoV-2 variants circulating within the Lebanese population.

2. Materials and Methods

2.1. Sample Collection

Nasopharyngeal and oropharyngeal swabs of confirmed COVID-19 patients were collected from the Lebanese American University (LAU) Medical Center-Rizk Hospital and Al-Hadi Medical Center. Twenty-eight samples were collected between 18 March and 1 December 2020 and labeled with the prefix LAU-R (24 strains) and LB-R (4 strains).

2.2. Sample Processing and SARS-CoV-2 Genome Sequencing

The total SARS-CoV-2 nucleic acid was extracted from 140 μL of each swab sample using a QIAmp viral RNA mini kit (cat. No. 52906; QIAGEN, Düsseldorf, Germany) according to the manufacturer’s instructions. Extracted viral RNA was stored at −80 °C until further use. The presence of SARS-CoV-2 was confirmed by real-time PCR using a Verso 1-step RT-qPCR Kit (cat. No. AB4100C; Thermo Fisher Scientific, Massachusetts, USA), primer F-S COV2-NRC (5′-TACCCATTGGTGCAGGTATATGC-3′), and primer R-S COV2-NRC (5′-GTGAGGCAATGATGGATTGACTA-3′) targeting the S and ORF genes in a 40 cycle run. Complementary DNA (cDNA) was generated using the SuperScript™ IV One-Step RT-PCR kit (cat. No. 12594025; Thermo Fisher Scientific, Massachusetts, USA) for genome sequencing, using the following protocol: reverse transcription for 10 min at 50 °C; initial PCR activation for 2 min at 98 °C, followed by standard initial denaturation for 10 sec at 98 °C; annealing for 10 sec at 42 °C; extension for 1 min at 72 °C, with repeating the last three steps 40 times; and then final extension for 5 min at 72 °C. Next generation sequencing (NGS) of the SARS-CoV-2 genome was then performed using a Nextera Illumina protocol as per manufacturer instructions at the Multi-Omics Laboratory at the Lebanese American University (for LAU-R strains) and Saint Jude Children’s Research Hospital, USA, for (LB-R strains).

2.3. Bioinformatic Analysis

The genome sequences of 30 previously published SARS-CoV-2 Lebanese strains were recovered from GISAID and GenBank, in addition to the 28 strains sequenced in this study. Their collection date was from 21 February 2020 to 8 January 2021. Those 58 genomes were assigned a clade and a Pango lineage based on marker mutations, via the built-in GISAID detector and the Pangolin COVID-19 Lineage Assigner (https://pangolin.cog-uk.io/, last accessed on 5 May 2021). Five randomly chosen “reference” strains from each of the eight clades designated by GISAID and the NCBI Refseq Wuhan-Hu-1 strain (Refseq accession number: NC_045512.2) were included in the phylogenetic analysis. In total, 99 sequences were aligned using the MUSCLE algorithm on MegaX (version 10.0.5; [28]). A phylogenetic tree was then constructed via a maximum likelihood approach with a bootstrap value of 500. The resulting tree was rooted with the branch holding the Wuhan-Hu-1 strain.
Furthermore, the mutations present in the 58 SARS-CoV-2 genomes of the Lebanese isolates were analyzed using “CoVsurver: The Mutation Analysis of hCoV-19” (https://www.gisaid.org/epiflu-applications/covsurver-mutations-app/; last accessed on 5 May 2021). The amino acid sequences of the translated ORF were compared to those of the reference strain Wuhan_WIV04 (GISAID accession number: EPI_ISL_402124; GISAID classification Clade: L; Pango Lineage: B). All mutations present in any of the viral protein were examined, particularly those present in the structural proteins S, N, M, and E. Mutations caused by the presence of “N”, hence resulting in an “X” amino acid, were disregarded. Then, detected amino acid substitutions that are considered uncommon according to CoVsurver were cross checked with NCBI’s “Mutations in SARS-CoV-2 SRA Data” (https://www.ncbi.nlm.nih.gov/labs/virus/vssi/#/scov2_snp; last accessed on 5 May 2021), which compares SARS-CoV-2 sequences to the NCBI reference strain Wuhan-Hu-1 (Refseq accession number: NC_045512.2).

3. Results

3.1. SARS-CoV-2 Lebanese Strains: Sequences and Classification

Fifty-eight Lebanese SARS-CoV-2 strains were sequenced and submitted to the GISAID and NCBI GenBank depositories, out of which 28 were sequenced within this study. The genome length of the latter strains ranged from 28,846 to 29,899 bp. As for the other 30 Lebanese strains, their genome was recovered from the GISAID and GenBank depositories, and its length ranged from 29,802 to 29,851 bp. In comparison to the GISAID reference sequence Wuhan_WIV04 (GISAID accession number: EPI_ISL_402124), the insertion or deletion of one or more nucleotides was seen in the genome of several strains (Supplementary Table S1). A frequently occurring deletion, detected in nine strains, consisted of a single nucleotide at position 9850 (position in the reference genome) within the ORF1ab gene, more precisely at the NS4, although no change was detected in direct link to this deletion. In other cases, the insertion of one or more nucleotides was also reported in comparison with Wuhan_WIV04, the longest one being an insertion of 12 nucleotides in strain LB-R12, but that also resulted in no change in the protein sequences (Supplementary Table S1).
Each strain was placed in a GISAID clade, a classification that allows grouping of strains based on marker mutations. For the analyzed strains, clade G included 48% of the strains (28 out of 58 strains). Strains associated with this clade were collected between 19 March 2020 and 1 December 2020. Clades GH and GRY had 21% (12 out of 58) and 17% (10 out of 58), respectively (Figure 1A).
A time course visualization of clade distribution of the 58 analyzed strains showed that clades G and GH overlapped and were the most prevalent in early and mid-2020. However, this was no longer the case in December 2020–January 2021, in which clade GRY took over (Figure 2). This corresponds to what was seen in the time course visualization for the available sequences in the GISAID database ([29]; last updated 5 May 2021).
Moreover, no strains were detected from clades GV, L, V, and S, whereas three strains, S3, S5, and S9, were detected from clade O, a less clearly defined group. Strains S5 and S9 were analyzed in a previous study and identified as belonging to the Nextstrain clade 19A, a group associated with the GISAID clades L or V. As for strain S3, it was excluded in the mentioned study, due to a poor sequence quality [27].
The three most prevalent Pango lineages in the analyzed strains were B.1.398 (45%), B.1 (17%), and B.1.1.7 (17%) (Figure 1B). These lineages are associated with clades G, GH, and GRY, respectively. The Lebanese SARS-CoV-2 strains analyzed in this study and their collection date, GISAID and/or GenBank accession numbers, and assigned GISAID clades and Pango lineages are summarized in Table 1.

3.2. Phylogenetic Relationships

The genomes of the 58 Lebanese SARS-CoV-2 strains and 41 reference strains from each of the eight GISAID clades (Table 2) were aligned and their phylogenetic relationship assessed via a maximum likelihood approach. The resulting tree was rooted with the branch carrying the original SARS-CoV-2 strain Wuhan-Hu-1.
As expected, strains belonging to the same clade clustered together (Figure 3). Nevertheless, a phylogenetic distance was seen in some cases, such as for strains belonging to clade GH, where a small distance existed between strains LAU-R70, -R110, -R113, and -R170 and the rest of this clade’s strains (Figure 3). Another case is that of clade G, where the majority of the strains clustered together, except for the three strains LB-R8, LBN-9-B.1.258 Delta, and Germany BY chVir 929 (GISAID accession number: EPI_ISL_406862). The strains S3, S5, and S9 clustered with the clade V GISAID strains, suggesting that they may in fact belong to this clade, despite their previous classification in clade O, the less clearly defined GSIAID clade. Moreover, seven Lebanese SARS-CoV-2 strains clustered together on a separate branch, within the supercluster encompassing the clade G.

3.3. Mutation Survey: A Focus on Structural Proteins

The classification of the strains and the phylogenetic distance observed in some clades was reflected by the presence of mutations which resulted in a change of the amino acid composition of the encoded proteins. Hence, in order to further investigate the SARS-CoV-2 variants in Lebanon, a mutation survey of the Lebanese SARS-CoV-2 strains showed several amino acid substitutions in both nonstructural and structural proteins. The most commonly found amino acid mutation, detected in 13 out of the 58 strains, was P323L in NS12, the RNA-dependent RNA polymerase [30]. Another mutation, located in 12 strains, was the Q58H located in the accessory NS3, a transmembrane protein that forms ion channels in the host membrane [31].
As for the structural proteins, the number of amino acid changes within the Lebanese strains ranged from 1 to 14 in comparison to the GISAID reference strain Wuhan_WIV04 (Supplementary Table S1). The nine strains FAS/2021 and LBN_1 to 8 presented the highest number of amino acid changes in the structural proteins, whereas the eight strains LAU-R242, S1, S2 and S7–11 presented only one amino acid change in the structural proteins.

3.3.1. Spike S Protein

For the S protein, 50 amino acid substitutions were detected (Figure 4a, Supplementary Table S1), with D614G being the most prevalent mutation within the analyzed strains (55 out of 58). This change is quite common and was previously reported in 96.73% of all GISAID Spike sequences in 170 countries. Furthermore, T95I was an amino acid change found in 19 of the Lebanese strains and was previously reported in 1.27% of all samples with S sequences in 79 countries. Moreover, T716I combined with nine other changes (S982A, D1118H, H69del, V70del, Y144del, N501Y, A580D, P681H, and D614G) marks the variant SARS-CoV-2 VUI 202012/0, also known as Pango lineage B.1.1.7 or clade GRY, which was first reported in the United Kingdom in December 2020 [32]. This combination of ten mutations was found in nine Lebanese strains collected between 24 December 2020 and 8 January 2021. Additionally, a tenth strain, Leb-UK-B, presented eight out of the ten substitutions, with D1118H and N501Y being absent.
While many of the detected mutations were known to be common, some having one million reports in GISAID at the time of the submission, others were less commonly found according to CoVsurver. Although some surveillance tools show the distribution of an amino acid substitution in accordance with when it was added to the database, the search for a substitution using CoVsurve was done based on the strain’s collection date, not the submission date of its genome. Various amino acid substitutions were detected once in the analyzed sequences and, according to GISAID’s CoVsurver, were either first or only reported in the corresponding strain (Table 3). Concerning the spike protein, a total of 15 amino acid changes were found in eight strains. While some were an “only report” by CoVsurver, this was in relation to the new amino acid, not the position of the change. This is the case of three S protein amino acid substitutions at positions G35, S31, and L5.
LAU-R186, collected on 1 August 2020, was one of the strains in which these uncommon mutations were found. In fact, LAU-R186 had ten amino acid substitutions in its S protein (Table 3; Supplementary Table S1). Three out of the ten changes are of particular interest: (I) N17S removes a potential N-glycosylation site at position 17 and has occurred 797 times in eight countries so far according to the GISAID’s CoVsurver, (II) F59Y has 106 reports until the time of the submission, and (III) Y38Stop can cause a 97.1% truncation of the protein. Another strain with an uncommon combination of amino acid changes is LAU-R45, collected on 1 June 2020. It has eight amino acid changes in its S protein, the most noteworthy being P39T, S31Y, and V11D (Table 3). S31Y is involved in the viral oligomerization interface, according to structural resolution of the S protein (PDB chain: 7cwm_B; data not shown). As for the substitution V11D, an equivalent mutation previously reported in the human MERS-CoV virus suggested an alteration in the antigenic drift [33]. The latter is also affected by a change in amino acid position L5, in which the change L5P was detected in strain LAU-R204, according to the analysis of the equivalent mutation L5F [18].

3.3.2. Nucleocapsid N Protein

In total, 23 different amino acid substitutions were reported for the N protein in 36 of the 58 analyzed genomes, the most prevalent being R03K and G204R, both found in 15 strains (Figure 4b). These substitutions were commonly found in SARS-CoV-2 strains and registered over 500,000 occurrences in the GISAID database. Nevertheless, five amino acid substitutions in the N protein were uncommon (Table 3), four of which, at positions E367, L352, L407, and S318, were previously flagged for changes, but for a different substituting amino acid. The first was found to be involved in an alteration of antigenic drift, [34], whereas the last might be implicated in the viral oligomerization interface, according to the structural resolution of the N protein (PDB chain: 6wzq_D; data not shown). For example, L407F and L407L were previously reported, whereas L407M was labeled as an only report in the Lebanese LAU-R43 strain by GISAID’s CoVsurver and NCBI’s mutations data. Amino acid changes in positions 361–370 were reported in two of the analyzed strains. Moreover, the two stop mutations, E367Stop and L352Stop, respectively found in strains LAU-R145 and LAU-R231, may result in a 12.6 and 16.2% truncation of the protein.

3.3.3. Membrane M and Envelope E proteins

Unlike S and N, the two remaining structural proteins presented only two amino acid substitutions, each detected in the Lebanese SARS-CoV-2 strains. For the M protein, the two amino acid changes were L16I, and A2V, previously reported 18 and 1866 times, respectively. As for the E protein, the two detected amino acid changes, C44S and V62D, were uncommon and detected in LAU-R110 and strain S6, respectively. The former might be implicated in the viral oligomerization interface, according to the structural resolution of the E protein (PDB chain: 5 × 30_A; data not shown).

4. Discussion

Lebanon, a small Mediterranean country with around 7 million inhabitants [35], presented more than 534,968 confirmed COVID-19 cases and 7569 deaths as of 14 May 2021 [25]. The current study forms an observational description of the 58 available genomes of Lebanese SARS-CoV-2 strains, their classification, phylogenetic relationships, and the mutations they hold. The Lebanese SARS-CoV-2 strains included the 28 sequenced for this study and 30 recovered from the GISAID and GenBank databases. Analysis was conducted on phylogenetic and mutational levels: each strain was assigned a GISAID clade and a Pango lineage, a phylogenetic analysis was done between the Lebanese strains and reference strains randomly chosen from the GISAID database, and the amino acid substitutions in the structural proteins were analyzed.
The classification according to the GISAID and Pango lineages revealed a prevalence of the clade G (48% of the strains) and its associated Pango lineage B.1.398 (45%). The main genetic marker of this clade was the amino acid substitution D614G, found in 94.8% (55 out of 58) of the analyzed strains. This substitution was shown to allow for an open conformational state of the S protein, thus enhancing the exposure of the receptor binding domain, which interacts with the ACE2 protein. This site is also the target for neutralizing antibodies. Moreover, this mutation was shown to increase the stability and infectivity rate of the virus in vitro [17]. This contributed to its rapid dissemination and its quick rise to dominant status, following its first appearance in January 2020, at the start of the pandemic [17]. The first confirmed COVID-19 incidence in Lebanon was on 21 February 2020, a time when the D614G mutation was already rapidly spreading, which explains its presence in 94.8% of the analyzed strains.
The clade GH, which differs from G by a marker mutation (NS3/ORF3a Q58H), included 21% of the strains along with its associated Pango lineage B.1 (17%). The clade GRY (18%) and its associated lineage B.1.1.7 included 17% of sequences. This clade was transmitted rapidly in December 2020 first through the UK, then Europe, then to more than 90 countries. It is considered a variant of concern by the Center for Disease Control and Prevention in the USA [36]. The series of amino acid changes in its S protein were behind an increase in transmission rates estimated to be 43 to 90% [37]. The phylogenetic analysis of the strains revealed a clustering mostly according to the GISAID clade and Pango lineage, with several subclades within each superclade, in addition to a number of phylogenetic distances that were noted between strains of the same clade. This might be explained by the subtle changes on the nucleotide level that do not necessarily translate into amino acid substitutions, but are taken into account on the DNA alignment level.
No strain belonged to clades GV, L, V, and S. For the first, it rapidly spread through Europe, mainly Italy, as of July 2020, reaching a peak in November 2020. However, according to the GISAID time course of clades distribution in collected sequences, the GV clade saw a steep decline as of January 2021 and is now being replaced by other clades [29]. As for clades L, V, and S, the first three clades designated at the beginning of the pandemic, detected in the GISAID sequence collection up to June 2020, their absence among the sequences of Lebanese SARS-CoV-2 strains could be explained by the fact that, up until June 2020, the community transmission of this virus in Lebanon was limited due to strict lockdown rules. The increase in the number of confirmed positive cases began in August 2020 [25], which coincides with a prevalence of clades G, GH, and GR, according to the GISAID database [29].
RNA viruses have a high mutation rate, a way to evolve and adapt to new hosts. In this study, we identified 77 amino acid substitutions in the structural proteins S, N, M, and E, 22 of which, located in LAU-R strains sequenced within this study, were labeled as uncommon (Table 3). Four of these substitutions resulted in a stop codon and the truncation of the corresponding protein in the range of 12 to 97%. However, mRNA in RNA viruses may be subject to a noncanonical translation, which includes a stop-codon readthrough, meaning that in a specific context, the resulting protein will be continued and not truncated at the site of the stop-codon [38,39].
The GISAID time course analysis of clade distribution in collected sequences shows that, as of late 2020, the GRY clade is taking over in lieu of G, GR, and GH, which were the dominant clades in mid-2020 [29]. When comparing this time course to the available SARS-CoV-2 Lebanese strains, we see that the G, GH, and GR could be gradually giving way to GRY in December 2020–January 2021.
Knowledge about the SARS-CoV-2 variants circulating among the Lebanese community is imperative and should be taken into account for further public health or epidemiological surveys and decisions. However, this study was limited by the small number of sequences analyzed. Additionally, among those sequences, the majority came from the same geographic area. This small number of sequences is due to the high cost of conducting full genome sequencing in Lebanon and complications associated with the political and economic situation of this country [22]. Hence, this may have affected the conclusions we are drawing in this study.

5. Conclusions

This study is the first to highlight 21 uncommon mutations in the structural proteins found in the newly sequenced LAU-R SARS-CoV-2 strains. Combined with other mutations found in these strains, LAU-R strains present an interesting combination of mutations in the spike protein. Moreover, this study was the first to classify 58 Lebanese SARS-CoV-2 strains in total, according to the GISAID clade classification and Pango lineage, and to compare this classification’s timeline with that of the GISAID sequence collection. Nevertheless, there is a need for a substantial improvement in the sequencing of Lebanese SARS-CoV-2 strains. Hence, future studies should follow a purposive sampling, bearing different sociodemographic statuses, such as a patient’s age, sex, and region of residence in addition to, most importantly, clinical presentation of disease.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biology10060531/s1, Table S1: The 58 analyzed Lebanese SARS-CoV-2 strains. This table summarizes the constellation of AA substitutions they carry as well as noteworthy remarks.

Author Contributions

Conceptualization, G.K., J.A. and J.G.; methodology, G.K., J.A., J.G., N.F. and W.A.H.; software, N.F.; validation, G.K., J.A., J.G. and N.F.; formal analysis, N.F.; investigation, A.K., M.N.K., R.E.-S. and W.A.H.; resources, J.M., W.A.H. and Y.M.; data curation, N.F.; writing—original draft preparation, G.K., J.A. and N.F.; writing—review and editing, G.K., J.A. and N.F.; visualization, N.F.; supervision, G.K., J.A. and J.G.; project administration, J.A. and J.G.; funding acquisition, G.K., J.A., J.M. and W.A.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Institute of Allergy and Infectious Diseases, National Institutes of Health, and US Department of Health and Human Services (under contract HHSN272201400006C), with the support of the National Council for Scientific Research in Lebanon CNRS-L (grant number CNRS-849).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of the Lebanese American University (protocol code LAU.SOP.JA1.15/Apr/2020).

Informed Consent Statement

Not applicable.

Data Availability Statement

SARS-CoV-2 sequences are available on GISAID [40] and NCBI’s GenBank [41].

Acknowledgments

The authors would like to acknowledge Rana Abou Salbi and Rebecca Badra for their administrative support and Hanna Semaan and Omar Bahlawan for their technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Z.; Xiao, K.; Zhang, X.; Roy, A.; Shen, Y. Emergence of SARS-like coronavirus in China: An update. J. Infect. 2020, 80, e28–e29. [Google Scholar] [CrossRef]
  2. World Health Organization. WHO Director General’s Opening Remarks at the Media Briefing on COVID-19. 11 March 2020. Available online: https://www.who.int/director-general/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19---11-march-2020 (accessed on 5 May 2021).
  3. Siddiqui, A.J.; Jahan, S.; Ashraf, S.A.; Alreshidi, M.; Ashraf, M.S.; Patel, M.; Snoussi, M.; Singh, R.; Adnan, M. Current Status and Strategic Possibilities on Potential Use of Combinational Drug Therapy against COVID-19 Caused by SARS-CoV-2. J. Biomol. Struct. Dyn. 2020. [Google Scholar] [CrossRef]
  4. Ciotti, M.; Angeletti, S.; Minieri, M.; Giovannetti, M.; Benvenuto, D.; Pascarella, S.; Sagnelli, C.; Bianchi, M.; Bernardini, S.; Ciccozzi, M. COVID-19 outbreak: An overview. Chemotherapy 2020, 64, 215–223. [Google Scholar] [CrossRef]
  5. Ellis, P.; Somogyvári, F.; Virok, D.P.; Noseda, M.; Mclean, G.R. Decoding COVID-19 with the SARS-CoV-2 genome. Curr. Genet. Med. Rep. 2021, 9, 1–12. [Google Scholar] [CrossRef]
  6. Gültekin, V.; Allmer, J. Novel perspectives for SARS-CoV-2 genome browsing. J. Integr. Bioinform. 2021, 18, 19–26. [Google Scholar] [CrossRef]
  7. Duffy, S. Why are RNA virus mutation rates so damn high? PLoS Biol. 2018, 16, e3000003. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, J.; Wang, R.; Wang, M.; Wei, G.W. Mutations Strengthened SARS-CoV-2 Infectivity. J. Mol. Biol. 2020, 432, 5212–5226. [Google Scholar] [CrossRef]
  9. Wang, R.; Hozumi, Y.; Yin, C.; Wei, G.W. Decoding SARS-CoV-2 Transmission and Evolution and Ramifications for COVID-19 Diagnosis, Vaccine, and Medicine. J. Chem. Inf. Model. 2020, 60, 5853–5865. [Google Scholar] [CrossRef]
  10. Shu, Y.; Mccauley, J. GISAID: Global initiative on sharing all influenza data-from vision to reality. Euro Surveill. 2017, 22, 30494. [Google Scholar] [CrossRef] [Green Version]
  11. 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] [Green Version]
  12. Alm, E.; Broberg, E.K.; Connor, T.; Hodcroft, E.B.; Komissarov, A.B.; Maurer-Stroh, S.; Melidou, A.; Neher, R.A.; O’Toole, Á.; Pereyaslov, D.; et al. Geographical and temporal distribution of SARS-CoV-2 clades in the WHO European Region, January to June 2020. Euro Surveill. 2020, 25, 1–8. [Google Scholar] [CrossRef]
  13. Mercatelli, D.; Giorgi, F.M. Geographic and genomic distribution of SARS-CoV-2 mutations. Front. Microbiol. 2020, 11, 1–13. [Google Scholar] [CrossRef]
  14. Sengupta, A.; Sarif, S.; Pal, P. Clade GR and clade GH isolates of SARS-CoV-2 in Asia show highest amount of SNPs. Infect. Genet. Evol. 2020, 89, 104724. [Google Scholar] [CrossRef] [PubMed]
  15. Hodcroft, E.B.; Zuber, M.; Nadeau, S.; Crawford, K.H.D.; Bloom, J.D.; Veesler, D.; Vaughan, T.G.; Comas, I.; Candelas, F.G.; Stadler, T.; et al. Emergence and spread of a SARS-CoV-2 variant through Europe in the summer of 2020. MedRxiv 2020. [Google Scholar] [CrossRef]
  16. Young, B.E.; Wei, W.E.; Fong, S.-W.; Mak, T.-M.; Anderson, D.E.; Chan, Y.-H.; Pung, R.; Heng, C.S.; Ang, L.W.; Zheng, A.K.E.; et al. Association of SARS-CoV-2 clades with clinical, inflammatory and virologic outcomes: An observational study. EBioMedicine 2021, 66, 103319. [Google Scholar] [CrossRef] [PubMed]
  17. Mansbach, R.A.; Chakraborty, S.; Nguyen, K.; Montefiori, D.C.; Korber, B.; Gnanakaran, G. The SARS-CoV-2 spike variant D614G favors an open conformational state. Biophys. J. 2021, 120, 298a. [Google Scholar] [CrossRef]
  18. Li, Q.; Wu, J.; Nie, J.; Zhang, L.; Hao, H.; Liu, S.; Zhao, C.; Zhang, Q.; Liu, H.; Nie, L.; et al. The impact of mutations in SARS-CoV-2 spike on viral infectivity and antigenicity. Cell 2020, 182, 1284–1294. [Google Scholar] [CrossRef] [PubMed]
  19. 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. [Google Scholar] [CrossRef]
  20. Rambaut, A.; Holmes, E.C.; O’Toole, Á.; Hill, V.; McCrone, J.T.; Ruis, C.; du Plessis, L.; Pybus, O.G. A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. Nat. Microbiol. 2020, 5, 1403–1407. [Google Scholar] [CrossRef]
  21. Lebanese Ministry of Public Health. Novel Coronavirus 2019-Lebanon. Available online: https://www.moph.gov.lb/en/Pages/2/24870/novel-coronavirus-2019 (accessed on 6 May 2021).
  22. Bizri, A.R.; Khachfe, H.H.; Fares, M.Y.; Musharrafieh, U. COVID-19 pandemic: An insult over injury for Lebanon. J. Community Health 2020. [Google Scholar] [CrossRef]
  23. El Deeb, O.; Jalloul, M. The dynamics of COVID-19 spread: Evidence from Lebanon. Math. Biosci. Eng. 2020, 17, 5618–5632. [Google Scholar] [CrossRef]
  24. Younis, N.K.; Rahm, M.; Bitar, F.; Arabi, M. COVID-19 in the MENA Region: Facts and findings. J. Infect. Dev. Ctries. 2021, 15, 342–349. [Google Scholar] [CrossRef]
  25. Lebanese Ministry of Public Health. Epidemiological Surveillance Program of COVID-19 in Lebanon. Available online: https://moph.gov.lb/en/Media/view/43750/monitoring-of-covid-19 (accessed on 7 May 2021).
  26. Abou-Hamdan, M.; Hamze, K.; Abdel, A.; Akl, H.; El-Zein, N.; Dandache, I.; Abdel-Sater, F. Variant analysis of the first Lebanese SARS-CoV-2 isolates. Genomics 2020, 113, 892–895. [Google Scholar] [CrossRef]
  27. Feghali, R.; Merhi, G.; Kwasiborski, A.; Hourdel, V.; Ghosn, N.; Tokajian, S. Genomic characterization and phylogenetic analysis of the first SARS-CoV-2 variants introduced in Lebanon. PeerJ 2021, 9, 1–15. [Google Scholar] [CrossRef]
  28. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  29. GISAID. hCoV-19 Latest Updates by GISAID; GISID: Munich, Germany, 2021. [Google Scholar]
  30. Astuti, I. Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2): An overview of viral structure and host response. Diabetes Metab. Syndr. Clin. Res. Rev. 2020, 14, 407–412. [Google Scholar] [CrossRef] [PubMed]
  31. Arya, R.; Kumari, S.; Pandey, B.; Mistry, H.; Bihani, S.C.; Das, A.; Prashar, V.; Gupta, G.D.; Panicker, L.; Kumar, M. Structural insights into SARS-CoV-2 proteins. J. Mol. Biol. 2021, 433. [Google Scholar] [CrossRef]
  32. European Centre for Disease Prevention and Control. Rapid Increase of a SARS-CoV-2 Variant with Multiple Spike Protein Mutations Observed in the United Kingdom. Available online: http://covid19-country-overviews.ecdc.europa.eu/#34_United_Kingdom (accessed on 6 May 2021).
  33. Zhou, H.; Chen, Y.; Zhang, S.; Niu, P.; Qin, K.; Jia, W.; Huang, B.; Zhang, S.; Lan, J.; Zhang, L.; et al. Structural definition of a neutralization epitope on the N-terminal domain of MERS-CoV spike glycoprotein. Nat. Commun. 2019, 10, 1–13. [Google Scholar] [CrossRef]
  34. Nelde, A.; Bilich, T.; Heitmann, J.S.; Maringer, Y.; Salih, H.R.; Roerden, M.; Lübke, M.; Bauer, J.; Rieth, J.; Wacker, M.; et al. SARS-CoV-2-derived peptides define heterologous and COVID-19-induced T cell recognition. Nat. Immunol. 2021, 22, 74–85. [Google Scholar] [CrossRef]
  35. The World Bank. Population, Total-Lebanon. Available online: https://data.worldbank.org/indicator/SP.POP.TOTL?locations=LB (accessed on 7 May 2021).
  36. Center for Disease Control and Prevention. Genomic Surveillance for SARS-CoV-2 Variants. Available online: https://www.cdc.gov/coronavirus/2019-ncov/cases-updates/variant-surveillance.html#print (accessed on 7 May 2021).
  37. Davies, N.G.; Abbott, S.; Barnard, R.C.; Jarvis, C.I.; Kucharski, A.J.; Munday, J.D.; Pearson, C.A.B.; Russell, T.W.; Tully, D.C.; Washburne, A.D.; et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Science 2021, 372, 3055. [Google Scholar] [CrossRef]
  38. Firth, A.E.; Brierley, I. Non-canonical translation in RNA viruses. J. Gen. Virol. 2012, 93, 1385–1409. [Google Scholar] [CrossRef] [PubMed]
  39. Namy, O.; Rousset, J.-P. Specification of Standard Amino Acids by Stop Codons. In Recording: Expansion of Decoding Rules Enriches Gene Expression; Springer: New York, NY, USA, 2010; pp. 79–100. ISBN 9780387893822. [Google Scholar]
  40. Lebanese SARS-CoV-2 sequences on GISAID. Available online: https://www.epicov.org/epi3/frontend#30f5b0 (accessed on 5 May 2021).
  41. Lebanese SARS-CoV-2 sequences on NCBI’s GenBank. Available online: https://www.ncbi.nlm.nih.gov/nuccore/?term=lebanon%20covid-19 (accessed on 5 May 2021).
Biology 10 00531 g001 550
Figure 1. Percentage pie chart distribution of (A) GISAID clades and (B) Pango lineages among the 58 analyzed Lebanese SARS-CoV-2 strains. Each (A) clade or (B) lineage, as well as the percentages, are shown either on or just outside the corresponding section.
Figure 1. Percentage pie chart distribution of (A) GISAID clades and (B) Pango lineages among the 58 analyzed Lebanese SARS-CoV-2 strains. Each (A) clade or (B) lineage, as well as the percentages, are shown either on or just outside the corresponding section.
Biology 10 00531 g001
Biology 10 00531 g002 550
Figure 2. Time course GISAID clade distribution of the analyzed Lebanese SARS-CoV-2 strains over the last year. A data table with legend keys is shown directly under the chart.
Figure 2. Time course GISAID clade distribution of the analyzed Lebanese SARS-CoV-2 strains over the last year. A data table with legend keys is shown directly under the chart.
Biology 10 00531 g002
Biology 10 00531 g003 550
Figure 3. Phylogenetic tree representing the relationship among the 58 Lebanese SARS-CoV-2 and 41 reference genomes: (a) Using MegaX, genome sequences were first aligned with MUSCLE, and the relationship tree was established via maximum likelihood, with a bootstrap value of 500, followed by a rooting on Wuhan-Hu-1; (b) SARS-CoV-2 strains grouping according to the GISAID clades; (c) the various Pango Lineages within a certain cluster: the number of each lineage over the total number of strains within a cluster is shown. Color key is shown in the bottom right corner.
Figure 3. Phylogenetic tree representing the relationship among the 58 Lebanese SARS-CoV-2 and 41 reference genomes: (a) Using MegaX, genome sequences were first aligned with MUSCLE, and the relationship tree was established via maximum likelihood, with a bootstrap value of 500, followed by a rooting on Wuhan-Hu-1; (b) SARS-CoV-2 strains grouping according to the GISAID clades; (c) the various Pango Lineages within a certain cluster: the number of each lineage over the total number of strains within a cluster is shown. Color key is shown in the bottom right corner.
Biology 10 00531 g003
Biology 10 00531 g004 550
Figure 4. Various amino acid substitutions detected in: (a) the spike and (b) the nucleocapsid proteins. The number of strains in which each substitution was detected is shown next to said substitution on the wheel.
Figure 4. Various amino acid substitutions detected in: (a) the spike and (b) the nucleocapsid proteins. The number of strains in which each substitution was detected is shown next to said substitution on the wheel.
Biology 10 00531 g004
Table
Table 1. Analyzed Lebanese SARS-CoV-2 strains.
Table 1. Analyzed Lebanese SARS-CoV-2 strains.
Virus NameGISAID Accession Number (Genbank if Available)Collection DateGISAID CladePango Lineage
FAS/2021EPI_ISL_11593754-January-2021GRYB.1.1.7
LAU1-53249EPI_ISL_63711013-August-2020GB.1.398
LAU2-53435EPI_ISL_637111August-2020GB.1.398
LAU3-53451EPI_ISL_63711214-August-2020GB.1.398
LAU4-53460EPI_ISL_63711314-August-2020GB.1.398
LAU-R110 1EPI_ISL_10090851-September-2020GHB.1
LAU-R113 1EPI_ISL_10091231-September-2020GHB.1.36
LAU-R115 1EPI_ISL_10091271-September-2020GB.1.398
LAU-R117 1EPI_ISL_10091281-September-2020GB.1.398
LAU-R129 1EPI_ISL_10091591-September-2020GB.1.398
LAU-R136 1EPI_ISL_1009160 1-October-2020GB.1.398
LAU-R145 1EPI_ISL_10092131-October-2020GB.1.398
LAU-R148 1EPI_ISL_10092141-October-2020GB.1.398
LAU-R151 1EPI_ISL_10096171-September-2020GB.1.398
LAU-R153 1EPI_ISL_1009626 1-September-2020GB.1.398
LAU-R170 1EPI_ISL_10096281-September-2020GHB.1.438
LAU-R186 1EPI_ISL_10096291-August-2020GRB.1.1.315
LAU-R204 1EPI_ISL_1009630 1-July 2020GB.1.398
LAU-R231 1EPI_ISL_10096311-August-2020GB.1.398
LAU-R242 1EPI_ISL_1009656 1-December-2020GB.1.398
LAU-R249 1EPI_ISL_10096581-December-2020GB.1.398
LAU-R43 1EPI_ISL_9822981-June-2020GB.1.398
LAU-R45 1EPI_ISL_1009011 1-June-2020GB.1.398
LAU-R59 1EPI_ISL_10090121-June-2020GB.1.398
LAU-R70 1EPI_ISL_10090131-July 2020GB.1.398
LAU-R73 1EPI_ISL_1009014 1-August-2020GB.1.398
LAU-R77 1EPI_ISL_10090151-August-2020GB.1.398
LAU-R87 1EPI_ISL_10090161-September-2020GHB.1.468
LAU-R92 1EPI_ISL_1009017 1-September-2020GB.1.398
LAU-uk1EPI_ISL_76875121-December-2020GB.1.398
LAU-uk2EPI_ISL_76879821-December-2020GB.1.398
LAU-uk3EPI_ISL_76879921-December-2020GB.1.398
LAU-uk4EPI_ISL_76874322-December-2020GRB.1.1
LB-R8 1EPI_ISL_498551 (MT801000)19-March-2020GB.1
LB-R9 1EPI_ISL_498552 (MT801001)18-March-2020GRB.1.1.106
LB-R10 1EPI_ISL_498554 (MT801002)20-March-2020GHB.1
LB-R12 1EPI_ISL_498556 (MT801003)20-March-2020GHB.1
LEB-UK-BEPI_ISL_107298530-Decmember-2020GRYB.1.1.7
S1_758EPI_ISL_4505084-March-2020GHB.1
S2_759EPI_ISL_45050915-March-2020GHB.1
S3_760EPI_ISL_4505104-March-2020OB.4
S4_761EPI_ISL_4505114-March-2020GRB.1.1.29
S5_762EPI_ISL_45051227-February-2020OB.4
S6_766EPI_ISL_4505139-March-2020GRB.1.1.29
S7_763EPI_ISL_45442021-February-2020GHB.1
S8_767EPI_ISL_45051413-March-2020GHB.1
S9_764EPI_ISL_45051511-March-2020OB.4
S10_768EPI_ISL_45051615-March-2020GHB.1
S11_765EPI_ISL_4505179-March-2020GHB.1
LBN_1-B.1.1.7_2021(MW686007)4-January-2021GRYB.1.1.7
LBN_2-B.1.1.7_2021(MW692113)4-January-2021GRYB.1.1.7
LBN_3-B.1.1.7_2021(MW692114)8-January-2021GRYB.1.1.7
LBN_4-B.1.1.7_2021(MW692115)8-January-2021GRYB.1.1.7
LBN_5-B.1.1.7_2020(MW692116)28-December-2020GRYB.1.1.7
LBN_6-B.1.1.7_2020(MW692117)28-December-2020GRYB.1.1.7
LBN_7-B.1.1.7_2020(MW692118)24-December-2020GRYB.1.1.7
LBN_8-B.1.1.7_2020(MW692119)24-December-2020GRYB.1.1.7
LBN_9-B.1.258Delta(MW720771)12-December-2020GB.1.258
1 Strains sequenced in this study.
Table
Table 2. Randomly selected strains from each GISAID clade.
Table 2. Randomly selected strains from each GISAID clade.
Virus Name 1GISAID Accession Number
(Genbank if Available)		Collection DateGISAID CladePango Lineage
England MLK-9E05B3EPI_ISL_60144320-September-2020GRYB.1.1.7
Ukraine Kyiv-PHC-0030EPI_ISL_149513818-March-2021GRYB.1.1.7
Poland 0252WSN2021_wsseolEPI_ISL_149334925-March-2021GRYB.1.1.7
Belgium ULG-13710EPI_ISL_149310419-March-2021GRYB.1.1.7
Sweden 19482581EPI_ISL_149259315-February-2021GRYB.1.1.7
England 20168037604EPI_ISL_46661516-February-2020GRB.1.1.1
Brazil MG-LBI172EPI_ISL_149497622-February-2021GRP.2
Poland 0221WSN2021_wsseolEPI_ISL_149331822-March-2021GRB.1.1
Czech Republic NRL-3602-21EPI_ISL_149219110-February-2021GRB.1.1
Russia SPE-RII-MH14837SEPI_ISL_14917329-February-2021GRB.1.1.397
Spain CT-ISCIII-2013597EPI_ISL_53954826-June-2020GVB.1.177
Spain PV-HUD-67010593EPI_ISL_14959781-February-2021GVB.1.177
Austria CeMM6068EPI_ISL_14958376-November-2020GVB.1.177
Portugal PT5231EPI_ISL_14948949-February-2021GVB.1.177
Wales PHWC-PYYSCAEPI_ISL_147648514-March-2021GVB.1.177
Canada ON-PHL-8751EPI_ISL_4183451-February-2020GHB.1.*
USA OR-TRACE-BENT-031521-431EPI_ISL_149493814-March-2021GHB.1.429
Portugal PT5252EPI_ISL_14949149-February-2021GHB.1.367
Colombia VAC-INS-VG-1482EPI_ISL_14949525-March-2021GHB.1.596
USA OR-TRACE-BENT-030821-421EPI_ISL_14949308-March-2021GHB.1.429
Wuhan WIV04 2019EPI_ISL_40212430-December-2019LB
England TFCI-2707EC6EPI_ISL_14768283-April-2020LB.3
Japan PG-24179EPI_ISL_142976119-February-2020LB.12
Kazakhstan 123631EPI_ISL_13411501-March-2020LB
Hong Kong HKPU-00440EPI_ISL_128943926-March-2020LB.3
Wuhan-Hu-1(NC_045512.2)1-December-2019LB
Italy LAZ-INMI-SPL1EPI_ISL_41297429-January-2020VB.2
South Korea KDCA2676EPI_ISL_149016718-February-2020VB.41
Canada QC-CHUM-2008201265AEPI_ISL_137897322-March-2020VB.40
Thailand CU-SI2011343-NTEPI_ISL_12964502-April-2020VB.40
Slovenia 21-2850EPI_ISL_128987430-March-2020VB
Guangdong 20SF012EPI_ISL_40393214-January-2021SA
France ARA-700000210773EPI_ISL_149023611-March-2021SA.27
Scotland QEUH-1409CFDEPI_ISL_138994718-March-2021SA
Switzerland BS-UHB-351261EPI_ISL_138812425-January-2021SA.27
Luxembourg LNS7034377EPI_ISL_13838551-February-2021SA.28
Germany BY-ChVir-929EPI_ISL_40686228-January-2020GB.1
Brazil MG-LBI57EPI_ISL_149496527-January-2021GP.2
USA IL-IDPH-WHI-C-001647EPI_ISL_149450314-July 2020GB.1.234
Australia NSW4430EPI_ISL_14947213-April-2021GB.1
Mexico SLP-UASLP-AH1COV2SS027_S16EPI_ISL_149472522-May-2020GB.1
1 The country of origin of each strain is indicated prior to its name.
Table
Table 3. List of uncommon AA substitutions detected in the structural proteins S, N, and E of the analyzed Lebanese SARS-CoV-2 strains.
Table 3. List of uncommon AA substitutions detected in the structural proteins S, N, and E of the analyzed Lebanese SARS-CoV-2 strains.
SARS-CoV-2 StrainAA SubstitutionPrevalence in Lebanese StrainsPrevalence between the Strains Reported on GISAID 1First/Only Report According to GISAID’s CoVsurverNumber of CountriesStrain Collection Date
Spike protein
LAU-R186spike F59Y1106First401-Aug-20
spike G35S16First2
spike N17S1797First8
spike R44S12First2
spike Y38stop11Only1
LAU-R45spike P39T18First601-Jun-20
spike S31Y11Only1
spike V11D11Only1
LAU-R204spike L10R11Only101-July 20
spike L5P12First2
LAU-R87spike D40V11Only101-Sep-20
LAU-R129spike L611R15First201-September-20
LAU-R151spike P82T13First201-September-20
LAU-R43spike Y91stop11Only101-Jun-20
spike F92L12Last2
Nucleocapsid protein
LAU-R145N E367stop13First301-Oct-20
LAU-R231N L352stop11Only101-Aug-20
LAU-R43N L407M11Only101-Jun-20
LAU-R153N M411K11Only101-September-20
LAU-R92N S318W11Only101-September-20
Envelope protein
LAU-R110E C44S16First301-September-20
S6E V62D11Only109-Mar-20
1 Up to 18 April 2021.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fayad, N.; Abi Habib, W.; Kandeil, A.; El-Shesheny, R.; Kamel, M.N.; Mourad, Y.; Mokhbat, J.; Kayali, G.; Goldstein, J.; Abdallah, J. SARS-CoV-2 Variants in Lebanon: Evolution and Current Situation. Biology 2021, 10, 531. https://doi.org/10.3390/biology10060531

AMA Style

Fayad N, Abi Habib W, Kandeil A, El-Shesheny R, Kamel MN, Mourad Y, Mokhbat J, Kayali G, Goldstein J, Abdallah J. SARS-CoV-2 Variants in Lebanon: Evolution and Current Situation. Biology. 2021; 10(6):531. https://doi.org/10.3390/biology10060531

Chicago/Turabian Style

Fayad, Nancy, Walid Abi Habib, Ahmed Kandeil, Rabeh El-Shesheny, Mina Nabil Kamel, Youmna Mourad, Jacques Mokhbat, Ghazi Kayali, Jimi Goldstein, and Jad Abdallah. 2021. "SARS-CoV-2 Variants in Lebanon: Evolution and Current Situation" Biology 10, no. 6: 531. https://doi.org/10.3390/biology10060531

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