A SARS-CoV-2 lineage A variant (A.23.1) with altered spike has emerged and is dominating the current Uganda epidemic

Status of the SARS-CoV-2 epidemic in Uganda

SARS-CoV-2 infection was first detected in Uganda in March 2020, initially among international travellers until passenger flights were stopped in late March 2020. A second route of virus entry with truck drivers from adjacent countries then became apparent (8). Since August 2020, community transmission dominated the Uganda case numbers. By March 2021 total cases in Uganda were 40,535, with 334 deaths attributed to the virus. We have continued our efforts to generate SARS-CoV-2 genomic sequence data to monitor virus movement and genetic changes and we report here on a novel sub-lineage A (A.23.1) that emerged and is dominating the local epidemic. The A.23.1 variant encodes multiple changes in the spike protein as well as in nsp6, ORF8 and ORF9, some predicted to be functionally similar to those observed in VOC in lineage B.

Changes in prevalence of lineage A viruses

The genomes generated here were classified into Pango lineages(4) using the Pangolin module pangoLEARN (https://github.com/cov-lineages/pangolin) and into NextStrain clades using NextClade (9) (https://clades.nextstrain.org/). The distribution of virus lineages circulating in Uganda changed dramatically over the course of the year. A clear feature of the earlier COVID-19 epidemic in the country was the diversity of viruses found throughout the country attributed to frequent flights into Uganda from Europe, UK, US and Asia; this is reflected in the 9 lineages seen from March to May 2020 with a mixture of both lineage A and B viruses (Figure 1, panel a). After passenger flights were limited in March 2020, the virus entered via land travel with truck drivers. Uganda is landlocked country, characterised by its important geographical position, i.e. the crossing of two main routes of the Trans-Africa Highway in East Africa. The essential nature of produce and goods transport allowed virus movement from/to Kenya, South Sudan, DRC, Rwanda and Tanzania. In the period of June to August 2020 lineage B.1 and B.1.393 strains were abundant, similar to patterns observed in Kenya (10) (Figure 1, panel b) although lineage A viruses did not decline as seen in US and Europe. Lineage A.23 strains were first observed in two prison outbreaks in Amuru and Kitgum, Uganda in August 2020 and by the September-November period, the A.23 was the major lineage circulating throughout the country (Figure 1, panel c). The A.23 virus continued to evolve into the lineage A.23.1, first observed in late October 2020. Given the diversity of virus lineages found in the country from March until November 2020, it was unexpected that by late December 2020 to January 2021, lineage A.23.1 viruses represented 90% (102 of 113 genomes) of all viruses observed in Uganda (Figure 1, panel d). In all time periods, the SARS-CoV-2 positive sample were obtained from multiple clinical and surveillance locations throughout Uganda indicating that the differences are unlikely to be due to sampling different subpopulations in the country at different times.

• Download figure
• Open in new tab

Figure 1 SARS-CoV-2 lineage diversity in Uganda.

All high-coverage complete from Uganda (n=322) were lineage typed using the pangolin resource (https://github.com/cov-lineages/pangolin) Lineage counts were stratified into four periods (panel a: March-May 2020, panel b:June-August 2020, panel c: September to November 2020 and panel d: December 2020 to January 2021). The percentage of each lineage within each set was plotted as a treemap using squarify (https://github.com/laserson/squarify) with the size of each sector proportional to the number of genomes, genomes numbers are listed with “n=”.

Virus sequence diversity including fatal cases

All newly and previously generated Uganda genomes that were complete and high-coverage (n=322) were used to construct a maximum-likelihood phylogenetic tree (Figure 2).

• Download figure
• Open in new tab

Figure 2. Maximum-likelihood phylogenetic tree comparing all available complete and high-coverage Uganda sequences (N=322).

Strain names are coloured according to the case profile (cases from the community: dark red, prison: orange, truck driver: light brown, return traveller: light blue). The case clusters from prisons in Kitgum and Amuru are highlighted in colour boxes in light yellow and light green, respectively. The lineages A.23 and A.23.1 are indicated. The tree was rooted where lineages A and B were split. The branch length is drawn to the scale of number of nucleotide substitutions per site, indicated in lower left, and only bootstrap values of major nodes were shown.

A number of A and B variant lineages were observed briefly at low frequencies and may have undergone extinction, similar to patterns observed in the UK (11) and Scotland (12). Genomes identified from a truck driver are often observed basal to community clusters (Figure 2), suggesting the importance of this route in the introduction and spread of the virus into Uganda. Most of genomes from truck drivers sampled at points of entry (POEs) bordering Kenya belonged to lineage B.1 and B.1.393 consistent with the pattern reported in Kenya (10). However, genomes identified from truck drivers from Tanzania, and from the Elegu POE bordering South Sudan, albeit small numbers, belonged to both A and B.1 lineages. Continued monitoring of truck drivers coming in and out of the Uganda provides a useful description of the inland circulation of strains in this part of world, where genomic surveillance is not as detailed as in other parts of the world.

Emergence of A.23 and A.23.1

Outbreaks of SARS-CoV-2 infections were reported in the Amuru and Kitgum prisons in August 2020 (13)(14). The SARS-CoV-2 genome sequences from individuals in the prisons were exclusively belonging to lineage A (Figure 2) with three amino acid (aa) changes encoded in the spike protein (F157L, V367F and Q613H, Figure 3) that now define lineage

• Download figure
• Open in new tab

Figure 3. Spike protein changes in lineage A.23 and A.23.1 relative to the SARS-CoV-2 reference strain (NC_045512) encoded protein are documented.

Panel a: The locations of important spike protein features are indicated. NTD: N-terminal domain, RBD: receptor-binding domain, S1: spike 1, S1: Spike 2, TM: transmembrane domain, HR1: helical repeat 1, HR2: helical repeat 2, NTD super: N-terminal domain supersite. Panel b: Each line represents the encoded spike protein sequence from a single genome, ordered by date of samples collection (bottom earliest, top most recent). Sequences from Amuru in August 2020, Kitgum in September 2020 and Uganda in October, November, December 2020/January 2021 are indicated. Markers indicating the positions of amino acid (aa) differences from the reference strain, changes observed in multiple genomes are annotated with the annotation (original aa position new aa). Panel c. Current global distribution of and A.23.1 All available SARS-CoV-2 complete genomes annotated as complete and lineage A from GISAID were retrieved on Feb 4 2021 and lineage typed using Pangolin(19). and confirmed as A.23 and A.23.1 by extracting examining the encoded spike protein. All novel Uganda A.23 and A.23.1 reported here were also included. Genomes were plotted by country and sample collection date.

A.23 (see below). By October 2020, lineage A.23 viruses were also found outside of the prisons in a community sample from Lira (a town 140 km from Amuru), in two samples from the Kitgum hospital, in several community samples from Kampala, Jinja, Mulago, Tororo, Soroti as well as in 2 truck drivers collected at POE bordering Kenya. By November 2020, the A.23 viruses spread further to northern Uganda in Gulu and Adjumani. Lineage A.23 viruses were not seen in Uganda (or anywhere in the world) before August 2020 (Figure 3 panel c), yet the A.23 viruses were attributed to 32% of the viruses in Uganda (Figure 1) from June to August 2020 and 50% of the observed viruses in September to November 2020. In late October, the A.23.1, a variant evolving from A.23, with additional change in the spike protein (P681R) was observed (Figure 3, panel b, c) and by December 2020-January 2021, 90% of identified genomes (102 out of 113) belonged to the new A.23.1 lineage (Figure 1 and 2). The mutations in A.23.1 were consistent with evolution from an original A.23 virus observed in Amuru/Kitgum cluster (Figure 2 and Supplemental Figure 1) as well as changes nsp6 and ORF9 (Supplemental Figure 2 and 4).

Important changes observed in the spike protein

The spike protein is crucial for virus entry into host cells, for tropism, and is a critical component of COVID-19 vaccine development and monitoring. The changes in spike protein observed in Uganda and global A.23 and A.23.1 viruses are shown in Figure 3 panel b.

Many amino acid (aa) changes were single events with no apparent transmission observed. However, the initial lineage A.23 genomes from Amuru and Kitgum encoded three amino acid changes in the exposed S1 domain of spike (F157L, V367F and Q613H, Figure 3 panel b). The V367F change is reported to modestly increase infectivity(15), the Q613H change may have similar consequences as the D614G change observed in the B.1 lineage found predominantly in Europe and USA; in particular, D614G was reported to increase infectivity, spike trimer stability and furin cleavage (15),(16),(17),(18). These changes were not observed in previously reported genomes from Uganda (8). Of some concern, the mutations E484K and N501Y amino acid changes in the receptor binding domain (RBD) were observed in the A.23 viruses identified in Adjumani cases on 9^th to 11^th November 2020 (Figure 3, panel b). These two amino acid changes are shown to substantially compromise the vaccine efficacy as well as antibody treatments.

Of concern, the recent Kampala and global A.23.1 virus sequences from December 2020-January 2021 now encoded 4 or 5 amino acid changes in the spike protein (now defining lineage A.23.1, see below) plus additional protein changes in nsp3, nsp6, ORF8 and ORF9 (Figure 3 panel b, Figure 4). The P681R spike change adds a basic amino acid adjacent to the spike furin cleavage site. This same change has been shown in vitro to enhance the fusion activity of the SARS-CoV-2 spike protein, likely due to increased cleavage by the cellular furin protease (20); importantly, a similar change (P681H) is encoded by the recently emerging VOC B.1.1.7 that is now spreading globally across 75 countries as of 5 February 2021 (5) (21). There are also changes in the spike N-terminal domain (NTD), a known target of immune selection, observed in samples from Kampala A.23.1 lineage, including P26S and R102I (Figure 3 panel b). Additionally and importantly, a A.23.1 strain identified in Kampala on 11^th December 2020 carried the E484K change in the RBD, which may add further concern of this particular variant as it gains higher transmissibility and enhanced resistance to vaccine and therapeutics. Outside of the spike protein, a single nucleotide change (G27870T) leading to early termination of the ORF7b (E39*) was observed in the A.23.1 from the community cases in Tororo in late December 2020. Although the clinical implication of this change is yet to be determined, it is important to document such change for further follow-up.

• Download figure
• Open in new tab

Figure 4. All protein changes across lineage new variants.

All forward open reading frames from the 35 early lineage B SARS-CoV-2 genomes were translated, and processed into 44 aa peptides (with 22 aa overlap), clustered at 0.65 identity using uclust (25), aligned with MAAFT (26) and converted into pHMMs using HMMER-3 (27). The presence of each domains and its bit-score (a measure of the similarity between the query sequence and the sequences used for the pHMM(27)) was sought in each set of SARS-CoV-2 VOC genomes and the 1-mean of the normalized domain bit-scores was plotted across the genome (e.g. 1 - the similarity of the identified query domain to the reference lineage B SARS-CoV-2 domain). Domains were coloured by the proteins from which they were derived with the colour code indicated below the figure. Panel a. Query set are 49 most Uganda lineage genomes. Panel b. All B.1.1.7 full genomes lacking Ns deposited in GISAID on Jan 26 2021, Panel c. All B.1.351 full genomes lacking Ns present in GISAID on Jan 26 2021, Panel d. All P.1 full genomes lacking Ns present in GISAID on Jan 26 2021.

New lineage A designations

The viruses detected in Amuru and Kitgum met the criteria for a new SARS-CoV-2 lineage (22)(23) by clustering together on a global phylogenetic tree, sharing epidemiological history and source from a single geographical origin, and encoding multiple defining SNPs. These features including especially the three spike changes F157L, Q613H and V367F define the new lineage A.23. Continued circulation and evolution of A.23 in Uganda was observed and two additional changes in spike R102I and P681R were observed in December 2020 in Kampala; these SNPs define the sub-lineage A.23.1. Additional changes in non-spike regions also define the A.23 and A.23.1, including nsp3: E95K, nsp6: M86I, L98F, ORF 8: L84S, E92K and ORF9 N: S202N, Q418H. These new lineages can be assigned since pangolin version v2.1.10 and pangoLEARN data release 2021-02-01.

Screening SARS-CoV-2 genomic data from GISAID (March 12, 2021), A.23 and A.23.1 viruses are now found in 26 countries outside of Uganda (Figure 3 panel c). The A.23 was first observed in Uganda in August 2020, subsequently in USA in October and Kenya and Rwanda in December (Figure 3 panel c). The A.23.1 was first seen in Uganda in the community cases in Mbale on 28^th October 2020 and in Jinja on 29^th October 2020, and soon spreading across the country in early November 2020. Outside of Uganda, the A.23.1 was found in England and Cambodia from the end of November, in Rwanda from the beginning of December. Of note, the international flights out of Uganda were restarted on 1 October 2020 with flights to Europe, Asia and USA. Phylogenetic analysis supports the close evolution of A.23 to A.23.1 (Supplementary Figure 1).

Additional changes in Ugandan A.23 and A.23.1 genomes compared to other VOC genomes

Although a main focus has been on spike protein changes, there are changes in other genomic regions of the SARS-CoV-2 virus accompanying the adaptation to human infection. We employed profile Hidden Markov Models (pHMMs) prepared from 44 amino acid peptides across the SARS-CoV-2 proteome (24) to detect and visualize protein changes from the early lineage B reference strain NC_045512. Measuring the identity score (bit-score) of each pHMMs across a query genome provides a measure of protein changes in 44 amino acid steps across the viral genome (Figure 4 panel a). This method applied to A.23 and A.23.1 genome sequences revealed the changes in spike (discussed above) as well as changes in the transmembrane protein nsp6 and the interferon modulators ORF8 and 9 (Figure 5 panel a).

We asked if a similar pattern of evolution was appearing in VOC as SARS-CoV-2 adapted to human infection. We gathered the sets of genomes described in the initial published descriptions of these VOC (B.1.1.7 (5), B.1.351 (28) or P.1 (7)) and applied the same profileHMM analysis. Similar to A.23/A.23.1, the B.1.1.7 lineage encodes nsp6, spike, ORF8 and 9 changes as well as changes in nsp3 and the RNA-dependent RNA polymerase (RRP, Figure 4 panel b). Lineage B.1.351 encodes nsp3, nsp6, RDRP, spike and ORF6 changes (Figure 4 panel c) and lineage P.1 encodes nsp3, nsp6, RDRP, nsp13, spike and ORF8 and 9 changes (Figure 4 panel d). Although the exact amino acid and positions of change within the proteins differ in each lineage, there are some striking similarities in the common proteins that have been altered. Of interest, the nsp6 change present in B.1.1.7, B.1.351 and P.1 is a 3 amino acid deletion (106, 107 and 108) in a protein loop of nsp6 predicted to be on exterior of the autophagy vesicles on which the protein accumulates (29).The three amino acid nsp6 changes of lineage A.23.1 are L98F in the same exterior loop region, and the M86l and M183I changes predicted to be in intramembrane regions but adjacent to where the protein exits the membrane (29), (Supplementary Figure 2). A compilation of the amino acid changes in A.23.1 and the VOC lineages is found in Supplementary Table 2 with proteins that are altered in all 4 lineages marked in red.

Sample collection, whole genome MinION sequencing and genome assembly

Residual nucleic acid extract from SARS-CoV-2 RT-PCR positive samples were obtained from Central Public Health Laboratory (Kampala, Uganda). The nucleic acid was converted to cDNA and amplified using SARS-CoV specific 1500bp-amplicon spanning the entire genome as previously described(33).The resulting DNA amplicons were used to prepare sequencing libraries, barcoded individually and then pooled to sequence on MinION R.9.4.1 flowcells, following the standard manufacturer’s protocol.

The genome assemblies were performed as previously described (8). Briefly, reads from fast5 files were base-called and demultiplexed using Guppy 3.6 running on the UMIC HPC. Adapters and primers sequences were removed using Porechop (https://github.com/rrwick/Porechop) and the resulting reads were mapped to the reference genome Wuhan-1 (GenBank NC_045512.2) using minimap2(34) and consensus genomes were generated in Geneious (Biomatters Ltd). Genome polishing was performed in Medaka, and SNPs and mismatches were checked and resolved by consulting raw reads.

Phylogenetic analyses

For the local Uganda virus comparison, all newly and previously generated genomes from Uganda (N=322) were aligned using MAFFT (26) and manually checked in AliView (35). The 5’ and 3’ untranslated regions (UTRs) were trimmed. Maximum-likelihood (ML) phylogenetic tree was constructed using RAxML-NG (36) under the GTR+I+G4 model as best-fitted substitution model according to Akaike Information Criterion (AIC) determined by ModelTest-NG (37) and run for 100 pseudo-replicates. Resulting tree was visualised in Figtree(38) and rooted at the point of splitting lineage And B.

For phylogenetic analyses of Uganda lineage A.23 and A.23.1 strains comparing to global A.23/A.23.1 strains, the global SARS-CoV-2 lineage A.23 (N=8) and A.23.1 (N=38) genomes were retrieved from GISAID on 12 March 2021. These global A.23/A.23.1 genomes combining with Ugandan A.23/A.23.1 genomes (N=191) were aligned using MAFFT and manually checked in AliView, followed by trimming 5’ and 3’ UTRs. The global and Ugandan A.23/A.23.1 genomes were used to construct a ML tree under the GTR+I+G4 model as best-fitted substitution model according to AIC determined by ModelTest-NG (37) and run for 100 pseudo-replicates using RAxML-NG. Resulting tree was visualised in Figtree and rooted using the A.23 lineage.

Profile Hidden Markov Model (profileHMM) domain analysis of A.23/A.23.1 and VOC genomes was performed as previously described (24) with some changes. A database of profileHMMs was generated from the first 65 lineage B SARS-CoV-2 genome sequences. All 3 forward open reading frames of each genome were translated computationally and then sliced into 44 amino acid segment with overlapping with 22 amino acids. All 44 amino acid query peptides were then clustered with uclust (25) and their original identity and coordinates determined by blastp search against a protein database made from the NC_045512 reference strain.

Query sets of genomes were processed to remove any genomes containing Ns (which disrupt the HMM scoring process). The hmmscan function from HMMER-3 (27) was used with the early B database. Query matches were identified using an E-value cutoff of 0.0001 and the bit-score values for each hit (a measure of the distance between the query 44 amino acid peptide and the B-lineage reference) was collected. Bit-scores for each domain were normalized by dividing each query score by the maximum score for that domain (x/x_max). In all analyses the original B lineage NC_045512 reference genome was included to define the maximum bit-score.