Zoosyst. Evol. 100 (3) 2024, 897-911 | DOI 10.3897/zse.100.125798 Ate BERLIN Revalidated after having been described more than a century ago: Calamaria berezowskii Giinther, 1896 (Squamata, Colubridae) from Sichuan, Southwestern China Ya-Ting Liang!*°, Zi-Dan Huang! *°, Li Ding*, Gernot Vogel®, Natalia B. Ananjeva®, Nikolai L. Orlov®, Sheng-Chao Shi’, Zheng-Jun Wu!*?, Ze-Ning Chen! 1 Key Laboratory of Ecology of Rare and Endangered Species and Environmental Protection (Guangxi Normal University), Ministry of Education, Guilin 541006, Guangxi, China Guangxi Key Laboratory of Rare and Endangered Animal Ecology, Guangxi Normal University, Guilin 541006, Guangxi, China School of Life Sciences, Guangxi Normal University, Guilin 541006, Guangxi, China Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu 610041, Sichuan, China Society for Southeast Asian Herpetology, D-69115 Heidelberg, Germany Division of Herpetology, Zoological Institute, Russian Academy of Sciences, St. Petersburg 199034, Russia N DO oO fF W PY Hubei Engineering Research Center for Protection and Utilization of Special Biological Resources in the Hanjiang River Basin, School of Life Science, Jianghan University, Wuhan 430056, China https://zoobank. org/453CA 1 DC-07E7-40E8-9979-243C9B9I48373 Corresponding authors: Natalia B. Ananjeva (nananjeva09@gmail.com); Ze-Ning Chen (chenzn@gxnu.edu.cn) Academic editor: Justin Bernstein # Received 19 April 2024 # Accepted 9 June 2024 Published 5 July 2024 Abstract The reed snakes of the genus Calamaria Boie, 1827 are one of the largest groups of Asian snakes, distributed from northeast India to the Maluku Islands of east Indonesia. Recent research on the genus in China has revealed that the species diversity of the group was underestimated. In this study, morphological comparisons and mitochondrial DNA analysis revealed that a junior synonym of C. pavimentata Dumeéril, Bibron & Dumeéril, 1854 — Calamaria berezowskii Gunther, 1896 is valid, hence we redescribed and re- covered the validity of C. berezowskii. This species can be distinguished from other congeners by the combination of the following characters: four supralabials; one preocular; rostral shield width larger than height; mental not touching anterior chin shields; eye diameter less than the distance from eye to mouth edge; less than 1/2 of the posterior chin shield meets in the midline; dorsal scales reduced to six rows at tail; indistinct light ring present in the nuchal region or a more or less distinct yellowish collar. Phylogeneti- cally, this species is sister to C. pavimentata, with significant genetic differences (0.190) on mitochondrial gene Cyt b. Key Words Calamaria berezowskii, Colubridae, morphology, phylogenetics, taxonomy Introduction 1854; C. pavimentata Duméril, Bibron & Duméril, 1854; C. lovii Boulenger, 1887; C. septentrionalis Boulenger, Calamaria Boie, 1827 is the largest group of the colubrid 1890; C. prakkei Lidth De Jeude, 1893; C. buchi Marx & subfamily Calamariinae (Reed Snakes), containing more than 68 species (Uetz et al. 2023). In mainland Southeast Asiaand China, 20 species of Calamaria are known at pres- ent, including C. /umbricoidea Boie, 1827; C. albiventer (Gray, 1834); C. schlegeli Duméril, Bibron & Dumeéril, Inger, 1955; C. yunnanensis Chernov, 1962 (Yeung et al. 2022; Uetz et al. 2023); C. thanhi Ziegler & Quyet, 2005; C. gialaiensis Ziegler, Nguyen & Nguyen, 2009; C. san- gi Nguyen, Koch et Ziegler, 2009; C. abramovi Orlov, 2009; C. concolor Orlov, Nguyen, Nguyen, Ananjeva & Copyright Liang, Y.-T. et al. 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. 898 Ho, 2010; C. andersoni Yang & Zheng, 2018; C. dominici Ziegler, Tran & Nguyen, 2019; C. strigiventris Poyarkov, Nguyen, Orlov & Vogel, 2019; C. nebulosa Lee, 2021; C. arcana Yeung, Lau & Yang, 2022; C. jinggangensis Cai, Jiang, Wu, Huang, Fei & Ding, 2023. There were only three species of genus Calamaria re- corded in China 26 years ago: C. pavimentata, C. septen- trionalis and C. yunnanensis (Zhao et al. 1998). Recently, three more species were described: C. andersoni, C. ar- cana, C. jinggangensis (Yang and Zheng 2018; Yeung et al. 2022; Cai et al. 2023; Uetz et al. 2023). The species diversity of the genus in China is suggested to be under- estimated, and those widely distributed species should be re-evaluated. For example, the widely recorded species C. pavimentata has a synonym, C. berezowskii Gunther 1896, which was described based on two specimens from Lun-ngan-fu (7% )ff Long’an Fu, now J 2244 Long’an Town of */ 70 Pingwu County) of Sze-chuen (Sichuan Province of China) (Gunther 1896). During the scientific expeditions in the Gongga Moun- tains, Luding County, Sichuan Province, between 2017 and 2022, we collected three Calamaria specimens. Mor- phological and molecular phylogenetic analyses showed that these specimens represent a species that differs from all currently recognized congeners of the genus. However, 100°0'0"E 35°0'0"N 30°0'0"N 100°0'0"E Liang, Y.-T. et al.: Calamaria berezowskii Gunther, 1896 the morphology of these specimens matches the original description of C. berezowskii. Besides, Luding County and Pingwu County both are located on the eastern slope of Qinghai-Tibet Plateau bordering the Sichuan Basin. Therefore, we identify the specimens as C. berezowshkii, and the species is revalidated and redescribed here. Materials and methods Sampling Three specimens of Calamaria were collected from east- ern slope of Mt. Gongga, Mox1 Town, Luding County, Ganzi Tibetan Autonomous Prefecture, Sichuan Province, China during 2017 to 2022 (Fig. 1). One adult female GXNU DLR195 (29.645105°N, 102.111076°E, 1736 m a.s.l.), one adult male GXNU DLR194 (29.615872°N, 102.109208°E, 1680 m as.l.), one juvenile female GXNU_ 20221215002 (29.641749°N, 102.110340°E, 1827 ma.s.l.), collected by Xu Zhang on 31 August 2018, Li Ding on 21 August 2017 and Congcong Du on 15 De- cember 2022. The three specimens were fixed and stored in 80% ethanol and deposited at the School of Life Sci- ences, Guangxi Normal University. 105°0'0"E 110°0'0"E 3 5°0'0"N ~ Low: -505 105°0'0"E 110°0'0"E Figure 1. Distribution map of Calamaria berezowskii. Red star indicates the type locality, purple triangle indicates the distribution found in this study. zse.pensoft.net Zoosyst. Evol. 100 (3) 2024, 897-911 Molecular phylogenetic analysis Genomic DNA was extracted from muscle or liver tissue taken from the collected specimens using the TIANamp Marine Animals DNA Kit (TIANGEN Biotech). We sequenced a fragment of the mitochondrial gene cyto- chrome 5 (Cyt 5) using the primer pair L14910/H16064 (Burbrink et al. 2000). PCR amplification was per- formed in 25 pl reactions containing 12.5 wl 2x SanTaq PCR Master Mix (with Blue Dye), 10 pl ddH,O, 1 pl F-primers, 1 ul R-primers and 0.5 ul DNA template. The PCR conditions were an initial denaturing step at 94 °C for 5 min, followed by 35 cycles: denaturing at 94 °C for 30 s, annealing at 48.5 °C for 45 s, an extension step at 72 °C for 35 s; and a final extension of 72 °C for 10 min. Finally, the PCR products were sent to Sangon Biotechnologies Co., Ltd. (Shanghai, China), where the purified PCR products were sequenced using the same forward and reverse primers. Sequences for comparison of available species were downloaded from GenBank (Table 1). We evaluated and manually corrected the sequencing peak maps of the three sequences obtained, downloaded 29 Cyt b sequences of 12 species of the genus Calamaria that have been reported in Southeast Asia and China in Table 1. DNA sequences used in this study. 899 GenBank, and selected three species Elaphe quatuorlin- eata, Lycodon rufozonatus and Orientocoluber spinalis as outgroups (Yang and Zheng 2018; Yeung et al. 2022), for a total of 35 sequences to be used for multiple sequence comparisons. All sequences were aligned with other re- trieved sequences in the same gene loci by using software MEGA 11 (Kumar et al. 2018). Phylogenetic trees were constructed based on mitochondrial gene Cyt 6b using maximum likelihood (ML) and Bayesian inference (BI). Maximum likelihood analysis was conducted in RAXxML v8.2.4 (Stamatakis 2014). Confidence intervals were de- termined with 1000 bootstrap replicates utilizing the rap- id bootstrap option under the GTR+GAMMA substitu- tion model. Bayesian inference analysis was conducted in MrBayes 3.2 (Ronquist et al. 2012). PartitionFinder 2.1.1 software was used to select the most suitable nucleotide substitution models for Cyt b sequence data based on the Bayesian information criterion (BIC): GTR+G, HKY+G and GTR+G (Lanfear et al. 2012; Lanfear et al. 2017). We ran our analyses for 20 million generations with the chains, sampling every 1000 generations using the Mar- kov Chain Monte Carlo (MCMC). After removing out- group taxa, MEGA11 (Kumar et al. 2018) was used to calculate uncorrected pairwise sequence divergence be- tween the Calamaria species. No. Taxa Voucher Locality Cyt b Reference 1 C. septentrionalis HSR19100 Mt. Huangshan, Anhui, China 0Q354842 Cai et al. 2023 2 C. septentrionalis HS11145 Mt. Nanling, Guangdong, China 0Q354840 Cai et al. 2023 3 C. septentrionalis DL2021610-1 Huangsha, Guangxi, China 0Q354838 Cai et al. 2023 4 C. septentrionalis KFBG14506 Hainan Island, China MH445956 Yang and Zheng 2018 5 C. septentrionalis ROM35605 Nguyén Binh, CaoBang, Vietnam AF471081 Lawson et al. 2005 6 C. septentrionalis ROM35597 Cao Bang, Vietnam KX694890 Alencar et al. 2016 7 C. pavimentata KFBG14507 Ningming, Guangxi, China MH445957 Yang and Zheng 2018 8 C. andersoni HSR20101 Dehong, Yunnan, China 00354844 Cai et al. 2023 9 C. andersoni HSR20181 Tengchong, Yunnan, China 0Q354845 Cai et al. 2023 10 ‘C. andersoni SYSr001699 Yingjiang, Yunnan, China MH445955 Yang and Zheng 2018 11 C. yunnanensis ROM41547 Simao, Yunnan, China KX694891 Zaher et al. 2009 12. C. yunnanensis YPx503 Unknown JQ598922 Grazziotin et al. 2012 13. -C. arcana HS17082 Mt. Dawu, Guangdong, China 0Q354835 Cai et al. 2023 14. ~C. arcana KFBG14611 Mt. Dadongshan, Guangdong, China ON482335 Yeung et al. 2022 15 C. arcana GP9975 Yongxing, Hunan, China OP980549 Cai et al. 2023 16 ~=~C. arcana DLR199 Mt. Wuyi, Fujian, China 0Q354834 Cai et al. 2023 17 _ C. jinggangensis DL20200725 Mt. Jinggangshan, Jiangxi, China 0Q354830 Cai et al. 2023 18 C. jinggangensis DL20200625-2 Mt. Jinggangshan, Jiangxi, China 0Q354831 Cai et al. 2023 19 _‘C. jinggangensis DL20200625-3 Mt. Jinggangshan, Jiangxi, China 0Q354832 Cai et al. 2023 20 —C. jinggangensis DL20200625-4 Mt. Jinggangshan, Jiangxi, China 0Q354833 Cai et al. 2023 21 ~~ ‘C. muelleri RMB1283 Gowa, South Sulawesi, Indonesia MT819391 Weinell et al. 2021 22 ~~ C. muelleri TNHC58955 Gowa, South Sulawesi, Indonesia MT819390 Weinell et al. 2021 23 ~~ ‘C. lumbricoidea USMHC1560 Penang, Peninsular, Malaysia MN338526 Quah et al. 2019 24 ~ C. palavanensis KU311411 Mt. Mantalingahan, Palawan, Philippine MT819387 Weinell et al. 2021 25 C. palavanensis KU309445 Barangay lrawan, Palawan, Philippine MT819386 Weinell et al. 2021 26 ~—sC..: gervaisii KU334485 Municipality, llocos Sur, Philippines MT819385 Weinell et al. 2021 27 ~—=«C.: gervaisii KU324661 Puguis, Benguet, Philippines MT819384 Weinell et al. 2021 28 ~=-C.. schlegeli LSUHC10278 Perak, Peninsular, Malaysia MN338525 Quah et al. 2019 29. ~=C. nebulosa FMNH258666 Phongsaly, Laos MN338524 Quah et al. 2019 30 ~~‘ C.. berezowskii GXNU DLR194 Mt. Gongga, Sichuan, China PP747047 This study 31. ‘C. berezowskii GXNU DLR195 Mt. Gongga, Sichuan, China PP747048 This study 32 C. berezowskii GXNU 20221215002 Mt. Gongga, Sichuan, China PP747049 This study 33. Elaphe quatuorlineata LSUMZ40626 Turkey, European Turkey AY486931 Nagy et al. 2004 34 — Lycodon rufozonatus LSUMZ44977 Unknown AF471063 Lawson et al. 2005 35 Orientocoluber spinalis MVZ211019 Yinnan, Ningxia, China AY486924 Nagy et al. 2004 zse.pensoft.net 900 Morphological analysis Terminology and descriptions follow the views of Inger and Marx (1965) and Ziegler et al. (2008). Body and tail length were measured with a tape ruler to the nearest 1 mm: total length (TL), from the tip of snout to the tip of tail; snout-vent length (SVL), from the tip of snout to posterior margin of cloaca; tail length (TaL), from posterior margin of cloaca to the tip of tail. Other measurements were con- ducted with a digital caliper to the nearest 0.1 mm: head length (HL), from the snout tip to the posterior margin of the mandible; head width (HW), measured at the widest part of the head on posterior side; head height (HH), at the maximal highest part of the head; the eye horizontal diameter (EyeD); and eye-mouth distance (Eye-Mouth D), measured from the anterior point of the eye to the mouth gap. Ventral scales (VEN) were counted according to Dowling (1951). The enlarged shield(s) anterior to the first ventral were regarded as preventral(s). The number of dorsal scale rows (DSR) are given at one head length be- hind head, at midbody, and at one head length before vent, respectively. The tail tip was not included in the number of subcaudal scales (SC). Sex was determined by examining the presence or absence of hemipenis. Type specimen (lectotype ZISP 8823) was examined for morphological comparisons. Relevant morphological data of other Calamaria species were obtained from ex- amined specimens (Appendix 1) and literatures. Institution acronyms FMNH: Field Museum of Natural History, Chicago. GXNU: Guangxi Normal University. KFBG: Kadoorie Farm and Botanic Garden. KU: University of Kansas Biodiversity Institute. LSUHC: La Sierra University Liang, Y.-T. et al.: Calamaria berezowskii Gunther, 1896 Herpetological Collections. LSUMZ: Louisiana State University Museum of Natural Science. MVZ: Muse- um of Vertebrate Zoology. RMB: Rafe M. Brown field tag (specimen deposited in Museum Zoologicum Bo- goriense, Indonesia). ROM: Royal Ontario Museum, Canada. SYS: Biological Museum of Sun Yat-sen Uni- versity, Guangzhou, China. TNHC: Texas Natural Histo- ry Collections. USMHC: Universiti Sains Malaysia Her- petological Collection, Malaysia. ZFMK: Zoologisches Forschungsmuseum Alexander Koenig, Bonn, Germany. ZISP: Zoological Institute, Russian Academy of Scienc- es, St. Petersburg, Russia. Molecular results The DNA dataset contains 35 samples with a total of 1105 base pairs. The phylogenetic topologies that resulted from Bayesian Inference (BI) and Maximum likelihood analy- sis (ML) are generally consistent in phylogenetic structure (Fig. 2). Phylogenetically, the species Calamaria from Mt. Gongga strongly clustered into a single lineage with high nodal support (Fig. 2, BPP 1.00 / UFB 100). Subsequently, it clustered with C. pavimentata KFBG14507 from Ning- ming, Guangxi (BPP 0.96 / UFB 94), and was clustered into the same clade with C. arcana, C. jinggangensis, C. septen- trionalis, C. pavimentata, C. andersoni and C. yunnanensis. The uncorrected pairwise p-distances of Cyt b sequenc- es between the specimen of Mt. Gongga and the other 12 congeners included in the study were 17.1—31.2%, with the minimum value observed in the comparison with sequences of C. arcana (P = 17.1%) (Table 2). This di- vergence is clearly among interspecies level since these levels of divergences are distinctly higher than those ob- served between two other well distinguished species of C. andersoni and C. yunananensis (P = 12.2—12.6%). The Table 2. Uncorrected p-distances between Calamaria species based on 1105 base pairs from the mitochondrial genes Cyt 6. The serial numbers in Table 2 are consistent with those in Table 1. No. Taxa 1-6 7 8-10 11-12 13-16 1-6 _C. septentrionalis 0.000-0.038 7 C. pavimentata 0.174-0.194 8-10 C. andersoni 0.130-0.148 0.246 0-0.019 11-12 C. yunnanensis 0.159-0.164 0.237 0.122-0.126 0.000 13-16 C. arcana 0.090-0.116 0.179- 0.181-0.206 0.183- 0.007- 0.199 0.209 0.034 17-20 C. jinggangensis 0.095-0.103 0.180 0.161-0.171 0.158 0.062- 0.074 21-22 C. muelleri 0.202-0.223 0.266 0.228-0.238 0.261- 0.178- 0.273 0.198 23 C. lumbricoidea 0.243-0.248 0.258 0.244-0.256 0.251 0.228- 0.239 24-25 C. palavanensis 0.187-0.222 0.231- 0.227-0.243 0.254- 0.207- 0.242 0.266 0.229 26-27 C. gervaisii 0.178-0.207 0.243- 0.202-0.212 0.243- 0.197- 0.258 0.244 0.223 28 C. schlegeli 0.232-0.243 0.263 0.258 0.271 0.217- 0.228 29 C. nebulosa 0.186-0.196 0.214 0.172 0.173 0.162- 0.182 30-32 C. berezowskii 0.187- 0.190 0.216- 0.208 0.171- 0.203 0.226 0.200 zse.pensoft.net 17-20. 21-22 23 24-25. 26-27 28 29 30-32 0.000 0.169- 0.007 0.178 0.228 0.164 0.173- O0.167- 0.166- 0.022 0.192 0182 0.181 0.159- 0173 0197- 0112 0.100 0.163 0.183 0.223 0.174 0.217 0.182- 0205 0170 0.186- 0.192 0.175 0.218 0.176 0.208 0.193 0.202- 0187- 0.197 0.219 0.207 0.234 0.202 0.223- 0.255 0.254 0.254 0.312 0.176 0.000 0.234 0.277. 0.271 Zoosyst. Evol. 100 (3) 2024, 897-911 901 1.00/100 Peps iane Pitted bee Atat i sapteveotss KFBG14611 aa bein | GP9975 C. arcana 1.00/99 ~ DLR199 1.00/100 eng ~ DL20200625-4 ~ DL20200625-3 C. jinggangensis 1.00/100 |. ~ DL20200625-2 0.99/97 | en DL20200725 Sas ttites PSE CAAA B83 tc 8b Res UE HSR19100 0.79/79 enn HS11145 C. septentrionalis 1.00/100 OPD2 Fteteacttenbiceth DL2021610-1 x posses ROMIZ5605 1.00/100 noi 1.00/99 poo HSR20101 con 1.00/100 ~ HSR20181 C. andersoni 1.00/100 om SYS1001699 ~ ROM41547 : 1.00/100 b----- YPx503 2 ann enen sts 1.00/1 00 iajhhs piste ee GXNU DLR1 94 1.00/100 0.64/88 Pe GXNU DLR195 C. berezowskii 0.96/94 ~ GXNU 2022121500 KFBG14507 | C. pavimentata eae FMNH258666 IC. nebulosa 10008 ers lc palavanensis 1.00/100 ~ KU309445 0.99/96 » KU334485 he Patten le gervaisil 0.99/99 1.00/100 ~ RMB1283 . a3 aint saoss | C. muelleri - USMHC1560 1 C. lumbricoidea ~ LSUHC10278 IC. schlegeli 0.96 ooo LSUMZ40626 | Elaphe quatuorlineata VUES OYSTELOT GPa | eee | EE ethics Pe tale CELL OPS nw Lk Bierseln Oe LSUM/744977 | Lycodon rufozonatus ~ MVZ211019 | Orientocoluber spinalis Figure 2. Phylogenetic relationships of Calamaria species based on Cyt 5. Bayesian posterior probabilities (BPP) and Ultrafast bootstrap approximation (UFB) were denoted beside each node (those lower than 75/0.7 were denoted as “-”). Calamaria specimens from Mt. Gongga, Luding, Sich- uan represent an evolutionary distinct lineage. Morphological results Morphology of the three Calamaria specimens from Luding County, Sichuan province matches the original description of C. berezowskii Gunther 1896 and the mor- phological characters of lectotype ZISP 8823 (Figs 3, 4, 5, 6). They share the following characters: one preocular and one postocular; presence of a more or less distinct yellowish collar; similar number of ventral scales (149-— 165 vs 155-167) and subcaudal scales (16—25 vs 14—22); presence of a faint dark line on each side of the back; tail blunt and conical at the end. Besides, Luding County and the type locality Pingwu County are both located on the eastern slope of Qinghai-Tibet Plateau bordering the Si- chuan Basin, with a geographic distance of about 350 km (Fig. 1). Hence, the Luding specimens were identified as conspecific with C. berezowskii. Morphological comparisons based on some characters between Calamaria berezowskii and other known species of the genus from mainland Southeast Asia and China were listed in Table 3. It can be distinguished by a com- bination of the following characters: four supralabials; preocular present; rostral shield width greater than height; mental not touching anterior chin shields; eye diameter less than distance from eye to mouth edge; less than 1/2 of the posterior chin shield meets in the midline; dorsal scales reduced to six rows at tail; presence of indistinct light rings in the nuchal region or a more or less distinct yellowish collar; and the two outermost dorsal scale rows light khaki, with upper margins partly dark pigmented. zse.pensoft.net 902 Liang, Y.-T. et al.: Calamaria berezowskii Gunther, 1896 Table 3. Morphological comparisons of species of Calamaria from mainland Southeast Asia and China. “1+***” indicates preven- trals + ventrals. Entries marked with ‘—’ are not available. Species and width Supralabials End of tail paraparietal only in anteriorly (0) and eye-mouth distance from (0) anterior chin shields Comparison between eye diameter Number of scales surrounding Mental touching (1) or separated Posterior chin shields meeting in midline (1), diverging or meeting Comparison between rostral height Preocular: present (1) or absent (0) C. berezowskii C. pavimentata Sharp point Ventrals Subcaudals Total length s =] o o ea] > —_ 0: x oO = zo} o = Zo} je) = = above position of subcaudal on tail w co © = © (ae) = 4 Obtuse point 24+149- | 153-— | 22-25 Not reduced | 248-290 | 123-305 155 2+165 Last-13" Dorsal scales reducing to four rows | Female Male Female 115485 subcaudal C. septentrionalis Broadly rounded 8-9 | 125-168 | 137-206} 3-33 | 8-20 8-9 | 148-166 | 168-188] 15-19] 6-11 Not reduced | 111-344 |117-384 C. yunnanensis H 1/4 Obtuse point 170-176 192 | 20-22 Not reduced | 144-303.2] 36.5 C.lumbricoidea |HW subcaudal C. albiventer H