Comparative Phylogeography of Forest-Dependent Mammals Reveals Paleo-Forest Corridors throughout Sundaland

Mitogenomic analysis of Thai Sunda pangolins reveals regional phylogeography and informs conservation management

Pangolins are considered the most trafficked mammals in the world with all eight species listed on CITES Appendix I. Despite this pervasive threat to their survival, there remains a limited understanding of genetic diversity and connectivity among populations of Asian pangolin species, hampering effective conservation management. We analysed mitogenome sequences of the Sunda pangolin (Manis javanica) from across their Southeast Asia continental distribution, as well as Borneo. Phylogenetic reconstruction revealed six lineages, with clear separation north and south of the Kangar-Pattani biogeographic line in southern Thailand, revealing clear differentiation between Sundaland and Indochinese Sunda pangolin lineages. Further divergence across an east-west divide was observed in central and northern Thailand, extending northwards towards China. Our results provide new insights into the evolutionary relationships among Sunda pangolin populations in Southeast Asia, building on other recent research in this field and helping to establish the species' baseline phylogeography. These inferences will aid conservation planning and support the genetic traceability of the illegal pangolin trade.

OrthoMaM v12: a database of curated single-copy ortholog alignments and trees to study mammalian evolutionary genomics

To date, the databases built to gather information on gene orthology do not provide end-users with descriptors of the molecular evolution information and phylogenetic pattern of these orthologues. In this context, we developed OrthoMaM, a database of ORTHOlogous MAmmalian Markers describing the evolutionary dynamics of coding sequences in mammalian genomes. OrthoMaM version 12 includes 15,868 alignments of orthologous coding sequences (CDS) from the 190 complete mammalian genomes currently available. All annotations and 1-to-1 orthology assignments are based on NCBI. Orthologous CDS can be mined for potential informative markers at the different taxonomic levels of the mammalian tree. To this end, several evolutionary descriptors of DNA sequences are provided for querying purposes (e.g. base composition and relative substitution rate). The graphical web interface allows the user to easily browse and sort the results of combined queries. The corresponding multiple sequence alignments and ML trees, inferred using state-of-the art approaches, are available for download both at the nucleotide and amino acid levels. OrthoMaM v12 can be used by researchers interested either in reconstructing the phylogenetic relationships of mammalian taxa or in understanding the evolutionary dynamics of coding sequences in their genomes. OrthoMaM is available for browsing, querying and complete or filtered download at https://orthomam.mbb.cnrs.fr/.

Genomic analysis reveals a cryptic pangolin species

Eight extant species of pangolins are currently recognized. Recent studies found that two mitochondrial haplotypes identified in confiscations in Hong Kong could not be assigned to any known pangolin species, implying the existence of a species. Here, we report that two additional mitochondrial haplotypes identified in independent confiscations from Yunnan align with the putative species haplotypes supporting the existence of this mysterious species/population. To verify the new species scenario we performed a comprehensive analysis of scale characteristics and 138 whole genomes representing all recognized pangolin species and the cryptic new species, 98 of which were generated here. Our morphometric results clearly attributed this cryptic species to Asian pangolins (Manis sp.) and the genomic data provide robust and compelling evidence that it is a pangolin species distinct from those recognized previously, which separated from the Philippine pangolin and Malayan pangolin over 5 Mya. Our study provides a solid genomic basis for its formal recognition as the ninth pangolin species or the fifth Asian one, supporting a new taxonomic classification of pangolins. The effects of glacial climate changes and recent anthropogenic activities driven by illegal trade are inferred to have caused its population decline with the genomic signatures showing low genetic diversity, a high level of inbreeding, and high genetic load. Our finding greatly expands current knowledge of pangolin diversity and evolution and has vital implications for conservation efforts to prevent the extinction of this enigmatic and endangered species from the wild.

Phylogeography of the Sunda pangolin, Manis javanica: Implications for taxonomy, conservation management and wildlife forensics

The Sunda pangolin (Manis javanica) is the most widely distributed Asian pangolin species, occurring across much of Southeast Asia and in southern China. It is classified as Critically Endangered and is one of the most trafficked mammals in the world, which not only negatively impacts wild Sunda pangolin populations but also poses a potential disease risk to other species, including humans and livestock. Here, we aimed to investigate the species' phylogeography across its distribution to improve our understanding of the species' evolutionary history, elucidate any taxonomic uncertainties and enhance the species' conservation genetic management and potential wildlife forensics applications. We sequenced mtDNA genomes from 23 wild Sunda pangolins of known provenance originating from Malaysia to fill sampling gaps in previous studies, particularly in Borneo. To conduct phylogenetic and population genetic analyses of Sunda pangolins across their range, we integrated these newly generated mitochondrial genomes with previously generated mtDNA and nuclear DNA data sets (RAD-seq SNP data). We identified an evolutionarily distinct mtDNA lineage in north Borneo, estimated to be ~1.6 million years divergent from lineages in west/south Borneo and the mainland, comparable to the divergence time from the Palawan pangolin. There appeared to be mitonuclear discordance, with no apparent genetic structure across Borneo based on analysis of nuclear SNPs. These findings are consistent with the 'out of Borneo hypothesis', whereby Sunda pangolins diversified in Borneo before subsequently migrating throughout Sundaland, and/or a secondary contact scenario between mainland and Borneo. We have elucidated possible taxonomic issues in the Sunda/Palawan pangolin complex and highlight the critical need for additional georeferenced samples to accurately apportion its range-wide genetic variation into appropriate taxonomic and conservation units. Additionally, these data have improved forensic identification testing involving these species and permit the implementation of geographic provenance testing in some scenarios.

Biogeography of the Sunda Shelf revisited: Insights from Macaranga section Pruinosae (Euphorbiaceae)

The Southeast Asian region of Sundaland is among the world’s major biodiversity hotspots. The region’s biodiversity coupled with its complex and dynamic geographic and climatic histories makes it an ideal region to study the various factors that determine the diversification and distribution patterns of tropical biota. Here we investigate the biogeographic patterns in the partly myrmecophytic Macaranga section Pruinosae to reveal some of the factors that play a role in determining the distribution of biota in Sundaland. We used single nucleotide polymorphisms (SNP) data derived from GBS, a next generation sequencing technique, in maximum likelihood and cluster analyses to determine phylogenetic relationships and population structures within this taxonomic section. Bayesian inference based on secondary calibration points and ancestral area reconstruction analyses were performed to infer spatial–temporal origins of the major lineages in the section. The results from these analyses were further substantiated using nuclear microsatellite data obtained from a broader sample set of two widely distributed species within the section: Macaranga gigantea and Macaranga pruinosa. Phylogenetic and cluster analyses reveal four well-defined, discrete species groups within section Pruinosae, all of which but one originated in Borneo with the crown node dated at 3.58 mya. Biogeographic patterns within the species reveal a biogeographic barrier between east and west Sundaland besides bringing to light the role played by various geological factors, especially the Crocker Range, on Borneo. Patterns also reveal a biogeographic barrier between the Bangka/Belitung islands and Sumatra for ant-free, swamp-adapted species. This study provides evidence that geographic barriers, edaphic conditions, and ecological adaptations are tightly linked and that their mutual interaction determines the diversification and distribution of species.

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.

Island Biogeography Revisited: Museomics Reveals Affinities of Shelf Island Birds Determined by Bathymetry and Paleo-Rivers, Not by Distance to Mainland

Island biogeography is one of the most powerful subdisciplines of ecology: its mathematical predictions that island size and distance to mainland determine diversity have withstood the test of time. A key question is whether these predictions follow at a population-genomic level. Using rigorous ancient-DNA protocols, we retrieved approximately 1,000 genomic markers from approximately 100 historic specimens of two Southeast Asian songbird complexes from across the Sunda Shelf archipelago collected 1893-1957. We show that the genetic affinities of populations on small shelf islands defy the predictions of geographic distance and appear governed by Earth-historic factors including the position of terrestrial barriers (paleo-rivers) and persistence of corridors (Quaternary land bridges). Our analyses suggest that classic island-biogeographic predictors may not hold well for population-genomic dynamics on the thousands of shelf islands across the globe, which are exposed to dynamic changes in land distribution during Quaternary climate change.

Divergence and introgression in small apes, the genus Hylobates, revealed by reduced representation sequencing

Gibbons of the genus Hylobates, which inhabit Southeast Asia, show great diversity and comprise seven to nine species. Natural hybridisation has been observed in several species contact zones, but the history and extent of hybridisation and introgression in possibly historical and the current contact zones remain unclear. To uncover Hylobates species phylogeny and the extent of introgression in their evolution, genotyping by random amplicon sequencing-direct (GRAS-Di) was applied to 47 gibbons, representing seven Hylobates species/subspecies and two outgroup gibbon species. Over 200,000 autosomal single-nucleotide variant sites were identified. The autosomal phylogeny supported that divergence from the mainland species began ~3.5 million years ago, and subsequently occurred among the Sundaic island species. Significant introgression signals were detected between H. lar and H. pileatus, H. lar and H. agilis and H. albibarbis and H. muelleri, which all are parapatric and form ongoing hybrid zones. Furthermore, the introgression signals were detected in every analysed individual of these species, indicating a relatively long history of hybridisation, which might have affected the entire gene pool. By contrast, signals of introgression were either not detected or doubtful in other species pairs living on different islands, indicating the rarity of hybridisation and introgression, even though the Sundaic islands were connected during the Pliocene and Pleistocene glacial events.

Evolutionary history of Sundaland shrews (Eulipotyphla: Soricidae: Crocidura) with a focus on Borneo

The hyperdiverse shrew genus Crocidura is one of few small mammal genera distributed across Sundaland and all of its boundaries. This represents a rare opportunity to study the geological history of this region through the evolutionary history of these shrews. We generate a phylogeny of all recognized species of Sundaland Crocidura and show that most speciation events took place during the Pleistocene, prior to the inundation of the Sunda Shelf around 400 000 years ago. We find east–west differentiation within two separate lineages on Borneo, and that the current taxonomy of its two endemic species does not reflect evolutionary history, but ecophenotypic variation of plastic traits related to elevation. Sulawesi shrews are monophyletic, with a single notable exception: the black-footed shrew (C. nigripes). We show that the black-footed shrew diverged from its relatives on Borneo recently, suggesting a human-assisted breach of Wallace’s line. Overall, the number of Crocidura species, especially on Borneo, probably remains an underestimate.

Brown banded bamboo shark (Chiloscyllium punctatum) shows high genetic diversity and differentiation in Malaysian waters

The demersal brown banded bamboo shark Chiloscyllium punctatum is a major component of sharks landed in Malaysia. However, little is known about their population structure and the effect of high fishing pressure on these weak swimming sharks. Both mitochondrial DNA control region (1072 bp) and NADH dehydrogenase subunit 2 (1044 bp) were used to elucidate the genetic structure and connectivity of C. punctatum among five major areas within the Sundaland region. Our findings revealed (i) strong genetic structure with little present day mixing between the major areas, (ii) high intra-population genetic diversity with unique haplotypes, (iii) significant correlation between genetic differentiation and geographical distance coupled with detectable presence of fine scale geographical barriers (i.e. the South China Sea), (iv) historical directional gene flow from the east coast of Peninsular Malaysia towards the west coast and Borneo, and (v) no detectable genetic differentiation along the coastline of east Peninsular Malaysia. Genetic patterns inferred from the mitochondrial DNA loci were consistent with the strong coastal shelf association in this species, the presence of contemporary barriers shaped by benthic features, and limited current-driven egg dispersal. Fine scale population structure of C. punctatum highlights the need to improve genetic understanding for fishery management and conservation of other small-sized sharks.

A new subspecies of Trypanosoma cyclops found in the Australian terrestrial leech Chtonobdella bilineata

Previously, it was suggested that haemadipsid leeches represent an important vector of trypanosomes amongst native animals in Australia. Consequently, Chtonobdella bilineata leeches were investigated for the presence of trypanosome species by polymerase chain reaction (PCR), DNA sequencing and in vitro isolation. Phylogenetic analysis ensued to further define the populations present. PCR targeting the 28S rDNA demonstrated that over 95% of C. bilineata contained trypanosomes; diversity profiling by deep amplicon sequencing of 18S rDNA indicated the presence of four different clusters related to the Trypanosoma (Megatrypanum) theileri. Novy–MacNeal–Nicolle slopes with liquid overlay were used to isolate trypanosomes into culture that proved similar in morphology to Trypanosoma cyclops in that they contained a large numbers of acidocalcisomes. Phylogeny of 18S rDNA/GAPDH/ND5 DNA sequences from primary cultures and subclones showed the trypanosomes were monophyletic, with T. cyclops as a sister group. Blood-meal analysis of leeches showed that leeches primarily contained blood from swamp wallaby (Wallabia bicolour), human (Homo sapiens) or horse (Equus sp.). The leech C. bilineata is a host for at least five lineages of Trypanosoma sp. and these are monophyletic with T. cyclops; we propose Trypanosoma cyclops australiensis as a subspecies of T. cyclops based on genetic similarity and biogeography considerations.

Gene flow creates a mirage of cryptic species in a Southeast Asian spotted stream frog complex

Most new cryptic species are described using conventional tree- and distance-based species delimitation methods (SDMs), which rely on phylogenetic arrangements and measures of genetic divergence. However, although numerous factors such as population structure and gene flow are known to confound phylogenetic inference and species delimitation, the influence of these processes is not frequently evaluated. Using large numbers of exons, introns, and ultraconserved elements obtained using the FrogCap sequence-capture protocol, we compared conventional SDMs with more robust genomic analyses that assess population structure and gene flow to characterize species boundaries in a Southeast Asian frog complex (Pulchrana picturata). Our results showed that gene flow and introgression can produce phylogenetic patterns and levels of divergence that resemble distinct species (up to 10% divergence in mitochondrial DNA). Hybrid populations were inferred as independent (singleton) clades that were highly divergent from adjacent populations (7%-10%) and unusually similar (<3%) to allopatric populations. Such anomalous patterns are not uncommon in Southeast Asian amphibians, which brings into question whether the high levels of cryptic diversity observed in other amphibian groups reflect distinct cryptic species-or, instead, highly admixed and structured metapopulation lineages. Our results also provide an alternative explanation to the conundrum of divergent (sometimes nonsister) sympatric lineages-a pattern that has been celebrated as indicative of true cryptic speciation. Based on these findings, we recommend that species delimitation of continuously distributed "cryptic" groups should not rely solely on conventional SDMs, but should necessarily examine population structure and gene flow to avoid taxonomic inflation.

Ancient Divergence Driven by Geographic Isolation and Ecological Adaptation in Forest Dependent Sundaland Tree Squirrels

A surprising amount of hidden phylogenetic diversity exists in the small to medium size, drab colored squirrels of the genus Sundasciurus. This genus is endemic to Sundaland and the Philippines, where it is widespread. An earlier revision of this genus found that the high elevation ‘populations’ of the widespread, lowland slender squirrel (S. tenuis) were different species. Previous phylogenies based on mitochondrial cytochrome b sequences also suggested that the widespread, lowland Low’s squirrel (S. lowii) and the narrow endemic Fraternal squirrel (S. fraterculus) are not reciprocally monophyletic. Additionally, deep divergences have been identified between lineages within Low’s squirrel that date to the early Pliocene. Here we focus on evaluating the relationships and differences within and between populations of these two nominal species using whole mitochondrial genome sequences, nuclear intron sequences, and morphology. We reassess the taxonomy of this group, revalidate the species status of Robinson’s squirrel (Sundasciurus robinsoniBonhote, 1903) support the species level recognition of the Natuna squirrel (Sundasciurus natunensisThomas, 1895) and identify three other lineages that require further study. We estimate times of divergence and integrate geologic history to find that most of the divergences are pre-Pleistocene, and thus predate the Pleistocene flooding of Sundaland. Biogeographic, and ecological factors may have played a more important role than climatic factors in generating these patterns. While divergence in allopatry seems to be the main process driving speciation in lowland Sundaland squirrels (Sundasciurus), ecomorphological and behavioral adaptations in this clade suggest an important role of niche divergence.

Quaternary land bridges have not been universal conduits of gene flow

Quaternary climate oscillations are a well-known driver of animal diversification, but their effects are most well studied in areas where glaciations lead to habitat fragmentation. In large areas of the planet, however, glaciations have had the opposite effect, but here their impacts are much less well understood. This is especially true in Southeast Asia, where cyclical changes in land distribution have generated enormous land expansions during glacial periods. In this study, we selected a panel of five songbird species complexes covering a range of ecological specificities to investigate the effects Quaternary land bridges have had on the connectivity of Southeast Asian forest biota. Specifically, we combined morphological and bioacoustic analysis with an arsenal of population genomic and modelling approaches applied to thousands of genome-wide DNA markers across a total of more than 100 individuals. Our analyses show that species dependent on forest understorey exhibit deep differentiation between Borneo and western Sundaland, with no evidence of gene flow during the land bridges accompanying the last 1-2 ice ages. In contrast, dispersive canopy species and habitat generalists have experienced more recent gene flow. Our results argue that there remains much cryptic species-level diversity to be discovered in Southeast Asia even in well-known animal groups such as birds, especially in nondispersive forest understorey inhabitants. We also demonstrate that Quaternary land bridges have not been equally suitable conduits of gene flow for all species complexes and that life history is a major factor in predicting relative population divergence time across Quaternary climate fluctuations.

Genetic diversity and structure of the binturong Arctictis binturong (Carnivora: Viverridae) – status of the elusive Palawan binturong and implications for conservation

The binturong or bearcat is a forest mesocarnivore ranging from Nepal to Indonesia and the Philippines. Several subspecies of binturongs are recognized but a revision is needed. The binturong from Palawan was described as a species and is now considered a subspecies, but its status has never been checked using molecular approaches. Owing to its restricted range and the pressure on its habitat, the Palawan binturong may be endangered. It is, therefore, of crucial importance to clarify its taxonomic status, particularly for the management of captive populations. We sequenced one nuclear and two mitochondrial markers for binturongs from locations across the species range and from zoos. Our results provide an assessment of the genetic polymorphism and structure within the binturong, resulting in two groups, corresponding to the Indochinese and the Sundaic regions. Within the latter were found the Palawan binturongs on one side, and an individual from Sulu archipelago (a locality not reported before) on the other side. The Palawan binturongs form a monophyletic group, genetically close to Bornean binturongs, which suggests that they may have dispersed from Borneo, and represents a lineage worth preserving, but which is not a separate species nor a separate subspecies.