| Citation: |
Herman L Mays Jr, Bailey D McKay, Dieter Thomas Tietze, Cheng-Te Yao, Lindsey N Miller, Kathleen N Moreland, Fumin Lei. 2015: A multilocus molecular phylogeny for the avian genus Liocichla (Passeriformes: Leiothrichidae: Liocichla). Avian Research, 6(1): 17. DOI: 10.1186/s40657-015-0025-y
|
Historically the babblers have been assigned to the family Timaliidae but several recent studies have attempted to rest the taxonomy of this diverse passerine assemblage on a more firm evolutionary footing. The result has been a major rearrangement of the group. A well-supported and comprehensive phylogeny for this widespread avian group is an important part of testing evolutionary and biogeographic hypotheses, especially in Asia where the babblers are a key component of many forest ecosystems. However, the genus Liocichla is poorly represented in these prior studies of babbler systematics.
We used a multilocus molecular genetic approach to generate a phylogenetic hypothesis for all five currently recognized species in the avian genus Liocichla. Multilocus DNA sequence data was used to construct individual gene trees using maximum likelihood and species trees were estimated from gene trees using Bayesian analyses. Divergence dates were obtained using a molecular clock approach.
Molecular data estimate a probable window of time for the origin for the Liocichla from the mid to late Miocene, between 5.55 and 12.87 Ma. Despite plumage similarities between the insular Taiwan endemic, L. steerii, and the continental L. bugunorum and L. omeiensis, molecular data suggest that L. steerii is the sister taxon to all continental Liocichla. The continental Liocichla are comprised of two lineages; a lineage containing L. omeiensis and L. bugunorum and a lineage comprised of L. phoenicea and L. ripponi. The comparatively early divergence of L. steerii within the Liocichla may be illusory due to extinct and therefore unsampled lineages. L. ripponi and L. phoenicea are parapatric with a Pleistocene split (0.07-1.88 Ma) occurring between an Eastern Himalayan L. phoenicea and a Northern Indochina distributed L. ripponi. L. bugunorum and L. omeiensis underwent a similar split between the Eastern Himalaya (L. bugunorum) and Central China (L. omeiensis) divided by the Hengduan Mountains.
This study supports an origin of the Liocichla occurring sometime prior to the Miocene-Pliocene boundary, a period of significant climatic upheaval in Asia. The biogeographical patterns within the Liocichla mirror those of other birds in the region and allude to common geological and climatic drivers of avian diversification in Asia.
The genus Glaucidium consists mainly of small owl species, which are also known as owlets or pygmy owls. Owing to inter-specific plumage similarities, the taxonomy of this genus is confused, with a variation of 26 to 35 species recognised across different taxonomic treatments. Presently, approximately six species of Glaucidium owlets are recognised across Asia (Ritschard and Schweizer 2007; Dickinson and Remsen 2013; del Hoyo and Collar 2014). Among these, the taxonomic status of the Collared Owlet (G. brodiei) species complex is one of the most debatable among the Old World Glaucidium owlets, with a discrepancy in treatments ranging from one (Dickinson and Remsen 2013; del Hoyo and Collar 2014; Gill and Donsker 2017; Clements et al. 2018) to two (Eaton et al. 2016) and possibly three species (Ritschard and Schweizer 2007). The species complex occurs in mostly montane and submontane forest (G. b. brodiei also occurs down to lowlands), and is widespread, with the taxon brodiei occurring from Afghanistan through the Himalayas to Central China and Southeast Asia, sylvaticum from Sumatra, borneense from Borneo and pardalotum from Taiwan Island (Fig. 1). Two other races, tubiger and garoense usually synonymised under the nominate subspecies, although borneense and sylvaticum have also been synonymised by König and Weick (2008).
Discrepancies in the taxonomy of the G. brodiei species complex can be largely attributed to the variation of plumage colours observed across individuals unrelated to their distribution (König and Weick 2008). Although these owlets were traditionally thought to exhibit true colour polymorphism (with grey and rufous morphs recognised), a study by Lin et al. (2014) suggests that this intra-specific variation in plumage colours may be age-dependent. At any rate, plumage comparison may be an unreliable species delineation tool for this group of owls and other species delimitation methods such as bioacoustics and genetics should be employed to accurately determine species boundaries within the G. brodiei complex. Furthermore, several studies have shown the pitfalls of relying solely on morphology in owl taxonomy, and the importance of bioacoustics in species delimitation of owls is well documented (König 1994; King 2002; Gwee et al. 2017).
Although regarded as a single species, vocal differences within the G. brodiei species complex have been noted by field observers. For example, Eaton et al. (2016) split the Sumatran and Bornean populations as a separate species from other members of the G. brodiei species complex based on notably different vocalisations, though quantitative analysis was unavailable. The pattern of species distribution proposed by Eaton et al. (2016) is quite unusual as both islands have been repeatedly connected to the geographically intervening Malay Peninsula throughout the Pleistocene, suggesting the incidence of a leapfrog pattern (Remsen 1984) in which the terminal taxa sylvaticum and borneense are divided from peninsular brodiei. In the present study, we investigated potential species boundaries among each member of the G. brodiei species complex using bioacoustics as well as plumage comparison of museum specimens. These character suites allowed us to examine species delimitation in the complex in the absence of modern DNA material for the insular populations.
We collected a total of 76 sound recordings of the G. brodiei species complex from the online sound library xeno-canto (https://www.xeno-canto.org) and from our personal collection (see Additional file 1: Table S1). In order to avoid duplicate recordings of the same individual, only one out of all the samples recorded by the same person at the same site and time was used for the study. Each sound recording was measured using the default settings on Raven Pro 1.5 (Bioacoustics Research Program, Cornell Laboratory of Ornithology, Ithaca, NY, USA). We measured a total of 8 vocal parameters: (1) number of elements per motif, (2) duration of a motif, (3) lowest frequency in a motif, (4) highest frequency in a motif, (5) bandwidth (= highest minus lowest frequency), (6) duration of first break within a motif, (7) duration of second break within a motif, and (8) duration of breaks between motifs. A motif is defined as a complete song that a bird usually repeats several times, while an element is defined as an individual note in the song.
Rstudio version 1.1.453 (https://www.rstudio.com) and R version 3.5.0 (R Core Team 2018) were used to conduct principal component analysis (PCA) on the vocal dataset to distinguish clinal bioacoustic variation from discrete variation. Pairwise comparison of each vocal parameter between two taxa were done using the criterion outlined by Isler et al. (1998), henceforth referred to as the Isler criterion. The Isler criterion is based on two conditions: (1) there must be no overlap between the ranges of measurements between the two taxa being compared, and (2) the means (x) and standard deviations (SD) of the taxon (t) with the smaller set of measurements (a) and the taxon with the larger set of measurements (b) have to meet the following requirement: xa + taSDa ≤ xb − tbSDb, where ti refers to the one-tailed t-score at the 97.5th percentile of the t distribution for n − 1 degrees of freedom. Although Isler et al. (1998) first applied this method for the species delimitation of suboscine antbirds, this criterion has also been employed across non-oscines such as pigeons (Rheindt et al. 2011; Ng et al. 2016; Ng and Rheindt, 2016), nightjars (Sangster and Rozendaal 2004), owls (Gwee et al. 2017), as well as oscines (Cros and Rheindt 2017; Prawiradilaga et al. 2017; Gwee et al. 2019).
A total of 13 adult specimens, including five mainland brodiei, six Bornean borneense, one Sumatran sylvaticum and one Taiwan pardalotum, were inspected at the Natural History Museum at Tring, UK (NHM Tring) and the Lee Kong Chian Natural History Museum, Singapore (LKCNHM). Photographic evidence of each specimen was taken and relevant plumage traits, such as the colouration of underparts, upperparts and neck collar, were visually assessed across specimens. Additionally, wing and tail measurements of the six specimens from the LKCNHM were obtained by C.Y. Gwee. Given the small sample size, biometric analysis was restricted to a mere comparison of ranges, with no significance testing.
Spectrograms of song examples from across the G. brodiei species complex reflect the stark vocal differences between populations from Borneo and Sumatra versus populations from mainland Asia and Taiwan (Fig. 2). The vocalisation of Bornean and Sumatran populations comprises seven notes per motif, in contrast to the four notes per motif vocalisation of mainland and Taiwan populations (Fig. 2).
PCA across eight vocal traits confirmed that Bornean borneense and Sumatran sylvaticum populations together form a cluster distinct from continental brodiei and Taiwan pardalotum populations (Fig. 3). This division was further corroborated by the Isler criterion results which showed at least two diagnosable vocal parameters between the Sundaic insular cluster and the other two taxa (Tables 1, 2). On the other hand, Bornean and Sumatran populations were vocally indistinguishable from each other (Fig. 3, Tables 1, 2). The mainland Asian and Taiwan populations were vocally indistinguishable, with large spatial overlap on PCA and an absence of any vocal parameter passing the threshold of diagnosability under the Isler criterions (Fig. 3, Table 2).
| Taxon | a | b | c | d | e | f | g | h |
| sylvaticum | 7.00 ± 0.00 | 3.05 ± 0.117 | 618 ± 39.6 | 959 ± 48.1 | 342 ± 60.1 | 0.382 ± 0.0481 | 0.115 ± 0.0158 | 17.8 ± 11.8 |
| borneense | 7.00 ± 0.00 | 2.83 ± 0.0771 | 614 ± 26.2 | 988 ± 76.8 | 374 ± 61.8 | 0.326 ± 0.0802 | 0.149 ± 0.0432 | 17.8 ± 13.8 |
| brodiei | 3.98 ± 0.139 | 1.41 ± 0.180 | 806 ± 53.0 | 1210 ± 56.7 | 400 ± 63.9 | 0.415 ± 0.0553 | 0.0953 ± 0.143 | 0.901 ± 0.223 |
| pardalotum | 4.00 ± 0.00 | 1.89 ± 0.00416 | 879 ± 42.3 | 1170 ± 97.6 | 290 ± 140 | 0.601 ± 0.0185 | 0.102 ± 0.0259 | 1.95 ± 0.116 |
| The following parameters were assessed: (a) number of elements per motif, (b) duration of a motif, (c) lowest frequency, (d) highest frequency, (e) bandwidth, (f) first break length within a motif, (g) second break length within a motif, and (h) break length between motifs | ||||||||
DownLoad:
CSV
| a | b | c | d | e | f | g | h | |
| sylvaticum vs borneense | ||||||||
| sylvaticum vs brodiei | X | X | X | |||||
| sylvaticum vs pardalotum | X | X | X | |||||
| borneense vs brodiei | X | X | X | |||||
| borneense vs pardalotum | X | X | ||||||
| brodiei vs pardalotum | ||||||||
| The following parameters were assessed: (a) number of elements per motif, (b) duration of a motif, (c) lowest frequency, (d) highest frequency, (e) bandwidth, (f) first break length within a motif, (g) second break length within a motif, and (h) break length between motifs | ||||||||
DownLoad:
CSV
Plumage inspection of the 13 specimens from the NHM Tring and the LKCNHM revealed continental brodiei and Taiwan pardalotum share a rufous neck collar, whereas Sumatran sylvaticum and Bornean borneense share a while neck collar (Fig. 4, Additional file 2: Fig. S1). Furthermore, biometric comparisons of specimens from the LKCNHM suggest brodiei has a slightly shorter wing (85 to 91 mm) than sylvaticum and borneense, though the latter two taxa were represented by only one specimen each (Table 3). On the other hand, the tail measurement of brodiei (44 to 57 mm) overlapped with that of borneense (52 mm). The tail length of the sylvaticum specimen was significantly shorter (29 mm), possibly due to tail moult. Given small sample sizes, no statistic comparison was attempted on the basis of our biometric measurements.
| Taxon | Voucher number | Sex | Locality | Wing measurement (mm) | Tail measurement (mm) |
| sylvaticum | ZRC3.8604 | Female | Northeast Sumatra, Indonesia | 105 | 29 |
| borneense | ZRC3.8603 | Male | Sarawak, Malaysia | 107 | 52 |
| borneense | 1900.2.14.11 | Male | Sarawak, Malaysia | – | – |
| borneense | 95.11.15.57 | Male | Sabah, Malaysia | – | – |
| borneense | 98.11.24.62 | Female | Sabah, Malaysia | – | – |
| borneense | 93.6.22.3 | Male | Sabah, Malaysia | – | – |
| borneense | 93.6.21.1 | Male | Sarawak, Malaysia | – | – |
| brodiei | ZRC3.8606 | Female | Chiang Mai, Thailand | 91 | 57 |
| brodiei | ZRC3.8607 | Male | Chiang Mai, Thailand | 88 | 50 |
| brodiei | ZRC3.8610 | Male | Pahang, Malaysia | 85 | 44 |
| brodiei | ZRC3.8612 | Female | Pahang, Malaysia | 91 | 53 |
| brodiei | 86.2.1.683 | Male | Sikkim, India | – | – |
| pardalotum | 1907.12.12.115 | Male | Mt Ho Ho, Taiwan, China | – | – |
DownLoad:
CSV
The innate vocalisations of owls are an important taxonomic indicator of species limits as they are inherited and subject to strong pressure for mate selection in nocturnal birds (King 2002; Gwee et al. 2017). Our bioacoustic data reflected a stark contrast between the well-known homogenous four-note song separated only by short pauses from other motifs of the Collared Owlet throughout mainland Asia and the seven-note song widely separated from other motifs uttered by insular populations in Sundaland (Fig. 2). Similarly, plumage differences in these owlets, particularly the neck collar colouration, seem to mirror the vocal pattern (Fig. 4). Although Sharpe (1875) commented sylvaticum differs from brodiei by its "rufous brown back", we failed to observe this in specimen ZRC3.8604 (Fig. 4c). Given potential age-related plumage variation in Glaucidium owlets, we caution that such plumage comparisons should be conservatively regarded (Ritschard and Schweizer 2007; Lin et al. 2014). Our biometric results seem to suggest sylvaticum and borneense have a greater wing length than brodiei regardless of their sex (Table 1). However, the measurements in Sharpe's (1875) record show the wing length of sylvaticum (3.8 in.) and borneense (3.65 in.) overlapped with the female brodiei (3.6 to 3.8 in.) which are greater than the male brodiei specimens (3.2 to 3.4 in.). Nevertheless, the distinct vocalisations and consistent neck collar colouration pattern of the Sunda Owlet support an elevation of G. sylvaticum to species status, with the junior Bornean taxon G. s. borneense to be reclassified as a subspecies of G. sylvaticum based on the Principle of Priority (ICZN 1999: Art. 23.1).
The Isthmus of Kra is a prominent biogeographical and avifaunal divide between northern Southeast Asian monsoon forests and equatorial Sundaic rainforests (Hughes et al. 2003). Numerous Oriental bird species pairs are characterised by a more monsoon-adapted northern and rainforest-adapted southern species that abut somewhere around this transition zone, including owl pairs such as Sunda Scops Owl (Otus lempiji) versus Collared Scops Owl (O. lettia), and Spot-bellied Eagle Owl (Bubo nipalensis) versus Barred Eagle Owl (B. sumatranus). However, in the case of the Collared Owlet complex, the Sundaic population from the Thai-Malay Peninsula is vocally and morphologically undifferentiated from the more northerly monsoon populations, generating an unusual division between continental (Peninsula Malaysia and other parts of mainland Asia) versus archipelagic populations (Sumatra and Borneo). Despite regular linkage through Quaternary land bridges via the geographically intervening Malay Peninsula (Bintanja et al. 2005), the islands of Sumatra and Borneo are known to share multiple montane and submontane species which are absent on the continental mainland. This pattern is observed in the Collared Owlet complex, as well as other avian species such as Rajah's Scops Owl (Otus brookii) and Black-capped White-eye (Zosterops atricapilla).
While our bioacoustic data did not support distinct divergence between Taiwan pardalotum and continental brodiei (Tables 1, 2), we note that one of the three sound recordings of pardalotum was quite different from continental brodiei (Fig. 3). More samples are required to ascertain whether this vocal difference is clinal as Taiwan forms a continuous landmass with the mainland during the Pleistocene glaciation, or that sample was an anomalous recording. Although the pardalotum specimen seems to have a darker brown upperpart than the brodiei specimens (Fig. 4), intra-specific plumage variation may be present (Ritschard and Schweizer 2007; Lin et al. 2014). Furthermore, Sharpe (1893) found the two taxa to be indistinguishable in plumage during specimen inspection. Therefore, we propose that Taiwan pardalotum be retained as a subspecies of G. brodiei.
Our study leads to a taxonomic division of an important and well-known Asian owl species complex into two resultant species on the basis of bioacoustic and morphological data. It is therefore one of numerous contributions that highlight the importance of bioacoustics as a tool for species delineation in Strigiformes, which are often challenging to identify by plumage due to age-related or ecomorphological variation (König and Weick 2008; Lin et al. 2014; Mikkola 2014; Sadanandan et al. 2015; Gwee et al. 2017). Modern genetic materials of owl species can be difficult to obtain due to their elusive behaviours and ancient genetic materials can be challenging to work with, thus bioacoustics serve as a reliable tool to scan for cryptic diversity in the absence of genetic data. Future studies can look into using molecular tools to investigate the level of genetic divergence within the G. brodiei species complex, as well as playback experiments to assess species recognition between the Sunda Owlet and Collared Owlet (Freeman and Montgomery 2017). In conclusion, we found vocal evidence further supported by plumage comparisons differentiating the insular Sumatran and Bornean taxa from the mainland and Taiwan taxa. We thereby propose the elevation of G. sylvaticum to species status under the common name Sunda Owlet.
Supplementary information accompanies this paper at https://doi.org/10.1186/s40657-019-0175-4.
Additional file 1: Table S1. Information of all recordings used in the study and measurements of all vocal parameters.
Additional file 2: Fig. S1. Photos of all borneense specimens inspected at the Natural History Museum at Tring, UK.
We are indebted to all sound recordists who deposited recordings on Xeno-Canto (https://www.xeno-canto.org), along with Roger McNeill and Oscar Johnson, both of whom deposited their sound recordings with Macaulay Library (https://www.macaulaylibrary.org) and were kind enough to share their recordings with us.
CYG analysed the vocal data and was the main contributor in writing the manuscript. JAE acquired and measured the sound recordings. EYXN interpreted the results. FER designed the research. All authors read and approved the final manuscript.
All data generated or analysed during this study are included in this published article [and its additional information files].
Not applicable.
Not applicable.
The authors declare that they have no competing interests.
|
Collar NJ, Robson C (2007) Family Timaliidae (Babblers). In: del Hoyo J, Elliott A, Christie D (eds) Handbook of birds of the World. Lynx Edicions, Barcelona, pp 70-291
|
|
Pfeifer R (2013) What is a babbler? Fascinating diversity of life strategies and challenge for systematists. Vogelwarte 51:117-126
|
|
Collar NJ (2006) A partial revision of the Asian babblers (Timaliidae). Forktail 22:85-112
|
|
Delacour J (1946) Les timaliinés. L'Oiseau 16:7-36
|
|
Collar NJ (2011) Taxonomic notes on some Asian babblers (Timaliidae). Forktail 27:100-102
|
|
Dickinson EC, Christidis L (2014) The Howard and Moore complete checklist of the birds of the world, vol 2: Passerines, 4th edn. Christopher Helm, London
|
|
Athreya R (2006) A new species of Liocichla (Aves: Timaliidae) from Eaglenest Wildlife Sanctuary, Arunachal Pradesh, India. Indian Birds 2:82-94
|
|
Peterson AT, Papeş M (2006) Potential geographic distribution of the Bugun Liocichla Liocichla bugunorum, a poorly-known species from north-eastern India. Indian Birds 2:146-149
|
|
Dickinson EC (2003) The Howard and Moore complete checklist of the birds of the world, 3rd edn. Christopher Helm, London
|
|
McKay BD (2011) A new timeframe for the diversification of Japan's mammals. J Biogeogr 39:1134-1143
|
|
Päckert M, Martens J, Severinghaus LL (2008) The Taiwan Firecrest (Regulus goodfellowi) belongs to the Goldcrest assemblage (Regulus regulus s. l.): evidence from mitochondrial DNA and the territorial song of the Regulidae. J Ornithol 150:205-220
|
| 1. | Sreehari K. Mohan, R. Anjitha, K. Maxim Rodrigues. The avian diversity of Chemmattamvayal Wetlands and adjacent areas of Kasaragod District, Kerala, India. Journal of Threatened Taxa, 2023, 15(4): 23047. DOI:10.11609/jott.5246.15.4.23047-23060 |
| 2. | Shuchan Zhou, Kan Wang, Beata Messyasz, et al. Functional and taxonomic beta diversity of butterfly assemblages in an archipelago: relative importance of island characteristics, climate, and spatial factors. Ecological Indicators, 2022, 142: 109191. DOI:10.1016/j.ecolind.2022.109191 |
| 3. | Asha Sohil, Neeraj Sharma. Assessing the bird guild patterns in heterogeneous land use types around Jammu, Jammu and Kashmir, India. Ecological Processes, 2020, 9(1) DOI:10.1186/s13717-020-00250-9 |
| 4. | Zhifeng Ding, Jianchao Liang, Yiming Hu, et al. Different responses of avian feeding guilds to spatial and environmental factors across an elevation gradient in the central Himalaya. Ecology and Evolution, 2019, 9(7): 4116. DOI:10.1002/ece3.5040 |
| 5. | Xingfeng Si, Marc W. Cadotte, Di Zeng, et al. Functional and phylogenetic structure of island bird communities. Journal of Animal Ecology, 2017, 86(3): 532. DOI:10.1111/1365-2656.12650 |
| 6. | Xingfeng Si, Andrés Baselga, Fabien Leprieur, et al. Selective extinction drives taxonomic and functional alpha and beta diversities in island bird assemblages. Journal of Animal Ecology, 2016, 85(2): 409. DOI:10.1111/1365-2656.12478 |
Figures(5) / Tables(3)
| Taxon | a | b | c | d | e | f | g | h |
| sylvaticum | 7.00 ± 0.00 | 3.05 ± 0.117 | 618 ± 39.6 | 959 ± 48.1 | 342 ± 60.1 | 0.382 ± 0.0481 | 0.115 ± 0.0158 | 17.8 ± 11.8 |
| borneense | 7.00 ± 0.00 | 2.83 ± 0.0771 | 614 ± 26.2 | 988 ± 76.8 | 374 ± 61.8 | 0.326 ± 0.0802 | 0.149 ± 0.0432 | 17.8 ± 13.8 |
| brodiei | 3.98 ± 0.139 | 1.41 ± 0.180 | 806 ± 53.0 | 1210 ± 56.7 | 400 ± 63.9 | 0.415 ± 0.0553 | 0.0953 ± 0.143 | 0.901 ± 0.223 |
| pardalotum | 4.00 ± 0.00 | 1.89 ± 0.00416 | 879 ± 42.3 | 1170 ± 97.6 | 290 ± 140 | 0.601 ± 0.0185 | 0.102 ± 0.0259 | 1.95 ± 0.116 |
| The following parameters were assessed: (a) number of elements per motif, (b) duration of a motif, (c) lowest frequency, (d) highest frequency, (e) bandwidth, (f) first break length within a motif, (g) second break length within a motif, and (h) break length between motifs | ||||||||
DownLoad:
CSV
| a | b | c | d | e | f | g | h | |
| sylvaticum vs borneense | ||||||||
| sylvaticum vs brodiei | X | X | X | |||||
| sylvaticum vs pardalotum | X | X | X | |||||
| borneense vs brodiei | X | X | X | |||||
| borneense vs pardalotum | X | X | ||||||
| brodiei vs pardalotum | ||||||||
| The following parameters were assessed: (a) number of elements per motif, (b) duration of a motif, (c) lowest frequency, (d) highest frequency, (e) bandwidth, (f) first break length within a motif, (g) second break length within a motif, and (h) break length between motifs | ||||||||
DownLoad:
CSV
| Taxon | Voucher number | Sex | Locality | Wing measurement (mm) | Tail measurement (mm) |
| sylvaticum | ZRC3.8604 | Female | Northeast Sumatra, Indonesia | 105 | 29 |
| borneense | ZRC3.8603 | Male | Sarawak, Malaysia | 107 | 52 |
| borneense | 1900.2.14.11 | Male | Sarawak, Malaysia | – | – |
| borneense | 95.11.15.57 | Male | Sabah, Malaysia | – | – |
| borneense | 98.11.24.62 | Female | Sabah, Malaysia | – | – |
| borneense | 93.6.22.3 | Male | Sabah, Malaysia | – | – |
| borneense | 93.6.21.1 | Male | Sarawak, Malaysia | – | – |
| brodiei | ZRC3.8606 | Female | Chiang Mai, Thailand | 91 | 57 |
| brodiei | ZRC3.8607 | Male | Chiang Mai, Thailand | 88 | 50 |
| brodiei | ZRC3.8610 | Male | Pahang, Malaysia | 85 | 44 |
| brodiei | ZRC3.8612 | Female | Pahang, Malaysia | 91 | 53 |
| brodiei | 86.2.1.683 | Male | Sikkim, India | – | – |
| pardalotum | 1907.12.12.115 | Male | Mt Ho Ho, Taiwan, China | – | – |
DownLoad:
CSV
Email Alerts
RSS Feeds