Skip to main content

Identification of severe fever with thrombocytopenia syndrome virus genotypes in patients and ticks in Liaoning Province, China

Abstract

Background

Severe fever with thrombocytopenia syndrome (SFTS), caused by the SFTS virus (SFTSV), is an acute infectious disease transmitted by ticks that has recently been identified. There are no reports of epidemic serotypes in Liaoning Province, PR China. The aim of this study was, therefore, to identify genotypes of SFTSV in this province.

Methods

In 2019, quantitative PCR testing was performed on 17 patients suspected of being infected with SFTS in Liaoning Province and on 492 ticks from the counties and cities surrounding the patients’ residences. Four samples were subjected to virus isolation and whole-genome amplification.

Results

Molecular diagnostic results confirmed SFTSV infection in five of the 17 suspected cases of SFTS and in 12 of the 492 ticks, with a prevalence of 2.4%. Four strains of SFTSV were successfully isolated from patients’ blood and ticks. Phylogenetic analysis after whole-genome amplification and sequencing showed that they all belonged to genotype A of SFTSV.

Conclusions

This study is the first to determine the genotype of SFTSV in patients and ticks in Liaoning Province, PR China. The results deepen our understanding of the SFTS epidemic and provide information on the variability in mortality rate among genotypes.

Graphical Abstract

Background

Severe fever with thrombocytopenia syndrome (SFTS) is a new acute infectious disease that was first reported in 2009 in rural areas of the Hubei and Henan provinces in PR China [1]. The pathogen of the disease was initially classified as genus Phlebovirus [1], but in 2020 it was designated as Dabie bandavirus (genus Bandavirus; family Phenuiviridae) by the International Committee on Taxonomy of Viruses (ICTV). The main clinical manifestations of SFTS are fever, anorexia, muscle pain, chills, lymphadenopathy, leukopenia, thrombocytopenia and multiple organ failure [2, 3]. Severe cases can even lead to death.

SFTS virus (SFTSV) is an arbovirus mainly carried by ticks [4]. Haemaphysalis longicornis, Amblyomma testudinarium, Ixodes nipponensis and Rhipicephalus microplus are the four main tick species that cause infections [5]. It can also be transmitted from person to person by direct contact with the blood or mucus of an infected person [6,7,8]. The general population is susceptible; therefore, its epidemic area tends to expand following initial infection [5]. Thousands of incidents of SFTS have been reported in more than 16 provinces and cities in China [9] as well as in the Republic of Korea [10], Japan [11] and the USA [12].

Liaoning Province is located in northeast China and borders North Korea, Russia, Japan and South Korea. Its eastern territory is mountainous and inhabited by ticks. SFTS is prevalent in Liaoning Province, with 183 cases confirmed in 2011–2014 [9]. Its prevalence [13], symptoms [2, 14] and trends [15] in Liaoning Province have been studied, but so far there are no reports of pathogen isolation and dominant epidemic strains.

The aim of this study was to identify genotypes of SFTSV in Liaoning Province. We therefore attempted to isolate the SFTS pathogen from the blood of individuals with suspected SFTS and from vector ticks. Our findings provide insights into the prevalence of SFTSV in Liaoning Province, as well as its differentiation.

Methods

Specimen collection

Between May and October 2019, 17 suspected cases of SFTS were reported by the Shenyang Infectious Disease Hospital, Shenyang, Liaoning Province. All 17 patients exhibited symptoms and were farmers, mainly from forest areas, who had been in contact with animals (rats) or been bitten by tick(s). Serum was collected from these patients and stored at − 80 °C until testing. A total of 492 adult tick samples (not collected in Beipiao City) were collected from patients’ residences and surrounding areas, including Dandong and Dalian cities, which are in forest areas of the Changbai Mountain range extending along the border area between China and North Korea (Fig. 1). After morphological identification, the tick samples were ground in liquid nitrogen. PBS was then added followed by centrifugation at 10,000 g for 10 min as preparation for nucleic acid extraction of the suspension.

Fig. 1
figure 1

Region of patients and ticks in the study. Number coding: 1‒6 locations where samples were taken from 17 patients, 1‒7 locations where ticks were collected, 4, 7 locations were SFTSV was isolated from ticks. SFTSV, Severe fever with thrombocytopenia syndrome

Nucleic acid extraction

Total nucleic acids were extracted from the sera of the 17 patients suspected of having SFTS and from the tick suspensions, using the EasyPure Viral DNA/RNA Kit (TransGen Biotech, Beijing, China) according to the manufacturer’s instructions.

Tick species identification

Tick species were first identified based on morphological characters, and then the 16S ribosomal (rRNA) of partial samples was amplified and sequenced following a previously reported method [16]. For a detailed analysis of Haemaphysalis, we analyzed the extended sequence by PCR-restriction fragment length polymorphism (RFLP) of the 16S rRNA gene [17].

SFTSV detection

SFTSV was detected using the following primer pairs (560 bp) targeting the M gene [18]: MF3 (5′-GATGAGATGGTCCATGCTGATTCT-3′) and MR2 (5′-CTCATGGGGTGGAATGTCCTCAC-3′). Reverse transcription and PCR amplification were performed using the HiScript II Q RT SuperMix for quantitative (q)PCR (+ gDNA wiper; Vazyme Biotech, Nanjing, China), according to the manufacturer's instructions, after DNase treatment of total nucleic acids. PCR conditions were: pre-denaturation for 5 min at 94 °C, followed by 40 cycles of 94 °C for 30 s, 55 °C for 30 s and 72 °C for 1 min.

Virus isolation

Tick suspensions and patient sera that were determined positive for SFTSV by RT-PCR [19] were filtered through a 0.22 µm bacterial filter and mixed with the cell maintenance solution (Dulbecco's modified Eagle medium containing 2% foetal bovine serum). The mixture was added to the Vero cell layer and cultured at 37 °C under 5% CO2. Cytopathic effects (CPEs) were observed daily. Cells exhibiting CPEs were collected and verified by RT-PCR to verify successful infection. The infected cells were stored at − 80 °C until analysis.

SFTSV whole-genome amplification

The SFTSV genome is composed of three independent gene segments: L segment (6368 bp; encodes the RNA-dependent RNA polymerase), M segment (3378 bp; encodes the membrane protein precursor), and the antisense S segment (1744 bp; encodes the non-structural S protein and nucleocapsid protein) [20]. Primers that were previously designed for the amplification of S (three primer sets), M (seven primer sets), and L (13 primer sets) segments were used for whole-genome sequencing [21] (Table 1). Amplification was performed on a thermal cycler at 94 °C for 30 s, 55 °C (63 °C) for 30 s, and 72 °C for 1 min for 35 cycles.

Table 1 SFTSV genome sequencing primer list

SFTSV whole-genome sequencing and analysis

The PCR products were verified by agarose gel electrophoresis and sent to Sangon Biotech (Shanghai) Co., Ltd. for sequencing. Sequences were assembled using the SeqMan program of DNASTAR 7.0.

Full-length sequences of 49 SFTSV isolates from China, Japan and the Republic of Korea, collected from 2010 to 2016, were downloaded from GenBank; these isolates cover the six known (A‒F) genotypes of SFTSV. MEGA6.0 software was used to conduct a comprehensive analysis of the SFTSV strain, with the aim to identify the genotype of the virus and analyze the genetic evolution of the strain. The maximum likelihood (ML) method was used for multiple alignments to build a phylogenetic tree based on the L, M and S sequences. A bootstrap method with 1000 replications was used to estimate the reliability of the results.

Results

Tick species identification

The 16S rRNA sequences of 59 samples taken for verification purposes had 89.33–100% identity with those of H. longicornis (Additional file 1: Dataset S1), and were found to be 87.56–100% identical. Based on PCR–RFLP analysis of the 16S rRNA gene, 389 tick samples were identified as H. longicornis, and three tick samples [P(7), P(65), P(a96)] were identified as Haemaphysalis juxtakochi.

SFTSV nucleic acid detection

Of the 17 patients with suspected SFTS, five (2 men, 3 women) were confirmed to be infected with SFTSV through molecular diagnostic tests. Of the 492 tick samples, 12 samples were positive for SFTSV nucleic acids, with a positivity rate of 2.4%.

Virus isolation results

Four SFTSV strains were successfully isolated from two positive patients and two tick samples [P(36) and P(47)]; these were designated LNHDD2019-9 (Dandong City, Liaoning Province), LNHDG2019-6 (Donggang City, Liaoning Province), LNTEDG2019-6 (Erdaogou, Dandong City, Liaoning Province) and LNTYJ2019-6 (Yanji City, Jilin Province) based on the geographic location of the sample source (Fig. 1).

Systematic evolutionary analysis of the SFTSV strain

The L, M and S fragments of the four SFTSV strains were amplified separately and assembled. The sequences of S, M and L fragments of SFTSV (LNHDD2019-9/Human/China/2019 in the evolutionary tree) were uploaded to the NCBI database with accession numbers MT232960, MT232961 and and MT232962, respectively.

The nucleotide (and deduced amino acid) sequence similarities of the isolates to other genetic A-type strains were 95.95‒100% (98.90‒100%), 97.12‒100% (98.32‒100%), and 95.08‒100% (99.59‒100%) in the L, M and S segments, respectively (Fig. 2), suggesting that the four SFTSV strains isolated are genotype A strains.

Fig. 2
figure 2

Comparison (% similarity) of nucleotide (top right) and amino acid (bottom left) sequences of L, M and S segments of the SFTSV genotype. The four SFTSV strains successfully isolated in this study are shown in the red boxes

The amino acid sequences of the M and S segments have been used to distinguish the six SFTSV genotypes [22]. Therefore, we used the ML method for phylogenetic analysis based on amino acid sequences. The results showed that the L, M and S segments of the four SFTSV isolates identified in this study clustered in a clade of genotype A (Fig. 3a–c) and were closest to the viruses isolated in Jiangsu, Zhejiang and other provinces.

Fig. 3
figure 3

Phylogenetic analyses. Phylogenetic analysis based on the complete nucleotide sequences of the L, M and S segments of the four SFTSV isolates identified (filled red circles). All four isolates are SFTSV genotype A strains.

The L and S segment sequences of SFTSV isolated from the ticks in Dandong City, Liaoning Province, and Yanji City, Jilin Province, were 100% identical to those of SFTSV strains isolated from patients. In contrast, the similarity of the M segment sequences was 99% (Table 2). These results indicated that the main transmission sources of SFTS in patients were ticks distributed throughout the local or nearby areas.

Table 2 Sequence information for SFTSV isolates

Discussion

In this study, the detection, isolation and evolution analysis of SFTSV isolates were carried out using serum collected from patients suspected of having SFTS and from ticks in the area surrounding their homes. The results show that five of 14 patients with suspectged SFTS and 12 of 492 ticks were positive for SFTSV. Whole-genome sequencing of the four isolated strains of SFTSV revealed that all were SFTSV genotype A. This is the first report of genotyping on the SFTSV genotype in Liaoning Province, and the findings have deepened our understanding of the SFTS epidemic.

To the best of our knowledge, this is the first work identifying SFTSV genotypes in patients with SFTS and ticks in Liaoning Province. Previous studies in Liaoning have focused on epidemiological characteristics and symptoms and mixed symptoms of SFTS. The SFTSV strains isolated in this study were all genotype A, consistent with previous reports that the dominant strain of SFTSV in China is genotype A.

Northeast China is characterized by a mountainous terrain with abundant natural resources, geographical complexity and significant biodiversity, all factors that provide an ecological and biological basis for the survival and reproduction of ticks. The tick species identified in this study are H. longicornis, Dermacentor silvarum, D. nuttalli and Ixodes persulcatus [19]. In a previous study, H. longicornis was the dominant species collected in northeast China and the only tick species collected in Liaoning province [19]. Our results showing that most of the ticks collected in the present study were H. longicornis verifies this earlier finding. However, three samples were confirmed to be H. juxtakochi by PCR–RFLP analysis of the 16S rRNA gene. This tick species has not been found in China earlier. Further investigation and verification are required in future studies.

In this study, the SFTSV-positive rate of ticks was 2.4%, which is lower than that (3.0%) reported in Jilin province [19]. However, the number of SFTS cases (708) in Liaoning Province in 2010–2019 was much higher than that (< 94) in Jilin Province [23], which may be related to the activities (working and social) of the people in the area, their awareness of potential illness associated with ticks and testing methods.

Reports from the East Asian region show that the average mortality rate of SFTSV varies greatly, from 5.3 to 16.2% in China [9] to 20% in Japan [24] and 23.3% in South Korea [25]; in comparison, mortality was 5.3% in 2011–2014 [9] and 3.2% in 2010–2019 [23] in Liaoning Province. Among the six known A-F genotypes of SFTSV [26], genotype B is associted with the highest morbidity and mortality rates. In one study, the incidence of the F genotype was lower than that of the B genotype, but the associated mortality rate of the former was higher, while genotype A had the lowest mortality rate [26]. The patients with SFTS in this study were infected from the end of March to the end of June. They all developed the most frequent symptoms of SFTS, such as fever, headache, dizziness, muscle joint pain, nausea and vomiting. After treatment, they all recovered. The identification of SFTSV genotype A in the present study explains the low mortality rate of patients with SFTS in Liaoning Province.

The most common genotypes of SFTSV in China are genotypes A, D and F, which are associated with relatively low mortality rates. In contrast, the most common genotype of SFTSV in South Korea is genotype B, with a high mortality rate [22, 26, 27]. Because of the geographical proximity between Liaoning and South Korea, genotype B may soon spread to China. Therefore, it is important to improve methods for detecting SFTSV genotypes to control and diagnose highly lethal SFTS.

In Liaoning Province, H. longicornis is the dominant tick species; the other species identified include Haemaphysalis concinna, H. japonica and I. persulcatus. In the present study, only H. longicornis-infected ticks were identified. Further studies are required to identify different tick species carrying other SFTSV. At the same time, this study only collected ticks in counties and cities of the patients. Tick species in other counties and towns and the SFTSV infection pattern must be further investigated.

Conclusion

The study identified SFTSV genotypes in patients with suspected SFTS and in ticks in Liaoning Province. The results enhance our understanding of the SFTS epidemic in Liaoning Province. Because the highly lethal genotype B is the main epidemic strain in South Korea, it is critical to improve methods for detecting SFTSV genotypes to prevent the occurrence and spread of highly lethal SFTSV.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article.

Abbreviations

CK:

Creatine kinase

ML:

Maximum likelihood

SFTS:

Severe fever with thrombocytopenia syndrome

SFTSV:

SFTS virus

PBS:

Phosphate-buffered saline

RFLP:

Restriction fragment length polymorphism

References

  1. Yu XJ, Liang MF, Zhang SY, Liu Y, Li JD, Sun YL, et al. Fever with thrombocytopenia associated with a novel bunyavirus in China. N Engl J Med. 2011;364:1523–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Bing T, Da Q, Akihito S, Jingyi C, Baocheng D. Acute pancreatitis in patients with severe fever with thrombocytopenia syndrome virus infection. Pancreatology. 2020;20:1631–6.

    Article  Google Scholar 

  3. Sun Y, Guo B, Yan H, Wu AL, Yao WW, Chen K, et al. Patient with severe fever with thrombocytopenia syndrome virus infection and central nervous system disturbance in Dongyang, Zhejiang Province, China, 2017. Virol J. 2019;16:129.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Liu Y, Li Q, Hu W, Wu J, Wang Y, Mei L, et al. Person-to-person transmission of severe fever with thrombocytopenia syndrome virus. Vector Borne Zoonotic Dis. 2012;12:156–60.

    Article  PubMed  Google Scholar 

  5. Seo JW, Kim D, Yun N, Kim DM. Clinical update of severe fever with thrombocytopenia syndrome. Viruses. 2021;13:1213.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Chen H, Hu K, Zou J, Xiao J. A cluster of cases of human-to-human transmission caused by severe fever with thrombocytopenia syndrome bunyavirus. Int J Infect Dis. 2013;17:e206–8.

    Article  PubMed  Google Scholar 

  7. Huang D, Jiang Y, Liu X, Wang B, Shi J, Su Z, et al. A cluster of symptomatic and asymptomatic infections of severe fever with thrombocytopenia syndrome caused by person-to-person transmission. Am J Trop Med Hyg. 2017;97:396–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Hu J, Li Z, Cai J, Liu D, Zhang X, Jiang R, et al. A cluster of bunyavirus-associated severe fever with thrombocytopenia syndrome cases in a coastal plain area in China, 2015: identification of a previously unidentified endemic region for severe fever with thrombocytopenia bunyavirus. Open Forum Infect Dis. 2019;6:ofz209. https://doi.org/10.1093/ofid/ofz209.

  9. Li Y, Zhou H, Mu D, Yin W, Yu H. Epidemiological analysis on severe fever with thrombocytopenia syndrome under the national surveillance data from 2011 to 2014, China. Zhonghua Liu Xing Bing Xue Za Zhi. 2015;36:598–602.

    PubMed  Google Scholar 

  10. Denic S, Janbeih J, Nair S, Conca W, Tariq WU, Al-Salam S. Acute thrombocytopenia, leucopenia, and multiorgan dysfunction: the first case of SFTS Bunyavirus outside China. Case Rep Infect Dis. 2011. https://doi.org/10.3390/v13061086.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Iwao K, Kawaguchi T, Kimura M, Iwao C, et al. Severe fever with thrombocytopenia syndrome accompanied by invasive pulmonary aspergillosis: an autopsy case. Viruses. 2021;13:1086.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Sharma D, Kamthania M. A new emerging pandemic of severe fever with thrombocytopenia syndrome (SFTS). Virusdisease. 2021;32:220–7.

    Article  PubMed Central  Google Scholar 

  13. Wang ZJ, Zhang J, Mao LL, Liu Y, Sun TT, Tian J, et al. Epidemiological characteristics of severe fever with thrombocytopenia syndrome in Liaoning province, 2010–2013. Zhonghua Yu Fang Yi Xue Za Zhi. 2016;50:362–3.

    CAS  PubMed  Google Scholar 

  14. Lu QB, Li H, Jiang F-C, Mao L-L, Liu X-S, Wang N, et al. The differential characteristics between severe fever with thrombocytopenia syndrome and hemorrhagic fever with renal syndrome in the endemic regions. Open Forum Infect Dis. 2019;6:0fz477.. https://doi.org/10.1093/ofid/ofz477.

  15. Sun J, Lu L, Wu H, Yang J, Liu K, Liu Q. Spatiotemporal patterns of severe fever with thrombocytopenia syndrome in China, 2011–2016. Ticks Tick Borne Dis. 2018;9:927–33.

    Article  PubMed  Google Scholar 

  16. Black WC 4th, Piesman J. Phylogeny of hard- and soft-tick taxa (Acari: Ixodida) based on mitochondrial 16S rDNA sequences. Proc Natl Acad Sci USA. 1994;91:10034–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Thompson AT, Dominguez K, Cleveland CA, Dergousoff SJ, Falco RC, Greay T, et al. Molecular characterization of haemaphysalis species and a molecular genetic key for the identification of Haemaphysalis of North America. Front Vet Sci. 2020;7:141.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Li J-X, Yun SM, Kim SY, Wen P, et al. Ticks and their potential infection with severe fever with thrombocytopenia syndromes virus in Yanbian, Jilin. Acta Parasitologica ET Medica Entomologica Sinica. 2017;24:132–40.

  19. Liu H, Li Z, Wang Z, He B, Wang S, Wei F, et al. The first molecular evidence of severe fever with thrombocytopenia syndrome virus in ticks in Jilin Northeastern China. Ticks Tick Borne Dis. 2016;7:1280–3.

    Article  PubMed  Google Scholar 

  20. Sun Y, Qi Y, Liu C, Chen P, Fu L, Peng B, et al. Nonmuscle myosin heavy chain IIA is a critical factor contributing to the efficiency of early infection of severe fever with thrombocytopenia syndrome virus. J Virol. 2014;88:237–48.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Shan J, Tang Z, Cui LB, Xian QI, et al. Detection of RNA and analysis of complete genome of severe fever with thrombocytopenia syndrome virus (SFTSV) in Jiangsu. Modern Prev Med. 2013;40:2686–96.

    CAS  Google Scholar 

  22. Won YJ, Kang LH, Lee SG, Park SW, Han JI, Paik SY. Molecular genomic characterization of severe fever with thrombocytopenia syndrome virus isolates from South Korea. J Microbiol. 2019;57:927–37.

    Article  CAS  PubMed  Google Scholar 

  23. Huang X, Li J, Li A, Wang S, Li D. Epidemiological characteristics of severe fever with thrombocytopenia syndrome from 2010 to 2019 in Mainland China. Int J Environ Res Public Health. 2021;18(6):3092. https://doi.org/10.3390/ijerph18063092.

  24. Yasukawa M. Human herpesvirus 6 in hematological diseases. Rinsho Ketsueki. 2008;49:247–56.

    PubMed  Google Scholar 

  25. Wu XB, Na RH, Wei SS, Zhu JS, Peng HJ. Distribution of tick-borne diseases in China. Parasit Vectors. 2013;6:119.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Yun SM, Park SJ, Kim YI, Park SW, Yu MA, Kwon HI, et al. Genetic and pathogenic diversity of severe fever with thrombocytopenia syndrome virus (SFTSV) in South Korea. JCI Insight. 2020;5:e129531.

    Article  PubMed Central  Google Scholar 

  27. Beom PK, Yul RJ, Yoon KS, Kim NY, Jo YH, Jeong RD, et al. Deep sequencing and phylogenetic analysis of severe fever with thrombocytopenia syndrome virus from the tick, Haemaphysalis longicornis, in Korea. Entomol Res. 2020;51:3–11.

    Google Scholar 

Download references

Acknowledgements

We thank The Sixth People's Hospital of Shenyang for providing patient samples.

Funding

This work was supported by the Scientific and Technological Assistance Projects for Developing Countries (KY201901014). Title: China-Mongolia "One Belt, One Road" Joint Laboratory of Zoonosis. This work was also supported by the National Important Scientific & Technology Project (2018ZX10101002-002). Title: Major vector-borne pathogens in the Belt and Road.

Author information

Authors and Affiliations

Authors

Contributions

ZC and XH conceived and designed the study. YM, HL, DL, YW, FJ, QG, FJ and XH performed the experiments and analyzed the data; FJ participated in sample collection; YM, GS, and ZC drafted and revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Bao-Shan Liu, Guo-Shun Shen or Ze-Liang Chen.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Shenyang Sixth People's Hospital (Letter Number: 2020-01-001-01). Written informed consent was obtained from 17 patients with suspected SFTS (tick-borne diseases).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: 16S.mas.

The alignment of the 16S rRNA sequences of 59 tick samples.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Han, XH., Ma, Y., Liu, HY. et al. Identification of severe fever with thrombocytopenia syndrome virus genotypes in patients and ticks in Liaoning Province, China. Parasites Vectors 15, 120 (2022). https://doi.org/10.1186/s13071-022-05237-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-022-05237-3

Keywords