Phylogeography of Cuban Rivulus: evidence for allopatric speciation and secondary dispersal across a marine barrier

A phylogeny of the genus Limia (Teleostei: Poeciliidae) suggests a single-lake radiation nested in a Caribbean-wide allopatric speciation scenario

OBJECTIVE

The Caribbean is an important global biodiversity hotspot. Adaptive radiations there lead to many speciation events within a limited period and hence are particularly prominent biodiversity generators. A prime example are freshwater fish of the genus Limia, endemic to the Greater Antilles. Within Hispaniola, nine species have been described from a single isolated site, Lake Miragoâne, pointing towards extraordinary sympatric speciation. This study examines the evolutionary history of the Limia species in Lake Miragoâne, relative to their congeners throughout the Caribbean.

RESULTS

For 12 Limia species, we obtained almost complete sequences of the mitochondrial cytochrome b gene, a well-established marker for lower-level taxonomic relationships. We included sequences of six further Limia species from GenBank (total N  = 18 species). Our phylogenies are in concordance with other published phylogenies of Limia. There is strong support that the species found in Lake Miragoâne in Haiti are monophyletic, confirming a recent local radiation. Within Lake Miragoâne, speciation is likely extremely recent, leading to incomplete lineage sorting in the mtDNA. Future studies using multiple unlinked genetic markers are needed to disentangle the relationships within the Lake Miragoâne clade.

Biogeography of the West Indies: A complex scenario for species radiations in terrestrial and aquatic habitats

Studies of the biogeography of the West Indies are numerous but not all taxonomic groups have received the same attention. Many of the contributions to this field have historically focused on terrestrial vertebrates from a perspective closely linked to the classical theory of island biogeography. However, some recent works have questioned whether some of the assumptions of this theory are too simplistic. In this review, we compiled information about the West Indies biogeography based on an extensive and rigorous literature search. While we offer some background of the main hypotheses that explain the origin of the Caribbean biota, our main purpose here is to highlight divergent diversification patterns observed in terrestrial versus aquatic groups of the West Indian biota and also to shed light on the unbalanced number of studies covering the biogeography of these groups of organisms. We use an objective method to compile existing information in the field and produce a rigorous literature review. Our results show that most of the relevant literature in the field is related to the study of terrestrial organisms (mainly vertebrates) and only a small portion covers aquatic groups. Specifically, livebearing fishes show interesting deviations from the species-area relationship predicted by classical island biogeography theory. We found that species richness on the Greater Antilles is positively correlated with island size but also with the presence of elevations showing that not only island area but also mountainous relief may be an important factor determining the number of freshwater species in the Greater Antilles. Our findings shed light on mechanisms that may differently drive speciation in aquatic versus terrestrial environments suggesting that ecological opportunity could outweigh the importance of island size in speciation. Investigations into freshwater fishes of the West Indies offer a promising avenue for understanding origins and subsequent diversification of the Caribbean biota.

Diversification of Neotropical Freshwater Fishes

Neotropical freshwater fishes (NFFs) constitute the most diverse continental vertebrate fauna on Earth, with more than 6,200 named species compressed into an aquatic footprint <0.5% of the total regional land-surface area and representing the greatest phenotypic disparity and functional diversity of any continental ichthyofauna. Data from the fossil record and time-calibrated molecular phylogenies indicate that most higher taxa (e.g., genera, families) diversified relatively continuously through the Cenozoic, across broad geographic ranges of the South American platform. Biodiversity data for most NFF clades support a model of continental radiation rather than adaptive radiation, in which speciation occurs mainly in allopatry, and speciation and adaptation are largely decoupled. These radiations occurred under the perennial influence of river capture and sea-level oscillations, which episodically fragmented and merged portions of adjacent river networks. The future of the NFF fauna into the Anthropocene is uncertain, facing numerous threats at local, regional, and continental scales.

Genetic diversity of a widespread annual killifish from coastal Tanzania

Background

African annual killifishes (Nothobranchius spp.) are adapted to seasonally desiccating habitats (ephemeral pools), surviving dry periods as dormant eggs. Given their peculiar life history, geographic aspects of their diversity uniquely combine patterns typical for freshwater taxa (river basin structure and elevation gradient) and terrestrial animals (rivers acting as major dispersal barriers). However, our current knowledge on fine-scale inter-specific and intra-specific genetic diversity of African annual fish is limited to a single, particularly dry region of their distribution (subtropical Mozambique). Using a widespread annual killifish from coastal Tanzania and Kenya, we tested whether the same pattern of genetic divergence pertains to a wet equatorial region in the centre of Nothobranchius distribution.

Results

In populations of Nothobranchius melanospilus species group across its range, we genotyped a part of mitochondrial cytochrome oxidase subunit 1 (COI) gene (83 individuals from 22 populations) and 10 nuclear microsatellite markers (251 individuals from 16 populations). We found five lineages with a clear phylogeographic structure but frequent secondary contact. Mitochondrial lineages were largely congruent with main population genetic clusters identified on microsatellite markers. In the upper Wami basin, populations are isolated as a putative Nothobranchius prognathus, but include also a population from a periphery of the middle Ruvu basin. Other four lineages (including putative Nothobranchius kwalensis) coexisted in secondary contact zones, but possessed clear spatial pattern. Main river channels did not form apparent barriers to dispersal. The most widespread lineage had strong signal of recent population expansion.

Conclusions

We conclude that dispersal of a Nothobranchius species from a wet part of the genus distribution (tropical lowland) is not constrained by main river channels and closely related lineages frequently coexist in secondary contact zones. We also demonstrate contemporary connection between the Ruvu and Rufiji river basins. Our data do not provide genetic support for existence of recently described cryptic species from N. melanospilus complex, but cannot resolve this issue.

Evidences for genetic differentiation within the highly endemic and endangered annual fish Austrolebias nigrofasciatus (Cyprinodontiformes: Rivulidae)

Samples of Austrolebias nigrofasciatus (n = 103), an endangered species of annual fish endemic to a small area of the Patos-Mirim lagoon system encompassing the São Gonçalo Channel lowlands, were collected from eight isolated temporary ponds, four located at the known distribution range of the species and four located along the Piratini River lowlands, where morphologically different individuals were found. In the laboratory, fragments of the mitochondrial cytochrome c oxidase I (coI), cytochrome b (cytb) and nuclear rhodopsin (rho) genes were amplified, purified and sequenced for 100, 99 and 58 of these individuals, respectively. Samples were further analysed using phylogenetic and phylogeographic methods to evaluate the patterns of genetic diversity and differentiation presented within and between populations, while assessing their evolutionary history, in order to guide the application of further conservation strategies. We found that the four new populations from the Piratini River lowlands encompass a different lineage of A. nigrofasciatus that diverged from that encountered in the São Gonçalo Channel at approximately 0.165 M years before present, during a population expansion and did not yet attain reciprocal monophyly. This divergence was associated with a glacial event that was preceded by an interglacial period putatively associated with the dispersal. Moreover, significant levels of genetic differentiation and a high number of exclusive haplotypes could be encountered even in micro-geographical scales, as in the comparisons between populations located within the same major lineage, indicating each of them may encompass independent management units. Conservation actions are certainly urgent, especially in the face of signs of a recent bottleneck.

Phylogeographic evidence that the distribution of cryptic euryhaline species in the Gambusia punctata species group in Cuba was shaped by the archipelago geological history

The main drivers of diversification of freshwater fishes in Cuba are not yet well understood. For example, salt tolerance was thought as the main factor involved in the diversification of Gambusia punctata species group in this archipelago. However, evidence from a recent DNA barcoding survey suggested the presence of cryptic species and no correlation between species delimitation and level of salinity. In this study, we analyzed the cryptic diversification of G. punctata species group in Cuba, based on a comprehensive sampling of its distribution and including habitats with different salinity levels. We evaluated the patterns of molecular divergence of the samples by sequencing a set of mitochondrial DNA (mtDNA) regions and genotyping nine nuclear microsatellite loci. We also used cytochrome b gene (cytb) partial sequences and these microsatellite loci to analyze population structure inside putative species. Five mtDNA well-differentiated haplogroups were found, four of them also identified by the analysis of the microsatellite polymorphism which corresponds to two already recognized species, G. punctata, and G. rhizophorae, and three putative new species. The extent of hybrid zones between these groups is also described. In each group, populations inhabiting environments with contrasting salinity levels were identified, indicating a generalized trait not specific to G. rhizophorae. The geographic distribution of the groups suggested a strong association with major relict territories of the Cuban Archipelago that was periodically joined or split-up by changes in seawater levels and land uplifts. Salinity tolerance might have facilitated sporadic and long-distance oversea dispersal but did not prevent speciation in the Cuban archipelago.

Poecilia vivipara Bloch & Schneider, 1801 (Cyprinodontiformes, Poeciliidae), a guppy in an oceanic archipelago: from where did it come?

Poecilia vivipara, a small euryhaline guppy is reported at the Maceió River micro-basin in the Fernando de Noronha oceanic archipelago, northeast Brazil. However, the origin (human-mediated or natural dispersal) of this insular population is still controversial. The present study investigates how this population is phylogenetically related to the surrounding continental populations using the cytochrome oxidase I mitochondrial gene from eleven river basins in South America. Our phylogenetic reconstruction showed a clear geographical distribution arrangement of P. vivipara lineages. The Fernando de Noronha haplotype fell within the 'north' clade, closely related to a shared haplotype between the Paraíba do Norte and Potengi basins; the geographically closest continental drainages. Our phylogenetic reconstruction also showed highly divergent lineages, suggesting that P. vivipara may represent a species complex along its wide distribution. Regarding to the insular population, P. vivipara may have been intentionally introduced to the archipelago for the purpose of mosquito larvae control during the occupation of a U.S. military base following World War II. However, given the euryhaline capacity of P. vivipara, a potential scenario of natural (passive or active) dispersal cannot be ruled out.

Phylogeny of the subfamily Stelliferinae suggests speciation in Ophioscion Gill, 1863 (Sciaenidae: Perciformes) in the western South Atlantic

Phylogenies based on morphological and molecular data confirm the monophyly of the subfamily Stelliferinae; however, there is no consensus on the intergeneric and interspecific relationships in the group. Previous studies suggested the non-monophyly of Ophioscion and Stellifer, and possible cryptic species in Ophioscion punctatissimus. Therefore, we used mitochondrial (16S rDNA and COI) and nuclear (Rhodopsin, EGR1, and RAG1) regions to examine phylogenetic relationships among species of this subfamily. Our results confirmed the monophyly of Stelliferinae and supports the close relationship among Bardiella, Corvula and Odontoscion, which form a sister group to Stellifer and Ophioscion. Notwithstanding, all the results support the non-monophyly of Stellifer and Ophioscion and we suggest that a taxonomic revision should consider Ophioscion as a junior synonym of Stellifer. Moreover, O. punctatissimus was grouped into two clades, with the O. punctatissimus lineage I (LI) being closer to O. scierus from the eastern Pacific than to the O. punctatissimus lineage II (LII). The most recent common ancestor (TMRCA) for the O. scierus and O. punctatissimus LI and O. punctatissimus LII clade dates from 7.2 (HPD: 4.3-10.5) Ma, whereas TMRCA for the O. scierus and O. punctatissimus LI clade dates from 5.3 (HPD: 2.4-8.6) Ma, indicating that speciation processes may be related to the rise of the Isthmus of Panama. Phylogeographic analyses corroborate the hypothesis of speciation in O. punctatissimus. These results suggest that lineages of O. punctatissimus originated from distinct ancestors and, by morphological similarity, were considered the same taxon. A taxonomic revision should be performed to validate the species status of such lineages.

Review of the family Rivulidae (Cyprinodontiformes, Aplocheiloidei) and a molecular and morphological phylogeny of the annual fish genus Austrolebias Costa 1998

ABSTRACT The family Rivulidae is the fourth most diverse clade of Neotropical fishes. Together with some genera of the related African family Nothobranchiidae, many rivulids exhibit a characteristic annual life cycle, with diapausing eggs and delayed embryonic development, which allows them to survive in the challenging seasonal ponds that they inhabit. Rivulidae also includes two species known as the only the self-fertilizing vertebrates and some species with internal fertilization. The first goal of this article is to review the systematics of the family considering phylogenetic relationships and synapomorphies of subfamilial clades, thus unifying information that is dispersed throughout the literature. From this revision, it is clear that phylogenetic relationships within Rivulidae are poorly resolved, especially in one of the large clades that compose it, the subfamily Rivulinae, where conflicting hypotheses of relationships of non-annual and annual genera are evident. The second goal of this work is to present an updated phylogenetic hypothesis (based on mitochondrial, nuclear, and morphological information) for one of the most speciose genus of Rivulidae, Austrolebias. Our results confirm the monophyly of the genus and of some subgeneric clades already diagnosed, but propose new relationships among them and their species composition, particularly in the subgenus Acrolebias.

Plio-Pleistocene phylogeography of the Southeast Asian Blue Panchax killifish, Aplocheilus panchax

The complex climatic and geological history of Southeast Asia has shaped this region’s high biodiversity. In particular, sea level fluctuations associated with repeated glacial cycles during the Pleistocene both facilitated, and limited, connectivity between populations. In this study, we used data from two mitochondrial and three anonymous nuclear markers to determine whether a fresh/brackish water killifish, Aplocheilus panchax, Hamilton, 1822, could be used to further understand how climatic oscillations and associated sea level fluctuations have shaped the distribution of biota within this region, and whether such patterns show evidence of isolation within palaeodrainage basins. Our analyses revealed three major mitochondrial clades within A. panchax. The basal divergence of A. panchax mitochondrial lineages was approximately 3.5 Ma, whilst the subsequent divergence timings of these clades occurred early Pleistocene (~2.6 Ma), proceeding through the Pleistocene. Continuous phylogeographic analysis showed a clear west-east dispersal followed by rapid radiation across Southeast Asia. Individuals from Krabi, just north of the Isthmus of Kra, were more closely related to the Indian lineages, providing further evidence for a freshwater faunal disjunction at the Isthmus of Kra biogeographic barrier. Our results suggest that Sulawesi, across the Wallace Line, was colonised relatively recently (~30 ka). Nuclear DNA is less geographically structured, although Mantel tests indicated that nuclear genetic distances were correlated with geographic proximity. Overall, these results imply that recent gene flow, as opposed to historical isolation, has been the key factor determining patterns of nuclear genetic variation in A. panchax, however, some evidence of historical isolation is retained within the mitochondrial genome. Our study further validates the existence of a major biogeographic boundary at the Kra Isthmus, and also demonstrates the use of widely distributed fresh/brackishwater species in phylogeographic studies, and their ability to disperse across major marine barriers in relatively recent time periods.

First complete mitochondrial genome of the South American annual fish Austrolebias charrua (Cyprinodontiformes: Rivulidae): peculiar features among cyprinodontiforms mitogenomes

Background

Among teleosts, the South American genus Austrolebias (Cyprinodontiformes: Rivulidae) includes 42 taxa of annual fishes divided into five different species groups. It is a monophyletic genus, but morphological and molecular data do not resolve the relationship among intrageneric clades and high rates of substitution have been previously described in some mitochondrial genes. In this work, the complete mitogenome of a species of the genus was determined for the first time. We determined its structure, gene order and evolutionary peculiar features, which will allow us to evaluate the performance of mitochondrial genes in the phylogenetic resolution at different taxonomic levels.

Results

Regarding gene content and order, the circular mitogenome of A. charrua (17,271 pb) presents the typical pattern of vertebrate mitogenomes. It contains the full complement of 13 proteins-coding genes, 22 tRNA, 2 rRNA and one non-coding control region. Notably, the tRNA-Cys was only 57 bp in length and lacks the D-loop arm. In three full sibling individuals, heteroplasmatic condition was detected due to a total of 12 variable sites in seven protein-coding genes. Among cyprinodontiforms, the mitogenome of A. charrua exhibits the lowest G+C content (37 %) and GCskew, as well as the highest strand asymmetry with a net difference of T over A at 1st and 3rd codon positions. Considering the 12 coding-genes of the H strand, correspondence analyses of nucleotide composition and codon usage show that A and T at 1st and 3rd codon positions have the highest weight in the first axis, and segregate annual species from the other cyprinodontiforms analyzed. Given the annual life-style, their mitogenomes could be under different selective pressures. All 13 protein-coding genes are under strong purifying selection and we did not find any significant evidence of nucleotide sites showing episodic selection (dN >dS) at annual lineages. When fast evolving third codon positions were removed from alignments, the “supergene” tree recovers our reference species phylogeny as well as the Cytb, ND4L and ND6 genes. Therefore, third codon positions seem to be saturated in the aforementioned coding regions at intergeneric Cyprinodontiformes comparisons.

Conclusions

The complete mitogenome obtained in present work, offers relevant data for further comparative studies on molecular phylogeny and systematics of this taxonomic controversial endemic genus of annual fishes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2090-3) contains supplementary material, which is available to authorized users.