Abstract
Anticoagulants are a major component of rodenticides used worldwide, which function by effectively blocking the vitamin K cycle in rodents. The rat Vitamin K epoxide Reductase Complex (VKORC) subunit 1 is the enzyme responsible for recycling vitamin K, and five substitution mutations (Tyr139Cys, Tyr139Ser, Tyr139Phe and Leu128Gln and Leu120Gln) located in the VKORC1 could result in resistance to anticoagulant rodenticides. This study carried out a VKORC1-based survey to estimate the anticoagulant rodenticide resistance in three Rattus species (R. losea, R. norvegicus, and R. tanezumi) collected in Hong Kong. A total of 202 rats captured in Hong Kong between 2017 and 2021 were analysed. Sequencing of molecular marker cytochrome c oxidase subunit 1 (COX1) was carried out to assist the species identification, and the identities of 52 lesser ricefield rats (R. losea), 81 common rats (R. norvegicus) and 69 house rats (R. tanezumi) were confirmed. Three VKORC1 exons were amplified from individuals by PCR followed by Sanger sequencing. A total of 47 R. tanezumi (68.1%) contained Tyr139Cys mutation in VKORC1 gene, and half of them were homozygous. None of the collected R. losea and R. norvegicus were detected with the five known substitutions leading to anticoagulant rodenticides resistance, and previously undescribed missense mutations were revealed in each species. Whole genome sequencing was further carried out on some individuals, and single nucleotide polymorphisms (SNPs) were also identified in the introns. This is the first study investigating the situation of anticoagulant rodenticide resistance in the rats collected in Hong Kong. Given that the efficacy of rodenticides is crucial for effective rodent management, regular genetic testing as well as population genomic analyses will be required to both monitor the situation and understand the adaption of different rat haplotypes for integrated pest management. Susceptibility tests for individual rodenticides should also be conducted regularly to assess their effectiveness on local species.
Similar content being viewed by others
Introduction
Rodents have been generally regarded as pests as they cause economic losses and transmit rodent-borne diseases1,2. In Hong Kong, eight species of rats and mice had been previously identified including Bandicota indica, Mus caroli, M. musculus, Niviventer fulvescens, Rattus norvegicus, R. rattus, R. tanezumi, and R. sikkimensis3. In recent, the first ever reported transmission of rat hepatitis E virus species C genotype 1 to human had also been identified in Hong Kong4. This emphasizes the importance of maintaining efficient rodent control in order to safeguard public health.
Anticoagulant pesticides are commonly used in agricultural and urban rodent controls since few decades ago5,6. The anticoagulant rodenticides including warfarin and coumarin derivatives function effectively via binding with the vitamin K epoxide reductase of the rodents7,8,9. Resistance to several anticoagulant rodenticides including has been reported worldwide since 1960s10,11,12,13, and the vitamin recycling gene Vitamin K epoxide reductase complex subunit 1 (VKORC1) is now known to associate with the anticoagulant rodenticides-resistance14,15,16.
Studying mutations of the exonic nucleotide composition or single nucleotide polymorphisms (SNPs) of VKORC1 gene provides crucial information on resistance to As rodenticides and efficacy of pest control17,18,19. For instance, ~ 70% of sampled common or Norwegian rats (R. norvegicus) in the United Kingdom carried one of the five known missense mutations (Tyr139Cys, Tyr139Ser, Tyr139Phe and Leu128Gln and Leu120Gln)20, while these mutations could confer certain level of resistance to both first and second generation of anticoagulant rodenticides16,17,19. In a recent VKORC1-based SNP survey in mice and rats in the United States, it has also been suggested that resistances detected in the 1980s were likely due to mutations of Leu128Ser and Tyr139Cys in house mice (M. musculus domesticus), Arg35Pro in common or Norwegian rats (R. norvegicus), and Tyr25Phe in roof rats (R. rattus)21. Nevertheless, limited information was obtained from Asia, including Hong Kong. We therefore collected rodents from Hong Kong and carried out a VKORC1-based survey to estimate the anticoagulant rodenticide resistance situation that could compromise pest management.
Materials and methods
Sampling and DNA extraction
A total of 202 tail samples from dead rodents were provided to The Chinese University of Hong Kong by the Food and Environmental Hygiene Department, The Government of the Hong Kong Special Administrative Region and the City University of Hong Kong. The rodents were captured using traps from different locations in Hong Kong between 2017 and 2021. Tail samples were stored at − 20 °C before further experimental procedures. Genomic DNA extraction was carried out using QIAamp DNA mini kit (QIAgen, Germany) following the manufacturer’s instructions. In brief, 0.02 g of tail tissue were homogenized and incubated with proteinase K at 55 °C for 2 h. The quantity and quality of DNA were determined by Nanodrop (Ratio of 260/280 ~ 1.8 and 260/230 ~ 1.8–2.0) and gel electrophoresis under Gel Doc™ EZ imager (Bio-Rad), respectively.
Species identification
Molecular identification was carried out via the polymerase chain reaction (PCR) of mitochondrial DNA cytochrome c oxidase subunit 1 (COX1) gene using a model of T100™ thermocycler (Bio-Rad). COX1 gene was amplified using rodent specific primer BatL5310 (5′‐CCT ACT CRG CCA TTT TAC CTA TG‐3′) and R6036R (5′-ACT TCT GGG TGT CCA AAG AAT CA‐3′)22 with following parameters: 3 min of denaturation at 95 °C,39 cycles of 30 s at 95 °C, 30 s at 57 °C, and 40 s at 72 °C; and 5 min of final extension at 72 °C. Each reaction consisted of DNA sample (~ 10–20 ng), 1 × PCR buffer, 0.8 mM of dNTPs, 1.5 mM of MgCl2, 0.4 μM of each forward and reverse primers, 11.2 μL of dd H2O and 1 unit of Taq DNA polymerase. The amplified products (762 bp) were confirmed on 1% agarose gel stained as well as Sanger sequencing (BGI Genomics Company Hong Kong). Obtained sequences were edited with software SnapGene Viewer, and aligned using MEGA X for phylogenetic analysis (Neighbour-joining method and 1000 bootstrap replications).
VKORC1 sequence analysis
All three exons of VKORC1 gene were amplified following a previous study using specific primers: (Exon1 forward: 5′-GTG GCG GGT TCT TCC CTC-3′; Exon 1 reverse: 5′-GAC TCC AAA ATC ATC TGG CAA CC-3′), (Exon 2 forward: 5′-AAG AGT AGG GGAC AAG GTG GC-3′; Exon 2 reverse: 5′-GGG TCA CCA AGA CAT GAG GTG-3′) and (Exon 3 forward: 5′-TTT CAC CAG AAG CAC CTG CTG CC-3′; Exon 3 reverse: 5′-ACA CTT GGG CAA GGC TCA TGT G-3′)13. The amplified products were confirmed on 2% agarose gel stained as well as Sanger sequencing (BGI Genomics Company Hong Kong). SNP of each exon sequence was compared to the available sequence from NCBI database (VKORC1 GenBank accession no. AY423047) with MEGA X software. BlastX searches with adjusted sequences were also carried out to locate any missense mutation. Homozygous and heterozygous genotypes of five published missense mutations on exon 3 were further confirmed on each chromatogram using SnapGene Viewer.
Genome sequencing of selected individuals
DNA of R. norvegicus and R. tanezumi from 4 localities including Yuen Long (YL_2, YL_3), Wan Chai (Wch_1, Wch_2), Kwun Tong (KTo_4, KTo_5) and Islands (Is_1, Is_6) were proceeded with low-coverage whole genome sequencing (Table 1). Raw sequenced reads were mapped to the R. norvegicus reference genome (GenBank assembly accession: GCF_000001895.5) and SNPs were called with Genome Analysis Toolkit (GATK)23. The SNP dataset was annotated with the gene models of the reference assembly using SnpEff24. The NGS data have been uploaded to NCBI under the BioProject accession number PRJNA723168.
Results
VKORC1 exon 3 of Rattus losea, R. norvegicus, and R. tanezumi in Hong Kong
In the 202 collected rats, 52, 81, and 69 of them were, R. losea, R. norvegicus, and Rattus tanezumi, respectively (Fig. 1). Greater genetic diversity was also observed in the COX1 of R. norvegicus than the two other captured species.
Among the five previous reported mutations in VKORC1 exon 3 reported elsewhere in the world (Tyr139Cys, Tyr139Ser, Tyr139Phe and Leu128Gln and Leu120Gln), only Tyr139Cys mutation was found in the R. tanezumi samples but not in the other collected species.
In the 69 collected R. tanezumi, 47 of them (68.1%) were found to carry Tyr139Cys mutations with 25 homozygotes and 22 heterozygotes. Details of their sampling locations and number of mutations are summarised in Table 2 and Fig. 2.
Other SNPs on the VKORC1 gene
Table 3 summarised all the located SNPs obtained from the selected R. tanezumi and R. norvegicus samples subjected to whole genome sequencing based on their geographical distributions.
In addition to the known Y139C mutation, nonsynonymous mutations were also found from one R. norvegicus sample and one R. losea sample, respectively. Further, six synonymous mutations were also found among three species. Details are provided in Supplementary information S1.
Besides the exons, a total of nine SNPs was revealed locating at the introns, with three coming from R. norvegicus and the other six from R. tanezumi (Table 4).
Ethics declaration
Animal ethics approval was granted by the Animal Research Ethics Sub-Committee of City University of Hong Kong. All methods were carried out in accordance with relevant guidelines and regulations. All methods are reported in accordance with ARRIVE guidelines.
Discussion
Efficacy of rodenticides is crucial for effective rodent management, and this study carried out the first VKORC1-based survey to estimate the anticoagulant rodenticide resistance situation. In contrast to the previous rodent species identification in Hong Kong revealing eight species of rats and mice, with Rattus norvegicus and R. rattus to be the dominant rat species in urban areas3. This study, nevertheless, identified three Rattus species including the report of the R. losea, R. norvegicus, and R. tanezumi based on molecular marker COX1. The number of captured rats has revealed the abundance of R. tanezumi and R. norgevicus, while the R. losea were captured from two locations only. Despite R. rattus and R. tanezumi were well known to be difficult to be morphologically differentiated from one another22,25, given the previous and present studies were carried out at different time (more than ten years) and places using different collection method, it is unclear whether the situation represents misidentification, distribution in different biotopes, different collection methods, or changes in dominant rodent species spatiotemporally.
In the limited studies carried out on anticoagulant rodenticide resistance in Asia, a relatively low warfarin-resistance rate (11%, 4 out of 36 samples) was determined by lethal feeding test in R. tanezumi collected from mainland China ten years ago26. It should be noticed that the use of anticoagulant rodenticides in China was believed to have started in the early 1980s27, which has a shorter history than other places in the world. A recent study also suggested a low anticoagulant rodenticide resistance rate in R. norvegicus collected from two cities in mainland China28. This study, based on the VKORC1 gene survey, discovered 68.1% of R. tanezumi in Hong Kong carried the Tyr139Cys mutation.
Previous studies suggested that the Tyr139Cys mutation could confer resistance to first- and second-generation anticoagulant rodenticides including bromadiolone and difenacoum in Norway rat and house mouse19. Given the relationships between anticoagulant rodenticide resistance and the Tyr139Cys mutation in R. tanezumi has not been tested, the cause and significance of such mutation being only observed in R. tanezumi but not in R. losea and R. norvegicus remains to be revealed. In case if the Tyr139Cys mutation in R. tanezumi also confer certain type of anticoagulant rodenticide resistance, other substances such as difethialone and flocoumafen could to be used29. Regarding to the rodent nuisance in Hong Kong, anticoagulant compound is more desirable and safer rodenticide for controlling rodents compared with acute poison within the densely populated urban area. Anticoagulant compound has been widely adopted by both private and public pest control operators. Currently, there is no statutory regulation to monitor the use of rodenticide from local pest control operators, however, the low efficacy of certain compound and good prevention practice should be aware in order to decrease the influence of rodent problems.
This study also revealed other SNP variants not documented previously, for instances, two synonymous SNPs and one nsSNP (Trp101Cys) in R. losea. It is also worth noting that no SNPs located in exon 1 identified from the nine R. losea samples were Arg58Gly mutation which confer anticoagulant rodenticide resistance30. These data bring up the issues that there are huge gaps in knowledge regarding the origin, introduction, genetic diversity, population connectivity of Rattus between different places in Asia, as well as the relationships of mutations brought in to VKORC1 genes and their anticoagulant rodenticide resistance of different Rattus population from different places in Asia.
Conclusion
This study provided the baseline information of rodenticide resistance status and distribution of 202 rodents belonging to 3 Rattus species in Hong Kong. The investigation indicates a distinctive anticoagulant rodenticide resistance pattern. The relatively high Tyr139Cys mutation found in VKORC1 gene of R. tanezumi suggested further susceptibility tests will be needed to reveal whether they are resistance to individual anticoagulant rodenticide and to ensure effectiveness on local species. Regular genetic testing and genomic analyses will also be required to understand the situations of rodent populations for integrated pest management.
Data availability
The raw reads generated in this study have been deposited to the NCBI database under the BioProject accession PRJNA723168.
References
Perry, R. D. & Fetherston, J. D. Yersinia pestis-etiologic agent of plague. Clin. Microbiol. Rev. 10, 35–66 (1997).
Meerburg, B. G., Singleton, G. R. & Kijlstra, A. Rodent-borne diseases and their risks for public health. Crit. Rev. Microbiol. 35, 221–270 (2009).
Chung, K. P. S. & Corlett, R. T. Rodent diversity in a highly degraded tropical landscape: Hong Kong, South China. Biodivers. Conserv. 15, 4521–4532 (2006).
Sridhar, S. et al. Transmission of rat hepatitis E virus infection to humans in Hong Kong: A clinical and epidemiological analysis. Hepatology 73, 10–22 (2020).
Hadler, M. & Buckle, A. Forty five years of anticoagulant rodenticides—Past, present and future trends. In Proceedings of the Fifteenth Vertebrate Pest Conference 1992, Vol. 36, 149–155 (1992).
Watt, B. E., Proudfoot, A. T., Bradberry, S. M. & Vale, J. A. Anticoagulant rodenticides. Toxicol. Rev. 24, 259–269 (2005).
Whitlon, D. S., Sadowski, J. A. & Suttie, J. W. Mechanism of coumarin action: Significance of vitamin K epoxide reductase inhibition. Biochemistry 17, 1371–1377 (1978).
Stafford, D. W. The vitamin K cycle. J. Thromb. Haemost. 3, 1873–1878 (2005).
Tie, J. K., Nicchitta, C., von Heijne, G. & Stafford, D. W. Membrane topology mapping of vitamin K epoxide reductase by in vitro translation/cotranslocation. J. Biol. Chem. 280, 16410–16416 (2005).
Boyle, C. M. Case of apparent resistance of Rattus norvegicus Berkenhout to anticoagulant poisons. Nature 188, 517 (1960).
Berny, P., Esther, A., Jacob, J. & Prescott, P. In Development of Resistance to Anticoagulant Rodenticides in Rodents, Chapter 10 in Anticoagulant Rodenticides and Wildlife Vol. 5 (eds van den Brink, N. W. et al.) 259–286 (Springer, 2018).
Hodroge, A., Longin-Sauvageon, C., Fourel, I., Benoit, E. & Lattard, V. Biochemical characterization of spontaneous mutants of rat VKORC1 involved in the resistance to antivitamin K anticoagulants. Arch. Biochem. Biophys. 515, 14–20 (2011).
Grandemange, A., Lasseur, R., Longin-Sauvageon, C., Benoit, E. & Berny, P. Distribution of VKORC1 single nucleotide polymorphism in wild Rattus norvegicus in France. Pest Manag. Sci. 66, 270–276 (2010).
Li, T. et al. Identification of the gene for vitamin K epoxide reductase. Nature 427, 541–544 (2004).
Rost, S. et al. Mutations in VKORC1 cause warfarin resistance and multiple coagulation factor deficiency type 2. Nature 427, 537–541 (2004).
Pelz, H. J. et al. Distribution and frequency of VKORC1 sequence variants conferring resistance to anticoagulants in Mus musculus. Pest Manag. Sci. 88, 254–259 (2012).
Rost, S. et al. Novel mutations in the VKORC1 gene of wild rats and mice—A response to 50 years of selection pressure by warfarin?. BMC Genet. 10, 4 (2009).
Pelz, H. J. et al. The genetic basis of resistance to anticoagulants in rodents. Genetics 170, 1839–1847 (2005).
McGee, C. F., McGilloway, D. A. & Buckle, A. P. Anticoagulant rodenticides and resistance development in rodent pest species: A comprehensive review. J. Stored Prod. Res. 88, 101688–101688 (2020).
Jones C, Talavera M, Buckle A and Prescott C, Anticoagulant resistance in rats and mice in the UK—Summary report with new data for 2019. Report from the Campaign for Responsible Rodenticide Use (CRRU) UK for the Government Oversight Group, Vertebrate Pests Unit. The University of Reading. Accessed 18 March 2021 https://www.thinkwildlife.org/downloads/ (2019).
Díaz, J. C. & Kohn, M. H. A VKORC1-based SNP survey of anticoagulant rodenticide resistance in the house mouse, Norway rat and roof rat in the USA. Pest Manag. Sci. 77, 234–242 (2020).
Robins, J., Hingston, M., Matisoo-Smith, E. & Ross, H. Identifying Rattus species using mitochondrial DNA. Mol. Ecol. Notes 7, 717–729 (2007).
DePristo, M. A. et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491–498 (2011).
Cingolani, P. et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 6, 80–92 (2012).
Aplin, K. P. et al. Multiple geographic origins of commensalism and complex dispersal history of black rats. PLoS ONE 6, 11 (2011).
Huang, B. H. et al. Warfarin resistance test and polymorphism screening in the VKORC1 gene in Rattus flavipectus. J. Pest Sci. 84, 87–92 (2011).
Liang, L. The Resistance of Rattus flavipectus and R. norvegicus to anticoagulant rodenticide in Zhanjiang Proper. Chin. J. Vector Biol. Control 16, 21–22 (2005) (in Chinese).
Ma, X. H. et al. Low warfarin resistance frequency in Norway rats in two cities in China after 30 years of usage of anticoagulant rodenticides. Pest Manag. Sci. 74, 2555–2560 (2018).
RRAG. Anticoagulant resistance in the Norway rat and guidelines for the management of resistant rat infestations in the UK. Rodenticide Resistance Action Group, UK. Revision Septemeber 2018. Accessed 4 June 2019 https://www.pestmagazine.co.uk/media/246897/management-of-resistant-norwayratinfestations-in-the-uk-rrag-june-2010.pdf (2018).
Wang, J. et al. Warfarin resistance in Rattus losea in Guangdong Province, China. Pest. Biochem. Physiol. 91, 90–95 (2008).
Acknowledgements
The authors would like to thank Susan Lu and Nicole Leong for carrying out some of the DNA extraction. The authors would like to thank K.L. Chan and Ming Wai Lee of Pest Control Advisory Section of the Food and Environmental Hygiene Department for their discussion in this study. This work was supported by The Chinese University of Hong Kong and Hong Kong Research Grant Council Collaborative Research Fund (C4015-20EF).
Author information
Authors and Affiliations
Contributions
J.H.L.H. conceived the study. J.H.L.H. supervised the study. E.Y.Y.H. carried out the PCR and sequence analyses. S.T.S.L. and W.N. conducted the population genomic analyses. H.Y.Y. provided the logistics support. TUA and IM collected the animals from field. E.Y.Y.H. and J.H.L.H. wrote the first draft of the manuscript. All authors (E.Y.Y.H., S.T.S.L., W.N., H.Y.Y., T.U.A., I.M., J.H.L.H.) approved the final version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare 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
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/.
About this article
Cite this article
Huang, E.Y.Y., Law, S.T.S., Nong, W. et al. The screening for anticoagulant rodenticide gene VKORC1 polymorphism in the rat Rattus norvegicus, Rattus tanezumi and Rattus losea in Hong Kong. Sci Rep 12, 12545 (2022). https://doi.org/10.1038/s41598-022-16550-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-16550-3
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.