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.
