- Research article
- Open Access
- Open Peer Review
Polymorphisms of Plasmodium falciparum k13-propeller gene among migrant workers returning to Henan Province, China from Africa
BMC Infectious Diseasesvolume 17, Article number: 560 (2017)
Henan Province has been in the malaria elimination stage, with all reports of the disease being imported since 2012 and over 90% coming from Africa. Surveillance and population studies are essential for the early detection and subsequent prevention of the spread of drug resistance. The K13-propeller gene was recently identified as a proposed molecular marker of artemisinin (ART) resistance. In this study, we detected mutations of the K13-propeller gene in samples taken from imported malaria cases in Henan Province from 2012 to 2015.
There were 483 samples that were obtained from Plasmodium falciparum–infected malaria migrant workers who returned to Henan Province from Africa between 2012 and 2015. The single nucleotide polymorphisms in the K13-propeller gene were assessed by nested PCR with DNA sequencing. Frequency and geographic difference of K13-propeller gene mutant types were analyzed.
Of 483 patients, 476 were cured and 7 died. There were no K13-propeller mutations in the blood samples from the 7 patients who died, but there were 23 different genotypes of the K13-propeller that were observed in 24 (4.97%) of the samples. C580Y, which was the predominant one in the resistance of ART, was not detected in the samples, but R539T and P574L which have also been associated with ART resistance, were observed in two samples from Angola and Equatorial Guinea. No mutations were detected in 11 samples from North Africa. The frequency of the K13-propeller was 6.50% (8/123) in Central Africa, followed by East Africa (1/19, 5.26%), West Africa (9/198, 4.55%) and South Africa (6/132, 4.55%). There was no significant difference among these four areas (P = 0.795).
R539T and P574L were found in migrant workers who traveled from Africa to Henan Province, although the frequency of the K13-propeller mutants was low. These data may enrich the molecular surveillance of antimalarial resistance and will be helpful for developing and updating the antimalarial policy in Henan Province.
Malaria is the most critical human parasitic disease caused by Plasmodium, it is a serious threat to human health, and there were about 214 million cases of malaria and 438,000 deaths from it in 2015, of which 90% occurred in Africa . In China, an action plan for malaria elimination was launched in 2010, and the Henan government carried out the action the same year. Indigenous P. falciparum was successfully eliminated in 1988 and no indigenous malaria cases have been reported in Henan Province since 2012 . However, imported malaria cases have increased year after year. P. falciparum which was imported from Africa, has been mainly responsible for the increasing impacts [3, 4].
The effectiveness of antimalarial medicines is a powerful guarantee for malaria control and elimination. Artemisinin-based combination therapies (ACTs) were recommended as the first-line treatment against uncomplicated P. falciparum infection in 2006 by the World Health Organization (WHO) . Artemisinin has proven to be the most successful antimalarial drug over 10 years, and the morbidity and mortality of malaria have fallen dramatically due to the use of ACT in highly endemic areas of the world . However, as with other drugs, the curative effect of ACT has declined gradually along with its use. The emergence of P. falciparum resistance to artemisinin and its derivatives was first reported in 2008 in Western Cambodia . In 2009, Dondorp et al.  also reported that the susceptibility of P. falciparum in vivo to ACT was reduced. Then, ACT resistance became prevalent in the Greater Mekong Subregion (GMS) [9,10,11]. K13-propeller, as a molecular marker of artemisinin-resistance, was identified by Ariey in 2014 . Since then, numerous studies have been conducted about K13-propeller mutations. The K13-propeller mutations associated with artemisinin-resistance were mainly found in Southeast Asia, with C580Y being the predominant one. At the same time, many studies were carried out in Africa. While few resistance-associated mutations have been found in Africa , there have been numerous mutations found in the K13-propeller. In 2015, Feng et al.  reported that C580Y and R539T were observed, and the most prevalent mutation was C580Y, with the rate of 2.7% in the samples from Ghana.
As a proposed molecular marker of ART resistance, the K13-propeller was recently identified and has been studied in an increasing number of reports. Mutations in the K13-propeller were first reported in 2008 in western Cambodia . More mutations in this gene have been also subsequently found in resistant parasites in samples from Vietnam and Burma. To date, more than 200 nonsynonymous mutations in the K13-propeller gene have been reported, of which mutations of N458Y, Y493H, R539T, I543T, R561H, and C580Y were validated to be associated with ART resistance, and M476I and M579I were also reported to be associated with resistance in vivo or in vitro tests [15, 16]. In this study, we evaluated the prevalence of polymorphisms in the K13-propeller gene from imported P. falciparum patients who returned from Africa to Henan Province during 2012- 2015.
Sample collection and DNA extraction
Blood samples were obtained from migrant workers with P. falciparum infection prior to treatment. All of the patients returned from Africa to Henan Province between 2012 and 2015. Thick and thin blood films were smeared and stained with 10% Giemsa for 20 min. The final diagnosis was confirmed by microscopic examination of Giemsa-stained thick blood films and nested PCR. For each patient, approximately 200 μl of finger-prick blood was spotted on Whatman No. 903 filter paper. All of the blood samples were labeled with study numbers, names, and dates and stored in individual plastic bags at −20 °C until use. DNA was extracted from the blood samples using the QIAamp DNA Mini kit (QIAGEN Inc., Germany) according to the manufacturer’s instructions and stored at −20 °C for use in PCR assays.
K13-propeller gene amplification and sequencing
The K13-propeller gene was amplified by the nested PCR method as previously described . Nested PCR was performed to amplify the fragments of the K13-propeller using a DNA thermal cycler (Mastercycler nexus, Eppendorf Ltd., Germany). The first-round PCR was performed in a 25-μl reaction volume containing 8.5 μl of ddH2O, 1.0 μl each of the primers Artinner F (GCCTTGTTGAAAGAAGCAGAA) and Artouter R (CGCCATTTTCTCCTCCTGTA) (10 μmol/l), 12.5 μl of 2 × Go Taq Green Master Mix (Promega, USA), and 3.0 μl of DNA template. The PCR reaction condition was 95 °C for 15 min, followed by 35 cycles at 94 °C for 1 min, 58 °C for 1 min, and 72 °C for 2 min, and a final extension at 72 °C for 5 min. The second round contained the primers ArtinnerF and ArtinnerR (GTGGCAGCTCCAAAATTCAT) and 5 μL of the first-round product as the template, with the same cycling conditions as the first round. The amplified products were sequenced by Sangon Biotech Co. Ltd. (Shanghai, China). The primers were also synthesized by Sangon Biotech Co. Ltd.
Sequencing alignments and data analysis
Sequences were analyzed with the BLAST program. Sequences were aligned to reference PF3D7_1343700 (http://www.plasmodb.org) using the BioEdit Sequence Alignment Editor. The K13-propeller allele frequency was calculated with Microsoft Excel to assess the differences by SPSS 17.0. Pearson’s chi-square test was used to determine the significance of the results. P values were calculated and considered to be statistically significant at <0.05.
Epidemiological and clinical characteristics of patients
A total of 483 samples were collected from migrant workers who had returned from 27 countries in Africa and were diagnosed with P. falciparum in Henan Province during 2012-2015. Of them, 89 were collected in 2012, 125 were obtained in 2013, 141 were taken in 2014, and 128 were collected in 2015. Of the 483 patients, the majority (40.99%, 198/483) had returned from West Africa, followed by South Africa (27.33%, 132/483) and Central Africa (25.47%, 123/483). Only a minority of the patients came back from East Africa (3.93%, 19/483) and North Africa (2.28%, 11/483) (Table 1). The male: female ratio was 79.5:1 (477/6), and the age range was 17 to 70 with the mean age of 37.76 ± 9.30 years.
The clinical characteristics, including fever, chills, sweating, headache, and diarrhea, were analyzed. There were 99.79% of the patients who had a fever and 21.95% (106/483) who had serious complications, including cerebral lesion, coma, and hemolysis (see Table 2).
After diagnosis, the patients with no complications were treated with Dihydroartemisinin and Piperaquine phosphate tablets, while the others were treated with injected Artesunate at the same time for symptomatic treatment. The results of the treatment were that most patients recovered, with the exception of the 7 patients who died; K13-propeller mutations were not observed in the samples from these 7 patients.
K13-propeller point mutations
A 751-bp fragment was amplified by nested PCR for each of the 483 samples. All of the products were sequenced successfully, and 4.97% (24/483) contained single nucleotide polymorphisms (SNPs) at 23 locations, of which 8 were described previously and 15 were unreported (Table 3). Notably, the R539T, and P574L mutations that were associated with ART resistance were observed in the samples from Angola and Equatorial Guinea. The M476I substitution was also present in one sample from Equatorial Guinea. The M476I mutation which was produced with artemisinin tolerance in vitro in a Tanzanian strain, was developed specifically to test ACT activity in a test tube. However, none of the I543T, Y493H, and C580Y mutations, which were reported to be directly associated with ART resistance, were observed in these samples.
Of the 23 mutations, there were 4 synonymous mutations and 19 nonsynonymous mutations. The mutations occurred in 1-3 clones (0.21-0.62%). Of the nonsynonymous mutations, N664D was the most prevalent (0.62%, n = 3), while R471R was the most prevalent (0.41%, n = 2) of the synonymous mutations. The rest were found in only one sample. Furthermore, six patients were found to have six mixed genotypes: I476I589, K556Y648, T578I579, T626A627, N648D664, and C662D664 (Table 3).
Distribution of K13-propeller point mutations
A total of 23 point mutations were detected in the isolates from 27 countries in Africa during 2012-2015 in Henan Province. The mutations were distributed in 11 countries. There were no mutations detected in 11 isolates from 2 countries in North Africa. Five polymorphisms (R471R, R539T, M579I, Q613E, and N664D) were detected in 132 isolates from 4 countries in South Africa, with 4.55% (6/132) prevalence. There was only one mutation, A578S, which was detected in 19 isolates from four countries in East Africa, with 5.26% (1/19) prevalence. There were 12 polymorphisms (D464N, G496G, E556K, A578S, K610R, A626T, A627A, N629S, D648N, D648Y, F662C, and N664D) that were found in 198 isolates from 10 countries in West Africa, with 4.55% (9/198) prevalence. There were 10 polymorphisms (D464N, M476I, P574L, A578T, M579I, V589I, Q613E, K658Q, L663 V, and N664 N) that were detected in 123 isolates from 7 countries in Central Africa, with 6.50% (8/123) prevalence (see Fig. 1, Tables 1 and 3). There was no significant difference among these 4 areas (P = 0.795).
In this study, the sequencing of the K13-propeller gene from 483 P. falciparum samples collected from Africa identified 23 point mutations in 24 samples. We observed 4 synonymous and 19 nonsynonymous mutations, of which 8 mutations had been reported before and 15 had yet to be described [18,19,20,21]. All of the samples were from 27 countries in Africa, and 24 samples with mutations were distributed in 11 countries. The K13-propeller mutations were most prevalent in Equatorial Guinea, Angola, and Nigeria; P574L and R539T were detected in the samples from Equatorial Guinea and Angola. In our samples, we could not detect C580Y in the 23 isolates from Ghana, a possible reason for this might be that the samples from Ghana were small. The results in our study were not consistent with the report by Feng et al.  and other reports with no mutations associated with artemisinin resistance in Africa [17, 22], but the low mutation rate agreed with other reports in samples from Africa [23, 24]. The findings in our study indicated that the mutants of the K13-propeller gene alleles exist in migrant workers returning from Africa, which required us to focus on these individuals and strengthen the early monitoring of the K13-propeller.
In this study, the total rate of mutation was 4.97% and the major mutations were from West Africa, South Africa, and Central Africa. The difference of the mutations in these regions was not significant. In the 102 samples from Angola, one K13-propeller resistance-confirmed mutation, R539T was observed in one sample and the mutation rate was 0.98%. We also detected another resistance-confirmed mutation, P574L, in one isolate from Equatorial Guinea and the rate was 1.54% (1/65). The other mutations associated with ART resistance were not found in our study, but the extensive distribution of the K13-propeller mutations and the emergence of the mutations (P574L and R539T) being directly related to ART resistance suggests that there are potential ART-resistant genes in Africa that could make ART-resistant mutations spread all over the world. Routine monitoring must continue in order to ensure that the recommended ACTs are effective, that timely changes to national treatment policies can be implemented, and that ART resistance can be detected early.
The mutations of the K13-propeller detected in our study mainly came from Equatorial Guinea, Angola, and other West African, South African, and Central African countries. Furthermore, the mutations corresponding to the most prevalent mutations in Cambodia were not observed in East Africa. Historically, the emergence of drug-resistant P. falciparum strains (to chloroquine and sulphadoxine-pyrimethamine) first occurred in Southeast Asia, spread to East Africa, and then moved outwards to the rest of the African continent . But in our study, this did not seem to be so. The mutations might be due to increased international travel and migration.
In China, ACTs have been used against uncomplicated P. falciparum infection since 2009 [26,27,28]. So far no evidence has indicated that there was ART resistance in China. However, in recent years, more than 90% of the registered malaria cases were imported cases in China [29, 30]. To understand the distribution of antimalarial drug resistance is essential for the prevention, control, and elimination of malaria in China. Early surveillance of the occurrence and spreading of antimalarial drug resistance would be greatly enhanced by the application of valid antimalarial resistance molecular markers .
Since 2010 Henan Province has been in the malaria elimination stage, with all reports of the disease being imported and over 90% coming from migrant workers returning from Africa after 2012. Angola, Nigeria, and Equatorial Guinea were the top three countries from which the patients returned [32, 33]. In our study, 483 patients returned from 27 countries in Africa, and half of them came back from: Angola (102), Nigeria (72), and Equatorial Guinea (65). Although C580Y was not observed and the two other mutations that are associated with ART resistance were at low frequencies, this study showed that the mutated K13-propeller gene was in the migrant workers who returned from Africa, requiring our attention to the emergence of the resistance to ACTs in Africa. The use of the molecular marker K13-propeller is fundamental for surveillance in malaria control programs, in order to prolong the life span of the ACTs in Africa. Surveillance and population studies are essential for the early detection and subsequent prevention of the spread of drug resistance . It is also important to change the treatment policy in a timely manner, and the information will be useful for developing and updating antimalarial guidance. Therefore, surveillance should continuously be strengthened in migrant workers returning to Henan Province, and the data may be used to enhance molecular surveillance of antimalarial resistance and informed decisions in the rational treatment policy in Henan Province.
We observed 23 point mutations, including 4 synonymous mutations and 19 nonsynonymous mutations, in the K13-propeller in migrant workers who had returned from Africa. C580Y was not detected, but R539T and P574L were found in 483 samples from Africa. It was meaningful to study artemisinin resistance in returning migrant workers in Henan Province, and these data should help provide a reasonable treatment policy in Henan Province.
Artemisinin-based Combination Therapies;
Greater Mekong Subregion
Single Nucleotide Polymorphisms
World Health Organization
World Health Organization. World malaria report. Geneva: WHO; 2016. http://www.who.int/malaria/publications/world_malaria_report_2016/en/
Liu Y, Zhang HW, Zhou RM, Yang CY, Qian D, Zhao YL, et al. First imported relapse case of Plasmodium vivax malaria and analysis of its origin CSP sequencing in Henan Province. China Malar J. 2014;13:448. doi:10.1186/1475-2875-13-448.
Feng J, Yan H, Feng XY, Zhang L, Li M, Xia ZG et al. Imported malaria in China, 2012. Emerg Infect Dis, 2014;20:1778-1780. doi:org/10.1186/1475-2875-13-29.
Zhang HW, Su YP, Zhao XD, Yan QY, Liu Y, Chen JS. [Imported falciparum malaria situation in Henan Province during 2005-2009] (In chinese). Zhongguo Ji Sheng Chong Xue Yu Ji Sheng Chong Bing Za Zhi. 2010;28: 476-477.
World Health Organization. Guidelines for the treatment of malaria. Geneva: WHO; 2006.
Wang LQ, Liu YH, Zhao SX, Yang ZQ. [advances in the study of artemisinin resistance in Plasmodium falciparum and methods of its detection] (in chinese). Journal of Pathogen Biology. 2014;9:1142–4. doi:10.13350/j.cjpb.141221.
Noedl H, Se Y, Schaecher K, Smith BL, Socheat D, Fukuda MM. Evidence of artemisinin-resistant malaria in western Cambodia. N Engl J Med. 2008;359:2619–20. doi:10.1056/NEJMc0805011.
Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, Lwin KM, Ariey F, Hanpithakpong W, Lee SJ, et al. Artemisinin Resistance in Plasmodium falciparum Malaria. N Engl J Med. 2009;361:455–67. doi:10.1056/NEJMoa0808859.
Muller O, Sie A, Meissner P, Schirmer RH, Kouyate B. Artemisinin resistance on the Thai-Cambodian border. Lancet. 2009;374:1419. doi:10.1016/S0140-6736(09)61857-2.
Wang ZL, Shrestha S, Li XL, Miao J, Yuan LL, Cabrera M, Grube C, Yang ZQ, Cui LW. Prevalence of K13-propeller polymorphisms in Plasmodium falciparum from China-Myanmar border in 2007–2012. Malar J. 2015;14:168. doi:10.1186/s12936-015-0672-9.
Takala-Harrison S, Jacob CG, Arze C, Cummings MP, Silva JC, Dondorp AM, Fukuda MM, Hien TT, Mayxay M, Noedl H, et al. Independent emergence of Artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia. J Infect Dis. 2014; 211:670–679. doi:https://doi.org/10.1093/infdis/jiu491.
Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois AC, Khim N, Kim S, Duru V, Bouchier C, Ma L, et al. A molecular marker of artemisininresistant Plasmodium falciparum malaria. Nature. 2014;505:50–5. doi:10.1038/nature12876.
Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, Sreng S, Anderson JM, Mao S, Sam B, et al. Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014;371:411–23. doi:10.1056/NEJMoa1314981.
Feng J, Li J, Yan H, Feng MY, Xia ZG. Evaluation of Antimalarial resistance marker polymorphism in returned migrant Workers in China. Antimicrob Agents Chemother. 2015;59:326–30. doi:10.1128/AAC.04144-14.
World Health Organization. Artemisinin and artemisinin-based combination therapy resistance. Geneva: WHO; http://apps.who.int/iris/bitstream/10665/250294/1/WHO-HTM-GMP-2016.11-eng. pdf.
Zhang MH, Lu F, Cao J, Gao Q. [progress in researches on molecular markers of Plasmodium falciparum drug resistance] (in Chinese). Chin J Schisto Control. 2015;27:323–7. doi:10.16250/j.32.1374.2014227.
Taylor SM, Parobek CM, DeConti DK, Kayentao K, Coulibaly SO, Greenwood BM, Tagbor H, Williams J, Bojang K, Njie F, et al. Absence of putative Artemisinin resistance mutations among Plasmodium falciparum in sub-Saharan Africa: a molecular epidemiologic study. J Infect Dis. 2015;211:680–8. doi:10.1093/infdis/jiu467.
Nyunt MH, Hlaing T, Oo HW, Tin-Oo LK, Phway HP, Wang B, Zaw NN, Han SS, Tun T, San KK, et al. Molecular assessment of artemisinin-resistance markers, polymorphisms in the K13 propeller and a multidrug-resistance gene, in eastern and western border areas of Myanmar. Clin Infect Dis. 2015;60:1208–15. doi:10.1093/cid/ciu1160.
Tun KM, Imwong M, Lwin KM, Win AA, Hlaing TM, Hlaing T, Lin K, Kyaw MP, Plewes K, Faiz MA, et al. Spread of artemisinin-resistant Plasmodium falciparum in Myanmar: a cross-sectional survey of the K13 molecular marker. Lancet Infect Dis. 2015;15:416-421. doi:org/10.1016/S1473-3099(15)70032-0.
Aye A. Win AA, Imwong M, Kyaw MP, Woodrow CJ, Chotivanich K, Hanboonkunupakarn B, Pukrittayakamee S. K13 mutations and pfmdr1 copy number variation in Plasmodium falciparum malaria in Myanmar. Malar J 2016; 15:110. doi: org/10.1093/cid/ciu1160.
Ménard D, Khim N, Beghain J, Adegnika AA, Shafiul-Alam M, Amodu O, et al. A worldwide map of Plasmodium falciparum K13-propeller polymorphisms. N Engl J Med. 2016;37:2453–64. doi:10.1056/NEJMoa1513137.
Kamau E, Campino S, Amenga-Etego L, Drury E, Ishengoma D, Johnson K, Mumba D, Kekre M, Yavo W, Mead D, et al. K13-propeller polymorphisms in Plasmodium falciparum parasites from sub-Saharan Africa. J Infect Dis. 2015;211:1352–5. doi:10.1093/infdis/jiu608.
Torrentino-Madamet M, Fall B, Benoit N, Camara C, Amalvict R, Fall M, Dionne P, Fall KB, Diatta B, Diemé Y, et al. Limited polymorphisms in k13 gene in Plasmodium falciparum isolates from Dakar, Senegal in 2012–2013. Malar J. 2014;13:472. doi:10.1186/1475-2875-13-472.
Sylla K, Abiola A, Kouly Tine RC, Faye B, Sow D, Ndiaye JL, Ndiaye M, LO AC, Folly K, Ndiaye LA, et al. Monitoring the efficacy and safety of three artemisinin based-combinations therapies in Senegal: results from two years surveillance. Malar J. 2013;13:598. doi:10.1186/1471-2334-13-598.
Wongsrichanalai C, Pickard AL, Wernsdorfer WH, Meshnick SR. Epidemiology of drug-resistant malaria. Lancet Infect Dis. 2002; 2:209–218. doi: http://dx.doi.org/10.1016/S1473-3099(02)00239-6.
WHO. The application of artemisinin monotherapies for malaria treatment was demanded to end by WHO (In Chinese). ADRJ.2006;8:152.
Antimalarial drug policy in China. http://www.nhfpc.gov.cn/zwgkzt/wsbysj/200907/41610.shtml.
Chinese Center for Disease Control and Prevention. China’s technical scheme of malaria elimination. Beijing: Chinese Center for Disease Control and Prevention; 2011.
Feng J, Xiao HH, Xia ZG, Zhang L, Xiao N. Analysis of malaria epidemiological characteristics in the People's Republic of China, 2004-2013. Am J Trop Med Hyg. 2015;93:293–9. doi:10.4269/ajtmh.14-0733.
Zhang L, Zhou SS, Feng J, Fang W, Xia ZG. [malaria situation in the people’ s republic of China in 2014] (in Chinese). Chin J Parasitol Parasit Dis. 2015;33:321–6.
Feng J, Zhou DL, Lin YX, Xiao HH, Yan H, Xia ZG. Amplification of pfmdr1, pfcrt, pvmdr1, and K13 propeller polymorphisms associated with Plasmodium falciparum and Plasmodium vivax isolates from the China-Myanmar border. Antimicrob Agents Chemother. 2015;59:2554–9. doi:10.1128/AAC.04843-14.
Zhou RM, Zhang HW, Yang CY, Liu Y, Zhao YL, Li SH, Qian D, Xu BL. Molecular mutation profile of pfcrt in plasmodium falciparum isolates imported from Africa in Henan province. Malar J. 2016;15:265. doi:10.1186/s12936-016-1306-6.
Yang CY, Qian D, Chen WQ, Liu Y, Zhou RM, Zhang HW, et al. Chin J Schisto Control. 2016;28:444–6. doi:10.16250/j.32.1374.2015200.
Mvumbi DM, Kayembe J, Situakibanza H, BobangaTL, Nsibu CN, Mvumbi GL, Melin P, Mol PD, Hayette M. Falciparum malaria molecular drug resistance in the Democratic Republic of Congo: a systematic review. Malar J. 2015; 14:354. doi: 10.1186/s12936-015-0892-z.
The authors thank ZQX and the workers from City level CDC, who participated in the study.
This work was supported by the Medical Science and Technique Project of Henan, China (No. 201304053), and the Special Funding of the Henan Health Science and Technology Innovation Talent Project (No.4045). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The datasets used and /or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was reviewed and approved by the Project of Medical Science and Technique of Henan province (No. 201304053) and the Special Funding of the Henan Health Science and Technology Innovation Talent Project (No.4045). Before sampling, the purpose was illustrated to participants and consent was obtained from those who agreed to participate in the survey.
Consent for publication
Written informed consent was obtained from the patients for the publication of this manuscript.
The authors declare that they have no interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.