Molecular phylogeny and phylogeography of the Australian freshwater fish genus Galaxiella, with an emphasis on dwarf galaxias (G. pusilla)

Cryptic diversity within two widespread diadromous freshwater fishes (Teleostei: Galaxiidae)

Identification of taxonomically cryptic species is essential for the effective conservation of biodiversity. Freshwater-limited organisms tend to be genetically isolated by drainage boundaries, and thus may be expected to show substantial cryptic phylogenetic and taxonomic diversity. By comparison, populations of diadromous taxa, that migrate between freshwater and marine environments, are expected to show less genetic differentiation. Here we test for cryptic diversity in Australasian populations (both diadromous and non-diadromous) of two widespread Southern Hemisphere fish species, Galaxias brevipinnis and Galaxias maculatus. Both mtDNA and nuclear markers reveal putative cryptic species within these taxa. The substantial diversity detected within G. brevipinnis may be explained by its strong climbing ability which allows it to form isolated inland populations. In island populations, G. brevipinnis similarly show deeper genetic divergence than those of G. maculatus, which may be explained by the greater abundance of G. maculatus larvae in the sea allowing more ongoing dispersal. Our study highlights that even widespread, 'high-dispersal' species can harbour substantial cryptic diversity and therefore warrant increased taxonomic and conservation attention.

Same process, different patterns: pervasive effect of evolutionary time on species richness in freshwater fishes

Tropical lands harbour the highest number of species, resulting in the ubiquitous latitudinal diversity gradient (LDG). However, exceptions to this pattern have been observed in some taxa, explained by the interaction between the evolutionary histories and environmental factors that constrain species' physiological and ecological requirements. Here, we applied a deconstruction approach to map the detailed species richness patterns of Actinopterygian freshwater fishes at the class and order levels and to disentangle their drivers using geographical ranges and a phylogeny, comprising 77% (12 557) of all described species. We jointly evaluated seven evolutionary and ecological hypotheses posited to explain the LDG: diversification rate, time for speciation, species-area relationship, environmental heterogeneity, energy, temperature seasonality and past temperature stability. We found distinct diversity gradients across orders, including expected, bimodal and inverse LDGs. Despite these differences, the positive effect of evolutionary time explained patterns for all orders, where species-rich regions are inhabited by older species compared to species-poor regions. Overall, the LDG of each order has been shaped by a unique combination of factors, highlighting the importance of performing a joint evaluation of evolutionary, historical and ecological factors at different taxonomic levels to reach a comprehensive understanding on the causes driving global species richness patterns.

Geographic range size and speciation in honeyeaters

Background

Darwin and others proposed that a species’ geographic range size positively influences speciation likelihood, with the relationship potentially dependent on the mode of speciation and other contributing factors, including geographic setting and species traits. Several alternative proposals for the influence of range size on speciation rate have also been made (e.g. negative or a unimodal relationship with speciation). To examine Darwin’s proposal, we use a range of phylogenetic comparative methods, focusing on a large Australasian bird clade, the honeyeaters (Aves: Meliphagidae).

Results

We consider the influence of range size, shape, and position (latitudinal and longitudinal midpoints, island or continental species), and consider two traits known to influence range size: dispersal ability and body size. Applying several analytical approaches, including phylogenetic Bayesian path analysis, spatiophylogenetic models, and state-dependent speciation and extinction models, we find support for both the positive relationship between range size and speciation rate and the influence of mode of speciation.

Conclusions

Honeyeater speciation rate differs considerably between islands and the continental setting across the clade’s distribution, with range size contributing positively in the continental setting, while dispersal ability influences speciation regardless of setting. These outcomes support Darwin’s original proposal for a positive relationship between range size and speciation likelihood, while extending the evidence for the contribution of dispersal ability to speciation.

Geography vs. past climate: the drivers of population genetic structure of the Himalayan langur

Background

Contemporary species distribution, genetic diversity and evolutionary history in many taxa are shaped by both historical and current climate as well as topography. The Himalayas show a huge variation in topography and climatic conditions across its entire range, and have experienced major climatic fluctuations in the past. However, very little is known regarding how this heterogenous landscape has moulded the distribution of Himalayan fauna. A recent study examined the effect of these historical events on the genetic diversity of the Himalayan langurs in Nepal Himalaya. However, this study did not include the samples from the Indian Himalayan region (IHR). Therefore, here we revisit the questions addressed in the previous study with a near complete sampling from the IHR, along with the samples from the Nepal Himalaya. We used the mitochondrial Cytochrome-b (Cyt-b, 746 bp) region combined with multiple phylogeographic analyses and palaeodistribution modelling.

Results

Our dataset contained 144 sequences from the IHR as well as the Nepal Himalaya. Phylogenetic analysis showed a low divergent western clade nested within high divergent group of eastern lineages and in the network analysis we identified 22 haplotypes over the entire distribution range of the Himalayan langurs. Samples from the Nepal Himalaya showed geographically structured haplotypes corresponding to different river barriers, whereas samples from IHR showed star-like topology with no structure. Our statistical phylogeography analysis using diyABC supported the model of east to west colonisation of these langurs with founder event during colonisation. Analysis of demographic history showed that the effective population size of the Himalayan langurs decreased at the onset of last glacial maximum (LGM) and started increasing post LGM. The palaeodistribution modelling showed that the extent of suitable habitat shifted from low elevation central Nepal, and adjoining parts of north India, during LGM to the western Himalaya at present.

Conclusion

The current genetic diversity and distribution of Himalayan langurs in the Nepal Himalaya has been shaped by river barriers, whereas the rivers in the IHR had relatively less time to act as a strong genetic barrier after the recent colonisation event. Further, the post LGM expansion could have had confounding effect on Himalayan langur population structure in both Nepal Himalaya and IHR.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02054-1.

Insight into the evolutionary history of Indoplanorbis exustus (Bulinidae: Gastropoda) at the scale of population and species

The history of a lineage is intertwined with the history of the landscape it inhabits. Here we showcase how the geo-tectonic and climatic evolution of South Asia and surrounding landmasses have shaped the biogeographical history of Indoplanorbis exustus, a tropical Asian, freshwater snail. We amplified partial COI gene fragments from all over India and combined this with a larger dataset from South and Southeast Asia to carry out phylogenetic reconstruction, species delimitation analysis and population genetic analyses. Two nuclear genes were also amplified from a few individuals per putative species to carry out divergence dating and ancestral area reconstruction analyses. The results suggest that I. exustus dispersed out of Africa into India during the Eocene. Furthermore, molecular data suggest I. exustus is a species complex consisting of multiple putative species. Primary diversification took place in the Northern Indian plains or in Northeast India. The speciation events appear to be primarily allopatric caused by a series of aridification events starting from the late Miocene to early Pleistocene. None of the species appears to have any underlying genetic structure suggestive of high vagility. All the species underwent population fluctuations during the Pleistocene, probably driven by the Quaternary climatic fluctuations.

Miocene biome turnover drove conservative body size evolution across Australian vertebrates

On deep time scales, changing climatic trends can have a predictable influence on macroevolution. From evidence of mass extinctions, we know that rapid climatic oscillations can indirectly open niche space and precipitate adaptive radiation, changing the course of ecological diversification. These dramatic shifts in the global climate, however, are rare events relative to extended periods of protracted climate change and biome turnover. It remains unclear whether during gradually changing periods, shifting habitats may instead promote non-adaptive speciation by facilitating allopatry and phenotypic conservatism. Using fossil-calibrated, species-level phylogenies for five Australian radiations comprising more than 800 species, we investigated temporal trends in biogeography and body size evolution. Here, we demonstrate that gradual Miocene cooling and aridification correlates with the restricted phenotypic diversification of multiple ecologically diverse vertebrate groups. This probably occurred as species ranges became fractured and isolated during continental biome restructuring, encouraging a shift towards conservatism in body size evolution. Our results provide further evidence that abiotic changes, not only biotic interactions, may act as selective forces influencing phenotypic macroevolution.

Molecular systematics, species delimitation and diversification patterns of the Phyllodactylus lanei complex (Gekkota: Phyllodactylidae) in Mexico

The description of cryptic gecko species worldwide has revealed both that many putative species are, in fact, conformed by a complex of morphologically conserved species that are genetically distinct and highly divergent, and that gecko species diversity could be underestimated. The taxonomy and species delimitation of geckos belonging to the genus Phyllodactylus is still controversial, 16 of which are distributed in Mexico and 13 are endemic. Although the large morphological variation shown by the Phyllodactylus species from Mexico has been amply documented, little is known about their genetic diversity and evolutionary relationships, and much less regarding cryptic speciation. Here, we included the most comprehensive sampling of populations and species of the Phyllodactylus lanei complex distributed in Mexico, and applied an analytical approach that included probabilistic phylogenetic analyses, jointly with species delimitation methods and Bayesian putative species validation analysis. Our results suggest the existence of 10 lineages within the complex, supporting the existence of cryptic species, and in great contrast with the current taxonomic proposal that includes only four subspecies. The most recent common ancestor (MRCA) for the P. lanei clade originated on the Early Eocene (∼54Mya), along the southern coasts of Mexico, followed by the highest diversification of the complex MRCA during the Eocene (34-56Mya). Lineages subsequently dispersed and diversified towards the northwest, and the diversification process ended with the most recent lineages inhabiting two islands on the coasts of Nayarit (Miocene; 5.5-23Mya). Our results highlight three vicariant events associated with the evolution of the lineages, two of them intimately related to the formation of the Sierra Madre del Sur and the Transmexican Volcanic Belt mountain ranges, main geographic barriers that isolated and facilitated the divergence and speciation in this group of geckos. Finally, we propose that there are 10 species in the P. lanei complex, from which four represent taxonomic changes and six are new species and require a formal description. We acknowledge that more analyses, including a detailed evaluation of morphological characters and use of more unlinked nuclear loci with enough variability, are needed to further support their taxonomic description.

Why Australasian vertebrate animals are so unique - A palaeontological perspective

Australasia has a unique fauna of living vertebrates, which include the oldest known species on the planet (the lungfish Neoceratodus) as well as many diverse, highly endemic families of fish, amphibians, reptiles, birds and mammals. The origin of most of the Australian vertebrate fauna has developed from two phases. Firstly, when Australia was subsumed within the greater Gondwana landmass, migration of animals from one region to another was possible by a land connection. Many of our most primitive forms of reptiles and mammals probably entered the country at this time, such as varanids, madtsooid snakes, monotremes and basal marsupials. Secondly, following the breakup of Gondwana, the isolation of Australia for its last 40 million years and subsequent changing climatic conditions drove the radiation of marsupial, reptile and amphibian families within the continent. The gradual aridification of central Australia further divided the landmass into discrete regional areas characterised by rainfall, vegetation, and climatic zones.

DNA barcoding and genetic diversity analyses of fishes of Kaladan River of Indo-Myanmar biodiversity hotspot

Species are considered as a fundamental unit of biodiversity. Therefore, the prerequisite for biodiversity management and conservation is to know the number of species one is dealing with. Consequently, the need of the present study was conceptualized, which dealt with the comprehensive molecular appraisal of Kaladan's Fish fauna. A total of 291 specimens representing 49 species, 28 genera, 11 families, and 4 orders, were collected from 11 sampling stations situated along the main Kaladan River and its four major tributaries, i.e. Tiau, Tuipui, Mat, and Tuichang, flowing in Mizoram state of India (part of Indo-Myanmar biodiversity hotspot) and COI sequences of all the 291 samples were generated. All the analyses conducted in the present study, i.e. K2P genetic divergence, bPTP and Neighbour-Joining suggest that DNA Barcoding is an efficient and reliable tool for species identification and deciding the species boundary. Most of the species of Kaladan showed the clear existence of barcode gap. However, the presence of intra-specific and inter-specific genetic distance overlap in two species revealed the existence of putative sibling species or hidden taxa. This study also revealed the presence of two cryptic species and putative previously unknown species of genus Garra and Schistura. The COI barcode database of Kaladan's fish fauna, established in the present study, may serve as a reference library for accurate identification of fishes and will help ichthyologist, researcher, students, biodiversity managers and policy makers in proper planning with regard to conservation and management of the resources.

Phylogeography, historical demography and habitat suitability modelling of freshwater fishes inhabiting seasonally fluctuating Mediterranean river systems: a case study using the Iberian cyprinid Squalius valentinus

The Mediterranean freshwater fish fauna has evolved under constraints imposed by the seasonal weather/hydrological patterns that define the Mediterranean climate. These conditions have influenced the genetic and demographic structure of aquatic communities since their origins in the Mid-Pliocene. Freshwater species in Mediterranean-type climates will likely constitute genetically well-differentiated populations, to varying extents depending on basin size, as a consequence of fragmentation resulting from drought/flood cycles. We developed an integrative framework to study the spatial patterns in genetic diversity, demographic trends, habitat suitability modelling and landscape genetics, to evaluate the evolutionary response of Mediterranean-type freshwater fish to seasonal fluctuations in weather. To test this evolutionary response, the model species used was Squalius valentinus, an endemic cyprinid of the Spanish Levantine area, where seasonal weather fluctuations are extreme, although our findings may be extrapolated to other Mediterranean-type species. Our results underscore the significant role of the Mediterranean climate, along with Pleistocene glaciations, in diversification of S. valentinus. We found higher nuclear diversity in larger drainage basins, but higher mitochondrial diversity correlated to habitat suitability rather than basin size. We also found strong correlation between genetic structure and climatic factors associated with Mediterranean seasonality. Demographic and migration analyses suggested population expansion during glacial periods that also contributed to the current genetic structure of S. valentinus populations. The inferred models support the significant contribution of precipitation and temperature to S. valentinus habitat suitability and allow recognizing areas of habitat stability. We highlight the importance of stable habitat conditions, fostered by typical karstic springs found on the Mediterranean littoral coasts, for the preservation of freshwater species inhabiting seasonally fluctuating river systems.

Assessing species boundaries using multilocus species delimitation in a morphologically conserved group of neotropical freshwater fishes, the Poecilia sphenops species complex (Poeciliidae)

Accurately delimiting species is fundamentally important for understanding species diversity and distributions and devising effective strategies to conserve biodiversity. However, species delimitation is problematic in many taxa, including 'non-adaptive radiations' containing morphologically cryptic lineages. Fortunately, coalescent-based species delimitation methods hold promise for objectively estimating species limits in such radiations, using multilocus genetic data. Using coalescent-based approaches, we delimit species and infer evolutionary relationships in a morphologically conserved group of Central American freshwater fishes, the Poecilia sphenops species complex. Phylogenetic analyses of multiple genetic markers (sequences of two mitochondrial DNA genes and five nuclear loci) from 10/15 species and genetic lineages recognized in the group support the P. sphenops species complex as monophyletic with respect to outgroups, with eight mitochondrial 'major-lineages' diverged by ≥2% pairwise genetic distances. From general mixed Yule-coalescent models, we discovered (conservatively) 10 species within our concatenated mitochondrial DNA dataset, 9 of which were strongly supported by subsequent multilocus Bayesian species delimitation and species tree analyses. Results suggested species-level diversity is underestimated or overestimated by at least ~15% in different lineages in the complex. Nonparametric statistics and coalescent simulations indicate genealogical discordance among our gene tree results has mainly derived from interspecific hybridization in the nuclear genome. However, mitochondrial DNA show little evidence for introgression, and our species delimitation results appear robust to effects of this process. Overall, our findings support the utility of combining multiple lines of genetic evidence and broad phylogeographical sampling to discover and validate species using coalescent-based methods. Our study also highlights the importance of testing for hybridization versus incomplete lineage sorting, which aids inference of not only species limits but also evolutionary processes influencing genetic diversity.

Catchment-scale conservation units identified for the threatened Yarra pygmy perch (Nannoperca obscura) in highly modified river systems

Habitat fragmentation caused by human activities alters metapopulation dynamics and decreases biological connectivity through reduced migration and gene flow, leading to lowered levels of population genetic diversity and to local extinctions. The threatened Yarra pygmy perch, Nannoperca obscura, is a poor disperser found in small, isolated populations in wetlands and streams of southeastern Australia. Modifications to natural flow regimes in anthropogenically-impacted river systems have recently reduced the amount of habitat for this species and likely further limited its opportunity to disperse. We employed highly resolving microsatellite DNA markers to assess genetic variation, population structure and the spatial scale that dispersal takes place across the distribution of this freshwater fish and used this information to identify conservation units for management. The levels of genetic variation found for N. obscura are amongst the lowest reported for a fish species (mean heterozygosity of 0.318 and mean allelic richness of 1.92). We identified very strong population genetic structure, nil to little evidence of recent migration among demes and a minimum of 11 units for conservation management, hierarchically nested within four major genetic lineages. A combination of spatial analytical methods revealed hierarchical genetic structure corresponding with catchment boundaries and also demonstrated significant isolation by riverine distance. Our findings have implications for the national recovery plan of this species by demonstrating that N. obscura populations should be managed at a catchment level and highlighting the need to restore habitat and avoid further alteration of the natural hydrology.

Paleoclimatic modeling and phylogeography of least killifish, Heterandria formosa: insights into Pleistocene expansion-contraction dynamics and evolutionary history of North American Coastal Plain freshwater biota

BACKGROUND

Climatic and sea-level fluctuations throughout the last Pleistocene glacial cycle (~130-0 ka) profoundly influenced present-day distributions and genetic diversity of Northern Hemisphere biotas by forcing range contractions in many species during the glacial advance and allowing expansion following glacial retreat ('expansion-contraction' model). Evidence for such range dynamics and refugia in the unglaciated Gulf-Atlantic Coastal Plain stems largely from terrestrial species, and aquatic species Pleistocene responses remain relatively uninvestigated. Heterandria formosa, a wide-ranging regional endemic, presents an ideal system to test the expansion-contraction model within this biota. By integrating ecological niche modeling and phylogeography, we infer the Pleistocene history of this livebearing fish (Poeciliidae) and test for several predicted distributional and genetic effects of the last glaciation.

RESULTS

Paleoclimatic models predicted range contraction to a single southwest Florida peninsula refugium during the Last Glacial Maximum, followed by northward expansion. We inferred spatial-population subdivision into four groups that reflect genetic barriers outside this refuge. Several other features of the genetic data were consistent with predictions derived from an expansion-contraction model: limited intraspecific divergence (e.g. mean mtDNA p-distance = 0.66%); a pattern of mtDNA diversity (mean Hd = 0.934; mean π = 0.007) consistent with rapid, recent population expansion; a lack of mtDNA isolation-by-distance; and clinal variation in allozyme diversity with higher diversity at lower latitudes near the predicted refugium. Statistical tests of mismatch distributions and coalescent simulations of the gene tree lent greater support to a scenario of post-glacial expansion and diversification from a single refugium than to any other model examined (e.g. multiple-refugia scenarios).

CONCLUSIONS

Congruent results from diverse data indicate H. formosa fits the classic Pleistocene expansion-contraction model, even as the genetic data suggest additional ecological influences on population structure. While evidence for Plio-Pleistocene Gulf Coast vicariance is well described for many freshwater species presently codistributed with H. formosa, this species demography and diversification departs notably from this pattern. Species-specific expansion-contraction dynamics may therefore have figured more prominently in shaping Coastal Plain evolutionary history than previously thought. Our findings bolster growing appreciation for the complexity of phylogeographical structuring within North America's southern refugia, including responses of Coastal Plain freshwater biota to Pleistocene climatic fluctuations.