Geographical variation in and evolutionary history of the Sunda clouded leopard (Neofelis diardi) (Mammalia: Carnivora: Felidae) with the description of a new subspecies from Borneo

First mitogenome phylogeny of the sun bear Helarctos malayanus reveals a deep split between Indochinese and Sundaic lineages

The sun bear Helarctos malayanus is one of the most endangered ursids, and to date classification of sun bear populations has been based almost exclusively on geographic distribution and morphology. The very few molecular studies focussing on this species were limited in geographic scope. Using archival and non-invasively collected sample material, we have added a substantial number of complete or near-complete mitochondrial genome sequences from sun bears of several range countries of the species' distribution. We here report 32 new mitogenome sequences representing sun bears from Cambodia, Thailand, Peninsular Malaysia, Sumatra, and Borneo. Reconstruction of phylogenetic relationships revealed two matrilines that diverged ~295 thousand years ago: one restricted to portions of mainland Indochina (China, Cambodia, Thailand; "Mainland clade"), and one comprising bears from Borneo, Sumatra, Peninsular Malaysia but also Thailand ("Sunda clade"). Generally recent coalescence times in the mitochondrial phylogeny suggest that recent or historical demographic processes have resulted in a loss of mtDNA variation. Additionally, analysis of our data in conjunction with shorter mtDNA sequences revealed that the Bornean sun bear, classified as a distinct subspecies (H. m. euryspilus), does not harbor a distinctive matriline. Further molecular studies of H. malayanus are needed, which should ideally include data from nuclear loci.

Whole genome analysis of clouded leopard species reveals an ancient divergence and distinct demographic histories

Similar to other apex predator species, populations of mainland (Neofelis nebulosa) and Sunda (Neofelis diardi) clouded leopards are declining. Understanding their patterns of genetic variation can provide critical insights on past genetic erosion and a baseline for understanding their long-term conservation needs. As a step toward this goal, we present draft genome assemblies for the two clouded leopard species to quantify their phylogenetic divergence, genome-wide diversity, and historical population trends. We estimate that the two species diverged 5.1 Mya, much earlier than previous estimates of 1.41 Mya and 2.86 Mya, suggesting they separated when Sundaland was becoming increasingly isolated from mainland Southeast Asia. The Sunda clouded leopard displays a distinct and reduced effective population size trajectory, consistent with a lower genome-wide heterozygosity and SNP density, relative to the mainland clouded leopard. Our results provide new insights into the evolutionary history and genetic health of this unique lineage of felids.

Taxonomic boundaries in Lesser Treeshrews (Scandentia, Tupaiidae: Tupaia minor )

The Lesser Treeshrew, Tupaia minor  Günther, 1876, is a small mammal from Southeast Asia with four currently recognized subspecies: T. m. minor from Borneo; T. m. malaccana from the Malay Peninsula; T. m. humeralis from Sumatra; and T. m. sincepis from Singkep Island and Lingga Island. A fifth subspecies, T. m. caedis, was previously synonymized with T. m. minor; it was thought to occur in northern Borneo and on the nearby islands of Banggi and Balambangan. These subspecies were originally differentiated based on pelage color, a plastic feature that has proven to be an unreliable indicator of taxonomic boundaries in treeshrews and other mammals. To explore infraspecific variation among T. minor populations across the Malay Peninsula, Borneo, Sumatra, and smaller islands, we conducted multivariate analyses of morphometric data collected from the hands and skulls of museum specimens. Principal component and discriminant function analyses reveal limited differentiation in manus and skull proportions among populations of T. minor from different islands. We find no morphometric support for the recognition of the four allopatric subspecies and no support for the recognition of T. m. caedis as a separate subspecies on Borneo.

The genome of Shorea leprosula (Dipterocarpaceae) highlights the ecological relevance of drought in aseasonal tropical rainforests

Hyperdiverse tropical rainforests, such as the aseasonal forests in Southeast Asia, are supported by high annual rainfall. Its canopy is dominated by the species-rich tree family of Dipterocarpaceae (Asian dipterocarps), which has both ecological (e.g., supports flora and fauna) and economical (e.g., timber production) importance. Recent ecological studies suggested that rare irregular drought events may be an environmental stress and signal for the tropical trees. We assembled the genome of a widespread but near threatened dipterocarp, Shorea leprosula, and analyzed the transcriptome sequences of ten dipterocarp species representing seven genera. Comparative genomic and molecular dating analyses suggested a whole-genome duplication close to the Cretaceous-Paleogene extinction event followed by the diversification of major dipterocarp lineages (i.e. Dipterocarpoideae). Interestingly, the retained duplicated genes were enriched for genes upregulated by no-irrigation treatment. These findings provide molecular support for the relevance of drought for tropical trees despite the lack of an annual dry season.

Ecogeographic variation and taxonomic boundaries in Large Treeshrews (Scandentia, Tupaiidae: Tupaia tana Raffles, 1821) from Southeast Asia

The Large Treeshrew, Tupaia tana Raffles, 1821, is a small mammal (~205 g) from Southeast Asia with a complicated taxonomic history. Currently, 15 subspecies are recognized from Borneo, Sumatra, and smaller islands, and many were originally differentiated based on minor pelage differences and small sample sizes. We explored intraspecific variation in T. tana using quantitative osteological data obtained from the hands and skulls of museum specimens. Multivariate analyses reveal extensive overlap among T. tana populations in morphospace, indicating that the majority of currently recognized subspecies are not morphometrically distinct. In contrast, the separation between Bornean and Sumatran populations of T. tana is sufficient to recognize them as different subspecies. Comparisons of Bornean specimens to those on small, offshore islands reveal that the latter average smaller body size. This pattern is inconsistent with Foster’s island rule, which predicts that island populations of small mammals (< 5 kg) will average larger body size relative to mainland forms. A similar lack of support for ecogeographic rules has been noted in T. glis (Diard, 1820), suggesting that these “rules” are poor predictors of geographic variation in treeshrews.

Mitogenome Phylogeny Including Data from Additional Subspecies Provides New Insights into the Historical Biogeography of the Eurasian lynx Lynx lynx

Previous molecular studies of the wide-ranging Eurasian lynx Lynx lynx focused mainly on its northern Palearctic populations, with the consequence that the reconstruction of this species’ evolutionary history did not include genetic variation present in its southern Palearctic distribution. We sampled a previously not considered Asian subspecies (L. l. dinniki), added published data from another Asian subspecies (L. l. isabellinus), and reassessed the Eurasian lynx mtDNA phylogeny along with previously published data from northern Palearctic populations. Our mitogenome-based analyses revealed the existence of three major clades (A: Central Asia, B: SE Europe/SW Asia, C: Europe and Northern Asia) and at least five lineages, with diversification in Lynx lynx commencing at least 28kyr earlier than hitherto estimated. The subspecies L. l. isabellinus harbors the most basal matriline, consistent with the origin of Lynx lynx in this subspecies’ current range. L. l. dinniki harbors the second most basal matriline, which is related to, and may be the source of, the mtDNA diversity of the critically endangered Balkan lynx L. l. balcanicus. Our results suggest that the Anatolian peninsula was a glacial refugium for Eurasian lynx, with previously unconsidered implications for the colonization of Europe by this species.

The illegal exploitation of the Javan Leopard (Panthera pardus melas) and Sunda Clouded Leopard (Neofelis diardi) in Indonesia

Indonesia is home to the Javan Leopard (Panthera pardus melas) and the Sunda Clouded Leopard (Neofelis diardi), both of which are threatened by habitat loss, human-wildlife conflict issues and the illegal wildlife trade. Leopards and clouded leopards are threatened by the illegal wildlife trade across their range, however, very little is known of the illegal trade in these two species in Indonesia, or of the efforts made to tackle this crime. Both the Javan Leopard and Sunda Clouded Leopard are protected species in Indonesia and both species are listed in Appendix I of the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES), meaning commercial international trade is generally prohibited. To better understand the trade, and efforts to end this trade, we collected records of seizures and prosecutions relating to Javan Leopards and Sunda Clouded Leopards in Indonesia for the period 2011–2019. Despite both species being protected by law, this study reveals the prevalence of extensive poaching, illegal domestic trade and international trafficking of both species. A total of 41 seizure records were obtained from 2011 to 2019, which was estimated to amount to approximately 83 animals, which likely represents only a fraction of the total number of cases and therefore the risk may be substantially greater. Approximately half of the cases resulted in successful prosecution and of these, the highest sentence given was 2 years in jail and a fine of IDR50mil (~USD3300). The majority of the penalties handed down for these crimes were far below the maximum potential penalties and are unlikely to be effective deterrents.

The Sicilian Wolf: Genetic Identity of a Recently Extinct Insular Population

Historically, many local grey wolf (Canis lupus) populations have undergone substantial reductions in size or become extinct. Among these, the wolf population once living in Sicily, the largest island in the Mediterranean Sea, was completely eradicated by human activity in the early decades of the 20th century. To gain a better understanding of the genetic identity of the Sicilian wolf, we used techniques for the study of ancient DNA to analyze the mitochondrial (mt) variability of six specimens stored in Italian museums. We were able to amplify a diagnostic mtDNA fragment of the control region (CR) in four of the samples. Two of the samples shared the same haplotype, differing by two substitutions from the currently most diffused Italian wolf haplotype (W14) and one substitution from the only other Italian haplotype (W16). The third sample showed a previously unreported wolf-like haplotype, and the fourth a haplotype commonly found in dogs. All of the wolf haplotypes analyzed in this study belonged to the mitochondrial haplogroup that includes haplotypes detected in all the known European Pleistocene wolves and in several modern southern European populations. Unfortunately, this endemic island population, which exhibited unique mtDNA variability, was definitively lost before it was possible to understand its taxonomic uniqueness and conservational value.

Possible Male Infanticide in Wild Orangutans and a Re-evaluation of Infanticide Risk

Infanticide as a male reproductive tactic is widespread across mammals, and is particularly prevalent in catarrhine primates. While it has never been observed in wild orangutans, infanticide by non-sire males has been predicted to occur due to their extremely long inter-birth intervals, semi-solitary social structure, and the presence of female counter-tactics to infanticide. Here, we report on the disappearance of a healthy four-month-old infant, along with a serious foot injury suffered by the primiparous mother. No other cases of infant mortality have been observed at this site in 30 years of study. Using photographic measurements of the injury, and information on the behavior and bite size of potential predators, we evaluate the possible causes of this injury. The context, including the behavior of the female and the presence of a new male at the time of the injury, lead us to conclude that the most likely cause of the infant loss and maternal injury was male infanticide. We suggest that in orangutans, and other species where nulliparous females are not preferred mates, these females may be less successful at using paternity confusion as an infanticide avoidance tactic, thus increasing the likelihood of infanticide of their first-born infants.

Phylogeography of the termite Macrotermes gilvus and insight into ancient dispersal corridors in Pleistocene Southeast Asia

Dispersal of soil-dwelling organisms via the repeatedly exposed Sunda shelf through much of the Pleistocene in Southeast Asia has not been studied extensively, especially for invertebrates. Here we investigated the phylogeography of an endemic termite species, Macrotermes gilvus (Hagen), to elucidate the spatiotemporal dynamics of dispersal routes of terrestrial fauna in Pleistocene Southeast Asia. We sampled 213 termite colonies from 66 localities throughout the region. Independently inherited microsatellites and mtDNA markers were used to infer the phylogeographic framework of M. gilvus. Discrete phylogeographic analysis and molecular dating based on fossil calibration were used to infer the dynamics of M. gilvus dispersal in time and space across Southeast Asia. We found that the termite dispersal events were consistently dated within the Pleistocene time frame. The dispersal pattern was multidirectional, radiating eastwards and southwards out of Indochina, which was identified as the origin for dispersal events. We found no direct dispersal events between Sumatra and Borneo despite the presence of a terrestrial connection between them during the Pleistocene. Instead, central Java served as an important link allowing termite colonies to be established in Borneo and Sumatra. Our findings support the hypothesis of a north-south dispersal corridor in Southeast Asia and suggest the presence of alternative dispersal routes across Sundaland during the Pleistocene. For the first time, we also propose that a west-east dispersal through over-water rafting likely occurred across the Pleistocene South China Sea. We found at least two independent entry routes for terrestrial species to infiltrate Sumatra and Borneo at different times.

Phylogeography of red muntjacs reveals three distinct mitochondrial lineages

BACKGROUND

The members of the genus Muntiacus are of particular interest to evolutionary biologists due to their extreme chromosomal rearrangements and the ongoing discussions about the number of living species. Red muntjacs have the largest distribution of all muntjacs and were formerly considered as one species. Karyotype differences led to the provisional split between the Southern Red Muntjac (Muntiacus muntjak) and the Northern Red Muntjac (M. vaginalis), but uncertainties remain as, so far, no phylogenetic study has been conducted. Here, we analysed whole mitochondrial genomes of 59 archival and 16 contemporaneous samples to resolve uncertainties about their taxonomy and used red muntjacs as model for understanding the evolutionary history of other species in Southeast Asia.

RESULTS

We found three distinct matrilineal groups of red muntjacs: Sri Lankan red muntjacs (including the Western Ghats) diverged first from other muntjacs about 1.5 Mya; later northern red muntjacs (including North India and Indochina) and southern red muntjacs (Sundaland) split around 1.12 Mya. The diversification of red muntjacs into these three main lineages was likely promoted by two Pleistocene barriers: one through the Indian subcontinent and one separating the Indochinese and Sundaic red muntjacs. Interestingly, we found a high level of gene flow within the populations of northern and southern red muntjacs, indicating gene flow between populations in Indochina and dispersal of red muntjacs over the exposed Sunda Shelf during the Last Glacial Maximum.

CONCLUSIONS

Our results provide new insights into the evolution of species in South and Southeast Asia as we found clear genetic differentiation in a widespread and generalist species, corresponding to two known biogeographical barriers: The Isthmus of Kra and the central Indian dry zone. In addition, our molecular data support either the delineation of three monotypic species or three subspecies, but more importantly these data highlight the conservation importance of the Sri Lankan/South Indian red muntjac.

Two species of Southeast Asian cats in the genus Catopuma with diverging histories: an island endemic forest specialist and a widespread habitat generalist

Background. The bay cat Catopuma badia is endemic to Borneo, whereas its sister species the Asian golden cat Catopuma temminckii is distributed from the Himalayas and southern China through Indochina, Peninsular Malaysia and Sumatra. Based on morphological data, up to five subspecies of the Asian golden cat have been recognized, but a taxonomic assessment, including molecular data and morphological characters, is still lacking. Results. We combined molecular data (whole mitochondrial genomes), morphological data (pelage) and species distribution projections (up to the Late Pleistocene) to infer how environmental changes may have influenced the distribution of these sister species over the past 120 000 years. The molecular analysis was based on sequenced mitogenomes of 3 bay cats and 40 Asian golden cats derived mainly from archival samples. Our molecular data suggested a time of split between the two species approximately 3.16 Ma and revealed very low nucleotide diversity within the Asian golden cat population, which supports recent expansion of the population. Discussion. The low nucleotide diversity suggested a population bottleneck in the Asian golden cat, possibly caused by the eruption of the Toba volcano in Northern Sumatra (approx. 74 kya), followed by a continuous population expansion in the Late Pleistocene/Early Holocene. Species distribution projections, the reconstruction of the demographic history, a genetic isolation-by-distance pattern and a gradual variation of pelage pattern support the hypothesis of a post-Toba population expansion of the Asian golden cat from south China/Indochina to Peninsular Malaysia and Sumatra. Our findings reject the current classification of five subspecies for the Asian golden cat, but instead support either a monotypic species or one comprising two subspecies: (i) the Sunda golden cat, distributed south of the Isthmus of Kra: C. t. temminckii and (ii) Indochinese, Indian, Himalayan and Chinese golden cats, occurring north of the Isthmus: C. t. moormensis.

Planning tiger recovery: Understanding intraspecific variation for effective conservation

Although significantly more money is spent on the conservation of tigers than on any other threatened species, today only 3200 to 3600 tigers roam the forests of Asia, occupying only 7% of their historical range. Despite the global significance of and interest in tiger conservation, global approaches to plan tiger recovery are partly impeded by the lack of a consensus on the number of tiger subspecies or management units, because a comprehensive analysis of tiger variation is lacking. We analyzed variation among all nine putative tiger subspecies, using extensive data sets of several traits [morphological (craniodental and pelage), ecological, molecular]. Our analyses revealed little variation and large overlaps in each trait among putative subspecies, and molecular data showed extremely low diversity because of a severe Late Pleistocene population decline. Our results support recognition of only two subspecies: the Sunda tiger, Panthera tigris sondaica, and the continental tiger, Panthera tigris tigris, which consists of two (northern and southern) management units. Conservation management programs, such as captive breeding, reintroduction initiatives, or trans-boundary projects, rely on a durable, consistent characterization of subspecies as taxonomic units, defined by robust multiple lines of scientific evidence rather than single traits or ad hoc descriptions of one or few specimens. Our multiple-trait data set supports a fundamental rethinking of the conventional tiger taxonomy paradigm, which will have profound implications for the management of in situ and ex situ tiger populations and boost conservation efforts by facilitating a pragmatic approach to tiger conservation management worldwide.

Is the clouded leopard Neofelis nebulosa extinct in Taiwan, and could it be reintroduced? An assessment of prey and habitat

During 1997–2012 we conducted a nationwide camera-trapping survey and assessed the availability of prey and habitat for the clouded leopard Neofelis nebulosa in Taiwan. We surveyed 1,249 camera-trap sites over 113,636 camera-trap days, from the seashore to an altitude of 3,796 m and covering various types of vegetation. No clouded leopards were photographed during 128,394 camera-trap days, including at 209 sites in other studies, confirming the presumed extinction of clouded leopards in Taiwan. Assessment of the prey base revealed altitudinal distribution patterns of prey species and prey biomass. Areas at lower altitudes and with less human encroachment and hunting supported a higher prey biomass and more of the typical prey species of clouded leopards. Habitat analysis revealed 8,523 km2 of suitable habitat but this was reduced to 6,734 km2 when adjacent areas of human encroachment were subtracted. In the absence of hunting and large mammalian carnivores the major prey of clouded leopards in Taiwan, such as Formosan macaques Macaca cyclopis, Reeves's muntjacs Muntiacus reevesi, Formosan serow Capricornis swinhoei and sambar Rusa unicolor, could become over-abundant. Thus, it is important to address the cascading effect of the disappearance of top-down predator control. Our assessment indicated that, with proper regulation of hunting, habitat restoration and corridor improvement, it may be possible to reintroduce the clouded leopard.

Repeated landmass reformation limits diversification in the widespread littoral zone mosquito Anopheles sundaicus sensu lato in the Indo-Oriental Region

Southeast Asia harbours abundant biodiversity, hypothesized to have been generated by Pliocene and Pleistocene climatic and environmental change. Vicariance between the island of Borneo, the remaining Indonesian archipelago and mainland Southeast Asia caused by elevated sea levels during interglacial periods has been proposed to lead to diversification in the littoral zone mosquito Anopheles (Cellia) sundaicus (Rodenwaldt) sensu lato. To test this biogeographical hypothesis, we inferred the population history and assessed gene flow of A. sundaicus s.l. sampled from 18 populations across its pan-Asian species range, using sequences from mitochondrial cytochrome c oxidase subunit 1 (CO1), the internal transcribed spacer 2 (ITS2) and the mannose phosphate isomerase (Mpi) gene. A hypothesis of ecological speciation for A. sundaicus involving divergent adaptation to brackish and freshwater larval habitats was also previously proposed, based on a deficiency of heterozygotes for Mpi allozyme alleles in sympatry. This hypothesis was not supported by Mpi sequence data, which exhibited no fixed differences between brackish and freshwater larval habitats. Mpi and CO1 supported the presence of up to eight genetically distinct population groupings. Counter to the hypothesis of three allopatric species, divergence was often no greater between Borneo, Sumatra/Java and the Southeast Asian mainland than it was between genetic groupings within these landmasses. An isolation-with-migration (IM) model indicates recurrent gene flow between the current major landmasses. Such gene flow would have been possible during glacial periods when the current landmasses merged, presenting opportunities for dispersal along expanding and contracting coastlines. Consequently, Pleistocene climatic variation has proved a homogenizing, rather than diversifying, force for A. sundaicus diversity.

Nuclear and chloroplast DNA phylogeography reveals Pleistocene divergence and subsequent secondary contact of two genetic lineages of the tropical rainforest tree species Shorea leprosula (Dipterocarpaceae) in South-East Asia

Tropical rainforests in South-East Asia have been affected by climatic fluctuations during past glacial eras. To examine how the accompanying changes in land areas and temperature have affected the genetic properties of rainforest trees in the region, we investigated the phylogeographic patterns of a widespread dipterocarp species, Shorea leprosula. Two types of DNA markers were used: expressed sequence tag-based simple sequence repeats and chloroplast DNA (cpDNA) sequence variations. Both sets of markers revealed clear genetic differentiation between populations in Borneo and those in the Malay Peninsula and Sumatra (Malay/Sumatra). However, in the south-western part of Borneo, genetic admixture of the lineages was observed in the two marker types. Coalescent simulation based on cpDNA sequence variation suggested that the two lineages arose 0.28-0.09 million years before present and that following their divergence migration from Malay/Sumatra to Borneo strongly exceeded migration in the opposite direction. We conclude that the genetic structure of S. leprosula was largely formed during the middle Pleistocene and was subsequently modified by eastward migration across the subaerially exposed Sunda Shelf.

Density of the Vulnerable Sunda clouded leopard Neofelis diardi in two commercial forest reserves in Sabah, Malaysian Borneo

Recently the Sunda clouded leopard Neofelis diardi was recognized as a separate species distinct from the clouded leopard Neofelis nebulosa of mainland Asia. Both species are categorized as Vulnerable on the IUCN Red List. Little is known about the newly identified species and, in particular, information from forests outside protected areas is scarce. Here we present one of the first density estimates calculated with spatial capture–recapture models using camera-trap data. In two commercial forest reserves in Sabah (both certified for their sustainable management practices) the density of the Sunda clouded leopard was estimated to be c. 1 per 100 km2 (0.84±SE 0.42 and 1.04±SE 0.58). The presence of the Sunda clouded leopard in such forests is encouraging for its conservation but additional studies from other areas, including protected forests, are needed to compare and evaluate these densities.