Biogeographic factors contributing to the diversification of Euphoniinae (Aves, Passeriformes, Fringillidae): a phylogenetic and ancestral areas analysis

Factors such as the Andean uplift, Isthmus of Panama, and climate changes have influenced bird diversity in the Neotropical region. Studying bird species that are widespread in Neotropical highlands and lowlands can help us understand the impact of these factors on taxa diversification. Our main objectives were to determine the biogeographic factors that contributed to the diversification of Euphoniinae and re-evaluate their phylogenetic relationships. The nextRAD and mitochondrial data were utilized to construct phylogenies. The ancestral distribution range was then estimated using a time-calibrated phylogeny, current species ranges, and neotropical regionalization. The phylogenies revealed two main Euphoniinae clades, Chlorophonia and Euphonia, similar to previous findings. Furthermore, each genus has distinctive subclades corresponding to morphology and geography. The biogeographic results suggest that the Andean uplift and the establishment of the western Amazon drove the vicariance of Chlorophonia and Euphonia during the Miocene. The Chlorophonia lineage originated in the Andes mountains and spread to Central America and the Mesoamerican highlands after the formation of the Isthmus of Panama. Meanwhile, the ancestral area of Euphonia was the Amazonas, from which it spread to trans-Andean areas during the Pliocene and Pleistocene due to the separation of the west lowlands from Amazonas due to the Northern Andean uplift. Chlorophonia and Euphonia species migrated to the Atlantic Forest during the Pleistocene through corridors from the East Andean Humid Forest and Amazonas. These two genera had Caribbean invasions with distinct geographic origins and ages. Finally, we suggested taxonomic changes in the genus Euphonia based on the study's phylogenetic, morphological, and biogeographic findings.

Evolution and environment of Caribbean coastal ecosystems

Isolation of the Caribbean Sea from the tropical Eastern Pacific by uplift of the Isthmus of Panama in the late Pliocene was associated with major, taxonomically variable, shifts in Caribbean biotic composition, and extinction, but inferred causes of these biological changes have remained elusive. We addressed this through falsifiable hypotheses about how independently determined historical changes in oceanographic conditions may have been responsible. The most striking environmental change was a sharp decline in upwelling intensity as measured from decreases in intra-annual fluctuations in temperature and consequently in planktonic productivity. We then hypothesized three general categories of biological response based upon observed differences in natural history between the oceans today. These include changes in feeding ecology, life histories, and habitats. As expected, suspension feeders and predators became rarer as upwelling declined. However, predicted increases in benthic productivity by reef corals, and benthic algae were drawn out over more than 1 Myr as seagrass and coral reef habitats proliferated; a shift that was itself driven by declining upwelling. Similar time lags occurred for predicted shifts in reproductive life history characteristics of bivalves, gastropods, and bryozoans. Examination of the spatial variability of biotic change helps to understand the time lags. Many older species characteristic of times before environmental conditions had changed tended to hang on in progressively smaller proportions of locations until they became extinct as expected from metapopulation theory and the concept of extinction debt. Faunal turnover may not occur until a million or more years after the environmental changes ultimately responsible.

Weaving Mitochondrial DNA and Y-Chromosome Variation in the Panamanian Genetic Canvas

The Isthmus of Panama was a crossroads between North and South America during the continent’s first peopling (and subsequent movements) also playing a pivotal role during European colonization and the African slave trade. Previous analyses of uniparental systems revealed significant sex biases in the genetic history of Panamanians, as testified by the high proportions of Indigenous and sub-Saharan mitochondrial DNAs (mtDNAs) and by the prevalence of Western European/northern African Y chromosomes. Those studies were conducted on the general population without considering any self-reported ethnic affiliations. Here, we compared the mtDNA and Y-chromosome lineages of a new sample collection from 431 individuals (301 males and 130 females) belonging to either the general population, mixed groups, or one of five Indigenous groups currently living in Panama. We found different proportions of paternal and maternal lineages in the Indigenous groups testifying to pre-contact demographic events and genetic inputs (some dated to Pleistocene times) that created genetic structure. Then, while the local mitochondrial gene pool was marginally involved in post-contact admixtures, the Indigenous Y chromosomes were differentially replaced, mostly by lineages of western Eurasian origin. Finally, our new estimates of the sub-Saharan contribution, on a more accurately defined general population, reduce an apparent divergence between genetic and historical data.

Base-substitution mutation rate across the nuclear genome of Alpheus snapping shrimp and the timing of isolation by the Isthmus of Panama

Background

The formation of the Isthmus of Panama and final closure of the Central American Seaway (CAS) provides an independent calibration point for examining the rate of DNA substitutions. This vicariant event has been widely used to estimate the substitution rate across mitochondrial genomes and to date evolutionary events in other taxonomic groups. Nuclear sequence data is increasingly being used to complement mitochondrial datasets for phylogenetic and evolutionary investigations; these studies would benefit from information regarding the rate and pattern of DNA substitutions derived from the nuclear genome.

Results

To estimate the genome-wide neutral mutation rate (µ), genotype-by-sequencing (GBS) datasets were generated for three transisthmian species pairs in Alpheus snapping shrimp. A range of bioinformatic filtering parameters were evaluated in order to minimize potential bias in mutation rate estimates that may result from SNP filtering. Using a Bayesian coalescent approach (G-PhoCS) applied to 44,960 GBS loci, we estimated µ to be 2.64E−9 substitutions/site/year, when calibrated with the closure of the CAS at 3 Ma. Post-divergence gene flow was detected in one species pair. Failure to account for this post-split migration inflates our substitution rate estimates, emphasizing the importance of demographic methods that can accommodate gene flow.

Conclusions

Results from our study, both parameter estimates and bioinformatic explorations, have broad-ranging implications for phylogeographic studies in other non-model taxa using reduced representation datasets. Our best estimate of µ that accounts for coalescent and demographic processes is remarkably similar to experimentally derived mutation rates in model arthropod systems. These results contradicted recent suggestions that the closure of the Isthmus was completed much earlier (around 10 Ma), as mutation rates based on an early calibration resulted in uncharacteristically low genomic mutation rates. Also, stricter filtering parameters resulted in biased datasets that generated lower mutation rate estimates and influenced demographic parameters, serving as a cautionary tale for the adherence to conservative bioinformatic strategies when generating reduced-representation datasets at the species level. To our knowledge this is the first use of transisthmian species pairs to calibrate the rate of molecular evolution from GBS data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01836-3.

A biogeographic framework of octopod species diversification: the role of the Isthmus of Panama

The uplift of the Isthmus of Panama (IP) created a land bridge between Central and South America and caused the separation of the Western Atlantic and Eastern Pacific oceans, resulting in profound changes in the environmental and oceanographic conditions. To evaluate how these changes have influenced speciation processes in octopods, fragments of two mitochondrial (Cytochrome oxidase subunit I, COI and 16S rDNA) and two nuclear (Rhodopsin and Elongation Factor-1α, EF-1α) genes were amplified from samples from the Atlantic and Pacific oceans. One biogeographical and four fossil calibration priors were used within a relaxed Bayesian phylogenetic analysis framework to estimate divergence times among cladogenic events. Reconstruction of the ancestral states in phylogenies was used to infer historical biogeography of the lineages and species dispersal routes. The results revealed three well-supported clades of transisthmian octopus sister species pair/complex (TSSP/TSSC) and two additional clades showing a low probability of species diversification, having been influenced by the IP. Divergence times estimated in the present study revealed that octopod TSSP/TSSC from the Atlantic and Pacific diverged between the Middle Miocene and Early Pliocene (mean range = 5-18 Ma). Given that oceanographic changes caused by the uplift of the IP were so strong as to affect the global climate, we suggest that octopod TSSP/TSSC diverged because of these physical and environmental barriers, even before the complete uplift of the IP 3 Ma, proposed by the Late Pliocene model. The results obtained in this phylogenetic reconstruction also indicate that the octopus species pairs in each ocean share a recent common ancestor from the Pacific Ocean.

Towards a synthesis of the Caribbean biogeography of terrestrial arthropods

BACKGROUND

The immense geologic and ecological complexity of the Caribbean has created a natural laboratory for interpreting when and how organisms disperse through time and space. However, competing hypotheses compounded with this complexity have resulted in a lack of unifying principles of biogeography for the region. Though new data concerning the timing of geologic events and dispersal events are emerging, powerful new analytical tools now allow for explicit hypothesis testing. Arthropods, with varying dispersal ability and high levels of endemism in the Caribbean, are an important, albeit understudied, biogeographic model system. Herein, we include a comprehensive analysis of every publicly available genetic dataset (at the time of writing) of terrestrial Caribbean arthropod groups using a statistically robust pipeline to explicitly test the current extent of biogeographic hypotheses for the region.

RESULTS

Our findings indicate several important biogeographic generalizations for the region: the South American continent is the predominant origin of Caribbean arthropod fauna; GAARlandia played a role for some taxa in aiding dispersal from South America to the Greater Antilles; founder event dispersal explains the majority of dispersal events by terrestrial arthropods, and distance between landmasses is important for dispersal; most dispersal events occurred via island hopping; there is evidence of 'reverse' dispersal from islands to the mainland; dispersal across the present-day Isthmus of Panama generally occurred prior to 3 mya; the Greater Antilles harbor more lineage diversity than the Lesser Antilles, and the larger Greater Antilles typically have greater lineage diversity than the smaller islands; basal Caribbean taxa are primarily distributed in the Greater Antilles, the basal-most being from Cuba, and derived taxa are mostly distributed in the Lesser Antilles; Jamaican taxa are usually endemic and monophyletic.

CONCLUSIONS

Given the diversity and deep history of terrestrial arthropods, incongruence of biogeographic patterns is expected, but focusing on both similarities and differences among divergent taxa with disparate life histories emphasizes the importance of particular qualities responsible for resulting diversification patterns. Furthermore, this study provides an analytical toolkit that can be used to guide researchers interested in answering questions pertaining to Caribbean biogeography using explicit hypothesis testing.

Army ant invasions reveal phylogeographic processes across the Isthmus of Panama

Female army ants cannot fly, making them very poor dispersers across water barriers. This dependence on terrestrial corridors motivated the investigation by Winston et al. (), published in this issue of Molecular Ecology, into the role of Panamanian isthmus formation in the diversification of Eciton army ants. Complete closure of this isthmus occurred around three million years ago (3 Ma), but it has also been hypothesized that earlier, temporary land connections facilitated additional colonization events between South and Central America over the past 13 million years or more. The phylogenomic and population genomic analyses by Winston et al. () uncovered multiple incursions of Eciton lineages into Central America between 4 and 7 Ma. Their study contributes to a growing body of evidence arguing that transitory land bridges predating 3 Ma supported substantial intercontinental biotic exchange.

Multilocus molecular systematics of the circumtropical reef-fish genus Abudefduf (Pomacentridae): history, geography and ecology of speciation

We investigated a pantropical sub-family and genus of damselfishes, the sergeant-majors (Pomacentridae: Abudefdufinae: Abudefduf), to identify the tempo and mechanisms of speciation in the lineage. We examined sequence capture data from 500 loci and 20 species, with multiple individuals sampled from across the geographic ranges of widespread species. Utilizing a maximum likelihood framework, as well as a time-calibrated Bayesian phylogeny, the following key questions are addressed: What is the historical tempo of speciation? What are the relative contributions of vicariant, peripatric and parapatric speciation to sergeant-major diversity? How is speciation related to major variation in trophic ecology? The approximately 20 species of sergeant-majors fall into three main lineages. The ancestral condition appears to be benthivory, which is predominant in two lineages comprising six species. The remaining species of sergeant-majors, of which there are at least 15, fall within a clade composed entirely of planktivores. This clade is sister to a benthivore clade that included one species, Abudefduf notatus, in transition to planktivory. Most speciation of sergeant-majors, which appeared ∼24 million years ago, occurred in the last 10 million years. Present distributional patterns indicate vicariant speciation precipitated by the closure of land barriers between both sides of the Atlantic and the Pacific, and the emergence of land between the Indian and Pacific Oceans. Within this backdrop, frequent oscillations in sea level over the last 10 million years also appear to have generated conditions suitable for both peripatric and vicariant speciation, and most speciation within the genus appears linked to these changes in sea level. Diversification within the genus has been concentrated in planktivorous seargeant-majors rather than benthivores. The root cause is unclear, but does not appear to be related to differences in dispersal potential, which is greater in the planktivorous species, due to the ability of their post-larval juveniles to raft with floating debris. This elevated speciation rate in planktivores and their propensity to form local endemics may reflect relaxation of selective pressures (e.g., on crypticity) that limit speciation in benthivorous sergeant-majors. Finally, our data allow us to clarify relationships of geminate sergeant-major species, indicating that there are subdivisions within the Atlantic for both benthivore and planktivore geminate pairs that may have misled previous studies.

Mitogenomic divergence between three pairs of putative geminate fishes from Panama

We sequenced the complete mitochondrial genomes of three pairs of congeneric peripheral fishes distributed on either side of the Isthmus of Panama in order to test their status as geminate species pairs. Our phylogenetic analysis did not support a sister relationship between Gobiomorus dormitor and G. maculatus and therefore they cannot be considered geminates. The average genetic distance of protein-coding genes between Sicydium altum and S. salvini was more than two times larger than between Atlantic and Pacific Awaous banana, suggesting different timings for their divergence across the Isthmus of Panama.

Patterns of divergence in fish species separated by the Isthmus of Panama

Background

The Pleistocene closure of Isthmus of Panama, separating the basins of the Eastern Pacific and the Caribbean Sea, created a unique natural experiment that reveals how marine faunas respond to environmental change. To explore how fishes have been affected by this tectonic event, I compare transisthmian patterns in phylogeny and morphology for geminate lineages in two families, Eleotridae (sleepers) and Apogonidae (cardinalfishes).

Results

Time-calibrated phylogenies for these families show different diversification patterns. In Eleotridae, several independent shallow instances of transisthmian divergences occur, with one or a few species on either side of the Isthmus. Among Apogonidae, a single clade of Eastern Pacific species is nested within a broad Caribbean radiation that also includes a species known from the Mediterranean. Divergence time estimates for taxa isolated by closure of the Isthmus are broadly congruent. Hypotheses dated with deeper, fossil-based legacy calibrations put the divergences in the Miocene at 7.4–15.1 Ma, while those estimated with a shallow biogeographic calibration of final Isthmus closure range from 5.1 to 9.9 Ma, in the late Miocene/early Pliocene. Eleotridae are more euryhaline than Apogonidae, but do not exhibit shallower transisthmian divergences. In both families, descendent lineages on either side of the Isthmus of Panama exhibit significant shape differences, although that distinction disappears for Apogonidae when I apply a correction for phylogenetic relationships. To evaluate the tempo and mode of continuous character evolution, I fit several single and multiple rate evolutionary models to morphometric data reconstructed on the Apogonidae phylogeny. I find that the most highly favored model, as estimated on both legacy and isthmus calibrated hypotheses, is a multiple rate Ornstein-Uhlbeck model, with a mosaic of rate shifts postulated for shape changes among fishes in the Caribbean and Eastern Pacific.

Conclusions

Although many transisthmian taxa have been compared and their phylogenies calibrated to estimate the dates associated with population sundering, few studies correlate these timing estimates with morphological change. I show that in transisthmian fish lineages, morphometric distinctions are detectable across the Isthmus, and that rates and patterns of shape change have also shifted, with variable manifestations across the body and between the Caribbean and Eastern Pacific.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0957-4) contains supplementary material, which is available to authorized users.

Timing and patterns of diversification in the Neotropical bat genus Pteronotus (Mormoopidae)

We investigate the biogeographic processes related to the origin and current patterns of distribution of the extant species of the genus Pteronotus. This clade of insectivorous bats is widely distributed in the Neotropical Region and has recently gone through a taxonomic update which increased more than twice its diversity. Using six molecular markers of 15 Pteronotus lineages ranging from Mexico to Central Brazil, we reconstruct a time-calibrated tree and infer molecular evolutionary rates for this bat genus. In addition, estimates of range evolution across phylogeny were obtained through statistical model testing among six different biogeographic models. The origin of the genus Pteronotus occurred approximately 16million years ago (Ma), with initial cladogenesis events being evenly distributed across the phylogeny. Divergence between most closely related species is recent, falling in the Pleistocene period less than 2.6Ma. Mainland lineages present congruent patterns of north versus south continent splitting while insular clades differ in their time of arrival in the Caribbean Islands. Temporal and geographic range estimates for early nodes of Pteronotus phylogeny suggest a central role of Neogene tectonic reorganizations of Central America in the group diversification process. Also, South American colonization by Pteronotus occurred early in the genus history. Founder-event speciation was an important mode of lineage splitting in Pteronotus, with two independent dispersal jumps having occurred to the Greater Antilles. Finally, Pleistocenic sea-level variation and climatic oscillations are possibly associated with divergence between sister-species and recent ages of MRCA for Pteronotus species.

Early and dynamic colonization of Central America drives speciation in Neotropical army ants

The emergence of the Isthmus of Panama is one of the most important events in recent geological history, yet its timing and role in fundamental evolutionary processes remain controversial. While the formation of the isthmus was complete around 3 million years ago (Ma), recent studies have suggested prior intercontinental biotic exchange. In particular, the possibility of early intermittent land bridges facilitating colonization constitutes a potential mechanism for speciation and colonization before full closure of the isthmus. To test this hypothesis, we employed genomic methods to study the biogeography of the army ant genus Eciton, a group of keystone arthropod predators in Neotropical rainforests. Army ant colonies are unable to disperse across water and are therefore ideally suited to study the biogeographic impact of land bridge formation. Using a reduced representation genome sequencing approach, we show that all strictly Central American lineages of Eciton diverged from their respective South American sister lineage between 4 and 7 Ma, significantly prior to the complete closure of the isthmus. Furthermore, three of the lineage pairs form extensive and coincident secondary contact zones in Costa Rica and Nicaragua, with no evidence of gene flow. Such a discrete and repeated biogeographic pattern indicates at least two waves of army ant dispersal into Central America that were separated by significant genetic divergence times. Thus, by integrating phylogenomic, population genomic and geographic evidence, we show that early colonization of Central America across the emerging Isthmus of Panamá drove parallel speciation in Eciton army ants.

Robust estimates of divergence times and selection with a poisson random field model: a case study of comparative phylogeographic data

Mutation frequencies can be modeled as a Poisson random field (PRF) to estimate speciation times and the degree of selection on newly arisen mutations. This approach provides a quantitative theory for comparing intraspecific polymorphism with interspecific divergence in the presence of selection and can be used to estimate population genetic parameters. Although the original PRF model has been extended to more general biological settings to make statistical inference about selection and divergence among model organisms, it has not been incorporated into phylogeographic studies that focus on estimating population genetic parameters for nonmodel organisms. Here, we modified a recently developed time-dependent PRF model to independently estimate genetic parameters from a nuclear and mitochondrial DNA data set of 22 sister pairs of birds that have diverged across a biogeographic barrier. We found that species that inhabit humid habitats had more recent divergence times and larger effective population sizes than those that inhabit drier habitats, and divergence time estimated from the PRF model were similar to estimates from a coalescent species-tree approach. Selection coefficients were higher in sister pairs that inhabited drier habitats than in those in humid habitats, but overall the mitochondrial DNA was under weak selection. Our study indicates that PRF models are useful for estimating various population genetic parameters and serve as a framework for incorporating estimates of selection into comparative phylogeographic studies.