The Roles of Introgression and Climate Change in the Rise to Dominance of Acropora Corals

Standing genetic variation and introgression shape the cryptic radiation of Aquilegia in the mountains of Southwest China

Cryptic diversity in evolutionary radiation offers an excellent system for investigating the intricacies of evolutionary progress. Understanding the evolution of cryptic diversity is imperative for unraveling the hidden complexities of biodiversity. However, empirical evidence elucidating the mechanisms behind cryptic radiation remains limited, particularly in plants. Here, we focus on a monophyletic group of Aquilegia species mainly distributed in the mountains of Southwest China, one of the world's biodiversity hotspots. Using whole-genome resequencing of 158 individuals from 23 natural populations, we identify three to four paraphyletic lineages within each morphological species. Our findings reveal that 39 out of 43 detected instances of introgression occurred post-lineage formation. Identifying shared genomic regions indicates that the divergence of fixed singletons in lineages from morphological species A. kansuensis and A. rockii predates lineage formation, supporting a scenario where incomplete lineage sorting of standing variation contributes to morphological parallelism. Furthermore, strong positive correlations among genomic differentiation, divergence, and introgression suggest that standing variations and introgression from non-sister lineages contribute to the rapid genetic divergence. Our study illuminates the important roles of standing variations and introgression in plant cryptic radiation, advancing our understanding of the complex mechanisms behind the evolution of biodiversity in recent radiation events.

Living Coral Displays, Research Laboratories, and Biobanks as Important Reservoirs of Chemodiversity with Potential for Biodiscovery

Over the last decades, bioprospecting of tropical corals has revealed numerous bioactive compounds with potential for biotechnological applications. However, this search involves sampling in natural reefs, and this is currently hampered by multiple ethical and technological constraints. Living coral displays, research laboratories, and biobanks currently offer an opportunity to continue to unravel coral chemodiversity, acting as "Noah's Arks" that may continue to support the bioprospecting of molecules of interest. This issue is even more relevant if one considers that tropical coral reefs currently face unprecedent threats and irreversible losses that may impair the biodiscovery of molecules with potential for new products, processes, and services. Living coral displays provide controlled environments for studying corals and producing both known and new metabolites under varied conditions, and they are not prone to common bottlenecks associated with bioprospecting in natural coral reefs, such as loss of the source and replicability. Research laboratories may focus on a particular coral species or bioactive compound using corals that were cultured ex situ, although they may differ from wild conspecifics in metabolite production both in quantitative and qualitative terms. Biobanks collect and preserve coral specimens, tissues, cells, and/or information (e.g., genes, associated microorganisms), which offers a plethora of data to support the study of bioactive compounds' mode of action without having to cope with issues related to access, standardization, and regulatory compliance. Bioprospecting in these settings faces several challenges and opportunities. On one hand, it is difficult to ensure the complexity of highly biodiverse ecosystems that shape the production and chemodiversity of corals. On the other hand, it is possible to maximize biomass production and fine tune the synthesis of metabolites of interest under highly controlled environments. Collaborative efforts are needed to overcome barriers and foster opportunities to fully harness the chemodiversity of tropical corals before in-depth knowledge of this pool of metabolites is irreversibly lost due to tropical coral reefs' degradation.

A genomic approach to Porites (Anthozoa: Scleractinia) megadiversity from the Indo-Pacific

Porites corals are vital components of tropical reef ecosystems worldwide, serving as ecosystem engineers and hubs of biodiversity in shallow water coral reefs. Despite their ecological significance and the widespread use of Porites spp. as models for research, the richness and evolutionary relationships of species within the genus remain elusive. In this study, we analyzed genomic data from 330 colonies of Porites from 17 localities across the Indo-Pacific region based on the reduced representation genomic approach ezRAD. We retrieved 25,163 SNPs and provided a phylogenomic hypothesis for 29 nominal species and 10 unknown morphologies, recovering 15 deeply rooted molecular clades. Among these, 12 clades included samples corresponding to single distinct morphospecies. One did not match any nominal species. The remaining two clades comprised species complexes, which included various massive and encrusting morphologies commonly used in experimental biology. Within these complexes, we observed additional geographic or morphological structure, indicating complex evolutionary dynamics, possibly reflecting distinct species, isolated populations or hybridization. Additionally, a series of divergent samples underscored the importance of more sampling to define species boundaries and refine phylogenomic relationships. We also integrated our findings with previous phylogenetic datasets and their respective sampling localities, challenging traditional notions about Porites species geographic distributions. Overall, our findings indicate a need to revise past synonymies and to formally establish new species. A precise understanding of Porites species and their diversity and distributions is necessary for effective reef conservation and management.

Photosymbiosis shaped animal genome architecture and gene evolution as revealed in giant clams

Symbioses are major drivers of organismal diversification and phenotypic innovation. However, how long-term symbioses shape whole genome evolution in metazoans is still underexplored. Here, we use a giant clam (Tridacna maxima) genome to demonstrate how symbiosis has left complex signatures in an animal's genome. Giant clams thrive in oligotrophic waters by forming a remarkable association with photosymbiotic dinoflagellate algae. Genome-based demographic inferences uncover a tight correlation between T. maxima global population change and major paleoclimate and habitat shifts, revealing how abiotic and biotic factors may dictate T. maxima microevolution. Comparative analyses reveal genomic features that may be symbiosis-driven, including expansion and contraction of immunity-related gene families and a large proportion of lineage-specific genes. Strikingly, about 70% of the genome is composed of repetitive elements, especially transposable elements, most likely resulting from a symbiosis-adapted immune system. This work greatly enhances our understanding of genomic drivers of symbiosis that underlie metazoan evolution and diversification.

Divergence in Regulatory Regions and Gene Duplications May Underlie Chronobiological Adaptation in Desert Tortoises

Many cellular processes and organismal behaviours are time-dependent, and asynchrony of these phenomena can facilitate speciation through reinforcement mechanisms. The Mojave and Sonoran desert tortoises (Gopherus agassizii and G. morafkai respectively) reside in adjoining deserts with distinct seasonal rainfall patterns and they exhibit asynchronous winter brumation and reproductive behaviours. We used whole genome sequencing of 21 individuals from the two tortoise species and an outgroup to understand genes potentially underlying these characteristics. Genes within the most diverged 1% of the genome (FST ≥ 0.63) with putatively functional variation showed extensive divergence in regulatory elements, particularly promoter regions. Such genes related to UV nucleotide excision repair, mitonuclear and homeostasis functions. Genes mediating chronobiological (cell cycle, circadian and circannual) processes were also among the most highly diverged regions (e.g., XPA and ZFHX3). Putative promoter variants had significant enrichment of genes related to regulatory machinery (ARC-Mediator complex), suggesting that transcriptional cascades driven by regulatory divergence may underlie the behavioural differences between these species, leading to asynchrony-based prezygotic isolation. Further investigation revealed extensive expansion of respiratory and intestinal mucins (MUC5B and MUC5AC) within Gopherus, particularly G. morafkai. This expansion could be a xeric-adaptation to water retention and/or contribute to differential Mycoplasma agassizii infection rates between the two species, as mucins help clear inhaled dust and bacterial. Overall, results highlight the diverse array of genetic changes underlying divergence, adaptation and reinforcement during speciation.

Ancient Hybridisation Fuelled Diversification in Acropora Corals

Introgression is the infiltration or flow of genes from one species to another through hybridisation followed by backcrossing. This may lead to incorrect phylogenetic reconstruction or divergence-time estimation. Acropora is a dominant genus of reef-building corals; however, whether this group has an introgression history before their diversification remains unclear, and previous divergence-time estimates of Acropora have not considered the impact of introgression. In this study, we broke through the limitation of a few genes and a few species and proved the existence of ancient introgression in the evolution of Acropora from whole-genome protein-coding sequences. We inferred 21.9% of all triplet loci (homologous loci from three different species) with a history of introgression and a series of introgression events with a genetic material contribution of up to 30.9% before diversification. Furthermore, 7756 nuclear loci were clustered into three groups using a multidimensional scaling algorithm, the heterogeneity of which resulted in different phylogenetic relationships. The diversification time of Acropora was estimated to be middle to late Miocene when we retained only the gene group with the lowest degree of introgression. The collision of Australia with the Pacific arcs and the Southeast Asian margin in the early Miocene, and a series of cooling events in the middle to late Miocene, may provide geographical and climatic conditions for the diversification of Acropora, respectively. Therefore, our results indicate that at the genome-wide level, ancient introgressive hybridisation may have promoted the radiation evolution of Acropora. Based on our results, the influence of introgression should be taken into account when reconstructing phylogenetic relationships and evaluating divergence time.

Interspecific hybridisation provides a low-risk option for increasing genetic diversity of reef-building corals

Interspecific hybridisation increases genetic diversity and has played a significant role in the evolution of corals in the genus Acropora. In vitro fertilisation can be used to increase the frequency of hybridisation among corals, potentially enhancing their ability to adapt to climate change. Here, we assessed the field performance of hybrids derived from the highly cross-fertile coral species Acropora sarmentosa and Acropora florida from the Great Barrier Reef. Following outplanting to an inshore reef environment, the 10-month survivorship of the hybrid offspring groups was intermediate between that of the purebred groups, although not all pairwise comparisons were statistically significant. The A. florida purebreds, which had the lowest survivorship, were significantly larger at 10 months post-deployment compared to the other three groups. The four offspring groups harboured the same intracellular photosymbiont communities (Symbiodiniaceae), indicating that observed performance differences were due to the coral host and not photosymbiont communities. The limited differences in the performance of the groups and the lack of outbreeding depression of the F1 hybrids in the field suggest that interspecific hybridisation may be a useful method to boost the genetic diversity, and as such increase the adaptive capacity, of coral stock for restoration of degraded and potentially genetically eroded populations.

Integrative taxonomic analyses reveal that rapid genetic divergence drives Acropora speciation

Reef-building corals provide the structural basis for one of Earth's most spectacular and diverse but increasingly threatened ecosystems. The reef-building coral genus Acropora may have undergone substantial speciation during the Pleistocene climate and sea-level changes. Here, we aimed to evaluate the speciation history of four morphologically similar tabular Acropora species (Acropora aff. hyacinthus, A. cf. bifurcata, A. cf. cytherea, and A. cf. subulata) using an integrative approach with morphology, genetic, and reproduction methodology. Extensive morphological analyses showed that these four species are distinct and exhibited high gamete incompatibility, preventing hybridization. Furthermore, population structure and principal component analyses with SNPs (>60,000) indicated that these species were genetically distinct, and the ABBA-BABA test did not support introgression among these species. Many of their coding and noncoding RNA sequences showed high genetic variance at loci with high Fst values along the genome. Comparison of these orthologs with those of other Acropora species suggested that many of these genes are under positive selection, which could be associated with spawning time, gamete, and morphological divergence. Our findings show that the speciation of tabular Acropora occurred without hybridization, and the divergence accompanying the rapid evolution of genes in species-rich Acropora could be associated with speciation.

Unravelling the relationships among Madrepora Linnaeus, 1758, Oculina Lamark, 1816 and Cladocora Ehrenberg, 1834 (Cnidaria: Anthozoa: Scleractinia)

Despite the widespread use of integrative taxonomic approaches, many scleractinian coral genera and species remain grouped in polyphyletic families, classified as incertae sedis or simply understudied. Oculinidae Gray, 1847 represents a family for which many taxonomic questions remain unresolved, particularly those related to some of the current genera, such as Oculina Lamark, 1816 or recently removed genera, including Cladocora Ehrenberg, 1834 and Madrepora Linnaeus, 1758. Cladocora is currently assigned to the family Cladocoridae Milne Edwards & Haime, 1857 and a new family, Bathyporidae Kitahara, Capel, Zilberberg & Cairns, 2024, was recently raised to accommodate Madrepora . However, the name Bathyporidae is not valid because this was not formed on the basis of a type genus name. To resolve taxonomic questions related to these three genera, the evolutionary relationships are explored through phylogenetic analyses of 18 molecular markers. The results of these analyses support a close relationship between the species Oculina patagonica and Cladocora caespitosa , indicating that these may belong to the same family (and possibly genus), and highlighting the need for detailed revisions of Oculina and Cladocora . By contrast, a distant relationship is found between these two species and Madrepora oculata , with the overall evidence supporting the placement of Madrepora in the resurrected family Madreporidae Ehrenberg, 1834. This study advances our knowledge of coral systematics and highlights the need for a comprehensive review of the genera Oculina , Cladocora and Madrepora .

Exploring coral speciation: Multiple sympatric Stylophora pistillata taxa along a divergence continuum on the Great Barrier Reef

Understanding how biodiversity originates and is maintained are fundamental challenge in evolutionary biology. Speciation is a continuous process and progression along this continuum depends on the interplay between evolutionary forces driving divergence and forces promoting genetic homogenisation. Coral reefs are broadly connected yet highly heterogeneous ecosystems, and divergence with gene flow at small spatial scales might therefore be common. Genomic studies are increasingly revealing the existence of closely related and sympatric taxa within taxonomic coral species, but the extent to which these taxa might still be exchanging genes and sharing environmental niches is unclear. In this study, we sampled extensively across diverse habitats at multiple reefs of the Great Barrier Reef (GBR) and comprehensively examined genome-wide diversity and divergence histories within and among taxa of the Stylophora pistillata species complex. S. pistillata is one of the most abundant and well-studied coral species, yet we discovered five distinct taxa, with wide geographic ranges and extensive sympatry. Demographic modelling showed that speciation events have occurred with gene flow and that taxa are at different stages along a divergence continuum. We found significant correlations between genetic divergence and specific environmental variables, suggesting that niche partitioning may have played a role in speciation and that S. pistillata taxa might be differentially adapted to different environments. Conservation actions rely on estimates of species richness, population sizes and species ranges, which are biased if divergent taxa are lumped together. As coral reefs are rapidly degrading due to climate change, our study highlights the importance of recognising evolutionarily distinct and differentially adapted coral taxa to improve conservation and restoration efforts aiming at protecting coral genetic diversity.

Extensive introgression among strongylocentrotid sea urchins revealed by phylogenomics

Gametic isolation is thought to play an important role in the evolution of reproductive isolation in broadcast-spawning marine invertebrates. However, it is unclear whether gametic isolation commonly evolves early in the speciation process or only accumulates after other reproductive barriers are already in place. It is also unknown whether gametic isolation is an effective barrier to introgression following speciation. Here, we used whole-genome sequencing data and multiple complementary phylogenomic approaches to test whether the well-documented gametic incompatibilities among the strongylocentrotid sea urchins have limited introgression. We quantified phylogenetic discordance, inferred reticulate phylogenetic networks, and applied the Δ statistic using gene tree topologies reconstructed from multiple sequence alignments of protein-coding single-copy orthologs. In addition, we conducted ABBA-BABA tests on genome-wide single nucleotide variants and reconstructed a phylogeny of mitochondrial genomes. Our results revealed strong mito-nuclear discordance and considerable nonrandom gene tree discordance that cannot be explained by incomplete lineage sorting alone. Eight of the nine species examined demonstrated a history of introgression with at least one other species or ancestral lineage, indicating that introgression was common during the diversification of the strongylocentrotid urchins. There was strong support for introgression between four extant species pairs (Strongylocentrotus pallidus ⇔ S. droebachiensis, S. intermedius ⇔ S. pallidus, S. purpuratus ⇔ S. fragilis, and Mesocentrotus franciscanus ⇔ Pseudocentrotus depressus) and additional evidence for introgression on internal branches of the phylogeny. Our results suggest that the existing gametic incompatibilities among the strongylocentrotid urchin species have not been a complete barrier to hybridization and introgression following speciation. Their continued divergence in the face of widespread introgression indicates that other reproductive isolating barriers likely exist and may have been more critical in establishing reproductive isolation early in speciation.

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.

Power of Bayesian and Heuristic Tests to Detect Cross-Species Introgression with Reference to Gene Flow in the Tamias quadrivittatus Group of North American Chipmunks

In the past two decades, genomic data have been widely used to detect historical gene flow between species in a variety of plants and animals. The Tamias quadrivittatus group of North America chipmunks, which originated through a series of rapid speciation events, are known to undergo massive amounts of mitochondrial introgression. Yet in a recent analysis of targeted nuclear loci from the group, no evidence for cross-species introgression was detected, indicating widespread cytonuclear discordance. The study used the heuristic method HYDE to detect gene flow, which may suffer from low power. Here we use the Bayesian method implemented in the program BPP to re-analyze these data. We develop a Bayesian test of introgression, calculating the Bayes factor via the Savage-Dickey density ratio using the Markov chain Monte Carlo (MCMC) sample under the model of introgression. We take a stepwise approach to constructing an introgression model by adding introgression events onto a well-supported binary species tree. The analysis detected robust evidence for multiple ancient introgression events affecting the nuclear genome, with introgression probabilities reaching 63%. We estimate population parameters and highlight the fact that species divergence times may be seriously underestimated if ancient cross-species gene flow is ignored in the analysis. We examine the assumptions and performance of HYDE and demonstrate that it lacks power if gene flow occurs between sister lineages or if the mode of gene flow does not match the assumed hybrid-speciation model with symmetrical population sizes. Our analyses highlight the power of likelihood-based inference of cross-species gene flow using genomic sequence data. [Bayesian test; BPP; chipmunks; introgression; MSci; multispecies coalescent; Savage-Dickey density ratio.].

Using genomics, morphometrics, and environmental niche modeling to test the validity of a narrow-range endemic snail, Pateranantahala (Gastropoda, Polygyridae)

Terrestrial gastropods are among the most imperiled groups of organisms on Earth. Many species have a complex taxonomic history, often including poorly defined subspecies, most of which have not been the focus of modern systematics research. Genomic tools, geometric morphometrics, and environmental niche modeling were used to assess the taxonomic status of Pateraclarkiinantahala (Clench & Banks, 1932), a subspecies of high conservation concern with a restricted range of approximately 3.3 km² in North Carolina, USA. A genome-scale dataset was generated that included individuals with morphologies matching P.c.nantahala, P.c.clarkii, and one individual with an intermediate form between P.c.nantahala and P.c.clarkii that was initially hypothesized as a potential hybrid. Mitochondrial phylogenetics, nuclear species tree inference, and phylogenetic networks were used to assess relationships and gene flow. Differences in shell shape via geometric morphometrics and whether the environmental niches of the two subspecies were significantly different were also examined. Molecular analyses indicated an absence of gene flow among lineages of P.clarkii sensu lato. Analyses rejected our hypothesis that the intermediate shelled form represented a hybrid, but instead indicated that it was a distinct lineage. Environmental niche models indicated significant differences in environmental niche between P.c.clarkii and P.c.nantahala, and geometric morphometrics indicated that P.c.nantahala had a significantly different shell shape. Given multiple lines of evidence, species-level recognition of P.nantahala is warranted.

Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef

Genomic studies are uncovering extensive cryptic diversity within reef-building corals, suggesting that evolutionarily and ecologically relevant diversity is highly underestimated in the very organisms that structure coral reefs. Furthermore, endosymbiotic algae within coral host species can confer adaptive responses to environmental stress and may represent additional axes of coral genetic variation that are not constrained by taxonomic divergence of the cnidarian host. Here, we examine genetic variation in a common and widespread, reef-building coral, Acropora tenuis, and its associated endosymbiotic algae along the entire expanse of the Great Barrier Reef (GBR). We use SNPs derived from genome-wide sequencing to characterize the cnidarian coral host and organelles from zooxanthellate endosymbionts (genus Cladocopium). We discover three distinct and sympatric genetic clusters of coral hosts, whose distributions appear associated with latitude and inshore-offshore reef position. Demographic modelling suggests that the divergence history of the three distinct host taxa ranges from 0.5 to 1.5 million years ago, preceding the GBR's formation, and has been characterized by low-to-moderate ongoing inter-taxon gene flow, consistent with occasional hybridization and introgression typifying coral evolution. Despite this differentiation in the cnidarian host, A. tenuis taxa share a common symbiont pool, dominated by the genus Cladocopium (Clade C). Cladocopium plastid diversity is not strongly associated with host identity but varies with reef location relative to shore: inshore colonies contain lower symbiont diversity on average but have greater differences between colonies as compared with symbiont communities from offshore colonies. Spatial genetic patterns of symbiont communities could reflect local selective pressures maintaining coral holobiont differentiation across an inshore-offshore environmental gradient. The strong influence of environment (but not host identity) on symbiont community composition supports the notion that symbiont community composition responds to habitat and may assist in the adaptation of corals to future environmental change.

Not out of the Mediterranean: Atlantic populations of the gorgonian Paramuricea clavata are a separate sister species under further lineage diversification

The accurate delimitation of species boundaries in nonbilaterian marine taxa is notoriously difficult, with consequences for many studies in ecology and evolution. Anthozoans are a diverse group of key structural organisms worldwide, but the lack of reliable morphological characters and informative genetic markers hampers our ability to understand species diversification. We investigated population differentiation and species limits in Atlantic (Iberian Peninsula) and Mediterranean lineages of the octocoral genus Paramuricea previously identified as P. clavata. We used a diverse set of molecular markers (microsatellites, RNA-seq derived single-copy orthologues [SCO] and mt-mutS [mitochondrial barcode]) at 49 locations. Clear segregation of Atlantic and Mediterranean lineages was found with all markers. Species-tree estimations based on SCO strongly supported these two clades as distinct, recently diverged sister species with incomplete lineage sorting, P. cf. grayi and P. clavata, respectively. Furthermore, a second putative (or ongoing) speciation event was detected in the Atlantic between two P. cf. grayi color morphotypes (yellow and purple) using SCO and supported by microsatellites. While segregating P. cf. grayi lineages showed considerable geographic structure, dominating circalittoral communities in southern (yellow) and western (purple) Portugal, their occurrence in sympatry at some localities suggests a degree of reproductive isolation. Overall, our results show that previous molecular and morphological studies have underestimated species diversity in Paramuricea occurring in the Iberian Peninsula, which has important implications for conservation planning. Finally, our findings validate the usefulness of phylotranscriptomics for resolving evolutionary relationships in octocorals.

A decade of population genetics studies of scleractinian corals: A systematic review

Coral reefs are the most diverse marine ecosystems. However, coral cover has decreased worldwide due to natural disturbances, climate change, and local anthropogenic drivers. In recent decades, various genetic methods and molecular markers have been developed to assess genetic diversity, structure, and connectivity in different coral species to determine the vulnerability of their populations. This review aims to identify population genetic studies of scleractinian corals in the last decade (2010-2020), and the techniques and molecular markers used. Bibliometric analysis was conducted to identify journals and authors working in this field. We then calculated the number of genetic studies by species and ecoregion based on data obtained from 178 studies found in Scopus and Web of Science. Coral Reefs and Molecular Ecology were the main journals published population genetics studies, and microsatellites are the most widely used molecular markers. The Caribbean, Australian Barrier Reef, and South Kuroshio in Japan are among the ecoregions with the most population genetics data. In contrast, we found limited information about the Coral Triangle, a region with the highest biodiversity and key to coral reef conservation. Notably, only 117 (out of 1500 described) scleractinian coral species have genetic studies. This review emphasizes which coral species have been studied and highlights remaining gaps and locations where such data is critical for coral conservation.

Evolutionary responses of a reef-building coral to climate change at the end of the last glacial maximum

Climate change threatens the survival of coral reefs on a global scale, primarily through mass bleaching and mortality as a result of marine heatwaves. While these short-term effects are clear, predicting the fate of coral reefs over the coming century is a major challenge. One way to understand the longer-term effects of rapid climate change is to examine the response of coral populations to past climate shifts. Coastal and shallow-water marine ecosystems such as coral reefs have been reshaped many times by sea-level changes during the Pleistocene, yet, few studies have directly linked this with its consequences on population demographics, dispersal, and adaptation. Here we use powerful analytical techniques, afforded by haplotype phased whole-genomes, to establish such links for the reef-building coral, Acropora digitifera. We show that three genetically distinct populations are present in northwestern Australia, and that their rapid divergence since the last glacial maximum (LGM) can be explained by a combination of founder-effects and restricted gene flow. Signatures of selective sweeps, too strong to be explained by demographic history, are present in all three populations and overlap with genes that show different patterns of functional enrichment between inshore and offshore habitats. In contrast to rapid divergence in the host, we find that photosymbiont communities are largely undifferentiated between corals from all three locations, spanning almost 1000 km, indicating that selection on host genes and not acquisition of novel symbionts, has been the primary driver of adaptation for this species in northwestern Australia.

A polar bear paleogenome reveals extensive ancient gene flow from polar bears into brown bears

Polar bears (Ursus maritimus) and brown bears (Ursus arctos) are sister species possessing distinct physiological and behavioural adaptations that evolved over the last 500,000 years. However, comparative and population genomics analyses have revealed that several extant and extinct brown bear populations have relatively recent polar bear ancestry, probably as the result of geographically localized instances of gene flow from polar bears into brown bears. Here, we generate and analyse an approximate 20X paleogenome from an approximately 100,000-year-old polar bear that reveals a massive prehistoric admixture event, which is evident in the genomes of all living brown bears. This ancient admixture event was not visible from genomic data derived from living polar bears. Like more recent events, this massive admixture event mainly involved unidirectional gene flow from polar bears into brown bears and occurred as climate changes caused overlap in the ranges of the two species. These findings highlight the complex reticulate paths that evolution can take within a regime of radically shifting climate.

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.

Spatially varying selection between habitats drives physiological shifts and local adaptation in a broadcast spawning coral on a remote atoll in Western Australia

At the Rowley Shoals in Western Australia, the prominent reef flat becomes exposed on low tide and the stagnant water in the shallow atoll lagoons heats up, creating a natural laboratory for characterizing the mechanisms of coral resilience to climate change. To explore these mechanisms in the reef coral Acropora tenuis, we collected samples from lagoon and reef slope habitats and combined whole-genome sequencing, ITS2 metabarcoding, experimental heat stress, and transcriptomics. Despite high gene flow across the atoll, we identified clear shifts in allele frequencies between habitats at relatively small linked genomic islands. Common garden heat stress assays showed corals from the lagoon to be more resistant to bleaching, and RNA sequencing revealed marked differences in baseline levels of gene expression between habitats. Our results provide new insight into the complex mechanisms of coral resilience to climate change and highlight the potential for spatially varying selection across complex coral reef seascapes to drive pronounced ecological divergence in climate-related traits.

First evidence for backcrossing of F1 hybrids in Acropora corals under sperm competition

Acropora is a species-rich genus of reef-building corals with highly diverse morphologies. Hybridization among intercrossing species potentially influences species diversity within Acropora. However, the mechanisms that allow hybridization/backcrossing remain unknown. Although we tested a limited number of species, we hypothesized that Acropora gametes in the Indo-Pacific may preferentially fertilize conspecific gametes despite their compatibility with heterospecific gametes, leading to infrequent hybridization between potentially intercrossing species. In this study, F₁ hybrids of Acropora florida and A. intermedia showed specific fertilization trends. For example, sperm had the ability to backcross with the parental species even in the presence of sperm from the parental species. Also, eggs of the hybrids produced from A. florida eggs and A. intermedia sperm (“FLOint”) exhibited self-fertilization. Since a low ratio of hybridization between A. florida and A. intermedia is predicted, the population size of hybrids should be small. Therefore, self-fertilization would facilitate reproduction of the hybrid in nature, while remaining sperm could outcompete parental species sperm to backcross with eggs. Although we succeeded in breeding two colonies of hybrids, it is reasonable to speculate that hybrids show a high tendency to choose the most efficient sexual reproduction tactics.

Genetic divergence and range expansion in a western North Pacific coral

Coral poleward range expansions have recently been observed in response to warming oceans. Range expansion can lead to reduced genetic diversity and increased frequency of deleterious mutations that were rare in core populations, potentially limiting the ability for adaptation and persistence in novel environments. Successful expansions that overcome these founder effects and colonize new habitat have been attributed to multiple introductions from different sources, hybridization with native populations, or rapid adaptive evolution. Here, we investigate population genomic patterns of the reef-building coral Acropora hyacinthus along a latitudinal cline that includes a well-established range expansion front in Japan using 2b-RAD sequencing. A total of 184 coral samples were collected across seven sites spanning from ~24°N to near its northern range front at ~33°N. We uncover the presence of three cryptic lineages of A. hyacinthus, which occupy discrete reefs within this region. Only one lineage is present along the expansion front and we find evidence for its historical occupation of marginal habitats. Within this lineage we also find evidence of bottleneck pressures associated with expansion events including higher clonality, increased linkage disequilibrium, and lower genetic diversity in range edge populations compared to core populations. Asymmetric migration between populations was also detected with lower migration from edge sites. Lastly, we describe genomic signatures of local adaptation potentially attributed to lower winter temperatures experienced at the more recently expanded northern populations. Together these data illuminate the genomic consequences of range expansion in a coral and highlight how adaptation to discrete environments along expansion fronts may facilitate further range expansion in this temperate coral lineage.

Deep connections: divergence histories with gene flow in mesophotic Agaricia corals

Largely understudied, mesophotic coral ecosystems lie below shallow reefs (at >30 m depth) and comprise ecologically distinct communities. Brooding reproductive modes appear to predominate among mesophotic‐specialist corals and may limit genetic connectivity among populations. Using reduced representation genomic sequencing, we assessed spatial population genetic structure at 50 m depth in an ecologically important mesophotic‐specialist species Agaricia grahamae, among locations in the Southern Caribbean. We also tested for hybridisation with the closely related (but depth‐generalist) species Agaricia lamarcki, within their sympatric depth zone (50 m). In contrast to our expectations, no spatial genetic structure was detected between the reefs of Curaçao and Bonaire (~40 km apart) within A. grahamae. However, cryptic taxa were discovered within both taxonomic species, with those in A. lamarcki (incompletely) partitioned by depth and those in A. grahamae occurring sympatrically (at the same depth). Hybrid analyses and demographic modelling identified contemporary and historical gene flow among cryptic taxa, both within and between A. grahamae and A. lamarcki. These results (1) indicate that spatial connectivity and subsequent replenishment may be possible between islands of moderate geographic distances for A. grahamae, an ecologically important mesophotic species, (2) that cryptic taxa occur in the mesophotic zone and environmental selection along shallow to mesophotic depth gradients may drive divergence in depth‐generalists such as A. lamarcki, and (3) highlight that gene flow links taxa within this relativity diverse Caribbean genus.

The Role of Photobionts as Drivers of Diversification in an Island Radiation of Lichen-Forming Fungi

Speciation in oceanic islands has attracted the interest of scientists since the 19th century. One of the most striking evolutionary phenomena that can be studied in islands is adaptive radiation, that is, when a lineage gives rise to different species by means of ecological speciation. Some of the best-known examples of adaptive radiation are charismatic organisms like the Darwin finches of the Galapagos and the cichlid fishes of the great African lakes. In these and many other examples, a segregation of the trophic niche has been shown to be an important diversification driver. Radiations are known in other groups of organisms, such as lichen-forming fungi. However, very few studies have investigated their adaptive nature, and none have focused on the trophic niche. In this study, we explore the role of the trophic niche in a putative radiation of endemic species from the Macaronesian Region, the Ramalina decipiens group. The photobiont diversity was studied by Illumina MiSeq sequencing of the ITS2 region of 197 specimens spanning the phylogenetic breadth and geographic range of the group. A total of 66 amplicon sequence variants belonging to the four main clades of the algal genus Trebouxia were found. Approximately half of the examined thalli showed algal coexistence, but in most of them, a single main photobiont amounted to more than 90% of the reads. However, there were no significant differences in photobiont identity and in the abundance of ITS2 reads across the species of the group. We conclude that a segregation of the trophic niche has not occurred in the R. decipiens radiation.

Divergence and hybridization in sea turtles: Inferences from genome data show evidence of ancient gene flow between species

Reconstructing past events of hybridization and population size changes are required to understand speciation mechanisms and current patterns of genetic diversity, and ultimately contribute to species' conservation. Sea turtles are ancient species currently facing anthropogenic threats including climate change, fisheries, and illegal hunting. Five of the seven extant sea turtle species are known to currently hybridize, especially along the Brazilian coast where some populations can have ~32%-42% of hybrids. Although frequently observed today, it is not clear what role hybridization plays in the evolutionary diversification of this group of reptiles. In this study, we generated whole genome resequencing data of the five globally distributed sea turtle species to estimate a calibrated phylogeny and the population size dynamics, and to understand the role of hybridization in shaping the genomes of these ancient species. Our results reveal discordant species divergence dates between mitochondrial and nuclear genomes, with a high frequency of conflicting trees throughout the nuclear genome suggesting that some sea turtle species frequently hybridized in the past. The reconstruction of the species' demography showed a general decline in effective population sizes with no signs of recovery, except for the leatherback sea turtle. Furthermore, we discuss the influence of reference bias in our estimates. We show long-lasting ancestral gene flow events within Chelonioidea that continued for millions of years after initial divergence. Speciation with gene flow is a common pattern in marine species, and it raises questions whether current hybridization events should be considered as a part of these species' evolutionary history or a conservation issue.

Genomic analysis of distinct bleaching tolerances among cryptic coral species

Reef-building coral species are experiencing an unprecedented decline owing to increasing frequency and intensity of marine heatwaves and associated bleaching-induced mortality. Closely related species from the Acropora hyacinthus species complex differ in heat tolerance and in their association with heat-tolerant symbionts. We used low-coverage full genome sequencing of 114 colonies monitored across the 2015 bleaching event in American Samoa to determine the genetic differences among four cryptic species (termed HA, HC, HD and HE) that have diverged in these species traits. Cryptic species differed strongly at thousands of single nucleotide polymorphisms across the genome which are enriched for amino acid changes in the bleaching-resistant species HE. In addition, HE also showed two particularly divergent regions with strong signals of differentiation. One approximately 220 kb locus, HES1, contained the majority of fixed differences in HE. A second locus, HES2, was fixed in HE but polymorphic in the other cryptic species. Surprisingly, non-HE individuals with HE-like haplotypes at HES2 were more likely to bleach. At both loci, HE showed particular sequence similarity to a congener, Acropora millepora. Overall, resilience to bleaching during the third global bleaching event was strongly structured by host cryptic species, buoyed by differences in symbiont associations between these species.

Elucidation of the speciation history of three sister species of crown-of-thorns starfish (Acanthaster spp.) based on genomic analysis

The crown-of-thorns starfish (COTS) is a coral predator that is widely distributed in Indo-Pacific Oceans. A previous phylogenetic study using partial mitochondrial sequences suggested that COTS had diverged into four distinct species, but a nuclear genome-based analysis to confirm this was not conducted. To address this, COTS species nuclear genome sequences were analysed here, sequencing Northern Indian Ocean (NIO) and Red Sea (RS) species genomes for the first time, followed by a comparative analysis with the Pacific Ocean (PO) species. Phylogenetic analysis and ADMIXTURE analysis revealed clear divergences between the three COTS species. Furthermore, within the PO species, the phylogenetic position of the Hawaiian sample was further away from the other Pacific-derived samples than expected based on the mitochondrial data, suggesting that it may be a PO subspecies. The pairwise sequentially Markovian coalescent model showed that the trajectories of the population size diverged by region during the Mid-Pleistocene transition when the sea-level was dramatically decreased, strongly suggesting that the three COTS species experienced allopatric speciation. Analysis of the orthologues indicated that there were remarkable genes with species-specific positive selection in the genomes of the PO and RS species, which suggested that there may be local adaptations in the COTS species.

Evolution and connectivity influence the persistence and recovery of coral reefs under climate change in the Caribbean, Southwest Pacific, and Coral Triangle

Corals are experiencing unprecedented decline from climate change-induced mass bleaching events. Dispersal not only contributes to coral reef persistence through demographic rescue but can also hinder or facilitate evolutionary adaptation. Locations of reefs that are likely to survive future warming therefore remain largely unknown, particularly within the context of both ecological and evolutionary processes across complex seascapes that differ in temperature range, strength of connectivity, network size, and other characteristics. Here, we used eco-evolutionary simulations to examine coral adaptation to warming across reef networks in the Caribbean, the Southwest Pacific, and the Coral Triangle. We assessed the factors associated with coral persistence in multiple reef systems to understand which results are general and which are sensitive to particular geographic contexts. We found that evolution can be critical in preventing extinction and facilitating the long-term recovery of coral communities in all regions. Furthermore, the strength of immigration to a reef (destination strength) and current sea surface temperature robustly predicted reef persistence across all reef networks and across temperature projections. However, we found higher initial coral cover, slower recovery, and more evolutionary lag in the Coral Triangle, which has a greater number of reefs and more larval settlement than the other regions. We also found the lowest projected future coral cover in the Caribbean. These findings suggest that coral reef persistence depends on ecology, evolution, and habitat network characteristics, and that, under an emissions stabilization scenario (RCP 4.5), recovery may be possible over multiple centuries.

Eighteen Coral Genomes Reveal the Evolutionary Origin of Acropora Strategies to Accommodate Environmental Changes

The genus Acropora comprises the most diverse and abundant scleractinian corals (Anthozoa, Cnidaria) in coral reefs, the most diverse marine ecosystems on Earth. However, the genetic basis for the success and wide distribution of Acropora are unknown. Here, we sequenced complete genomes of 15 Acropora species and 3 other acroporid taxa belonging to the genera Montipora and Astreopora to examine genomic novelties that explain their evolutionary success. We successfully obtained reasonable draft genomes of all 18 species. Molecular dating indicates that the Acropora ancestor survived warm periods without sea ice from the mid or late Cretaceous to the Early Eocene and that diversification of Acropora may have been enhanced by subsequent cooling periods. In general, the scleractinian gene repertoire is highly conserved; however, coral- or cnidarian-specific possible stress response genes are tandemly duplicated in Acropora. Enzymes that cleave dimethlysulfonioproprionate into dimethyl sulfide, which promotes cloud formation and combats greenhouse gasses, are the most duplicated genes in the Acropora ancestor. These may have been acquired by horizontal gene transfer from algal symbionts belonging to the family Symbiodiniaceae, or from coccolithophores, suggesting that although functions of this enzyme in Acropora are unclear, Acropora may have survived warmer marine environments in the past by enhancing cloud formation. In addition, possible antimicrobial peptides and symbiosis-related genes are under positive selection in Acropora, perhaps enabling adaptation to diverse environments. Our results suggest unique Acropora adaptations to ancient, warm marine environments and provide insights into its capacity to adjust to rising seawater temperatures.

Genomic variation of an endosymbiotic dinoflagellate (Symbiodinium 'fitti') among closely related coral hosts

Mutualisms where hosts are coupled metabolically to their symbionts often exhibit high partner fidelity. Most reef-building coral species form obligate symbioses with a specific species of photosymbionts, dinoflagellates in the family Symbiodiniaceae, despite needing to acquire symbionts early in their development from environmental sources. Three Caribbean acroporids (Acropora palmata, A. cervicornis and their F₁ hybrid) are sympatric across much of their range, but often occupy different depth and light habitats. Throughout this range, both species and their hybrid associate with the endosymbiotic dinoflagellate Symbiodinium 'fitti'. Because light (and therefore depth) influences the physiology of dinoflagellates, we investigated whether S. 'fitti' populations from each host taxon were differentiated genetically. Single nucleotide polymorphisms (SNPs) among S. 'fitti' strains were identified by aligning shallow metagenomic sequences of acroporid colonies sampled from across the Caribbean to a ~600-Mb draft assembly of the S. 'fitti' genome (from the CFL14120 A. cervicornis metagenome). Phylogenomic and multivariate analyses revealed that genomic variation among S. 'fitti' strains partitioned to each host taxon rather than by biogeographical origin. This is particularly noteworthy because the hybrid has a sparse fossil record and may be of relatively recent origin. A subset (37.6%) of the SNPs putatively under selection were nonsynonymous mutations predicted to alter protein efficiency. Differences in genomic variation of S. 'fitti' strains from each host taxon may reflect the unique selection pressures created by the microenvironments associated with each host. The nonrandom sorting among S. 'fitti' strains to different hosts could be the basis for lineage diversification via disruptive selection, leading to ecological specialization and ultimately speciation.

Morphological stasis masks ecologically divergent coral species on tropical reefs

Coral reefs are the epitome of species diversity, yet the number of described scleractinian coral species, the framework-builders of coral reefs, remains moderate by comparison. DNA sequencing studies are rapidly challenging this notion by exposing a wealth of undescribed diversity, but the evolutionary and ecological significance of this diversity remains largely unclear. Here, we present an annotated genome for one of the most ubiquitous corals in the Indo-Pacific (Pachyseris speciosa) and uncover, through a comprehensive genomic and phenotypic assessment, that it comprises morphologically indistinguishable but ecologically divergent lineages. Demographic modeling based on whole-genome resequencing indicated that morphological crypsis (across micro- and macromorphological traits) was due to ancient morphological stasis rather than recent divergence. Although the lineages occur sympatrically across shallow and mesophotic habitats, extensive genotyping using a rapid molecular assay revealed differentiation of their ecological distributions. Leveraging "common garden" conditions facilitated by the overlapping distributions, we assessed physiological and quantitative skeletal traits and demonstrated concurrent phenotypic differentiation. Lastly, spawning observations of genotyped colonies highlighted the potential role of temporal reproductive isolation in the limited admixture, with consistent genomic signatures in genes related to morphogenesis and reproduction. Overall, our findings demonstrate the presence of ecologically and phenotypically divergent coral species without substantial morphological differentiation and provide new leads into the potential mechanisms facilitating such divergence. More broadly, they indicate that our current taxonomic framework for reef-building corals may be scratching the surface of the ecologically relevant diversity on coral reefs, consequently limiting our ability to protect or restore this diversity effectively.

A Bayesian Implementation of the Multispecies Coalescent Model with Introgression for Phylogenomic Analysis

Recent analyses suggest that cross-species gene flow or introgression is common in nature, especially during species divergences. Genomic sequence data can be used to infer introgression events and to estimate the timing and intensity of introgression, providing an important means to advance our understanding of the role of gene flow in speciation. Here, we implement the multispecies-coalescent-with-introgression model, an extension of the multispecies-coalescent model to incorporate introgression, in our Bayesian Markov chain Monte Carlo program Bpp. The multispecies-coalescent-with-introgression model accommodates deep coalescence (or incomplete lineage sorting) and introgression and provides a natural framework for inference using genomic sequence data. Computer simulation confirms the good statistical properties of the method, although hundreds or thousands of loci are typically needed to estimate introgression probabilities reliably. Reanalysis of data sets from the purple cone spruce confirms the hypothesis of homoploid hybrid speciation. We estimated the introgression probability using the genomic sequence data from six mosquito species in the Anopheles gambiae species complex, which varies considerably across the genome, likely driven by differential selection against introgressed alleles.

Ecomorphological divergence and trophic resource partitioning in 15 syntopic Indo-Pacific parrotfishes (Labridae: Scarini)

Adaptive diversification is a product of both phylogenetic constraint and ecological opportunity. The species-rich parrotfish genera Scarus and Chlorurus display considerable variation in trophic cranial morphology, but these parrotfishes are often described as generalist herbivores. Recent work has suggested that parrotfish partition trophic resources at very fine spatial scales, raising the question of whether interspecific differences in cranial morphology reflect trophic partitioning. We tested this hypothesis by comparing targeted feeding substrata with a previously published dataset of nine cranial morphological traits. We sampled feeding substrata of 15 parrotfish species at Lizard Island, Great Barrier Reef, Australia, by following individuals until focused biting was observed, then extracting a bite core 22 mm in diameter. Three indices were parameterized for each bite core: substratum taphonomy, maximum turf height and cover of crustose coralline algae. Parrotfish species were spread along a single axis of variation in feeding substrata: successional status of the substratum taphonomy and epilithic and endolithic biota. This axis of trophic variation was significantly correlated with cranial morphology, indicating that morphological disparity within this clade is associated with interspecific partitioning of feeding substrata. Phylogenetic signal and phylomorphospace analyses revealed that the evolution of this clade involved a hitherto-unrecognized level of trophic diversification.

Adaptive venom evolution and toxicity in octopods is driven by extensive novel gene formation, expansion, and loss

BACKGROUND

Cephalopods represent a rich system for investigating the genetic basis underlying organismal novelties. This diverse group of specialized predators has evolved many adaptations including proteinaceous venom. Of particular interest is the blue-ringed octopus genus (Hapalochlaena), which are the only octopods known to store large quantities of the potent neurotoxin, tetrodotoxin, within their tissues and venom gland.

FINDINGS

To reveal genomic correlates of organismal novelties, we conducted a comparative study of 3 octopod genomes, including the Southern blue-ringed octopus (Hapalochlaena maculosa). We present the genome of this species and reveal highly dynamic evolutionary patterns at both non-coding and coding organizational levels. Gene family expansions previously reported in Octopus bimaculoides (e.g., zinc finger and cadherins, both associated with neural functions), as well as formation of novel gene families, dominate the genomic landscape in all octopods. Examination of tissue-specific genes in the posterior salivary gland revealed that expression was dominated by serine proteases in non-tetrodotoxin-bearing octopods, while this family was a minor component in H. maculosa. Moreover, voltage-gated sodium channels in H. maculosa contain a resistance mutation found in pufferfish and garter snakes, which is exclusive to the genus. Analysis of the posterior salivary gland microbiome revealed a diverse array of bacterial species, including genera that can produce tetrodotoxin, suggestive of a possible production source.

CONCLUSIONS

We present the first tetrodotoxin-bearing octopod genome H. maculosa, which displays lineage-specific adaptations to tetrodotoxin acquisition. This genome, along with other recently published cephalopod genomes, represents a valuable resource from which future work could advance our understanding of the evolution of genomic novelty in this family.

Genomic signatures in the coral holobiont reveal host adaptations driven by Holocene climate change and reef specific symbionts

Genetic signatures caused by demographic and adaptive processes during past climatic shifts can inform predictions of species' responses to anthropogenic climate change. To identify these signatures in Acropora tenuis, a reef-building coral threatened by global warming, we first assembled the genome from long reads and then used shallow whole-genome resequencing of 150 colonies from the central inshore Great Barrier Reef to inform population genomic analyses. We identify population structure in the host that reflects a Pleistocene split, whereas photosymbiont differences between reefs most likely reflect contemporary (Holocene) conditions. Signatures of selection in the host were associated with genes linked to diverse processes including osmotic regulation, skeletal development, and the establishment and maintenance of symbiosis. Our results suggest that adaptation to post-glacial climate change in A. tenuis has involved selection on many genes, while differences in symbiont specificity between reefs appear to be unrelated to host population structure.

AMPA receptor auxiliary subunits emerged during early vertebrate evolution by neo/subfunctionalization of unrelated proteins

In mammalian synapses, the function of ionotropic glutamate receptors is critically modulated by auxiliary subunits. Most of these specifically regulate the synaptic localization and electrophysiological properties of AMPA-type glutamate receptors (AMPARs). Here, we comprehensively investigated the animal evolution of the protein families that contain AMPAR auxiliary subunits (ARASs). We observed that, on average, vertebrates have four times more ARASs than other animal species. We also demonstrated that ARASs belong to four unrelated protein families: CACNG-GSG1, cornichon, shisa and Dispanin C. Our study demonstrates that, despite the ancient origin of these four protein families, the majority of ARASs emerged during vertebrate evolution by independent but convergent processes of neo/subfunctionalization that resulted in the multiple ARASs found in present vertebrate genomes. Importantly, although AMPARs appeared and diversified in the ancestor of bilateral animals, the ARAS expansion did not occur until much later, in early vertebrate evolution. We propose that the surge in ARASs and consequent increase in AMPAR functionalities, contributed to the increased complexity of vertebrate brains and cognitive functions.

Population genetics of the coral Acropora millepora: Toward genomic prediction of bleaching

Although reef-building corals are declining worldwide, responses to bleaching vary within and across species and are partly heritable. Toward predicting bleaching response from genomic data, we generated a chromosome-scale genome assembly for the coral Acropora millepora We obtained whole-genome sequences for 237 phenotyped samples collected at 12 reefs along the Great Barrier Reef, among which we inferred little population structure. Scanning the genome for evidence of local adaptation, we detected signatures of long-term balancing selection in the heat-shock co-chaperone sacsin We conducted a genome-wide association study of visual bleaching score for 213 samples, incorporating the polygenic score derived from it into a predictive model for bleaching in the wild. These results set the stage for genomics-based approaches in conservation strategies.

Defining Species When There is Gene Flow

Whatever one's definition of species, it is generally expected that individuals of the same species should be genetically more similar to each other than they are to individuals of another species. Here, we show that in the presence of cross-species gene flow, this expectation may be incorrect. We use the multispecies coalescent model with continuous-time migration or episodic introgression to study the impact of gene flow on genetic differences within and between species and highlight a surprising but plausible scenario in which different population sizes and asymmetrical migration rates cause a genetic sequence to be on average more closely related to a sequence from another species than to a sequence from the same species. Our results highlight the extraordinary impact that even a small amount of gene flow may have on the genetic history of the species. We suggest that contrasting long-term migration rate and short-term hybridization rate, both of which can be estimated using genetic data, may be a powerful approach to detecting the presence of reproductive barriers and to define species boundaries.[Gene flow; introgression; migration; multispecies coalescent; species concept; species delimitation.].

TREEasy: An automated workflow to infer gene trees, species trees, and phylogenetic networks from multilocus data

Multilocus genomic data sets can be used to infer a rich set of information about the evolutionary history of a lineage, including gene trees, species trees, and phylogenetic networks. However, user-friendly tools to run such integrated analyses are lacking, and workflows often require tedious reformatting and handling time to shepherd data through a series of individual programs. Here, we present a tool written in Python-TREEasy-that performs automated sequence alignment (with MAFFT), gene tree inference (with IQ-Tree), species inference from concatenated data (with IQ-Tree and RaxML-NG), species tree inference from gene trees (with ASTRAL, MP-EST, and STELLS2), and phylogenetic network inference (with SNaQ and PhyloNet). The tool only requires FASTA files and nine parameters as inputs. The tool can be run as command line or through a Graphical User Interface (GUI). As examples, we reproduced a recent analysis of staghorn coral evolution, and performed a new analysis on the evolution of the "WGD clade" of yeast. The latter revealed novel patterns that were not identified by previous analyses. TREEasy represents a reliable and simple tool to accelerate research in systematic biology (https://github.com/MaoYafei/TREEasy).

The Impact of Cross-Species Gene Flow on Species Tree Estimation

Recent analyses of genomic sequence data suggest cross-species gene flow is common in both plants and animals, posing challenges to species tree estimation. We examine the levels of gene flow needed to mislead species tree estimation with three species and either episodic introgressive hybridization or continuous migration between an outgroup and one ingroup species. Several species tree estimation methods are examined, including the majority-vote method based on the most common gene tree topology (with either the true or reconstructed gene trees used), the UPGMA method based on the average sequence distances (or average coalescent times) between species, and the full-likelihood method based on multilocus sequence data. Our results suggest that the majority-vote method based on gene tree topologies is more robust to gene flow than the UPGMA method based on coalescent times and both are more robust than likelihood assuming a multispecies coalescent (MSC) model with no cross-species gene flow. Comparison of the continuous migration model with the episodic introgression model suggests that a small amount of gene flow per generation can cause drastic changes to the genetic history of the species and mislead species tree methods, especially if the species diverged through radiative speciation events. Estimates of parameters under the MSC with gene flow suggest that African mosquito species in the Anopheles gambiae species complex constitute such an example of extreme impact of gene flow on species phylogeny. [IM; introgression; migration; MSci; multispecies coalescent; species tree.]

A Likely Ancient Genome Duplication in the Speciose Reef-Building Coral Genus, Acropora

Whole-genome duplication (WGD) has been recognized as a significant evolutionary force in the origin and diversification of multiple organisms. Acropora, a speciose reef-building coral genus, is suspected to have originated by polyploidy. Yet, there is no genetic evidence to support this hypothesis. Using comprehensive phylogenomic and comparative genomic approaches, we analyzed six Acroporid genomes and found that a WGD event likely occurred ∼31 million years ago in the most recent common ancestor of Acropora, concurrent with a worldwide coral extinction. We found that duplicated genes were highly enriched in gene regulation functions, including those of stress responses. The functional clusters of duplicated genes are related to the divergence of gene expression patterns during development. Some proteinaceous toxins were generated by WGD in Acropora compared with other cnidarian species. Collectively, this study provides evidence for an ancient WGD event in corals, which helps explain the origin and diversification of Acropora.

•

An ancient genome duplication occurred in the most recent common ancestor of Acropora

•

This WGD event likely occurred between 28 and 36 mya in Acropora

•

The WGD event potentially contributes to the origin and diversification of Acropora

•

Duplications of toxic proteins were found in Acropora following the WGD