cov-lineages / pango-designation

Repository for suggesting new lineages that should be added to the current scheme
Other
1.04k stars 97 forks source link

New potential recombinant BA.1.1/BA.2.9 mainly circulanting in Panama and other countries (53 genomes) #759

Closed Vjimenez-vasquez closed 2 years ago

Vjimenez-vasquez commented 2 years ago

Description New potential recombinant BA.1.1 (lacking G8393A) / BA.2.9 mainly circulanting in Panama identified by the group of Genomic Surveillance of INS-Peru. This potential recombinant presents a breakpoint in the middle of ORF1a. This genomes cannot be assigned by PANGO and are misidentified (unassigned flag) as belong to XQ by NEXTCLADE, however their mutational profile are no compatible with this recombinant lineage. We have recently identified 9 additional genomes from Peru.

Posible recombination squeme recombination_breaks pink box indicates BA.1.1 genome contribution and green box indicates BA.2.9 contribution.

Mutations mutations

Genomes genomes.txt

sc2rf simplified output

sc2rf_output

sc2rf ouput indicates one breakpoint and identifies BA.1 and BA.2 as parental genomes.

Earliest sequence: 2022-03-19 hCoV-19/Panama/M221120-GMI/2022 (EPI_ISL_11859557)

Most recent sequence: 2022-05-30 hCoV-19/USA/NC-CORVASEQ-1085-4923/2022 (EPI_ISL_13163618)

Countries circulating Chile(1), Colombia(2), Denmark(2), Mexico(1), Panama(41) and USA (7)

Vjimenez-vasquez commented 2 years ago

Just an observation. I have noted that parental BA.2.9 strain for this recombinant has been first detected in Costa Rica (Central America) by the end of january 2022 and contains the additional mutation A24433G. Then dispersed to Canada and USA and finally to Panama where the recombination event probably took place.

Mutations

BA 2 9_mutations

Geographic contribution

Important representation in both Costa Rica and Panama.

geographic_distrution

International comparison

Costa Rica and Pamana have currently deposited 12k and 4k total genomes, respectively. It means the atribution of this BA.2.9 strain to COVID-19 cases in both countries, should bot be neglilible.

international_comparison

Relative growth advantage

CovSPECTRUM estimates 27% of relative growth rate

relative_growth_advantage

Wen1953 commented 2 years ago

Good observation @Vjimenez-vasquez.It seems to be a new recombinant. @InfrPopGen @corneliusroemer @chrisruis could it be?

agolsby commented 2 years ago

Excellent work in tracking down this apparent recombinant!

I'm concerned about the growth advantage, as CoV-Spectrum seems to have generated a curve ill-fitted to the underlying data. With full caution that single-chart epidemiology is highly prone to overinterpretation, I'll carefully mention the following early impressions from these data:

My concern is that its failure to continue amplifying indicates inability to compete with the likes of BA.2.12.1 and BA.5. However, there are plenty of other possible explanations. By facilitating tracking of this lineage, a designation could help answer questions about its epidemiologic significance in the current environment.

c19850727 commented 2 years ago

@Vjimenez-vasquez it popped out of nowhere after nearly a month, now in the UK, and with a bunch of additional amino changes (and reversions!). image

c19850727 commented 2 years ago

BTW, this potential recombinant has 141 genomes as of today.

agolsby commented 2 years ago

254 as of today when using the CoV-Spectrum query below. 286 on its accompanying UShER tree, where this clade is labelled "proposed759".

It seems concentrated mainly in Panama and the United States. Exportations include Israel, Denmark, Colombia, Peru, and Mexico.

Cladistic evidence in the UShER tree matching this query reveals multiple weeks of transmission in the USA with three multi-state clusters: one bearing C4246T, one with G8371A, and one with G13948T (which is also present in Panama). (I can't entirely rule out that sequencing artifacts created false subclades.)

https://cov-spectrum.org/explore/World/AllSamples/Past6M/variants?nucMutations=C9344T%2CA9424G%2CC9534T%2CC9866T%2CC11514T%2CC15714T%2CC17410T%2CC25624T%2CA24433G%2CC241T%2CC2470T%2CA2832G%2CC3037T%2CT5386G

Unlike the above-discussed parent clade, this recombinant does not show signs of having lost its growth advantage. The following two images show today's estimates. First is from all samples, second is from surveillance samples only.

0627 proposed759 growth 0627 proposed759 growth surveillance only

The new sequences strengthen the case that an epidemiologic event has occurred.

c19850727 commented 2 years ago

5 sequences from the UK now, all with additional AA changes mentioned above. (I am tempted to propose a sublineage of this current proposal if the growth continues)

Vjimenez-vasquez commented 2 years ago

Hi there, just an update of this proposal. There are currently 316 identified genomes by USHER.

proposition_759

InfrPopGen commented 2 years ago

Thanks for submitting. We've added recombinant lineage XAM with 341 newly designated sequences, and 3 updated designations.