A phylogeny for the pomatiopsidae (Gastropoda: Rissooidea): a resource for taxonomic, parasitological and biodiversity studies
© Liu et al.; licensee BioMed Central Ltd. 2014
Received: 30 September 2013
Accepted: 7 February 2014
Published: 18 February 2014
The Pomatiopsidae are reported from northern India into southern China and Southeast Asia, with two sub-families, the Pomatiopsinae (which include freshwater, amphibious, terrestrial and marine species) and the freshwater Triculinae. Both include species acting as intermediate host for species of the blood-fluke Schistosoma which cause a public health problem in East Asia. Also, with around 120 species, triculine biodiversity exceeds that of any other endemic freshwater molluscan fauna. Nevertheless, the origins of the Pomatiopsidae, the factors driving such a diverse radiation and aspects of their co-evolution with Schistosoma are not fully understood. Many taxonomic questions remain; there are problems identifying medically relevant species. The predicted range is mostly unsurveyed and the true biodiversity of the family is underestimated. Consequently, the aim of the study was to collect DNA-sequence data for as many pomatiopsid taxa as possible, as a first step in providing a resource for identification of epidemiologically significant species (by non-malacologists), for use in resolving taxonomic confusion and for testing phylogeographical hypotheses.
The evolutionary radiation of the Triculinae was shown to have been rapid and mostly post late Miocene. Molecular dating indicated that the radiation of these snails was driven first by the uplift of the Himalaya and onset of a monsoon system, and then by late-Pliocene global warming. The status of Erhaia as Anmicolidae is supported. The genera Tricula and Neotricula are shown to be non-monophyletic and the tribe Jullieniini may be polyphyletic (based on convergent characters). Triculinae from northern Vietnam could be derived from Gammatricula of Fujian/Yunnan, China.
The molecular dates and phylogenetic estimates in this study are consistent with an Australasian origin for the Pomatiopsidae and an East to West radiation via Oligocene Borneo-Philippines island hopping to Japan and then China (Triculinae arising mid-Miocene in Southeast China), and less so with a triculine origin in Tibet. The lack of monophyly in the medically important genera and indications of taxonomic inaccuracies, call for further work to identify epidemiologically significant taxa (e.g., Halewisia may be potential hosts for Schistosoma mekongi) and highlight the need for surveys to determine the true biodiversity of the Triculinae.
KeywordsPomatiopsidae Triculinae Neotricula Schistosomiasis Phylogeny Phylogeography Taxonomy Amnicolidae Biodiversity
The importance of the pomatiopsidae
In addition to their parasitological importance, the Triculinae are of great interest owing to their biodiversity, conservation status, and the extent of their adaptive radiation. Radiations of Triculinae occur in Yunnan/Sichuan (Southwest China), 11 species; Hunan (Southeast China), 13 species (mainly Neotricula) [8, 9]; and an endemic fauna that includes over 90 species occurs along a 300 km stretch of the lower Mekong river in Thailand and Laos . The biodiversity of the Triculinae exceeds that of any other endemic freshwater molluscan fauna. If the endemic gastropods are counted for the ancient lakes, Baikal, Biwa, Ohrid, Tanganyika, and Titicaca, the species of the Mekong river Triculinae outnumber the entire endemic fauna in any one lake and show more than double the number of species found in any one endemic family in a lake. Gastropod radiations in lotic systems are generally more speciose than those in lakes and the Alabama River system (Alabama, Cahaba, Coosa, and Mobile rivers) radiation of the Southeast USA is considered to be a remarkable radiation of endemic freshwater molluscs ; however, the lower Mekong river triculine radiation exceeds even this . In addition, the low vagility, high levels of endemism, strict habitat specialisation and sensitivity to environmental disturbance make triculine snails good indicators of near-pristine or long-term undisturbed habitats. All Pomatiopsidae are calciphilic and many are associated with endangered karst habitats; this further enhances their conservation significance.
In spite of the importance of this group many questions remain regarding their biology. For example, the origins of the Pomatiopsidae, the factors driving such a diverse radiation and aspects of their co-evolution with Schistosoma are not fully understood. Many taxonomic questions remain and there are still problems assigning known taxa to genera, further, with much of their predicted range unsurveyed (most triculines are highland species), the true biodiversity of the family is certainly greatly underestimated. Public health surveillance is inhibited by confusion between snail intermediate hosts of trematodes and similar looking sympatric species; there is also confusion as to which snails are acting as hosts for Schistosoma. The problems of small size (most Triculinae are < 3 mm in shell height), for many species an inability to breed in the laboratory, remoteness and inaccessibility of habitats, a high prevalence of convergent evolution and a relatively high degree of intraspecific variation among anatomical characters, all hamper studies of the Triculinae.
Questions remain to be addressed at all taxonomic levels concerning the Triculinae. The currently accepted taxonomy is based on that of Brandt  which has been extensively revised by G.M. Davis (e.g., ). The Triculinae were originally divided into three tribes based on a cladistic analysis which used morphological characters (namely, the Jullieniini, Lacunopsini and Triculini) . In 1990 these clades were revised, a new tribe, Pachydrobiini was defined, thus dividing the Triculini sensu Davis 1979 into the Triculini (Tricula, Delavaya, Fenouilia, Lacunopsis) and the Pachydrobiini (Gammatricula, Guoia, Halewisia, Jinghongia, Neotricula, Pachydrobia, Robertsiella, Wuconchona) . The Pomatiopsinae were also revised following the discovery of a new genus, Erhaia, in Yunnan, which was provisionally assigned to a new tribe of the Pomatiopsidae: Pomatiopsinae, the Erhaiini. Genera originally placed in the Pomatiopsinae (Pomatiopsis, Oncomelania, Cecina, et al.) were placed in a new tribe, the Pomatiopsini . Further revisions were necessary in 1992, with the creation of a new tribe, the Pseudobythinellini, to accommodate the genera Akiyoshia and Pseudobythinella which had been recognized as Pomatiopsinae; these included some taxa originally referred to the Erhaiini and so Erhaia was placed in synonymy with Pseudobythinella. More recently Erhaia (and other Pseudobythinellini such as Akiyoshia) were transferred to the Amnicolidae: Amnicolinae [16, 17]. Similar taxonomic problems concern the genus Manningiella, with most (maybe all) of its species re-assigned to the genera Halewisia and Hubendickia; the fact that these genera belong to different tribes further illustrates the problems encountered in triculine taxonomy. The revision of the Triculinae is a work in progress.
The current systematics of the Pomatiopsidae is largely based on morphological characters. Clearly this is a situation where DNA-sequence based phylogenies would be useful in resolving taxonomic questions. Several authors have published such phylogenies involving Pomatiopsidae, but the focus has been on relationships within medically important genera [3–5, 18–21] or the study was addressing questions on a different family  or higher divisions in the Rissooidea . A few DNA based phylogenies have been published addressing questions at higher taxonomic levels, but each has been generally restricted to single radiation. For example, a DNA based phylogeny has been published for Chinese Triculinae , but this study omitted the Mekong radiation and included only one pomatiopsine snail. The study was also adopted a disputed taxonomy (that of [25–27], which is based solely on characters known to be unreliable such as those of shell, radula and head-foot; the present study uses the better documented taxonomy of Davis  and later works by that author. The taxonomy of Davis is based on descriptions of taxa involving more reliable characters and more complete evidence (e.g., more extensive photomicrography, SEM and fine-scale dissection work). The study of  found two main clades of Triculinae, with the three taxa they sampled from Fujian, Guangxi and Zhejiang provinces in one clade (together with N. aperta), and those from Hunan, Sichuan and Hubei provinces in the other clade; this the authors related to the greater distance of Fujian, Guangxi and Zhejiang from the Yangtze river. These authors also found Tricula to be paraphyletic; however, the taxa in question were T. wumingensis and T. fujianensis, which are taxa described by Liu et al.  and Hu et al.  based on superficial characters only. The phylogeny published by Kameda and Kato  covers a broader range of taxa than that of the 2008 study and is based on both nuclear (18S and 28S rDNA) and mitochondrial (16S rDNA and cox1) DNA sequences; however, it is focused on the Pomatiopsinae and includes Triculinae sequences from GenBank for comparative purposes (these authors followed the triculine taxonomy of Liu et al. ; Tang et al. ; Hu et al., . The 2011 study provided new sequence data for Japanese Akiyoshia confirming that it is not Pomatiopsidae. Finally, there are reports, based on morphology, of Pachydrobiini from Vietnam that are described as Vietricula[29, 30]; however, these snails resemble Pachydrobia and not any kind of Tricula.
Questions regarding the origin and evolution of the pomatiopsidae
Tibet (west–east) radiation hypotheses
Phylogeographical events and divergence time predictions used to construct phylogenies corresponding to, and in accordance with, alternative hypotheses for the origin and radiation of the Triculinae and related taxa
Hypothesis … Event
Predicted date [source]
Tibet (West–east) radiation
1. Lithoglyphidae diverge from Hydrobiidae
Pangaea/fossil record 
305-169 Ma 
2. Amnicolidae and proto-Pomatiopsidae diverge
Pangaea after initial breakup into northern and Gondwanaland continents/amnicolids common and diverse in North, but less so in Asia, therefore it is assumed that they crossed into Southwest China before the major Himalayan uplift) 
200-190 Ma 
3. Pomatiopsidae arise after break up of East Gondwana
Gondwana/Fossil record indicates no Triassic rafting of pomatiopsids on West Burma block and Cimmeria to Asia
<210-150 Ma 
4. Pomatiopsinae and Triculinae diverge
Indian Craton - after separation from Gondwana (as no Triculinae in South Africa) 
5. The tribes of the Triculinae diverge
Yunnan/Prior to closure of the Brahmaputra-Irrawaddy-Mekong corridor (because all three tribes of the Triculinae are found in lower Mekong) 
c.a. 18 Ma 
6. Lacunopsis clade diverged from Triculini
Yunnan/Fenouilia of Yunnan, like Lacunopsis, is derived from proto-Tricula bollingi therefore, date lies between divergence of Triculini (18 Ma) and major regional uplift
18-7 Ma 
7. Lower Mekong triculine radiation in Laos
Triculine taxa become more derived as they radiate out from Dali/evolution occurs in concert with the evolution of the rivers as they cut southwards towards the sea 
Hunan (East to West) radiation
1. Lithoglyphidae diverge from Hydrobiidae
Pangaea/fossil record 
305-169 Ma 
2. Amnicolidae and proto-Pomatiopsidae diverge
East Gondwana/After separation from India-Madagascar (because no confirmed reports of pomatiopsids in India, either extant or fossil)
3. Pomatiopsinae diverge in marine habitats
Eastern Gondwana (Australia)/inundation by high sea levels, increased coastal habitats
c.a. 112 Ma 
4. Proto-Oncomelania arises
Tertiary island hopping along extensive island complex (Borneo-Philippines) 
40-20 Ma 
5. Radiation of terrestrial/amphibious pomatiopsines
Japan/Miocene orogeny in Japan plus local and global climate change 
6. Oncomelania enters China and diverges
On Yangtze plain/Prior to opening of Sea of Japan 
c.a. 15 Ma 
7. Triculinae diverge from proto-Oncomelania
West China/Triculinae adapt to new conditions from Pliocene major uplift of Himalaya 
c.a. 7 Ma 
8. Hunan, northern Lao & Vietnam Pachydrobiini isolate
Indosinia/Red river corridor between southern China (Hunan) and Sundaland broken 
c.a. 4.4 Ma 
9. Yunnan & Sichuan Triculinae isolate in uplifting terrane
Southwest China/Uplift of Hengduan mountains and associated ranges
c.a. 3.4 Ma 
10. Lacunopsis clade diverged from Triculini
Yunnan/Prior to Hengduan orogeny
c.a. 3.4 Ma 
11. Northern Thailand Triculini isolated from China/Sundaland taxa
12. Lower Mekong extensive radiation of Pachydrobbiini and other Triculinae
Climatic change and range contraction/fragmentation
c.a. 2.5 Ma 
13. Robertsiella diverges from Mekong Pachydrobiini in Pahang river drainage, Malaysia
2.6-0.8 Ma 
14. Divergence of Erhai & Dianchi Basin taxa in Yunnan
Yunnan/Stage 3–2 of Sanjiang orogeny in Yunnan and associated tectonic events
0.9 Ma 
15. Neotricula aperta strains diverge
Final surge of Himalayan uplift, tilting of Khorat Basin, Mul-river flow reversal, and volcanism in southern Laos , drive further divergence in the lower Mekong
0.8-0.9 Ma 
Very little fossil evidence has been presented in the literature in support of the India rafting hypothesis and Miocene entry of Pomatiopsidae into eastern Asia, from India, via Tibet. Fossils resembling Tomichia have been reported from the upper Cretaceous of South Africa ; however, this is after the currently accepted date for the separation of Africa and India, and a fossil (being shell only) provides insufficient characters for the identification of hydrobiid genera. More recently Neubauer and colleagues  reported what they described as a triculine genus from the upper most Eocene to lower Oligocene of nothern Vietnam, namely Bacbotricula gen. nov. (Pomatiopsidae: Triculinae). The fossil was assigned to the Triculinae because, despite the name Bacbotricula, it resembled jullieniinine taxa, such as Hydrorissoia, Karelainia, Neoprososthenia and Paraprososthenia, in having a turreted shell with spiral sculpture; however, no extant Triculinae shows such a marked keel and angulated parietal wall as “Bacbotricula”. On the basis of the shell alone, this taxon could be Rissooidea or small Cerithioidea, and the authors did note that their taxonomy was tentative. In the present study fossil evidence is used to date only the deepest cladogenic events, and such data are not suitable for application at higher taxonomic levels. Further, if we date the find as Priabonian-Rupelian (40–28 Ma), Bacbotricula’s arrival in Vietnam appears to pre-date the evolution of the rivers of Yunnan, which are proposed to have driven and shaped phylogenesis among the Pomatiopsidae under the Tibet hypothesis. Nevertheless, Neubauer et al., have provided data consistent with the idea that there could have been Triculinae in Oligocene Vietnam and an insight into their freshwater ecology; these authors also proposed the Red river as a possible route for triculine snails to have entered Southeast Asia from Tibet/Yunnan, thus echoing the “Red river hypothesis”  of earlier papers (see following sub-section).
Hunan (east to west) radiation hypothesis
As mentioned above the Tibet hypothesis fails to explain the current deployment of triculine taxa. For example, the lack of Tricula spp. (and Oncomelania) in India, the absence of Neotricula from Yunnan and its mono-specific radiation in Southeast Asia. Estimated dates of divergence for Asian Schistosoma spp. also conflict with the Tibet hypothesis. DNA sequence data supported a 3–6 Ma age for the Schistosoma sinensium group (the clade to which S. malayensis, S. mekongi, and S. japonicum belong ), with the divergence of S. japonicum at around 4 Ma and S. malayensis/S. mekongi at 2.5 Ma . An alternative hypothesis was proposed which explained the deployment of the Triculinae as a mid-Pliocene radiation from Hunan in southern China, westwards into Southeast Asia [9, 20]. At least 5 species of Neotricula are known from Hunan and the closest relatives of Robertsiella (Guoia and Wuconchona) are found in Hunan and not in Yunnan or Southeast Asia . The apparent “southern continental” distribution of the phylogenetically basal Pomatiopsinae can be as readily explained by correlation with karst or limestone distribution as by a Gondwanan origin  (see Figure 2B).
The 900 km long valley of the Red River fault was proposed as the mid-Pliocene dispersal corridor for proto-Neotricula through the Mesozoic Annam mountain chain which forms a barrier between Hunan and northern Laos/Vietnam . In the past the Red River fault ran 400 km further West of its current path, and therefore closer to Laos . In addition, the Yangtze river is believed to have flowed along a common course with the Red River . The Red-River radiation into Southeast Asia was dated at c.a. 4 Ma, just prior to major tectonic events in the region which may have led to a Pliocene switch from previously southwards courses to the eastwards courses seen in the major Chinese rivers today ; this change in the drainage pattern would have been a significant isolating event in the region . The Pliocene Dong-Ngai-Mekong River, which once flowed North to South across the Sunda shelf (now off shore of Vietnam’s East coast) to the East of the Annam mountains , could have introduced proto-Neotricula to Southeast Asia from Hunan. The ancestral condition of Triculinae is spring or primary stream dwelling, and this is seen in almost all Triculinae. Neotricula aperta and other lower Mekong triculines are unusual in that they have become adapted to life in large rivers such as the Mekong. This adaptation may have been a response to climatic change and range contraction/fragmentation occurring 2.5 Ma following a marked intensification of monsoon winds affecting rainfall patterns in the region . Cambodia would have been the most likely point of entry into Southeast Asia from the vast plain of wetlands and rivers that existed off the present day coast of Vietnam, with a similar plain linking Cambodia and West Malaysia. The uplift of volcanic highlands in Southeast Cambodia at 0.8 Ma  probably cut off the Vietnam-Cambodia corridor just South of the Kontum range. The Pleistocene Mekong river flowed West, along the present day border between Thailand and Cambodia, then southwards along the Tonle Sap river of today, from where it flowed across the exposed Sunda shelf towards Malaysia . The lower Mekong River is unlikely to have occupied its present course until 5–6 ka (thousand years ago) . Consequently, the Triculinae of the lower Mekong river most likely radiated in different habitats and river systems until the Holocene. Under the East to West hypothesis, the Triculinae are assumed to have diverged from the Pomatiopsinae as these taxa migrated westwards and were caught up in the rapidly uplifting terrain of the Shan and north Indosinian mountains (Yunnan/northern Laos). Tricula (and other Triculini) was proposed to have diverged in the higher elevations of Yunnan and Neotricula (and other Pachydrobiini) in the small streams of the lesser mountains of Hunan. In this model Oncomelania is regarded as diverging from other Pomatiopsidae, which are ancestrally marine or brackish water, in Australasia and colonizing China from Japan and eastern China (Figure 2B), rather than from India via Tibet and Yunnan .
The Triculinae are proposed to have arisen from a common ancestor with Pomatiopsinae which arose in the Cretaceous on eastern Gondwana (Australia-Antarctica-Papua New Guinea-New Zealand etc.) after separation of India/Madagascar (because Pomatiopsidae are absent from India and Africa). Rather than rafting via India, the proto-Pomatiopsidae are considered to have rafted from Northwest Australia on small continental fragments of Gondwanan origin which today form parts of Borneo and eastern Indonesia . Onward dispersal to the southern Philippines, then in a more southerly location than at present, would have been via island hopping across the extensive island complex that formed Oligocene North Borneo (Figure 2B). In addition, northwards rafting and rotation brought southerly Gondwanan terranes to abut the Philippines during the Mio-Pliocene. Colonization of Japan would have occurred about mid-Miocene, followed by invasion of the Yangtze plain of China from its East coast. Allozyme and DNA-sequence based phylogenies for Oncomelania show the “southern continental” pomatiopsines and Japanese O. minima as basal to all other Oncomelania[67, 68], suggesting that Oncomelania sensu stricto arose in Japan, with O. hupensis then colonizing China and back-tracking into the Philippines and Sulawesi (probably displacing its ancestral form). The proto-pomatiopsid leaving Gondwana most likely held the ancestral condition and was marine; the dispersal to south Africa and South America (originally in much closer proximity) was therefore not unlikely. Indeed, continental connections or close proximity of southern Africa and South America may have persisted until c.a. 105 Ma . The March 2011 Pacific tsunami demonstrated that large aggregates of material may cross the ocean in less than 15 months and that these may harbour viable communities of exotic aquatic organisms (including molluscs) . This hypotheses not only explains why genera such as Tomichia show relatively little diversity, but also why the basal genera in the pomatiopsid phylogeny are not those best adapted to mountainous conditions but are those adapted to marine or freshwater wet-land conditions. The phylogeny, dispersal tracts and time markers associated with this hypotheses, together with additional references and details are further explained by Table 1 and Figures 2B and 3B.
Aims of the study
In view of the many questions regarding the taxonomy of the Triculinae the objective of the present study was to collect samples of as many relevant taxa as possible and to obtain DNA sequence data for representatives of all genera which could be reliably identified by morphological criteria. These data would then be combined with existing published data and phylogenies then estimated by several independent methods so as to obtain a well supported evolutionary history for these taxa. The resulting phylogeny would be used to assess the monophyly of clades within the Triculinae. In addition, alternative hypotheses regarding the origin of the Triculinae were to be tested, and the timing and dispersal tracts of their radiation estimated. The findings would then be used to provide guidance for species identification in the Triculinae, information of use in studies of host range in Schistosoma and Paragonimus and of general investigations into host-parasite co-evolution, clarification of higher taxonomic groupings in the Pomatiopsidae, and to shed light on the nature and triggers of rapid and highly speciose radiations of freshwater invertebrate taxa. Finally, the work was expected to facilitate conservation and biodiversity studies through clarification of the numbers of taxa present, their habitat and geographical ranges and through aiding identification.
It should be noted that, although the Tibet hypothesis is based on that of Davis , the original hypothesis published by Davis was not presented in a form which could be tested relative to a Hunan hypothesis. Consequently, the Tibet hypothesis is an adaptation of the original ideas presented in Davis  and other papers. The Hunan hypothesis is also not exclusive of the ideas presented by Davis, for example it looks to Davis  for the origin of the Jullieniini and Lacunopsini and also begins with the break up of Gondwana. In view of this, the subject of this paper is not a comparison of a Davis hypothesis with that of other authors, it is more accurately a comparison of two hypotheses, both influenced by the work of Davis, albeit to different degrees.
In this study we estimate phylogenies from DNA-sequence data available at homologous loci for as comprehensive a sample of Pomatiopsidae as possible (16 out of the 32 validated pomatiopsid genera are sampled, including at least 3 taxa which are apparently extinct). A relaxed molecular clock is applied to the phylogenies and key divergence dates estimated by a Bayesian approach. The divergence dates and phylogeny are shown to agree with predictions for the Hunan hypothesis as well, or even better than, those of the traditional Tibet hypothesis. Major cladogenic events are linked to palaeo-climate and geographical changes in the region which could have driven divergence. Phylogenetic affinities are linked to dispersal and divergence along past dispersal corridors. The medically important genera Neotricula and Tricula were not found to be monophyletic and suggestions are made for their revision.
The use of vicariant calibrations and panbiogeography
The approaches used in phylogeography have been the subject of much controversy, with objections raised regarding the assumption of vicariance , dating using molecular clock rates , and the correlation of the age of taxa with the dates of palaeogeographical events . For example, a correlation between continental breakup and phylogeny may only occur where the continental biota was homogeneous prior to the tectonic events , and multiple events, occurring at different times, can generate the appearance of false and misleading vicariance (“pseudo-congruence” . The present study has attempted to avoid these problems by adopting a more panbiogeographical approach. The aim is not to use the data to infer or prove dispersal and vicariance events for the phylogeny of the Pomatiopsidae, instead the approach is to compare patterns in the phylogeny with the two main hypotheses for the historical biogeography and evolution of these snails and to assess their relative degrees of fit. Neither of these hypotheses was formulated in the present study, and the events or molecular clock rates used in this study are for the most part independent of those used to derive these hypotheses. The earliest divergences in the hypotheses do assume that cladogenesis is due to plate movement and not post-breakup dispersal or divergence in allopatry prior to isolation by continental fragmentation. Such assumptions have been heavily criticised [71, 75]; however, these criticisms often refer to speciation events. In the present study such assumptions are applied mainly to deep cladogenic events and at continental scales, where entire stem lineages fail to enter certain regions (e.g., proto-Pomatiopsidae into India). The events hypothesized are also deep in time (e.g., Cretaceous) so that the gradual nature of tectonic upheavals have less of an effect when treated as events rather than episodes, a potential problem noted by critics . In addition, the more recent dating events in the hypotheses are often seen as triggers of cladogenic eposides (e.g., climate changes after major Himalayan uplift driving triculine radiation) rather than the simple isolating events in vicariance-dispersal approaches. A further criticism of phylogeographical studies is the use of present day regions (e.g. Borneo, China, Southeast Asia) as biogeographical units, and the failure to recognize that these regions can have composite, independent, affinities , such as the past connections between parts of Borneo and China or Australia and Indonesia, recognized by the present study, or be influenced by extinct landforms such as the Oligocene island chain north of Australia. Even critics of the use of vicariant calibrations concede that their use can be appropriate in certain cases such as the evolution of freshwater taxa of low vagility . Nevertheless, some assumptions do remain in the present study, such as the assumption that the ecological habit of the Triculinae has only changed once in their history (from marshland to spring or primary stream dwelling), that the freshwater Pomatiopsinae could not cross open seas (whereas their marine ancestors presumably did in the Cretaceous) and that taxa dispersed always and whenever a route became open to them. In view of these considerations, as with all phylogeographical reconstructions, it is necessary to consider the hypotheses presented in this study as possible scenarios rather than definitive histories. Such histories provide frameworks for further enquiry and testing.
The use of molecular clock rates as calibration data, rather than the fossil record, has also been criticised [72, 73]. Numerous examples have been cited where the use of molecular clock rates for fishes led to estimated divergence dates either pre-dated the event attributed to their isolation [76–78] or the existence of their ancestral habitat [79, 80]. Such criticisms are now rather dated, and all the papers just cited used either a non-parametric method or a penalized likelihood approach (e.g., as used in the software r8S ), as was common at the time; the present study uses a Bayesian approach, which not only allows calibration dates to be incorporated as statistical distributions (reflecting the uncertainty of fossil dates or the gradual processes of continental breakup), but also gives more reliable results where a small number of loci are used  (see Phylogenetic methods: parameters and model priors). Further, most of the examples of problematic analyses cited in critiques of phylogenetic dating [72, 73] are of papers which used either clock rates, or fossils or palaeogeographical events, but not all three sources (as in the present study). Nevertheless, there are important caveats for the use of molecular clock rates. For example, it is highly problematic to use a clock rate from a study which employed a calibration date, e.g., a dated vicariance event, in a subsequent study which aims to date same vicariance event as used in the earlier study . Consequently, in this study care has been taken to choose calibration sources which are as independent as possible from the hypotheses being tested, and to give full details of the priors (calibrations) being used and, where secondary calibrations are used (e.g., clock rates from other studies), of the calibrations used in the primary studies. Further details of calibration points and sources of error in using fossil calibrations are given in Section Hypotheses testing.
The PCR primers used in this study amplified a 598 bp stretch of legible DNA sequence at the cox1 locus and 541 bp at the rrnL locus. The haplotype diversity was 1.000 for the cox1 data but only 0.999 for the rrnL data (the DNA sequences observed at this locus for Manningiella conica and Manningiella polita were identical). The combined data set was 1139 bp (1061 bp excluding sites with gaps/missing data) and comprised of 1177 characters including indels (coded by SIC) and alignment gaps/missing data.
Tests for deviations from neutrality were not significant: cox1 Fu* Li’s D* -0.37712 (P > 0.10), F* -0.73853 (P > 0.10); rrnL D* -1.46885 (P > 0.10), F* -1.82856 (P > 0.10), all calculations were based on total number of mutations and the outgroup was specified. Tajima’s D: cox1 -1.39619 (P > 0.10); rrnL −1.63067 (P > 0.05). The test of Xia et al.  did not indicate significant saturation in the sequence data (Iss 0.295, ISS.cSym 0.752, P < 0.0001, a lack of statistical significance here would imply a poor phylogenetic signal). Also, plots of transitions and transversions against distance suggested negligible levels of saturation at both loci.
The ML tree was also estimated for the data with the indels treated as missing data (not coded) (Additional file 1: Figure S1), for comparison with Figure 4. The trees differ in that the Amnicolinae are not monophyletic, Hubendickia spiralis is not basal to the other two Lao Julliennini, and the two Pachydrobia species are in slightly different locations on the tree. As the clade support values were similarly low and the Amnicolinae were polyphyletic, the inclusion of the indels was considered to be non-detrimental to the estimation of the phylogeny of the taxa in the present study.
Bayesian without clock model
The harmonic mean of the marginal likelihoods for the Beast runs was −11231.23; ASDSF was 0.001015 and PSRF values for all parameters were 1.000. Clade credibility values ranged from 32 to 100 percent (only 8 nodes had values < 50%). The topology of the tree representing the phylogeny estimated by MrBayes differed from that estimated by Beast in the following aspects. In the MrBayes tree the Jullieniini are basal at the root of the Pomatiopsidae (not derived from the same ancestor as the West China Triculini) and the Pomatiopsinae are found in the same clade as all Triculinae except for those of West China. In addition the two Oncomelania species are not monophyletic in the MrBayes phylogeny and the Gammatricula clade is not well resolved. The topology of the MrBayes tree is identical to that of the best-scoring tree found by RAxML. The Beast tree is described below.
Bayesian with clock model and date priors
Palaeogeographical/fossil record (FR) calibration points used as priors molecular dating
Prior code/taxa diverging
Geological events involved/biogeography
Root/Lithoglyphidae diverge from other hydrobioids
Rissooidea in Carboniferous – Hydrobiidae widespread by upper Jurassic 
237 ± 34.7
Core group/Proto-Pomatiopsidae and Amnicolidae
161 - 83
122 ± 19.9
Malaya/Robertsiella diverges from Sundaland Neotricula clade
2.0 ± 1.98
2.0 ± 1.01
Khorat/β and γ strains of N. aperta diverge in N E Thailand
Volcanism on eastern margin of Khorat Basin and river flow reversal 
0.8 ± 0.10
Parsimony method with POY
POY5 was used to implement a parsimony method using dynamic homology characters with direct optimization (i.e., no a priori multiple sequence alignment was required). In the initial sensitivity analyses the maximum MRI value (0.97990256) was obtained for a scheme with CP112 and rrnL partitions and cost weightings 6481; however, this was little different from the simpler scheme of cox1 and rrnL as one single partition and cost scheme of 6411 (MRI 0.96813765). Consequently, the data were used unpartitioned, Ts and Tv costs were both 1, and gap cost was 64. The MP analysis was found to be much more sensitive to cost weighting than to partitioning scheme.
The estimated phylogeny failed to recover any of the well recognized hydrobioid taxonomic groups (Additional file 3: Figure S3). Neither the Amnicolidae, Pomatiopsidae, Pomatiopsinae, Triculinae nor any of the three tribes of Triculinae appear monophyletic on the tree. The results of the POY analysis failed to shed any light on the evolutionary relationships of these taxa.
Shimodaira-hasegawa test (SH-test) and topological method tests
Constraint trees were drawn to represent the two main hypotheses regarding the evolution of the Triculinae, the Tibet origin with a West to East radiation and the Hunan origin with an East to West radiation (Figure 3). The log likelihood of the best-scoring tree found with a heuristic unconstrained search with PAUP* was −11378.017 and those for the same search constrained by the Tibet hypothesis and the Hunan hypotheses were −11449.67498 and −11447.08823, respectively. These findings can be used in LRTs as a crude comparison of the fit of the hypotheses to the data. In both cases the likelihood under the hypothesis was significantly lower than that of the unconstrained search (P < 0.0001); however, there was no significant difference when the likelihoods under the two models were compared with each other (P = 0.9999). The SH-test suggested that the constrained tree for the Hunan hypothesis was significantly worse than the unconstrained tree (P = 0.0056), as were 5/14 of the MP trees found by PAUP*. For the Tibet hypotheses, again the constrained tree was significantly worse than the unconstrained tree (P = 0.002409) as were the same five MP trees.
The topology based test (using slightly modified versions of the hypothesis trees, see Hypotheses testing in Methods Section) showed that 0/200100 trees in the posterior distribution were compatible with the Tibet hypothesis (P = 0.0000) and 35/200100 with the Hunan hypothesis (P = 0.0002). Consequently, all of the tree based hypothesis tests indicate that the Hunan is a slightly better fit to the data, although this difference is not significant and neither hypothesis agrees well with the data (both show a significantly poorer fit than the unconstrained tree).
Estimation of divergence times and substitution rates with Beast
Results of a Bayesian estimation of nucleotide substitution rates
Mean ± S.E.
95% HPD (lower – upper)
-ln (likelihood) Deep divergences
−10921.7241 ± 0.201
−10921.3095 – –10902.1831
-ln (likelihood) Most recent divergences
−10924.1375 ± 0.1375
−10944.6406 – –10905.5067
Mean substitution rate (cox1) First codons
0.70605 ± 0.0031356
0.37368 – 1.1606
Mean substitution rate (cox1) First codons
0.14204 ± 0.0014602
0.057033 – 0.2659
Mean substitution rate (cox1) Second codons
0.087271 ± 0.00046307
0.036775 - 0.14957
Mean substitution rate (cox1) Second codons
0.017227 ± 0.00016606
0.0055299 – 0.033703
Mean substitution rate (cox1) Third codons
25.54 ± 0.06864
17.16 – 34.53
Mean substitution rate (cox1) Third codons
7.2543 ± 0.079013
2.7901 – 13.05
Mean substitution rate (cox1) All Sites
2.1253 ± 0.033205
0.8405 – 3.9242
Mean substitution rate (rrnL)
1.2869 ± 0.0038937
0.79066 – 1.854
Mean substitution rate (rrnL)
0.27031 ± 0.0021339
0.12723 – 0.47463
Mean substitution rate indels
0.73603 ± 0.0033437
0.37765 – 1.094
Mean substitution rate indels
0.14923 ± 0.0012716
0.06371 – 0.26997
Results of a Bayesian estimation of divergence times
Parameter [node on tree]
Mean ± S.E.
95% HPD (lower – upper)
tmrca (ingroup) 
114.533 ± 0.6187
52.5946 – 176.4489
≤ H1, T1
tmrca (Pomatiopsidae) []
110.760 ± 0.5864
51.8809 – 170.102
tmrca (Pomatiopsinae) []
59.0213 ± 0.5583
16.9746 – 111.8963
tmrca ( Oncomelania ) []
22.4588 ± 0.2018
15.4905 – 34.7171
tmrca (Sunda) 
7.925 ± 0.048796
3.9652 – 13.1893
tmrca (Thai) 
5.963 ± 0.020944
3.621 – 8.469
≥ H11, ≤T7
tmrca (North_Sunda) 
5.1588 ± 0.01804
3.1019 – 7.1554
H8, T5, T7
tmrca (Yunnan) 
4.1509 ± 0.01424
2.4844 – 5.9712
tmrca (Malaya) 
3.5555 ± 0.01417
2.3954 – 5.0629
≥ H13, ≤ T7
tmrca (Mekong) 
2.3658 ± 0.00790
1.4563 – 3.3614
tmrca (Khorat) 
0.8679 ± 0.00140
0.6755 – 1.0775
tmrca (Jullieniini) 
0.7847 ± 0.00363
0.3492 – 1.3228
tmrca (Lao) 
0.3561 ± 0.00154
0.181 – 0.569
Origin of the pomatiopsidae
Table 4 indicates that the Pomatiopsidae diverged from the Amnicolidae around 111 Ma; this is consistent with the break up of the minor Gondwanan continent (which included Australasia) around 132 Ma and rafting of pomatiopsine-like proto-Pomatiopsidae on these Gondwanan fragments, with subsequent island hopping northwards (and some (tectonic) rafting) and radiation of marine pomatiopsines along a Borneo-Philippines tract and also to the southern continents (South Africa and South America). The estimated date is less consistent with a divergence prior to the break up of Pangaea and rafting of proto-pomatiopsids (resembling Triculinae) and Asian amnicolids on the Indian craton (around 180 Ma); however, the upper bound of the HPD for this estimate was 170.102 Ma, which is close to 180 Ma. Consequently, this date alone is not sufficient to rule out any hypothesis.
The common ancestor of the Triculinae and Pomatiopsinae was estimated to have occurred around 59 Ma (Table 4). This taxon may have inhabited Eastern Gondwana (Australia) after the loss of India/Madagascar, with later island hopping from Borneo and Sulawesi along the extensive island complex existing in the region during the Eocene-Oligocene. These snails could have then given rise to the antecedent Triculinae when they reached China, and Pomatiopsinae on the Japanese islands. The Tibet hypothesis would require the Pomatiopsinae to be derived from the Triculinae or the former to have crossed the uplifting Himalaya (something which for the extant Pomatiopsinae show little adaptation); the southern margin of the Qinghai-Tibet Plateau achieved present day levels 22–15 Ma. In addition, extensive and semi-continuous volcanism, 45–26 Ma and around 20 Ma , in the Himalayas makes trans-Himalayan dispersal tracts even more unlikely. Volcanism along the Central Volcanic Arc also would have restricted migration across Myanmar even until the Pleistocene .
Radiation of the pomatiopsinae
The TMRCA estimates next suggested very little divergence among Pomatiopsidae (at least among lineages with extant representatives) until around 22 Ma, with high rates of cladogenesis 8 – 2 Ma (Table 4). During the early Miocene the Japan Sea opened and the Japanese islands separated from China. Many of the Japanese islands underwent a cycle of emergence, submergence and re-emergence (18–16 Ma) . The regional climate at that time also became warm and humid [52, 96]. These events could have driven the divergence of amphibious and terrestrial Pomatiopsinae in Japan. Figure 6 and Table 4 suggest that this occurred around 22 Ma, with Oncomelania shown to have appeared prior to the other Japanese Pomatiopsinae; this would have allowed ancestors of Chinese Oncomelania (O. hupensis) to have colonized mainland China before land connections were broken (around 15 Ma ). O. minima and the terrestrial pomatiopsines are not found in China. The Pomatiopsinae are also absent from Korea (which is unexpected if they colonized Japan from China), possibly because deep marine basins have isolated the Korean peninsula since the Palaeogene .
Divergence of the triculinae
The divergence of the Pachydrobiini from a clade comprising the Triculini and Jullieniini is estimated to have occurred around 8 Ma (Table 4, Sunda). This date coincides with the attainment of present day elevations across the whole Qinghai-Tibet Plateau  and global changes in climate, with the Himalaya obstructing global air circulation and the predominance of a monsoon climate in the region . Most Triculini are restricted to primary streams (usually from springs) in highland areas, whereas almost all Pachydrobiini live in larger streams and rivers. The isolation of spring-dwelling highland taxa in uplifting regions of Southwest China and the more arid climate could have triggered the divergence of the proto-Triculinae, with the invasion of larger rivers and streams. From this point onwards a rapid acceleration in triculine radiation and cladogenesis is observed (Figure 6); the acceleration among the Triculini could have been driven by deepening of the valleys in which they live (i.e., increased isolation) and that among the Pachydrobiini by the greater variety of new habitats and area covered following the invasion of rivers. The next divergence event is that of the Tricula of northern Thailand (including T. hortensis of Sichuan) from a clade comprised of Jullieniini and the Yunnan Triculini; this event is dated at around 6 Ma (Thai, Table 4). Table 1 links this event (under the Hunan hypothesis) to Late-Pliocene (3.4 Ma) block faulting which may have induced several course changes along the proto-Mekong river  and the break up of connexions between Yunnan/Guangxi China and northern Thailand via the extended Mekong-Ping river, and flow reversals  in Northwest Sundaland. The estimated date of 6 Ma, however, is more closely aligned with the second major uplift of the Himalaya (7.5-6.5 Ma, ), which would have also led to uplift along the Shan formation and other ranges acting as barriers between China and Sundaland. Uplift of northern Thailand would also have accelerated and would have isolated taxa there from those radiating in the lower Mekong river. The next divergence, dated at around 5 Ma, concerns the Pachydrobiini, with all other Pachydrobiini (mostly Sundaland and Mekong river taxa) arising from a lineage now represented by Gammatricula (of Guangxi/Fujian China) (North Sunda, Table 4). Table 1 (Hunan hypothesis) associates this event with the loss of the Red river corridor between southern China (Hunan/Guangxi) and Northeast Sundaland  at around 4.4 Ma .
The sundaland radiation
The next major event dated was the divergence of Manningiella and Neotricula in the lower Mekong Basin, which was estimated at around 2.4 Ma (Table 4, Mekong). This event coincides with the end of a global cooling period (2.7-2.5 Ma) , and a significant Asian monsoon intensification occurred prior to the major glaciation at around 2.7 Ma, which triggered an extremely arid 1.4 Ma interval . Such climatic events could have driven the divergence of proto-Mekong river Triculinae as they re-encountered more shallow water, stream-like habitats. Finally, the divergence of the beta and gamma strains of Neotricula aperta is dated at around 0.87 Ma (Khorat, Table 4); this date agrees well with the predictions of earlier studies (Table 1, H15). The divergence of N. aperta was considered to be driven by events linked to the final surge of the Himalayan uplift , such as the tilting of the Khorat Basin, reversal of the flow of the Mul-river, and volcanism in southern Laos which temporarily separated the lower Mul river from the Mekong .
Erhaiini and jullieniini
The phylogenies in Figures 4 and 5 both indicate similar phylogenetic affinities, which help to clarify the taxonomy of the Asian “hydrobioids”. The taxon resembling Erhaia, which was collected in Guangxi Province, China, is seen to form an amnicolid clade together with Akiyoshia of Japan and the European Bythinella austriaca; this further supports the decision of other authors [17, 18] to transfer all Pseudobythinellini/Erhaiini to the Amnicolidae. The position of Erhaia was uncertain because although this taxon possesses a Hydrobia-like central tooth (as in Amnicolidae) it shows closure of the ventral channel to form a spermathecal-duct like structure (as in Pomatiopsidae) and the penial hold-fast is absent . Erhaia is known from Fujian, Hubei, Hunan, Sichuan and to a lesser extent, Yunnan, China; the finding of Erhaia in Guangxi, though expected (Guangxi lies between Fujian and Yunnan), extends the known range of this genus.
The position of Jullienia rolfbrandti at the base of the pomatiopsid clade (Figures 4, 5 and6) is unexpected. In the tree inferred by Beast (Figures 5 and 6) all other Jullieniini form a distinct clade within the Triculinae and arise from the base of the West China Triculini lineage. In the phylogeny estimated by RAxML, the Jullieniini do occur on (monophyletic) adjacent branches with Jullienia the most basal, and all Jullieniini being clustered at the base of the pomatiopsid phylogeny. Consequently, the unexpected position of Jullienia on the Beast tree may be a failure of the analysis to land upon a correctly resolved tree; however, it is equally likely that the Jullieniini is an artificial group. The Jullieniini were envisaged as having diverged from the lineage leading to the Pachydrobiini, Triculini and Lacunopsis during the Miocene; they were considered to be united by a modification of the central tooth and adaptations of the female reproductive system (e.g., the oviduct circle complex), and another derived character state of the Triculinae, that of large shell size . The misidentification of specimens or mis-labelling of DNA extracts in the laboratory is unlikely because there were only samples of Triculinae in the laboratory at that time and no other genus from the collecting site has a shell, radula and gonad resembling a Jullienia sp. that is not Triculinae.
The Triculini and their relationships with other triculine taxa
Turning next to the Triculini, it can be seen that Tricula is monophyletic in neither the Bayesian (Beast) nor ML (RAxML) phylogenies (Figures 4 and 6). The northern Thai Tricula and Tricula hortensis are found at the root of the Triculini. This Thai clade forms a bifurcation with a clade comprising the Jullieniini (at the root) and the West China (Yunnan) Triculini, which include Tricula bambooensis, T. hudiequanensis, T. ludongbini and Delavaya, with Lacunopsis at the base. The remaining Tricula species, T. xiaolongmenensis, is found at the root of the Triculini + Jullieniini clade. Even allowing the possibility of a small degree of systematic error affecting the exact distribution among clades, the genus Tricula is clearly not monophyletic (the credibility of the aforementioned clade is 0.87, which is greater than the 0.81 for the monophyly of the Triculinae, Figure 5). Earlier morphological studies placed the Chinese Tricula in a clade closely associated with Delavaya and Lacunopsis, with the remaining Tricula (including the type species T. montana) in a clade at the base of the Triculini and closest (among Triculini) to Hubendickia. The present results appear to support this division which was recognized in a series of remarkably detailed earlier morphological studies . Although it is not possible to obtain a sample of Tricula montana (due to India’s bio-resources protection policy), the basal position of T. montana in the cladogram estimated by Davis et al., , suggests that the northern Thai species (including T. bollingi), and possibly also T. hortensis, and T. xiaolongmenensis are Tricula sensu stricto, whereas, the Yunnan species should not be regarded as Tricula. The Yunnan “Tricula” are congeners of Delavaya, whether or not they can all be classed as Delavaya spp. would depend on a detailed re-examination of the anatomy of T. hudiequanensis. In Figure 4, the ML phylogeny, T. hortensis is found at the base of a poorly supported (8%) clade including the Pachydrobiini and Pomatiopsinae! This is most likely a failure of phylogenetic estimation by direct ML; however, it is worth considering that previous authors noted the difficulties of determining the affinities of this taxon – which was described as difficult to place within the Pomatiopsidae , or as showing character states common to both Triculinae and Pomatiopsinae. The inclusion of Lacunopsis at the root of the Yunnan Triculini clade agrees with the conclusions of morphological studies that Lacunopsis arose from the Triculini and is close to Delavaya. Finally, the relationship needs to be considered between T. bollingi and a species of Tricula found in a spring only 11.2 km distant. The new taxon is most likely a species distinct from T. bollingi because, when compared with other sister species of Tricula (e.g., T. bambooensis and T. ludongbini) the divergence between this new taxon and T. bollingi is much greater. The two taxa also appear to have diverged around 5 Ma; this suggests that the two have very different histories, for example, the new taxon may have originated locally, whereas T. bollingi could have been introduced from China or Myanmar (for example on religious artifacts transported by monks – the collecting site for T. bollingi was downstream of a large temple). Tricula bollingi is also known from Yunnan, China . Additional file 4: Table S1 lists taxa highlighted for revision and suggested changes.
Considering next the relationships of the remaining Jullieniini, the direct ML trees (PAUP* Additional file 1: Figure S1 and RaXML Figure 4) and the Bayesian estimated phylogenies all show these taxa as arising at the root of the clade leading to the Triculini; however, in all these trees these jullieniinines do not occur in a cohesive monophyletic clade and the Triculini are at the base of the Pomatiopsidae, whereas in the Beast estimated tree, the Triculini form a bifurcation with the Pachydrobiini in a monophyletic Triculinae clade, and these Jullieniini form a monophyletic clade at the base of the West China Triculini (Figure 5). Hubendickia spiralis lies at the root of this clade, with Hubendickia schuetti and Paraprososthenia levayi as sister taxa; this is well supported (credibilities 0.99 and 1.0, respectively). The implication is that Hubendickia schuetti is Paraprososthenia. Indeed, this taxon was first described as Paraprososthenia and was later transferred to Hubendickia because of its ovate-pyriform aperture. The tribe Jullieniini is characterized by a specialized central tooth (excluding some Hubendickia) and elongation (Hubendickia) to extreme elongation (Paraprososthenia, Jullienia, Hydrorissoia, and Pachydrobiella) of the seminal receptacle. Increasing degrees of elongation, reflexion or coiling of different parts of the vas deferens are also seen as a characteristic of increasingly derived jullieniinine genera, the extreme being seen in Jullienia and Hydrorissoia. These adaptations of the reproductive system have been considered as responses to mating in the early flood period of the lower Mekong river, which is the time when these taxa complete their life-cycles [10, 43, 103]. In view of this it is possible that the character states serving as synapomorphies uniting the Jullieniini are simply due to convergence; thus, the apparent affinity of schuetti with species of Hubendickia may simply reflect a similar microhabitat and rapid radiation in the evolving Pleistocene lower Mekong river (this divergence was dated to around 0.8 Ma, Jullieniini, Table 4). Convergence may also explain the failure of Jullienia rolfbrandti to cluster with the other Jullieniini in this study.
Taxonomy of the pachydrobiini
The Pachydrobiini show a near isochronous burst of divergence events with that seen among the Triculini, suggesting that both lineages responded to the same event (probably the second major uplift of the Himalaya and Sibamasu orogeny, see above). In the phylogenies estimated by ML and Bayesian approaches Gammatricula is seen to occur at the root of the Pachydrobiini. In the Beast (and PAUP*) phylogeny the two Gammatricula spp. form a monophyletic clade, in the RAxML tree they co-occur with Jinghongia and the Vietnamese taxon resembling Pachydrobia. The observation is interesting because Gammatricula was described as a genus showing character states that were overall derived among Pachydrobiini . The Beast phylogeny shows that the taxon from northern Vietnam (Pachydrobia sp.) lies at the root of the Pachydrobiini, just above the Gammatricula clade (or inside it according to ML), and therefore is unlikely to be Pachydrobia. Snails described as Vietricula species have been described from Cao Bang in northern Vietnam , from a site approximately 200 km East of the collection locality for the taxon collected in the present study. The descriptions in the 2012 study appear to be based on shell character states and so the possibility that these snails were derived from early Pliocene Chinese Gammatricula in the nearby provinces of Yunnan and Guangxi cannot be excluded.
Jinhongia, which is found in the middle Mekong river in Yunnan, China, is found at the base of the Sundaland Pachydrobiini clade; this sugests that around 4 Ma taxa from Guangxi/Yunnan entered Sundaland from the North and began a rapid radiation southwards in the lower Mekong river. The southwards radiation is in accordance with the Tibet hypothesis as originally proposed , although the time scale is different. The next extant taxon in this southwards radiation is Neotricula burchi. Regardless of the phylogenetic method, the present data indicated that Neotricula was not monophyletic because species of Manningiella, Pachydrobia and Robertsiella are found with the Neotricula clade. The type for the genus Neotricula is N. aperta (assignation ), so that Neotricula burchi should be placed in a new genus (or in Jinhongia – additional samples of Robertsiella are needed to decide this). The three species of Manningiella sampled form a well supported (Figure 5, 0.98) clade with N. aperta; therefore it is possible that these are in fact species of Neotricula. Additional, sampling (of loci and species) and morphological work is required on Manningiella and the two taxa from northern Thailand and central Laos (here denoted as Neotricula sp. and Neotricula sp. 2) to determine if these taxa are indeed all Neotricula spp. The question is one of particular importance because Neotricula aperta acts as intermediate host for Schistosoma mekongi and it is clearly very closely related to taxa described as Manningiella. It should be mentioned that most, perhaps all, species originally assigned to Manningiella by  have been transferred to other genera (even including non-pachydrobiinine genera, such as Hubendickia), e.g., M. conica and M. expansa were assigned to a new genus, Halewisia, on the observation that the spermathecal duct runs directly into the duct of the bursa copulatrix in an anterior position (in other Triculini it enters in a posterior position), the sperm duct does not form a common sperm duct with the duct of the seminal receptacle (so resembling Pachydrobia and Tricula), and shells are smooth, ovate-conic, with expanded inner lip . In the present study we used the name Manningiella because the detailed descriptions of taxa such as polita and velimirovici refer to Manningiella, because Halewisia is not in common usage, and because there are no detailed accounts of the anatomies of several species assigned to the nominal taxon Halewisia. Nevertheless, it should be noted that the species polita and conica listed here as Manningiella have been referred to Halewisia in some works .
The present study has demonstrated that the idea of an Australasian origin for the Pomatiopsidae, followed by an East to West radiation is a plausible alternative to a hypothesis of rafting of Pomatiopsidae on the Indian craton, with divergence in the Himalaya and a westwards colonization of Sundaland and China in the early Miocence. The findings provide support for a Hunan origin for the Triculinae, no later than around 8 Ma. The radiation the Pomatiopsinae is seen to begin some 14.5 Ma prior to that of the Triculinae; this suggests that the former diverged first (possibly in Japan) and the Triculinae diverged in China much later and probably from pomatiopsine stock. The scenario that both sub-families diverged in Tibet as they colonized the upper Yangtze and Mekong rivers, respectively, is less probable in the light of the current data. The results confirm the rapid radiation and burst of cladogenesis in the Triculinae, as postulated by previous authors [15, 21, 64], and suggest that the acceleration began around 6 Ma. The molecular dating indicates that the radiation of these snails was driven first by the uplift of the Himalaya and onset of a monsoon system, and subsequently by late-pliocene global warming.
The analyses have helped confirm the status of Erhaia (formerly tribe Erhaiini) as Amnicolidae rather than Pomatiopsinae: Pseudobythinellini as previously thought and have extended the known range of Erhaia in China. The phylogenetic estimations demonstrated that the genera Tricula and Neotricula are not monophyletic and that the Jullieniini maybe an artificial taxon based on convergent responses of the hydrobioid female reproductive anatomy to life in the faster and more extreme flow regime of the lower Mekong river. The study showed that several species assigned to Halewisia/Manningiella are very closely related to Neotricula aperta and this is important for the transmission of Mekong schistosomiasis because N. aperta is the only known intermediate host for Schistosoma mekongi. The intermediate host of Schistosoma malayensis, Robertsiella, was also shown to have a possibly independent evolutionary history from Neotricula aperta, with a more direct entry into Cambodia via now extinct and submarine river systems linking Hunan-Vietnam and Cambodia; this also has implications for the radiation of Schistosoma. The work also supported the conclusions of George Davis, based on intricate and remarkably detailed anatomical studies of the Triculinae, that there were divisions within the genus Tricula and of the origin of the Lacunopsini from the stem of the Triculini clade.
The present study has also raised several new questions which require further study. These include additional work on Pachydrobia, Robertsiella and Jinghongia with the inclusion of more species (except Jinghongia) and more loci, so that the generic status of Neotricula burchi can be elucidated. Surveys in Fujian and Guangxi, China, are required to detect further species of Gammatricula or Vietricula, to determine the genera of these taxa and their phylogeographical history. Much more anatomical work is required for Halewisia, to determine which species should be assigned to this genus and the extent of affinities with Neotricula aperta. Renewed parasitological investigations are also required in this context, past investigations suggested that miracidia of S. mekongi can locate and penetrate Manningiella expansa but that the snails died before the sporocysts would have become patent . In later experiments 12% of Manningiella conica exposed to S. mekongi were observed to shed cercariae, which is a rate very close to that of similarly exposed N. aperta. Aside from this, the transmission of Schistosoma appears limited to 3 clades, scattered across the phylogeny, with a most recent common ancestor around the Paleocene; this suggests that the evolutionary relationship between these snails and Schistosoma is more like loose phylogenetic tracking, rather than co-evolutionary. The finding has implications for future studies of these host-parasite relationships and prediction of the long-term impacts of snail control in regional schistosomiasis. The present study included 16 genera of Pomatiopsidae but this represents only 4 of 9 Pomatiopsinae, 3 of 5 Triculini, 3 of 9 Jullieniini, and 6 of 9 Pachydrobiini genera (this excludes those genera of uncertain affinity or data deficient, e.g. Saduniella, Vietricula). Clearly there is much scope for future studies with more extensive taxon sampling and, if taxa can be re-discovered (where necessary) and re-sampled, the inclusion of additional loci in the phylogenetic analyses. It is hoped that the present study has provided the first step in solving the taxonomic confusion and lack of data surrounding this at risk, highly diverse and medically important group of freshwater snails.
Taxonomic and collection details of the taxa included in the phylogenetic reconstructions
GenBank: AB611823, GenBank: AB611822
GenBank: AF445333, GenBank: AF445342
GenBank: KC832701, GenBank: KC832722
Müggelsee (near Berlin)
GenBank: AF445334, GenBank: AF445344
Hydrobia acuta/Hydrobia sp
Ban Xieng-Wang Village
GenBank: KC832688, GenBank: KC832709
Ban Xieng-Wang Village
GenBank: KC832689, GenBank: KC832710
Ban Xieng-Wang Village
GenBank: KC832697, GenBank: KC832718
GenBank: KC832698, GenBank: KC832719
GenBank: KC832694, GenBank: KC832715
GenBank: KC832695, GenBank: KC832716
GenBank: KC832687, GenBank: KC832708
Spree (near Berlin)
GenBank: AF445332, GenBank: AF445341
GenBank: AF213342, GenBank: AF212896
GenBank: KC832703, GenBank: KC832724
GenBank: KC832707, GenBank: KC832728
GenBank: AF531542, GenBank: AF531544
GenBank: KC832704, GenBank: KC832725
GenBank: KC832706, GenBank: KC832727
GenBank: AF531541, GenBank: AF531556
GenBank: KC832700, GenBank: KC832721
Nam Sai River
Van Ban District
GenBank: KC832690, GenBank: KC832711
GenBank: AF531550, GenBank: AF531548
Blanfordia j. japonica
GenBank: AB611727, GenBank: AB611726
GenBank: AB611747, GenBank: AB611746
Fukuia kurodai ooyagii
GenBank: AB611783, GenBank: AB611782
Oncomelania h. robertsoni*
GenBank: AB531547, GenBank: AB531545
GenBank: AB611795, GenBank: AB611794
GenBank: KC832692, GenBank: KC832713
GenBank: KC832705, GenBank: KC832726
Hei Qiong Duan
GenBank: KC832699, GenBank: KC832720
GenBank: AF531553, GenBank: AF531551
Han Wang Zhen
GenBank: AB531539, GenBank: AB531554
GenBank: KC832691, GenBank: KC832712
Hei Long Tan
GenBank: KC832696, GenBank: KC832717
Mae Fang River
GenBank: KC832702, GenBank: KC832723
GenBank: KC832693, GenBank: KC832714
DNA amplification and sequencing
The shells were gently cracked and the bodies removed. The gut and digestive gland were discarded and DNA extracted from the remainder by standard methods . A section of the coding region of the mitochondrial (mt) cytochrome-c oxidase subunit I (CO1) gene (cox1) was amplified by PCR with the HCO-2198 and LCO-1490 primer pair developed by Folmer et al. (1994) , and part of the gene coding for the mt 16S rRNA (rrnL) using the primers of . The cox1 gene was chosen because it is a popular locus in barcoding projects  and showed adequate levels of interspecific variation in past studies of triculine phylogenetics (33% of sites polymorphic across Pomatiopsidae ). The mt 16S locus had previously been used successfully to investigate within family relationships among lymnaeids  and Pomatiopsidae [3, 18, 20]. The rrnL data were checked for artefacts such as chimeras using Pintail 1.1 , with the sequences for Tricula bollingi and Neotricula aperta used as subjects and T. hortensis and N. burchi used to verify integrity of the subjects, respectively (these taxa were chosen because they had been sequenced using several pairs of PCR and sequencing primers and had been used in phylogenetic studies of these taxa for many years without any anomalous results). The two loci were also chosen for compatibility with sequence data available for taxa which are no longer known to science as extant (i.e., all known populations are extinct, as in the case of many snail populations around Dianchi Lake) but for which our laboratory, or the GenBank, had cox1 and rrnL sequences on record.
Total genomic DNA was used as a template for PCR amplification on a Progene thermal cycler (MWG) employing standard PCR conditions . Unincorporated primers and nucleotides were removed from PCR products using the QIAQuick PCR purification kit (QIAGEN). Sequences were determined bidirectionally, directly from the products by thermal-cycle-sequencing using Big Dye fluorescent dye terminators and an ABI 377 automated sequencer (Perkin-Elmer), following procedures recommended by the manufacturers. Sequences were assembled and aligned using Sequencher (version 3.1 Gene Codes Corp. Ann Arbor, Michigan). DNA sequences for both strands were aligned and compared to verify accuracy. Controls without DNA template were included in all PCR runs to exclude any cross-over contamination. A multiple sequence alignment was generated for all taxa for each locus and the cox1 and rrnL sequences were then concatenated to form a combined data set. Potential outgroup sequences were obtained from the GenBank; these included two taxa of the Bithyniidae/Hydrobiidae and one Lithoglyphinae (Bithynia tentaculata, Hydrobia sp. and Lithoglyphus naticoides, respectively, see Table 5). In addition to one taxon sampled in this study (a previously unknown species of Erhaia from Guangxi, China), sequence data for two Amnicolidae were also taken from the GenBank. A variety of outgroup taxa was made available to the analysis in order to help detect any long branch attraction (LBA), which may affect the estimate of phylogeny, through the commonly applied method of taxon addition/removal and varying outgroup [123, 124].
Initial handling of data and selection of partitioning scheme and substitution models
Tajima’s test (based on the total number of mutations) and Fu and Li’s test for neutrality was performed using DNAsp (v. 5.10.01) . The data were tested for substitution saturation using plots of the numbers of transitions (Ts) and transversions (Tv) against the F81 genetic distance . The indications of these plots were further evaluated using the entropy-based test  as found in the DAMBE (v. 4.5.20) software package , which provides a statistical test for saturation. For all phylogenetic reconstructions (except for those using POY, which employed dynamic homology) indels were coded using Simple Indel Coding .
The incongruence length-difference (ILD) test  as implemented in PAUP* (v. 4.0b10; ) with 5000 replicates, was used to test for homogeneity between the cox1 and rrnL data partitions; the test was applied to informative sites only . A suitable substitution model was selected for each data partition using an hierarchical test of alternative models as implemented in jModelTest (v. 2.1.2) [132, 133], with default settings. Models were compared by AICC. The starting parameters of the chosen models were then “optimized” using a Maximum Likelihood (ML) method with the Brent Powell algorithm in the phylogenetics software suite P4 (v. 0.90) . Molecular evolution of indels was modeled by a F81 model, which is a monotone model fitting the similarly monotone substitution of a binary character evolving along a phylogeny ; however, in the analyses using Beast (v. 1.7.5) , a stochastic Dollo model  was also evaluated. It was considered that the Stochastic Dollo might provide a more realistic modelling of the indel process for our data if each indel is a unique, independent, event in the phylogeny. The Stochastic Dollo model here assumes that an indel is gained and then differentially lost in descendent lineages. No backward mutation is assumed. Each new indel is assumed independent of any pre-existing indel. The model has two rate parameters, an indel may arise at Poisson rate λ and is lost at constant rate μ . The model was implemented both with and without correction for ascertainment bias (which may occur because the invariant sites have been excluded and are thus unobservable). The correction involves re-scaling the finite-time transition probabilities .
Phylogenetic methods: parameters and model priors
Phylogenetic estimation used three different approaches; this strategy was adopted to look for resilience of the hypothesis to changes in method and associated assumptions. The rationale was that any clade that was represented in phylogenies found by all methods (and well supported) would be a reliable reconstruction of phylogenetic history for these taxa. In addition, the strategy increases the chance that artefacts peculiar to a specific class of methods will become apparent. For example, parsimony methods are perhaps more susceptible to the effects of long-branch attraction (LBA) than probabilistic methods [123, 139] (although ML methods can also show LBA problems). Similarly, Bayesian methods have been known to erroneously converge on Long-Tree (LT) solutions in cases where partitioned data are described by complex models with many parameters close to nonidentifiability [140, 141].
Likelihood based methods
A standard maximum likelihood method was applied using RAxML (v. 7.3.2) . RAxML was chosen for this purpose because of its ability to handle partitioned data, whilst at the same time perform rapid bootstrapping and make use of parallel computing resources. RAxML also speeds ML search times using “Lazy Subtree Rearrangements” for branch swapping, and thus allows more replicates/multiple runs within the same time frame. RAxML also uses a simulated annealing approach, which incorporates a cooling schedule and allows “backward steps” during the hill-climbing process.
Bayesian phylogenetic analyses were implemented separately through two different packages, which were MrBayes (v. 3.2.1)  and Beast (v. 1.7.5) . Bayesian approach (BA) was included in the analyses because these are not only faster in terms of computing time (for analyses with an equivalent level of confidence), but are also statistically superior to a method based on unintegrated likelihoods . For example, such methods do not assume approximate normality or large sample sizes as would general ML methods . BAs differ from direct ML methods in that the former consider the posterior probability of the model (with parameters) and tree after observing the data. The posterior probability of a hypothesis is proportional to the product of its prior probability and the probability of observing the dataset given the hypothesis (i.e., the likelihood of the hypothesis). Consequently, unlike direct ML, a BA allows incorporation of prior information about the phylogenetic process into the analysis. The analyses with Beast differed from those with MrBayes in that, with Beast, a molecular clock was applied to the tree and specified calibration point (date) priors were used in addition to substitution model parameter priors. Beast was also used to simultaneously estimate the divergence times for phylogenetic events relevant to phylogeographical hypotheses. The Bayesian approach implemented in Beast is currently considered superior to other approaches (e.g., non-parametric methods such as NPRS  or penalized likelihood methods , particularly for divergences with a low time depth, because it can allow for uncertainty in dates assigned to calibration points and does not require untested assumptions about the pattern of clock rate variation among lineages . The greatest benefit of using a Bayesian method for dating is that the specification of prior distributions can be used to ensure that the analysis realistically incorporates the uncertainty associated with the calibration points used . In the present study the indels were modelled by the TKF91 model (the default in Beast), which was found to perform better than the Stochastic Dollo. The settings particular to each phylogenetic method, and methods used in their selection, are described in detail in in Additional file 5 (Detailed description of methods). Details of priors used for dating in the Beast analyses are given below in the “Hypotheses testing” section.
Parsimony method using POY
POY (v. 5.0.0)  (parallelized build) was used to implement a Maximum Parsimony approach. The use of MP afforded an analysis free of the assumptions underlying ML methods, and the use of POY with its dynamic homology approach (where characters (transformation series) are inferred during the process of phylogenetic reconstruction) frees the analysis of any effects particular to the alignment inferred by Sequencher . The choice of weighting schemes and run settings used followed published procedures .
Topology based method
All non-Pomatiopsidae (including Jullienia rolfbrandti) were restricted to an unresolved (polytomic) clade at the root of the tree (after the outgroup)
Tricula xiaolongmenensis at root of the pomatiopsine clade
Tricula hortensis at the root of all Triculinae, except for the West China Triculinae clade
Gammatricula sp. and then (G. fujianensis, J. jinghongensis) at the root of all Sundaland Triculinae.
Shimodaira-hasegawa (SH) test
SH-tests  were performed as implemented in PAUP* with RELL optimization and 1000000 bootstrap replicates. The test compares between and among topologies from constrained and unconstrained searches; the set of 14 “best” trees from a MP analysis using PAUP* and the ML tree from an unconstrained PAUP* search were used as the set of trees for this comparison. The trees from Figure 3 (unedited) were used as the constraints for the two hypotheses under test. Although these two trees were written a priori, the use of the Kishino-Hasegawa test was not considered appropriate because the authors could not rule out the influence of earlier studies (some of which involved a sub-set of the present data) on the design of the trees. PAUP* was run with a CP123 and rrnL partitioning scheme, site specific rates, and the TPM1uf + I + G model with parameters as chosen by jModelTest for all sites combined. Gaps were treated as missing data.
Molecular dating method with beast
All Beast runs were performed with up to four date priors as calibration data. In reflection of the lack of certainty surrounding fossil dates, for example, divergence of taxa often pre-dates their first appearance in the fossil record , fossils maybe incorrectly assigned to extant taxa , and strata in which fossils formed may be eroded and subsequently overlain (leading to incorrect dating [155, 156]), a wide time interval was used for each prior. Fossil dates are also subject to revision as new “earliest” fossils are discovered . Molecular clock based calibrations often use clock rates set by fossil calibrations or palaeogeographical events (as in this study) and are, therefore, no more reliable than those based directly on fossils. Indeed molecular calibrations may suffer from the additional complications of rate variation among lineages or over evolutionary time . As detailed below, measures were taken in the present study to minimize these problems. In addition to specifying wide ranging date priors, the priors were allowed to vary over a normal distribution such that the prior range (in most cases) corresponded to a 95% confidence interval. The HPD traces for parameter estimates were examined following each run to ensure that the estimates were not being truncated or overly skewed by the prior. Consequently, the prior effectively acts as a guide for the analysis, with the final estimates representing the result of an optimal combination of substitution model parameters, data, phylogeny and priors (this is a strength of all Bayesian methods). Care was also taken to achieve independence of calibration prior and hypothesis being tested.
Beast was run with the priors Root and Core-Group (Deep divergence runs) or Root, Core-Group Malaya and Khorat (Recent divergence runs). The “Root” prior provides a probability distribution for the TMRCA of the ingroup and Lithoglyphopsis (the root). The fossil record suggests that the Lithoglyphidae diverged from the Bithyniidae and Amnicolidae (i.e., from other hydrobioids) in the Jurassic . In addition, fossil putative Rissooidea are known from the Carboniferous [157, 158] and fossil hydrobioids are widespread by the upper Jurassic . Consequently, the root of the phylogeny was dated as upper Carboniferous to mid-Jurassic (305–169 Ma). This time interval is also consistent with estimates based on a mt DNA clock (calibrated using minimum divergence dates from the littorinid fossil record) . The Core-Group prior refers to the TMRCA of all the taxa except Lithoglyphopsis and Hydrobia, and is effectively the divergence date for the proto-Pomatiopsidae and the Amnicolidae. Fossil Bithyniidae are known from the Purbeckian , and the assimineid-pomatiopsid-truncatellid radiation is also thought to be Purbeckian . The Purbeckian is roughly equivalent to the Berriasian . Consequently, the upper range for this divergence was taken to be upper Jurassic. Hydrobiid taxa are seen to be well diverged (diverse) by the Santonian-Campanian, consequently the mid-Cretaceous was taken as the latest date for this range (i.e., 161–83 Ma, see Table 2). The prior Malaya (for the divergence of Robertsiella from the Sundaland Neotricula clade) was obtained from the divergence time estimated by a mt DNA molecular clock for Schistosoma malayensis and S. mekongi; this clock was in turn calibrated using DNA substitution rates from the S. incognitum and S. mansoni divergence , which was based on hypotheses unrelated to those of the present study. Nevertheless, a normal prior (5.0 ± 1.0 Ma) was applied to the TMRCA of the S. sinensium group in the initial paper . The prior was based on a date for the second major Himalayan orogeny; however, this was referring to events isolating Central Asia from the Orient , i.e., the acceleration of Tian Shan uplift at 5 Ma . The prior is therefore not directly related to Himalayan events hypothesized in the Tibet or Hunan hypotheses, which mostly concern uplift in Southwest China.
Considering the more recent priors, the Malaya prior is based on Schistosoma and so assumes correspondence between the divergence of parasites and snails. It is not necessary to assume co-evolution, but merely that both the snails and parasites responded in the same way to the same events (i.e., Pleistocene changes in river courses where the snails live and eustatics). Also, although host-switching is not uncommon in Schistosoma, this tends to occur at higher taxonomic levels (of the parasites). As Robertsiella is the only pomatiopsid lineage found in peninsular Malaysia, and the entire S. sinensium-group is restricted to Pomatiopsidae, it is not unreasonable to assume that the ancestor of S. malayensis arrived in Malaysia through dispersal in the ancestor of Robertsiella. The prior Khorat is based on volcanism on the eastern margin of the Khorat Basin and flow reversal of the Mul river (c.a. 0.8 Ma), which are likely to have had an impact on the divergence of Neotricula aperta. The precise dates of the more recent geological events are likely to be revised and refined in future geological studies and, as with the Malaya prior, appropriately wide confidence intervals were applied to this prior.
The calibration priors are summarised in Table 2 and were applied as normally distributed calibration dates to guide the inference of other divergence dates by Beast. The SD of each prior was chosen so that the 95% CI spanned the range of dates (obtained from the literature) given in column two of Table 2. Beast was then set to estimate the divergence times or TMRCAs of key clades in order to help test phylogeographical hypotheses.
Figure 6 shows the phylogenetic tree estimated by Beast (from recent dating runs) and the locations of each date prior and TMRCA estimated. The TMRCAs for the more recent clades were estimated separately from those of deeper divergences because Beast is currently unable to allow the clock rates to vary with time along a lineage. It is likely that clock rates will vary over long time scales, and attempts to estimate TMRCAs for divergences believed to have occurred 400–20 Ma simultaneously with those for divergences on the 15–0.5 Ma time scale, are unlikely to estimate the correct dates (at either scale) [147, 164]. In general the molecular clock rate tends to be greater during evolutionary bursts or radiations (if these involve frequent population subdivision ) and along more recent lineages, than those in the deeper past, because in more recent lineages mutations may be present which, given sufficient time, will be removed by drift [166, 167]. At the same time the effects of saturation accumulate in older divergences. Consequently, those TMRCAs expected to be Miocene or later were estimated separately from expectedly older divergences, and with the inclusion of more recent calibration date priors. Separating the estimates also leads to greater computational speed and replication when performing analyses on multiple machines/cpus.
A Jeffrey’s prior was placed on the clock rate for the indels (following the indications of test runs and BF comparisons). For the remaining partitions an informative prior was applied based on rates published in the literature for these loci. For the cox1 partitions a 1.75% clock rate was applied, with SD set to allow the rate to vary from 0 to 3.5% (at 95% confidence); the rate was allowed to vary to 35% for the third codon positions. For the rrnL data the rate prior was 1.097% (again with 95% CI covering 0-3%). A percent clock rate is equivalent to substitutions per site per year × 10-8. The rate priors were chosen with reference to the following sources. For rrnL, triculine snail divergences over the early Pleistocene , early Pliocene Hydrobiidae , and other divergences of Hydrobiidae . The choice of rate for the cox1 data was guided by rates for Mio-Pliocene freshwater pulmonate snail divergence in Southeast Asia . All runs were also performed with sampling from the prior only (no data) to ensure that the topology, clock rates, range of parameter estimates and tmrcas were not being directed by the priors over the data. Particular care was taken where a TMRCA was estimated for a divergence upon which a prior had also been specified (see initial paragraph of this sub-section).
Availability of supporting data
DNA sequence data generated during this study are deposited in GenBank under accession numbers KC832687-KC832728 (http://www.ncbi.nih.gov/entrez/query.fcgi?db=Nucleotide&cmd=search&term=KC832687). All other data supporting the results of this article are included within the article and its additional files.
Highest posterior density
Long branch attraction
Likelihood ratio test
Megaannum or million years
Markov chain Monte Carlo
ribosomal DNA (ribosomal RNA gene cluster or operon)
Scanning electron Microscopy
Simple indel coding
Time to the most recent common ancestor.
Sampling for this work was supported by Wellcome Trust project grant No. 068706 (2003–2006) and China NSF grant No. 30972565 (2009–2012) to SWA. Thanks are due to the Natural History Museum (London) for providing access to library facilities and to the Mollusca collection. The funding bodies had no role in design, in the collection, analysis, and interpretation of data; nor in the writing of the manuscript; and in the decision to submit the manuscript for publication.
- Davis GM: Evolution of prosobranch snails transmitting Asian Schistosoma; coevolution with Schistosoma: a review. Progr Clin Parasit. 1992, 3: 145-204.Google Scholar
- Attwood SW: Schistosomiasis in the Mekong region: epidemiology and phylogeography. Adv Parasit. 2001, 50: 87-152.Google Scholar
- Attwood SW, Upatham ES, Zhang Y-P, Yang Z-Q, Southgate VR: A DNA-sequence based phylogeny for triculine snails (Gastropoda: Pomatiopsidae: Triculinae), intermediate hosts for Schistosoma (Trematoda: Digenea): phylogeography and the origin of Neotricula. J Zool. 2004, 262: 47-56. 10.1017/S0952836903004424.Google Scholar
- Blair D, Davis GM, Wu B: Evolutionary relationships between trematodes and snails emphasizing schistosomes and paragonimids. Parasitology. 2001, 123: S229-S243.PubMedGoogle Scholar
- Spolsky C, Davis GM, Yi Z: Sequencing methodology and phylogenetic analysis: cytochrome b gene sequence reveals significant diversity in Chinese populations of Oncomelania (Gastropda: Pomatiopsidae). Malacologia. 1996, 38: 213-221.Google Scholar
- Becker RA, Wilks AR, Minka TP, Brownrigg R: Maps: Draw Geographical Maps. 2013, Wellington, New Zealand: Victoria University of WellingtonGoogle Scholar
- Becker RA, Wilks AR, Minka TP, Brownrigg R: Mapdata: Extra Map Databases. R Package. 2013, Wellington, New Zealand: Victoria University of WellingtonGoogle Scholar
- Davis GM, Chen C-E, Wu C, Kuang T-F, Xing X-G, Li L, Liu W-J, Yan Y-L: The Pomatiopsidae of Hunan, China (Gastropoda: Rissoacea). Malacologia. 1992, 34: 143-342.Google Scholar
- Attwood SW: Mekong Schistosomiasis: where Did It Come from and Where Is It Going?. The Mekong: Biophysical Environment of an International River Basin. Edited by: Campbell IC. 2009, New York: Academic Press, 464-Google Scholar
- Davis GM: The origin and evolution of the gastropod family Pomatiopsidae, with emphasis on the Mekong river Triculinae. Acad Nat Sci Philad, Monogr. 1979, 20: 1-120.Google Scholar
- Whelan NV, Johnson PD, Harris PM: Rediscovery of Leptoxis compacta (Anthony, 1854) (Gastropoda: Cerithioidea: Pleuroceridae). PLoS ONE. 2012, 7: e42499-10.1371/journal.pone.0042499.PubMedPubMed CentralGoogle Scholar
- Brandt RA: Descriptions of new non-marine mollusks from Asia. Archiv für Molluskenkunde. 1968, 98: 213-289.Google Scholar
- Davis GM, Kitikoon V, Temcharoen P: Monograph on “Lithoglyphopsis” aperta, the snail host of Mekong river schistosomiasis. Malacologia. 1976, 15: 241-287.PubMedGoogle Scholar
- Davis GM, Kang Z-B: The genus Wuconchona of China (Gastropoda: Pomatiopsidae: Triculinae): anatomy, systematics, cladistics, and transmission of Schistosoma. Proc of the Acad Nat Sci Philad. 1990, 142: 119-142.Google Scholar
- Davis GM, Kuo Y-H, Hoagland KE, Chen P-L, Yang H-M, Chen D-J: Erhaia, a new genus and new species of Pomatiopsidae from China (Gastropoda: Rissoacea). Proc of the Acad Nat Sci Philad. 1985, 137: 48-78.Google Scholar
- Wilke T, Davis GM, Chen C-E, Zhou XN, Zeng XP, Zhang Y, Spolsky C: Oncomelania hupensis (Gastropoda: Rissooidea) in eastern China: molecular phylogeny, population structure, and ecology. Acta Trop. 2000, 77: 215-227. 10.1016/S0001-706X(00)00143-1.PubMedGoogle Scholar
- Wilke T, Davis GM, Falniowski A, Giusti F, Bodon M, Szarowska M: Molecular systematics of Hydrobiidae (Mollusca: Gastropoda: Rissooidea): testing monophyly and phylogenetic relationships. Proc of the Acad Nat Sci Philad. 2001, 151: 1-21. 10.1635/0097-3157(2001)151[0001:MSOHMG]2.0.CO;2.Google Scholar
- Wilke T, Davis GM, Gong X, Liu H-X: Erhaia (Gastropoda: Rissooidea): Phylogenetic relationships and the question of Paragonimus coevolution in Asia. Amer J Trop Med & Hyg. 2000, 62: 453-459.Google Scholar
- Attwood SW, Johnston DA: Nucleotide sequence differences reveal genetic variation in Neotricula aperta (Gastropoda: Pomatiopsidae), the snail host of schistosomiasis in the lower Mekong Basin. Biol J Linn Soc. 2001, 73: 23-41. 10.1111/j.1095-8312.2001.tb01344.x.Google Scholar
- Attwood SW, Ambu S, Meng XH, Upatham ES, Xu F-S, Southgate VR: The phylogenetics of triculine snails (Rissooidea: Pomatiopsidae) from south-east Asia and southern China: historical biogeography and the transmission of human schistosomiasis. J Mollus Stud. 2003, 69: 263-271. 10.1093/mollus/69.3.263.Google Scholar
- Attwood SW, Brown DS, Meng XH, Southgate VR: A new species of Tricula (Pomatiopsidae: Triculinae) from Sichuan Province, PR China: intermediate host of Schistosoma sinensium. Syst & Biodiv. 2003, 1: 109-116. 10.1017/S147720000300104X.Google Scholar
- Perez KE, Ponder WF, Colgan DJ, Clark SA, Lydeard C: Molecular phylogeny and biogeography of spring-associated hydrobiid snails of the Great Artesian Basin, Australia. Mol Phylogenet & Evol. 2005, 34: 545-556. 10.1016/j.ympev.2004.11.020.Google Scholar
- Davis GM, Wilke T, Spolsky C, Zhang Y, Xia M-Y, Rosenberg G: Cytochrome Oxidase I-based phylogenetic relationships among the Hydrobiidae, Pomatiopsidae, Rissoidae, and Truncatellidae (Gastropoda: Prosobranchia: Rissoacea). Malacologia. 1998, 40: 251-266.Google Scholar
- Guan F, Niu AO, Attwood SW, Zhu YH YLLBZ: Molecular phylogenetics of Triculine snails (Gastropoda: Pomatiopsidae) from southern China. Mol Phylogenet Evol. 2008, 48: 702-707. 10.1016/j.ympev.2008.04.021.PubMedGoogle Scholar
- Liu Y-Y, Zhang W-Z, Wang Y-X: Studies on Tricula (Prosobranchia: hydrobiidae) From China [in Chinese with English abstract]. Act Zootax Sin. 1983, 8: 135-140.Google Scholar
- Tang C, Liu GY, Xiong TS, Lian ZY, Li MB: A newly recorded host of Pagumogonimus skrjabini: Tricula hongshanensis sp. nov [in Chinese with English abstract]. Act Hydrob Sin. 1986, 10: 286-289.Google Scholar
- Hu W, Zhou SH, Long ZP: Study on the laboratory raising of a tricula snail [in Chinese with English abstract]. Chin J Parasitol & Parasit Dis. 1994, 3: 216-218.Google Scholar
- Kameda Y, Kato M: Terrestrial invasion of pomatiopsid gastropods in the heavy-snow region of the Japanese Archipelago. BMC Evol Biol. 2011, 11: 118-10.1186/1471-2148-11-118.PubMedPubMed CentralGoogle Scholar
- Hai H-T, Thanh DN: Contribution to the Study on Triculinae snail group (Pomatiopsidae - Mollusca) in Tay Nguyen highland, Vietnam [in Vietnamese]. Tạp Chí Sinh Học (J Biol Viet). 2010, 32: 14-18.Google Scholar
- Hai HT, Thanh DN: Two new species of snail Triculinae group (Pomatiopsidae - Prosobranchia) from Nothern Vietnam [in Vietnamese]. Tạp Chí Sinh Học (J Biol Viet). 2012, 33: 17-23.Google Scholar
- Hutchinson CS: Geological Evolution of South-East Asia. 1989, Oxford: Clarendon PressGoogle Scholar
- Yamada K, Tanaka Y, Irizuki T: Paleoceanographic shifts and global events recorded in late Pliocene shallow marine deposits (2.80–2.55 Ma) of the Sea of Japan. Palaeogeogr, Palaeoclim, Palaeoecol. 2005, 220: 255-271. 10.1016/j.palaeo.2005.01.006.Google Scholar
- Dey AK: Geology of India. 1968, National Book Trust: New DelhiGoogle Scholar
- Hall R: The plate tectonics of cenozoic Asia and the distribution of land and sea. Biogeography and Geological Evolution of SE Asia. Edited by: Hall R, Holloway JD. 1998, Leiden: Backbuys, 99-131.Google Scholar
- Hall R: Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer-based reconstructions, model and animations. J Asian Earth Sci. 2002, 20: 353-431. 10.1016/S1367-9120(01)00069-4.Google Scholar
- Hartt F: Geology and Physical Geography of Brazil. 1870, Cambridge: Boston, Fields, Osgood & Co.Google Scholar
- Ministry of Geology A: Carta Geologico-Economica de la Republica Argentina. 1947, Buenos Aires: Direccion General de Minas y GeologiaGoogle Scholar
- Van Royen W: Mineral Resources of the World. 1952, New York: Prentice-Hall, 2-Google Scholar
- Yamada N, Teraoka Y, Hata M: Geological Atlas of Japan. 1987, Tokyo: Geological Survey of JapanGoogle Scholar
- Abbott RT: Handbook of medically important mollusks of the Orient and the western Pacific. Bull Mus Comp Zoo. 1948, 100: 269-328.Google Scholar
- Nesemann H, Sharma S, Sharma G, Khanal SN, Pradhan B, Shah DN, Tachamo RD: Aquatic Invertebrates of the Ganga River System: Mollusca, Annelida, Crustacea. 2007, Nesemann: Kathmandu HGoogle Scholar
- Davis GM: Snail Hosts of Asian Schistosoma Infecting man: evolution and Coevolution. In The Mekong Schistosome, Malacological Review, Suppl. 2. Edited by Bruce JI, Sornmani S, Asch HL, Crawford KA. Ann Arbor, Michigan: The University of Michigan; 1980:195–238.Google Scholar
- Davis GM: Different modes of evolution and adaptive radiation in the Pomatiopsidae (Prosobranchia: Mesogastropoda). Malacologia. 1981, 21: 209-262.Google Scholar
- Hausdorf B, Röpstorf P, Riedel F: Relationships and origin of endemic Lake Baikal gastropods (Caenogastropoda: Rissooidea) based on mitochondrial DNA sequences. Mol Phylogen Evol. 2003, 26: 435-443. 10.1016/S1055-7903(02)00365-2.Google Scholar
- Coltice N, Seton M, Rolf T, Müller RD, Tackley PJ: Convergence of tectonic reconstructions and mantle convection models for significant fluctuations in seafloor spreading. Earth & Planet Sci Let. 2013, 383: 92-100.Google Scholar
- McLoughlin S: The breakup history of Gondwana and its impact on pre-Cenozoic floristic provincialism. Aust J Bot. 2001, 49: 271-300. 10.1071/BT00023.Google Scholar
- Seton M, Müller RD, Zahirovic S, Gaina C, Torsvik T, Shephard G, Talsma A, Gurnis M, Turner M, Maus S, Chandler M: Global continental and ocean basin reconstructions since 200 Ma. Earth Sci Rev. 2012, 113: 212-270. 10.1016/j.earscirev.2012.03.002.Google Scholar
- Davis GM, Kuo Y-H, Hoagland KE, Chen P-L, Yang H-M, Chen D-J: Advances in the systematics of the Triculinae (Gastropoda: Pomatiopsidae): the genus Fenouilia of Yunnan, China. Proc Acad Nat Sci Philad. 1983, 135: 177-199.Google Scholar
- Amano K, Taira A: Two-phase uplift of Higher Himalayas since 17 Ma. Geology. 1992, 20: 391-394. 10.1130/0091-7613(1992)020<0391:TPUOHH>2.3.CO;2.Google Scholar
- Brookfield ME: The evolution of the great river systems of southern Asia during the Cenozoic India-Asia collision: rivers draining southwards. Geomorphology. 1998, 22: 285-312. 10.1016/S0169-555X(97)00082-2.Google Scholar
- Veevers JJ: Gondwanaland from 650–500 Ma assembly through 320 Ma merger in Pangea to 185–100 Ma breakup: supercontinental tectonics via stratigraphy and radiometric dating. Earth Sci Rev. 2004, 68: 1-132. 10.1016/j.earscirev.2004.05.002.Google Scholar
- Liu L-P, Eronen JT, Fortelius M: Significant mid-latitude aridity in the middle Miocene of East Asia. Palaeogeogr, Palaeoclim, Palaeoecol. 2009, 279: 201-206. 10.1016/j.palaeo.2009.05.014.Google Scholar
- Attwood SW: Studies on the parasitology, phylogeography and the evolution of host-parasite interactions for the snail intermediate hosts of medically important trematode genera in Southeast Asia. Adv Parasitol. 2010, 73: 405-440.PubMedGoogle Scholar
- Chinzei K: Late Cenozoic zoogeography of the Sea of Japan. Episodes. 1991, 14: 231-235.Google Scholar
- Allen CR, Gillespie AR, Yuan H, Sieh KE, Buchun Z, Chengnan Z: Red river and associated faults Yunnan Province, China: Quaternary geology, slip rates and seismic hazard. Bull Geol Soc Amer. 1984, 95: 686-700. 10.1130/0016-7606(1984)95<686:RRAAFY>2.0.CO;2.Google Scholar
- Chen F-B: Hengduan event: An important tectonic event of the late Cenozoic in eastern Asia [in Chinese with English abstract]. Mount Res. 1992, 10: 195-202.Google Scholar
- Métivier F, Gaudemer Y: Stability of output fluxes of large rivers in South and East Asia during the last 2 million years: implications on floodplain processes. Basin Res. 1999, 11: 293-303. 10.1046/j.1365-2117.1999.00101.x.Google Scholar
- Xiao J, An Z: Three large shifts in East Asian monsoon circulation indicated by loess-paleosol sequences in China and late Cenozoic deposits in Japan. Palaeogeogr, Palaeoclim, Palaeoecol. 1999, 154: 179-189. 10.1016/S0031-0182(99)00110-8.Google Scholar
- Attwood SW, Lokman HS, Ong KY: Robertsiella silvicola, a new species of triculine snail (Caenogastropoda: Pomatiopsidae) from peninsular Malaysia, intermediate host of Schistosoma malayensis (Trematoda: Digenea). J Mollus Stud. 2005, 71: 379-391. 10.1093/mollus/eyi040.Google Scholar
- Attwood SW, Fatih FA, Upatham ES: A study of DNA-sequence variation among Schistosoma mekongi (Trematoda: Digenea) populations and related taxa; phylogeography and the current distribution of Asian schistosomiasis. PLoS Negl Trop Dis. 2008, 2: 1-10.Google Scholar
- Workman DR: Geology of Laos, Cambodia, south Vietnam and the eastern part of Thailand: a review. 1972, London: Institute of Geological Science Overseas Division Report, 19-Google Scholar
- Newton RB: On some freshwater fossils from central South Africa. Annal & Mag Nat Hist. 1920, 9: 241-249.Google Scholar
- Neubauer TA, Schneider S, Böhme M, Prieto J: First records of freshwater rissooidean gastropods from the Palaeogene of Southeast Asia. J Mollus Stud. 2012, 78: 275-282. 10.1093/mollus/eys011.Google Scholar
- Attwood SW, Upatham ES, Meng XH, Qiu D-C, Southgate VR: The phylogeography of Asian Schistosoma (Trematoda: Schistosomatidae). Parasitology. 2002, 125: 1-13.Google Scholar
- Lacassin R, Leloup PH, Tapponier P: Bounds on strain in large Tertiary shear zones of SE Asia from boudinage restoration. J Struct Geol. 1993, 15: 677-692. 10.1016/0191-8141(93)90055-F.Google Scholar
- Workman DR: Geology of Laos, Cambodia, South Vietnam and the Eastern Part of Thailand. 1977, London: Institute of Geological SciencesGoogle Scholar
- Okamoto M, Lo C-T, Tiu WU, Qiu D, Hadidjaja P, Upatham S, Sugiyama H, Taguchi T, Hirai H, Saitoh Y, Habe S, Kawanaka M, Hirata M, Agatsuma T: Phylogenetic relationships of snails of the genera Oncomelania and Tricula inferred from the mitochondrial 12S rRNA gene. Jap J Trop Med & Hyg. 2003, 31: 5-10. 10.2149/tmh1973.31.5.Google Scholar
- Woodruff DS, Carpenter MP, Upatham ES, Viyanant V: Molecular phylogeography of Oncomelania lindoensis (Gastropoda: Pomatiopsidae), the intermediate host of Schistosoma japonicum in Sulawesi. J Mollus Stud. 1999, 65: 21-32. 10.1093/mollus/65.1.21.Google Scholar
- Barron EJ: Cretaceous plate tectonic reconstructions. Palaeogeogr, Palaeoclim, Palaeoecol. 1987, 59: 3-29.Google Scholar
- Bagulayan A, Bartlett-Roa J, Carter A, Inman B, Keen E, Orenstein E, Patin N, Sato K, Sibert E, Simonis A, Van Cise A, Franks P: Journey to the Center of the Gyre: the fate of the Tohoku tsunami debris field. Oceanography. 2012, 25: 200-207. 10.5670/oceanog.2012.55.Google Scholar
- Kodandaramaiah U: Tectonic calibrations in molecular dating. Curr Zool. 2011, 57: 116-124.Google Scholar
- Waters JM, Rowe DL, Apte S, King TM, Wallis GP, Anderson L, Norris RJ, Craw D, Burridge CP: Geological dates and molecular rates: rapid divergence of rivers and their biotas. Syst Biol. 2007, 56: 271-282. 10.1080/10635150701313855.PubMedGoogle Scholar
- Heads M: Dating nodes on molecular phylogenies: a critique of molecular biogeography. Cladistics. 2005, 21: 62-78. 10.1111/j.1096-0031.2005.00052.x.Google Scholar
- Donoghue MJ, Moore BR: Toward an integrative historical biogeography. Integr Comp Biol. 2003, 43: 261-270. 10.1093/icb/43.2.261.PubMedGoogle Scholar
- Upchurch P: Gondwanan break-up: legacies of a lost world?. Trends Ecol Evol (Amst). 2008, 23: 229-236. 10.1016/j.tree.2007.11.006.Google Scholar
- Strange RM, Burr BM: Intraspecific phylogeography of North American highland fishes: a test of the Pleistocene vicariance hypothesis. Evolution. 1997, 51: 885-897. 10.2307/2411163.Google Scholar
- Burridge CP: Biogeographic history of geminate cirrhitoids (Perciformes: Cirrhitoidea) with east–west allopatric distributions across southern Australia, based on molecular data. Glob Ecol Biogeogr. 2000, 9: 517-525. 10.1046/j.1365-2699.2000.00204.x.Google Scholar
- Perdices A, Doadrio I: The molecular systematics and biogeography of the European cobitids based on mitochondrial DNA sequences. Mol Phylogenet Evol. 2001, 19: 468-478. 10.1006/mpev.2000.0900.PubMedGoogle Scholar
- Ovenden JR, White RWG, Adams M: Mitochondrial and allozyme genetics of two Tasmanian galaxiids (Galaxias auratus and G. tanycephalus, Pisces: Galaxiidae) with restricted lacustrine distributions. Heredity. 1993, 70: 223-230. 10.1038/hdy.1993.33.Google Scholar
- Waters JM, López JA, Wallis GP: Molecular phylogenetics and biogeography of galaxiid fishes (Osteichthyes: Galaxiidae): dispersal, vicariance, and the position of Lepidogalaxias salamandroides. Syst Biol. 2000, 49: 777-795. 10.1080/106351500750049824.PubMedGoogle Scholar
- Sanderson MJ: A nonparametric approach to estimating divergence times in the absence of rate constancy. Mol Biol & Evol. 1997, 14: 1218-1231. 10.1093/oxfordjournals.molbev.a025731.Google Scholar
- Mulcahy DG, Noonan BP, Moss T, Townsend TM, Reeder TW, Sites JW, Wiens JJ: Estimating divergence dates and evaluating dating methods using phylogenomic and mitochondrial data in squamate reptiles. Mol Phylogenet Evol. 2012, 65: 974-991. 10.1016/j.ympev.2012.08.018.PubMedGoogle Scholar
- Arbogast BS, Edwards SV, Wakeley J, Beerli P, Slowinski JB: Estimating divergence times from molecular data on phylogenetic and population genetic timescales. Ann Rev Ecol & System. 2002, 33: 707-740. 10.1146/annurev.ecolsys.33.010802.150500.Google Scholar
- Xia X, Xie Z, Salemi M, Chen L, Wang Y: An index of substitution saturation and its application. Mol Phylogenet Evol. 2003, 26: 1-7. 10.1016/S1055-7903(02)00326-3.PubMedGoogle Scholar
- Ponder WF: Prosobranch phylogeny: the truncatelloidean (= Rissoacean) Radiation - a Preliminary Phylogeny. Malacological Review, Suppl. Edited by: Ponder WF, Ann A. 1988, Michigan: The University of Michigan, 129-164.Google Scholar
- Bandel K: Caenogastropoda during Mesozoic times. Script Geolog, Spec. 1993, 2: 7-56.Google Scholar
- Bandel K, Riedel F: The Late Cretaceous gastropod fauna from Ajka (Bakony Mountains, Hungary): a revision. Ann Naturhist Mus Wien. 1994, 96A: 1-65.Google Scholar
- Métivier F, Gaudemer Y, Tapponnier P, Klein M: Mass accumulation rates in Asia during the Cenozoic. Geophys J Inter. 1999, 137: 280-318.Google Scholar
- Marko PB: Fossil calibration of molecular clocks and the divergence times of geminate species pairs separated by the isthmus of panama. Mol Biol Evol. 2002, 19: 2005-2021. 10.1093/oxfordjournals.molbev.a004024.PubMedGoogle Scholar
- Schultheiß R, Bocxlaer BV, Wilke T, Albrecht C: Old fossils–young species: evolutionary history of an endemic gastropod assemblage in Lake Malawi. Proc R Soc B. 2009, 276: 2837-2846. 10.1098/rspb.2009.0467.PubMedPubMed CentralGoogle Scholar
- Attwood SW, Fatih FA, Campbell IC, Upatham ES: The distribution of Mekong schistosomiasis, past and future: preliminary indications from an analysis of genetic variation in the intermediate host. Parasitol Int. 2008, 57: 256-270. 10.1016/j.parint.2008.04.003.PubMedGoogle Scholar
- Cunningham CW, Blackstone NW, Buss LW: Evolution of king crabs from hermit crab ancestors. Nature. 1992, 355: 539-542. 10.1038/355539a0.PubMedGoogle Scholar
- Xia L, Li X, Ma Z, Xu X, Xia Z: Cenozoic volcanism and tectonic evolution of the Tibetan plateau. Gondwana Res. 2011, 19: 850-866. 10.1016/j.gr.2010.09.005.Google Scholar
- Stephenson D, Marshall TR: The petrology and mineralogy of Mt. Popa Volcano and the nature of the late-Cenozoic Burma Volcanic Arc. J Geol Soc. 1984, 141: 747-762. 10.1144/gsjgs.141.4.0747.Google Scholar
- Kano K, Yoshikawa T, Yanagisawa Y, Ogasawara K, Danhara T: An unconformity in the early Miocene syn-rifting succession, northern Noto Peninsula, Japan: Evidence for short-term uplifting precedent to the rapid opening of the Japan Sea. Island Arc. 2002, 11: 170-184. 10.1046/j.1440-1738.2002.00363.x.Google Scholar
- Fortelius M, Eronen J, Jernvall J, Liu L, Pushkina D, Rinne J, Tesakov A, Vislobokova I: Fossil mammals resolve regional patterns of Eurasian climate change over 20 million years. Europe. 2002, 4: 1005-1016.Google Scholar
- Itoh Y, Nagasaki Y: Crustal shortening of Southwest Japan in the Late Miocene. Island Arc. 1996, 5: 337-353. 10.1111/j.1440-1738.1996.tb00035.x.Google Scholar
- Harrison TM, Copeland P, Kidd WSF, Yin A: Raising Tibet. Science. 1992, 255: 1663-1670. 10.1126/science.255.5052.1663.PubMedGoogle Scholar
- Zhisheng A, Kutzbach JE, Prell WL, Porter SC: Evolution of Asian monsoons and phased uplift of the Himalaya-Tibetan plateau since Late Miocene times. Nature. 2001, 411: 62-66. 10.1038/35075035.PubMedGoogle Scholar
- Zhu H-H: Formation of fault lakes in Yunnan province, and its sediment and evolution in late Cenozoic. Ocean Limnol Sin. 1991, 22: 509-516.Google Scholar
- Brigham-Grette J, Melles M, Minyuk P, Andreev A, Tarasov P, DeConto R, Koenig S, Nowaczyk N, Wennrich V, Rosén P, Haltia E, Cook T, Gebhardt C, Meyer-Jacob C, Snyder J, Herzschuh U: Pliocene warmth, polar amplification, and stepped Pleistocene cooling recorded in NE arctic Russia. Science. 2013, 340: 1421-1427. 10.1126/science.1233137.PubMedGoogle Scholar
- Zhang YG, Ji J, Balsam W, Liu L, Chen J: Mid-Pliocene Asian monsoon intensification and the onset of Northern Hemisphere glaciation. Geology. 2009, 37: 599-602. 10.1130/G25670A.1.Google Scholar
- Attwood SW: A demographic analysis of y-Neotricula aperta (Gastropoda: Pomatiopsidae) populations in Thailand and southern Laos, in relation to the transmission of schistosomiasis. J Mollus Stud. 1995, 61: 29-42. 10.1093/mollus/61.1.29.Google Scholar
- Davis GM, Subba Rao NV, Hoagland KE: In search of Tricula (Gastropoda: Pomatiopsidae): Tricula defined, and a new genus described. Proc Acad Nat Sci Philad. 1986, 138: 426-442.Google Scholar
- Brandt RA: New freshwater gastropods from the Mekong. Archiv für Molluskenkunde. 1970, 100: 183-205.Google Scholar
- Brandt RA, Temcharoen P: The molluscan fauna of the Mekong at the foci of schistosomiasis in South Laos and Cambodia. Archiv für Molluskenkunde. 1971, 101: 111-140.Google Scholar
- Kitikoon V: Studies on Tricula aperta and related taxa, the snail intermediate hosts of Schistosoma mekongi. IIII. Susceptibility studies. Malacol Rev. 1981, 14: 37-42.Google Scholar
- Lydeard C, Holznagel WE, Glaubrecht M, Ponder WF: Molecular phylogeny of a circum-global, diverse gastropod superfamily (Cerithioidea: mollusca: Caenogastropoda): pushing the deepest phylogenetic limits of mitochondrial LSU rDNA sequences. Mol Phylogenet Evol. 2002, 22: 399-406. 10.1006/mpev.2001.1072.PubMedGoogle Scholar
- Wilke T, Davis GM: Infraspecific mitochondrial sequence diversity in Hydrobia ulvae and Hydrobia ventrosa (Hydrobiidae: Rissooidea: Gastropoda): do their different life histories affect biogeographic patterns and gene flow?. Biol J Linn Soc. 2000, 70: 89-105. 10.1111/j.1095-8312.2000.tb00202.x.Google Scholar
- Brandt RA: The non-marine mollusca of Thailand. Archiv für Molluskenkunde. 1974, 105: 1-423.Google Scholar
- Davis GM, Liu YY, Chen YG: New genus of Triculinae (Prosobranchia: Pomatiopsidae) from China: phylogenetic relationships. Proc Acad Nat Sci Philad. 1990, 142: 143-165.Google Scholar
- Davis GM, Guo YH, Hoagland KE, Chen P-L, Zheng LC, Yang H-M, Chen D-J, Zhou YF: Anatomy and systematics of Triculini (Prosobranchia: Pomatiopsidae: Triculinae), freshwater snails from Yunnan, China, with descriptions of new species. Proc Acad Nat Sci Philad. 1986, 138: 466-575.Google Scholar
- Davis GM: New Tricula from Thailand. Archiv für Molluskenkunde. 1968, 98: 291-317.Google Scholar
- Davis GM, Kang Z-B: Advances in the systematics of Erhaia (Gastropoda: Pomatiopsidae) from the People’s Republic of China. Proc Acad Nat Sci Philad. 1995, 146: 391-427.Google Scholar
- Davis GM, Subba Rao NV: Discovery of Erhaia (Gastropoda: Pomatiopsidae) in northern India with description of new genus of Erhaiini from China. Proc Acad Nat Sci Philad. 1997, 148: 273-299.Google Scholar
- Winnepenninck B, Backeljau T, De Wachter R: Extraction of high molecular weight DNA from molluscs. Trend Genet. 1993, 9: 407-Google Scholar
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol & Biotech. 1994, 3: 294-299.Google Scholar
- Palumbi SR, Martin A, Romano S, McMillan WO, Stice L, Grabowski G: The Simple fool’s Guide to PCR. Version 2.0. 1991, Honolulu: University of HawaiiGoogle Scholar
- Hebert PDN, Cywinska A, Ball SL, deWaard JR: Biological identifications through DNA barcodes. Proc Biol Sci. 2003, 270: 313-321. 10.1098/rspb.2002.2218.PubMedPubMed CentralGoogle Scholar
- Remigio EA, Blair D: Molecular systematics of the freshwater snail family Lymnaeidae (Pulmonata: Basommatophora) utilising mitochondrial ribosomal DNA sequences. J Mollus Stud. 1997, 63: 173-185. 10.1093/mollus/63.2.173.Google Scholar
- Ashelford KE, Chuzhanova NA, Fry JC, Jones AJ, Weightman AJ: At least 1 in 20 16S rRNA sequence records currently held in public repositories is estimated to contain substantial anomalies. Appl Environ Microbiol. 2005, 71: 7724-7736. 10.1128/AEM.71.12.7724-7736.2005.PubMedPubMed CentralGoogle Scholar
- Clackson T, Güssow D, Jones PT: General Applications of PCR to Gene Cloning and Manipulation. PCR: A Practical Approach, Volume 1. Edited by: McPherson MJ, Quirke P, Taylor GR. 1991, Oxford: OUP, 187-214.Google Scholar
- Bergsten J: A review of long-branch attraction. Cladistics. 2005, 21: 163-193. 10.1111/j.1096-0031.2005.00059.x.Google Scholar
- O’Connor T, Sundberg K, Carroll H, Clement M, Snell Q: Analysis of long branch extraction and long branch shortening. BMC Genomics. 2010, 11: S14-PubMedPubMed CentralGoogle Scholar
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.PubMedGoogle Scholar
- DeSalle R, Freedman T, Prager EM, Wilson AC: Tempo and mode of sequence evolution in mitochondrial DNA of Hawaiian Drosophila. J Mol Evol. 1987, 26: 157-164. 10.1007/BF02111289.PubMedGoogle Scholar
- Xia X: DAMBE (Data Analysis in Molecular Biology and Evolution) Software, Version 3.7.48. 1999, Hong kong: University of Hong kong (Department of Ecology and Biodiversity)Google Scholar
- Simmons MP, Ochoterena H: Gaps as characters in sequence-based phylogenetic analyses. Syst Biol. 2000, 49: 369-381. 10.1093/sysbio/49.2.369.PubMedGoogle Scholar
- Farris JS, Kallersjo M, Kluge AG, Bult C: Constructing a significance test for incongruence. Syst Biol. 1995, 44: 570-572.Google Scholar
- Swofford D: PAUP: Phylogenetic Analysis Using Parsimony (and Other Methods). Version 4.0 Beta Version. 2002, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
- Lee MS: Uninformative characters and apparent conflict between molecules and morphology. Mol Biol Evol. 2001, 18: 676-680. 10.1093/oxfordjournals.molbev.a003848.PubMedGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D: jModelTest 2: more models, new heuristics and parallel computing. Nat Meth. 2012, 9: 772-772.Google Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.PubMedGoogle Scholar
- Foster PG: Modeling compositional heterogeneity. Syst Biol. 2004, 53: 485-495. 10.1080/10635150490445779.PubMedGoogle Scholar
- Gascuel O, Steel M: Inferring ancestral sequences in taxon-rich phylogenies. Math Biosci. 2010, 227: 125-135. 10.1016/j.mbs.2010.07.002.PubMedGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.PubMedPubMed CentralGoogle Scholar
- Alekseyenko AV, Lee CJ, Suchard MA: Wagner and Dollo: a stochastic duet by composing two parsimonious solos. Syst Biol. 2008, 57: 772-784. 10.1080/10635150802434394.PubMedPubMed CentralGoogle Scholar
- Bouckaert R, Lemey P, Dunn M, Greenhill SJ, Alekseyenko AV, Drummond AJ, Gray RD, Suchard MA, Atkinson QD: Mapping the origins and expansion of the Indo-European language family. Science. 2012, 337: 957-960. 10.1126/science.1219669.PubMedPubMed CentralGoogle Scholar
- Felsenstein J: Cases in which parsimony or compatibility methods will be positively misleading. Syst Zool. 1978, 27: 401-410. 10.2307/2412923.Google Scholar
- Brown JM, Hedtke SM, Lemmon AR, Lemmon EM: When trees grow too long: investigating the causes of highly inaccurate Bayesian branch-length estimates. Syst Biol. 2010, 59: 145-161. 10.1093/sysbio/syp081.PubMedGoogle Scholar
- Marshall DC: Cryptic failure of partitioned Bayesian phylogenetic analyses: lost in the land of long trees. Syst Biol. 2010, 59: 108-117. 10.1093/sysbio/syp080.PubMedGoogle Scholar
- Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.PubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.PubMedGoogle Scholar
- Holder M, Lewis PO: Phylogeny estimation: traditional and Bayesian approaches. Nat Rev Genet. 2003, 4: 275-284. 10.1038/nrg1044.PubMedGoogle Scholar
- Van Dongen S: Prior specification in Bayesian statistics: three cautionary tales. J Theor Biol. 2006, 242: 90-100. 10.1016/j.jtbi.2006.02.002.PubMedGoogle Scholar
- Sanderson MJ: Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Mol Biol Evol. 2002, 19: 101-109. 10.1093/oxfordjournals.molbev.a003974.PubMedGoogle Scholar
- Ho SYW, Phillips MJ, Cooper A, Drummond AJ: Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. Mol Biol Evol. 2005, 22: 1561-1568. 10.1093/molbev/msi145.PubMedGoogle Scholar
- Pulquério MJ, Nichols RA: Dates from the molecular clock: how wrong can we be?. Trends Ecol Evol. 2007, 22: 180-184. 10.1016/j.tree.2006.11.013.PubMedGoogle Scholar
- Andrés Varón LSV: POY version 4: phylogenetic analysis using dynamic homologies. Cladistics. 2010, 26: 72-85. 10.1111/j.1096-0031.2009.00282.x.Google Scholar
- Wheeler W: Fixed character states and the optimization of molecular sequence data. Cladistics. 1999, 15: 379-385. 10.1111/j.1096-0031.1999.tb00274.x.Google Scholar
- Liu L, Mondal M, Idris M, Lokman H, Rajapakse PJ, Satrija F, Diaz J, Upatham ES, Attwood S: The phylogeography of Indoplanorbis exustus (Gastropoda: Planorbidae) in Asia. Parasites Vectors. 2010, 3: 57-10.1186/1756-3305-3-57.PubMedPubMed CentralGoogle Scholar
- Shimodaira H, Hasegawa M: Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol Biol Evol. 1999, 16: 1114-1116. 10.1093/oxfordjournals.molbev.a026201.Google Scholar
- De Jong R, Congdon TCE: The Montane butterflies of the Eastern Afrotropics. Biogeography and Ecology of the Rain Forests of Eastern Africa. Edited by: Lovett JC, Wasser SK. 1993, Cambridge, UK: Cambridge University Press, 133-172.Google Scholar
- Budd GE: The Cambrian fossil record and the origin of the phyla. Integr Comp Biol. 2003, 43: 157-165. 10.1093/icb/43.1.157.PubMedGoogle Scholar
- Heads M: Integrating earth and life sciences in New Zealand natural history: the parallel arcs model. NZJ Zoo. 1989, 16: 549-585. 10.1080/03014223.1989.10422922.Google Scholar
- Craw RC, Grehan JR, Heads MJ: Panbiogeography: Tracking the History of Life: Tracking the History of Life. 1999, Oxford, UK: Oxford University PressGoogle Scholar
- Kabat AR, Hershler R: The prosobranch snail family hydrobiidae (Gastropoda: Rissooidea): review of classification and supraspecific taxa. Smithson Contrib Zool. 1993, 547: 1-94.Google Scholar
- Tracey S, Todd JA, Erwin DH: Mollusca: Gastropoda. The Fossil Record 2. Edited by: Benton MJ. 1993, London: Chapman and Hall, 131-167.Google Scholar
- Wenz W: Gastropoda extramarina tertiaria. Foss Cat I Anim. 1926, 32: 1863-2230.Google Scholar
- Bandel K: Gastropods from brackish and fresh water of the Jurassic-Cretaceous transition (a systematic reevaluation). Berl Geowiss Abh. 1991, 134: 9-55.Google Scholar
- Bell DG: Norwegian North Sea Well 31/2-2. 1980, Trondheim, Norway: Continental Shelf Institute (IKU)Google Scholar
- Dejong RJ, Morgan JA, Wilson WD, Al-Jaser MH, Appleton CC, COULIBALY G, Doenhoff MJ, Haas W, IDRIS MA, Magalhaes LA, MONE H, MOUAHID G, Mubila L, Pointier J-P, Webster JP, Zanotti-Magalhaes EM, Paraense WL, Mkoji GM, Loker ES: Phylogeography of Biomphalaria glabrata and B. pfeifferi, the two most important hosts of Schistosoma mansoni in the New and Old world tropics. Mol Ecol. 2003, 12: 3041-3056. 10.1046/j.1365-294X.2003.01977.x.PubMedGoogle Scholar
- Attwood SW, Fatih FA, Mondal MM, Alim MA, Fadjar S, Rajapakse RPV, Rollinson D: A DNA-sequence based study of the Schistosoma indicum (Trematoda: Digenea) group; population phylogeny, taxonomy and historical biogeography. Parasitology. 2007, 134: 2009-2020.PubMedGoogle Scholar
- Ho SYW, Larson G: Molecular clocks: when times are a-changin’. Trends Genet. 2006, 22: 79-83. 10.1016/j.tig.2005.11.006.PubMedGoogle Scholar
- Bromham L, Cardillo M: Testing the link between the latitudinal gradient in species richness and rates of molecular evolution. J Evol Biol. 2003, 16: 200-207. 10.1046/j.1420-9101.2003.00526.x.PubMedGoogle Scholar
- MacEachern S, McEwan J, Mather A, McCulloch A, Sunnucks P, Goddard M: Testing the neutral theory of molecular evolution using genomic data: a comparison of the human and bovine transcriptome. Genet Sel Evol. 2006, 38: 321-341. 10.1186/1297-9686-38-3-321.PubMedPubMed CentralGoogle Scholar
- Kivisild T, Shen P, Wall DP, Do B, Sung R, Davis K, Passarino G, Underhill PA, Scharfe C, Torroni A, Scozzari R, Modiano D, Coppa A, de Knijff P, Feldman M, Cavalli-Sforza LL, Oefner PJ: The role of selection in the evolution of human mitochondrial genomes. Genetics. 2006, 172: 373-387.PubMedPubMed CentralGoogle Scholar
- Wilke T: Salenthydrobia gen. nov. (Rissooidea: Hydrobiidae): a potential relict of the Messinian salinity crisis. Zool J Linn Soc. 2003, 137: 319-336. 10.1046/j.1096-3642.2003.00049.x.Google Scholar
- Liu H-P, Hershler R: A test of the vicariance hypothesis of western North American freshwater biogeography. J Biogeog. 2007, 34: 534-548. 10.1111/j.1365-2699.2006.01611.x.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.