Next Article in Journal
Molecular Analysis of SARS-CoV-2 Lineages in Armenia
Next Article in Special Issue
Non-Lethal Detection of Frog Virus 3-Like (RUK13) and Common Midwife Toad Virus-Like (PDE18) Ranaviruses in Two UK-Native Amphibian Species
Previous Article in Journal
Fecal Feline Coronavirus RNA Shedding and Spike Gene Mutations in Cats with Feline Infectious Peritonitis Treated with GS-441524
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 5
10.3390/v14051073
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prevalence of Ranavirus Infection in Three Anuran Species across South Korea

by
Namho Roh
1,†,
Jaejin Park
2,†,
Jongsun Kim
2,
Hyerim Kwon
2 and
Daesik Park
2,*
1
Department of Biological Sciences, Kangwon National University, Chuncheon 24341, Kangwon, Korea
2
Division of Science Education, Kangwon National University, Chuncheon 24341, Kangwon, Korea
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2022, 14(5), 1073; https://0-doi-org.brum.beds.ac.uk/10.3390/v14051073
Submission received: 25 April 2022 / Revised: 15 May 2022 / Accepted: 16 May 2022 / Published: 17 May 2022
(This article belongs to the Special Issue Recent Advances in Virus Infections of Lower Vertebrates)
Browse Figures
Review Reports Versions Notes

Abstract

:
To cope with amphibian die-offs caused by ranavirus, it is important to know the underlying ranavirus prevalence in a region. We studied the ranavirus prevalence in tadpoles of two native and one introduced anuran species inhabiting agricultural and surrounding areas at 49 locations across eight provinces of South Korea by applying qPCR. The local ranavirus prevalence and the individual infection rates at infected locations were 32.6% and 16.1%, respectively, for Dryophytes japonicus (Japanese tree frog); 25.6% and 26.1% for Pelophylax nigromaculatus (Black-spotted pond frog); and 30.5% and 50.0% for Lithobates catesbeianus (American bullfrog). The individual infection rate of L. catesbeianus was significantly greater than that of D. japonicus. The individual infection rate of P. nigromaculatus was related to the site-specific precipitation and air temperature. The individual infection rate gradually increased from Gosner development stage 39, and intermittent infection was confirmed in the early and middle developmental stages. Our results show that ranavirus is widespread among wild amphibians living in agricultural areas of South Korea, and mass die-offs by ranavirus could occur at any time.

1. Introduction

Ranavirus is an infectious agent that affects ectothermic organisms such as fish, amphibians, and reptiles [1,2,3]. Ranavirus is known as one of the main causes of the global decline in amphibians, in addition to chytrid fungi [4,5,6,7]. As a major threat to various amphibians, ranavirus, for example, has caused mass die-offs of the endangered giant salamander in China [8], the larval group of endangered gold-spotted pond frogs in Korea [9], and the abundance of common frogs in England [10]. Despite being a continuous threat, ranaviruses are much less well studied than the amphibian chytrid fungus [11,12,13]. Since ranavirus was first discovered in the United States [14], it has spread and is now found on five continents [3,7,15]. Recently, the number of cases of discovery and mass die-offs due to ranavirus infection has also increased in Asia.
Among Asian countries, China has had the highest number of published studies, with 153 registered publications in the Web of Science database between 2010 and 2019 [10]. In addition to China, cases of ranavirus infection or mass die-offs due to ranavirus infection have been reported in eight Asian countries and regions, including Taiwan, Hong Kong, Japan, Korea, Thailand, Singapore, Malaysia, and India [16]. To date, ranavirus infection or mass mortality reports have been made for five amphibian species in South Korea, including larval Rana huanrenensis [17], adult R. uenoi [18], adult Kaloula borealis [19], larval P. chosenicus [12], and larval D. japonicus [19]. To predict the possibility of ranavirus-caused mass die-offs according to species and region, it is important to determine the ranavirus prevalence in various amphibian species in a wide area. While ranavirus-caused mass die-offs have been frequently reported in various countries [20,21], few studies have screened ranavirus prevalence on a regional scale [21,22,23]. Even within China, where the most research has been performed, screening was performed only in Heilongjiang Province and some areas of northeastern China, and the studies included only a few species, such as R. dybowskii and R amurensis [24,25]. Such background studies have not been conducted in most Asian countries.
To determine the ranavirus prevalence in amphibian species inhabiting agricultural and surrounding areas, such as wetlands, ponds, and reservoirs, we screened the ranavirus prevalence in three anuran species: the introduced L. catesbeianus and the native D. japonicus and P. nigromaculatus. Our study covered 49 locations across eight provinces in South Korea. In addition, the relationship between the individual ranavirus infection rate and habitat characteristics was explored.

2. Materials and Methods

2.1. Sampling Locations and Tissue Collection

The selected three species (D. japonicus, P. nigromaculatus, and L. catesbeianus) are commonly distributed in agricultural areas in South Korea [26]. We selected L. catesbeianus because it is an introduced species widely inhabiting the central and southern parts of the Korean Peninsula [27], and there have been many reports of ranavirus infection in this species in many countries [28,29,30]. Dryophytes japonicus and P. nigromaculatus largely share distribution ranges and habitats with L. catesbeianus, and they are also representative anuran species in rice paddies [26]. To identify sampling locations, we first selected the areas where all three species or at least two native species have been reported within a 3 km radius based on the results of the third (2006–2012) and fourth (2013–2018) “National Ecosystem Surveys (NESs)”. The South Korean Ministry of Environment has conducted three NESs since 1986, and the data are available from EcoBank [31]. After the confirmation of species presence by directly visiting the locations, we finally selected 49 locations across eight provinces (Figure 1; Table 1). Considering the effective collection of samples over a wide range, we collected tadpoles of the three species using a dip net, seine or trap between June 2020 and October 2021. To prevent potential transmission of ranavirus between sampling locations, we sanitized all collecting equipment using the following protocol: brushing off mud and vegetation, spraying 10% bleach (Yuhan Clorox, Seoul, South Korea) and cleaning with tap water, spraying 70% ETOH, and completely air-drying all field equipment after each collection [32]. We individually preserved sampled tadpoles (Voucher number G01621DJ-G03995LC) in 99% ETOH after euthanasia by submerging them for more than 15 min in 0.5% MS222 [33]. In the laboratory, we randomly subsampled a maximum of eight tadpoles of each species from each location, determined their Gosner developmental stage [34], extracted liver tissues, and preserved the tissues at −80 °C.

2.2. DNA Extraction and qPCR

DNA extraction from liver tissue was performed with the DNeasy Blood and Tissue kit (Qiagen, Hilden, Germany) following the manufacturer’s protocol. The extracted DNA was quantified with a Qubit3 Fluorometer (Invitrogen, Waltham, MA, USA) using the Qubit 1X dsDNA HS Assay Kit (Invitrogen, Waltham, MA, USA) and then stored at −80 °C until qPCR experiments. We used RVMCPKim3_F and RVMCPKim3_R qPCR primers to determine ranavirus infection [35]. The composition of the amplification reaction solution for qPCR included 10 µL of power SYBR green PCR master mix (Applied Biosystems, Waltham, MA, USA), 0.5 µL of forward primer, 0.5 µL of reverse primer, and 4 ng of DNA. We adjusted the volume to a final value of 20 µL using molecular biology grade water. qPCR was performed in QuantStudio 1 (Applied Biosystems, Waltham, MA, USA) at 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 62.5 °C for 20 s. qPCR of all samples was performed in triplicate with a negative control (sterile, molecular grade water) and a positive control. The positive control was PCR-amplified ranavirus MCP DNA from Kaloula borealis, and this sample was obtained and verified for ranavirus infection in 2016 [19]. If a positive reaction was detected in two or more wells of the three replicate samples, the melting temperature (Tm) value of the melting curve coincided with the positive control, and the cycle threshold (CT) value was 35 or less, then we considered the tadpole to be infected with ranavirus (Figure S1A–H). If a positive reaction was confirmed in only one well, the test was rerun, and only when the above conditions were satisfied was it judged to be an infected tadpole.

2.3. Weather Factors and Habitat Characteristics

To investigate the relationship between the ranavirus infection rate and habitat characteristics, we analyzed 16 habitat characteristics at each sampling location. The land cover rate (urban area, agricultural area, forest area, grass area, wetland area, bare area, hydrosphere area) within a radius of 3 km2 from the mid-sampling point was calculated on the 2020 South Korean land cover map [36]. From the land cover data, greenhouse cultivation areas within agricultural areas, salt fields and tidal flat areas within wetland areas, and sea areas within hydrosphere areas were excluded during data handling because these areas are inhabitable for amphibians. Additionally, the shortest distance from the mid-sampling point to agricultural land, mountains, water bodies, and urban areas was calculated in units of 1 m. When a sampling point was located within a specific coverage area, the distance was set to 2 m. All land cover rates, distance data, and altitude values of each sampling point were calculated using QGIS (ver. 3.4.7, QGIS.org 2021; https://www.qgis.org/ko, accessed on 27 July 2021). Furthermore, the average air temperature, lowest air temperature, highest air temperature, and average precipitation for the immediately preceding quarter, based on the sampling date of the tadpoles used in the qPCR experiment, were obtained from the local meteorological station closest to the sampling point.

2.4. Data Analyses

For analyses, we log-transformed 16 habitat and climate characteristics. As most of the data did not present a normal distribution even after transformation, nonparametric statistics were used for statistical analyses. We used the chi-square test to determine the difference in the prevalence of infected locations (local ranavirus prevalence) among the three species using the Kruskal–Wallis test, and the Dunn-Bonferroni post hoc test was conducted following the provided option of the test. The correlations of the individual ranavirus infection rates within the three species and between the three species and the 16 habitat characteristics were verified by Spearman correlation analysis. All statistics were performed in SPSS (Version 26, IBM). To determine the change in the individual ranavirus infection rate according to the Gosner developmental stage, the number of ranavirus-infected individuals was calculated at each stage and then compared to the number of individuals, which was qPCR-tested at each stage for each species. Afterward, the average ranavirus infection rate of the three species based on each developmental stage was obtained, and trends were presented. All data in the text are presented as the mean ± 1 standard error unless otherwise noted.

3. Results

Out of 49 populations, we collected larval D. japonicus from 43 locations (22.1 ± 4.2 individuals per location, n = 43), larval P. nigromaculatus from 43 locations (21.0 ± 4.0, n = 43), and larval L. catesbeianus from 26 locations (19.5 ± 10.2, n = 26) (Figure 1). The Gosner developmental stage of the collected larvae ranged from 25 to 46. In the qPCR experiments, we subsampled and used 344 larval D. japonicus (8 tadpoles per location), 342 larval P. nigromaculatus (8 per location, except for 6 in one location), and 187 larval L. catesbeianus (7.2 ± 2.1 SD per location, ranging from 1–8 individuals).

3.1. Rates of Local and Individual Ranavirus Infections

We found at least one ranavirus-infected individual at 28 of the 49 studied locations (57.1%) (Figure 1; Table 1). The local ranavirus prevalence was 32.6% (14 out of 43 locations), 25.6% (11 out of 43 locations), and 30.8% (8 out of 26 locations), and the individual infection rates within the infected locations were 16.1 ± 2.0% (12.5–37.5%, n = 14), 26.1 ± 5.7% (12.5–75.0%, n = 11), and 50.0% 50.0 ± 9.7% (12.5–100%, n = 8) for D. japonicus, P. nigromaculatus, and L. catesbeianus, respectively (Table 1). The local ranavirus prevalences were not significantly different between the three species (p = 0.77), but the individual ranavirus infection rates at the infected locations were significantly different (H (2) = 12.47, p = 0.002) (Figure 2). In particular, the rate of L. catesbeianus was greater than that of D. japonicus (p = 0.001). In contrast, the rates were not significantly different between L. catesbeianus and P. nigromaculatus or between P. nigromaculatus and D. japonicus (ps > 0.05).

3.2. Relationship with Weather Conditions, Habitat Characteristics, and Developmental Stages

The individual infection rate was positively correlated with grass area (r = 0.579, n = 14, p = 0.030) in D. japonicus and negatively correlated with distance to mountains (r = −0.604, n = 11, p = 0.049) in P. nigromaculatus and with forest area (r = −0.781, n = 8, p = 0.022) in L. catesbeianus. In P. nigromaculatus, the rate was also positively related to the precipitation (r = 0.634, n = 11, p = 0.036), average air temperature (r = 0.837, n = 11, p = 0.001), and lowest air temperature (r = 0.721, n = 11, p = 0.012) in the first quarter prior to the sampling date. The individual ranavirus infection rates of the three species were not related (ps > 0.05). Additionally, all other remaining relationships were not significant (ps > 0.05). The individual infection rate of the three species showed a gradual increase from Gosner developmental stage 39. In the early and middle developmental stages, there were intermittent cases of ranavirus infection, such as at stages 32 and 37 (Figure 3).

4. Discussion

Our study shows that ranavirus is widely spread in amphibians inhabiting agricultural and surrounding areas in South Korea. Ranavirus-infected tadpoles were identified in over 25.6% of the surveyed locations, and the average individual infection rate in the infected area was as high as 50% depending on the species. Amphibian larvae often have a higher rate of ranavirus infection than adults [37] due to factors such as immunologically vulnerable metamorphosis and high population density [38,39]. To date, although ranavirus infection in larval D. japonicus has been confirmed in agricultural areas [19], there has been no report of mass die-offs of either larval or adult D. japonicus or P. nigromaculatus in South Korea. Our results show that given the high 50% individual infection rate of L. catesbeianus, amphibian mass die-offs due to ranavirus could occur in agricultural areas whenever outbreak conditions are present.
The high individual infection rate of L. catesbeianus tadpoles may be due to the habitat and ecological characteristics of the species. Ranavirus can infect various ectothermic vertebrates through direct individual contact or indirect water matrices [40,41,42]. In this study, the individual infection rate of L. catesbeianus was 50.0%, which was higher than the rates of 16.1% and 26.1% for D. japonicus and P. nigromaculatus, respectively. Several factors might have affected this result. First, while the larvae of D. japonicus and P. nigromaculatus generally spend approximately one month in water before metamorphosis [26], the larvae of L. catesbeianus often spend more than one year in water [43]. This difference may increase the likelihood of being directly or indirectly infected by ranavirus. Second, considering that the main habitats of L. catesbeianus are ponds and reservoirs adjacent to farmland and villages, there is a high possibility of long-term exposure to anthropogenic pollutants. Such pollutants are often related to ranavirus infection [44,45]. Third, adult L. catesbeianus live in ponds throughout the year. In contrast, adult D. japonicus and P. nigromaculatus live in water only during the breeding season [26]. Therefore, adult L. catesbeianus are highly likely to be infected with ranavirus, which may have subsequently caused the high rate of larval infection. In future studies, screening the degree of ranavirus infection in adult L. catesbeianus is necessary.
In this study, we could not find any association of ranavirus infection between introduced L. catesbeianus and two native species (D. japonicus and P. nigromaculatus). Despite screening 49 locations, we confirmed only four locations with more than two species infected by ranavirus. It has been reported that L. catesbeianus functions as a ranavirus carrier and transmission vector [46,47]. Lithobates catesbeianus was introduced to Korea approximately 50 years ago and is widely distributed in the central and southern regions of South Korea [27]. In addition, considering that ranavirus has been previously detected in bullfrogs in many countries [28,29,30], L. catesbeianus may be responsible, at least in part, for the wide ranavirus prevalence in South Korea. However, in this study, joint ranavirus infections of both L. catesbeianus and P. nigromaculatus were confirmed at only one location. Therefore, we could not appropriately test the hypothesis that L. catesbeianus may serve as a ranavirus vector for native anuran species. In future studies, screening for ranavirus and comparing the strains through sequencing in amphibians, fish, and insects that share a water source with or cohabit with L. catesbeianus are necessary.
The relationships between the individual ranavirus infection rate and several habitat characteristics were identified. In contrast with the cases in L. catesbeianus, the infection rate of P. nigromaculatus showed a high positive correlation with precipitation, average air temperature, and lowest air temperature in the first quarter prior to the collection date. This result was consistent with those of previous studies showing that the rate of ranavirus infection was related to habitat temperature [2]. In addition, in P. nigromaculatus, the greater the distance from the mountain was, the higher the infection rate, implying that the infection rate might be low when the anthropogenic factors are low near mountainous or forest areas. However, in L. catesbeianus, there was no correlation with climatic factors, and only the more forested areas had a lower infection rate. This result suggests that ranavirus infection in the introduced L. catesbeianus might be related to anthropogenic factors such as various pollutants rather than to climatic factors. In D. japonicus, the infection rate at 12 of the locations was the same, at 12.5%, so it might be statistically difficult to detect any significant relationships with habitat characteristics. Many amphibians can survive in a state of ranavirus infection [15,46], and many field studies have conducted post hoc mass mortality analysis [17,18,28]. Thus, the direct causal relationship between population habitat characteristics and ranavirus-caused mass mortality has not been clearly established. Further studies are needed to link the rate of infection with habitat characteristics.
There were two trends in the infection rate according to the developmental stage of tadpoles. The individual ranavirus infection rate gradually increased after Gosner developmental stage 39, and intermittent ranavirus infection occurred in the early to mid-developmental stages. Among the factors influencing the ranavirus infection rate in amphibians, the developmental stage is a well-known factor [2,3,4,5]. The ranavirus infection rate is often high at Gosner stages 44–46, which is when metamorphosis occurs and often leads to die-off events [38,48,49]. There was also a case report in South Korea where many larval R. huanrenensis died during metamorphosis [17]. Although ranavirus-caused die-off events of amphibian tadpoles are often found close to metamorphosis [50,51], tadpoles at lower developmental stages are also susceptible to ranavirus infection [52]. Our results imply that anuran tadpoles can be infected by ranavirus regardless of the specific developmental stage and that ranavirus-caused mass die-offs can occur whenever environmental triggers form suitable conditions for an outbreak. In the paddy fields in which rice is cultivated, water flooding and draining are repeated to help the successful growth of rice, and this cycle greatly affects amphibian activities [39,53]. Often, several thousand dead tadpoles at various developmental stages are found in small hollow patches in water-drained rice paddies in summer [54]. Our findings suggest a potential link between these deaths and ranavirus, suggesting the urgent need for further studies.

5. Conclusions

Our study is of great significance in that we, for the first time, determined the ranavirus prevalence in representative amphibian species living in agricultural and surrounding areas across South Korea. Our results show that mass die-offs of amphibians due to ranavirus could occur at any time in such areas. In particular, the further screening of adult frogs and the investigation of possible ranavirus causation on large die-offs of tadpoles in rice fields are urgently necessary.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/v14051073/s1, Figure S1: The representative qPCR results of D. japonicus (A,B) collected at Haenam, P. nigromaculatus (C,D) collected at Pohang, and L. catesbeianus (E,F) collected at Gimchoen, and the photographs of representative gel electrophoresis of the MCP fragment (~520 bp), which confirmed the qPCR result again by sequencing the fragment (G,H) in the further study.

Author Contributions

Conceptualization, N.R., J.P. and D.P.; methodology, J.P., J.K. and D.P.; formal analysis, J.P. and D.P.; investigation, N.R., J.P., H.K. and D.P.; data curation, J.P., J.K. and H.K.; writing—original draft preparation, N.R., J.P. and J.K.; writing—review and editing, J.K. and D.P.; supervision, D.P.; funding acquisition, J.P. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (No. 2020R1A6A3A13060949).

Institutional Review Board Statement

This study was reviewed and approved by the Institutional Animal Care and Use Committee of Kangwon National University (KW-200618-3).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data on the weather and habitat characteristics, collected from 49 sampling sites and the representative qPCR results and the representative photographs of the gel electrophoresis of MCP fragment (~520 bp), which were taken during further sequencing procedures, are openly available in Dryad Digital Repository at https://datadryad.org/stash/dataset/doi:10.5061/dryad.0zpc8670q (accessed on 11 May 2021).

Acknowledgments

We thank Woo-Jin Choi, Il-Kook Park, Ji-Ho Park for their help during field sampling and MinWoo Park, Ho-Jun Jung for their help during the dissection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hoverman, J.T.; Gray, M.J.; Haislip, N.A.; Miller, D.L. Phylogeny, life history, and ecology contribute to differences in amphibian susceptibility to ranaviruses. EcoHealth 2011, 8, 301–319. [Google Scholar] [CrossRef] [PubMed]
  2. Brunner, J.L.; Storfer, A.; GrayJason, M.J.; Hoverman, T. Ranavirus ecology and evolution: From epidemiology to extinction. In Ranaviruses: Lethal Pathogens of Ectothermic Vertebrates; Gray, M.J., Chinchar, V.G., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 71–104. ISBN 978-3-319-13754-4. [Google Scholar]
  3. Duffus, A.L.; Waltzek, T.B.; Stöhr, A.C.; Allender, M.C.; Gotesman, M.; Whittington, R.J.; Hick, P.; Hines, M.K.; Marschang, R.E. Distribution and host range of ranaviruses. In Ranaviruses: Lethal Pathogens of Ectothermic Vertebrates; Gray, M.J., Chinchar, V.G., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 9–57. ISBN 978-3-319-13754-4. [Google Scholar]
  4. Rosa, I.D.; Simoncelli, F.; Fagotti, A.; Pascolini, R. The proximate cause of frog declines? Nature 2007, 447, E4–E5. [Google Scholar] [CrossRef] [PubMed]
  5. Robert, J. Emerging ranaviral infectious diseases and amphibian decline. Diversity 2010, 2, 314–330. [Google Scholar] [CrossRef]
  6. Price, S.J.; Garner, T.W.; Nichols, R.A.; Balloux, F.; Ayres, C.; de Alba, A.M.C.; Bosch, J. Collapse of amphibian communities due to an introduced ranavirus. Curr. Biol. 2014, 24, 2586–2591. [Google Scholar] [CrossRef] [Green Version]
  7. Brunner, J.L.; Olson, D.H.; Gray, M.J.; Miller, D.L.; Duffus, A.L. Global patterns of ranavirus detections. Facets 2021, 6, 912–924. [Google Scholar] [CrossRef]
  8. O’hanlon, S.J.; Rieux, A.; Farrer, R.A.; Rosa, G.M.; Waldman, B.; Bataille, A.; Kosch, T.A.; Murray, K.A.; Brankovics, B.; Fumagalli, M.; et al. Recent Asian origin of chytrid fungi causing global amphibian declines. Science 2018, 360, 621–627. [Google Scholar] [CrossRef] [Green Version]
  9. Fisher, M.C.; Pasmans, F.; Martel, A. Virulence and pathogenicity of chytrid fungi causing amphibian extinctions. Annu. Rev. Microbiol. 2021, 75, 673–693. [Google Scholar] [CrossRef]
  10. Wirth, W.; Lesbarrères, D.; Ariel, E. Ten years of ranavirus research (2010–2019): An analysis of global research trends. Facets 2021, 6, 44–57. [Google Scholar] [CrossRef]
  11. Geng, Y.; Wang, K.Y.; Zhou, Z.Y.; Li, C.W.; Wang, J.; He, M.; Yin, Z.Q.; Lai, W.M. First report of a ranavirus associated with morbidity and mortality in farmed Chinese giant salamanders (Andrias davidianus). J. Comp. Pathol. 2011, 145, 95–102. [Google Scholar] [CrossRef]
  12. Kim, S.; Sim, M.Y.; Eom, A.H.; Park, D.S.; Ra, N.Y. PCR detection of ranavirus in gold-spotted pond frogs (Rana plancyi chosenica) from Korea. Korean J. Environ. Biol. 2009, 27, 110–113. [Google Scholar]
  13. Teacher, A.G.F.; Cunningham, A.A.; Garner, T.W.J. Assessing the long-term impact of ranavirus infection in wild common frog populations. Anim. Conserv. 2010, 13, 514–522. [Google Scholar] [CrossRef]
  14. Granoff, A.; Came, P.E.; Rafferty, K.A., Jr. The isolation and properties of viruses from Rana pipiens: Their possible relationship to the renal adenocarcinoma of the leopard frog. Ann. N. Y. Acad. Sci. 1965, 126, 237–255. [Google Scholar] [CrossRef]
  15. Miller, D.L.; Gray, M.J.; Rajeev, S.; Schmutzer, A.C.; Burton, E.C.; Merrill, A.; Baldwin, C.A. Pathologic findings in larval and juvenile anurans inhabiting farm ponds in Tennessee, USA. J. Wildl. Dis. 2009, 45, 314–324. [Google Scholar] [CrossRef] [Green Version]
  16. Herath, J.; Ellepola, G.; Meegaskumbura, M. Patterns of infection, origins, and transmission of ranaviruses among the ectothermic vertebrates of Asia. Ecol. Evol. 2021, 11, 15498–15519. [Google Scholar] [CrossRef]
  17. Kwon, S.; Park, J.; Choi, W.J.; Koo, K.S.; Lee, J.G.; Park, D. First case of ranavirus-associated mass mortality in a natural population of the Huanren frog (Rana huanrenensis) tadpoles in South Korea. Anim. Cells Syst. 2017, 21, 358–364. [Google Scholar] [CrossRef] [Green Version]
  18. Park, J.; Grajal-Puche, A.; Roh, N.H.; Park, I.K.; Ra, N.Y.; Park, D. First detection of ranavirus in a wild population of Dybowski’s brown frog (Rana dybowskii) in South Korea. J. Ecol. Environ. 2021, 45, 2. [Google Scholar] [CrossRef]
  19. Park, I.K.; Koo, K.S.; Moon, K.Y.; Lee, J.G.; Park, D. PCR detection of ranavirus from dead Kaloula borealis and sick Hyla japonica tadpoles in the wild. Korean J. Herpetol. 2017, 8, 10–14. [Google Scholar]
  20. Mu, W.H.; Geng, Y.; Yu, Z.H.; Wang, K.Y.; Huang, X.L.; Ou, Y.P.; Chen, D.F.; He, C.L.; Zhong, Z.J.; Yang, Z.X.; et al. FV3-like ranavirus infection outbreak in black-spotted pond frogs (Rana nigromaculata) in China. Microb. Pathog. 2018, 123, 111–114. [Google Scholar] [CrossRef]
  21. Box, E.K.; Cleveland, C.A.; Subramaniam, K.; Waltzek, T.B.; Yabsley, M.J. Molecular confirmation of ranavirus infection in amphibians from Chad, Africa. Front. Vet. Sci. 2021, 8, 733939. [Google Scholar] [CrossRef]
  22. Davis, D.R.; Farkas, J.K.; Kruisselbrink, T.R.; Watters, J.L.; Ellsworth, E.D.; Kerby, J.L.; Siler, C.D. Prevalence and distribution of ranavirus in amphibians from southeastern Oklahoma, USA. Herpetol. Conserv. Biol. 2019, 14, 360–369. [Google Scholar]
  23. Galt, N.; Atkinson, M.S.; Glorioso, B.; Waddle, H.; Litton, M.; Savage, A.E. Widespread ranavirus and perkinsea infections in Cuban treefrogs (Osteopilus septentrionalis) invading New Orleans, USA. Herpetol. Conserv. Biol. 2021, 16, 17–29. [Google Scholar]
  24. Xu, K.; Zhu, D.Z.; Wei, Y.; Schloegel, L.M.; Chen, X.F.; Wang, X.L. Broad distribution of ranavirus in free-ranging Rana dybowskii in Heilongjiang, China. EcoHealth 2010, 7, 18–23. [Google Scholar] [CrossRef]
  25. Zhu, Y.Q.; Wang, X.L. Genetic diversity of ranaviruses in amphibians in China: 10 new isolates and their implications. Pak. J. Zool. 2016, 48, 107–114. [Google Scholar]
  26. Jang, H.J.; Suh, J.H. Distribution of amphibian species in South Korea. Korean J. Herpetol. 2010, 2, 45–51. [Google Scholar]
  27. Kang, H.J.; Koo, K.S.; Sung, H.C. Current distribution of American bullfrog Rana catesbeiana Shaw, 1802 in the Republic of Korea. Bioinvasions Rec. 2019, 8, 942–946. [Google Scholar] [CrossRef]
  28. Une, Y.; Sakuma, A.; Matsueda, H.; Nakai, K.; Murakami, M. Ranavirus outbreak in North American bullfrogs (Rana catesbeiana), Japan, 2008. Emerg. Infect. Dis. 2009, 15, 1146. [Google Scholar] [CrossRef]
  29. Saucedo, B.; Serrano, J.M.; Jacinto-Maldonado, M.; Leuven, R.S.; Rocha García, A.A.; Méndez Bernal, A.; Gröne, A.; van Beurden, S.J.; Escobedo-Bonilla, C.M. Pathogen risk analysis for wild amphibian populations following the first report of a ranavirus outbreak in farmed American bullfrogs (Lithobates catesbeianus) from Northern Mexico. Viruses 2019, 11, 26. [Google Scholar] [CrossRef] [Green Version]
  30. Hsieh, C.Y.; Rairat, T.; Chou, C.C. Detection of ranavirus by PCR and in situ hybridization in the American bullfrog (Rana catesbeiana) in Taiwan. Aquaculture 2021, 543, 736955. [Google Scholar] [CrossRef]
  31. ECO Bank. Available online: https://www.nie-ecobank.kr (accessed on 20 February 2020).
  32. NEPARC. Disinfection of Field Equipment to Minimize Risk of Spread of Chytridiomycosis and Ranavirus. Northeast Partners in Amphibian and Reptile Conservation. Available online: https://www.northeastparc.org/products/pdfs/NEPARC-Pub-2014-02_Disinfection_Protocol.pdf (accessed on 10 May 2021).
  33. Galex, I.A.; Gallant, C.M.; D’Avignon, N.; Kuchenbrod, L.M.; Fletcher, C.A.; Rogala, A.R. Evaluation of effective and practical euthanasia methods for larval African clawed frogs (Xenopus laevis). J. Am. Assoc. Lab. Anim. Sci. 2020, 59, 269–274. [Google Scholar] [CrossRef]
  34. Gosner, K.L. A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologica 1960, 16, 183–190. [Google Scholar]
  35. Kimble, S.J.; Karna, A.K.; Johnson, A.J.; Hoverman, J.T.; Williams, R.N. Mosquitoes as a potential vector of ranavirus transmission in terrestrial turtles. EcoHealth 2015, 12, 334–338. [Google Scholar] [CrossRef]
  36. Environmental Geographic Information Service (EGIS). Available online: http://egis.me.go.kr (accessed on 12 December 2021).
  37. Cullen, B.R.; Owens, L.; Whittington, R.J. Experimental infection of Australian anurans (Limnodynastes terraereginae and Litoria latopalmata) with Bohle iridovirus. Dis. Aquat. Org. 1995, 23, 83–92. [Google Scholar] [CrossRef]
  38. Carey, C.; Cohen, N.; Rollins-Smith, L. Amphibian declines: An immunological perspective. Dev. Comp. Immunol. 1999, 23, 459–472. [Google Scholar] [CrossRef]
  39. Ruso, G.E.; Hogan, N.S.; Sheedy, C.; Gallant, M.J.; Jardine, T.D. Effects of agricultural stressors on growth and an immune status indicator in wood frog (Lithobates sylvaticus) tadpoles and metamorphs. Environ. Toxicol. Chem. 2021, 40, 2269–2281. [Google Scholar] [CrossRef]
  40. Brunner, J.L.; Schock, D.M.; Collins, J.P. Transmission dynamics of the amphibian ranavirus Ambystoma tigrinum virus. Dis. Aquat. Org. 2007, 77, 87–95. [Google Scholar] [CrossRef]
  41. Cunningham, A.A.; Hyatt, A.D.; Russell, P.; Bennett, P.M. Emerging epidemic diseases of frogs in Britain are dependent on the source of ranavirus agent and the route of exposure. Epidemiol. Infect. 2007, 135, 1200–1212. [Google Scholar] [CrossRef]
  42. Brenes, R.; Gray, M.J.; Waltzek, T.B.; Wilkes, R.P.; Miller, D.L. Transmission of ranavirus between ectothermic vertebrate hosts. PLoS ONE 2014, 9, e92476. [Google Scholar] [CrossRef]
  43. Gray, M.J.; Miller, D.L.; Schmutzer, A.C.; Baldwin, C.A. Frog virus 3 prevalence in tadpole populations inhabiting cattle-access and non-access wetlands in Tennessee, USA. Dis. Aquat. Org. 2007, 77, 97–103. [Google Scholar] [CrossRef]
  44. Daszak, P.; Cunningham, A.A.; Hyatt, A.D. Anthropogenic environmental change and the emergence of infectious diseases in wildlife. Acta Trop. 2001, 78, 103–116. [Google Scholar] [CrossRef]
  45. Reeve, B.C.; Crespi, E.J.; Whipps, C.M.; Brunner, J.L. Natural stressors and ranavirus susceptibility in larval wood frogs (Rana sylvatica). EcoHealth 2013, 10, 190–200. [Google Scholar] [CrossRef]
  46. Gray, M.J.; Miller, D.L.; Hoverman, J.T. First report of ranavirus infecting lungless salamanders. Herpetol. Rev. 2009, 40, 316. [Google Scholar]
  47. Brunner, J.L.; Olson, A.D.; Rice, J.G.; Meiners, S.E.; Le Sage, M.J.; Cundiff, J.A.; Goldberg, C.S.; Pessier, A.P. Ranavirus infection dynamics and shedding in American bullfrogs: Consequences for spread and detection in trade. Dis. Aquat. Org. 2019, 135, 135–150. [Google Scholar] [CrossRef] [PubMed]
  48. Green, D.E.; Converse, K.A.; Schrader, A.K. Epizootiology of sixty-four amphibian morbidity and mortality events in the USA, 1996–2001. Ann. N. Y. Acad. Sci. 2002, 969, 323–339. [Google Scholar] [CrossRef] [PubMed]
  49. Warne, R.W.; LaBumbard, B.; LaGrange, S.; Vredenburg, V.T.; Catenazzi, A. Co-infection by chytrid fungus and ranaviruses in wild and harvested frogs in the tropical Andes. PLoS ONE 2016, 11, e0145864. [Google Scholar] [CrossRef]
  50. Green, D.E.; Converse, K.A. Diseases of frogs and toads. In Wildlife Diseases: Landscape Epidemiology, Spatial Distribution and Utilization of Remote Sensing Technology; Majumdar, S.K., Huffman, J.E., Brenner, F.J., Panah, A.I., Eds.; Pennsylvania Academy of Science: Easton, PA, USA, 2005; pp. 89–117. ISBN 978-0-9458-0919-7. [Google Scholar]
  51. Greer, A.L.; Berrill, M.; Wilson, P.J. Five amphibian mortality events associated with ranavirus infection in south central Ontario, Canada. Dis. Aquat. Org. 2005, 67, 9–14. [Google Scholar] [CrossRef] [Green Version]
  52. Haislip, N.A.; Gray, M.J.; Hoverman, J.T.; Miller, D.L. Development and disease: How susceptibility to an emerging pathogen changes through anuran development. PLoS ONE 2011, 6, e22307. [Google Scholar] [CrossRef]
  53. Borzée, A.; Heo, K.; Jang, Y. Relationship between agro-environmental variables and breeding Hylids in rice paddies. Sci. Rep. 2018, 8, 8049. [Google Scholar] [CrossRef]
  54. Park, D.; (Division of Science Education, Kangwon National University, Kangwon, Korea). Personal communication, 2022.
Viruses 14 01073 g001 550
Figure 1. Local ranavirus infection in larvae of three anuran species (Dryophytes japonicus, green circles; Pelophylax nigromaculatus, blue circles; Lithobates catesbeianus, red circles) at 49 locations across eight provinces in South Korea. Filled circles indicate the confirmation of at least one ranavirus-infected tadpole of the species at the location.
Figure 1. Local ranavirus infection in larvae of three anuran species (Dryophytes japonicus, green circles; Pelophylax nigromaculatus, blue circles; Lithobates catesbeianus, red circles) at 49 locations across eight provinces in South Korea. Filled circles indicate the confirmation of at least one ranavirus-infected tadpole of the species at the location.
Viruses 14 01073 g001
Viruses 14 01073 g002 550
Figure 2. Box plots (medians ± quartiles) of the ranavirus infection rates among larvae of three anuran species (Dryophytes japonicus, Pelophylax nigromaculatus, Lithobates catesbeianus). The different characters on the bars between species indicate significantly different rates using the Kruskal–Wallis test with the Dunn-Bonferroni post hoc test (p = 0.001).
Figure 2. Box plots (medians ± quartiles) of the ranavirus infection rates among larvae of three anuran species (Dryophytes japonicus, Pelophylax nigromaculatus, Lithobates catesbeianus). The different characters on the bars between species indicate significantly different rates using the Kruskal–Wallis test with the Dunn-Bonferroni post hoc test (p = 0.001).
Viruses 14 01073 g002
Viruses 14 01073 g003 550
Figure 3. Changes in the ranavirus infection rate according to the Gosner developmental stage, based on combined data from larvae of three anuran species (Dryophytes japonicus, Pelophylax nigromaculatus, Lithobates catesbeianus).
Figure 3. Changes in the ranavirus infection rate according to the Gosner developmental stage, based on combined data from larvae of three anuran species (Dryophytes japonicus, Pelophylax nigromaculatus, Lithobates catesbeianus).
Viruses 14 01073 g003
Table
Table 1. Sampling locations and ranavirus infection rates in larvae of three anuran species (Dryophytes japonicus, Pelophylax nigromaculatus, Lithobates catesbeianus) across 49 sampling locations in South Korea.
Table 1. Sampling locations and ranavirus infection rates in larvae of three anuran species (Dryophytes japonicus, Pelophylax nigromaculatus, Lithobates catesbeianus) across 49 sampling locations in South Korea.
ProvinceSampling
Location
(Abbreviation)		Infection Rate (%) (No. of Infected/Tested Tadpoles)
D. japonicusP. nigromaculatusL. catesbeianus
GangwonSokcho (SC)NoneNone-
Gangneung (GaN)12.5 (1/8)None-
Inje (InJ)12.5 (1/8)None-
Yanggu (YG)None12.5 (1/8)-
Chuncheon1 (CC1)12.5 (1/8)None-
Chuncheon2 (CC2)NoneNone-
Hongcheon1 (HC1)NoneNone-
Hongcheon2 (HC2)NoneNone-
Hongcheon3 (HC3)25.0 (2/8)75.0 (6/8)-
GyeonggiYeoju (YJ)NoneNone-
Goyang (GY)None--
ChungbukCheongju1 (CJ1)12.5 (1/8)None-
Cheongju2 (CJ2)None--
Okcheon (OkC)NoneNone-
ChungnamYesan (YS)--None
Seosan (SS)12.5 (1/8)NoneNone
Taean (TA)37.5 (3/8)NoneNone
Sejong (SeJ)NoneNone50.0 (1/2)
Cheongyang (CY)NoneNone25.0 (2/8)
Seocheon (SuCh)None25.0 (2/8)None
GyeongbukMungyeong (MK)NoneNoneNone
Sangju (SJ)NoneNone-
Uisung (UiS)NoneNone-
Gimcheon (KiC)NoneNone100.0 (1/1)
Seongju (SuJ)None12.5 (1/8)12.5 (1/8)
Daegu (DG)12.5 (1/8)12.5 (1/8)None
Pohang (PH)None37.5 (3/8)-
GyeongnamUlsan (UlS)NoneNone-
Busan1 (BS1)--None
Busan2 (BS2)12.5 (1/8)None-
Sacheon (SaC)None37.5 (3/8)None
Jinju (JJ)NoneNoneNone
JeonbukIksan (IS)12.5 (1/8)NoneNone
Gimjae (KJ)NoneNone75.0 (6/8)
Jeongup (JuU)12.5 (1/8)NoneNone
Imsil (ImS)None12.5 (1/8)-
Sunchang (SuCa)NoneNone-
JeonnamJangseong (JS)12.5 (1/8)NoneNone
Gokseong (GoS)NoneNone-
Muan (MA)12.5 (1/8)25.0 (2/8)None
Imja (IJ)-None50.0 (4/8)
Sinan1 (SAB)NoneNone-
Sinan2 (SABC)None12.5 (1/8)None
Sinan3 (SAGI)--None
Sinan4 (SAJD)-None37.5 (3/8)
Sinan5 (SAJU)--None
Yeongam (YA)None25.0 (2/8)None
Haenam (HN)25.0 (2/8)NoneNone
Goheung (GoH)NoneNone50.0 (4/8)
Average
(No. of infected/tested tadpoles)		49 locations16.1%
(18/344)		26.1%
(23/342)		50.0%
(22/187)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Roh, N.; Park, J.; Kim, J.; Kwon, H.; Park, D. Prevalence of Ranavirus Infection in Three Anuran Species across South Korea. Viruses 2022, 14, 1073. https://0-doi-org.brum.beds.ac.uk/10.3390/v14051073

AMA Style

Roh N, Park J, Kim J, Kwon H, Park D. Prevalence of Ranavirus Infection in Three Anuran Species across South Korea. Viruses. 2022; 14(5):1073. https://0-doi-org.brum.beds.ac.uk/10.3390/v14051073

Chicago/Turabian Style

Roh, Namho, Jaejin Park, Jongsun Kim, Hyerim Kwon, and Daesik Park. 2022. "Prevalence of Ranavirus Infection in Three Anuran Species across South Korea" Viruses 14, no. 5: 1073. https://0-doi-org.brum.beds.ac.uk/10.3390/v14051073

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop