The ancient lakes of Indonesia: towards integrated research on speciation

Ancient lakes have provided considerable insights into the drivers of speciation and adaptive radiation in aquatic organisms. Most studies of species-flocks, however, focus only on a single group of organisms, and few have attempted to integrate geological, limnological, ecological, and genetic drivers of speciation on multiple species-flocks at various trophic levels. As such, there is a need for a comprehensive model system for research on speciation in aquatic environments where multiple radiations are investigated at various levels of biological organization (e.g., individual, population, and ecosystem) and placed in light of geographical and geological setting. The ancient Malili Lakes of Sulawesi, Indonesia, are ideal candidates for such a model, and represent the only hydrologically connected ancient lakes in the world. These lakes are characterized by ultra-oligotrophy and unique physicochemical conditions that govern the composition and production of planktonic communities. At higher trophic levels, there are three recurring trends: (1) low taxonomic richness and simple community structures, (2) adaptive radiations with trophic specialization, and (3) remarkably high endemism with evolutionary innovations throughout the lakes and species-flocks. Furthermore, the restricted geographic distributions of species-flocks within the Malili Lakes indicate that each lake constitutes a unique environment, and dispersal among lakes is limited, despite close contemporary connectivity. These observations suggest that ecological and evolutionary processes are regulated from the bottom up, and speciation is primarily facilitated by interspecific and intraspecific competition for limited resources. The Malili Lakes represent an outstanding natural model for integrative research into speciation as they offer the opportunity to explore the roles of geography, dispersal, and selection in the radiation of aquatic organisms.     

Ancient lakes as hotspots of diversity: a morphological review of an endemic species flock of <Emphasis Type="Italic">Tylomelania</Emphasis> (Gastropoda: Cerithioidea: Pachychilidae) in the Malili lake system on Sulawesi,Indonesia

The viviparous freshwater gastropod Tylomelania (Caenogastropoda: Cerithioidea: Pachychilidae) endemic to the Indonesian island Sulawesi has radiated extensively in two ancient lake systems. We here present the first systematic species-level review of taxa in the five lakes of the Malili lake system, which contains the most diverse and best studied freshwater fauna on Sulawesi. Our results indicate a significantly higher diversity of Tylomelania in these lakes than previously perceived based on morphological evidence for delimiting the taxa. We describe nine new species, thus increasing the number of taxa known from the Malili lakes to 25. Tylomelania species are inhabiting all available substrates in the lakes, and the diversity of habitats is reflected in an unparalleled range of radula types in this closely related group. Several species show a high intraspecific variability in some characters, and their closer investigation will probably lead to the discovery of more cryptic species. As it is, this species flock on Sulawesi is among the largest freshwater mollusc radiations known. Since the Malili lake system also contains other large endemic species flocks of e.g. crustaceans and fishes, it is a major hotspot of freshwater biodiversity in Asia to become a conservation priority. Handling editor: K. Martens     

ADAPTIVE RADIATION AND ECOLOGICAL DIVERSIFICATION OF SULAWESI'S ANCIENT LAKE SHRIMPS

Ancient lakes are natural laboratories for the study of adaptive radiation. Recently, two lake systems on the Indonesian island Sulawesi have emerged as promising new model systems. A species flock of atyid freshwater shrimps in the Malili lake system comprises 15 colorful endemic taxa. Mitochondrial DNA data suggest two independent colonizations by riverine ancestors. Only one colonization event led to subsequent radiation into 14 species, while the second clade comprises just one species. The vast majority of species (n= 12) are habitat specialists, which are confined to the larger Malili clade and include all taxa with species‐specific color patterns and a restricted distribution within the five connected Malili lakes. Morphological, genetic, and ecological data are consistent with the existence of an adaptive radiation in the Malili lakes, involving the habitat‐specific diversification of trophic morphology. In addition to testing criteria for the recognition of an adaptive radiation, an ancestral state reconstruction reveals an equal probability for either a riverine generalist or a lacustrine specialist as ancestor of the large Malili clade, which is interpreted as indicative of an early stage of habitat specialization within this radiation. Finally, our results suggest that species diversification may have been primarily driven by ecological specialization and allopatric speciation.     

Historic hybridization and introgression between two iconic Australian anemonefish and contemporary patterns of population connectivity

Endemic species on islands are considered at risk of extinction for several reasons, including limited dispersal abilities, small population sizes, and low genetic diversity. We used mitochondrial DNA (D-Loop) and 17 microsatellite loci to explore the evolutionary relationship between an endemic anemonefish, Amphiprion mccullochi (restricted to isolated locations in subtropical eastern Australia) and its more widespread sister species, A. akindynos. A mitochondrial DNA (mtDNA) phylogram showed reciprocal monophyly was lacking for the two species, with two supported groups, each containing representatives of both species, but no shared haplotypes and up to 12 species, but not location-specific management units (MUs). Population genetic analyses suggested evolutionary connectivity among samples of each species (mtDNA), while ecological connectivity was only evident among populations of the endemic, A. mccullochi. This suggests higher dispersal between endemic anemonefish populations at both evolutionary and ecological timeframes, despite separation by hundreds of kilometers. The complex mtDNA structure results from historical hybridization and introgression in the evolutionary past of these species, validated by msat analyses (NEWHYBRIDS, STRUCTURE, and DAPC). Both species had high genetic diversities (mtDNA h > 0.90, π = 4.0%; msat genetic diversity, gd > 0.670). While high gd and connectivity reduce extinction risk, identifying and protecting populations implicated in generating reticulate structure among these species should be a conservation priority.     

Genetic affinities of an introduced predator: Nile perch in Lake Victoria, East Africa

Although the introduction of Nile perch, Lates niloticus , to Lake Victoria has received intense global attention, especially in relation to its impact on endemic cichlid species and on fishery yields, fundamental information on its taxonomy and population genetics is lacking. Most importantly, the introduced fish originated from two lakes (Lakes Albert and Turkana) containing three Lates species, and it has never been entirely clear which of these became established in Lake Victoria, or indeed whether the Lake Victoria population is derived from hybridization between Lates species. In addition, genetic drift caused by the relatively small founder population (≈ 400), the initially slow population increase followed by a period of explosive population growth, and selection pressures in the new environment may have resulted in substantial genetic changes. Allozyme data indicated that the introduced Nile perch of Lake Victoria were mainly L. niloticus from Lake Albert, although maximum likelihood estimates of stock contributions (GSI) suggested the presence of L. macrophthalmus. In contrast, introduced Nile perch in adjacent smaller lakes (Lakes Kyoga and Nabugabo) appeared to be entirely L. niloticus . The effect of the introductions on allozyme diversity varied among lakes and appeared to be uncorrelated to the number of fish introduced.     

Fitness and genetic variation of Viola calaminaria, an endemic metallophyte: implications of population structure and history

We investigated variations in genetic diversity and plant fitness in a rare endemic metallophyte of calamine soils, Viola calaminaria, in relation to population size, population connectivity and population history in order to evaluate and discuss potential conservation strategies for the species. Mean population genetic diversity (H(s) = 0.25) of V. calaminaria was similar to endemic non-metallophyte taxa. Twenty-one per cent of the genetic variation was partitioned among populations and a low (9%) but significant differentiation was found among geographical regions. Our results did not support the hypothesis that the acquisition of metal tolerance may result in reduced genetic diversity, and suggested that strict metallophytes do not exhibit higher inter-population differentiation resulting from scattered habitats. There were no relationships between population genetic diversity and population size. Significant correlations were found between plant fitness and (i) population size and (ii) connectivity index. Recently-founded populations exhibited the same level of genetic diversity as ancient populations and also possessed higher plant fitness. There was no indication of strong founder effects in recently-established populations. The results suggest that the creation of habitats through human activities could provide new opportunities for conservation of this species.     

Predicting global population connectivity and targeting conservation action for snow leopard across its range

Movements of individuals within and among populations help to maintain genetic variability and population viability. Therefore, understanding landscape connectivity is vital for effective species conservation. The snow leopard is endemic to mountainous areas of central Asia and occurs within 12 countries. We assess potential connectivity across the species’ range to highlight corridors for dispersal and genetic flow between populations, prioritizing research and conservation action for this wide‐ranging, endangered top‐predator. We used resistant kernel modeling to assess snow leopard population connectivity across its global range. We developed an expert‐based resistance surface that predicted cost of movement as functions of topographical complexity and land cover. The distribution of individuals was simulated as a uniform density of points throughout the currently accepted global range. We modeled population connectivity from these source points across the resistance surface using three different dispersal scenarios that likely bracket the lifetime movements of individual snow leopard: 100 km, 500 km and 1000 km. The resistant kernel models produced predictive surfaces of dispersal frequency across the snow leopard range for each distance scenario. We evaluated the pattern of connectivity in each of these scenarios and identified potentially important movement corridors and areas where connectivity might be impeded. The models predicted two regional populations, in the north and south of the species range respectively, and revealed a number of potentially important connecting areas. Discrepancies between model outputs and observations highlight unsurveyed areas of connected habitat that urgently require surveying to improve understanding of the global distribution and ecology of snow leopard, and target land management actions to prevent population isolation. The connectivity maps provide a strong basis for directed research and conservation action, and usefully direct the attention of policy makers.     

Watershed characteristics shape the landscape genetics of brook stickleback (Culaea inconstans) in shallow prairie lakes

Investigating the consequences of landscape features on population genetic patterns is increasingly important to elucidate the ecological factors governing connectivity between populations and predicting the evolutionary consequences of landscapes. Small prairie lakes in Alberta, Canada, and the brook stickleback (Culaea inconstans) that inhabit them, provide a unique aquatic system whereby populations are highly isolated from one another. These heterogeneous and extreme environments are prone to winterkills, an event whereby most of the fish die and frequent bottlenecks occur. In this study, we characterized the genetic population structure of brook stickleback among several lakes, finding that the species is hierarchically influenced by within‐lake characteristics in small‐scale watersheds. Landscape genetic analyses of the role of spatial features found support for basin characteristics associated with genetic diversity and bottlenecks in 20% of the sampled lakes. These results suggest that brook stickleback population genetic patterns may be driven, at least in part, by ecological processes that accelerate genetic drift and landscape patterns associated with reduced dispersal. Collectively, these results reinforce the potential importance of connectivity in the maintenance of genetic diversity, especially in fragmented landscapes.     

Montane refugia predict population genetic structure in the Large‐blotched Ensatina salamander

Understanding the biotic consequences of Pleistocene range shifts and fragmentation remains a fundamental goal in historical biogeography and evolutionary biology. Here, we combine species distribution models (SDM) from the present and two late Quaternary time periods with multilocus genetic data (mitochondrial DNA and microsatellites) to evaluate the effect of climate‐induced habitat shifts on population genetic structure in the Large‐blotched Ensatina (Ensatina eschscholtzii klauberi), a plethodontid salamander endemic to middle and high‐elevation conifer forest in the Transverse and Peninsular Ranges of southern California and northern Baja California. A composite SDM representing the range through time predicts two disjunct refugia, one in southern California encompassing the core of the species range and the other in the Sierra San Pedro Mártir of northern Baja California at the southern limit of the species range. Based on our spatial model, we would expect a pattern of high connectivity among populations within the northern refugium and, conversely, a pattern of isolation due to long‐term persistence of the Sierra San Pedro Mártir population. Our genetic results are consistent with these predictions based on the hypothetical refugia in that (i) historical measures of population connectivity among stable areas are correlated with gene flow estimates; and (ii) there is strong geographical structure between separate refugia. These results provide evidence for the role of recent climatic change in shaping patterns of population persistence and connectivity within the Transverse and Peninsular Ranges, an evolutionary hotspot.     

Directional gene flow patterns in disjunct populations of the black ratsnake (Pantheropis obsoletus) and the Blanding’s turtle (Emydoidea blandingii)

The estimation and maintenance of connectivity among local populations is an important conservation goal for many species at risk. We used Bayesian statistics and coalescent theory to estimate short- and long-term directional gene flow among subpopulations for two reptiles that occur in Canada as peripheral populations that are geographically disjunct from the core of their respective species’ ranges: the black ratsnake and the Blanding’s turtle. Estimates of directional gene flow were used to examine population connectivity and potential genetic source-sink dynamics. For both species, our estimates of directional short- and long-term gene flow were consistently lower than estimates inferred previously from F _ST measures. Short- and long-term gene flow estimates were discordant in both species, suggesting that population dynamics have varied temporally in both species. These estimates of directional gene flow were used to identify specific subpopulations in both species that may be of high conservation value because they are net exporters of individuals to other subpopulations. Overall, our results show that the use of more sophisticated methods to evaluate population genetic data can provide valuable information for the conservation of species at risk, including bidirectional estimates of subpopulation connectivity that rely on fewer assumptions than more traditional analyses. Such information can be used by conservation practitioners to better understand the geographic scope required to maintain a functional metapopulation, determine which habitat corridors within a working landscape may be most important to maintain connectivity among subpopulations, and to prioritize subpopulations with respect to their potential to act as genetic sources within the metapopulation.     

A complex history of introgression and vicariance in a threatened montane skink (<Emphasis Type="Italic">Pseudemoia cryodroma</Emphasis>) across an Australian sky island system

Species endemic to sky island systems are isolated to mountain peaks and high elevation plateaux both geographically and ecologically, making them particularly vulnerable to the effects of climate change. Pressures associated with climate change have already been linked to local extinctions of montane species, emphasizing the importance of understanding the genetic diversity and population connectivity within sky islands systems for the conservation management of remaining populations. Our study focuses on the endangered alpine skink Pseudemoia cryodroma, which is endemic to the Victorian Alps in south-eastern Australia, and has a disjunct distribution in montane habitats above 1100 m a.s.l. Using mitochondrial DNA (mtDNA) and microsatellite loci, we investigated species delimitation, genetic connectivity and population genetic structure across the geographic range of this species. We found discordance between genetic markers, indicating historical mtDNA introgression at one of the study sites between P. cryodroma and the closely related, syntopic P. entrecasteauxii. Molecular diversity was positively associated with site elevation and extent of suitable habitat, with inbreeding detected in three of the five populations. These results demonstrate the complex interaction between geography and habitat in shaping the population structure and genetic diversity of P. cryodroma, and highlight the importance of minimising future habitat loss and fragmentation for the long-term persistence of this species.     

The icthyofauna of the Mhlathuze coastal lakes: some preliminary results

The fish faunas of the four Mhlathuze coastal lakes and the lower river comprise a diverse assortment of over fifty marine, estuarine and freshwater species. Three freshwater species are endemic to KwaZulu-Natal and nine estuarine species are endemic to southern Africa. Five species are of conservation significance. Species numbers in Lakes Mzingazi and Cubhu are similar historically and both lakes served as secondary nursery habitats for estuarine associated fishes. This role has been impacted by the construction of weirs at their outlets which prevent successful recruitment of estuarine species, especially during drought years when lake water levels are low. The fish faunas of Lakes Nsezi and Mangeza are depauperate and lack marine or estuarine components. In order that these systems fulfil their potential function as secondary nursery habitats to many estuarine fish species, minimum lake water levels must be set to ensure sufficient outflows at proposed fish ladders during critical spawning and migration times.     

Genetic connectivity and diversity in inselberg populations of Acacia woodmaniorum,a rare endemic of the Yilgarn Craton banded iron formations

Historically rare plant species with disjunct population distributions and small population sizes might be expected to show significant genetic structure and low levels of genetic diversity because of the effects of inbreeding and genetic drift. Across the globe, terrestrial inselbergs are habitat for rich, often rare and endemic flora and are valuable systems for investigating evolutionary processes that shape patterns of genetic structure and levels of genetic diversity at the landscape scale. We assessed genetic structure and levels of genetic diversity across the range of the historically rare inselberg endemic Acacia woodmaniorum. Phylogeographic and genetic structure indicates that connectivity is not sufficient to produce a panmictic population across the limited geographic range of the species. However, historical levels of gene flow are sufficient to maintain a high degree of adaptive connectivity across the landscape. Genetic diversity indicates gene flow is sufficient to largely counteract any negative genetic effects of inbreeding and random genetic drift in even the most disjunct or smallest populations. Phylogeographic and genetic structure, a signal of isolation by distance and a lack of evidence of recent genetic bottlenecks suggest long-term stability of contemporary population distributions and population sizes. There is some evidence that genetic connectivity among disjunct outcrops may be facilitated by the occasional long distance dispersal of Acacia polyads carried by insect pollinators moved by prevailing winds.     

The species flocks of lacustrine gastropods: <Emphasis Type="Italic">Tylomelania</Emphasis> on Sulawesi as models in speciation and adaptive radiation

Endemic radiations provide splendid opportunities for studies in evolutionary biology. Species flocks in ancient lakes, such as in Tanganyika, Malawi or Baikal, have featured prominently in evolutionary biology, viewing these “evolutionary theatres” as hotspots of diversification. However, following a century of neglect, the endemic evolution of limnic cerithioidean gastropods in the two central lake systems on the Indonesian island of Sulawesi (i.e. Lake Poso and the lakes of the Malili system, e.g. Danau Matano, Mahalona and Towuti) also provide instructive model cases for the study of speciation mechanisms, adaptive radiation and annidation (i.e. niche exploitation). We here discuss the evolutionary and taxonomic implications of the lacustrine species flocks in Tylomelania from these lakes in Sulawesi as an exceptional endemic assemblage of morphologically distinct viviparous pachychilid gastropods. This first comprehensive compilation of data on both ancient lake systems, Poso and Malili, offers a new perspective on ecological differentiation in this radiation. Presented here within the framework of the theory of evolutionary ecology it provides a research program for acquiring a synthetical perspective that includes morphology, molecular genetics, ecology and biogeography. In this context, it will be possible to compare the species flocks of these truly “Darwinian snails” on Sulawesi with the long enigmatic, so-called thalassoid (i.e. marine-like) gastropod radiation in East African’s Lake Tanganyika.     

Gene flow at the margin of Lake Matano’s adaptive sailfin silverside radiation: Telmatherinidae of River Petea in Sulawesi

Classical speciation concepts focus almost exclusively on the evolution of strict reproductive isolation as a prerequisite for speciation. However, there is a growing body of evidence indicating that speciation is possible despite or even triggered by gene flow among populations or species. Previous findings indicate that introgressive hybridization is a dominant phenomenon in the adaptive radiation of sailfin silversides (Telmatherinidae) endemic to Lake Matano (Sulawesi). In this study, we investigate patterns of genotypic and phenotypic variation of “sharpfin” sailfin silversides in the outlet area of L. Matano and six locations along River Petea, which is the only connection between L. Matano and other lakes and streams of the Malili Lakes system. Fieldwork revealed no hints for a previously cited major waterfall in River Petea, which was thought to separate L. Matano’s sailfin silverside radiation from the diversity of the downstream lake drainages. Likewise, genomic (AFLP) and morphometric data suggest high levels of gene flow between upper and lower stretches of this river, as well as between riverine Petea and lacustrine Matano populations. Increasing levels of genotypic and phenotypic dissimilarity are correlated with distance over a remarkably short geographic range.     

Brood care among basommatophorans: a unique reproductive strategy in the freshwater limpet snail Protancylus (Heterobranchia: Protancylidae), endemic to ancient lakes on Sulawesi, Indonesia

In molluscan taxa inhabiting marine environments oviparity and reproduction via planktonic larvae is predominant while incubation and viviparity is most frequently found in taxa inhabiting brackish or freshwater aquatic habitats. Brooding has evolved repeatedly and independently in several limnic taxa among Bivalvia and Gastropoda. However, among basommatophoran gastropods no such cases were yet known. We here report on a unique reproductive strategy involving brood care in the lacustrine freshwater limpet genus Protancylus, endemic to the ancient lakes on central Sulawesi (former Celebes), Indonesia, namely the Lake Poso and the Malili lake system, because this constitutes the first known case of this behaviour among the Basommatophora. Protancylus live exclusively as epizoans on those pachychilid gastropods of the viviparous genus Tylomelania, also a Sulawesi endemic species, that inhabit mostly soft substrates. We found that the two known species Protancylus pileolus from Lake Poso and P. adhaerens from the Malili lake system both retain gelatinous egg strings underneath their outer mantle, where up to 15 (mostly eight or nine) shelled juveniles are brooded. Nourishment is provided within the egg capsule only. Thus, brood care in Protancylus resembles the reproductive strategy found recently among pachychilid gastropods Jagora from the Philippines, but differs from euviviparous (i.e. matrotrophic) incubation among thiarid gastropods possessing a brood pouch with juveniles being nourished via a ‘pseudoplacenta’ in several taxa.     

Migratory behaviour shapes spatial genetic structure of cyprinid fishes within the Lake Malawi catchment

• Genetic differences among freshwater fish populations are dependent on life‐history characteristics of the species, including the range of adult dispersal and the extent of homing to natal breeding grounds. However, the effects of variation in such characteristics on population genetic connectivity are rarely studied comparatively among closely related species.
• We studied population genetic structure within three congeneric cyprinid species from the Lake Malawi catchment that differ substantially in life‐history traits and conservation status, using a combination of microsatellite and mitochondrial DNA markers. Mpasa (Opsaridium microlepis) is a large (70 cm total length) migratory species that spawns in rivers, but as an adult is exclusively known from the main lake body. Sanjika (Opsaridium microcephalum), is a medium size (30 cm total length) species that exists in lake breeding, river‐lake migratory and apparently landlocked populations. Dwarf sanjika (Opsaridium tweddleorum) is a small non‐migratory species (15 cm total length) that persists in small tributaries surrounding the main lake and adjoining rivers.
• The results revealed striking differences among the three species in spatial genetic structuring. The river‐lake migratory mpasa showed only weak yet significant population genetic structure within the main Lake Malawi catchment, suggesting that there is no strong natal homing. The habitat‐generalist sanjika showed only weak spatial genetic differentiation at microsatellite loci within the Lake Malawi catchment, but moderate structure in mitochondrial DNA, potentially reflecting male‐biased dispersal. The river‐restricted dwarf sanjika showed strong genetic structure in both microsatellite and mitochondrial DNA, suggesting strictly limited dispersal at both adult and juvenile stages.
• We conclude that contrasting migration life histories have resulted in dramatically different patterns of population genetic structure among these congeneric species. The observed patterns demonstrate how divergent life‐history evolution may strongly influence broader patterns of population genetic connectivity in freshwater fish, with consequences for management and conservation. Specifically the results suggesting gene flow among Lake Malawi populations of mpasa, an IUCN red‐listed ‘Endangered’ species endemic to the lake catchment, imply that conservation initiatives operating at both local and catchment scales are needed to reverse local population decline.

     

Water availability strongly impacts population genetic patterns of an imperiled Great Plains endemic fish

Genetic, demographic, and environmental processes affect natural populations synergistically, and understanding their interplay is crucial for the conservation of biodiversity. Stream fishes in metapopulations are particularly sensitive to habitat fragmentation because persistence depends on dispersal and colonization of new habitat but dispersal is constrained to stream networks. Great Plains streams are increasingly fragmented by water diversion and climate change, threatening connectivity of fish populations in this ecosystem. We used seven microsatellite loci to describe population and landscape genetic patterns across 614 individuals from 12 remaining populations of Arkansas darter (Etheostoma cragini) in Colorado, a candidate species for listing under the U.S. Endangered Species Act. We found small effective population sizes, low levels of genetic diversity within populations, and high levels of genetic structure, especially among basins. Both at- and between-site landscape features were associated with genetic diversity and connectivity, respectively. Available stream habitat and amount of continuous wetted area were positively associated with genetic diversity within a site, while stream distance and intermittency were the best predictors of genetic divergence among sites. We found little genetic contribution from historic supplementation efforts, and we provide a set of management recommendations for this species that incorporate a conservation genetics perspective.