The origin and evolution of coral species richness in a marine biodiversity hotspot

Earliest known funerary rites in Wallacea after the last glacial maximum

The insular region of Wallacea has become a focal point for studying Pleistocene human ecological and cultural adaptations in island environments, however, little is understood about early burial traditions during the Pleistocene. Here we investigate maritime interactions and burial practices at Ratu Mali 2, an elevated coastal cave site on the small island of Kisar in the Lesser Sunda Islands of eastern Indonesia dated to 15,500-3700 cal. BP. This multidisciplinary study demonstrates extreme marine dietary adaptations, engagement with an extensive exchange network across open seas, and early mortuary practices. A flexed male and a female, interred in a single grave with abundant shellfish and obsidian at Ratu Mali 2 by 14.7 ka are the oldest known human burials in Wallacea with established funerary rites. These findings highlight the impressive flexibility of our species in marginal environments and provide insight into the earliest known ritualised treatment of the dead in Wallacea.

Seq' and ARMS shall find: DNA (meta)barcoding of Autonomous Reef Monitoring Structures across the tree of life uncovers hidden cryptobiome of tropical urban coral reefs

Coral reefs are among the richest marine ecosystems on Earth, but there remains much diversity hidden within cavities of complex reef structures awaiting discovery. While the abundance of corals and other macroinvertebrates are known to influence the diversity of other reef-associated organisms, much remains unknown on the drivers of cryptobenthic diversity. A combination of standardized sampling with 12 units of the Autonomous Reef Monitoring Structure (ARMS) and high-throughput sequencing was utilized to uncover reef cryptobiome diversity across the equatorial reefs in Singapore. DNA barcoding and metabarcoding of mitochondrial cytochrome c oxidase subunit I, nuclear 18S and bacterial 16S rRNA genes revealed the taxonomic composition of the reef cryptobiome, comprising 15,356 microbial ASVs from over 50 bacterial phyla, and 971 MOTUs across 15 metazoan and 19 non-metazoan eukaryote phyla. Environmental factors across different sites were tested for relationships with ARMS diversity. Differences among reefs in diversity patterns of metazoans and other eukaryotes, but not microbial communities, were associated with biotic (coral cover) and abiotic (distance, temperature and sediment) environmental variables. In particular, ARMS deployed at reefs with higher coral cover had greater metazoan diversity and encrusting plate cover, with larger-sized non-coral invertebrates influencing spatial patterns among sites. Our study showed that DNA barcoding and metabarcoding of ARMS constitute a valuable tool for quantifying cryptobenthic diversity patterns and can provide critical information for the effective management of coral reef ecosystems.

A hybrid-capture approach to reconstruct the phylogeny of Scleractinia (Cnidaria: Hexacorallia)

A well-supported evolutionary tree representing most major lineages of scleractinian corals is in sight with the development and application of phylogenomic approaches. Specifically, hybrid-capture techniques are shedding light on the evolution and systematics of corals. Here, we reconstructed a broad phylogeny of Scleractinia to test previous phylogenetic hypotheses inferred from a few molecular markers, in particular, the relationships among major scleractinian families and genera, and to identify clades that require further research. We analysed 449 nuclear loci from 422 corals, comprising 266 species spanning 26 families, combining data across whole genomes, transcriptomes, hybrid capture and low-coverage sequencing to reconstruct the largest phylogenomic tree of scleractinians to date. Due to the large number of loci and data completeness (<38% missing data), node supports were high across shallow and deep nodes with incongruences observed in only a few shallow nodes. The "Robust" and "Complex" clades were recovered unequivocally, and our analyses confirmed that Micrabaciidae Vaughan, 1905 is sister to the "Robust" clade, transforming our understanding of the "Basal" clade. Several families remain polyphyletic in our phylogeny, including Deltocyathiidae Kitahara, Cairns, Stolarski & Miller, 2012, Caryophylliidae Dana, 1846, and Coscinaraeidae Benzoni, Arrigoni, Stefani & Stolarski, 2012, and we hereby formally proposed the family name Pachyseridae Benzoni & Hoeksema to accommodate Pachyseris Milne Edwards & Haime, 1849, which is phylogenetically distinct from Agariciidae Gray, 1847. Results also revealed species misidentifications and inconsistencies within morphologically complex clades, such as Acropora Oken, 1815 and Platygyra Ehrenberg, 1834, underscoring the need for reference skeletal material and topotypes, as well as the importance of detailed taxonomic work. The approach and findings here provide much promise for further stabilising the topology of the scleractinian tree of life and advancing our understanding of coral evolution.

Fast-growing species shape the evolution of reef corals

Ecological interactions are ubiquitous on tropical coral reefs, where sessile organisms coexist in limited space. Within these high-diversity systems, reef-building scleractinian corals form an intricate interaction network. The role of biotic interactions among reef corals is well established on ecological timescales. However, its potential effect on macroevolutionary patterns remains unclear. By analysing the rich fossil record of Scleractinia, we show that reef coral biodiversity experienced marked evolutionary rate shifts in the last 3 million years, possibly driven by biotic interactions. Our models suggest that there was an overwhelming effect of staghorn corals (family Acroporidae) on the fossil diversity trajectories of other coral groups. Staghorn corals showed an unparalleled spike in diversification during the Pleistocene. But surprisingly, their expansion was linked with increases in both extinction and speciation rates in other coral families, driving a nine-fold increase in lineage turnover. These results reveal a double-edged effect of diversity dependency on reef evolution. Given their fast growth, staghorn corals may have increased extinction rates via competitive interactions, while promoting speciation through their role as ecosystem engineers. This suggests that recent widespread human-mediated reductions in staghorn coral cover, may be disrupting the key macroevolutionary processes that established modern coral reef ecosystems.

The effect of biotic interactions among reef corals on macroevolutionary patterns is unclear. Here, the authors study the rich coral fossil record, finding that reef coral diversity experienced potentially biotic interaction-driven evolutionary rate changes, and that Staghorn corals affected fossil diversity trajectories of other coral groups.

Coral community composition and carbonate production in an urbanized seascape

Coastal urbanization causes environmental modifications that directly and indirectly influence the distribution and functioning of coral reefs. However, the capacity of urban infrastructure to support corals and vertically accrete is less understood. Here, we investigated if coral communities on reefs and seawalls in Singapore are distinct, and examined the environmental variables influencing coral carbonate production. Surveys at 22 sites yielded 134 coral species, with richness significantly higher on reefs. Coral cover and Shannon index did not differ between habitat types. Community composition was distinct between habitat types, with seawalls supporting a higher proportion of massive and thick-plating species. 'Distance from mainland' was the single most important variable influencing normalized carbonate production rates (a function of species-specific linear extension rate and skeletal bulk density and site coral cover), which were higher further from the mainland where human activity and development pressures were greater. Our results indicate that environmental filtering strongly shapes coral communities and may influence ecosystem functioning in Singapore's urbanized reef system. The findings will guide the management of reefs on increasingly urbanized coastlines.

Evolutionary History Drives Biogeographic Patterns of Coral Reef Resilience

Modern-day Indo-Pacific coral reefs are characterized by rapid recovery driven by pulses of coral recruitment, but Caribbean reefs exhibit low rates of recruitment and poor recovery following a wide range of disturbance events. The contrasting evolutionary history of coral taxa offers key insight into biogeographic patterns of coral resilience. Following the closure of the Isthmus of Panama approximately 2.8 million years ago, widespread extinction of Caribbean corals led to an evolutionary bottleneck that favored large and long-lived species with a relatively high reliance on asexual versus sexual reproduction. In contrast, adaptive radiation led to the evolution of superrecruiting tabular, digitate, and corymbose corals that drive the rapid recovery of modern-day Indo-Pacific reefs following disturbance. The dominance of branching growth forms and evolutionary absence of superrecruiting growth forms throughout the entire evolutionary history of the Caribbean (approximately 38 million years ago to present) may explain the exceptionally high recruitment rates on modern-day Indo-Pacific reefs and low historical recruitment on Caribbean reefs. The evolutionary history of the Caribbean coral reef-building taxa implies that, even with a reversal of ecosystem state, widespread recovery of Caribbean reefs may be limited.

Historical biogeography of the widespread macroalga Sargassum (Fucales, Phaeophyceae)

Sargassum is a cosmopolitan brown algal genus spanning the three ocean basins of the Atlantic, Pacific and Indian Oceans, inhabiting temperate, subtropical and tropical habitats. Sargassum has been postulated to have originated in the Oligocene epoch approximately 30 mya according to a broad phylogenetic analysis of brown macroalgae, but its diversification to become one of the most widespread and speciose macroalgal genera remains unclear. Here, we present a Bayesian molecular clock study, which analyzed data from the order Fucales of the brown algal crown radiation (BACR) group to reconstruct a time-calibrated phylogeny of the Sargassum clade. Our phylogeny included a total of 120 taxa with 99 Sargassum species sampled for three molecular markers - ITS-2, cox3 and rbcLS - calibrated with an unambiguous Sargassaceae fossil from between the lower and middle Miocene. The analysis revealed a much later origin of Sargassum than expected at about 6.7 mya, with the genus diversifying since approximately 4.3 mya. Current geographic distributions of Sargassum species were then analyzed in conjunction with the time-calibrated phylogeny using the dispersal-extinction-cladogenesis (DEC) model to estimate ancestral ranges of clades in the genus. Results strongly support origination of Sargassum in the Central Indo-Pacific (CIP) region with subsequent independent dispersal events into other marine realms. The longer history of diversification in the ancestral CIP range could explain the much greater diversity there relative to other marine areas today. Analyses of these dynamic processes, when fine-tuned to a higher spatial resolution, enable the identification of evolutionary hotspots and provide insights into long-term dispersal patterns.

Genomic signatures of host-associated divergence and adaptation in a coral-eating snail, Coralliophila violacea (Kiener, 1836)

The fluid nature of the ocean, combined with planktonic dispersal of marine larvae, lowers physical barriers to gene flow. However, divergence can still occur despite gene flow if strong selection acts on populations occupying different ecological niches. Here, we examined the population genomics of an ectoparasitic snail, Coralliophila violacea (Kiener 1836), that specializes on Porites corals in the Indo-Pacific. Previous genetic analyses revealed two sympatric lineages associated with different coral hosts. In this study, we examined the mechanisms promoting and maintaining the snails' adaptation to their coral hosts. Genome-wide single nucleotide polymorphism (SNP) data from type II restriction site-associated DNA (2b-RAD) sequencing revealed two differentiated clusters of C. violacea that were largely concordant with coral host, consistent with previous genetic results. However, the presence of some admixed genotypes indicates gene flow from one lineage to the other. Combined, these results suggest that differentiation between host-associated lineages of C. violacea is occurring in the face of ongoing gene flow, requiring strong selection. Indeed, 2.7% of all SNP loci were outlier loci (73/2,718), indicative of divergence with gene flow, driven by adaptation of each C. violacea lineage to their specific coral hosts.

From marine park to future genomic observatory? Enhancing marine biodiversity assessments using a biocode approach

Few tropical marine sites have been thoroughly characterised for their animal species, even though they constitute the largest proportion of multicellular diversity. A number of focused biodiversity sampling programmes have amassed immense collections to address this shortfall, but obstacles remain due to the lack of identification tools and large proportion of undescribed species globally. These problems can be partially addressed with DNA barcodes ("biocodes"), which have the potential to facilitate the estimation of species diversity and identify animals to named species via barcode databases. Here, we present the first results of what is intended to be a sustained, systematic study of the marine fauna of Singapore's first marine park, reporting more than 365 animal species, determined based on DNA barcodes and/or morphology represented by 931 specimens (367 zooplankton, 564 macrofauna including 36 fish). Due to the lack of morphological and molecular identification tools, only a small proportion could be identified to species solely based on either morphology (24.5%) or barcodes (24.6%). Estimation of species numbers for some taxa was difficult because of the lack of sufficiently clear barcoding gaps. The specimens were imaged and added to "Biodiversity of Singapore" (http://singapore.biodiversity.online), which now contains images for > 13,000 species occurring in the country.

The Origin and Diversification of the Hyperdiverse Flora in the Chocó Biogeographic Region

Extremely high levels of plant diversity in the American tropics are derived from multiple interactions between biotic and abiotic factors. Previous studies have focused on macro-evolutionary dynamics of the Tropical Andes, Amazonia, and Brazil's Cerrado and Atlantic forests during the last decade. Yet, other equally important Neotropical biodiversity hotspots have been severely neglected. This is particularly true for the Chocó region on the north-western coast of South and Central America. This geologically complex region is Earth's ninth most biodiverse hotspot, hosting approximately 3% of all known plant species. Here, we test Gentry's [1982a,b] hypothesis of a northern Andean-Central American Pleistocene origin of the Chocoan flora using phylogenetic reconstructions of representative plant lineages in the American tropics. We show that plant diversity in the Chocó is derived mostly from Andean immigrants. Contributions from more distant biogeographical areas also exist but are fewer. We also identify a strong floristic connection between the Chocó and Central America, revealed by multiple migrations into the Chocó during the last 5 Ma. The dated phylogenetic reconstructions suggest a Plio-Pleistocene onset of the extant Chocó flora. Taken together, these results support to a limited extend Gentry's hypothesis of a Pleistocene origin and of a compound assembly of the Chocoan biodiversity hotspot. Strong Central American-Chocoan floristic affinity may be partly explained by the accretion of a land mass derived from the Caribbean plate to north-western South America. Additional densely sampled phylogenies of Chocoan lineages also well represented across the Neotropics could enlighten the role of land mass movements through time in the assembly of floras in Neotropical biodiversity hotspots.

Interaction among ploidy, breeding system and lineage diversification

If particular traits consistently affect rates of speciation and extinction, broad macroevolutionary patterns can be interpreted as consequences of selection at high levels of the biological hierarchy. Identifying traits associated with diversification rates is difficult because of the wide variety of characters under consideration and the statistical challenges of testing for associations from comparative phylogenetic data. Ploidy (diploid vs polyploid states) and breeding system (self-incompatible vs self-compatible states) are both thought to be drivers of differential diversification in angiosperms. We fit 29 diversification models to extensive trait and phylogenetic data in Solanaceae and investigate how speciation and extinction rate differences are associated with ploidy, breeding system, and the interaction between these traits. We show that diversification patterns in Solanaceae are better explained by breeding system and an additional unobserved factor, rather than by ploidy. We also find that the most common evolutionary pathway to polyploidy in Solanaceae occurs via direct breakdown of self-incompatibility by whole genome duplication, rather than indirectly via breakdown followed by polyploidization. Comparing multiple stochastic diversification models that include complex trait interactions alongside hidden states enhances our understanding of the macroevolutionary patterns in plant phylogenies.

Hierarchical partitioning of fish diversity and scale-dependent environmental effects in tropical coastal ecosystems

The spatial structure of the fish diversity and site-scale and landscape-scale environmental effects were investigated across hierarchical levels in tropical coastal ecosystems. Total diversity (γ) was hierarchically partitioned into α and β components using both the additive and multiplicative methods. A model selection based on the AICc was applied to generalized linear mixed models relating diversity measures to environmental variables and including random effects for hierarchical levels and season. Short-term seasonal effects were negligible. Spatial effects were more relevant at the site level and negligible at the subregion level, due to the high spatial heterogeneity and the natural pooling of ecosystems, respectively. Site-scale environmental effects were more relevant at the subregion level, with eutrophic conditions (continental influence) favoring the species richness (α and γ) and higher absence of species (β^A) in oligotrophic conditions (marine influence). At the system level, the positive effect of the distance from the ocean on γ and higher β^A in oligotrophic conditions reinforced the positive continental influence on fish diversity. Environmental homogenization processes were most likely associated with the negative effect of the pasture cover on α at the system level, and γ and β^A at the site level. The negative effect of the forest cover on the later diversity measure evidenced its relevance to maintain richer but more similar assemblages, whereas the positive continental influence was most likely due to the loss of stenohaline marine species. This study evidenced that disentangling spatial, land use, and marine vs. continental effects on diversity components is critical to understand the primary determinants of the fish diversity in tropical coastal ecosystems.

Live fast, diversify non-adaptively: evolutionary diversification of exceptionally short-lived annual killifishes

Background

Adaptive radiations are triggered by ecological opportunity – the access to novel niche domains with abundant available resources that facilitate the formation of new ecologically divergent species. Therefore, as new species saturate niche space, clades experience a diversity-dependent slowdown of diversification over time. At the other extreme of the radiation continuum, non-adaptively radiating lineages undergo diversification with minimal niche differentiation when ‘spatial opportunity’ (i.e. areas with suitable ‘ancestral’ ecological conditions) is available. Traditionally, most research has focused on adaptive radiations, while empirical studies on non-adaptive radiations remain lagging behind. A prolific clade of African fish with extremely short lifespan (Nothobranchius killifish), show the key evolutionary features of a candidate non-adaptive radiation – primarily allopatric species with minimal niche and phenotypic divergence. Here, we test the hypothesis that Nothobranchius killifish have non-adaptively diversified. We employ phylogenetic modelling to investigate the tempo and mode of macroevolutionary diversification of these organisms.

Results

Nothobranchius diversification has proceeded with minor niche differentiation and minimal morphological disparity among allopatric species. Additionally, we failed to identify evidence for a role of body size or biogeography in influencing diversification rates. Diversification has been homogeneous within this genus, with the only hotspot of species-richness not resulting from rapid diversification. However, species in sympatry show higher disparity, which may have been caused by character displacement among coexisting species.

Conclusions

Nothobranchius killifish have proliferated following the tempo and mode of a non-adaptive radiation. Our study confirms that this exceptionally short-lived group have diversified with minimal divergent niche adaptation, while one group of coexisting species seems to have facilitated spatial overlap among these taxa via the evolution of ecological character displacement.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1344-0) contains supplementary material, which is available to authorized users.

Cryptic lineages in the Wolf Cardinalfish living in sympatry on remote coral atolls

Coral reef health and biodiversity is under threat worldwide due to rapid climate change. However, much of the inter- and intra-specific diversity of coral reefs are undescribed even in well studied taxa such as fish. Delimiting previously unrecognised diversity is important for understanding the processes that generate and sustain biodiversity in coral reef ecosystems and informing strategies for their conservation and management. Many taxa that inhabit geographically isolated coral reefs rely on self-recruitment for population persistence, providing the opportunity for the evolution of unique genetic lineages through divergent selection and reproductive isolation. Many such lineages in corals and fish are morphologically similar or indistinguishable. Here, we report the discovery and characterisation of cryptic lineages of the Wolf Cardinalfish, Cheilodipterus artus, from the coral atolls of northwest Australia using multiple molecular markers from mitochondrial (CO1 and D-loop) and nuclear (microsatellites) DNA. Concordant results from all markers identified two highly divergent lineages that are morphologically cryptic and reproductively isolated. These lineages co-occurred at daytime resting sites, but the relative abundance of each lineage was strongly correlated with wave exposure. It appears, therefore, that fish from each lineage are better adapted to different microhabitats. Such cryptic and ecologically based diversity appears to be common in these atolls and may well aid resilience of these systems. Our results also highlight that underwater surveys based on visual identification clearly underestimate biodiversity, and that a taxonomic revision of the Cheilodipterus genus is necessary.

Explaining the ocean's richest biodiversity hotspot and global patterns of fish diversity

For most marine organisms, species richness peaks in the Central Indo-Pacific region and declines longitudinally, a striking pattern that remains poorly understood. Here, we used phylogenetic approaches to address the causes of richness patterns among global marine regions, comparing the relative importance of colonization time, number of colonization events, and diversification rates (speciation minus extinction). We estimated regional richness using distributional data for almost all percomorph fishes (17 435 species total, including approximately 72% of all marine fishes and approximately 33% of all freshwater fishes). The high diversity of the Central Indo-Pacific was explained by its colonization by many lineages 5.3-34 million years ago. These relatively old colonizations allowed more time for richness to build up through in situ diversification compared to other warm-marine regions. Surprisingly, diversification rates were decoupled from marine richness patterns, with clades in low-richness cold-marine habitats having the highest rates. Unlike marine richness, freshwater diversity was largely derived from a few ancient colonizations, coupled with high diversification rates. Our results are congruent with the geological history of the marine tropics, and thus may apply to many other organisms. Beyond marine biogeography, we add to the growing number of cases where colonization and time-for-speciation explain large-scale richness patterns instead of diversification rates.