Research Article |
Corresponding author: Rajeev Raghavan ( rajeevraq@hotmail.com ) Academic editor: Ralf Britz
© 2023 Neelesh Dahanukar, Remya L. Sundar, Duwaki Rangad, Graham Proudlove, Rajeev Raghavan.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Dahanukar N, Sundar RL, Rangad D, Proudlove G, Raghavan R (2023) The world’s largest cave fish from Meghalaya, Northeast India, is a new species, Neolissochilus pnar (Cyprinidae, Torinae). Vertebrate Zoology 73: 141-152. https://doi.org/10.3897/vz.73.e101011
|
The world’s largest subterranean fish was discovered in 2019, and was tentatively identified as a troglomorphic form of the golden mahseer, Tor putitora. Detailed analyses of its morphometric and meristic data, and results from molecular analyses now reveal that it is a new species of the genus Neolissochilus, the sister taxon of Tor. We formally describe the new species as Neolissochilus pnar, honouring the tribal communities of East Jaintia hills in Meghalaya, Northeast India, from where it was discovered. Neolissochilus pnar possesses a number of characters unique among species of Neolissochilus, with the exception of the similarly subterranean N. subterraneus from Thailand. The unique characters that diagnose N. pnar from all epigean congeners comprise highly reduced eye size to complete absence of externally visible eyes, complete lack of pigmentation, long maxillary barbels, long pectoral-fin rays, and scalation pattern. Neolissochilus pnar is distinguished from the hypogean N. subterraneus, the type locality of which is a limestone cave ~2000 kms away in Central Thailand, by a lesser pre-pelvic length (47.8–49.4 vs. 50.5–55.3 %SL), a shorter caudal peduncle (16.1–16.8 vs. 17.8–23.7 %SL), and shorter dorsal fin (17.4–20.8 vs. 21.5–26.3 %SL). In addition, Neolissochilus pnar is also genetically and morphologically distinct from its close congeners with a raw genetic divergence of 1.1–2.7% in the COI gene with putative topotype of N. hexastichus and 2.1–2.6% with putative topotype of N. hexagonolepis.
Eastern Himalaya, limestone cave, mahseer, new species, subterranean fishes
Roughly 1.6% (293 species) of all known (~18,000) freshwater fish species live their whole lives either in caves, or in groundwater aquifers (
Most subterranean fishes have evolved a small-sized body plan to meet the limitations in food resources, light availability and space in underground habitats. The mean size of subterranean fish species is 85.5 mm, with most species below 130 mm (
List of subterranean fish species with adult sizes in excess of 200 mm standard length (SL).
Family/Species | Country | Maximum SL (mm) |
Cyprinidae | ||
Sinocyclocheilus guanyangensis | China | 202 |
Sinocyclocheilus hugeibarbus | China | 217 |
Neolissochilus subterraneus | Thailand | 217 |
Heptapteridae | ||
Rhamdia enfurnada | Brazil | 218 |
Synbranchidae | ||
Typhlosynbranchus luticolus | Cameroon | 209 |
Rakthamichthys digressus | India | 242 |
Ophisternon infernale | Mexico | 325 |
Ophisternon candidum | Australia | 356 |
Stories of a ‘white cavefish’ from the Siju Caves in the Garo Hills of Meghalaya, Northeast India have been documented for 100 years, but were suggested to be slightly decolorized specimens of Neolissochilus hexastichus (M’Clelland) that appeared almost white when observed inside the water, under the light of a torch (
The availability of additional fresh specimens of this unique cyprinid fish has now enabled us to study its morphological characters in more detail and to include it in a molecular genetic analysis. This combined evidence reveals that the world’s largest cavefish is an undescribed species of the cyprinid genus Neolissochilus, for which we make a name available below.
The specimens of our study were collected from the Krem Um Ladaw and the Krem Chympe caves in Meghalaya, Northeast India in 2019 and 2020. All the specimens were fixed in 4% formaldehyde, after preserving pectoral fin-clips in absolute ethanol for DNA analysis. Specimens are deposited in the museum collection of the Kerala University of Fisheries and Ocean Studies (KUFOS), Kochi, India.
Characterization and analysis of morphometric and meristic information was carried out in line with previous studies on members of the subfamily Torinae (
DNA was extracted from alcohol preserved fin-clips of the three specimens in the type series using QIAamp® DNA Mini Kit (Qiagen, Germany) following manufacturer’s protocol. Three mitochondrial genes, i.e., cytochrome oxidase subunit 1 (COI), cytochrome b (cytb) and large subunit ribosomal ribonucleic acid (16S) were amplified, purified and sequenced following published protocols (
A total of nine sequences were generated for the three genes (COI, cytb and 16S) from the holotype and two paratypes. GenBank accession numbers and GenSeq nomenclature (Chakrabarty et al. 2013) for sequences generated in the current study are provided in Table
GenBank accession numbers and GenSeq nomenclature for sequences generated in the current study.
Species | Locality | Voucher | COI | cytb | 16S | GenSeq |
Neolissochilus pnar | Krem Um Ladaw | KUFOS.F.2022.701 | OQ351360 | OQ349705 | OQ357607 | genseq-1 COI, cytb, 16S |
Neolissochilus pnar | Krem Um Ladaw | KUFOS.F.2022.702 | OQ351361 | OQ349706 | OQ357608 | genseq-2 COI, cytb, 16S |
Neolissochilus pnar | Krem Chympe | KUFOS.F.2022.703 | OQ351362 | OQ349707 | OQ357609 | genseq-2 COI, cytb, 16S |
Because COI sequences were available for a larger dataset of Neolissochilus, we performed a separate ML analysis (as described above) including all available COI sequences of Neolissochilus in GenBank (Table S3). The best partition scheme and nucleotide substitution analysis for COI dataset is provided in Table S4. We performed molecular species delimitation using Assemble Species by Automatic Partitioning (ASAP), employing uncorrected genetic distances, for barcode gap analysis and species delimitation (
KUFOS.F.2022.701, 329.2 mm SL, 92 m below the surface in Krem Um Ladaw, Meghalaya, India; collected 7 Jan 2020.
KUFOS.F.2022.702, 179.7 mm, same locality as holotype, collected 21 Feb 2019; KUFOS.F.2022.703, 208.9 mm SL, Krem Chympe cave, Meghalaya, India, collected 7 Jan 2020.
The species name pnar, honours the ‘pnar’, the sub-tribal group of the Khasi people in the state of Meghalaya, India.
Neolissochilus pnar is distinguished from all its congeners by mandibular barbel long, reaching anterior margin of opercle (vs. short, not reaching margin of opercle). It is further distinguished from all epigean congeners by atrophied eyes, highly reduced in size in juveniles and small-adults and absence of externally visible eyes in adults (vs. presence of well-developed eyes in all life-stages); complete absence of pigmentation (vs. presence); long pectoral-fin reaching anterior base of pelvic fin (vs. short, not reaching anterior base of pelvic fin); and distinct scalation pattern with 28+2 (2) or 31+1 (1) lateral line scales, 8 scales in transverse series with 4 above the lateral line and 3 below the lateral line. Neolissochilus pnar is distinguished from the only other subterranean congener, N. subterraneus by shorter pre-pelvic length (47.8–49.4 vs. 50.5–55.3 %SL), shorter caudal-peduncle length (16.1–16.8 vs. 17.8–23.7 %SL) and shorter dorsal-fin length (17.4–20.8 vs. 21.5–26.3 %SL).
General appearance as in Fig.
Morphometric and meristic data of the holotype and two paratypes of Neolissochilus pnar.
Holotype | Paratype | Paratype | |
KUFOS.F.2022.701 | KUFOS.F.2022.702 | KUFOS.F.2022.703 | |
Morphometric information | |||
Total length (TL, mm) | 409.9 | 216.1 | 255.2 |
Standard length (SL, mm) | 329.2 | 179.7 | 208.9 |
Head length (HL, mm) | 98.9 | 47.8 | 60.8 |
% SL | |||
Head Length | 30.0 | 26.6 | 29.1 |
Snout Length | 12.6 | 10.6 | 11.1 |
Pre-dorsal length | 49.8 | 47.4 | 48.6 |
Pre-pectoral length | 30.8 | 26.6 | 29.3 |
Pre-pelvic length | 49.4 | 47.8 | 48.9 |
Pre-anal length | 73.5 | 71.2 | 75.1 |
Caudal peduncle length | 16.6 | 16.1 | 16.8 |
Dorsal-fin length | 17.4 | 20.8 | 19.4 |
Dorsal-fin base length | 16.1 | 15.8 | 16.4 |
Pectoral-fin length | 22.9 | 22.2 | 23.4 |
Pectoral-fin base length | 5.3 | 4.6 | 5.8 |
Pelvic-fin length | 19.8 | 18.0 | 16.4 |
Pelvic-fin base length | 4.7 | 5.3 | 5.0 |
Anal-fin length | 19.5 | 15.8 | 16.5 |
Anal-fin base length | 6.5 | 8.5 | 7.7 |
Caudal-fin length | 24.7 | 25.9 | 23.6 |
Caudal-fin base length | 10.7 | 12.5 | 11.6 |
Body depth at dorsal fin | 25.1 | 24.3 | 23.9 |
Body depth at anal fin | 13.9 | 15.7 | 16.4 |
Body width at dorsal fin | 12.6 | 13.5 | 14.3 |
Body width at anal fin | 7.3 | 8.3 | 7.2 |
Caudal-peduncle depth | 9.0 | 10.1 | 9.6 |
% HL | |||
Snout length | 42.0 | 39.9 | 38.3 |
Inter-orbital length | 38.5 | 36.8 | 37.7 |
Maxillary barbel length | 45.9 | 43.6 | 40.1 |
Rostral barbel length | 44.5 | 44.7 | 43.6 |
Meristic information | |||
Dorsal-fin rays | iv, 9 | iv, 9 | iv, 9 |
Pectoral-fin rays | i, 15 | i, 15 | i, 15 |
Pelvic-fin rays | i, 8 | i, 8 | i, 8 |
Anal-fin rays | iii, 5 | iii, 5 | iii, 5 |
Principal caudal-fin rays | 19 | 19 | 19 |
Perforated lateral-line scales | 31+1 | 28+2 | 28+2 |
Dorsal fin with 13 soft rays (iv+9), its origin almost midway between tip of snout and end of caudal peduncle, or slightly in advance. Posterior margin of adpressed dorsal fin reaching anal-fin origin. Pectoral fin with 16 rays (i+15), its length shorter than head length. Adpressed pectoral fin reaching vertical at dorsal-fin origin, and almost reaching pelvic-fin origin. Pelvic-fin with 9 rays (i+8), its origin slightly posterior to vertical at dorsal-fin origin. Anal fin with 8 rays (iii+5). Caudal fin forked with 19 principal caudal rays. Caudal peduncle 2–2.3 times as long as deep.
Body lateral line continuous, with 28–31 perforated scales, and an additional 1–2 on caudal-fin base. Transverse series with 8 scale rows, 4 scale rows between dorsal-fin origin and row of lateral line scales, 3 scale rows between row of lateral line scales and pelvic-fin origin. Pre-dorsal scales 9.
In life (Fig.
The species is known from the caves at Krem Um Ladaw, and the adjacent Krem Chympe in Jaintia Hills, Meghalaya, India, which drain into the Meghna River System (Fig.
The entrance to the cave in Krem Um Ladaw is in the form of a large open pitch head, lies in a large, rocky, seasonally dry streambed within a forest. The entrance series is predominantly vertical with some short (<20 m) horizontal to steeply sloping sections. After descending for just over 100 m, the entrance series drops into a horizontal and relatively narrow (3–4 m) streamway, the floor of which has several pools of standing water. The cave floor is predominantly rocky with areas of bedrock, boulders and coarse gravel (Fig.
The fish reside in small-sized (~3m x 4m) to large (>10m x 10m) pools. Although the invertebrate community in the cave is plentiful, it is not noticeably more abundant than that of many caves in Meghalaya. Amongst the terrestrial invertebrates were brown crickets (Eutachycines sp.), cellar spiders (Pholcidae) and fungus gnat larvae (Keroplatidae). Isopods were also frequently encountered including Cubaris sp. and Philoscia sp. Aquatic invertebrates included shrimp (Macrobrachium cf. cavernicola), snails (Paludomus sp.), pond skaters (Gerridae), and a few tadpoles. No significant bat roosts were encountered, and therefore no guano deposits or other obvious sources of nutrients were observed within the cave. It is conceivable that seasonal flood debris (bamboo, tree branches and leaf litter) carried into the cave from the surrounding forest provides the primary food source for the fish population. There is no plant growth in the caves and in the absence of bat guano, there is probably no other primary energy source in the habitat.
Unlike Um Ladaw, the Krem Chympe, where one of the paratypes were collected, is a broadly horizontal river cave, with a massive tunnel of deep water, and various small waterfalls/dams inside. Neolissochilus pnar occurs here in pools in a side passage. The biodiversity in this cave comprises of fish (Garra sp.), shrimps (Macrobrachium sp.), and tadpoles. Further details and photographs of both Um Ladaw and Chympe caves are available from
Phylogenetic analysis based on ML analysis revealed that the new species forms a distinct clade, and the sister taxon to a clade containing two other species of Neolissochilus, namely N. hexagonolepis (M’Clelland) and N. hexastichus, both from the Brahmaputra River system of northeast India (Fig.
Maximum likelihood analysis of Neolissochilus, Tor and Naziritor based on concatenated mitochondrial COI, cytb and 16S sequences. Bootstrap values based on 1000 iterations are shown along the nodes. Garra species are used as outgroup. Live specimen of Neolissochilus pnar is shown in inset (specimen not collected) (Photo: Uros Aksamovic).
Maximum likelihood analysis of available COI sequences for Neolissochilus with Tor putitora as outgroup. Species delimitation based on ASAP is shown as a dashed circle surrounding the phylogenetic tree. Clade containing Neolissochilus pnar is highlighted in blue. Clades containing morphologically identified putative topotypes of Neolissochilus hexagonolepis and N. hexastichus sensu
The limestone caves of Meghalaya, in northeastern India harbour a remarkable diversity of subterranean taxa (
The genus Neolissochilus represents a poorly-known group of medium- to large-sized cyprinids, with currently 31 species, distributed across South and Southeast Asia (Fricke et al. 2023). Despite its cultural and commercial importance, there have been no comprehensive studies on the taxonomy or systematics of this group since the description of the genus (
The sister taxon of Neolissochilus pnar (Fig.
One of the paratypes (KUFOS.F.2022.703) was collected from about 1.2 km inside the adjacent Krem Chympe cave, had a slightly different appearance including distinctly larger eyes and scalation pattern than fish of a similar size from those in the Krem Um Ladaw. Although the paratype of Neolissochilus pnar from the Krem Chympe cave was identical to the COI barcoding region and partial 16S genes of holotype and paratype from Krem Um Ladaw, there is a 1.4% raw genetic distance in the cyt b gene between the two populations.
The eye size in individuals of the Krem Um Ladaw population reduces as fish size increases. The smallest individuals have distinct, but atrophied eyes, which then become less distinct in larger individuals, and appear to be entirely absent in the largest individuals. This pattern appears consistent over all individuals photographed and videoed in the Krem Um Ladaw.
Genetic data for Neolissochilus subterraneus, the only other known subterranean species of Neolissochilus described from Tham Phra Wang Daeng cave in Thailand (
Principal component analysis biplot of factor scores and factor loadings of morphometric data of the two subterranean species of Neolissochilus. Factor scores are shown as scatter of points and factor loadings are shown as arrows. Percentage variation, out of total variation in the data, explained by each principal component, are provided in parenthesis.
Neolissochilus pnar is the largest known troglobitic species by a considerable margin. There has been the view that troglobitic adaptations are a consequence of the limited food availability in cave habitats. The need to locate sparse food reserves is thought to drive the development of the enhanced chemosensory capabilities typical of troglobites (
Despite the ichthyofaunal richness in aquifers and caves on the Indian subcontinent, there is only a limited number of studies dealing with their diversity and distribution. Recent descriptions of not only new species (Choudhury et al. 2017;
The authors thank Daniel Harries, Heriot Watt University, Edinburgh, UK for his generous support and help for this study. Studies of the Meghalayan cave biota would not have been possible without the support of the Meghalayan Adventurers Association (led by Brian Kharpran-Daly), which has supported international teams of cave explorers under the banner of the ‘Caving in the Abode of the Clouds Project’ for the last thirty years. All participants in this project deserve acknowledgement, but are too numerous to mention. Key individuals include Thomas Arbenz (Switzerland) and Simon Brooks (UK). Duwaki Rangad is grateful to the Shri S.M. Sahai IFS, Principal Chief Conservator of Forests (Bio-Diversity & Wildlife) and Chief Wildlife Warden, Government of Meghalaya, Shillong for necessary permissions for survey and collection of specimens (FWC/Research/128). Rajeev Raghavan thanks Mithun Sukumaran (Department of Aquatic Biology, University of Kerala, Thiruvananthapuram, India), Arya Sidharthan (KUFOS, Kochi) and Arjun C.P (MARC, Kannur) for their help and support during the study; Rohan Pethiyagoda (Australian Museum, Sydney) and Ralf Britz (Senckenberg Collections, Dresden, Germany) for useful discussions, guidance and support. Lukas Rüber (Natural History Museum, Bern, Switzerland), an anonymous reviewer, and the handling editor provided constructive comments and useful suggestions that greatly improved the paper.
Supplementary informations
Data type: .docx
Explanation notes: Tables S1. GenBank details for sequences used for Figure 3. — Table S2. Statistics for partition scheme and substitutional model analysis for maximum likelihood analysis provided in Figure 3. — Table S3. GenBank details for COI sequences used for Figure 4. — Table S4. Statistics for partition scheme and substitutional model analysis for maximum likelihood analysis provided in Figure 4. — Table S5. Statistics of barcode gap analysis using ASAP. Best partition is shown in the first row highlighted in grey.