Genomic surveillance of SARS-CoV-2 in COVID-19 vaccinated healthcare workers in Lebanon

Background The emergence of SARS-CoV-2 variants including the Delta and Omicron along with waning of vaccine-induced immunity over time contributed to increased rates of breakthrough infection specifically among healthcare workers (HCWs). SARS-CoV-2 genomic surveillance is an important tool for timely detection and characterization of circulating variants as well as monitoring the emergence of new strains. Our study is the first national SARS-CoV-2 genomic surveillance among HCWs in Lebanon. Methods We collected 250 nasopharyngeal swabs from HCWs across Lebanon between December 2021 and January 2022. Data on the date of positive PCR, vaccination status, specific occupation, and hospitalization status of participants were collected. Extracted viral RNA from nasopharyngeal swabs was converted to cDNA, library prepped using the coronaHIT method, followed by whole genome sequencing on the Illumina NextSeq 500 platform. Results A total of 133 (57.1%) samples belonging to the Omicron (BA.1.1) sub-lineage were identified, as well as 44 (18.9%) samples belonging to the BA.1 sub-lineage, 28 (12%) belonging to the BA.2 sub-lineage, and only 15 (6.6%) samples belonging to the Delta variant sub-lineage B.1.617.2. These results show that Lebanon followed the global trend in terms of circulating SARS-CoV-2 variants with Delta rapidly replaced by the Omicron variant. Conclusion This study underscores the importance of continuous genomic surveillance programs in Lebanon for the timely detection and characterization of circulating variants. The latter is critical to guide public health policy making and to timely implement public health interventions. Supplementary Information The online version contains supplementary material available at 10.1186/s12920-023-01443-9.


Background
Since its emergence in December 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains a global public health threat. According to the World Health Organization (WHO), excess mortality attributed to COVID-19 had reached 14.9 million between January 2020 and December 2021, with middle-income countries accounting for 81% of the death toll [1]. In Lebanon, more than 1.1 million confirmed cases and 10,457 deaths have been reported [2]. Vaccine development against SARS-CoV-2 proceeded in an unprecedented pace with 11 vaccines granted emergency use listing (EUL) by the World Health Organization (WHO) as of July 25, 2022 [3]. More than 11.4 billion vaccine doses have been administered worldwide since the start of COVID-19 vaccine rollout in December 2020 [4]. Despite the availability of effective vaccines against SARS-CoV-2, according to a modeling study, it is generally accepted that vaccine-induced immunity wanes overtime leading to an increased incidence of breakthrough infections, especially with the emergence of new variants [5].
Since the emergence of SARS-CoV-2, new variants have evolved from the original SARS-CoV-2 strain (Wuhan 19 strain (WA1/2020)). These variants are classified into variants under monitoring (VUM), variants of interest (VOI), and variants of concern (VOC) [6,7]. VUM are variants with genetic changes that are suspected to affect virus characteristics and may pose future risk but with yet no clear evidence of phenotypic or epidemiological impact. VOI are SARS-CoV-2 variants possessing predicted or known genetics changes that affect the characteristics of the virus (transmissibility, disease severity, immune escape, therapeutic escape) and known to cause significant community transmission or multiple COVID-19 clusters in multiple countries. VOC are SARS-CoV-2 variants that meet the definition of a VOI and are associated with increased transmissibility detrimental change in COVID- 19 [6]. The WHO is continuously monitoring and assessing the evolution of SARS-CoV-2 and the emergence of new variants with increased risk to global public health.
The B.1.617.2 lineage along with its sub-lineages made up the Delta variant that was responsible for the COVID-19 surge in India, eventually spreading and dominating globally [8]. Despite the high replicative efficiency, reduced sensitivity to host immune responses, and high transmissibility of the Delta variant compared to previous VOCs, vaccine effectiveness was sustained [8]. In November 2021, the surge of cases in South Africa marked the identification of a new VOC named Omicron (B.1.1529). Omicron replaced the Delta variant and was characterized by a higher number of amino acid substitutions, higher transmissibility and partial resistance to vaccine induced immunity compared to previous VOCs [9][10][11]. Studies showed that although Omicron had higher rates of reinfection, it was clinically less severe compared to the Delta variant, possibly due to prior infections and T cell mediated immune responses [8]. Recent studies showed that vaccine effectiveness against the Omicron variant (B.1.1.529) was lower than the Delta variant (B.1.617.2) after primary immunization with 2 doses of the ChAdOx1 nCoV-19, BNT162b2 or mRNA-1273 vaccines with significant reduction in vaccine effectiveness against the two variants ≥ 25 weeks following the second dose [12]. Nevertheless, vaccine effectiveness against symptomatic disease was restored following a booster shot, underscoring the importance of a third dose [13,14].
In Lebanon, 43.6% of the population were fully vaccinated (2 doses) as of June 2022 [2]. Healthcare workers (HCWs) are at increased risk of SARS-CoV-2 infection compared to the general population and have been prioritized in COVID-19 vaccine rollout worldwide [15] as per the WHO Strategic Advisory Group of Experts framework for the allocation and prioritization of COVID-19 vaccination [16] and the Advisory Committee on Immunization Practices [17]. Consequently, HCWs were among the priority groups to receive the BNT162b2 vaccine in Lebanon. Vaccine rollouts started in mid-February 2021 and the second dose was administered 21 days following the first dose. In October 2021, the third dose (i.e. booster shot) was introduced [18]. The use of personal protective equipment (PPE) including mask mandates (for patients and visitors of healthcare facilities) significantly reduced the occupational risk of acquiring COVID-19 by HCWs [19]. However, waning immunity and the emergence of antigenically drifted VOCs meant that HCWs are at risk of breakthrough infections [20][21][22][23][24][25][26][27].
Genomic surveillance of SARS-CoV-2 was first initiated in Lebanon in 2020 with direct support from the WHO [28]. Lebanon was the first in the Eastern Mediterranean Region (EMR) to identify the Delta variant in a timely manner. However, the sampling was done randomly and mainly on travelers. In collaboration with the Ministry of Public Health (MoPH) and support from the WHO, this study allowed for the establishment of the capacity to perform an ongoing genomic surveillance of SARS-CoV-2 among HCW's to support the timely implementation of public health measures to control the spread and emergence of new SARS-CoV-2 variants and to inform global surveillance. Here, we report the first genomic data from this surveillance effort focusing primarily on breakthrough infections detected among HCWs, a group with the highest vaccination coverage among the Lebanese population.

Study design, population and data collection
This study is part of a national surveillance program in collaboration with the Epidemiological Surveillance Unit Aliquots of the collected samples were stored at − 80 °C until processed. The date of positive PCR, vaccination status, specific occupation, and hospitalization status of participants were collected.

RNA extraction and whole genome sequencing (WGS)
Aliquots of the nasopharyngeal swabs (140 µl) were used to extract total RNA following manufacturer's instructions (QIAamp Viral RNA mini Kit, QIAGEN, Hilden, Germany, Cat. 52906). Aliquots were eluted in 30 µL of Buffer AVE. Both the concentration and quality of RNA samples were measured and checked with the Denovix Blue DS-11 Spectrophotometer. Viral RNA extracts were sequenced at the Quadram Institute Bioscience, UK. Briefly, viral RNA was converted to cDNA then amplified using the ARTIC protocol v3 (LoCost) [29] and using V4 of the primer set, with sequencing libraries prepared using CoronaHiT as previously described [30]. Genome sequencing was performed using the Illumina NextSeq 500 platform (Illumina, CA, USA) with one positive control and one negative control per 94 samples. The raw reads were demultiplexed using bcl2fastq (v2.20). The reads were used to generate a consensus sequence using the ARTIC bioinformatic pipeline [31]. Briefly, the reads had adapters trimmed with TrimGalore [32] and were aligned to the WuhanHu-1 reference genome (accession MN908947.3) using BWA-MEM (v0.7.17) [33]. The ARTIC amplicons were trimmed and a consensus was built using iVAR (v.1.3.1) [34]. PANGO lineages were assigned using Pangolin (v3.1.20) [35] and PangoLEARN model dated 2022-02-02 [36]. In this manuscript, we used the Pango lineage designation system.

Demographic characteristics of study participants
Our study included a total of 250 HCWs testing positive for SARS-CoV-2 by RT-PCR between December 1, 2021 (n = 27), and January 31, 2022 (n = 223). Data on the occupation and COVID-19 vaccine status were available for only 175 HCWs working at AUBMC. Among those, the majority were females (64%) and received at least two doses of the BNT162b2 COVID-19 vaccine (95.4%). Nurses (39.4%), medical doctors (14.8%) and technicians (10.3%) accounted for the majority of samples ( Table 1). None of the SARS-CoV-2-positive HCWs were hospitalized.

SARS-CoV-2 phylogenetic analysis
Lineage analysis was performed using the Pangolin COVID-19 lineage assigner. Overall, 10 lineages were identified among HCWs (Table 2; Additional file 1:  [41]. These mutations have been shown to provide potential immune escape characteristics and higher transmission than BA.2 [41,44]. In this study, we found that Omicron BA.1.1 and BA.1 lineages were the predominant circulating lineages in our cohort of HCWs between December 2021 and January 2022. The Delta variant was detected in only 6% of our samples suggesting the replacement of the Delta variant with Omicron as the predominant circulating VOC which is consistent with global trends observed during the same period [8]. We did not detect BA.3 lineage or any recombinant lineages in our sequenced samples. The former does not have specific mutations in the spike protein but rather a combination of mutations from BA.1 and BA.2 [45]. The rate of spread of the three Omicron lineages (BA.1, BA.2 and BA.3) differs with BA.1 and BA.2 being the predominant lineages. Between December 2021 and January 2022, BA.1 lineage accounted for 78% of sequenced samples submitted to the GISAID database compared to 16% of BA.2 [46]. This is similar to our findings reflecting the predominance of BA.1 over BA.2.
The subvariant BA.2 shares 32 mutations with BA.1 but has 28 distinct mutations, four unique ones in the RBD region alone, which according to a deep learning algorithm, made it far more likely than other lineages to be the next dominant subvariant [40]. Indeed, BA.2 had already become the dominant variant in multiple countries such as Denmark and UK in February 2022 [47]. Moreover, as of June 21, 2022, 52% of sequenced samples submitted to GISAID database were BA.2 followed by 19% BA.2.12.1, 18% BA.5, 11% BA.4 and < 1% BA.1 [46]. As we continue our national genomic surveillance beyond January 2022, we expect a shift in dominance in favor of the highly transmissible BA.2 subvariant as well as an expected detection of other Omicron subvariants. Following the Delta variant wave in August 2021, the number of confirmed cases in Lebanon was relatively stable with an average of 4000 cases reported weekly. On December 22, the number of cases doubled with a weekly change rate of 69.51%, marking the beginning of a new wave of SARS-CoV-2. The numbers kept increasing gradually reaching a peak of 57,800 weekly cases on January 31. All the samples included in our study were collected during this period, and 88% of them belonged to the Omicron variant. According to the WHO COVID-19 weekly update, a similar trend was observed globally during the week of 27 December until January 2, where the total number of new cases exhibited a sharp increase of 71% as compared to previous weeks. In addition, GISAID variant tracking revealed a predominant global circulation of the Omicron variants in more than 75 countries. Data on vaccines effectiveness and breakthrough infections against SARS-CoV-2 variants are lacking in the Eastern Mediterranean region. However, a study conducted in Qatar predicted the cumulative incidence of breakthrough infections to be 0.59% and 0.84% 6 months following receiving the second dose of the mRNA-1273 vaccine and the BNT162b2 vaccine, respectively [48].
A key limitation of our study is that it did not include samples from all regions in Lebanon and thus is not fully representative of the situation in Lebanon. However, given that Beirut is the capital of this small country with a lot of population movement during the week, and particularly on the weekends when its residents travel to their villages across Lebanon, we believe that data from Beirut are to some extent representative of the whole country. We were also unable to gather clinical data of HCWs, which hampered our ability to assess risk factors associated with Omicron breakthrough infections. Moreover, data on receiving the date of the second and booster shots were unavailable and thus we were unable to estimate vaccine effectiveness between the date of receiving the booster and the date of breakthrough infection.

Conclusion
Our findings underscore the importance of continuing genomic surveillance in Lebanon in order to monitor virus evolution and the emergence of novel SARS-CoV-2 variants. This is particularly important in HCWs as they are more likely to be exposed to emerging variants and can act as an advanced warning proxy to the wider community. More recently, two Omicron lineages (BA.4 and BA.5) have been identified in South Africa before being detected in several countries worldwide including Botswana, Belgium, Denmark, the United Kingdom, France, Germany, Portugal and China [42,49,50]. Therefore, continuing genomic surveillance will help assessing the characteristics and the public health implications of these lineages and other variants that might emerge and contribute to more informed public health intervention strategies.