Genome-wide SNP analysis explains coral diversity and recovery in the Ryukyu Archipelago

Chromosomal inversions harbour excess mutational load in the coral, Acropora kenti, on the Great Barrier Reef

The future survival of coral reefs in the Anthropocene depends on the capacity of corals to adapt as oceans warm and extreme weather events become more frequent. Targeted interventions designed to assist evolutionary processes in corals require a comprehensive understanding of the distribution and structure of standing variation, however, efforts to map genomic variation in corals have so far focussed almost exclusively on SNPs, overlooking structural variants that have been shown to drive adaptive processes in other taxa. Here, we show that the reef-building coral, Acropora kenti, harbours at least five large, highly polymorphic structural variants, all of which exhibit signatures of strongly suppressed recombination in heterokaryotypes, a feature commonly associated with chromosomal inversions. Based on their high minor allele frequency, uniform distribution across habitats and elevated genetic load, we propose that these inversions in A. kenti are likely to be under balancing selection. An excess of SNPs with high impact on protein-coding genes within these loci elevates their importance both as potential targets for adaptive selection and as contributors to genetic decline if coral populations become fragmented or inbred in future.

Genome and tissue-specific transcriptomes of the large-polyp coral, Fimbriaphyllia (Euphyllia) ancora: a recipe for a coral polyp

Coral polyps are composed of four tissues; however, their characteristics are largely unexplored. Here we report biological characteristics of tentacles (Te), mesenterial filaments (Me), body wall (Bo), and mouth with pharynx (MP), using comparative genomic, morpho-histological, and transcriptomic analyses of the large-polyp coral, Fimbriaphyllia ancora. A draft F. ancora genome assembly of 434 Mbp was created. Morpho-histological and transcriptomic characterization of the four tissues showed that they have distinct differences in structure, primary cellular composition, and transcriptional profiles. Tissue-specific, highly expressed genes (HEGs) of Te are related to biological defense, predation, and coral-algal symbiosis. Me expresses multiple digestive enzymes, whereas Bo expresses innate immunity and biomineralization-related molecules. Many receptors for neuropeptides and neurotransmitters are expressed in MP. This dataset and new insights into tissue functions will facilitate a deeper understanding of symbiotic biology, immunology, biomineralization, digestive biology, and neurobiology in corals.

Genome-wide analysis to uncover how Pocillopora acuta survives the challenging intertidal environment

Characterisation of genomic variation among corals can help uncover variants underlying trait differences and contribute towards genotype prioritisation in coastal restoration projects. For example, there is growing interest in identifying resilient genotypes for transplantation, and to better understand the genetic processes that allow some individuals to survive in specific conditions better than others. The coral species Pocillopora acuta is known to survive in a wide range of habitats, from reefs artificial coastal defences, suggesting its potential use as a starter species for ecological engineering efforts involving coral transplantation onto intertidal seawalls. However, the intertidal section of coastal armour is a challenging environment for corals, with conditions during periods of emersion being particularly stressful. Here, we scanned the entire genome of P. acuta corals to identify the regions harbouring single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) that separate intertidal colonies (n = 18) from those found in subtidal areas (n = 21). Findings revealed 74,391 high quality SNPs distributed across 386 regions of the P. acuta genome. While the majority of the detected SNPs were in non-coding regions, 12% were identified in exons (i.e. coding regions). Functional SNPs that were significantly associated with intertidal colonies were found in overrepresented genomic regions linked to cellular homeostasis, metabolism, and signalling processes, which may represent local environmental adaptation in the intertidal. Interestingly, regions that exhibited CNVs were also associated with metabolic and signalling processes, suggesting P. acuta corals living in the intertidal have a high capacity to perform biological functions critical for survival in extreme environments.

Multiomics data integration, limitations, and prospects to reveal the metabolic activity of the coral holobiont

Since their radiation in the Middle Triassic period ∼240 million years ago, stony corals have survived past climate fluctuations and five mass extinctions. Their long-term survival underscores the inherent resilience of corals, particularly when considering the nutrient-poor marine environments in which they have thrived. However, coral bleaching has emerged as a global threat to coral survival, requiring rapid advancements in coral research to understand holobiont stress responses and allow for interventions before extensive bleaching occurs. This review encompasses the potential, as well as the limits, of multiomics data applications when applied to the coral holobiont. Synopses for how different omics tools have been applied to date and their current restrictions are discussed, in addition to ways these restrictions may be overcome, such as recruiting new technology to studies, utilizing novel bioinformatics approaches, and generally integrating omics data. Lastly, this review presents considerations for the design of holobiont multiomics studies to support lab-to-field advancements of coral stress marker monitoring systems. Although much of the bleaching mechanism has eluded investigation to date, multiomic studies have already produced key findings regarding the holobiont's stress response, and have the potential to advance the field further.

Genomic Data Reveal Diverse Biological Characteristics of Scleractinian Corals and Promote Effective Coral Reef Conservation

Reef-building corals (Scleractinia, Anthozoa, Cnidaria) are the keystone organisms of coral reefs, which constitute the most diverse marine ecosystems. Since the first decoded coral genome reported in 2011, about 40 reference genomes are registered as of 2023. Comparative genomic analyses of coral genomes have revealed genomic characters that may underlie unique biological characteristics and coral diversification. These include existence of genes for biosynthesis of mycosporine-like amino acids, loss of an enzyme necessary for cysteine biosynthesis in family Acroporidae, and lineage-specific gene expansions of DMSP lyase-like genes in the genus Acropora. While symbiosis with endosymbiotic photosynthetic dinoflagellates is a common biological feature among reef-building corals, genes associated with the intricate symbiotic relationship encompass not only those shared by many coral species, but also genes that were uniquely duplicated in each coral lineage, suggesting diversified molecular mechanisms of coral-algal symbiosis. Coral genomic data have also enabled detection of hidden, complex population structures of corals, indicating the need for species-specific, local-scale, carefully considered conservation policies for effective maintenance of corals. Consequently, accumulating coral genomic data from a wide range of taxa and from individuals of a species not only promotes deeper understanding of coral reef biodiversity, but also promotes appropriate and effective coral reef conservation. Considering the diverse biological traits of different coral species and accurately understanding population structure and genetic diversity revealed by coral genomic analyses during coral reef restoration planning could enable us to "archive" coral reef environments that are nearly identical to natural coral reefs.

Coral larvae suppress heat stress response during the onset of symbiosis decreasing their odds of survival

The endosymbiosis between most corals and their photosynthetic dinoflagellate partners begins early in the host life history, when corals are larvae or juvenile polyps. The capacity of coral larvae to buffer climate-induced stress while in the process of symbiont acquisition could come with physiological trade-offs that alter behaviour, development, settlement and survivorship. Here we examined the joint effects of thermal stress and symbiosis onset on colonization dynamics, survival, metamorphosis and host gene expression of Acropora digitifera larvae. We found that thermal stress decreased symbiont colonization of hosts by 50% and symbiont density by 98.5% over 2 weeks. Temperature and colonization also influenced larval survival and metamorphosis in an additive manner, where colonized larvae fared worse or prematurely metamorphosed more often than noncolonized larvae under thermal stress. Transcriptomic responses to colonization and thermal stress treatments were largely independent, while the interaction of these treatments revealed contrasting expression profiles of genes that function in the stress response, immunity, inflammation and cell cycle regulation. The combined treatment either cancelled or lowered the magnitude of expression of heat-stress responsive genes in the presence of symbionts, revealing a physiological cost to acquiring symbionts at the larval stage with elevated temperatures. In addition, host immune suppression, a hallmark of symbiosis onset under ambient temperature, turned to immune activation under heat stress. Thus, by integrating the physical environment and biotic pressures that mediate presettlement event in corals, our results suggest that colonization may hinder larval survival and recruitment under projected climate scenarios.

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.

Genomic analysis of a reef-building coral, Acropora digitifera, reveals complex population structure and a migration network in the Nansei Islands, Japan

Understanding the structure and connectivity of coral populations is fundamental for developing marine conservation policies, especially in patchy environments such as archipelagos. The Nansei Islands, extending more than 1000 km in southwestern Japan, are characterized by high levels of biodiversity and endemism, supported by coral reefs, which make this region ideal for assessing genetic attributes of coral populations. In this study, we conducted population genomic analyses based on genome-wide, single-nucleotide polymorphisms (SNPs) of Acropora digitifera, a common species in the Nansei Islands. By merging newly obtained genome resequencing data with previously published data, we identified more than 4 million genome-wide SNPs in 303 colonies collected at 22 locations, with sequencing coverage ranging from 3.91× to 27.41×. While population structure analyses revealed genetic similarities between the southernmost and northernmost locations, separated by >1000 km, several subpopulations in intermediate locations suggested limited genetic admixture, indicating conflicting migration tendencies in the Nansei Islands. Although migration networks revealed a general tendency of northward migration along the Kuroshio Current, a substantial amount of southward migration was also detected, indicating important contributions of minor ocean currents to coral larval dispersal. Moreover, heterogeneity in the transition of effective population sizes among locations suggests different histories for individual subpopulations. The unexpected complexity of both past and present population dynamics in the Nansei Islands implies that heterogeneity of ocean currents and local environments, past and present, have influenced the population structure of this species, and similar unexpected population complexities may be expected for other marine species with similar reproductive modes.

Population genomics for symbiotic anthozoans: can reduced representation approaches be used for taxa without reference genomes?

Population genetic studies of symbiotic anthozoans have been historically challenging because their endosymbioses with dinoflagellates have impeded marker development. Genomic approaches like reduced representation sequencing alleviate marker development issues but produce anonymous loci, and without a reference genome, it is unknown which organism is contributing to the observed patterns. Alternative methods such as bait-capture sequencing targeting Ultra-Conserved Elements are now possible but costly. Thus, RADseq remains attractive, but how useful are these methods for symbiotic anthozoan taxa without a reference genome to separate anthozoan from algal sequences? We explore this through a case-study using a double-digest RADseq dataset for the sea anemone Bartholomea annulata. We assembled a holobiont dataset (3854 loci) for 101 individuals, then used a reference genome to create an aposymbiotic dataset (1402 loci). For both datasets, we investigated population structure and used coalescent simulations to estimate demography and population parameters. We demonstrate complete overlap in the spatial patterns of genetic diversity, demographic histories, and population parameter estimates for holobiont and aposymbiotic datasets. We hypothesize that the unique combination of anthozoan biology, diversity of the endosymbionts, and the manner in which assembly programs identify orthologous loci alleviates the need for reference genomes in some circumstances. We explore this hypothesis by assembling an additional 21 datasets using the assembly programs pyRAD and Stacks. We conclude that RADseq methods are more tractable for symbiotic anthozoans without reference genomes than previously realized.

Signatures of selection underpinning rapid coral adaptation to the world's warmest reefs

Coral populations in the world’s warmest reefs, the Persian/Arabian Gulf (PAG), represent an ideal model system to understand the evolutionary response of coral populations to past and present environmental change and to identify genomic loci that contribute to elevated thermal tolerance. Here, we use population genomics of the brain coral Platygyra daedalea to show that corals in the PAG represent a distinct subpopulation that was established during the Holocene marine transgression, and identify selective sweeps in their genomes associated with thermal adaptation. We demonstrate the presence of positive and disruptive selection and provide evidence for selection of differentially methylated haplotypes. While demographic analyses suggest limited potential for genetic rescue of neighboring Indian Ocean reefs, the presence of putative targets of selection in corals outside of the PAG offers hope that loci associated with thermal tolerance may be present in the standing genetic variation.

Two divergent haplogroups of a sacsin-like gene in Acropora corals

Reef-building corals are declining due to environmental changes. Sacsin is a member of the heat shock proteins and has been reported as a candidate protein associated with the stress response in Acropora corals. Recently, high nucleotide diversity and the persistence of two divergent haplogroups of sacsin-like genes in Acropora millepora have been reported. While it was not clear when the two haplogroups have split and whether the haplogroups have persisted in only A. millepora or the other lineages in the genus Acropora. In this study, we analyzed a genomic region containing a sacsin-like gene from Acropora and Montipora species. Higher nucleotide diversity in the sacsin-like gene compared with that of surrounding regions was also observed in A. digitifera. This nucleotide diversity is derived from two divergent haplogroups of a sacsin-like gene, which are present in at least three Acropora species. The origin of these two haplogroups can be traced back before the divergence of Acropora and Montipora (119 Ma). Although the link between exceptionally high genetic variation in sacsin-like genes and functional differences in sacsin-like proteins is not clear, the divergent haplogroups may respond differently to envionmental stressors and serve in the adaptive phsiological ecology of these keystone species.

Barriers and corridors of gene flow in an urbanized tropical reef system

Information about the distribution of alleles among marine populations is critical for determining patterns of genetic connectivity that are essential in modern conservation planning. To estimate population connectivity in Singapore's urbanized equatorial reef system, we analysed single nucleotide polymorphisms (SNPs) from two species of reef-building corals with distinct life histories. For Porites sp., a broadcast-spawning coral, we found cryptic lineages that were differentially distributed at inshore and central-offshore sites that could be attributed to contemporary surface current regimes. Near panmixia was observed for Pocillopora acuta with differentiation of colonies at the farthest site from mainland Singapore, a possible consequence of the brooding nature and relatively long pelagic larval duration of the species. Furthermore, analysis of recent gene flow showed that 60-80% of colonies in each population were nonmigrants, underscoring self-recruitment as an important demographic process in this reef system. Apart from helping to enhance the management of Singapore's coral reef ecosystems, findings here pave the way for better understanding of the evolution of marine populations in South-East Asia.

Higher Genetic Diversity of the Common Sea Cucumber Holothuria (Halodeima) atra in Marine Protected Areas of the Central and Southern Ryukyu Islands

Recently, sea cucumbers (Echinodermata: Holothuroidea) have been over-exploited in many areas of the world, including in the Ryukyu Islands, southern Japan, due to increases in their economic importance. Nevertheless, management and protection of sea cucumbers are insufficient worldwide. The black sea cucumber Holothuria (Halodeima) atra Jaeger, 1833, inhabits a large range across the Indo-West Pacific Ocean and is a widely harvested species. Here we conducted population genetic analyses on H. atra using partial mitochondrial DNA sequences of cytochrome c oxidase subunit I (COI) and 16S ribosomal RNA (16S) to examine 11 different populations around three island groups in the middle Ryukyus; Okinawajima Island, the Kerama Islands, and the Sakishima Islands, all within Okinawa Prefecture. We found 27 haplotypes for COI and 16 haplotypes for 16S. Locations within national and quasi-national parks (Zamami Island, Keramas, and Manza, Okinawajima; managed by the national Ministry of Environment and Okinawa Prefecture, respectively) had the highest number of haplotypes, whereas locations with less management and more anthropogenic pressure had lower numbers The mean of all samples' genetic diversity indices was moderate with regards to both haplotype and nucleotide diversity. According to our results, Zamami Ama was the most genetically diverse location based on both markers used, likely because it is located within Kerama-Shoto National Park with comparatively stricter regulations than most other locations. Based on our COI sequences, three-quarters of the locations with the highest haplotype diversity were found to be distant from Okinawajima Island, indicating that the genetic diversity of H. atra was reduced around Okinawajima Island. Our results possibly reflect negative impacts from anthropogenic pressures such as over-harvesting and coastal development, although future comprehensive research including sequences of nuclear loci is needed to confirm this hypothesis.

Correlation between Organelle Genetic Variation and RNA Editing in Dinoflagellates Associated with the Coral Acropora digitifera

In order to develop successful strategies for coral reef preservation, it is critical that the biology of both host corals and symbiotic algae are investigated. In the Ryukyu Archipelago, which encompasses many islands spread over ∼500 km of the Pacific Ocean, four major populations of the coral Acropora digitifera have been studied using whole-genome shotgun (WGS) sequence analysis (Shinzato C, Mungpakdee S, Arakaki N, Satoh N. 2015. Genome-wide single-nucleotide polymorphism (SNP) analysis explains coral diversity and recovery in the Ryukyu Archipelago. Sci Rep. 5:18211.). In contrast, the diversity of the symbiotic dinoflagellates associated with these A. digitifera populations is unknown. It is therefore unclear if these two core components of the coral holobiont share a common evolutionary history. This issue can be addressed for the symbiotic algal populations by studying the organelle genomes of their mitochondria and plastids. Here, we analyzed WGS data from ∼150 adult A. digitifera, and by mapping reads to the available reference genome sequences, we extracted 2,250 sequences representing 15 organelle genes of Symbiodiniaceae. Molecular phylogenetic analyses of these mitochondrial and plastid gene sets revealed that A. digitifera from the southern Yaeyama islands harbor a different Symbiodiniaceae population than the islands of Okinawa and Kerama in the north, indicating that the distribution of symbiont populations partially matches that of the four host populations. Interestingly, we found that numerous SNPs correspond to known RNA-edited sites in 14 of the Symbiodiniaceae organelle genes, with mitochondrial genes showing a stronger correspondence than plastid genes. These results suggest a possible correlation between RNA editing and SNPs in the two organelle genomes of symbiotic dinoflagellates.

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.

De novo transcriptome assembly from the gonads of a scleractinian coral, Euphyllia ancora: molecular mechanisms underlying scleractinian gametogenesis

BACKGROUND

Sexual reproduction of scleractinians has captured the attention of researchers and the general public for decades. Although extensive ecological data has been acquired, underlying molecular and cellular mechanisms remain largely unknown. In this study, to better understand mechanisms underlying gametogenesis, we isolated ovaries and testes at different developmental phases from a gonochoric coral, Euphyllia ancora, and adopted a transcriptomic approach to reveal sex- and phase-specific gene expression profiles. In particular, we explored genes associated with oocyte development and maturation, spermiogenesis, sperm motility / capacitation, and fertilization.

RESULTS

1.6 billion raw reads were obtained from 24 gonadal samples. De novo assembly of trimmed reads, and elimination of contigs derived from symbiotic dinoflagellates (Symbiodiniaceae) and other organisms yielded a reference E. ancora gonadal transcriptome of 35,802 contigs. Analysis of 4 developmental phases identified 2023 genes that were differentially expressed during oogenesis and 678 during spermatogenesis. In premature/mature ovaries, 631 genes were specifically upregulated, with 538 in mature testes. Upregulated genes included those involved in gametogenesis, gamete maturation, sperm motility / capacitation, and fertilization in other metazoans, including humans. Meanwhile, a large number of genes without homology to sequences in the SWISS-PROT database were also observed among upregulated genes in premature / mature ovaries and mature testes.

CONCLUSIONS

Our findings show that scleractinian gametogenesis shares many molecular characteristics with that of other metazoans, but it also possesses unique characteristics developed during cnidarian and/or scleractinian evolution. To the best of our knowledge, this study is the first to create a gonadal transcriptome assembly from any scleractinian. This study and associated datasets provide a foundation for future studies regarding gametogenesis and differences between male and female colonies from molecular and cellular perspectives. Furthermore, our transcriptome assembly will be a useful reference for future development of sex-specific and/or stage-specific germ cell markers that can be used in coral aquaculture and ecological studies.

Seascape genomics as a new tool to empower coral reef conservation strategies: An example on north-western Pacific Acropora digitifera

Coral reefs are suffering a major decline due to the environmental constraints imposed by climate change. Over the last 20 years, three major coral bleaching events occurred in concomitance with anomalous heatwaves, provoking a severe loss of coral cover worldwide. The conservation strategies for preserving reefs, as they are implemented now, cannot cope with global climatic shifts. Consequently, researchers are advocating for preservation networks to be set-up to reinforce coral adaptive potential. However, the main obstacle to this implementation is that studies on coral adaption are usually hard to generalize at the scale of a reef system. Here, we study the relationships between genotype frequencies and environmental characteristics of the sea (seascape genomics), in combination with connectivity analysis, to investigate the adaptive potential of a flagship coral species of the Ryukyu Archipelago (Japan). By associating genotype frequencies with descriptors of historical environmental conditions, we discovered six genomic regions hosting polymorphisms that might promote resistance against heat stress. Remarkably, annotations of genes in these regions were consistent with molecular roles associated with heat responses. Furthermore, we combined information on genetic and spatial distances between reefs to predict connectivity at a regional scale. The combination of these results portrayed the adaptive potential of this population: we were able to identify reefs carrying potential heat stress adapted genotypes and to understand how they disperse to neighbouring reefs. This information was summarized by objective, quantifiable and mappable indices covering the whole region, which can be extremely useful for future prioritization of reefs in conservation planning. This framework is transferable to any coral species on any reef system and therefore represents a valuable tool for empowering preservation efforts dedicated to the protection of coral reefs in warming oceans.

An Investigation into the Genetic History of Japanese Populations of Three Starfish, Acanthaster planci, Linckia laevigata, and Asterias amurensis, Based on Complete Mitochondrial DNA Sequences

Crown-of-thorns starfish, Acanthaster planci (COTS), are common in coral reefs of Indo-Pacific Ocean. Since they are highly fecund predators of corals, periodic outbreaks of COTS cause substantial loss of healthy coral reefs. Using complete mitochondrial DNA sequences, we here examined how COTS outbreaks in the Ryukyu Archipelago, Japan are reflected by the profile of their population genetics. Population genetics of the blue starfish, Linckia laevigata, which lives in the Ryukyu Archipelago, but not break out and the northern Pacific sea star, Asterias amurensis, which lives in colder seawater around the main Islands of Japan, were also examined as controls. Our results showed that As. amurensis has at least two local populations that diverged approximately 4.7 million years ago (MYA), and no genetic exchanges have occurred between the populations since then. Linckia laevigata shows two major populations in the Ryukyu Archipelago that likely diverged ∼6.8 MYA. The two populations, each comprised of individuals collected from coast of the Okinawa Island and those from the Ishigaki Island, suggest the presence of two cryptic species in the Ryukyu Archipelago. On the other hand, population genetics of COTS showed a profile quite different from those of Asterias and Linckia. At least five lineages of COTS have arisen since their divergence ∼0.7 MYA, and each of the lineages is present at the Okinawa Island, Miyako Island, and Ishigaki Island. These results suggest that COTS have experienced repeated genetic bottlenecks that may be associated with or caused by repeated outbreaks.

Parentage analyses identify local dispersal events and sibling aggregations in a natural population of Millepora hydrocorals, a free-spawning marine invertebrate

Dispersal is a critical process for the persistence and productivity of marine populations. For many reef species, there is increasing evidence that local demography and self-recruitment have major consequences on their genetic diversity and adaptation to environmental change. Yet empirical data of dispersal patterns in reef-building species remain scarce. Here, we document the first genetic estimates of self-recruitment and dispersal distances in a free-spawning marine invertebrate, the hydrocoral Millepora cf. platyphylla. Using twelve microsatellite markers, we gathered genotypic information from 3,160 georeferenced colonies collected over 27,000 m² of a single reef in three adjacent habitats in Moorea, French Polynesia; the mid slope, upper slope, and back reef. Although the adult population was predominantly clonal (85% were clones), our parentage analysis revealed a moderate self-recruitment rate with a minimum of 8% of sexual propagules produced locally. Assigned offspring often settled at <10 m from their parents and dispersal events decrease with increasing geographic distance. There were no discrepancies between the dispersal distances of offspring assigned to parents belonging to clonal versus nonclonal genotypes. Interhabitat dispersal events via cross-reef transport were also detected for sexual and asexual propagules. Sibship analysis showed that full siblings recruit nearby on the reef (more than 40% settled at <30 m), resulting in sibling aggregations. Our findings highlight the importance of self-recruitment together with clonality in stabilizing population dynamics, which may ultimately enhance local sustainability and resilience to disturbance.

Small-scale population genetic structure of the sand bubbler crab Scopimera ryukyuensis in the Ryukyu Islands, Japan

Generally, the gene flow of marine organisms is well maintained, but some local populations of coastal species are genetically differentiated even on a small scale (genetic patchiness). Small-scale isolation can be crucial for understanding genetic diversity within a species. The present study examined the population genetic structure of the sand bubbler crab Scopimera ryukyuensis, which is endemic to the Ryukyu Islands in the northwestern Pacific. A total of 52 haplotypes of mitochondrial cytochrome c oxidase subunit I were recovered from 197 specimens collected from four islands. The haplotype and nucleotide diversities were relatively high in the central Ryukyus (Amami-Oshima and Okinawa Islands) with some exceptions but were low at the southern edge of the geographical distribution of the species, i.e., the southern Ryukyus (Ishigaki and Iriomote Islands). Pairwise F_ST analysis suggested that the gene flow of S. ryukyuensis was largely restricted. The local populations of the species are differentiated among islands, except for stations on Ishigaki Island and a station on Iriomote Island. Moreover, a clear intra-island population genetic structure was observed within Amami-Oshima and Iriomote Islands, e.g., only 20 km between stations. Small-scale isolation among local populations may be a common tendency for coastal species in the Ryukyu Islands, considering the results of previous studies on corals.

Presence-absence polymorphisms of single-copy genes in the stony coral Acropora digitifera

Background

Despite the importance of characterizing genetic variation among coral individuals for understanding phenotypic variation, the correlation between coral genomic diversity and phenotypic expression is still poorly understood.

Results

In this study, we detected a high frequency of genes showing presence–absence polymorphisms (PAPs) for single-copy genes in Acropora digitifera. Among 10,455 single-copy genes, 516 (5%) exhibited PAPs, including 32 transposable element (TE)-related genes. Five hundred sixteen genes exhibited a homozygous absence in one (102) or more than one (414) individuals (n = 33), indicating that most of the absent alleles were not rare variants. Among genes showing PAPs (PAP genes), roughly half were expressed in adults and/or larvae, and the PAP status was associated with differential expression among individuals. Although 85% of PAP genes were uncharacterized or had ambiguous annotations, 70% of these genes were specifically distributed in cnidarian lineages in eumetazoa, suggesting that these genes have functional roles related to traits related to cnidarians or the family Acroporidae or the genus Acropora. Indeed, four of these genes encoded toxins that are usually components of venom in cnidarian-specific cnidocytes. At least 17% of A. digitifera PAP genes were also PAPs in A. tenuis, the basal lineage in the genus Acropora, indicating that PAPs were shared among species in Acropora.

Conclusions

Expression differences caused by a high frequency of PAP genes may be a novel genomic feature in the genus Acropora; these findings will contribute to improve our understanding of correlation between genetic and phenotypic variation in corals.

The active spread of adaptive variation for reef resilience

The speed at which species adapt depends partly on the rates of beneficial adaptation generation and how quickly they spread within and among populations. Natural rates of adaptation of corals may not be able to keep pace with climate warming. Several interventions have been proposed to fast-track thermal adaptation, including the intentional translocation of warm-adapted adults or their offspring (assisted gene flow, AGF) and the ex situ crossing of warm-adapted corals with conspecifics from cooler reefs (hybridization or selective breeding) and field deployment of those offspring. The introgression of temperature tolerance loci into the genomic background of cooler-environment corals aims to facilitate adaptation to warming while maintaining fitness under local conditions. Here we use research on selective sweeps and connectivity to understand the spread of adaptive variants as it applies to AGF on the Great Barrier Reef (GBR), focusing on the genus Acropora. Using larval biophysical dispersal modeling, we estimate levels of natural connectivity in warm-adapted northern corals. We then model the spread of adaptive variants from single and multiple reefs and assess if the natural and assisted spread of adaptive variants will occur fast enough to prepare receiving central and southern populations given current rates of warming. We also estimate fixation rates and spatial extent of fixation under multiple release scenarios to inform intervention design. Our results suggest that thermal tolerance is unlikely to spread beyond northern reefs to the central and southern GBR without intervention, and if it does, 30+ generations are needed for adaptive gene variants to reach fixation even under multiple release scenarios. We argue that if translocation, breeding, and reseeding risks are managed, AGF using multiple release reefs can be beneficial for the restoration of coral populations. These interventions should be considered in addition to conventional management and accompanied by strong mitigation of CO2 emissions.

Population Connectivity in the Common Reef Zoantharian Zoanthus sansibaricus (Anthozoa: Hexacorallia) in Southern Japan

Tropical and subtropical shallow benthic marine communities are highly diverse and balanced systems that constitute an important natural resource. Knowledge of the genetic diversity, connectivity and reproduction mode of each population is critical to understanding the fate of whole assemblages in times of disturbances. Importantly, the capability of populations to adapt to environmental challenges will be crucial to determining their survival. Here, we report on the population structure of the common reef zoantharian Zoanthus sansibaricus in the northwestern Pacific, by examining populations at three different locations in southern Japan using five highly variable microsatellite markers. Analyses of a population at the species' northern distribution limit combined with analyses of two subtropical populations suggest that habitat characteristics and ocean currents influence the connectivity and genetic diversity of this species. Our findings emphasize the adaptive ability of Z. sansibaricus to different environmental conditions and may help explain the wide distribution and generalist nature of this species.

Genome analysis of the rice coral Montipora capitata

Corals comprise a biomineralizing cnidarian, dinoflagellate algal symbionts, and associated microbiome of prokaryotes and viruses. Ongoing efforts to conserve coral reefs by identifying the major stress response pathways and thereby laying the foundation to select resistant genotypes rely on a robust genomic foundation. Here we generated and analyzed a high quality long-read based ~886 Mbp nuclear genome assembly and transcriptome data from the dominant rice coral, Montipora capitata from Hawai’i. Our work provides insights into the architecture of coral genomes and shows how they differ in size and gene inventory, putatively due to population size variation. We describe a recent example of foreign gene acquisition via a bacterial gene transfer agent and illustrate the major pathways of stress response that can be used to predict regulatory components of the transcriptional networks in M. capitata. These genomic resources provide insights into the adaptive potential of these sessile, long-lived species in both natural and human influenced environments and facilitate functional and population genomic studies aimed at Hawaiian reef restoration and conservation.

Binning enables efficient host genome reconstruction in cnidarian holobionts

Background

Many cnidarians, including stony corals, engage in complex symbiotic associations, comprising the eukaryotic host, photosynthetic algae, and highly diverse microbial communities—together referred to as holobiont. This taxonomic complexity makes sequencing and assembling coral host genomes extremely challenging. Therefore, previous cnidarian genomic projects were based on symbiont-free tissue samples. However, this approach may not be applicable to the majority of cnidarian species for ecological reasons. We therefore evaluated the performance of an alternative method based on sequence binning for reconstructing the genome of the stony coral Porites rus from a hologenomic sample and compared it to traditional approaches.

Results

Our results demonstrate that binning performs well for hologenomic data, producing sufficient reads for assembling the draft genome of P. rus. An assembly evaluation based on operational criteria showed results that were comparable to symbiont-free approaches in terms of completeness and usefulness, despite a high degree of fragmentation in our assembly. In addition, we found that binning provides sufficient data for exploratory k-mer estimation of genomic features, such as genome size and heterozygosity.

Conclusions

Binning constitutes a powerful approach for disentangling taxonomically complex coral hologenomes. Considering the recent decline of coral reefs on the one hand and previous limitations to coral genome sequencing on the other hand, binning may facilitate rapid and reliable genome assembly. This study also provides an important milestone in advancing binning from the metagenomic to the hologenomic and from the prokaryotic to the eukaryotic level.

Presence-Absence Polymorphisms of Highly Expressed FP Sequences Contribute to Fluorescent Polymorphisms in Acropora digitifera

Despite many hypotheses regarding the roles of fluorescent proteins (FPs), their biological roles and the genetic basis of FP-mediated color polymorphisms in Acropora remain unclear. In this study, we determined the genetic mechanism underlying fluorescent polymorphisms in A. digitifera. Using a high-throughput sequencing approach, we found that FP gene sequences in FP multigene family exhibit presence-absence polymorphism among individuals. A few particular sequences in short-to-middle wavelength emission and middle-to-long wavelength emission clades were highly expressed in adults, and different sequences were highly expressed in larvae. These highly expressed sequences were absent in the genomes of individuals with low total FP gene expression. In adults, presence-absence differences of the highly expressed FP sequences were consistent with measurements of emission spectra of corals, suggesting that presence-absence polymorphisms of these FP sequences contributed to the fluorescent polymorphisms. The functions of recombinant FPs encoded by highly expressed sequences in adult and larval stages were different, suggesting that expression of FP sequences with different functions may depend on the life-stage of A. digitifera. Highly expressed FP sequences exhibited presence-absence polymorphisms in subpopulations of A. digitifera, suggesting that presence-absence status is maintained during the evolution of A. digitifera subpopulations. The difference in FPs between adults and larvae and the polymorphisms of highly expressed FP genes may provide key insight into the biological roles of FPs in corals.

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

Reef-building corals provide the structural basis for one of Earth's most spectacular and diverse-but increasingly threatened-ecosystems. Modern Indo-Pacific reefs are dominated by species of the staghorn coral genus Acropora, but the evolutionary and ecological factors associated with their diversification and rise to dominance are unclear. Recent work on evolutionary radiations has demonstrated the importance of introgression and ecological opportunity in promoting diversification and ecological success. Here, we analyze the genomes of five staghorn coral species to examine the roles of introgression and ecological opportunity in the rise to dominance of Acropora. We found evidence for a history marked by a major introgression event as well as recurrent gene flow across species. In addition, we found that genes with topologies mismatching the species tree are evolving faster, which is suggestive of a role for introgression in spreading adaptive genetic variation. Demographic analysis showed that Acropora lineages profited from climate-driven mass extinctions in the Plio-Pleistocene, indicating that Acropora exploited ecological opportunity opened by a new climatic regime favoring species that could cope with rapid sea-level changes. Collectively, the genomes of reef-building corals have recorded an evolutionary history shaped by introgression and climate change, suggesting that Acropora-among most vulnerable corals to stressors-may be critical for understanding how reefs track the impending rapid sea-level changes of the Anthropocene.

Molecular signatures of host specificity linked to habitat specialization in Exaiptasia sea anemones

Rising ocean temperatures associated with global climate change induce breakdown of the symbiosis between coelenterates and photosynthetic microalgae of the genus Symbiodinium. Association with more thermotolerant partners could contribute to resilience, but the genetic mechanisms controlling specificity of hosts for particular Symbiodinium types are poorly known. Here, we characterize wild populations of a sea anemone laboratory model system for anthozoan symbiosis, from contrasting environments in Caribbean Panama. Patterns of anemone abundance and symbiont diversity were consistent with specialization of holobionts for particular habitats, with Exaiptasia pallida/S. minutum (ITS2 type B1) abundant on vertical substrate in thermally stable, shaded environments but E. brasiliensis/Symbiodinium sp. (ITS2 clade A) more common in shallow areas subject to high temperature and irradiance. Population genomic sequencing revealed a novel E. pallida population from the Bocas del Toro Archipelago that only harbors S. minutum. Loci most strongly associated with divergence of the Bocas-specific population were enriched in genes with putative roles in cnidarian symbiosis, including activators of the complement pathway of the innate immune system, thrombospondin-type-1 repeat domain proteins, and coordinators of endocytic recycling. Our findings underscore the importance of unmasking cryptic diversity in natural populations and the role of long-term evolutionary history in mediating interactions with Symbiodinium.

Genome-wide survey of single-nucleotide polymorphisms reveals fine-scale population structure and signs of selection in the threatened Caribbean elkhorn coral, Acropora palmata

The advent of next-generation sequencing tools has made it possible to conduct fine-scale surveys of population differentiation and genome-wide scans for signatures of selection in non-model organisms. Such surveys are of particular importance in sharply declining coral species, since knowledge of population boundaries and signs of local adaptation can inform restoration and conservation efforts. Here, we use genome-wide surveys of single-nucleotide polymorphisms in the threatened Caribbean elkhorn coral, Acropora palmata, to reveal fine-scale population structure and infer the major barrier to gene flow that separates the eastern and western Caribbean populations between the Bahamas and Puerto Rico. The exact location of this break had been subject to discussion because two previous studies based on microsatellite data had come to differing conclusions. We investigate this contradiction by analyzing an extended set of 11 microsatellite markers including the five previously employed and discovered that one of the original microsatellite loci is apparently under selection. Exclusion of this locus reconciles the results from the SNP and the microsatellite datasets. Scans for outlier loci in the SNP data detected 13 candidate loci under positive selection, however there was no correlation between available environmental parameters and genetic distance. Together, these results suggest that reef restoration efforts should use local sources and utilize existing functional variation among geographic regions in ex situ crossing experiments to improve stress resistance of this species.

A preliminary survey of zoantharian endosymbionts shows high genetic variation over small geographic scales on Okinawa-jima Island, Japan

Symbiotic dinoflagellates (genus Symbiodinium) shape the responses of their host reef organisms to environmental variability and climate change. To date, the biogeography of Symbiodinium has been investigated primarily through phylogenetic analyses of the ribosomal internal transcribed spacer 2 region. Although the marker can approximate species-level diversity, recent work has demonstrated that faster-evolving genes can resolve otherwise hidden species and population lineages, and that this diversity is often distributed over much finer geographical and environmental scales than previously recognized. Here, we use the noncoding region of the chloroplast psbA gene (psbA^ncr) to examine genetic diversity among clade C Symbiodinium associating with the common reef zoantharian Palythoa tuberculosa on Okinawa-jima Island, Japan. We identify four closely related Symbiodinium psbA^ncr lineages including one common generalist and two potential specialists that appear to be associated with particular microhabitats. The sea surface temperature differences that distinguish these habitats are smaller than those usually investigated, suggesting that future biogeographic surveys of Symbiodinium should incorporate fine scale environmental information as well as fine scale molecular data to accurately determine species diversity and their distributions.

Shifting communities after-- typhoon damage on an upper mesophotic reef in Okinawa, Japan

Very few studies have been conducted on the long-term effects of typhoon damage on mesophotic coral reefs. This study investigates the long-term community dynamics of damage from Typhoon 17 (Jelawat) in 2012 on the coral community of the upper mesophotic Ryugu Reef in Okinawa, Japan. A shift from foliose to bushy coral morphologies between December 2012 and August 2015 was documented, especially on the area of the reef that was previously recorded to be poor in scleractinian genera diversity and dominated by foliose corals. Comparatively, an area with higher diversity of scleractinian coral genera was observed to be less affected by typhoon damage with more stable community structure due to less change in dominant coral morphologies. Despite some changes in the composition of dominant genera, the generally high coverage of the mesophotic coral community is facilitating the recovery of Ryugu Reef after typhoon damage.

Contrasting patterns of population structure and gene flow facilitate exploration of connectivity in two widely distributed temperate octocorals

Connectivity is an important component of metapopulation dynamics in marine systems and can influence population persistence, migration rates and conservation decisions associated with Marine Protected Areas (MPAs). In this study, we compared the genetic diversity, gene flow and population structure of two octocoral species, Eunicella verrucosa and Alcyonium digitatum, in the northeast Atlantic (ranging from the northwest of Ireland and the southern North Sea, to southern Portugal), using two panels of 13 and 8 microsatellite loci, respectively. Our results identified regional genetic structure in E. verrucosa partitioned between populations from southern Portugal, northwest Ireland and Britain/France; subsequent hierarchical analysis of population structure also indicated reduced gene flow between southwest Britain and northwest France. However, over a similar geographical area, A. digitatum showed little evidence of population structure, suggesting high gene flow and/or a large effective population size; indeed, the only significant genetic differentiation detected in A. digitatum occurred between North Sea samples and those from the English Channel/northeast Atlantic. In both species the vast majority of gene flow originated from sample sites within regions, with populations in southwest Britain being the predominant source of contemporary exogenous genetic variants for the populations studied. Overall, historical patterns of gene flow appeared more complex, though again southwest Britain appeared to be an important source of genetic variation for both species. Our findings have major conservation implications, particularly for E. verrucosa, a protected species in UK waters and listed by the IUCN as ‘Vulnerable’, and for the designation and management of European MPAs.

Extremely low genetic variability within and among locations of the greenfish holothurian Stichopus chloronotus Brandt, 1835 in Okinawa, Japan

The greenfish sea cucumber Stichopus chloronotus is an economically and ecologically important sea cucumber species throughout its range. This species is widely distributed, inhabiting coral reefs of the Indo-Pacific Ocean. Our study evaluated population genetic structure and levels of genetic diversity in southern Japan. A total of 180 individuals were collected from eight locations from Okinawa and Okinoerabu Islands and sequenced using mitochondrial 16S ribosomal DNA (16S) and nuclear histone H3 (H3) gene. Only three 16S haplotypes were detected (518 bp) with haplotype diversity ranging from 0 to 0.56 and nucleotide diversity from 0 to 0.1%. H3 showed no variation among the studied locations. It is plausible that such results could be due to a shift to asexual reproduction. Additionally, the presence of the species on the east coast of Okinawa could only be detected in one location and all individuals consisted of a single haplotype. Genetic differences between the east and west coasts of Okinawa have been noticed in other coral reef organisms, and attributed to either ecological or biogeographical historical differences between the coasts due to differing levels of isolation during Pleistocene ice ages. Results from the present study should inform management and conservation policies of S. chloronotus in southern Japan.

The Mesoderm-Forming Gene brachyury Regulates Ectoderm-Endoderm Demarcation in the Coral Acropora digitifera

Blastoporal expression of the T-box transcription factor gene brachyury is conserved in most metazoans [1, 2]. Its role in mesoderm formation has been intensively studied in vertebrates [3-6]. However, its fundamental function near the blastopore is poorly understood in other phyla. Cnidarians are basal metazoans that are important for understanding evolution of metazoan body plans [7, 8]. Because they lack mesoderm, they have been used to investigate the evolutionary origins of mesoderm [1, 9-11]. Here, we focus on corals, a primitive clade of cnidarians that diverged from sea anemones ∼500 mya [12]. We developed a microinjection method for coral eggs to examine Brachyury functions during embryogenesis of the scleractinian coral, Acropora digitifera. Because Acropora embryos undergo pharynx formation after the blastopore closes completely [13-15], they are useful to understand Brachyury functions in gastrulation movement and pharynx formation. We show that blastoporal expression of brachyury is directly activated by Wnt/β-catenin signaling in the ectoderm of coral embryos, indicating that the regulatory axis from Wnt/β-catenin signaling to brachyury is highly conserved among eumetazoans. Loss-of-function analysis demonstrated that Brachyury is required for pharynx formation but not for gastrulation movement. Genome-wide transcriptome analysis demonstrated that genes positively regulated by Brachyury are expressed in the ectoderm of Acropora gastrulae, while negatively regulated genes are in endoderm. Therefore, germ layer demarcation around the blastopore appears to be the evolutionarily conserved role of Brachyury during gastrulation. Compared with Brachyury functions in vertebrate mesoderm-ectoderm and mesoderm-endoderm demarcation [4-6], our results suggest that the vertebrate-type mesoderm may have originated from brachyury-expressing ectoderm adjacent to endoderm.

Unexpectedly complex gradation of coral population structure in the Nansei Islands, Japan

To establish effective locations and sizes of potential protected areas for reef ecosystems, detailed information about source and sink relationships between populations is critical, especially in archipelagic regions. Therefore, we assessed population structure and genetic diversity of Acropora tenuis, one of the dominant stony coral species in the Pacific, using 13 microsatellite markers to investigate 298 colonies from 15 locations across the Nansei Islands in southwestern Japan. Genetic diversity was not significant among sampling locations, even in possibly peripheral locations. In addition, our results showed that there are at least two populations of A. tenuis in the study area. The level of genetic differentiation between these populations was relatively low, but significant between many pairs of sampling locations. Directions of gene flow, which were estimated using a coalescence-based approach, suggest that gene flow not only occurs from south to north, but also from north to south in various locations. Consequently, the Yaeyama Islands and the Amami Islands are potential northern and southern sources of corals. On the other hand, the Miyako Islands and west central Okinawa Island are potential sink populations. The Kerama Islands and the vicinity of Taketomi Island are potential contact points of genetic subdivision of coral populations in the Nansei Islands. We found that genetic population structure of A. tenuis in the Nansei Islands is more complex than previously thought. These cryptic populations are very important for preserving genetic diversity and should be maintained.