Phylogeography of the Atlanto-Mediterranean sea cucumber Holothuria (Holothuria) mammata: the combined effects of historical processes and current oceanographical pattern

We assessed the genetic structure of populations of the widely distributed sea cucumber Holothuria (Holothuria) mammata Grube, 1840, and investigated the effects of marine barriers to gene flow and historical processes. Several potential genetic breaks were considered, which would separate the Atlantic and Mediterranean basins, the isolated Macaronesian Islands from the other locations analysed, and the Western Mediterranean and Aegean Sea (Eastern Mediterranean). We analysed mitochondrial 16S and COI gene sequences from 177 individuals from four Atlantic locations and four Mediterranean locations. Haplotype diversity was high (H = 0.9307 for 16S and 0.9203 for COI), and the haplotypes were closely related (p = 0.0058 for 16S and 0.0071 for COI). The lowest genetic diversities were found in the Aegean Sea population. Our results showed that the COI gene was more variable and more useful for the detection of population structure than the 16S gene. The distribution of mtDNA haplotypes, the pairwise FST values and the results of exact tests and AMOVA revealed: (i) a significant genetic break between the population in the Aegean Sea and those in the other locations, as supported by both mitochondrial genes, and (ii) weak differentiation of the Canary and Azores Islands from the other populations; however, the populations from the Macaronesian Islands, Algarve and West Mediterranean could be considered to be a panmictic metapopulation. Isolation by distance was not identified in H. (H.) mammata. Historical events behind the observed findings, together with the current oceanographic patterns, were proposed and discussed as the main factors that determine the population structure and genetic signature of H. (H.) mammata

Comparative phylogeography of alpine invertebrates indicates deep lineage diversification and historical refugia in the Australian Alps

Aim: Comparative phylogeographic analyses of alpine biota from the Northern Hemisphere have linked patterns of genetic diversification to glacial expansion and contraction events in the Pliocene and Pleistocene. Furthermore, the extent of diversification across species groups appears to be associated with vagility. In this study we test whether these patterns apply to a geologically stable system from eastern Australia with comparatively shallow elevational gradients and minimal influence from historical glacial activity. Location: The Australian Alps, Victoria, eastern Australia. Methods: We considered phylogeographic patterns across five alpine invertebrate species based on mitochondrial and nuclear DNA sequence data. Bayesian inference methods were used to estimate species phylogenies and divergence times among lineages. GIS tools were used to map interpopulation genetic divergence and intrapopulation genetic diversity estimates and to visualize spatial patterns across species, providing insights into patterns of endemism and demographic history. Results: Phylogeographic patterns and the timing of lineage diversification were consistent across taxonomic groups. Mountain summits harbour highly differentiated haplogroups, including summits connected by high-elevational plateaus, pointing to diversifications being maintained since the early to mid-Pleistocene. These findings are consistent with previous studies of alpine mammals and reptiles, demonstrating a high degree of endemism in this region, regardless of species vagility. Main conclusions: The fine spatial scales at which deep genetic differentiation among alpine communities was observed in this study are unprecedented. This suggests that glacial periods have had less of an impact on species distributions and genetic diversity than they have in alpine systems in the Northern Hemisphere. Historical gene flow among sky-island populations has been limited despite connecting snowlines during glacial periods, suggesting that factors other than snow cover have influenced patterns of gene flow in this region. These findings emphasize the unique phylogeographic history affecting Victorian alpine biodiversity, and the importance of conserving biodiversity from multiple mountain summits in this region of high endemism.

Phylogeography of the antilopine wallaroo (Macropus antilopinus) across tropical northern Australia

The distribution of antilopine wallaroo, Macropus antilopinus, is marked by a break in the species’ range between Queensland and the Northern Territory, coinciding with the Carpentarian barrier. Previous work on M. antilopinus revealed limited genetic differentiation between the Northern Territory and Queensland M. antilopinus populations across this barrier. The study also identified a number of divergent lineages in the Northern Territory, but was unable to elucidate any geographic structure. Here, we re-examine these results to (1) determine phylogeographic patterns across the range of M. antilopinus and (2) infer the biogeographic barriers associated with these patterns. The tropical savannahs of northern Australia: from the Cape York Peninsula in the east, to the Kimberley in the west. We examined phylogeographic patterns in M. antilopinus using a larger number of samples and three mtDNA genes: NADH dehydrogenase subunit 2, cytochrome b, and the control region. Two datasets were generated and analyzed: (1) a subset of samples with all three mtDNA regions concatenated together and (2) all samples for just control region sequences that included samples from the previous study. Analysis included generating phylogenetic trees based on Bayesian analysis and intraspecific median-joining networks. The contemporary spatial structure of M. antilopinus mtDNA lineages revealed five shallow clades and a sixth, divergent lineage. The genetic differences that we found between Queensland and Northern Territory M. antilopinus samples confirmed the split in the geographic distribution of the species. We also found weak genetic differentiation between Northern Territory samples and those from the Kimberley region of Western Australia, possibly due to the Kimberley Plateau–Arnhem Land barrier. Within the Northern Territory, two clades appear to be parapatric in the west, while another two clades are broadly sympatric across the Northern Territory. MtDNA diversity of M. antilopinus revealed an unexpectedly complex evolutionary history involving multiple sympatric and parapatric mtDNA clades across northern Australia. These phylogeographic patterns highlight the importance of investigating genetic variation across distributions of species and integrating this information into biodiversity conservation.

Concordance in phylogeography and ecological niche modelling identify dispersal corridors for reptiles in arid Australia

Aim: Using the rock-specialist agamid Ctenophorus caudicinctus as a model, we test hypothesized biogeographical dispersal corridors for lizards in the Australian arid zone (across the western sand deserts), and assess how these dispersal routes have shaped phylogeographical structuring. Location: Arid and semi-arid Australia. Methods: We sequenced a c. 1400 bp fragment of mtDNA (ND2) for 134 individuals of C. caudicinctus as well as a subset of each of the mtDNA clades for five nuclear loci (BDNF, BACH1, GAPD, NTF3, and PRLR). We used phylogenetic methods to assess biogeographical patterns within C. caudicinctus, including relaxed molecular clock analyses to estimate divergence times. Ecological niche modelling (Maxent) was employed to estimate the current distribution of suitable climatic envelopes for each lineage. Results: Phylogenetic analyses identified two deeply divergent mtDNA clades within C. caudicinctus - an eastern and western clade - separated by the Western Australian sand deserts. However, divergences pre-date the Pleistocene sand deserts. Phylogenetic analyses of the nuclear DNA data sets generally support major mtDNA clades, suggesting past connections between the western C. c. caudicinctus populations in far eastern Pilbara (EP) and the lineages to the east of the sand deserts. Ecological niche modelling supports the continued suitability of climatic conditions between the Central Ranges and the far EP for C. c. graafi. Main conclusions: Estimates of lineage ages provide evidence of divergence between eastern and western clades during the Miocene with subsequent secondary contact during the Pliocene. Our results suggest that this secondary contact occurred via dispersal between the Central Ranges and the far EP, rather than the more southerly Giles Corridor. These events precede the origins of the western sand deserts and divergence patterns instead appear associated with Miocene and Pliocene climate change.

Deeply divergent mitochondrial lineages reveal patterns of local endemism in chironomids of the Australian Wet Tropics

The Wet Tropics bioregion of north-eastern Australia has been subject to extensive fluctuations in climate throughout the late Pliocene and Pleistocene. Cycles of rainforest contraction and expansion of dry sclerophyll forest associated with such climatic fluctuations are postulated to have played a major role in driving geographical endemism in terrestrial rainforest taxa. Consequences for the distributions of aquatic organisms, however, are poorly understood.The Australian non-biting midge species Echinocladius martini Cranston (Diptera: Chironomidae), although restricted to cool, well-forested freshwater streams, has been considered to be able to disperse among populations located in isolated rainforest pockets during periods of sclerophyllous forest expansion, potentially limiting the effect of climatic fluctuations on patterns of endemism. In this study, mitochondrial COI and 16S data were analysed for E. martini collected from eight sites spanning theWet Tropics bioregion to assess the scale and extent of phylogeographic structure. Analyses of genetic structure showed several highly divergent cryptic lineages with restricted geographical distributions. Within one of the identified lineages, strong genetic structure implied that dispersal among proximate (<1 km apart) streams was extremely restricted. The results suggest that vicariant processes, most likely due to the systemic drying of the Australian continent during the Plio-Pleistocene, might have fragmented historical E. martini populations and, hence, promoted divergence in allopatry.

Influence of historical landscapes, drainage evolution and ecological traits on patterns of genetic diversity in Southeast Asian freshwater snakehead fishes

Snakehead fishes in the family Channidae are obligate freshwater fishes represented by two extant genera, the African Parachannna and the Asian Channa. These species prefer still or slow flowing water bodies, where they are top predators that exercise high levels of parental care, have the ability to breathe air, can tolerate poor water quality, and interestingly, can aestivate or traverse terrestrial habitat in response to seasonal changes in freshwater habitat availability. These attributes suggest that snakehead fishes may possess high dispersal potential, irrespective of the terrestrial barriers that would otherwise constrain the distribution of most freshwater fishes. A number of biogeographical hypotheses have been developed to account for the modern distributions of snakehead fishes across two continents, including ancient vicariance during Gondwanan break-up, or recent colonisation tracking the formation of suitable climatic conditions. Taxonomic uncertainty also surrounds some members of the Channa genus, as geographical distributions for some taxa across southern and Southeast (SE) Asia are very large, and in one case is highly disjunct. The current study adopted a molecular genetics approach to gain an understanding of the evolution of this group of fishes, and in particular how the phylogeography of two Asian species may have been influenced by contemporary versus historical levels of dispersal and vicariance. First, a molecular phylogeny was constructed based on multiple DNA loci and calibrated with fossil evidence to provide a dated chronology of divergence events among extant species, and also within species with widespread geographical distributions. The data provide strong evidence that trans-continental distribution of the Channidae arose as a result of dispersal out of Asia and into Africa in the mid–Eocene. Among Asian Channa, deep divergence among lineages indicates that the Oligocene-Miocene boundary was a time of significant species radiation, potentially associated with historical changes in climate and drainage geomorphology. Mid-Miocene divergence among lineages suggests that a taxonomic revision is warranted for two taxa. Deep intra-specific divergence (~8Mya) was also detected between C. striata lineages that occur sympatrically in the Mekong River Basin. The study then examined the phylogeography and population structure of two major taxa, Channa striata (the chevron snakehead) and the C. micropeltes (the giant snakehead), across SE Asia. Species specific microsatellite loci were developed and used in addition to a mitochondrial DNA marker (Cyt b) to screen neutral genetic variation within and among wild populations. C. striata individuals were sampled across SE Asia (n=988), with the major focus being the Mekong Basin, which is the largest drainage basin in the region. The distributions of two divergent lineages were identified and admixture analysis showed that where they co-occur they are interbreeding, indicating that after long periods of evolution in isolation, divergence has not resulted in reproductive isolation. One lineage is predominantly confined to upland areas of northern Lao PDR to the north of the Khorat Plateau, while the other, which is more closely related to individuals from southern India, has a widespread distribution across mainland SE Asian and Sumatra. The phylogeographical pattern recovered is associated with past river networks, and high diversity and divergence among all populations sampled reveal that contemporary dispersal is very low for this taxon, even where populations occur in contiguous freshwater habitats. C. micropeltes (n=280) were also sampled from across the Mekong River Basin, focusing on the lower basin where it constitutes an important wild fishery resource. In comparison with C. striata, allelic diversity and genetic divergence among populations were extremely low, suggesting very recent colonisation of the greater Mekong region. Populations were significantly structured into at least three discrete populations in the lower Mekong. Results of this study have implications for establishing effective conservation plans for managing both species, that represent economically important wild fishery resources for the region. For C. micropeltes, it is likely that a single fisheries stock in the Tonle Sap Great Lake is being exploited by multiple fisheries operations, and future management initiatives for this species in this region will need to account for this. For C. striata, conservation of natural levels of genetic variation will require management initiatives designed to promote population persistence at very localised spatial scales, as the high level of population structuring uncovered for this species indicates that significant unique diversity is present at this fine spatial scale.

Evolutionary biology of Gondwanan non-biting midges (Diptera: Chironomidae)

The potential restriction to effective dispersal and gene flow caused by habitat fragmentation can apply to multiple levels of evolutionary scale; from the fragmentation of ancient supercontinents driving diversification and speciation on disjunct landmasses, to the isolation of proximate populations as a result of their inability to cross intervening unsuitable habitat. Investigating the role of habitat fragmentation in driving diversity within and among taxa can thus include inferences of phylogenetic relationships among taxa, assessments of intraspecific phylogeographic structure and analyses of gene flow among neighbouring populations. The proposed Gondwanan clade within the chironomid (non-biting midge) subfamily Orthocladiinae (Diptera: Chironomidae) represents a model system for investigating the role that population fragmentation and isolation has played at different evolutionary scales. A pilot study by Krosch et al (2009) indentified several highly divergent lineages restricted to ancient rainforest refugia and limited gene flow among proximate sites within a refuge for one member of this clade, Echinocladius martini Cranston. This study provided a framework for investigating the evolutionary history of this taxon and its relatives more thoroughly. Populations of E. martini were sampled in the Paluma bioregion of northeast Queensland to investigate patterns of fine-scale within- and among-stream dispersal and gene flow within a refuge more rigorously. Data was incorporated from Krosch et al (2009) and additional sites were sampled up- and downstream of the original sites. Analyses of genetic structure revealed strong natal site fidelity and high genetic structure among geographically proximate streams. Little evidence was found for regular headwater exchange among upstream sites, but there was distinct evidence for rare adult flight among sites on separate stream reaches. Overall, however, the distribution of shared haplotypes implied that both larval and adult dispersal was largely limited to the natal stream channel. Patterns of regional phylogeographic structure were examined in two related austral orthoclad taxa – Naonella forsythi Boothroyd from New Zealand and Ferringtonia patagonica Sæther and Andersen from southern South America – to provide a comparison with patterns revealed in their close relative E. martini. Both taxa inhabit tectonically active areas of the southern hemisphere that have also experienced several glaciation events throughout the Plio-Pleistocene that are thought to have affected population structure dramatically in many taxa. Four highly divergent lineages estimated to have diverged since the late Miocene were revealed in each taxon, mirroring patterns in E. martini; however, there was no evidence for local geographical endemism, implying substantial range expansion post-diversification. The differences in pattern evident among the three related taxa were suggested to have been influenced by variation in the responses of closed forest habitat to climatic fluctuations during interglacial periods across the three landmasses. Phylogeographic structure in E. martini was resolved at a continental scale by expanding upon the sampling design of Krosch et al (2009) to encompass populations in southeast Queensland, New South Wales and Victoria. Patterns of phylogeographic structure were consistent with expectations and several previously unrecognised lineages were revealed from central- and southern Australia that were geographically endemic to closed forest refugia. Estimated divergence times were congruent with the timing of Plio-Pleistocene rainforest contractions across the east coast of Australia. This suggested that dispersal and gene flow of E. martini among isolated refugia was highly restricted and that this taxon was susceptible to the impacts of habitat change. Broader phylogenetic relationships among taxa considered to be members of this Gondwanan orthoclad group were resolved in order to test expected patterns of evolutionary affinities across the austral continents. The inferred phylogeny and estimated divergence times did not accord with expected patterns based on the geological sequence of break-up of the Gondwanan supercontinent and implied instead several transoceanic dispersal events post-vicariance. Difficulties in appropriate taxonomic sampling and accurate calibration of molecular phylogenies notwithstanding, the sampling regime implemented in the current study has been the most intensive yet performed for austral members of the Orthocladiinae and unsurprisingly has revealed both novel taxa and phylogenetic relationships within and among described genera. Several novel associations between life stages are made here for both described and previously unknown taxa. Investigating evolutionary relationships within and among members of this clade of proposed Gondwanan orthoclad taxa has demonstrated that a complex interaction between historical population fragmentation and dispersal at several levels of evolutionary scale has been important in driving diversification in this group. While interruptions to migration, colonisation and gene flow driven by population fragmentation have clearly contributed to the development and maintenance of much of the diversity present in this group, long-distance dispersal has also played a role in influencing diversification of continental biotas and facilitating gene flow among disjunct populations.

Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer

Background Phylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction. Results Bayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallace's Line, a biogeographic boundary between Southeast Asia and Australia. Conclusion We describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.

A new geological model of Cenozoic New Zealand suggests a complex phylogeographic history for New Zealand terrestrial taxa

The new model of North Island Cenozoic palaeogeography developed by Kamp et al. has a range of important implications for the evolution of New Zealand terrestrial taxa over the past 30 Ma. Key aspects include the prolonged isolation of the biota on the North Island landmass from the larger and more diverse greater South Island, and the founding of North Island taxa from the potentially unusual ecosystem of a small island around Northland. The prolonged period of isolation is expected to have generated deep phylogenetic splits within taxa present on both islands, and an important current aim should be to identify such signals in surviving endemics to start building a picture of the historical phylogeography, and inferred ecology of both islands through the Cenozoic. Given the potential differences in founding terrestrial species and climatic conditions, it seems likely that the ecology may have been very diferent between the North and South Islands. New genetic data from the 10 or so species of extinct moa suggest that the radiation of moa was much more recent than previously suggested, and reveals a complex pattern that is inferred to result from the interplay of the Cenozoic biogeography, marine barriers, and glacial cycles.

Taxonomy and palaeoecology of Quaternary faunas from caves in eastern tropical Queensland : a record of broad-scale environmental change

Analysis of fossils from cave deposits at Mount Etna (eastern-central Queensland) has established that a species-rich rainforest palaeoenvironment existed in that area during the middle Pleistocene. This unexpected finding has implications for several fields (e.g., biogeography/phylogeography of rainforest-adapted taxa, and the impact of climate change on rainforest communities), but it was unknown whether the Mount Etna sites represented a small refugial patch of rainforest or was more widespread. In this study numerous bone deposits in caves in north-east Queensland are analysed to reconstruct the environmental history of the area during the late Quaternary. Study sites are in the Chillagoe/Mitchell Palmer and Broken River/Christmas Creek areas. The cave fossil records in these study areas are compared with dated (middle Pleistocene-Holocene) cave sites in the Mount Etna area. Substantial taxonomic work on the Mount Etna faunas (particularly dasyurid marsupials and murine rodents) is also presented as a prerequisite for meaningful comparison with the study sites further north. Middle Pleistocene sites at Mount Etna contain species indicative of a rainforest palaeoenvironment. Small mammal assemblages in the Mount Etna rainforest sites (>500-280 ka) are unexpectedly diverse and composed almost entirely of new species. Included in the rainforest assemblages are lineages with no extant representatives in rainforest (e.g., Leggadina), one genus previously known only from New Guinea (Abeomelomys), and forms that appear to bridge gaps between related but morphologically-divergent extant taxa ('B-rat' and 'Pseudomys C'). Curiously, some taxa (e.g., Melomys spp.) are notable for their absence from the Mount Etna rainforest sites. After 280 ka the rainforest faunas are replaced by species adapted to open, dry habitats. At that time the extinct ‘rainforest’ dasyurids and rodents are replaced by species that are either extant or recently extant. By the late Pleistocene all ‘rainforest’ and several ‘dry’ taxa are locally or completely extinct, and the small mammal fauna resembles that found in the area today. The faunal/environmental changes recorded in the Mount Etna sites were interpreted by previous workers as the result of shifts in climate during the Pleistocene. Many samples from caves in the Chillagoe/Mitchell-Palmer and Broken River/Christmas Creek areas are held in the Queensland Museum’s collection. These, supplemented with additional samples collected in the field as well as samples supplied by other workers, were systematically and palaeoecologically analysed for the first time. Palaeoecological interpretation of the faunal assemblages in the sites suggests that they encompass a similar array of palaeoenvironments as the Mount Etna sites. ‘Rainforest’ sites at the Broken River are here interpreted as being of similar age to those at Mount Etna, suggesting the possibility of extensive rainforest coverage in eastern tropical Queensland during part of the Pleistocene. Likewise, faunas suggesting open, dry palaeoenvironments are found at Chillagoe, the Broken River and Mount Etna, and may be of similar age. The 'dry' faunal assemblage at Mount Etna (Elephant hole Cave) dates to 205-170 ka. Dating of one of the Chillagoe sites (QML1067) produced a maximum age for the deposit of approximately 200 ka, and the site is interpreted as being close to that age, supporting the interpretation of roughly contemporaneous deposition at Mount Etna and Chillagoe. Finally, study sites interpreted as being of late Pleistocene-Holocene age show faunal similarities to sites of that age near Mount Etna. This study has several important implications for the biogeography and phylogeography of murine rodents, and represents a major advance in the study of the Australian murine fossil record. Likewise the survey of the northern study areas is the first systematic analysis of multiple sites in those areas, and is thus a major contribution to knowledge of tropical Australian faunas during the Quaternary. This analysis suggests that climatic changes during the Pleistocene affected a large area of eastern tropical Queensland in similar ways. Further fieldwork and dating is required to properly analyse the geographical extent and timing of faunal change in eastern tropical Queensland.