Middle Miocene closure of the Central American Seaway

A global blueberry phylogeny: Evolution, diversification, and biogeography of Vaccinieae (Ericaceae)

Vaccinieae is a morphologically diverse and species-rich (∼1430 species) tribe in Ericaceae. Although the majority of diversity is tropical, Vaccinieae are best known for temperate crops (i.e., blueberries, cranberries, and lingonberries) in Vaccinium. Vaccinium itself (∼500 species) has been previously suggested as highly polyphyletic and taxonomic boundaries among many of the other genera in the tribe remain uncertain. We assessed the evolutionary history of Vaccinieae with phylogenomic analyses based on a target-enrichment dataset containing 256 low-copy nuclear loci and 210 species representing 30 of the 35 genera in the tribe and 25 of the 29 sections of Vaccinium. We conducted time-calibrated biogeographic analyses and diversification analyses to explore the area of origin and global dispersal history of the tribe. The analysis recovered a temperate North American origin for Vaccinieae approximately 30 million years ago. Tropical diversity of Vaccinieae was inferred to result from multiple, independent movements into the tropics from north-temperate ancestors. Diversification rate increases corresponded to radiation into the Andes and SE Asia. The pseudo-10-locular ovary evolved once in the tribe from the five-locular state, coinciding with the diversification of a major clade that includes most Asian Vaccinium and the group from which commercial blueberries are derived (V. sect. Cyanococcus). A reconstruction from available chromosome counts suggests that a major polyploid event predated the evolution of nearly half the diversity of Vaccinieae. The extent of polyphyly in Vaccinium documented here supports the need for a generic reclassification of the tribe.

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae)

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

Phylogeographic patterns of Cyphocharax from trans-Andean rivers and northward expansion to lower Central America (Teleostei, Curimatidae)

Phylogenetic analyses of mitochondrial and nuclear data of 31 specimens of Cyphocharax from trans-Andean rivers support the presence of one lineage of Cyphocharax aspilos in Lago Maracaibo and three cryptic lineages of Cyphocharax magdalenae: (1) Cauca-Magdalena and Ranchería, (2) León and Atrato, and (3) Chucunaque-Tuira, Santa María, and Chiriquí basins of Central America. Results suggest that the Serranía del Perijá facilitated Late Miocene cladogenetic events, whereas post-Isthmian C. magdalenae expansion was enabled by gene flow across the lower Magdalena valley and Central American lowlands. Time-calibrated phylogenetics indicate that the C. magdalenae colonized lower Central America in the Pliocene (3.7 MYA; Ma), the divergence Atrato-Magdalena occurred in Late Pliocene (3.0 Ma) and the split Ranchería-Magdalena during the Middle Pleistocene (1.3 Ma). Updated geographic distribution data support the hypothesis that the Cordillera de Talamanca functions as a barrier to northward expansion of C. magdalenae in Central America.

Complete phylogeny of Micrathena spiders suggests multiple dispersal events among Neotropical rainforests, islands and landmasses, and indicates that Andean orogeny promotes speciation

The Neotropical region is the most diverse on the planet, largely owing to its mosaic of tropical rainforests. Multiple tectonic and climatic processes have been hypothesized to contribute to generating this diversity, including Andean orogeny, the closure of the Isthmus of Panama, the GAARlandia land bridge and historical connections among currently isolated forests. Micrathena spiders are diverse and widespread in the region, and thus a complete phylogeny of this genus allows the testing of hypotheses at multiple scales. We estimated a complete, dated phylogeny using morphological data for 117 Micrathena species and molecular data of up to five genes for a subset of 79 species. Employing eventc-based approaches and biogeographic stochastic mapping while considering phylogenetic uncertainty, we estimated ancestral distributions, the timing and direction of dispersal events and diversification rates among areas. The phylogeny is generally robust, with uncertainty in the position of some of the species lacking sequences. Micrathena started diversifying around 25 Ma. Andean cloud forests show the highest in-situ speciation, while the Amazon is the major dispersal source for adjacent areas. The Dry Diagonal generated few species and is a sink of diversity. Species exchange between Central and South America involved approximately 23 dispersal events and started ~20 Ma, which is consistent with a Miocene age for the Isthmus of Panama closure. We inferred four dispersal events from Central America to the Antilles in the last 20 Myr, indicating the spiders did not reach the islands through the GAARlandia land bridge. We identified important species exchange routes among the Amazon, Andean cloud forests and Atlantic forests during the Plio-Pleistocene. Sampling all species of the genus was fundamental to the conclusions above, especially in identifying the Andean forests as the area that generated the majority of species. This highlights the importance of complete taxonomic sampling in biogeographic studies.

Biogeography of Greater Antillean freshwater fishes, with a review of competing hypotheses

In biogeography, vicariance and long-distance dispersal are often characterised as competing scenarios. However, they are related concepts, both relying on collective geological, ecological, and phylogenetic evidence. This is illustrated by freshwater fishes, which may immigrate to islands either when freshwater connections are temporarily present and later severed (vicariance), or by unusual means when ocean gaps are crossed (long-distance dispersal). Marine barriers have a strong filtering effect on freshwater fishes, limiting immigrants to those most capable of oceanic dispersal. The roles of vicariance and dispersal are debated for freshwater fishes of the Greater Antilles. We review three active hypotheses [Cretaceous vicariance, Greater Antilles-Aves Ridge (GAARlandia), long-distance dispersal] and propose long-distance dispersal to be an appropriate model due to limited support for freshwater fish use of landspans. Greater Antillean freshwater fishes have six potential source bioregions (defined from faunal similarity): Northern Gulf of México, Western Gulf of México, Maya Terrane, Chortís Block, Eastern Panamá, and Northern South America. Faunas of the Greater Antilles are composed of taxa immigrating from many of these bioregions, but there is strong compositional disharmony between island and mainland fish faunas (>90% of Antillean species are cyprinodontiforms, compared to <10% in Northern Gulf of México and Northern South America, and ≤50% elsewhere), consistent with a hypothesis of long-distance dispersal. Ancestral-area reconstruction analysis indicates there were 16 or 17 immigration events over the last 51 million years, 14 or 15 of these by cyprinodontiforms. Published divergence estimates and evidence available for each immigration event suggests they occurred at different times and by different pathways, possibly with rafts of vegetation discharged from rivers or washed to sea during storms. If so, ocean currents likely provide critical pathways for immigration when flowing from one landmass to another. On the other hand, currents create dispersal barriers when flowing perpendicularly between landmasses. In addition to high salinity tolerance, cyprinodontiforms collectively display a variety of adaptations that could enhance their ability to live with rafts (small body size, viviparity, low metabolism, amphibiousness, diapause, self-fertilisation). These adaptations likely also helped immigrants establish island populations after arrival and to persist long term thereafter. Cichlids may have used a pseudo bridge (Nicaragua Rise) to reach the Greater Antilles. Gars (Lepisosteidae) may have crossed the Straits of Florida to Cuba, a relatively short crossing that is not a barrier to gene flow for several cyprinodontiform immigrants. Indeed, widespread distributions of Quaternary migrants (Cyprinodon, Gambusia, Kryptolebias), within the Greater Antilles and among neighbouring bioregions, imply that long-distance dispersal is not necessarily inhibitory for well-adapted species, even though it appears to be virtually impossible for all other freshwater fishes.

New World dung beetle (Coleoptera: Scarabaeidae: Scarabaeinae) colonization of a recent Miocene insular territory: The case of Costa Rica

Costa Rica emerged from the seas as a new geological territory during the Miocene as an insular archipelago. It later became part of a continental area once it became a segment of Central America. Two dung beetle genera that colonized this new territory from South and North America, Canthidium and Onthophagus (Coleoptera: Scarabaeidae: Scarabaeinae), are here studied, in the first analysis of a volcanic paleo-archipelago, colonized from its emergence, and then later becoming terra firma. To assess their biodiversity distribution patterns, we analyzed the effect of biogeography, ecosystem origins, and body size on their altitudinal distribution patterns in three geographic basins of Costa Rica. Based on 32 years of collecting representing more than 158,000 specimens from 1017 localities, we undertook Generalized Linear Models of the two dung beetle genera to assess the effects of biodiversity and biogeographical distribution patterns. Canthidium and Onthophagus species ranged from 0 to 3000 m a.s.l., with an abrupt diversity decline at altitudes above 1500 m. Endemic species tended to show a higher altitudinal mean with a narrow altitudinal band distribution than non-endemic dung beetle species. Although there was a trend of decreasing species body size with the increase in altitude, such a trend depended on the distribution pattern of the species group. This possible insular-mediated endemicity mechanism has generated baffling biodiversity levels, considered the highest worldwide per unit area. Costa Rica is an expanse represented by a geographic overlap of two or more temporally disjunct areas and is not part of a natural transition zone. The effect of the insular Miocene origin of Costa Rica still pervades today, reflected by different insular syndromes shown by the dung beetle fauna. The importance of geological origins in generating biodiversity seems to have been an underrated criterion for conservation biology practices and should be considered ex officio.

Molecular phylogeny and divergence time of Harpalyce (Leguminosae, Papilionoideae), a lineage with amphitropical diversification in seasonally dry forests and savannas

Our knowledge of the systematics of the papilionoid legume tribe Brongniartieae has greatly benefitted from recent advances in molecular phylogenetics. The tribe was initially described to include species marked by a strongly bilabiate calyx and an embryo with a straight radicle, but recent research has placed taxa from the distantly related core Sophoreae and Millettieae within it. Despite these advances, the most species-rich genera within the Brongniartieae are still not well studied, and their morphological and biogeographical evolution remains poorly understood. Comprising 35 species, Harpalyce is one of these poorly studied genera. In this study, we present a comprehensive, multi-locus molecular phylogeny of the Brongniartieae, with an increased sampling of Harpalyce, to investigate morphological and biogeographical evolution within the group. Our results confirm the monophyly of Harpalyce and indicate that peltate glandular trichomes and a strongly bilabiate calyx with a carinal lip and three fused lobes are synapomorphies for the genus, which is internally divided into three distinct ecologically and geographically divergent lineages, corresponding to the previously recognized sections. Our biogeographical reconstructions demonstrate that Brongniartieae originated in South America during the Eocene, with subsequent pulses of diversification in South America, Mesoamerica, and Australia. Harpalyce also originated in South America during the Miocene at around 20 Ma, with almost synchronous later diversification in South America and Mexico/Mesoamerica beginning 10 Ma, but mostly during the Pliocene. Migration of Harpalyce from South to North America was accompanied by a biome and ecological shift from savanna to seasonally dry forest.

Marine introgressions and Andean uplift have driven diversification in neotropical Monkey tree frogs (Anura, Phyllomedusinae)

The species richness in the Neotropics has been linked to environmental heterogeneity and a complex geological history. We evaluated which biogeographic processes were associated with the diversification of Monkey tree frogs, an endemic clade from the Neotropics. We tested two competing hypotheses: the diversification of Phyllomedusinae occurred either in a "south-north" or a "north-south" direction in the Neotropics. We also hypothesized that marine introgressions and Andean uplift had a crucial role in promoting their diversification. We used 13 molecular markers in a Bayesian analysis to infer phylogenetic relationships among 57 species of Phyllomedusinae and to estimate their divergence times. We estimated ancestral ranges based on 12 biogeographic units considering the landscape modifications of the Neotropical region. We found that the Phyllomedusinae hypothetical ancestor range was probably widespread throughout South America, from Western Amazon to Southern Atlantic Forest, at 29.5 Mya. The Phyllomedusines' ancestor must have initially diverged through vicariance, generally followed by jump-dispersals and sympatric speciation. Dispersal among areas occurred mostly from Western Amazonia towards Northern Andes and the South American diagonal of dry landscapes, a divergent pattern from both "south-north" and "north-south" diversification hypotheses. Our results revealed a complex diversification process of Monkey tree frogs, occurring simultaneously with the orogeny of Northern Andes and the South American marine introgressions in the last 30 million years.

Botany and geogenomics: Constraining geological hypotheses in the neotropics with large-scale genetic data derived from plants

Decades of empirical research have revealed how the geological history of our planet shaped plant evolution by establishing well-known patterns (e.g., how mountain uplift resulted in high rates of diversification and replicate radiations in montane plant taxa). This follows a traditional approach where botanical data are interpreted in light of geological events. In this synthesis, I instead describe how by integrating natural history, phylogenetics, and population genetics, botanical research can be applied alongside geology and paleontology to inform our understanding of past geological and climatic processes. This conceptual shift aligns with the goals of the emerging field of geogenomics. In the neotropics, plant geogenomics is a powerful tool for the reciprocal exploration of two long standing questions in biology and geology: how the dynamic landscape of the region came to be and how it shaped the evolution of the richest flora. Current challenges that are specific to analytical approaches for plant geogenomics are discussed. I describe the scale at which various geological questions can be addressed from biological data and what makes some groups of plants excellent model systems for geogenomics research. Although plant geogenomics is discussed with reference to the neotropics, the recommendations given here for approaches to plant geogenomics can and should be expanded to exploring long-standing questions on how the earth evolved with the use of plant DNA.

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.

Unraveling the complex phylogeographic history of freshwater fishes in Lower Central America: A study of the electric fish Brachyhypopomus occidentalis

Lower Central America (LCA) has a complex biogeographic history shaped by the rise of the Isthmus of Panama and the global climatic oscillations of the Pleistocene. These events have been crucial in structuring biodiversity in LCA, but their consequences for the distribution and partitions of genetic diversity across the region remain to be elucidated. We combined complete mitochondrial genomes and nuclear ultraconserved elements (UCEs) to study the phylogeographic history and population genetic structure of the electric fish Brachyhypopomus occidentalis in LCA. Our results are consistent with the known phylogeographic history of B. occidentalis in LCA, but we update this history in several important ways that help illuminate the phylogeographic history of freshwater fishes in the region. We provide: i) support for three waves of colonization, two of which occurred prior to the final closure of the Panama Isthmus; ii) a more precise understanding of each colonization event, with evidence for a larger footprint of the first event, as well as genetic exchange across the continental divide in subsequent events; and iii) evidence for high levels of previously unrecognized population genetic structure across LCA. This updated model of colonization and diversification of B. occidentalis consists of three waves of dispersal and colonization, which triggered the evolution of geographic breaks in both nuclear and mitochondrial genomes across LCA. These processes are tightly linked to the dynamic uplift of the Isthmus, recent volcanic activity in the region, and the sea-level oscillations of the Pleistocene. These results improve previous phylogeographic inferences regarding the distribution and diversification of freshwater fishes in LCA, and generate testable hypotheses to guide future research exploring the factors shaping biodiversity in the region.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

Extinction debt and functional traits mediate community saturation over large spatiotemporal scales

Determining if ecological communities are saturated (have a limit to the number of species they can support) has important implications for understanding community assembly, species invasions, and climate change. However, previous studies have generally been limited to short time frames that overlook extinction debt and have not explicitly considered how functional trait diversity may mediate patterns of community saturation. Here, we combine data from biodiversity surveys with functional and phylogenetic data to explore if the colonisation events after the Great American Biotic Interchange (closure of the Panamanian Isthmus) resulted in increases in species richness of communities of the snake family Dipsadidae. We determined the number and the direction of dispersal events between Central and South America by estimating ancestral areas based on a Bayesian time-calibrated phylogenetic analysis. We then evaluated whether variation in community saturation was mediated by the functional similarity of six traits for the resident and colonizing snakes and/or local environmental conditions. We found that colonised communities did not support more species than those that were not colonised. Moreover, we did not find an association between the functional diversity across sites and whether they were colonised by members from the lineages dispersing across the Isthmus or not. Instead, variation in species richness was predicted best by covariates such as time since colonisation and local environment. Taken together, our results suggest that snake communities of the Dipsadidae across the neotropics are saturated. Moreover, our research highlights two important factors to consider in studies of community saturation: extinction debt and the functional differences and similarities in species' ecological roles.

Phylogenetics and Population Genetics of the Petrolisthes lamarckii-P. haswelli Complex in China: Old Lineage and New Species

Petrolisthes lamarckii (Leach, 1821) and P. haswelli Miers, 1884 are a pair of sister species of porcelain crabs, both of which are common in the intertidal zone of southern China, typically found under rocks and in the crevices of coral reefs. However, the distribution, genetic relationship and diversity of the two species in China have not been rigorously studied. Meanwhile, P. lamarckii is considered as a complex of cryptic species due to their diverse morphological features. In this study, we identified 127 specimens of the P. lamarckii-P. haswelli complex (LH complex) and recognised a new species through morphological and molecular analysis. Furthermore, we constructed a time-calibrated phylogeny of the LH complex using three mitochondrial and two nuclear genes from all three species, finding that the divergence of the LH complex can be traced back to the Miocene epoch, and that the genetic diversity increased during the Mid-Pleistocene transition period. Glacial refugia formed during the Pleistocene climatic oscillations has been regarded as one of the contributing factors to the diversification of marine organisms in the north-western Pacific. Petrolisthes haswelli demonstrates a wide distribution along the southern coast of China, while other lineages display more restricted distributions. The research on the demographic history and gene flow of P. haswelli revealed that the Chinese coastal populations experienced an expansion event approximately 12.5 thousand years ago (Kya) and the asymmetrical gene flows were observed between the two sides of the Taiwan Strait and Qiongzhou Strait, respectively, which is likely influenced by the restriction of ocean currents.

Historical Assembly of Andean Tree Communities

Patterns of species diversity have been associated with changes in climate across latitude and elevation. However, the ecological and evolutionary mechanisms underlying these relationships are still actively debated. Here, we present a complementary view of the well-known tropical niche conservatism (TNC) hypothesis, termed the multiple zones of origin (MZO) hypothesis, to explore mechanisms underlying latitudinal and elevational gradients of phylogenetic diversity in tree communities. The TNC hypothesis posits that most lineages originate in warmer, wetter, and less seasonal environments in the tropics and rarely colonize colder, drier, and more seasonal environments outside of the tropical lowlands, leading to higher phylogenetic diversity at lower latitudes and elevations. In contrast, the MZO hypothesis posits that lineages also originate in temperate environments and readily colonize similar environments in the tropical highlands, leading to lower phylogenetic diversity at lower latitudes and elevations. We tested these phylogenetic predictions using a combination of computer simulations and empirical analyses of tree communities in 245 forest plots located in six countries across the tropical and subtropical Andes. We estimated the phylogenetic diversity for each plot and regressed it against elevation and latitude. Our simulated and empirical results provide strong support for the MZO hypothesis. Phylogenetic diversity among co-occurring tree species increased with both latitude and elevation, suggesting an important influence on the historical dispersal of lineages with temperate origins into the tropical highlands. The mixing of different floras was likely favored by the formation of climatically suitable corridors for plant migration due to the Andean uplift. Accounting for the evolutionary history of plant communities helps to advance our knowledge of the drivers of tree community assembly along complex climatic gradients, and thus their likely responses to modern anthropogenic climate change.

An elevational shift facilitated the Mesoamerican diversification of Azure-hooded Jays (Cyanolyca cucullata) during the Great American Biotic Interchange

The Great American Biotic Interchange (GABI) was a key biogeographic event in the history of the Americas. The rising of the Panamanian land bridge ended the isolation of South America and ushered in a period of dispersal, mass extinction, and new community assemblages, which sparked competition, adaptation, and speciation. Diversification across many bird groups, and the elevational zonation of others, ties back to events triggered by the GABI. But the exact timing of these events is still being revealed, with recent studies suggesting a much earlier time window for faunal exchange, perhaps as early as 20 million years ago (Mya). Using a time-calibrated phylogenetic tree, we show that the jay genus Cyanolyca is emblematic of bird dispersal trends, with an early, pre-land bridge dispersal from Mesoamerica to South America 6.3-7.3 Mya, followed by a back-colonization of C. cucullata to Mesoamerica 2.3-4.8 Mya, likely after the land bridge was complete. As Cyanolyca species came into contact in Mesoamerica, they avoided competition due to a prior shift to lower elevation in the ancestor of C. cucullata. This shift allowed C. cucullata to integrate itself into the Mesoamerican highland avifauna, which our time-calibrated phylogeny suggests was already populated by higher-elevation, congeneric dwarf-jays (C. argentigula, C. pumilo, C. mirabilis, and C. nanus). The outcome of these events and fortuitous elevational zonation was that C. cucullata could continue colonizing new highland areas farther north during the Pleistocene. Resultingly, four C. cucullata lineages became isolated in allopatric, highland regions from Panama to Mexico, diverging in genetics, morphology, plumage, and vocalizations. At least two of these lineages are best described as species (C. mitrata and C. cucullata). Continued study will further document the influence of the GABI and help clarify how dispersal and vicariance shaped modern-day species assemblages in the Americas.

Bridging two oceans: small toothed cetaceans (Odontoceti) from the Late Miocene Chagres Formation, eastern Caribbean (Colon, Panama)

Fossil cetaceans are often found in Miocene marine outcrops across the globe. However, because this record is not homogeneous, the dissimilar increase in occurrences, along with the sampling bias has created regions with extensive records and others with great scarcity. Among these, the Caribbean has remained enigmatic due to the lack of well-preserved cetacean fossils. Here, we report new Caribbean fossil cetaceans from the Upper Miocene Chagres Formation exposed along Piña beach, Eastern Panama, including a scaphokogiine kogiid, an Acrophyseter-like physeteroid and the phocoenid Piscolithax. Along with previous records of the iniid Isthminia panamensis and the kogiine Nanokogia isthmia, the Chagres cetacean fauna shows some similarities with other Late Miocene cetacean communities such as the Californias in the North Pacific, although their closest affinities lie with the eastern South Pacific Pisco Formation, Peru. Such findings indicate that though deep and intermediate Caribbean-Pacific water interchange was reduced by the Middle Miocene due to the shallowing of the Central American Seaway, shallow waters marine connection that persisted until the Pliocene might have facilitated the dispersal of coastal species across both sides of the Isthmus.

Diversification trajectories and paleobiogeography of Neogene chondrichthyans from Europe

Despite the rich fossil record of Neogene chondrichthyans (chimaeras, sharks, rays, and skates) from Europe, little is known about the macroevolutionary processes that generated their current diversity and geographical distribution. We compiled 4368 Neogene occurrences comprising 102 genera, 41 families, and 12 orders from four European regions (Atlantic, Mediterranean, North Sea, and Paratethys) and evaluated their diversification trajectories and paleobiogeographic patterns. In all regions analyzed, we found that generic richness increased during the early Miocene, then decreased sharply during the middle Miocene in the Paratethys, and moderately during the late Miocene and Pliocene in the Mediterranean and North Seas. Origination rates display the most significant pulses in the early Miocene in all regions. Extinction rate pulses varied across regions, with the Paratethys displaying the most significant pulses during the late Miocene and the Mediterranean and North Seas during the late Miocene and early Pliocene. Overall, up to 27% and 56% of the European Neogene genera are now globally and regionally extinct, respectively. The observed pulses of origination and extinction in the different regions coincide with warming and cooling events that occurred during the Neogene globally and regionally. Our study reveals complex diversity dynamics of Neogene chondrichthyans from Europe and their distinct biogeographic composition despite the multiple marine passages that connected the different marine regions during this time.

An updated phylogeny, biogeography, and PhyloCode-based classification of Cornaceae based on three sets of genomic data

PREMISE

A major goal of systematic biology is to uncover the evolutionary history of organisms and translate that knowledge into stable classification systems. Here, we integrate three sets of genome-wide data to resolve phylogenetic relationships in Cornaceae (containing only Cornus s.l.), reconstruct the biogeographic history of the clade, and provide a revised classification using the PhyloCode to stabilize names for this taxonomically controversial group.

METHODS

We conducted phylogenetic analyses using 312 single-copy nuclear genes and 70 plastid genes from Angiosperms353 Hyb-Seq, plus numerous loci from RAD-Seq. We integrated fossils using morphological data and produced a dated phylogeny for biogeographical analysis.

RESULTS

A well-resolved, strongly supported, comprehensive phylogeny was obtained. Biogeographic analyses support an origin and rapid diversification of Cornus into four morphologically distinct major clades in the Northern Hemisphere (with an eastern Asian ancestor) during the late Cretaceous. Dispersal into Africa from eastern Asia likely occurred along the Tethys Seaway during the Paleogene, whereas dispersal into South America likely occurred during the Neogene. Diversification within the northern hemisphere likely involved repeated independent colonization of new areas during the Paleogene and Neogene along the Bering Land Bridge, the North Atlantic Land Bridge, and the Tethys Seaway. Thirteen strongly supported clades were named following rules of the PhyloCode.

CONCLUSIONS

Our study provides an example of integrating genomic and morphological data to produce a robust, explicit species phylogeny that includes fossil taxa, which we translate into an updated classification scheme using the PhyloCode to stabilize names.

Environmental DNA (eDNA) reveals potential for interoceanic fish invasions across the Panama Canal

Interoceanic canals can facilitate biological invasions as they connect the world's oceans and remove dispersal barriers between bioregions. As a consequence, multiple opportunities for biotic exchange arise and the resulting establishment of migrant species often causes adverse ecological and economic impacts. The Panama Canal is a key region for biotic exchange as it connects the Pacific and Atlantic Oceans in Central America. In this study, we used two complementary methods (environmental DNA (eDNA) metabarcoding and gillnetting) to survey fish communities in this unique waterway. Using COI (cytochrome oxidase subunit I) metabarcoding, we detected a total of 142 fish species, including evidence for the presence of sixteen Atlantic and eight Pacific marine fish in different freshwater sections of the Canal. Of these, nine are potentially new records. Molecular data did not capture all species caught with gillnets, but generally provided a more complete image of the known fish fauna as more small-bodied fish species were detected. Diversity indices based on eDNA surveys revealed significant differences across different sections of the Canal reflecting in part the prevailing environmental conditions. The observed increase in the presence of marine fish species in the Canal indicates a growing potential for interoceanic fish invasions. The potential ecological and evolutionary consequences of this increase in marine fishes are not only restricted to the fish fauna in the Canal as they could also impact adjacent ecosystems in the Pacific and Atlantic Oceans.

Cryptic species diversity in a widespread neotropical tree genus: The case of Cedrela odorata

PREMISE

Reconciling the use of taxonomy to partition morphological variation and describe genetic divergence within and among closely related species is a persistent challenge in phylogenetics. We reconstructed phylogenetic relationships among Cedrela odorata (Meliaceae) and five closely allied species to test the genetic basis for the current model of species delimitation in this economically valuable and threatened genus.

METHODS

We prepared a nuclear species tree with the program SNPhylo and 16,000 single-nucleotide polymorphisms from 168 Cedrela specimens. Based on clades present and ancestral patterns ADMIXTURE, we designed nine species delimitation models and compared each model to current taxonomy with Bayes factor delimitation. Timing of major lineage divergences was estimated with the program SNAPP.

RESULTS

The resulting analysis revealed that modern C. odorata evolved from two genetically distinct ancestral sources. All species delimitation models tested better fit the data than the model representing current taxonomic delimitation. Models with the greatest marginal likelihoods separated Mesoamerican C. odorata and South American C. odorata into two species and lumped C. angustifolia and C. montana as a single species. We estimated that Cedrela diversified in South America within the last 19 million years following one or more dispersal events from Mesoamerican lineages.

CONCLUSIONS

Our analyses show that the present taxonomic understanding within the genus obscures divergent lineages in C. odorata due in part to morphological differentiation and taxonomic distinctions that are not predictably associated with genetic divergence. A more accurate application of taxonomy to C. odorata and related species may aid in its conservation, management, and restoration efforts.

Reinhold Fritzgaertner y los reportes de fósiles de Honduras a finales del siglo XIX

Durante el siglo XIX el Gobierno de Honduras propició la llegada de académicos extranjeros, siendo el geólogo y paleontólogo alemán Reinhold Fritzgaertner uno de los casos más notorios por su influencia en el impulso de la actividad minera en Honduras, sin embargo, se ha descrito muy poco su rol en la historia de la paleontología de América Central. A través de una extensa revisión documental en bibliotecas virtuales de Europa y América, en el repositorio virtual Tzibalnaah-UNAH y en el Archivo Nacional de Honduras, se buscó determinar los vínculos académicos de Reinhold Fritzgaertner y su aporte en la investigación de la paleontología de Honduras en el siglo XIX. Fritzgaertner llega a Honduras después de completar sus estudios de doctorado bajo la asesoría del renombrado geólogo y paleontólogo alemán Friedrich August von Quenstedt y de fungir como asistente del renombrado geólogo y paleontólogo norteamericano James Hall, manteniendo vínculos con importantes instituciones académicas de la época. Su presencia en Honduras trascendió la minería y desarrolló un papel notorio en la historia de la ciencia de este país, en una época en que la evidencia paleontológica de Honduras adquirió importancia en el mundo académico. Fritzgaertner puede ser considerado el académico que tuvo el mayor conocimiento de la paleontología de Honduras en el siglo XIX, y mediante su participación en las diferentes exposiciones mundiales, los reportes en medios escritos de divulgación y sus vínculos académicos internacionales, propició el reporte de fósiles de Honduras a finales del siglo XIX. Palabras clave: Paleontología, América Central, Historia de Honduras, Exposiciones mundiales.

Molecular phylogeny and historical biogeography of marine palaemonid shrimps (Palaemonidae: Palaemonella-Cuapetes group)

Palaemonidae is the most speciose shrimp family within the infraorder Caridea, composed predominately of freshwater species and marine symbiotic species. The subject of this study is a clade of mainly free-living marine taxa representing a basally separated lineage from most of the symbiotic marine palaemonid genera. Phylogenetic and biogeographic relationships were explored by analysing sequence data from two mitochondrial and four nuclear markers. Maximum likelihood and Bayesian analyses, based on sequences from 52 species of 11 genera, provided similar tree topologies revealing the genera Palaemonella, Cuapetes and Eupontonia as non-monophyletic groups. Divergence time and S-DIVA analyses reveals that the focal clade originated during the Late Cretaceous in the Paleotethys region respective to the present Indo-West Pacific area, a minor part of which spread out to the eastern Pacific during the Paleocene, followed by further migration into the Atlantic (before the closure of the Panama Isthmus). The ancestral state reconstruction of host associations revealed eight independent symbiotic lineages originating from free-living ancestors, entering primary symbioses. The first associations with Cnidaria are estimated to have evolved in the Eocene. This study points to the need of taxonomic revisions of the non-monophyletic genera concerned.

Evolution of andrenine bees reveals a long and complex history of faunal interchanges through the Americas during the Mesozoic and Cenozoic

Bees are presumed to have arisen in the early to mid-Cretaceous coincident with the fragmentation of the southern continents and concurrently with the early diversification of the flowering plants. Here, we apply DNA sequences from multiple genes to recover a dated phylogeny and historical biogeographic of andrenine bees, a large group of 3000 species mainly distributed in arid areas of North America, South America, and the Palearctic region. Our results corroborate the monophyly of Andreninae and points toward a South America origin for the group during the Late Cretaceous. Overall, we provide strong evidence of amphitropical distributional pattern currently observed in the American continent as result of faunal interchange in at least three historical periods, much prior to the Panama Isthmus closure. The Palearctic diversity is shown to have arisen from North America during the Eocene and Miocene, and the Afrotropical lineages likely originated from the Palearctic region in the Miocene when the Sahara Desert was mostly vegetated. The incursions from South to North America and then onto the Old World are chronological congruent with periods when open-vegetation habitats were available for trans-continental dispersal and at the times when aridification and temperature decline offered favorable circumstances for bee diversification.

Mitochondrial phylogeny of the brittle star genus Ophioderma

We reconstructed the mitochondrial phylogeny of the species of the brittle star genus Ophioderma, using sequences of the Cytochrome Oxidase I gene (COI) to address four questions: (i) Are the species of Ophioderma described on morphological evidence reflected in mitochondrial genealogy? (ii) Which species separated from which? (iii) When did speciation events occur? (iv) What is the rate of COI evolution in ophiuroids? We found that most of the 22 described species we sampled coincide with monophyletic clusters of COI sequences, but there are exceptions, most notably in the eastern Pacific, in which three undescribed species were indicated. The COI phylogeny lacks resolution in the deeper nodes, but it does show that there are four species pairs, the members of which are found on either side of the central American Isthmus. Two pairs with a genetic distance of ~ 4% between Atlantic and Pacific members were probably split during the final stages of Isthmus completion roughly 3 million years ago. The rate of divergence provided by these pairs allowed the calibration of a relaxed molecular clock. Estimated dates of divergence indicate that the lineages leading to extant species coalesce at times much older than congeneric species in other classes of echinoderms, suggesting that low extinction rates may be one of the reasons that ophiuroids are species-rich. The mean rate of COI substitution in Ophioderma is three times slower than that of echinoids. Conclusions of previous mitochondrial DNA studies of ophiuroids that relied on echinoid calibrations to determine divergence times need to be revised.

Biogeography and Diversification of Bumblebees (Hymenoptera: Apidae), with Emphasis on Neotropical Species

A detailed phylogeny of bumblebees is urgently needed to understand speciation and biogeographic diversification in the Neotropical region. We sequenced autosomal and mtDNA loci from nine Brazilian bumblebee species and compiled it with the data already available to obtain highly resolved phylogenetic trees with fossil-calibrated dates. The ancestral Bombus lineage was estimated to diversify between 47.08 and 34.27 million years ago (Ma) in the Holarctic region, but largely restricted to the eastern Old World. The Neotropical region was initially colonized in the Late Miocene, where bumblebee diversification was shown to be consistent with geologic and climatic events of the Late Cenozoic. Neotropical bumblebees likely originated from Nearctic lineages, which dispersed towards South America after 29 Ma.

Comparative phylogeography among eight Neotropical wild cat species: no single evolutionary pattern

The felid species of South America are thought to have arrived on the continent during the Great American Biotic Interchange (GABI) in the Pleistocene. However, molecular and palaeontological data do not agree on how this event affected speciation in felids. Here, we determine both the number of colonization events and the period when felines first migrated from North America to South America. In addition, we evaluate whether similar evolutionary events could have affected the eight Neotropical cat species in their levels of genetic diversity, spatial genetic structure and demographic changes. We analysed four concatenated mitochondrial genes of the jaguar, ocelot, margay, tigrina, pampas cat, Andean cat, puma and jaguarundi. The samples were representative of a wide distribution of these species in Central and South America. Our analysis suggests either three or four colonization events from North America to South America over the past 3 Myr, followed by subsequent speciation events and the attainment of high or very high genetic diversity levels for seven of the species. No unique evolutionary process was detected for any of the current Neotropical cat species.

The multicausal twilight of South American native mammalian predators (Metatheria, Sparassodonta)

Sparassodonts were the apex mammalian predators of South America throughout most of the Cenozoic, diversifying into a wide array of niches including fox-like and even saber-toothed forms. Their extinction is still controversial, with different authors suggesting competition with other predators (placental carnivorans, terror birds, and carnivorous opossums), extinction of prey, and climate change as causal explanations. Here, we analyse these hypotheses using a novel approach implicating Bayesian analyses. We find that speciation and extinction rates of sparassodonts can be correlated with (i) intrinsic biotic factors such as changes in body mass and diversity of sparassodonts, (ii) extrinsic biotic factors such as potential prey diversity, and iii) extrinsic abiotic factors like the atmospheric CO₂, sea level, temperature, and uplift of the Andes. Thus, sparassodonts are a good example of a multilevel mixed model of evolution, where various factors drove the evolutionary history of this clade in a pluralistic way. There is no evidence for competition between Sparassodonta and others predators, and the effect of competition in the face of extinctions of fossil species should be tested and not assumed. Furthermore, we propose a novel approach for evaluating the fossil record when performing macroevolutionary analyses.

Historical biogeographical analysis of the Udoteaceae (Bryopsidales, Chlorophyta) elucidates origins of high species diversity in the Central Indo-Pacific, Western Indian Ocean and Greater Caribbean regions

There is a growing interest in elucidating the biogeographical processes underlying biodiversity patterns of seaweeds, with recent studies largely focusing on red and brown macroalgae. This study focuses on the siphonous green algal family Udoteaceae, which is diverse and globally distributed in tropical to warm-temperate seas, and includes species that form important components of tropical reefs. We explored the historical processes that have shaped current biodiversity patterns in the family by analyzing a comprehensive dataset of 568 specimens sampled across its geographical range, and including 45 species, corresponding to 59% of the known diversity. Historical biogeographical analysis was based on a three-locus time-calibrated phylogeny, and probabilistic modeling of geographical range evolution. Many species were found to have restricted ranges, indicative of low dispersal capacity. Our analysis points toward a Western Tethys origin and early diversification of the Udoteaceae in the Triassic period. Three centers of diversity were identified, which are, in order of highest species richness, the Central Indo-Pacific, the Western Indian Ocean, and the Greater Caribbean. Different drivers have likely played a role in shaping these diversity centres. Species richness in the Central Indo-Pacific likely resulted from speciation within the region, as well as recolonization from neighbouring regions, and overlap of some wider ranged species, corroborating the "biodiversity feedback" model. Species richness in the Western Indian Ocean can be explained by ancient and more recent diversification within the region, and dispersal from the Central Indo-Pacific. The Greater Caribbean region was colonized more recently, followed by diversification within the region.

Divergence times of the Rhoadsia clade (Characiformes: Characidae)

Abstract The family Characidae is the most diverse group of fishes in the Neotropics with challenging systematics. The three genera Carlana, Parastremma, and Rhoadsia, formerly considered the subfamily Rhoadsiinae, are now included in the subfamily Stethaprioninae. Previous phylogenetic analyses did not include all genera of Rhoadsiinae, specifically Parastremma. Here, we estimated the phylogenetic relationships and divergence times of the genera of Rhoadsiinae (the Rhoadsia clade) relative to the most representative genera of the Characidae. We used six molecular markers from the mitochondrial and nuclear genome to estimate the phylogeny and divergence times. We confirmed the monophyly of the Rhoadsia clade. Furthermore, we estimated that the Central American genus Carlana and the western Colombian genus Parastremma diverged approximately 13 Mya (95% HPD 8.36–18.11), consistent with the early-closure estimates of the Isthmus of Panama (~15 Mya). The genus Rhoadsia, endemic to Western Ecuador and Northern Peru, was estimated to originate at around 20 Mya (95% HPD 14.35–25.43), consistent with the Andean uplift (~20 Mya).

Testing spatial and environmental factors to explain body shape variation in the widespread Central American Blackbelt cichlid Vieja maculicauda (Teleostei: Cichlidae)

Abstract Body-shape variability may respond to the interaction between history and environment, particularly in species whose range includes known biogeographical barriers. Central America has traditionally been regarded to have functioned as a complete land-bridge connecting two continents since the Plio-Pleistocene and as an incomplete one since much earlier. This history has helped shape species diversification and distributions. Vieja maculicauda is a widely distributed cichlid, found throughout most of the Atlantic slope of Central America, across three of the four geological blocks that make up this region and whose borders are known dispersal barriers. We asked whether there is an effect of geologic blocks on body shape of V. maculicauda as determined by geometric morphometric assessment of 151 individuals from across its range. We asked how much variance in body shape is explained by environmental variables. We used a CVA and a MANCOVA to determine whether a relationship between body shape and geologic blocks exists. For testing body shape correlation against geographic and environmental variables, we used a Mantel test. Our results suggest that body shape in V. maculicauda can be segregated by geologic blocks, while the correlation results showed no strong correlation between our sets of variables.

Phylogeny indicates polyphyly in Cnodocentron (Trichoptera: Xiphocentronidae): biogeography and revision of New World species (Caenocentron)

Cnodocentron has a disjunct Laurasian distribution, with species in South-East Asia and the New World. It is divided into two subgenera: Cnodocentron and Caenocentron. Here, we infer the evolutionary history of the genus through phylogenetic and biogeographic data combining COI and 46 morphological characters. Phylogenetic relationships and divergence-time estimation were simultaneously inferred through Bayesian inference. The dating analysis was performed through relaxed morphological and molecular clocks. The historical biogeography was investigated using the dispersal–extinction–cladogenesis model. Our results indicate the polyphyly of the genus Cnodocentron, with each subgenus being more related to other genera than to each other. Therefore, the subgenus Caenocentron is elevated to genus status here. The biogeographical analysis showed that the Oriental Cnodocentron diverged in the Indian subcontinent during the Middle Eocene, while the New World Caenocentron stat. nov. originated in the Chortis Block in the Late Eocene. The dispersal of Caenocentron to South America occurred only after the Late Miocene, around 10 Mya. Additionally, we provide a revision of Caenocentron, with an identification key and description of the male and female of two new species from Costa Rica: Caenocentron carlosdelarosai sp. nov. and Caenocentron rafamoralesi sp. nov..

Geodispersal of the Typical Neotropical cenocron from South America to the Mexican Transition Zone: a cladistic biogeographical test

Analysis of the biotic assembly of the Mexican Transition Zone (MTZ) is challenging because of the convergence of taxa with different origins and biogeographical histories. The typical Neotropical cenocron corresponds to genera widespread in South America that might have dispersed to the MTZ in the Pliocene–Holocene. It has been exemplified with distributional patterns of several plant and animal taxa; however, their historical congruence has not been tested. We examined the historical relationships among the areas of endemism where Typical Neotropical taxa are found through a cladistic biogeographical analysis to determine whether they have dispersed simultaneously, showing a geodispersal pattern. We searched for Neotropical clades that dispersed to the MTZ from the Pliocene to the present. The distribution of nine animal clades in seven areas of endemism was analysed by cladistic biogeography. Redundant areas and widespread taxa were treated with paralogy-free and transparent methods. Patterns of area relationships were searched using a parsimony analysis of paralogy-free subtrees. The parsimony analysis found a single resolved tree showing a general pattern of area relationships: (Chacoan ((South Brazilian (Boreal Brazilian + Lesser Antilles)), (South American Pacific (Mesoamerican-Central America + Western Mexico)))). Distributional and phylogenetic information on the groups analysed contrast with a previous cladistic biogeographical analysis that has shown different area relationships, thus supporting the general hypothesis of a geodispersal event of the Typical Neotropical cenocron in the MTZ. Additionally, the analysis provided evidence of vicariant events related to the lifting of the Northern Andes and climatic changes during the Pleistocene.

The Yorktown Formation: Improved Stratigraphy, Chronology, and Paleoclimate Interpretations from the U.S. Mid-Atlantic Coastal Plain

The Yorktown Formation records paleoclimate conditions along the mid-Atlantic Coastal Plain during the mid-Piacenzian Warm Period (3.264 to 3.025 Ma), a climate interval of the Pliocene in some ways analogous to near future climate projections. To gain insight into potential near future changes, we investigated Yorktown Formation outcrops and cores in southeastern Virginia, refining the stratigraphic framework. We analyzed 485 samples for alkenone-based sea surface temperature (SST) and productivity estimates from the Holland and Dory cores, an outcrop at Morgarts Beach, Virginia, and the lectostratotype of the Yorktown Formation at Rushmere, Virginia, and analyzed planktonic foraminferal assemblage data from the type section. Using the structure of the SST record, we improved the chronology of the Yorktown Formation by establishing the maximum age ranges of the Rushmere (3.3–3.2 Ma) and Morgarts Beach (3.2–3.15 Ma) Members. SST values for these members average ~26 °C, corroborating existing sclerochronological data. Increasing planktonic foraminifer abundance, productivity, and species diversity parallel increasing SST over the MIS M2/M1 transition. These records constitute the greatest temporal concentration of paleoecological estimates within the Yorktown Formation, aiding our understanding of western North Atlantic temperature patterns, seasonality and ocean circulation during this interval. We provide a chronologic framework for future studies analyzing ecological responses to profound climate change.

The evolutionary history of the Caribbean magnolias (Magnoliaceae): Testing species delimitations and biogeographical hypotheses using molecular data

The Caribbean islands provide an ideal setting for studying biodiversity, given their complex geological and environmental history, and their historical and current geographical proximity to the American mainland. Magnolia, a flagship tree genus that has 15 endemic and threatened taxa (12 species and 3 subspecies) on the Caribbean islands, offers an excellent case study to empirically test Caribbean biogeographical hypotheses. We constructed phylogenetic hypotheses to: (1) reveal their evolutionary history, (2) test the current largely morphology-based classification and assess species limits, and (3) investigate major biogeographic hypotheses proposed for the region. Nuclear and chloroplast DNA sequence data of all 15 Caribbean Magnolia taxa are included, supplemented by a selection of American mainland species, and species representing most major clades of the Magnoliaceae family. We constructed phylogenetic hypotheses in a time-calibrated Bayesian framework, supplemented with haplotype network analyses and ancestral range estimations. Genetic synapomorphies in the studied markers confirm the species limits of 14 out of 15 morphologically recognizable Caribbean Magnolia taxa. There is evidence for four colonization events of Magnolia into the Caribbean from the American mainland, which most likely occurred by overwater dispersal, given age estimates of maximum 16 mya for their presence on the Caribbean islands.

Around the world in 40 million years: Phylogeny and biogeography of Tecomeae (Bignoniaceae)

Intercontinental disjunct distributions can arise from vicariance, long distance dispersal, or both. Tecomeae (Bignoniaceae) are a nearly cosmopolitan clade of flowering plants providing us with an excellent opportunity to investigate global distribution patterns. While the tribe contains only about 57 species, it has achieved a distribution that is not only pantropical, but also extends into the temperate zones in both the Northern and Southern hemispheres. This distribution is similar to the distribution of its sister group, a clade of about 750 spp. that includes most remaining taxa in Bignoniaceae. To infer temporal and spatial patterns of dispersal, we generated a phylogeny of Tecomeae by gathering sequence data from chloroplast and nuclear markers for 41 taxa. Fossil calibrations were used to determine divergence times, and ancestral states were reconstructed to infer its biogeographic history. We found support for a South American origin and a crown age of the tribe estimated at ca. 40 Ma. Two dispersal events seem to have happened during the Eocene-Oligocene, one from South America to the Old World, and another from South America to North America. Furthermore, two other dispersal events seem to have taken place during the Miocene, one from North America to Asia, and another from Australia to South America. We suggest that intercontinental dispersal via land bridges and island hopping, as well as sweepstakes of long distance dispersal from the Eocene to the present explain the global distribution of Tecomeae.

Global biogeography of living brachiopods: Bioregionalization patterns and possible controls

The global distribution patterns of 14918 geo-referenced occurrences from 394 living brachiopod species were mapped in 5° grid cells, which enabled the visualization and delineation of distinct bioregions and biodiversity hotspots. Further investigation using cluster and network analyses allowed us to propose the first systematically and quantitatively recognized global bioregionalization framework for living brachiopods, consisting of five bioregions and thirteen bioprovinces. No single environmental or ecological variable is accountable for the newly proposed global bioregionalization patterns of living brachiopods. Instead, the combined effects of large-scale ocean gyres, climatic zonation as well as some geohistorical factors (e.g., formation of land bridges and geological recent closure of ancient seaways) are considered as the main drivers at the global scale. At the regional scale, however, the faunal composition, diversity and biogeographical differentiation appear to be mainly controlled by seawater temperature variation, regional ocean currents and coastal upwelling systems.

When adaptive radiations collide: Different evolutionary trajectories between and within island and mainland lizard clades

Oceanic islands are known as test tubes of evolution. Isolated and colonized by relatively few species, islands are home to many of nature's most renowned radiations from the finches of the Galápagos to the silverswords of the Hawaiian Islands. Despite the evolutionary exuberance of insular life, island occupation has long been thought to be irreversible. In particular, the presumed much tougher competitive and predatory milieu in continental settings prevents colonization, much less evolutionary diversification, from islands back to mainlands. To test these predictions, we examined the ecological and morphological diversity of neotropical Anolis lizards, which originated in South America, colonized and radiated on various islands in the Caribbean, and then returned and diversified on the mainland. We focus in particular on what happens when mainland and island evolutionary radiations collide. We show that extensive continental radiations can result from island ancestors and that the incumbent and invading mainland clades achieve their ecological and morphological disparity in very different ways. Moreover, we show that when a mainland radiation derived from island ancestors comes into contact with an incumbent mainland radiation the ensuing interactions favor the island-derived clade.

Phylogenetic analysis of a new morphological dataset elucidates the evolutionary history of Crocodylia and resolves the long-standing gharial problem

First appearing in the latest Cretaceous, Crocodylia is a clade of semi-aquatic, predatory reptiles, defined by the last common ancestor of extant alligators, caimans, crocodiles, and gharials. Despite large strides in resolving crocodylian interrelationships over the last three decades, several outstanding problems persist in crocodylian systematics. Most notably, there has been persistent discordance between morphological and molecular datasets surrounding the affinities of the extant gharials, Gavialis gangeticus and Tomistoma schlegelii. Whereas molecular data consistently support a sister taxon relationship, in which they are more closely related to crocodylids than to alligatorids, morphological data indicate that Gavialis is the sister taxon to all other extant crocodylians. Here we present a new morphological dataset for Crocodylia based on a critical reappraisal of published crocodylian character data matrices and extensive firsthand observations of a global sample of crocodylians. This comprises the most taxonomically comprehensive crocodylian dataset to date (144 OTUs scored for 330 characters) and includes a new, illustrated character list with modifications to the construction and scoring of characters, and 46 novel characters. Under a maximum parsimony framework, our analyses robustly recover Gavialis as more closely related to Tomistoma than to other extant crocodylians for the first time based on morphology alone. This result is recovered regardless of the weighting strategy and treatment of quantitative characters. However, analyses using continuous characters and extended implied weighting (with high k-values) produced the most resolved, well-supported, and stratigraphically congruent topologies overall. Resolution of the gharial problem reveals that: (1) several gavialoids lack plesiomorphic features that formerly drew them towards the stem of Crocodylia; and (2) more widespread similarities occur between species traditionally divided into tomistomines and gavialoids, with these interpreted here as homology rather than homoplasy. There remains significant temporal incongruence regarding the inferred divergence timing of the extant gharials, indicating that several putative gavialids ('thoracosaurs') are incorrectly placed and require future re-appraisal. New alligatoroid interrelationships include: (1) support for a North American origin of Caimaninae in the latest Cretaceous; (2) the recovery of the early Paleogene South American taxon Eocaiman as a 'basal' alligatoroid; and (3) the paraphyly of the Cenozoic European taxon Diplocynodon. Among crocodyloids, notable results include modifications to the taxonomic content of Mekosuchinae, including biogeographic affinities of this clade with latest Cretaceous-early Paleogene Asian crocodyloids. In light of our new results, we provide a comprehensive review of the evolutionary and biogeographic history of Crocodylia, which included multiple instances of transoceanic and continental dispersal.

Origins of Galápagos’ land-locked vertebrates: what, whence, when, how?

Based on a synthesis of new molecular phylogenetic data, a detailed review is presented for the origins of the Galápagos’ native land-locked vertebrates [42 species; 11 clades: geckos (3), lava lizards (2), giant tortoises (1), iguanas (1), racer snakes (1) and oryzomyine rodents (3)]. Nine groups have roots in coastal Ecuador and Peru and would have been transported to the archipelago on rafts, many on the Humboldt Current. Inferring the sources of the giant tortoises, which probably floated over unaided, and the iguanas is more challenging because their closest living relatives occupy ground remote from the Pacific. Acknowledging uncertainties with the age-dating of both the phylogenetic tree nodes and the landmass emergences, seven, probably eight, of the colonizations likely involved beachings on the modern-day islands within the last 4 Myr. Three, possibly four, of the earlier arrivals may have been on now-submerged landmasses that were created by the Galápagos volcanic hotspot. Alternatively, the true sister taxa of the Galápagos species could be extinct and these colonizations, too, are more recent. This is likely for the giant tortoises. The assembled data set hints at the oldest/youngest clades showing the highest/lowest levels of diversification, although other factors also exert an influence.

Asymmetry of marine invasions across tropical oceans

Understanding the mechanisms of spatial variation of biological invasions, across local-to-global scales, has been a major challenge. The importance of evolutionary history for invasion dynamics was noted by Darwin, and several studies have since considered how biodiversity of source and recipient regions can influence the probability of invasions. For over a century, the Panama Canal has connected water bodies and biotas with different evolutionary histories, and created a global shipping hot spot, providing unique opportunities to test mechanisms that affect invasion patterns. Here, we test for asymmetry in both the extent of invasions and predation effects, a possible mechanism of biotic resistance, between two tropical oceans at similar latitudes. We estimated nonnative species (NNS) richness for sessile marine invertebrates, using standardized field surveys and literature synthesis, to examine whether invasions are asymmetrical, with more NNS present in the less diverse Pacific compared to the Atlantic. We also experimentally tested whether predation differentially limits the abundance and distribution of these invertebrates between oceans. In standardized surveys, observed total NNS richness was higher in the Pacific (18 NNS, 30% of all Pacific species) than the Atlantic (11 NNS, 13% of all Atlantic species). Similarly, literature-based records also display this asymmetry between coasts. When considering only the reciprocal exchange of NNS between Atlantic and Pacific biotas, NNS exchange from Atlantic to Pacific was eightfold higher than the opposite direction, exceeding the asymmetry predicted by random exchange based simply on differences of overall diversity per region. Predation substantially reduced biomass and changed NNS composition in the Pacific, but no such effects were detected on the Atlantic coast. Specifically, some dominant NNS were particularly susceptible to predation in the Pacific, supporting the hypothesis that predation may reduce the abundance of certain NNS here. These results are consistent with predictions that high diversity in source regions, and species interactions in recipient regions, shape marine invasion patterns. Our comparisons and experiments across two tropical ocean basins, suggest that global invasion dynamics are likely driven by both ecological and evolutionary factors that shape susceptibility to and directionality of invasions across biogeographic scales.

Autosomal, sex-linked and mitochondrial loci resolve evolutionary relationships among wrens in the genus Campylorhynchus

Although there is general consensus that sampling of multiple genetic loci is critical in accurate reconstruction of species trees, the exact numbers and the best types of molecular markers remain an open question. In particular, the phylogenetic utility of sex-linked loci is underexplored. Here, we sample all species and 70% of the named diversity of the New World wren genus Campylorhynchus using sequences from 23 loci, to evaluate the effects of linkage on efficiency in recovering a well-supported tree for the group. At a tree-wide level, we found that most loci supported fewer than half the possible clades and that sex-linked loci produced similar resolution to slower-coalescing autosomal markers, controlling for locus length. By contrast, we did find evidence that linkage affected the efficiency of recovery of individual relationships; as few as two sex-linked loci were necessary to resolve a selection of clades with long to medium subtending branches, whereas 4-6 autosomal loci were necessary to achieve comparable results. These results support an expanded role for sampling of the avian Z chromosome in phylogenetic studies, including target enrichment approaches. Our concatenated and species tree analyses represent significant improvements in our understanding of diversification in Campylorhynchus, and suggest a relatively complex scenario for its radiation across the Miocene/Pliocene boundary, with multiple invasions of South America.

Tectonic and climatic drivers of Asian monsoon evolution

Asian Monsoon rainfall supports the livelihood of billions of people, yet the relative importance of different drivers remains an issue of great debate. Here, we present 30 million-year model-based reconstructions of Indian summer monsoon and South East Asian monsoon rainfall at millennial resolution. We show that precession is the dominant direct driver of orbital variability, although variability on obliquity timescales is driven through the ice sheets. Orographic development dominated the evolution of the South East Asian monsoon, but Indian summer monsoon evolution involved a complex mix of contributions from orography (39%), precession (25%), atmospheric CO2 (21%), ice-sheet state (5%) and ocean gateways (5%). Prior to 15 Ma, the Indian summer monsoon was broadly stable, albeit with substantial orbital variability. From 15 Ma to 5 Ma, strengthening was driven by a combination of orography and glaciation, while closure of the Panama gateway provided the prerequisite for the modern Indian summer monsoon state through a strengthened Atlantic meridional overturning circulation.

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.

Phylogenomics and Historical Biogeography of Seahorses, Dragonets, Goatfishes, and Allies (Teleostei: Syngnatharia): Assessing Factors Driving Uncertainty in Biogeographic Inferences

The charismatic trumpetfishes, goatfishes, dragonets, flying gurnards, seahorses, and pipefishes encompass a recently defined yet extraordinarily diverse clade of percomorph fishes-the series Syngnatharia. This group is widely distributed in tropical and warm-temperate regions, with a great proportion of its extant diversity occurring in the Indo-Pacific. Because most syngnatharians feature long-range dispersal capabilities, tracing their biogeographic origins is challenging. Here, we applied an integrative phylogenomic approach to elucidate the evolutionary biogeography of syngnatharians. We built upon a recently published phylogenomic study that examined ultraconserved elements by adding 62 species (total 169 species) and one family (Draconettidae), to cover ca. 25% of the species diversity and all 10 families in the group. We inferred a set of time-calibrated trees and conducted ancestral range estimations. We also examined the sensitivity of these analyses to phylogenetic uncertainty (estimated from multiple genomic subsets), area delimitation, and biogeographic models that include or exclude the jump-dispersal parameter (j). Of the three factors examined, we found that the j parameter has the strongest effect in ancestral range estimates, followed by number of areas defined, and tree topology and divergence times. After accounting for these uncertainties, our results reveal that syngnatharians originated in the ancient Tethys Sea ca. 87 Ma (84-94 Ma; Late Cretaceous) and subsequently occupied the Indo-Pacific. Throughout syngnatharian history, multiple independent lineages colonized the eastern Pacific (6-8 times) and the Atlantic (6-14 times) from their center of origin, with most events taking place following an east-to-west route prior to the closure of the Tethys Seaway ca. 12-18 Ma. Ultimately, our study highlights the importance of accounting for different factors generating uncertainty in macroevolutionary and biogeographic inferences.

Mesoamerica is a cradle and the Atlantic Forest is a museum of Neotropical butterfly diversity: insights from the evolution and biogeography of Brassolini (Lepidoptera: Nymphalidae)

Regional species diversity is explained ultimately by speciation, extinction and dispersal. Here, we estimate dispersal and speciation rates of Neotropical butterflies to propose an explanation for the distribution and diversity of extant species. We focused on the tribe Brassolini (owl butterflies and allies), a Neotropical group that comprises 17 genera and 108 species, most of them endemic to rainforest biomes. We inferred a robust species tree using the multispecies coalescent framework and a dataset including molecular and morphological characters. This formed the basis for three changes in Brassolini classification: (1) Naropina syn. nov. is subsumed within Brassolina; (2) Aponarope syn. nov. is subsumed within Narope; and (3) Selenophanes orgetorix comb. nov. is reassigned from Catoblepia to Selenophanes. By applying biogeographical stochastic mapping, we found contrasting species diversification and dispersal dynamics across rainforest biomes, which might be explained, in part, by the geological and environmental history of each bioregion. Our results revealed a mosaic of biome-specific evolutionary histories within the Neotropics, where butterfly species have diversified rapidly (cradles: Mesoamerica), have accumulated gradually (museums: Atlantic Forest) or have diversified and accumulated alternately (Amazonia). Our study contributes evidence from a major butterfly lineage that the Neotropics are a museum and a cradle of species diversity.

Molecular Phylogeny and Evolution of Amazon Parrots in the Greater Antilles

Amazon parrots (Amazona spp.) colonized the islands of the Greater Antilles from the Central American mainland, but there has not been a consensus as to how and when this happened. Today, most of the five remaining island species are listed as endangered, threatened, or vulnerable as a consequence of human activity. We sequenced and annotated full mitochondrial genomes of all the extant Amazon parrot species from the Greater Antillean (A. leucocephala (Cuba), A. agilis, A. collaria (both from Jamaica), A. ventralis (Hispaniola), and A. vittata (Puerto Rico)), A. albifrons from mainland Central America, and A. rhodocorytha from the Atlantic Forest in Brazil. The assembled and annotated mitogenome maps provide information on sequence organization, variation, population diversity, and evolutionary history for the Caribbean species including the critically endangered A. vittata. Despite the larger number of available samples from the Puerto Rican Parrot Recovery Program, the sequence diversity of the A. vittata population in Puerto Rico was the lowest among all parrot species analyzed. Our data support the stepping-stone dispersal and speciation hypothesis that has started approximately 3.47 MYA when the ancestral population arrived from mainland Central America and led to diversification across the Greater Antilles, ultimately reaching the island of Puerto Rico 0.67 MYA. The results are presented and discussed in light of the geological history of the Caribbean and in the context of recent parrot evolution, island biogeography, and conservation. This analysis contributes to understating evolutionary history and empowers subsequent assessments of sequence variation and helps design future conservation efforts in the Caribbean.