755
Views
0
CrossRef citations to date
0
Altmetric
Rapid Communications

DNA barcoding indicates the range extension in an endemic frog Nyctibatrachus jog, from the Western Ghats, India

ORCID Icon, & ORCID Icon
Pages 2468-2474 | Received 15 Apr 2021, Accepted 12 Jul 2021, Published online: 26 Jul 2021

Abstract

The frogs of genus Nyctibatrachus from the Western Ghats are endemic, with some taxa showing a narrow distribution range. Nyctibatrachus jog was known only from the type locality, Jog falls from Sharavathi river basin suggesting a restricted distribution. In this study, using DNA barcoding, we studied the distribution patterns of N. jog by sampling multiple river basins. 16S rRNA and Cytochrome b genes were used to distinguish N. jog from its congeners as well as to infer intra-species relationships. The results from the 16S rRNA gene showed 99% similarity of the collected individuals with the type specimen from Jog Falls confirming the identity of N. jog. The results indicate that N. jog has wide distribution extending its range in multiple river basins in the Western Ghats, India. This study also provides the Area of Occurrence and Extent of Occurrence of N. jog which could help in developing strategies for its conservation.

Introduction

DNA barcoding has been routinely used to identify species, contributing to biodiversity research (Hebert et al. Citation2003; Hajibabaei et al. Citation2007). It has also been useful to understand the geographical distributions of populations, providing insights into the patterns of geneflow within a species (Hajibabaei et al. Citation2007). Using DNA barcoding, many studies have identified the new extent of the species thus expanding their distribution range (Seshadri et al. Citation2016; Kundu et al. Citation2020; Lim et al. Citation2020).

The Western Ghats of India is known for the high diversity of amphibians with more than 90% of the species endemic to the region. There are nine families of anurans found here, of which three families, Micrixalidae Dubois et al. (Citation2001), Ranixalidae Dubois (Citation1987) and Nasikabatrachidae Biju and Bossuyt (Citation2003) are endemic to the Western Ghats while family Nyctibatrachidae Blommers-Schlösser (Citation1993) has genera Nyctibatrachus and Astrobatrachus distributed in the Western Ghats and genus Lankanectus endemic to Sri Lanka. Some of these taxa have narrow distribution ranges, making them of high conservation concern (Vijayakumar et al. Citation2019).

The genus Nyctibatrachus is endemic to the Western Ghats, India. There are 36 species at present (Biju et al. Citation2011; Gururaja et al. Citation2014; Garg et al. Citation2017; Krutha et al. Citation2017). The characteristics of Nyctibatrachus frogs include wrinkled skin, rhomboid pupils, pointed vomerine teeth, subocular glands and a notched tongue (Biju et al. Citation2011). The larger Nyctibatrachus species are associated with torrential streams whereas the smaller Nyctibatrachus species are associated with marshy pools. (Biju et al. Citation2011). A study by Van Bocxlaer et al. (Citation2012) revealed that the species show mountain-associated clade level endemism with a narrow distribution range. The authors mentioned that torrentially adapted Nyctibatrachus species disperse along the river systems but may have restricted distribution range due to the east-west orientation of the river system, thus limiting their expansion from one river system to another. However, fine-scale studies understanding the distribution range for many Nyctibatrachus species is lacking. Knowledge of the distribution extent can be helpful for their conservation.

Nyctibatrachus jog Biju et al. (Citation2011) is an endemic frog from central Western Ghats India. The species was discovered from Jog falls, Karnataka at an elevation of 600 m by Biju et al. (Citation2011). The species is adapted to fast-flowing streams and torrents (Biju et al. Citation2011). The authors mentioned that N. jog was found only at the type locality, Jog Falls, central Western Ghats, India. However, the exact distribution range of N. jog is not known. In this study, we provide additional distribution records of N. jog, identified using DNA barcoding methods. We also provide the Extent of Occurrence (EOO) and Area of Occurrence (AOO) for N. jog.

Material and methods

As a part of a study on the ecology and genetics of Nyctibatrachus species, we sampled the streams from Bedti, Aghanashini and Sharavathi () river basins from the central Western Ghats of India. These rivers originate in the Western Ghats and flow westwards to join the Arabian sea. The river basins are differentiated based on their catchment areas and delineated using Digital Elevation models (DEM). 43 streams were sampled from the three river basins. The streams were sampled for both adults and tadpoles. The tadpoles were sampled using dipnets (15 cm by 15 cm). To distinguish N. jog from its congeners, individuals were identified by their calls, size, sex, clutch size following Biju et al. (Citation2011), Gururaja et al. (Citation2014) and by DNA barcoding. A small tadpole tail tip was excised and preserved in 80% ethanol for molecular work and the tadpoles were released back into the streams.

Figure 1. Current distribution records of Nyctibatrachus jog along the Western Ghats. The type locality of Jog Falls is given in rhombus. Triangle indicates the type locality as per Biju et al. (Citation2011) latitude and longitude value, which is 150 km south of the actual type locality. Jog Falls. Yellow circles represent sites sampled in this study. Blue circle represents the presence of N. jog not sampled in this study.

Figure 1. Current distribution records of Nyctibatrachus jog along the Western Ghats. The type locality of Jog Falls is given in rhombus. Triangle indicates the type locality as per Biju et al. (Citation2011) latitude and longitude value, which is 150 km south of the actual type locality. Jog Falls. Yellow circles represent sites sampled in this study. Blue circle represents the presence of N. jog not sampled in this study.

Molecular methods

The DNA was extracted using the salt extraction method as mentioned in Vences et al. (Citation2012). A 530 bp of 16S rRNA gene (Palumbi et al. Citation2002) was used for identification as it has been widely used for DNA barcoding of amphibians (Vences et al. Citation2005). We also used 427 bp Cytochrome b (Cyt-b) gene (mcb 398: TACCATGAGGACAAATATCATTCTG3 mcb 869: CCTCCTAGTTTGTTAGGGATTGATCG3, Verma and Singh Citation2002). 23 individuals were sequenced for 16S rRNA and 26 individuals were sequenced for Cyt-b. The sequences were checked for identity using nucleotide BLAST (https://blast.ncbi.nlm.nih.gov/) and deposited in GenBank (16S Accession numbers MW081889–MW081909, Cyt-b Accession numbers MW762640–MW762665). Two 16S sequences viz adult N. jog (Voucher number BNHS5457) with accession number JN644900 (Van Bocxlaer et al. Citation2012) and tadpole N. jog (Voucher number BNHS5900) with accession number KP317819 (Priti et al. Citation2015) were downloaded from GenBank to confirm the identity of genetic sequences generated in this study. The monophyly of N. jog individuals were analyzed by building phylogenetic tree for 16S rRNA gene using Bayesian inference in Mr Bayes (Ronquist and Huelsenbeck Citation2003) and Maximum Likelihood methods in RAxMLGUI (Silvestro and Michalak Citation2012). GTR + I + G was used as substitution model as selected by Partition Finder v 1.1.1 (Lanfear et al. Citation2012). Lankanectes corrugatus (Accession number AF215393) was used as an outgroup. The genetic distance was calculated using Kimura 2 parameter method in MEGA 7 (Kumar et al. Citation2016). To understand the distribution range of N. jog, the Area of Occurrence (AOO) and the Extent of Occurrence (EEO) were estimated using GeoCat (Bachman et al. Citation2011).

Results

The new localities where N. jog individuals were observed are summarized in . The phylogenetic tree showed that the N. jog individuals formed a strong monophyly (Bootstrap support = 94, Posterior Probability = 1, ). The generated sequences showed 99.2% similarity with the published 16S rRNA sequence from the type locality, Jog falls mentioned in Biju et al. (Citation2011), Van Bocxlaer et al. (Citation2012) respectively (Accession number JN644900). The individuals collected from multiple localities of Sharavathi river basin clustered with the adult N. jog (Accession number JN644900) as well as with the tadpole N. jog (Accession number KP317819) confirming the identity of N. jog (). Although the genetic data of Biju et al. (Citation2011) matches with our study, the location mentioned as Jog falls (13.4421 N, 75.1704E) by Biju et al. (Citation2011) falls in Tunga river basin which is 150 km away from Jog falls. This could be an error in coordinates provided by the authors. The actual co-ordinates of Jog Falls are 14.2294°N, 74.74679°E. From our study, the occurrence of the N. jog now extends into multiple localities in Aghanashini, Bedti and Sharavathi river basins. The maximum interspecies genetic distance was with N. periyar (13%) while minimum interspecific distance was with N. petraeus (1%). Although the interspecies distance between N. jog and N. petraeus is less for 16S r RNA gene, the genetic distance between the two species for ND1 gene was 4.5% along with the difference in morphology, based on which, they were distinguished as distinct species (Biju et al. Citation2011). The intraspecific genetic distance among N. jog individuals for 16S rRNA gene ranged from 0.00 to 0.8% () while for the Cyt-b gene, the intraspecific genetic distance ranged between 0.00 to 1.95% (). The intraspecies genetic relationship as shown by Maximum likelihood and Bayesian methods showed populations of Sharavathi river basin as a separate clade (BS > 95, PP = 1, ) from the populations of Aghanashini and Bedti river basin whereas there was some admixture between populations of Aghanashini and Bedti river basins. Further studies understanding intra-species genetic diversity and gene flow in N. jog is needed. The area of occupancy (AOO) for N. jog is 32 km2 and the extent of occupancy (EOO) is 461.420 km2. Besides these sites, personal observations by naturalist Manu Nakathaya from Reserve Forest in Mala (Suvarna river basin) near Kudremukh also suggests the presence of N. jog based on morphology and call records, which is further south by 110 km from the type locality. If these points are included, the area of occupancy (AOO) for N. jog increases to 36 km2 and the extent of occupancy (EOO) increases to 1629.464 km2.

Figure 2. Maximum Likelihood phylogenetic tree depicting the inferred relationships among members of the genus Nyctibatrachus based on 530 base pairs of the mitochondrial 16S gene. The number above the branches denotes bootstrap support. The numbers below the branches denote Posterior Probability values. Blue vertical line represents Sharavathi river basin, Yellow Line represents Aghanashini river basin and Red line indicates Bedti river basin.

Figure 2. Maximum Likelihood phylogenetic tree depicting the inferred relationships among members of the genus Nyctibatrachus based on 530 base pairs of the mitochondrial 16S gene. The number above the branches denotes bootstrap support. The numbers below the branches denote Posterior Probability values. Blue vertical line represents Sharavathi river basin, Yellow Line represents Aghanashini river basin and Red line indicates Bedti river basin.

Table 1. Sampling localities of N. jog along with GenBank accession numbers.

Table 2. Kimura2 parameter pairwise genetic distance (%) among 36 Nyctibatrachus species from the Western Ghats based on 530 bp of mitochondrial 16S rRNA gene (1–23: GenBank Accession numbers of N.jog).

Table 3. Kimura2 parameter pairwise genetic distance (%) between individuals of Nyctibatrachus jog from the Western Ghats based on 427 bp of mitochondrial Cyt-b gene.

Discussion

The range extension of N. jog suggests either lack of fine-scale sampling by previous studies or the difficulty in sighting this species due to its specific microhabitat preferences. N. jog prefers fast flowing streams and waterfalls which may not be easily accessible for sampling leading to few observations. Even in this study, the sightings were only in 9 streams out of the 43 streams sampled. Seasonality could be another reason as N. jog calls are not heard post monsoon making it difficult to identify in the field. Also, Nyctibatrachus genus harbors cryptic species, hence previous studies based on morphological observations may have misidentified N. jog for another Nyctibatrachus species as it was earlier known to have restriction distribution. This DNA barcoding study now indicates the wide distribution of N. jog as against the previous notion that it is distributed only in the type locality of the Sharavathi river basin (Biju et al. Citation2011).

The habitats of N. jog are perennial fast-flowing streams and torrents. The major threats to its habitats are the conversion of the streams into areca plantations and paddy fields. Since N. jog inhabits torrents, there is a threat to their habitat by tourists who use the torrents for bathing and defecating. According to the International Union of Conservation of Nature (IUCN) Global Redlist assessment, the status of N. jog is not evaluated (Dandekar et al. Citation2020). From this study, it is evident that the EOO of N. jog has increased, thus it is no longer restricted to its type locality, Jog Falls of the Sharavathi river basin. Since this species appears to have wide distribution, efforts are needed for extensive fieldwork locating new populations. To increase the conservation efforts, multiple stakeholders must be involved that includes educating local people and forest department officials. Besides that, it is also essential to study population dynamics and population genetic studies for developing conservation strategies.

Acknowledgements

We would like to thank the Karnataka Forest Department for permission (No.PCCF (WL)/E2/CR-49/2018-19). Priti H and Gururaja KV are thankful to Ashok Hegde, Mavingundi; Suhas Hegde, Drongo Homestay Sirsi and Vivek Divekar, Satwik Homestay Sirsi for the logistical support. Priti H and Anisha A are thankful to Kartik Sunagar and Evolutionary Venomics Lab for providing the laboratory space.

Disclosure statement

The authors report no conflict of interest

Data availability statement

The data that support the findings of this study are openly available in [NCBI] at https://www.ncbi.nlm.nih.gov/, reference numbers MW081889 to MW081909, MW762640 to MW762665. The voucher specimens mentioned in the study are deposited in Bombay Natural History Society (BNHS), Mumbai, email: [email protected].

Additional information

Funding

This work was supported by Department of Science and Technology (DST), India to Priti Hebbar under the INSPIRE Faculty Fellowship [DST/INSPIRE/04/2017/003152].

References

  • Bachman S, Moat J, Hill AW, De La Torre J, Scott B. 2011. Supporting Red List threat assessments with GeoCAT: geospatial conservation assessment tool. Zookeys. 150:117–126.
  • Biju SD, Bossuyt F. 2003. New frog family from India reveals an ancient biogeographical link with the Seychelles. Nature. 425(6959):711–714.
  • Biju SD, Van Bocxlaer I, Mahony S, Dinesh KP, Radhakrishnan C, Zachariah A, Giri V, Bossuyt F. 2011. A taxonomic review of the Night Frog genus Nyctibatrachus Boulenger, 1882 in the Western Ghats, India (Anura: Nyctibatrachidae) with description of twelve new species. Zootaxa. 3029(1):1–96.
  • Blommers-Schlösser RM. 1993. Systematic relationships of the Mantellinae Laurenst 1946 (Anura Ranoidea). Ethol Ecol Evol. 5(2):199–218.
  • Dubois A, Ohler A, Biju SD. 2001. A new genus and species of Ranidae (Amphibia, Anura) from south-western India. Alytes. 19(2–4):53–79.
  • Dubois A. 1987. Miscellanea taxinomica batrachologica (I). Alytes. 5:7–95.
  • Garg S, Suyesh R, Sukesan S, Biju SD. 2017. Seven new species of Night Frogs (Anura, Nyctibatrachidae) from the Western Ghats Biodiversity Hotspot of India, with remarkably high diversity of diminutive forms. Peer J. 5:e3007.
  • Gururaja KV, Dinesh KP, Priti H, Ravikanth G. 2014. Mud-packing frog: a novel breeding behavior and parental care in a stream dwelling new species of Nyctibatrachus (Amphibia, Anura, Nyctibatrachidae). Zootaxa. 3796(1):33–61.
  • Hajibabaei M, Singer GA, Hebert PD, Hickey DA. 2007. DNA barcoding: how it complements taxonomy, molecular phylogenetics and population genetics. Trends Genet. 23(4):167–172.
  • Hebert PD, Cywinska A, Ball SL, Dewaard JR. 2003. Biological identifications through DNA barcodes. Proc R Soc Lond B. 270(1512):313–321.
  • Krutha K, Dahanukar N, Molur S. 2017. Nyctibatrachus mewasinghi, a new species of night frog (Amphibia: Nyctibatrachidae) from Western Ghats of Kerala. J Threat Taxa. 9(12):10985–10997.
  • Kumar S, Stecher G, Tamura K. 2016. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33(7):870–1874.
  • Kundu S, Lalremsanga HT, Purkayastha J, Biakzuala L, Chandra K, Kumar V. 2020. DNA barcoding elucidates the new altitude record and range-extension of lesser-known bullfrog (Hoplobatrachus litoralis) in northeast India. Mitochondrial DNA Part B. 5(3):2668–2672.
  • Lanfear R, Calcott B, Ho SYW, Guindon S. 2012. Partitionfinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 29(6):1695–1701.
  • Lim BK, Loureiro LO, Garbino GS. 2020. Cryptic diversity and range extension in the big-eyed bat genus Chiroderma (Chiroptera, Phyllostomidae). Zookeys. 918:41–63.
  • Dandekar N, Gururaja KV, Priti H, Hegde A, Komanduri K, Modak N, Padhye A, Patel H, Sreedharan G. 2020. Conservation needs assessment for Nyctibatrachus jog, India. AArk/ASG India Assessment Workshop 2020; [accessed 2021 Mar 23]. https://www.conservationneeds.org/Assessment/AssessmentResults?assessmentId=5606&countryId=146&speciesId=6808.
  • Palumbi S, Martin A, Romano S, McMillan WO, Stice L, Grabowski G. 2002. The simple fool’s guide to PCR. Version 2.0. Department of Zoology, Honolulu, 45. http://palumbi.stanford.edu/SimpleFoolsMaster.pdf.
  • Priti H, Gururaja KV, Ravikanth G. 2015. Morphology, natural history and molecular identification of tadpoles of three endemic frog species of Nyctibatrachus Boulenger, 1882 (Anura: Nyctibatrachidae) from Central Western Ghats, India. J Nat Hist. 49(43–44):2667–2681.
  • Ronquist F, Huelsenbeck JP. 2003. MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 19(12):1572–1574.
  • Seshadri KS, Priti H, Ravikanth G, Vidisha MK, Vineeth KK, Singal R, Sarma RR, Aravind NA, Gururaja KV. 2016. Redescription and range extension of Microhyla sholigari Dutta & Ray (Amphibia: Anura: Microhylidae) from South West India. Zootaxa. 4208(6):547.
  • Silvestro D, Michalak I. 2012. RAxMLGUI: a graphical front-end for RAxML. Org Divers Evol. 12(4):335–337.
  • Van Bocxlaer I, Biju SD, Willaert B, Giri VB, Shouche YS, Bossuyt F. 2012. Mountain-associated clade endemism in an ancient frog family (Nyctibatrachidae) on the Indian subcontinent. Mol Phylogenet Evol. 62(3):839–847.
  • Vences M, Thomas M, Van der Meijden A, Chiari Y, Vieites DR. 2005. Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians. Front Zool. 2(1):5–12.
  • Vences M, Nagy ZT, Sonet G, Verheyen E. 2012. DNA barcoding amphibians and reptiles. Methods Mol Biol. 858:79–107.
  • Verma SK, Singh L. 2002. Novel universal primers establish the identity of an enormous number of animal species for forensic application. Mol. Ecol. Notes. 3(1):28–31.
  • Vijayakumar SP, Pyron RA, Dinesh KP, Torsekar VR, Srikanthan AN, Swamy P, Stanley EL, Blackburn DC, Shanker K. 2019. A new ancient lineage of frog (Anura: Nyctibatrachidae: Astrobatrachinae subfam. nov.) endemic to the Western Ghats of Peninsular India. Peer J. 7:e6457.