Next Article in Journal
Histovariability and Palaeobiological Implications of the Bone Histology of the Dromornithid, Genyornis newtoni
Next Article in Special Issue
Bat Diversity in Cat Ba Biosphere Reserve, Northeastern Vietnam: A Review with New Records from Mangrove Ecosystem
Previous Article in Journal
Traits to Differentiate Lineages and Subspecies of Aegilops tauschii, the D Genome Progenitor Species of Bread Wheat
Previous Article in Special Issue
Genetic Introgression and Morphological Variation in Naked-Back Bats (Chiroptera: Mormoopidae: Pteronotus Species) along Their Contact Zone in Central America
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 13
Issue 5
10.3390/d13050218
diversity-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

Not the Cryptic Species: Diversity of Hipposideros gentilis (Chiroptera: Hipposideridae) in Indochina

by
Alexander P. Yuzefovich
1,
Ilya V. Artyushin
1,* and
Sergei V. Kruskop
2,*
1
Department of Vertebrate Zoology, Moscow State University, Leninskiye Gory 1 bld. 12, 119234 Moscow, Russia
2
Zoological Museum, Moscow State University, Bolshaya Nikitskaya 2, 125009 Moscow, Russia
*
Authors to whom correspondence should be addressed.
Diversity 2021, 13(5), 218; https://0-doi-org.brum.beds.ac.uk/10.3390/d13050218
Submission received: 9 April 2021 / Revised: 12 May 2021 / Accepted: 12 May 2021 / Published: 19 May 2021
(This article belongs to the Special Issue Diversity Aspects in Bats: Genetics, Morphology, Community Structure)
Browse Figures
Review Reports Versions Notes

Abstract

:
We present here the result of phylogenetic analysis for Vietnamese Hipposideros gentilis specimens using 7 nuclear genes and one mitochondrial gene. The complex distribution of divergent mitochondrial DNA lineages contradicts, at least in part, nuclear and morphological data. The most likely explanation for this discordance is the historical hybridization between ancestral populations of H. gentilis and H. rotalis/H. khaokhouayensis. Our data supports the species status of H. gentilis, while only partially corroborating its previously proposed subspecies delimitation. We suggest the lowland forest populations from south Vietnam may correspond to their own subspecies. At the same time, the close phylogenetic relationship and morphological similarity of mountain forms from south and central Vietnam to the north Vietnamese populations make doubtful the subspecies status of H. gentilis sinensis.

Graphical Abstract

1. Introduction

Hipposideros is the core genus to Hipposideridae (Old World leaf-nosed bats), a bat taxon prominent in Southeast Asia. It is also one of the largest mammalian genera, retaining a highly complicated taxonomy even after having lost two of its former species groups, namely “commersoni” and “cyclops” [1,2]. Within the genus, Hipposideros bicolor is a group that comprises the most forms. Many of them are similar in morphology, and their taxonomic placement is often uncertain. In the last decade, significant efforts were dedicated to solving this puzzle (e.g., [3,4,5]).
One of the common Southeast Asian species of leaf-nosed bats from the “bicolor” group was widely referred to as Hipposideros pomona by Andersen in 1918 [6,7,8]. Morphologically, this form was once treated as a subspecies or group of subspecies within H. bicolor (Temminck, 1834) [9]. Thereafter H. pomona was considered a distinct species, including the forms gentilis (Andersen, 1918) and sinensis (Andersen, 1918) as partial synonyms [10]. Paracoelops megalotis (Dorst, 1947) was initially described as a separate species and genus from central Vietnam [11]. However, a review of the single type specimen concluded that it belongs to H. pomona [12]. Molecular studies have revealed that several geographically limited haplogroups are present in H. pomona populations. Distances between these haplogroups are high enough to assume species-level divergence. Moreover, the monophyly of H. pomona was questioned [13,14]. In the Indochina peninsula, four main lineages have been discovered, namely northern, central, and two in the south. One of the southern clades takes a sister position to the lineage which includes H. khaokhouayensis (Guillen-Servent, Francis, 2006) and H. rotalis (Francis, Kock, Habersetzer, 1999) [13]. Some of these haplogroups were treated as subspecies in an earlier study of H. pomona genetic variability in China [5]. The clades SM and SY were assigned to the subspecies sinensis and gentilis, respectively, with no suggestions for the “central” clade (SH). The southern haplogroups were not discussed in that paper.
On the other hand, the shape of baculum (os penis) is very similar in different genetic lineages of H. pomona [15] (see also [16]). The baculum is greatly reduced and stick-like, about 0.5 mm in length. It is strikingly different from that of H. bicolor, H. cineraceus (Blyth, 1853) and other related species [4]. At the same time H. pomona s. str. from southwest India is remarkably distinct from all other populations in bacular and cranial morphology [17]. Hipposideros gentilis is the next senior name, and in this case should be accepted as valid for all populations under current study.
Similarity in baculum shape, as well as the absence of obvious differences between populations in the skull morphometry and the structure of the nasal leaves (orig., unpubl.), cast a shadow of doubt upon the putative species status of the H. gentilis haplogroups. Here, we apply a multi-locus analysis of nuclear genes to address this subject.

2. Materials and Methods

2.1. Analyzed Specimens

Our research is based on the Zoological Museum of Moscow State University (ZMMU) collections. Tissue samples (muscle fragments) were taken from specimens preserved in ethanol. In total, 14 specimens of H. gentilis (H. “pomona”) were used, together with a few other species of the genus Hipposideros. Aselliscus stoliczkanus was taken as outgroup. A full list of specimens and obtained sequences is provided in Table 1 and Table 2. Also, additional sequences of cytb and four nuclear genes were downloaded from GenBank (see Table 2 and Appendix A).

2.2. Morphometric Analysis.

To understand how the observed genetic diversity is reflected in morphology, we performed a morphometric analysis using cranial and dental measurements.
A total of 102 specimens of the Hipposideros pomona/gentilis complex (dry or alcohol preserved skins with extracted skulls) were examined for morphometric comparison. We excluded several specimens due to missing measurements. Only 99 were included in the final analysis. The full list of these specimens is provided in Appendix C. Localities are shown in Figure 1. Acronyms of the processed collection repositories are: MNH—British Museum of Natural History (London, UK); MNHN—National Museum of Natural History (Paris, France); ROM—Royal Ontario Museum (Toronto, Canada); NHW—Vienna Museum of Natural History (Vienna, Austria); ZMMU—Zoological Museum of Moscow State University (Moscow, Russia); ZMB—Berlin Zoological Museum (Berlin, Germany).
Cranial and dental measurements were made under stereo microscope with digital calipers, rounding to the nearest 0.01 mm. The following cranial measurements were taken: greatest skull length (TL), condylo-canine length (CCL), skull width at the mastoid (MW), brain case width above the mastoids (BCW), occiput height (height from the lower margins of the occipital condyles to the highest point just above them; OH), zygomatic width (ZW), width of the postorbital constriction (POC), rostral width at the level of the anterorbital foramina (RW), length of rostrum in front of the anterorbital foramen (RL), width across the upper canines (CC), width across posterior upper molars (MM), length of the upper tooth row (CM), length of the upper molariform row, distance from P4 to the posterior molar (PM), longitudinal length of the upper canine (C), width of nasal opening (NO), lower tooth row length (cm), articular length of the mandible (MdL), and height of the mandible (MdH).
To assess the variation pattern of quantitative characters, Principal Component (PC) and Discriminant Function (DF) analyses were performed for the 20 craniodental measurements, using STATISTICA for Windows version 9.0 (StatSoft, currently TIBCO Software Inc., Palo Alto, CA, USA). DF analysis was used to calculate squared Mahalanobis distances between groups and the significance of inter-group differences. The training set for calculating the squared Mahalanobis distances for the DF analysis included six samples, namely “Myanmar”, “Central Thailand”, “NW Indochina”, “NE Indochina and SE China”, “C Indochina” and “S Indochina”. The third and fourth training sets were selected based roughly on the two haplogroup (SY and SM) ranges as they are described in [5]. Indian material, due to low specimens’ availability, was not treated as a training set, but was included in the analysis as “unidentified”.
We were not able to genotype all specimens used in the morphometric study. However, all geographic samples from Indochina used in both the PS and DF analyzes include specimens genotyped for one or both mitochondrial genes and for nuclear genes (see Appendix C and appropriate Figures).

2.3. Molecular Analysis

Genomic DNA was extracted from ethanol-preserved muscle tissue by standard method of phenol–chloroform deproteinization [19]. One mitochondrial gene cytochrome b (cytb; 1137 bp) and fragments of seven nuclear genes (ABHD11, 460 bp; ACOX2, 597 bp; COPS, 744 bp; RAG2, 1035 bp; ROGDI2, 509 bp; SORBS2, 569 bp; and THY, 565 bp) were sequenced. Primers were taken from previously published papers [20,21,22,23]. The thermal profile amplification reaction included 35 cycles and was set up as follows: denaturation for 30 s at 94 °C, annealing for 1 min (at 55 °C for cytb and 60°C for the nuclear genes), and elongation for 1 min at 72 °C. Predenaturation lasted 3 min at 94 °C, and the final elongation was 6 min at 72 °C. Automatic sequencing was performed on an ABI PRISM 3500xl sequencer (Applied Biosystems (ABI), Thermo Fisher Scientific, Waltham, Massachusetts, United States) using ABI PRISM® BigDyeTM Terminator v. 3.1 reagent kits (Applied Biosystems, United States) in the laboratory of Eurogen Joint Stok Company (Moscow, Russia). To obtain additional mitochondrial topology, a total of 117 specimens, including 95 sequences of H. gentilis and seven other Asian Hipposideros species were analyzed for patterns of genetic divergence in the COI region. These data were taken from publicly accessible projects housed by the Barcode of Life Data System (BOLD, Guelph, Ontario, Canada) [24]. A list of BOLD Process ID numbers is provided in Appendix B.
Obtained sequences were aligned by the MAFFT v7.307 algorithm [25], then trimmed and corrected manually, after which heterozygous sites were marked by IUPAC ambiguity symbols for tree construction and phased in DNA 6.12.03 [26] for network building. Indels were replaced by “N.” The reconstruction of phylogenetic trees was carried out using maximum likelihood (ML) in IQ-TREE v2.0.6 [27,28]. The built-in algorithm was applied to select optimal models and partitioning scheme [29]. Bayesian analysis (BA) was also performed using the MrBayes 3.2.6 program [30,31,32]. Partitioning schemes and models for BA were determined in Partition-Finder v. 2.1.1 [33] using a BIC criterion with a greedy search method. The initial partition scheme included separate partitions for each gene. For mitochondrial genes, codon positions were used as the initial partitioning scheme.
Supertree was computed using the MRP method (matrix representation with parsimony) [34,35,36]. Using the alignment of individual nuclear genes, bootstrap samples of 1000 ML trees were constructed. The resulting topologies were combined using the maximum parsimony method, implemented in the phytools v0.7.70 [37] package in the R v 3.5.2 [38] environment. The reliability of the topology was calculated by applying a bootstrap procedure to total tree sample.

3. Results

3.1. Pattern of Variation in the mtDNA Marker

Phylogeny reconstructed from cytb sequences mostly corroborates previously published COI gene data [13]. Basal nodes and relationships between different species groups have low support. The monophyly of the species groups generally is well supported. Topologies of Bayes and ML trees didn’t contradict each other at highly supported nodes (Figure 2).
Hipposideros gentilis is non-monophyletic by this gene and is divided into two well-supported clades (Figure 2, Figure S1; Table 3, Table S1). The first clade is by itself paraphyletic with respect to the H. khaokhouayensis/rotalis lineage and includes samples from the southern lowland forests of Vietnam (Dong Nai, Tuy Phong, and Bihn Chau provinces) and from the southern part of the Vietnamese central highlands (Gia Lai province). Within the second clade of H. gentilis, specimens from Vietnam, Laos and China are represented. The Haplogroup from the northern half of central Vietnam forms a sister clade to the animals from Laos, Hainan (SH clade in [5]) and one specimen from Yunnan province, China. Animals from north Vietnam (Gulf of Tonkin) are grouped with those from Guangdong (SM clade). Such relations somewhat contradict the results of [5] and suggest a possible overlap in distribution of haplogroups. Both the central Vietnam and Gulf of Tonkin clades are monophyletic with high support.

3.2. Pattern of Variation in nuDNA Markers

Analyses of individual nuclear genes mainly demonstrate insufficient phylogenetic signal (Figure S2, Tables S2, S3). Three of seven genes (RAG2, ROGDI, THY) yielded monophyletic clades for the analyzed samples of H. gentilis. However, only the THY gene had a good bootstrap support value for this branch (87%). It is unclear which of the two species (H. ater, H. cineraceus) is closest to H. gentilis. The other three genes (ABHD11, ACOX2, SORBS2) did not support H. gentilis monophyly, and we placed H. khaokhouayensis, H. bicolor or H. cineraceus respectively within this clade. Not enough sequences were obtained for the COPS gene to draw any meaningful conclusions.
Low diversity levels and dataset incompleteness precluded reliable distance computation, although the substitution count between geographic samples of H. gentilis is generally lower than those between species in the “bicolor” group. We also found heterozygous SNPs in sites which bear differences between the three Vietnamese H. gentilis lineages.
These data definitely contradict to the results obtained from mitochondrial genes and generally indicate closer relationships of H. gentilis populations.
However, the tree built with the MRP algorithm on the combined data for all the seven nuclear genes supports the monophyly of H. gentilis and provides well-resolved topology within its clade (Figure 3). The position of this clade relative to H. bicolor, H. ater and H. khaokhouayensis is unresolved, and H. cineraceus is a sister branch to all of the above. Contrary to mitochondrial data, animals from the south Vietnam lowlands form their own branch, which is the most basal within the H. gentilis radiation.
Specimens from the Vietnamese central highlands are separated from the south Vietnam lowland populations and tend to aggregate with the specimens from other mountainous populations from the Dalat Plateau and central Vietnam. Animals from the north-east of Vietnam (Gulf of Tonkin) are paraphyletic with respect to the mountain populations.

3.3. Patterns of Morphological Variation

The morphometric data are moderately factorized: the total contribution of the first four factors to the total variance is about 70%, which means relatively low overall variation of cranial proportions in H. pomona s. lato. In the space of the first two factors, most geographic samples overlap highly, except for the Indian specimens, representing H. pomona s.str. (Figure 4). This agrees with the species-level dissimilarity between H. pomona and H. gentilis. The third factor greatly reduces the overlap between the samples from northern Indochina (both western and eastern) and the other analyzed samples. Then, the fourth factor separates the sample from southern Indochina from all the others (Figure 5). The most significant contribution to the first PC was made by CM, PM, C and cm variables (i.e., by the tooth row length); by MW, BCW, and CCL (skull width and overall size) to the second; and by CC and, to a lesser extent, other measurements related to skull width and mandible size (Table 4) to the fourth.
The results of the discriminant function analysis clearly show the high significance of the differences among all the training sets, with one exception: animals from western and eastern north Indochina are separated insignificantly, suggesting their taxonomic equality (Figure 6; Table 5). Surprisingly, the specimens set from Myanmar turned out to be one of the most distinct. The training set from southern Indochina is well-separated from all the rest by the values of the first two canonical variables. It shows some degree of similarity with the set from central Thailand and overlaps slightly with it, but not the others. The specimens from the Dalat Plateau and Gia Lai province, which are similar in mitochondrial sequences to the south Vietnamese ones, show similarities with different training sets (more often with central Indochina), but not with southern Indochina.

4. Discussion

Taxonomical issues regarding the status and delimitation of morphologically similar tropical forms are always quite complicated and cannot be solved by considering one particular trait alone. According to our results, a simultaneous analysis of a few (here seven) nuclear genes can provide a highly supported and reliable tree topology that cannot be obtained from any of the markers separately. The resulting phylogenetic relationships within the H. gentilis clade and between this species and other congenerics differs substantially from the previously known mtDNA tree topology. The first remarkable feature of the mitochondrial phylogeny is paraphyletic position of the different H. gentilis lineages, which suggests the assignment of species status to at least some of them. However, our nuclear data definitely support the monophyly of H. gentilis, which diminishes the likelihood of finding cryptic taxa within this complex. Such discordance may imply historical introgression events.
Specimens from the mountain and lowland regions of southern Vietnam arranged themselves into distinct, but related, haplogroups (also related to the H. khaokhouayensis/rotalis lineage) in the mtDNA analysis. However, the nuDNA data is more likely to group animals from different parts of the Vietnamese mountain areas (Dalat plateau, Gialai–Kontum plateaus, and the northern part of central Vietnam) and from northern Vietnam together, suggesting a common origin of these populations. At the same time, the southern lowland branch has full support and is a sister clade to all other H. gentilis specimens. Given such topology discordance, we may suppose that the H. gentilis populations of southern Vietnam once obtained mtDNA from the common ancestor of the H. khaokhouayensis/rotalis lineage (with which they have no direct contact in present). It is difficult to pinpoint the exact time frame of these introgressive events, but they might have happened approximately on the late Pliocene and early Pleistocene boundary (based on the estimated divergence time between H. gentilis and H. khaokhouayensis: [1]). Noteworthily, the H. gentilis population co-occurring with H. khaokhouayensis in sympatry on Catba island does not show any traces of such introgression. One may suggest other scenarios leading to the discordance between mitochondrial and nuclear data, for example, ancestral polymorphism, but they seem less likely to us than introgression. The verification of one hypothesis or another requires more material than we currently have, and, probably the employment of additional markers (e.g., microsatellites).
Hipposideros gentilis was described from the Thayet (Thayetmyo) district of central Myanmar (Burma) [9,38]. Since it was for long time considered a subspecies of H. pomona, all specimens from Myanmar, northeast India, and Thailand were allocated to this form [4,39]. The occurrence of the form gentilis in Yunnan and northern Laos was also presumed, mainly from the molecular studies [5].
The form sinensis was described from Fujian as a subspecies of H. gentilis [4,40]. According to a published study [5] and our original results, this name could be applied to animals from southeastern China and, at least partially, to this species in north Vietnam (e.g., Gulf of Tonkin). Both lineages are quite close to each other, which also agrees with their morphological similarity. It should be noted, however, that the genetically studied material presented in [5] was collected in northeastern Myanmar, far from the terra typica of H. gentilis. Specimens which belong to the lines SY and SM (sensu Zhao et al. [5]) are morphologically similar to each other, at the same time the morphologically distinct series from the vicinity of Mandalay demonstrates some difference in skull morphometry (Table 3). The type territory of H. gentilis is located even further south (Figure 1). Thus, the interpretation of the SY line as gentilis s. str. can be questioned. Perhaps the entire “northern” clade (from south China, north Vietnam, Laos, and the adjoining parts of Myanmar) should be designated as H. g. sinensis. The relationship between the later and the nominotypical form requires further clarification.
Animals from central Vietnam (SH clade in [5]) could be regarded as a distinct subspecies, though their distribution limits and relationships with the northern mitochondrial lineages require clarification through more extensive material. The only other named form, undoubtedly associated with H. gentilis, is megalotis, which is described from the vicinity of Vinh [11,12]. However, it cannot be used as a valid for this putative subspecies since it represents a junior homonym to the African Hipposideros megalotis (Heuglin, 1861).
We can conclude from our data that Hipposideros gentilis from the lowland forests of southern Vietnam (and probably from the adjacent areas of Cambodia) is quite a distinct form, both genetically and morphologically. It could be, according to nuclear markers, a sister group for all other H. gentilis populations. Meanwhile, no named forms are associated with the “H. pomona” (H. gentilis in current interpretation) populations from the southern part of the species range [41]. This suggests that the “southern” H. gentilis is an undescribed taxon. Its status is not entirely clear, given the common origin and probable past hybridization with other genetic lineages attributed to H. gentilis; and it requires further, more comprehensive, study. Therefore, here we refrain from its formal description.
On the whole, we can conclude that H. gentilis is, contrary to the previous opinion, a monophyletic taxon, presumably divided into two or three subspecies. Its morphological diversity is to some extent in agreement with the obtained genetic lineages, but the decision regarding the taxonomic rank of these lineages can vary greatly depending on the dataset used (mitochondrial or nuclear markers).

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/d13050218/s1, Table S1: Pairwise distances calculated in cytb gene analysis, Figure S1: Phylogenetic relationships of selected Hipposideros species inferred from cytb sequences, Figure S2: Median-joining networks showing the relationships among the alleles of the individual nuclear genes of Hipposideros gentilis and some related forms, Table S2: Correspondence between specimen IDs and haplotypes as they designated on Figure S1, Table S3: Correspondence between nodes of the networks on Figure S1 and additional haplotypes, belonging to the same nodes.

Author Contributions

Conceptualization, A.P.Y., I.V.A. and S.V.K.; methodology, I.V.A. and S.V.K.; validation, A.P.Y., I.V.A. and S.V.K.; formal analysis, A.P.Y., I.V.A. and S.V.K.; investigation, A.P.Y.; data curation, I.V.A. and S.V.K.; writing—original draft preparation, S.V.K.; writing—review and editing, A.P.Y., I.V.A. and S.V.K.; visualization, A.P.Y., I.V.A. and S.V.K.; supervision, S.V.K.; project administration, I.V.A. and S.V.K.; funding acquisition, S.V.K. All authors have read and agreed to the published version of the manuscript.

Funding

Study of collections was supported by the Russian Foundation for Basic Research (grant No. 17-04-00689a). The whole study was performed in line with the stated theme of scientific work of the ZMMU (“Taxonomic and chorological analysis of the animal world, as a ground for study and conservation of the biological diversity”, 121032300105-0).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Genetic data: data available in a publicly accessible repository: GenBank and BOLD (see numbers in Table 2 and Appendix B); morphometric data: data available on request from the corresponding author due to privacy reasons.

Acknowledgments

We would like to express our thanks to all colleagues who provided their priceless help on different stages of our study. Comparative collection materials were studied by SVK in the Royal Ontario Museum, Canada, due to kind permission of J. L. Eger and B. Lim; in the Museum of Natural History, Great Britain, under the support of R. Portela Miguez; in the Berlin Zoological Museum due to permission and help of Frieder Mayer and Nora Lange; in the Hungarian Natural History Museum, Hungary, under the support of G. Csorba; and in the Natural History Museum of Vienna, Austria, under support of Frank Zachos and Alexander Bibl. Study of the collection materials was done in the Zoological Museum of Moscow University, using the collection facilities, with the support of its director, M.V. Kalyakin. Molecular genetic studies were performed on the facilities of the Vertebrate Zoology Department of Moscow University, with invaluable support from A.A. Bannikova. Obtaining materials from Vietnam became possible through collaboration with the Joint Vietnamese-Russian Tropical Research and Technological Centre, due to the support of Nguyen Dang Hoi and A.N. Kuznetsov.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A. List of the GenBank Sequences of the Cytb Mitochondrial Gene, Used in the Analysis

Hipposideros gentilis: Laos: DQ054810.1; China (unspecified): KP876550.1; China, Yunnan: DQ888671.1, KJ623705.1, KP336273.1, KP336274.1; China, Guangdong: EU434950.1, KJ619513.1; China, Hainan: KJ623703.1, KJ623704.1; NE Myanmar: MK064112.1, MK064113.1, MK410336.1, MK410337.1, MK410340.1, MK462234.1; Vietnam, Ha Giang: MK091940.1; Vietnam, Quang Tri: MK091943.1; Vietnam, Tuyen Quang: MK091946.1; Vietnam, Phu Tho: MK091947.1, MK091949.1; Vietnam, Ninh Binh: MK410344.1; Vietnam, Binh Thuan: MK410350.1, MK410351.1, MK430028.1; Vietnam, Vinh Phuc: MK430029.1;
Hipposideros ruber: EU934477.1, EU934485.1, FJ347985.1, FJ347991.1; H. aff. ruber: HQ343240.1, HQ343255.1, HQ343258.1; H. alongensis alongensis: JN247006.1, JN247007.1; H. alongensis sungi: JN247009.1, JN247012.1; H. armiger: JN247034.1; H. armiger terasensis: JN247045.1; H. ater: DQ054807.1; H. bicolor: DQ054808.1, MT149741.1, MT149742.1; H. caffer: FJ347978.1; H. calcaratus: DQ054806.1; H. cervinus: DQ054805.1; H. cf. bicolor: MT149813.1; H. cineraceus: DQ054809.1, KX458067.1KX467584.1, LC406452.1, LC406453.1, LC406454.1, LC406456.1, MK091936.1, MK410352.1; H. diadema: DQ219421.1; H. durgadasi: KY176014.1; H. fuliginosus: EU934467.1, EU934468.1; H. griffini: JN247040.1, JX849199.1; H. khaokhouayensis: DQ054815.1, DQ054816.1; H. larvatus s. lato: DQ888672.1, JN247026.1, JN247027.1; H. lylei: JN247043.1, KR908662.1; H. swinhoei: KJ094477.1, KR908659.1; H. ridleyi: DQ054812.1; H. rotalis: DQ054813.1, DQ054814.1; H. sp. (China): EU434947.1, EU434948.1; H. pendelburyi: JN247028.1; H. turpis: JN247046.1.

Appendix B. List of the BOLD Process ID Numbers for Specimens Used in the Analysis

Hipposideros gentilis: China, Guangxi: ABCMA599-07, ABCMA601-07, ABCMA609-07; China, Guizhou: ABCMA695-07, ABCMA755-07; Laos: ABBM105-05, ABBM131-05, ABBM210-05, ABRLA158-06, BM160-03, ABBM256-05, ABBM277-05, ABRLA034-06, ABRLA087-06, ABRLA088-06, ABRLA130-06, ABRLA163-06, BM087-03, ABBM336-05, ABBM237-05, ABBM360-05, ABBM375-05, ABRLA066-06, BM036-03; Myanmar: ABBSI116-08, ABBM465-05, BM329-03, BM330-03, BM336-03; Vietnam, Ba Ria-Vung Tau: ABBSI346-11, SKMZM1111-12; Vietnam, Dong Nai: ABBSI259-10, ABBSI404-11, BM618-04; Vietnam, Gia Lai: SKZMR126-19, SKZMR127-19, SKZMR128-19; Vietnam, Hai Phong: SKMZM1097-12, SKMZM1098-12; Vietnam, Kien Giang: BM681-05; Vietnam, Lam Dong: SKMZM1199-13; Vietnam, Quang Nam: ABRVN520-06, ABRVN521-06, ABRVN522-06, ABRVN523-06, ABRVN524-06, ABRVN525-06, ABRVN526-06, ABRVN527-06, ABRVN528-06, ABRVN529-06, ABRVN530-06, ABRVN531-06, ABRVN532-06, ABRVN534-06, ABRVN535-06, ABRVN536-06, ABRVN537-06, ABRVN538-06, ABRVN539-06, ABRVN540-06, ABRVN548-06, ABRVN549-06, ABRVN560-06, ABRVN561-06, ABRVN562-06, ABRVN563-06, ABRVN564-06, ABRVN565-06, ABRVN566-06, ABRVN568-06, ABRVN569-06, ABRVN573-06; Vietnam, Quang Ninh: BM659-05; Vietnam, Tuyen Quang: ABRVN141-06, ABRVN175-06, BM328-04; Vietnam, Vinh Phuc: ABRVN039-06, ABRVN040-06, ABRVN041-06, ABRVN042-06, ABRVN043-06, ABRVN044-06, ABRVN045-06, ABRVN046-06;
Hipposideros cf. bicolor: Malaysia: BM453-04, BM452-04; Thailand ABBM049-05, ABBM050-05; H. cineraceus: Vietnam: ABBSI305-11, ABBSI306-11, ABBSI307-11, ABBSI308-11, ABBSI309-11, ABBSI262-10, ABBSI264-10, ABBSI265-10, ABBSI407-11, BM633-04, BM660-05; H. galeritus: Vietnam: ABBSI310-11, ABBSI311-11, ABBSI312-11, ABBSI313-11, ABBSI314-11, SKMZM1126-12, SKMZM1130-12, ABBSI236-10; H. halophyllus: Malaysia: ABBSI019-04; H. khaokhouayensis: Laos: ABBM376-05, ABBM382-05; Vietnam: SKMZM1108-12; H. ridleyi: Malaysia: BM470-04, BM471-04; H. rotalis: Laos: BM059-03, ABRLA167-06, ABRLA169-06, ABRLA170-06.

Appendix C. List of Specimens Used in the Morphometric Analysis. Specimens IDs Genotyped by at Least One Mitochondrial Gene Are in Italic, Specimens Genotyped by Nuclear Genes Are in Bold

India, Haleri: MNH 18.8.3.4 ? (type of H. pomona); India, SW India: MNH 2003.397 m; MNH 2003.398 m; MNH 2003.399 m; India, Darjeeling: MNH 21.1.17.78 m; MNH 21.1.17.79 f; MNH 21.1.17.87 m; Myanmar: ZMB 49464 m; ZMB 49465 m; ZMB 49466 m; ZMB 49467 m; ZMB 49471 m; ZMB 49474 m; ZMB 49480 f; ZMB 49482 m; ZMB 49484 m; ZMB 49486 f; ZMB 49487 f; ZMB 49488 f; ZMB 49489 f; ZMB 49491 f; ZMB 49492 f; ZMB 49496 f; China, Guangxi: ROM MAM 116072 m; ROM MAM 116077 m; China: Guizhou: ROM MAM 118538 f; ROM MAM 118549 f; China, Fujian: MNH 4.12.2.7 ?; China, Hong Kong: NMNS 06651 m; N. Thailand: ZRC 4.6714 ?; ZRC 4.6715 ?; Thailand, Uthai Thani: NMW 65447 m; NMW 65448 m; NMW 65449 m; NMW 65450 m; NMW 65451 f; NMW 65452 f; NMW 65453 m; NMW 65454 m; NMW 65455 m; NMW 65456 m; NMW 65457 m; NMW 65458 m; NMW 65459 m; Cambodia: HNHM 2005.81.11 f; HNHM 2005.81.13 m; Laos, Khammouan: HNHM 2005.82.46 f; Laos, Tham Phakeo: MNHN 2006-0075 m; MNHN 2006-0076 f; MNHN 2006-0077 m; MNHN 2006-0078 m; Vietnam, Lai Chau: MNH 1997.387 ?; Vietnam, Cao Bang: MNH 1997.326 ?; Vietnam, Tuyen Quang: MNH 1997.383 ?; ROM MAM 107660 m; ROM MAM 107700 m; Vietnam, Ha Noi: HNHM 88.24.2 f; Vietnam, Cat Ba Island: HNHM 98.90.3 f; ZMMU S-190298 f; ZMMU S-190299 m; Vietnam, Fuong Vong Islands ZMMU S-144297 f; Vietnam, Vinh Phuoc: ROM MAM 107545 m; ROM MAM 107546 f; ROM MAM 107547 m; ROM MAM 107548 m; Vietnam, Ninh Binh: HNHM 88.26.1 m; Vietnam, Thanh Hoa: HNHM 88.27.2 m; HNHM 88.28.2 m; HNHM 88.28.3 m; HNHM 88.28.5 m; HNHM 88.29.1 m; HNHM 88.29.2 f; HNHM 88.29.3 f; HNHM 88.29.4 f; Vietnam, Nghe An: MNHN 1947-644 m (type of “P.” megalotis); Vietnam, Ha Tinh: ZMMU S-164988 m; Vietnam, Quang Binh: ZMMU S-167172 m; ZMMU S-167174 f; ZMMU S-167175 f; Vietnam, Quang Nam: ROM MAM 111342 m; ROM MAM 111343 f; ROM MAM 111344 f; ROM MAM 111345 m; ROM MAM 111346 m; ROM MAM 111347 m; ROM MAM 111348 m; ROM MAM 111352 f; ROM MAM 111353 f; ROM MAM 111356 m; ROM MAM 111357 m; ROM MAM 111371 f; ROM MAM 111403 m; Vietnam, Gia Lai: ZMMU S-198153 m; ZMMU S-198154 m; ZMMU S-198155 m; Vietnam, Dak Lak: ZMMU S-190724 m; Vietnam, Dong Nai: ZMMU S-175109 f; ZMMU S-181871 f; ZMMU S-186562 m; ZMMU S-194706 m; Vietnam, Phu Quoc Island: ZMMU S-175402 f; ZMMU S-175403 f; Vietnam, Ba Ria-Vung Tau: ZMMU S-188169 m; ZMMU S-190302 m.

References

  1. Foley, N.M.; Goodman, S.M.; Whelan, C.V.; Puechmaille, S.J.; Teeling, E. Towards navigating the Minotaur’s labyrinth: Cryptic diversity and taxonomic revision within the speciose genus Hipposideros (Hipposideridae). Acta Chiropterol. 2017, 19, 1–18. [Google Scholar] [CrossRef]
  2. Monadjem, A.; Soisook, P.; Thong, V.D.; Kingston, T. Family Hipposideridae (Old World Leaf-nosed Bats). In Handbook of the Mammals of the World; Wilson, D.E., Mittermeier, R.A., Eds.; Lynx Edictions: Barcelona, Spain, 2019; Volume 9, pp. 210–259. [Google Scholar]
  3. Douangboubpha, B.; Bumrungsri, S.; Soisook, P.; Murray, S.W.; Puechmaille, S.J.; Satasook, C.; Bu, S.S.H.; Harrison, D.L.; Bates, P.J.J. A taxonomic review of Hipposideros halophyllus, with additional information on H. ater and H. cineraceus (Chiroptera: Hipposideridae) from Thailand and Myanmar. Acta Chiropterol. 2010, 12, 29–50. [Google Scholar] [CrossRef]
  4. Douangboubpha, B.; Bumrungsri, S.; Soisook, P.; Satasook, C.; Thomas, N.M.; Bates, P.J.J. A taxonomic review of the Hipposideros bicolor species complex and H. pomona (Chiroptera: Hipposideridae) in Thailand. Acta Chiropterol. 2010, 12, 415–438. [Google Scholar] [CrossRef]
  5. Zhao, L.-Z.; Bu, Y.-Z.; Zhou, H.-X.; Zhou, H.-W.; Zhang, Z.-X.; Niu, H.-X. Differences in Hipposideros pomona from three geographical regions in China based on morphology and molecular sequences data. J. Mammal. 2015, 96, 1305–1316. [Google Scholar] [CrossRef] [Green Version]
  6. Corbet, G.B.; Hill, J.E. The Mammals of the Indomalayan Region: A Systematic Review; Oxford University Press: Oxford, UK, 1992; pp. 1–488. [Google Scholar]
  7. Francis, C.M. A Field Guide to the Mammals of South-East Asia; New Holland: London, UK, 2008; pp. 1–392. [Google Scholar]
  8. Kruskop, S.V. Bats of Vietnam. Checklist and an Identification Manual, 2nd ed.; KMK Ltd.: Moscow, Russia, 2013; pp. 1–300. [Google Scholar]
  9. Hill, J.E. A revision of the genus Hipposideros. Bull. Brit. Mus. Nat. Hist. Zool. 1963, 11, 1–129. [Google Scholar] [CrossRef]
  10. Hill, J.E.; Zubaid, A.; Davison, G.W.H. The taxonomy of leaf-nosed bats of the Hipposideros bicolor group (Chiroptera: Hipposideridae) from southeastern Asia. Mammalia 1986, 50, 536–540. [Google Scholar] [CrossRef]
  11. Dorst, J. Une nouvelle chauve-souris de l’Indochine française Paracoelops megalotis. Bull. Paris Mus. Nat. Hist., ser. 2 1947, 19, 436–437. [Google Scholar]
  12. Thong, V.D.; Dietz, C.; Denzinger, A.; Bates, P.J.J.; Puechmaille, S.J.; Callou, C.; Schnitzler, H.-U. Resolving a mammal mystery: The identity of Paracoelops megalotis (Chiroptera: Hipposideridae). Zootaxa 2012, 3505, 75–85. [Google Scholar] [CrossRef] [Green Version]
  13. Francis, C.M.; Borisenko, A.V.; Ivanova, N.V.; Eger, J.L.; Lim, B.K.; Guillén-Servent, A.; Kruskop, S.V.; Mackie, I.; Hebert, P.D.N. The Role of DNA Barcodes in Understanding and Conservation of Mammal Diversity in Southeast Asia. PLoS ONE 2010, 5, e12575. [Google Scholar] [CrossRef] [Green Version]
  14. Murray, S.W.; Campell, P.; Kingston, T.; Zubaid, A.; Francis, C.M.; Kunz, T.H. Molecular phylogeny of hipposiderid bats from Southeast Asia and evidence of cryptic diversity. Mol. Phylogenetics Evol. 2012, 62, 597–611. [Google Scholar] [CrossRef]
  15. Kruskop, S.V. The bacula of some bat species from Indo-China: Rhinolophoids (Chiroptera: Rhinolophidae, Hipposideridae). Plecotus et al. 2014, 17, 3–17. [Google Scholar]
  16. Topal, G. Bacula of some Old World leaf-nosed bats (Rhinolophidae and Hipposideridae, Chiroptera: Mammalia). Vertebr. Hung. 1975, 16, 21–53. [Google Scholar]
  17. Shrinivasulu, B.; Shrinivasulu, C. In plain sight: Bacular and noseleaf morphology supports distinct specific status of Roundleaf Bats Hipposideros pomona Andersen, 1918 and Hipposideros gentilis Andersen, 1918 (Chiroptera: Hipposideridae). J. Threat. Taxa 2018, 10, 12018–12026. [Google Scholar] [CrossRef]
  18. Srinivasulu, C.; Srinivasulu, B.; Srinivasulu, A. Hipposideros gentilis. The IUCN Red List of Threatened Species 2020. Available online: https://0-dx-doi-org.brum.beds.ac.uk/10.2305/IUCN.UK.2020-3.RLTS.T180991219A180991293.en (accessed on 25 February 2021).
  19. Sambrook, J.; Fritsch, E.F.; Maniatis, T. Molecular Cloning: A Laboratory Manual, 2nd ed.; Cold Spring Harbor Laboratory Press: New York, NY, USA, 1989; pp. 1–1546. [Google Scholar]
  20. Matthee, C.A.; Burzlaff, J.D.; Taylor, J.F.; Davis, S.K. Mining the mammalian genome for artiodactyl systematics. Syst. Biol. 2001, 50, 1–24. [Google Scholar] [CrossRef]
  21. Artyushin, I.V.; Bannikova, A.A.; Lebedev, V.S.; Kruskop, S.V. Mitochondrial DNA relationships among North Palaearctic Eptesicus (Vespertilionidae, Chiroptera) and past hybridization between Common Serotine and Northern Bat. Zootaxa 2009, 2262, 40–52. [Google Scholar] [CrossRef]
  22. Igea, J.; Juste, J.; Castresana, J. Novel intron markers to study the phylogeny of closely related mammalian species. BMC Evol. Biol. 2010, 10, 369. [Google Scholar] [CrossRef] [Green Version]
  23. Yusefovich, A.P.; Artyushin, I.V.; Raspopova, A.A.; Bannikova, A.A.; Kruskop, S.V. An attempt to reconstruct the phylogeny of the Hipposideros leaf-nosed bats based on nuclear gene markers. Dokl. Biol. Sci. 2020, 493, 136–140. [Google Scholar] [CrossRef]
  24. Ratnasingham, S.; Hebert, P.D.N. BOLD: The Barcode of Life Data System (www.barcodinglife.org). Mol. Ecol. Notes 2007, 7, 355–364. [Google Scholar] [CrossRef] [Green Version]
  25. Katoh, T.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  26. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Datasets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  27. Chernomor, O.; von Haeseler, A.; Minh, B.Q. Terrace aware data structure for phylogenomic inference from supermatrices. Syst. Biol. 2016, 65, 997–1008. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; von Haeseler, A.; Lanfear, R. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Huelsenbeck, J.P.; Ronquist, F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 2001, 17, 754–755. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Ronquist, F.; Huelsenbeck, J.P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19, 1572–1574. [Google Scholar] [CrossRef] [Green Version]
  32. Altekar, G.; Dwarkadas, S.; Huelsenbeck, J.P.; Ronquist, F. Parallel Metropolis coupled Markov chain Monte Carlo for Bayesian phylogenetic inference. Bioinformatics 2004, 20, 407–415. [Google Scholar] [CrossRef] [PubMed]
  33. Lanfear, R.; Frandsen, P.B.; Wright, A.M.; Senfeld, T.; Calcott, B. PartitionFinder 2: New methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol. Biol. Evol. 2017, 34, 772–773. [Google Scholar] [CrossRef] [Green Version]
  34. Baum, B.R. Combining trees as a way of combining datasets for phylogenetic inference, and the desirability of combining gene trees. Taxon 1992, 41, 3–10. [Google Scholar] [CrossRef]
  35. Ragan, M.A. Phylogenetic inference based on matrix representation of trees. Mol. Phyl. Evol. 1992, 1, 53–58. [Google Scholar] [CrossRef]
  36. Ragan, M.A. Matrix representation in reconstructing phylogenetic relationships among the eukaryotes. Bio Syst. 1992, 28, 47–55. [Google Scholar] [CrossRef]
  37. Revell, L.J. Phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol. 2012, 3, 217–223. [Google Scholar] [CrossRef]
  38. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  39. Bates, P.J.J.; Harrison, D.L. Bats of the Indian Subcontinent; Harrison Zoological Museum: Sevenoaks, UK, 1997; pp. 1–258. [Google Scholar]
  40. Andersen, K. Diagnoses of new bats of the families Rhinolophidae and Megadermatidae. Ann. Mag. Nat. Hist. 1918, 2, 374–384. [Google Scholar] [CrossRef] [Green Version]
  41. Murray, S.W.; Khan, F.A.; Kingston, T.; Zubaid, A.; Campbell, P. A new species in the Hipposideros bicolor group (Chiroptera: Hipposideridae) from Peninsular Malaysia. Acta Chiropterol. 2018, 20, 1–29. [Google Scholar] [CrossRef]
Diversity 13 00218 g001 550
Figure 1. Distribution range of Hipposideros gentilis (adopted from [18]). Localities of morphologically studied material are marked with full circles; localities for original genetic samples (cytb and nuclear genes) are marked with empty circles. Type localities of named forms are numbered as follows: 1—gentilis Andersen, 2—sinensis Andersen, 3—megalotis Dorst.
Figure 1. Distribution range of Hipposideros gentilis (adopted from [18]). Localities of morphologically studied material are marked with full circles; localities for original genetic samples (cytb and nuclear genes) are marked with empty circles. Type localities of named forms are numbered as follows: 1—gentilis Andersen, 2—sinensis Andersen, 3—megalotis Dorst.
Diversity 13 00218 g001
Diversity 13 00218 g002 550
Figure 2. Phylogenetic relationships of selected Hipposideros species inferred from cytb (left) and COI (right) sequences. Topologies for both genes are constructed using maximum likelihood (ML). Bayesian posterior probabilities/maximum likelihood bootstrap supports are shown at the tree nodes. Only bootstrap support values over 70 are shown; maximum supports (Bayesian probability 0.99 and higher, 100% ML) are denoted by asterisks. Dashes denote nodes, which didn’t appear on MB tree. Scale bars indicate genetic distances estimated by the ML method. The bootstrap values are derived from 1000 replications.
Figure 2. Phylogenetic relationships of selected Hipposideros species inferred from cytb (left) and COI (right) sequences. Topologies for both genes are constructed using maximum likelihood (ML). Bayesian posterior probabilities/maximum likelihood bootstrap supports are shown at the tree nodes. Only bootstrap support values over 70 are shown; maximum supports (Bayesian probability 0.99 and higher, 100% ML) are denoted by asterisks. Dashes denote nodes, which didn’t appear on MB tree. Scale bars indicate genetic distances estimated by the ML method. The bootstrap values are derived from 1000 replications.
Diversity 13 00218 g002
Diversity 13 00218 g003 550
Figure 3. Phylogenetic tree based on seven nuclear gene topologies combined by MRP method. Bayesian/maximum likelihood/MRP supports are shown at the tree nodes. Maximum supports are denoted with asterisks. Dash symbols denote nodes, which didn’t appear on respective tree.
Figure 3. Phylogenetic tree based on seven nuclear gene topologies combined by MRP method. Bayesian/maximum likelihood/MRP supports are shown at the tree nodes. Maximum supports are denoted with asterisks. Dash symbols denote nodes, which didn’t appear on respective tree.
Diversity 13 00218 g003
Diversity 13 00218 g004 550
Figure 4. Bivariate scatter plot for the first two factors of a Principal Component analysis based on 18 cranial and dental measurements of 99 specimens of Hipposideros pomona s. lato. Specimens marked by closed circles were genotyped with at least one mitochondrial gene.
Figure 4. Bivariate scatter plot for the first two factors of a Principal Component analysis based on 18 cranial and dental measurements of 99 specimens of Hipposideros pomona s. lato. Specimens marked by closed circles were genotyped with at least one mitochondrial gene.
Diversity 13 00218 g004
Diversity 13 00218 g005 550
Figure 5. Bivariate scatter plot for the third and fourth factors of a Principal Component analysis based on 18 cranial and dental measurements of 99 specimens of Hipposideros pomona s. lato. Specimens marked by closed circles were genotyped with at least one mitochondrial gene.
Figure 5. Bivariate scatter plot for the third and fourth factors of a Principal Component analysis based on 18 cranial and dental measurements of 99 specimens of Hipposideros pomona s. lato. Specimens marked by closed circles were genotyped with at least one mitochondrial gene.
Diversity 13 00218 g005
Diversity 13 00218 g006 550
Figure 6. Bivariate scatter plot for the first and second canonical variances of a Discriminant Function analysis based on 18 cranial and dental measurements of six training sets of Hipposideros pomona s. lato. All Indian specimens were included as “undetermined.” Specimens marked by closed diamonds were genotyped with at least one mitochondrial gene.
Figure 6. Bivariate scatter plot for the first and second canonical variances of a Discriminant Function analysis based on 18 cranial and dental measurements of six training sets of Hipposideros pomona s. lato. All Indian specimens were included as “undetermined.” Specimens marked by closed diamonds were genotyped with at least one mitochondrial gene.
Diversity 13 00218 g006
Table
Table 1. List of specimens from ZMMU used in the cytb analysis.
Table 1. List of specimens from ZMMU used in the cytb analysis.
Mus. IDSpeciesLocalityAccession
S-190280A. stoliczkanusN. VietnamCat BaMZ219226
S-167159H. cf. armigerC. VietnamQuan BinhMZ219210
S-195483H. cf. griffiniS. VietnamDong NaiMZ219216
S-186724H. cineraceusS. VietnamCon DaoMZ219220
S-186725H. cineraceusS. VietnamCon DaoMZ219221
S-186730H. cineraceusS. VietnamCon DaoMZ219222
S-191867H. cineraceusS. VietnamDong NaiMZ219218
S-195484H. cineraceusS. VietnamDong NaiMZ219219
S-186567H. diademaS. VietnamBinh PhuocMZ219211
S-191868H. galeritusS. VietnamDong NaiMZ219225
S-167170H. gentilisC. VietnamQuan BinhMZ219215
S-167171H. gentilisC. VietnamQuan BinhMZ219214
S-167172H. gentilisC. VietnamQuan BinhMZ219223
S-167173H. gentilisC. VietnamQuan BinhMZ219224
S-190298H. gentilisN. VietnamCat BaMZ219212
S-190301H. gentilisN. VietnamCat BaMZ219213
S-190302H. gentilisS. VietnamBihn ChauMZ219227
S-191870H. gentilisS. VietnamDong NaiMZ219217
S-198154H. gentilisSC. VietnamGia LaiMZ219228
S-189221H. grandisS. VietnamDong NaiMZ219209
Table
Table 2. List of specimens used in the nuclear genes analysis.
Table 2. List of specimens used in the nuclear genes analysis.
speciesID NoReferenceLocalityTHYABHDROGDIRAG2ACOXCOPSSORBS
H. gentilisS-167174originalC Vietnam MZ219110 MZ219140 MZ219182
H. gentilisS-190298originalN. VietnamMZ219096 MZ219134
H. gentilisS-190299originalN. VietnamMZ219097 MZ219135
H. gentilisS-190301originalN. Vietnam MZ219118MZ219199 MZ219133 MZ219179
H. gentilisS-190302originalS VietnamMZ219098MZ219119MZ219200MZ219170MZ219137MZ219160MZ219180
H. gentilisS-191870originalS VietnamMZ219099MZ219122MZ219203MZ219172MZ219138MZ219159MZ219181
H. gentilisS-191871originalS Vietnam MZ219123MZ219204 MZ219139
H. gentilisS-191909originalS Vietnam MZ219124MZ219205MZ219173MZ219136
H. gentilisS-198253originalSC VietnamMZ219107 MZ219175MZ219152 MZ219191
H. gentilisS-198254originalSC VietnamMZ219108MZ219129MZ219207MZ219176MZ219153 MZ219192
H. gentilisS-198255originalSC VietnamMZ219109 MZ219208MZ219177MZ219154 MZ219193
H. abaeML162-210211-
HIPABA		[1]MaliKP176357KP176214KP176320KP176020
H. armigerS-195483originalS VietnamMZ219094MZ219127MZ219206MZ219174MZ219132MZ219157MZ219190
H. armigerT-171109-1[1]VietnamKP176354KP176210KP176317KP176016
H. centralisS-192897originalEthiopiaMZ219105MZ219125 MZ219149
H. cf. abaeS-189528originalEthiopiaMZ219104MZ219116 MZ219150
H. cf. grandisS-195421originalSC VietnamMZ219093MZ219126 MZ219131MZ219156MZ219189
H. cf. ruberML29-310111
-HIPCAF		[1]MaliKP176358KP176215KP176321KP176021
H. cineraceusS-167179originalC Vietnam MZ219111 MZ219146
H. cineraceusS-186724originalS VietnamMZ219100 MZ219165MZ219143MZ219161MZ219185
H. cineraceusS-186725originalS VietnamMZ219101MZ219113 MZ219166MZ219144 MZ219186
H. cineraceusS-186730originalS VietnamMZ219102MZ219114MZ219195MZ219167MZ219145
H. cineraceusS-191867originalS VietnamMZ219103MZ219120MZ219201MZ219171MZ219147MZ219162MZ219183
H. cineraceusS-195484originalS Vietnam MZ219128 MZ219142 MZ219184
H. diadema10.0011[1]ThailandKY552688KY552683KY552686
H. diademaS-186567originalS. VietnamMZ219095MZ219112MZ219194MZ219164 MZ219158
H. galeritusS-191868originalS. Vietnam MZ219121MZ219202 MZ219148 MZ219187
H. galeritusT-090708-6[1]VietnamKP176355KP176211 KP176017
H. grandisS-189221originalS VietnamMZ219092MZ219115MZ219196MZ219168MZ219130MZ219155MZ219178
H. halophyllusHhal2[1]ThailandKP176359KP176216KP176322KP176022
H. jonesiML155BIS-160211-HIPJON[1]Mali KP176213KP176319KP176019
H. khaokhouayensisS-190294originalN Vietnam MZ219117MZ219198 MZ219141MZ219163
H. khaokhouayensisT-070108-1[1]VietnamKY552735KY552730KY552733
H. larvatusHlar29[1]ThailandKP176360KP176217KP176323KP176023
H. “ater”*T-250608-2[1]VietnamKY552740KY552736
H. “pomona”*T-180809-3[1]VietnamKP176356KP176212KP176318KP176018
A. stoliczkanusS-190280originalN VietnamMZ219106 MZ219197MZ219169MZ219151 MZ219188
* Names of sequences taken from Genbank follow those used by the original authors.
Table
Table 3. Between-group pairwise p-distances for H. gentilis and related forms. Distances are shown below the diagonal, and standard error estimates are shown above the diagonal. Groups correspond to the cytb tree branches on Figure 2.
Table 3. Between-group pairwise p-distances for H. gentilis and related forms. Distances are shown below the diagonal, and standard error estimates are shown above the diagonal. Groups correspond to the cytb tree branches on Figure 2.
H. rotalisH. khaokhouayensisH. cineraceus
C. Vietnam		H. cineraceus
S. Vietnam		H. cineraceus
Malaysia		H. gentilis
S. Vietnam		H. gentilis
S-C. Vietnam		H. gentilis
C. Vietnam		H. gentilis
N. Vietnam		H. gentilis
Yunnan		H. gentilis
Hainan		H. gentilis
Laos
H. rotalis 0.00540.00840.00930.00800.00720.00820.00860.00880.00820.00860.0085
H. kha0.0402 0.00870.00940.00810.00810.00870.00900.00880.00820.00910.0089
H. cin C. Vietnam0.10210.0968 0.00760.00710.00880.00930.00730.00830.00730.00790.0073
H. cin S. Vietnam0.09880.09770.0803 0.00760.00910.00940.00890.00870.00770.00860.0087
H. cin Malaysia0.10920.10680.09650.0989 0.00790.00870.00750.00790.00700.00800.0080
H. gen S. Vietnam0.07530.07710.09530.09910.1028 0.00950.00890.00880.00800.00890.0091
H. gen S-C. Vietnam0.06830.07260.09850.09640.10370.0874 0.00920.00930.00890.00920.0094
H. gen C. Vietnam0.10070.09600.08860.09330.09880.09570.0985 0.00760.00650.00560.0047
H. gen N. Vietnam0.09800.09610.09040.09440.10180.09210.09640.0739 0.00610.00820.0079
H. gen Yunnan0.09940.09490.08960.08740.09810.09330.09340.06630.0583 0.00690.0069
H. gen Hainan0.10390.09890.08910.09110.10330.09770.09510.04170.07810.0769 0.0058
H. gen Laos0.09760.09370.08620.09040.10400.09730.10200.02660.07970.07220.0412
Table
Table 4. Results of PCA applied to specimens of the Hipposideros pomona/gentilis species complex based on 18 cranial measurements: factor loadings, eigenvalues, and percentage of total variance. See text for measurement abbreviations. Factor loadings over 0.7 are highlighted in boldface type.
Table 4. Results of PCA applied to specimens of the Hipposideros pomona/gentilis species complex based on 18 cranial measurements: factor loadings, eigenvalues, and percentage of total variance. See text for measurement abbreviations. Factor loadings over 0.7 are highlighted in boldface type.
PC IPC IIPC IIIPC IV
TL0.38740.57000.23480.5591
CCL0.47560.60680.17810.4944
MW0.01970.85540.00300.2787
BCW0.17660.8479−0.2165−0.0034
OH−0.02870.30390.03920.6229
ZW0.13890.3070−0.33320.6959
POC0.11690.1589−0.77890.0429
RW0.29310.25310.06220.6742
RL0.42990.32680.5776−0.1496
CC0.40590.0960−0.02710.7712
MM0.37450.0619−0.42800.6985
CM0.73340.33660.02060.4562
PM0.71230.1227−0.12820.4650
C0.8001−0.1206−0.07970.0967
NO−0.10330.52490.46120.2852
CM0.62730.22830.02390.3402
MDL0.46170.37840.03900.6815
MDH0.1897−0.06610.02500.6867
Eigenvalue8.08502.070051.48101.1372
% Total44.916811.50038.22766.3176
Cumulative variance44.916856.417164.644770.9622
Table
Table 5. Phenetic (squared Mahalanobis) distances between training sets (above diagonal) and significance of difference between them (p-level, below diagonal) calculated in a forward stepwise Discriminant Function analysis. Cases of insignificant difference re highlighted in boldface type.
Table 5. Phenetic (squared Mahalanobis) distances between training sets (above diagonal) and significance of difference between them (p-level, below diagonal) calculated in a forward stepwise Discriminant Function analysis. Cases of insignificant difference re highlighted in boldface type.
Training Sets123456
1. Myanmar 21.339628.999431.069819.368125.8214
2. C Thailand0.0000 13.544815.634910.288615.3410
3. NW Indochina0.00000.0000 3.19067.093234.5136
4. NE Indochina/SE China0.00000.00000.4175 10.614140.3524
5. C China0.00000.00020.00680.0000 20.9887
6. S Indochina0.00000.00020.00000.00000,0000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yuzefovich, A.P.; Artyushin, I.V.; Kruskop, S.V. Not the Cryptic Species: Diversity of Hipposideros gentilis (Chiroptera: Hipposideridae) in Indochina. Diversity 2021, 13, 218. https://0-doi-org.brum.beds.ac.uk/10.3390/d13050218

AMA Style

Yuzefovich AP, Artyushin IV, Kruskop SV. Not the Cryptic Species: Diversity of Hipposideros gentilis (Chiroptera: Hipposideridae) in Indochina. Diversity. 2021; 13(5):218. https://0-doi-org.brum.beds.ac.uk/10.3390/d13050218

Chicago/Turabian Style

Yuzefovich, Alexander P., Ilya V. Artyushin, and Sergei V. Kruskop. 2021. "Not the Cryptic Species: Diversity of Hipposideros gentilis (Chiroptera: Hipposideridae) in Indochina" Diversity 13, no. 5: 218. https://0-doi-org.brum.beds.ac.uk/10.3390/d13050218

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