A genomic glance through the fog of plasticity and diversification in Pocillopora

Niche breadth and divergence in sympatric cryptic coral species (Pocillopora spp.) across habitats within reefs and among algal symbionts

While the presence of morphologically cryptic species is increasingly recognized, we still lack a useful understanding of what causes and maintains co-occurring cryptic species and its consequences for the ecology, evolution, and conservation of communities. We sampled 724 Pocillopora corals from five habitat zones (the fringing reef, back reef, and fore reef at 5, 10, and 20 m) at four sites around the island of Moorea, French Polynesia. Using validated genetic markers, we identified six sympatric species of Pocillopora, most of which cannot be reliably identified based on morphology: P. meandrina (42.9%), P. tuahiniensis (25.1%), P. verrucosa (12.2%), P. acuta (10.4%), P. grandis (7.73%), and P. cf. effusa (2.76%). For 423 colonies (58% of the genetically identified hosts), we also used psbA ncr or ITS2 markers to identify symbiont species (Symbiodiniaceae). The relative abundance of Pocillopora species differed across habitats within the reef. Sister taxa P. verrucosa and P. tuahiniensis had similar niche breadths and hosted the same specialist symbiont species (mostly Cladocopium pacificum) but the former was more common in the back reef and the latter more common deeper on the fore reef. In contrast, sister taxa P. meandrina and P. grandis had the highest niche breadths and overlaps and tended to host the same specialist symbiont species (mostly C. latusorum). Pocillopora acuta had the narrowest niche breadth and hosted the generalist, and more thermally tolerant, Durusdinium gynnii. Overall, there was a positive correlation between reef habitat niche breadth and symbiont niche breadth-Pocillopora species with a broader habitat niche also had a broader symbiont niche. Our results show how fine-scale variation within reefs plays an important role in the generation and coexistence of cryptic species. The results also have important implications for how niche differences affect community resilience, and for the success of coral restoration practices, in ways not previously appreciated.

Latitudinal variation in thermal performance of the common coral Pocillopora spp

Understanding how tropical corals respond to temperatures is important to evaluating their capacity to persist in a warmer future. We studied the common Pacific coral Pocillopora over 44° of latitude, and used populations at three islands with different thermal regimes to compare their responses to temperature using thermal performance curves (TPCs) for respiration and gross photosynthesis. Corals were sampled in the local autumn from Moorea, Guam and Okinawa, where mean±s.d. annual seawater temperature is 28.0±0.9°C, 28.9±0.7°C and 25.1±3.4°C, respectively. TPCs for respiration were similar among latitudes, the thermal optimum (Topt) was above the local maximum temperature at all three islands, and maximum respiration was lowest at Okinawa. TPCs for gross photosynthesis were wider, implying greater thermal eurytopy, with a higher Topt in Moorea versus Guam and Okinawa. Topt was above the maximum temperature in Moorea, but was similar to daily temperatures over 13% of the year in Okinawa and 53% of the year in Guam. There was greater annual variation in daily temperatures in Okinawa than Guam or Moorea, which translated to large variation in the supply of metabolic energy and photosynthetically fixed carbon at higher latitudes. Despite these trends, the differences in TPCs for Pocillopora spp. were not profoundly different across latitudes, reducing the likelihood that populations of these corals could better match their phenotypes to future more extreme temperatures through migration. Any such response would place a premium on high metabolic plasticity and tolerance of large seasonal variations in energy budgets.

Disparate genetic divergence patterns in three corals across a pan-Pacific environmental gradient highlight species-specific adaptation

Tropical coral reefs are among the most affected ecosystems by climate change and face increasing loss in the coming decades. Effective conservation strategies that maximize ecosystem resilience must be informed by the accurate characterization of extant genetic diversity and population structure together with an understanding of the adaptive potential of keystone species. Here we analyzed samples from the Tara Pacific Expedition (2016-2018) that completed an 18,000 km longitudinal transect of the Pacific Ocean sampling three widespread corals-Pocillopora meandrina, Porites lobata, and Millepora cf. platyphylla-across 33 sites from 11 islands. Using deep metagenomic sequencing of 269 colonies in conjunction with morphological analyses and climate variability data, we can show that despite a targeted sampling the transect encompasses multiple cryptic species. These species exhibit disparate biogeographic patterns and, most importantly, distinct evolutionary patterns in identical environmental regimes. Our findings demonstrate on a basin scale that evolutionary trajectories are species-specific and can only in part be predicted from the environment. This highlights that conservation strategies must integrate multi-species investigations to discern the distinct genomic footprints shaped by selection as well as the genetic potential for adaptive change.

From genomics to integrative species delimitation? The case study of the Indo-Pacific Pocillopora corals

With the advent of genomics, sequencing thousands of loci from hundreds of individuals now appears feasible at reasonable costs, allowing complex phylogenies to be resolved. This is particularly relevant for cnidarians, for which insufficient data is available due to the small number of currently available markers and obscures species boundaries. Difficulties in inferring gene trees and morphological incongruences further blur the study and conservation of these organisms. Yet, can genomics alone be used to delimit species? Here, focusing on the coral genus Pocillopora, whose colonies play key roles in Indo-Pacific reef ecosystems but have challenged taxonomists for decades, we explored and discussed the usefulness of multiple criteria (genetics, morphology, biogeography and symbiosis ecology) to delimit species of this genus. Phylogenetic inferences, clustering approaches and species delimitation methods based on genome-wide single-nucleotide polymorphisms (SNP) were first used to resolve Pocillopora phylogeny and propose genomic species hypotheses from 356 colonies sampled across the Indo-Pacific (western Indian Ocean, tropical southwestern Pacific and south-east Polynesia). These species hypotheses were then compared to other lines of evidence based on genetic, morphology, biogeography and symbiont associations. Out of 21 species hypotheses delimited by genomics, 13 were strongly supported by all approaches, while six could represent either undescribed species or nominal species that have been synonymised incorrectly. Altogether, our results support (1) the obsolescence of macromorphology (i.e., overall colony and branches shape) but the relevance of micromorphology (i.e., corallite structures) to refine Pocillopora species boundaries, (2) the relevance of the mtORF (coupled with other markers in some cases) as a diagnostic marker of most species, (3) the requirement of molecular identification when species identity of colonies is absolutely necessary to interpret results, as morphology can blur species identification in the field, and (4) the need for a taxonomic revision of the genus Pocillopora. These results give new insights into the usefulness of multiple criteria for resolving Pocillopora, and more widely, scleractinian species boundaries, and will ultimately contribute to the taxonomic revision of this genus and the conservation of its species.

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.

The effects of marine heatwaves on acute heat tolerance in corals

Scleractinian coral populations are increasingly exposed to conditions above their upper thermal limits due to marine heatwaves, contributing to global declines of coral reef ecosystem health. However, historic mass bleaching events indicate there is considerable inter- and intra-specific variation in thermal tolerance whereby species, individual coral colonies and populations show differential susceptibility to exposure to elevated temperatures. Despite this, we lack a clear understanding of how heat tolerance varies across large contemporary and historical environmental gradients, or the selective pressures that underpin this variation. Here we conducted standardised acute heat stress experiments to identify variation in heat tolerance among species and isolated reefs spanning a large environmental gradient across the Coral Sea Marine Park. We quantified the photochemical yield (F_v /F_m ) of coral samples in three coral species, Acropora cf humilis, Pocillopora meandrina, and Pocillopora verrucosa, following exposure to four temperature treatments (local ambient temperatures, and + 3°C, +6°C and + 9°C above local maximum monthly mean). We quantified the temperature at which F_v /F_m decreased by 50% (termed ED50) and used derived values to directly compare acute heat tolerance across reefs and species. The ED50 for Acropora was 0.4-0.7°C lower than either Pocillopora species, with a 0.3°C difference between the two Pocillopora species. We also recorded 0.9°C to 1.9°C phenotypic variation in heat tolerance among reefs within species, indicating spatial heterogeneity in heat tolerance across broad environmental gradients. Acute heat tolerance had a strong positive relationship to mild heatwave exposure over the past 35 years (since 1986) but was negatively related to recent severe heatwaves (2016-2020). Phenotypic variation associated with mild thermal history in local environments provides supportive evidence that marine heatwaves are selecting for tolerant individuals and populations; however, this adaptive potential may be compromised by the exposure to recent severe heatwaves.

High-quality genome assembles from key Hawaiian coral species

Background

Coral reefs house about 25% of marine biodiversity and are critical for the livelihood of many communities by providing food, tourism revenue, and protection from wave surge. These magnificent ecosystems are under existential threat from anthropogenic climate change. Whereas extensive ecological and physiological studies have addressed coral response to environmental stress, high-quality reference genome data are lacking for many of these species. The latter issue hinders efforts to understand the genetic basis of stress resistance and to design informed coral conservation strategies.

Results

We report genome assemblies from 4 key Hawaiian coral species, Montipora capitata, Pocillopora acuta, Pocillopora meandrina, and Porites compressa. These species, or members of these genera, are distributed worldwide and therefore of broad scientific and ecological importance. For M. capitata, an initial assembly was generated from short-read Illumina and long-read PacBio data, which was then scaffolded into 14 putative chromosomes using Omni-C sequencing. For P. acuta, P. meandrina, and P. compressa, high-quality assemblies were generated using short-read Illumina and long-read PacBio data. The P. acuta assembly is from a triploid individual, making it the first reference genome of a nondiploid coral animal.

Conclusions

These assemblies are significant improvements over available data and provide invaluable resources for supporting multiomics studies into coral biology, not just in Hawaiʻi but also in other regions, where related species exist. The P. acuta assembly provides a platform for studying polyploidy in corals and its role in genome evolution and stress adaptation in these organisms.

Cophylogeny and specificity between cryptic coral species (Pocillopora spp.) at Mo'orea and their symbionts (Symbiodiniaceae)

The congruence between phylogenies of tightly associated groups of organisms (cophylogeny) reflects evolutionary links between ecologically important interactions. However, despite being a classic example of an obligate symbiosis, tests of cophylogeny between scleractinian corals and their photosynthetic algal symbionts have been hampered in the past because both corals and algae contain genetically unresolved and morphologically cryptic species. Here, we studied co-occurring, cryptic Pocillopora species from Mo'orea, French Polynesia, that differ in their relative abundance across depth. We constructed new phylogenies of the host Pocillopora (using complete mitochondrial genomes, genomic loci, and thousands of single nucleotide polymorphisms) and their Symbiodiniaceae symbionts (using ITS2 and psbA^ncr markers) and tested for cophylogeny. The analysis supported the presence of five Pocillopora species on the fore reef at Mo'orea that mostly hosted either Cladocopium latusorum or C. pacificum. Only Pocillopora species hosting C. latusorum also hosted taxa from Symbiodinium and Durusdinium. In general, the Cladocopium phylogeny mirrored the Pocillopora phylogeny. Within Cladocopium species, lineages also differed in their associations with Pocillopora haplotypes, except those showing evidence of nuclear introgression, and with depth in the two most common Pocillopora species. We also found evidence for a new Pocillopora species (haplotype 10), that has so far only been sampled from French Polynesia, that warrants formal identification. The linked phylogenies of these Pocillopora and Cladocopium species and lineages suggest that symbiont speciation is driven by niche diversification in the host, but there is still evidence for symbiont flexibility in some cases.

Full-Length Transcriptome Maps of Reef-Building Coral Illuminate the Molecular Basis of Calcification, Symbiosis, and Circa-Dian Genes

Coral transcriptomic data largely rely on short-read sequencing, which severely limits the understanding of coral molecular mechanisms and leaves many important biological questions unresolved. Here, we sequence the full-length transcriptomes of four common and frequently dominant reef-building corals using the PacBio Sequel II platform. We obtain information on reported gene functions, structures, and expression profiles. Among them, a comparative analysis of biomineralization-related genes provides insights into the molecular basis of coral skeletal density. The gene expression profiles of the symbiont Symbiodiniaceae are also isolated and annotated from the holobiont sequence data. Finally, a phylogenetic analysis of key circadian clock genes among 40 evolutionarily representative species indicates that there are four key members in early metazoans, including cry genes; Clock or Npas2; cyc or Arntl; and tim, while per, as the fifth member, occurs in Bilateria. In summary, this work provides a foundation for further work on the manipulation of skeleton production or symbiosis to promote the survival of these important organisms.

Species and population genomic differentiation in Pocillopora corals (Cnidaria, Hexacorallia)

Correctly delimiting species and populations is a prerequisite for studies of connectivity, adaptation and conservation. Genomic data are particularly useful to test species differentiation for organisms with few informative morphological characters or low discrimination of cytoplasmic markers, as in Scleractinians. Here we applied Restriction site Associated DNA sequencing (RAD-sequencing) to the study of species differentiation and genetic structure in populations of Pocillopora spp. from Oman and French Polynesia, with the objectives to test species hypotheses, and to study the genetic structure among sampling sites within species. We focused here on coral colonies morphologically similar to P. acuta (damicornis type β). We tested the impact of different filtering strategies on the stability of the results. The main genetic differentiation was observed between samples from Oman and French Polynesia. These samples corresponded to different previously defined primary species hypotheses (PSH), i.e., PSHs 12 and 13 in Oman, and PSH 5 in French Polynesia. In Oman, we did not observe any clear differentiation between the two putative species PSH 12 and 13, nor between sampling sites. In French Polynesia, where a single species hypothesis was studied, there was no differentiation between sites. Our analyses allowed the identification of clonal lineages in Oman and French Polynesia. The impact of clonality on genetic diversity is discussed in light of individual-based simulations.

Thermal stress triggers productive viral infection of a key coral reef symbiont

Climate change-driven ocean warming is increasing the frequency and severity of bleaching events, in which corals appear whitened after losing their dinoflagellate endosymbionts (family Symbiodiniaceae). Viral infections of Symbiodiniaceae may contribute to some bleaching signs, but little empirical evidence exists to support this hypothesis. We present the first temporal analysis of a lineage of Symbiodiniaceae-infecting positive-sense single-stranded RNA viruses ("dinoRNAVs") in coral colonies, which were exposed to a 5-day heat treatment (+2.1 °C). A total of 124 dinoRNAV major capsid protein gene "aminotypes" (unique amino acid sequences) were detected from five colonies of two closely related Pocillopora-Cladocopium (coral-symbiont) combinations in the experiment; most dinoRNAV aminotypes were shared between the two coral-symbiont combinations (64%) and among multiple colonies (82%). Throughout the experiment, seventeen dinoRNAV aminotypes were found only in heat-treated fragments, and 22 aminotypes were detected at higher relative abundances in heat-treated fragments. DinoRNAVs in fragments of some colonies exhibited higher alpha diversity and dispersion under heat stress. Together, these findings provide the first empirical evidence that exposure to high temperatures triggers some dinoRNAVs to switch from a persistent to a productive infection mode within heat-stressed corals. Over extended time frames, we hypothesize that cumulative dinoRNAV production in the Pocillopora-Cladocopium system could affect colony symbiotic status, for example, by decreasing Symbiodiniaceae densities within corals. This study sets the stage for reef-scale investigations of dinoRNAV dynamics during bleaching events.

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.

Coral Genus Differentiation Based on Direct Analysis in Real Time-High Resolution Mass Spectrometry-Derived Chemical Fingerprints

Coral reefs are one of the most biologically diverse ecosystems, and the accurate identification of the species is essential for diversity assessment and conservation. Current genus determination approaches are time-consuming and resource-intensive and can be highly subjective. To explore the hypothesis that the small-molecule profiles of coral are genus-specific and can be used as a rapid tool to catalogue and distinguish between coral genera, the small-molecule chemical fingerprints of the species Acanthastrea echinata, Catalaphyllia jardinei, Duncanopsammia axifuga, Echinopora lamellosa, Euphyllia divisa, Euphyllia paraancora, Euphyllia paradivisa, Galaxea fascicularis, Herpolitha limax, Montipora confusa, Monitpora digitata, Montipora setosa, Pachyseris rugosa, Pavona cactus, Plerogyra sinuosa, Pocillopora acuta, Seriatopora hystrix, Sinularia dura, Turbinaria peltata, Turbinaria reniformis, Xenia elongata, and Xenia umbellata were generated using direct analysis in real time-high resolution mass spectrometry (DART-HRMS). It is demonstrated here that the mass spectrum-derived small-molecule profiles for coral of different genera are distinct. Multivariate statistical analysis processing of the DART-HRMS data enabled rapid genus-level differentiation based on the chemical composition of the coral. Coral samples were analyzed with no sample preparation required, making the approach rapid and efficient. The resulting spectra were subjected to kernel discriminant analysis (KDA), which furnished accurate genus differentiation of the coral. Leave-one-out cross-validation (LOOCV) was carried out to determine the classification accuracy of each model and confirm that this approach can be used for coral genus attribution with prediction accuracies ranging from 86.67 to 97.33%. The advantages and application of the statistical analysis to DART-HRMS-derived coral chemical signatures for genus-level differentiation are discussed.

Mutualistic microalgae co-diversify with reef corals that acquire symbionts during egg development

The application of molecular genetics has reinvigorated and improved how species are defined and investigated scientifically, especially for morphologically cryptic micro-organisms. Here we show how species recognition improves understanding of the ecology and evolution of mutualisms between reef-building corals and their mutualistic dinoflagellates (i.e. Symbiodiniaceae). A combination of genetic, ecological, and morphological evidence defines two sibling species of Cladocopium (formerly Symbiodinium Clade C), specific only to host corals in the common genus Pocillopora, which transmit their obligate symbionts during oogenesis. Cladocopium latusorum sp. nov. is symbiotic with P. grandis/meandrina while the smaller-celled C. pacificum sp. nov. associates with P. verrucosa. Both symbiont species form mutualisms with Pocillopora that brood their young. Populations of each species, like their hosts, are genetically well connected across the tropical and subtropical Pacific Ocean, indicating a capacity for long-range dispersal. A molecular clock approximates their speciation during the late Pliocene or early Pleistocene as Earth underwent cycles of precipitous cooling and warming; and corresponds to when their hosts were also diversifying. The long temporal and spatial maintenance of high host fidelity, as well as genetic connectivity across thousands of kilometers, indicates that distinct ecological attributes and close evolutionary histories will restrain the adaptive responses of corals and their specialized symbionts to rapid climate warming.

Characterizing the Genetic Architecture of Parkinson's Disease in Latinos

OBJECTIVE

This work was undertaken in order to identify Parkinson's disease (PD) risk variants in a Latino cohort, to describe the overlap in the genetic architecture of PD in Latinos compared to European-ancestry subjects, and to increase the diversity in PD genome-wide association (GWAS) data.

METHODS

We genotyped and imputed 1,497 PD cases and controls recruited from nine clinical sites across South America. We performed a GWAS using logistic mixed models; variants with a p-value <1 × 10^-5 were tested in a replication cohort of 1,234 self-reported Latino PD cases and 439,522 Latino controls from 23andMe, Inc. We also performed an admixture mapping analysis where local ancestry blocks were tested for association with PD status.

RESULTS

One locus, SNCA, achieved genome-wide significance (p-value <5 × 10^-8 ); rs356182 achieved genome-wide significance in both the discovery and the replication cohorts (discovery, G allele: 1.58 OR, 95% CI 1.35-1.86, p-value 2.48 × 10^-8 ; 23andMe, G allele: 1.26 OR, 95% CI 1.16-1.37, p-value 4.55 × 10^-8 ). In our admixture mapping analysis, a locus on chromosome 14, containing the gene STXBP6, achieved significance in a joint test of ancestries and in the Native American single-ancestry test (p-value <5 × 10^-5 ). A second locus on chromosome 6, containing the gene RPS6KA2, achieved significance in the African single-ancestry test (p-value <5 × 10^-5 ).

INTERPRETATION

This study demonstrated the importance of the SNCA locus for the etiology of PD in Latinos. By leveraging the demographic history of our cohort via admixture mapping, we identified two potential PD risk loci that merit further study. ANN NEUROL 2021.

Among-genotype responses of the coral Pocillopora acuta to emersion: implications for the ecological engineering of artificial coastal defences

Stony corals are promising transplant candidates for the ecological engineering of artificial coastal defences such as seawalls as they attract and host numerous other organisms. However, seawalls are exposed to a wide range of environmental stressors associated with periods of emersion during low tide such as desiccation and changes in salinity, temperature, and solar irradiance. All of these variables have known deleterious effects on coral physiology, growth, and fitness. In this study, we performed parallel experiments (in situ and ex situ) to examine among-genotype responses of Pocillopora acuta to emersion by quantifying growth, photophysiological metrics (F_v/F_m, non-photochemical quenching [NPQ], endosymbiont density, and chlorophyll [chl] a concentration) and survival, following different emersion periods. Results showed that coral fragments emersed for longer durations (>2 h) exhibited reduced growth and survival. Endosymbiont density and NPQ, but not F_v/F_m and chl a concentration, varied significantly among genotypes across different durations of emersion. Overall, the ability of P. acuta to tolerate emersion for up to 2 h suggests its potential to serve as a 'starter species' for transplantation efforts on seawalls. Further, careful characterisation and selection of genotypes with a high capacity to withstand emersion can help maximise the efficacy of ecological engineering using coral transplants.

Response diversity in corals: hidden differences in bleaching mortality among cryptic Pocillopora species

Variation among functionally similar species in their response to environmental stress buffers ecosystems from changing states. Functionally similar species may often be cryptic species representing evolutionarily distinct genetic lineages that are morphologically indistinguishable. However, the extent to which cryptic species differ in their response to stress, and could therefore provide a source of response diversity, remains unclear because they are often not identified or are assumed to be ecologically equivalent. Here, we uncover differences in the bleaching response between sympatric cryptic species of the common Indo-Pacific coral, Pocillopora. In April 2019, prolonged ocean heating occurred at Moorea, French Polynesia. 72% of pocilloporid colonies bleached after 22 d of severe heating (>8o C-days) at 10 m depth on the north shore fore reef. Colony mortality ranged from 11% to 42% around the island four months after heating subsided. The majority (86%) of pocilloporids that died from bleaching belonged to a single haplotype, despite twelve haplotypes, representing at least five species, being sampled. Mitochondrial (open reading frame) sequence variation was greater between the haplotypes that experienced mortality versus haplotypes that all survived than it was between nominal species that all survived. Colonies > 30 cm in diameter were identified as the haplotype experiencing the most mortality, and in 1125 colonies that were not genetically identified, bleaching and mortality increased with colony size. Mortality did not increase with colony size within the haplotype suffering the highest mortality, suggesting that size-dependent bleaching and mortality at the genus level was caused instead by differences among cryptic species. The relative abundance of haplotypes shifted between February and August, driven by declines in the same common haplotype for which mortality was estimated directly, at sites where heat accumulation was greatest, and where larger colony sizes occurred. The identification of morphologically indistinguishable species that differ in their response to thermal stress, but share a similar ecological function in terms of maintaining a coral-dominated state, has important consequences for uncovering response diversity that drives resilience, especially in systems with low or declining functional diversity.

Applying model approaches in non-model systems: A review and case study on coral cell culture

Model systems approaches search for commonality in patterns underlying biological diversity and complexity led by common evolutionary paths. The success of the approach does not rest on the species chosen but on the scalability of the model and methods used to develop the model and engage research. Fine-tuning approaches to improve coral cell cultures will provide a robust platform for studying symbiosis breakdown, the calcification mechanism and its disruption, protein interactions, micronutrient transport/exchange, and the toxicity of nanoparticles, among other key biological aspects, with the added advantage of minimizing the ethical conundrum of repeated testing on ecologically threatened organisms. The work presented here aimed to lay the foundation towards development of effective methods to sort and culture reef-building coral cells with the ultimate goal of obtaining immortal cell lines for the study of bleaching, disease and toxicity at the cellular and polyp levels. To achieve this objective, the team conducted a thorough review and tested the available methods (i.e. cell dissociation, isolation, sorting, attachment and proliferation). The most effective and reproducible techniques were combined to consolidate culture methods and generate uncontaminated coral cell cultures for ~7 days (10 days maximum). The tests were conducted on scleractinian corals Pocillopora acuta of the same genotype to harmonize results and reduce variation linked to genetic diversity. The development of cell separation and identification methods in conjunction with further investigations into coral cell-type specific metabolic requirements will allow us to tailor growth media for optimized monocultures as a tool for studying essential reef-building coral traits such as symbiosis, wound healing and calcification at multiple scales.

Phylogenomics of Porites from the Arabian Peninsula

The advent of high throughput sequencing technologies provides an opportunity to resolve phylogenetic relationships among closely related species. By incorporating hundreds to thousands of unlinked loci and single nucleotide polymorphisms (SNPs), phylogenomic analyses have a far greater potential to resolve species boundaries than approaches that rely on only a few markers. Scleractinian taxa have proved challenging to identify using traditional morphological approaches and many groups lack an adequate set of molecular markers to investigate their phylogenies. Here, we examine the potential of Restriction-site Associated DNA sequencing (RADseq) to investigate phylogenetic relationships and species limits within the scleractinian coral genus Porites. A total of 126 colonies were collected from 16 localities in the seas surrounding the Arabian Peninsula and ascribed to 12 nominal and two unknown species based on their morphology. Reference mapping was used to retrieve and compare nearly complete mitochondrial genomes, ribosomal DNA, and histone loci. De novo assembly and reference mapping to the P. lobata coral transcriptome were compared and used to obtain thousands of genome-wide loci and SNPs. A suite of species discovery methods (phylogenetic, ordination, and clustering analyses) and species delimitation approaches (coalescent-based, species tree, and Bayesian Factor delimitation) suggested the presence of eight molecular lineages, one of which included six morphospecies. Our phylogenomic approach provided a fully supported phylogeny of Porites from the Arabian Peninsula, suggesting the power of RADseq data to solve the species delineation problem in this speciose coral genus.

Seascape genomics reveals candidate molecular targets of heat stress adaptation in three coral species

Anomalous heat waves are causing a major decline of hard corals around the world and threatening the persistence of coral reefs. There are, however, reefs that have been exposed to recurrent thermal stress over the years and whose corals appear to have been tolerant against heat. One of the mechanisms that could explain this phenomenon is local adaptation, but the underlying molecular mechanisms are poorly known. In this work, we applied a seascape genomics approach to study heat stress adaptation in three coral species of New Caledonia (southwestern Pacific) and to uncover the molecular actors potentially involved. We used remote sensing data to characterize the environmental trends across the reef system, and sampled corals living at the most contrasted sites. These samples underwent next generation sequencing to reveal single nucleotide polymorphisms (SNPs), frequencies of which were associated with heat stress gradients. As these SNPs might underpin an adaptive role, we characterized the functional roles of the genes located in their genomic region. In each of the studied species, we found heat stress-associated SNPs located in proximity of genes involved in pathways well known to contribute to the cellular responses against heat, such as protein folding, oxidative stress homeostasis, inflammatory and apoptotic pathways, and DNA damage-repair. In some cases, the same candidate molecular targets of heat stress adaptation recurred among species. Together, these results underline the relevance and the power of the seascape genomics approach for the discovery of adaptive traits that could allow corals to persist across wider thermal ranges.

Functional diversity of microboring Ostreobium algae isolated from corals

The filamentous chlorophyte Ostreobium sp. dominates shallow marine carbonate microboring communities, and is one of the major agents of reef bioerosion. While its large genetic diversity has emerged, its physiology remains little known, with unexplored relationship between genotypes and phenotypes (endolithic versus free-living growth forms). Here, we isolated nine strains affiliated to two lineages of Ostreobium (>8% sequence divergence of the plastid gene rbcL), one of which was assigned to the family Odoaceae, from the fast-growing coral host Pocillopora acuta Lamarck 1816. Free-living isolates maintained their bioerosive potential, colonizing pre-bleached coral carbonate skeletons. We compared phenotypes, highlighting shifts in pigment and fatty acid compositions, carbon to nitrogen ratios and stable isotope compositions (δ¹³ C and δ¹⁵ N). Our data show a pattern of higher chlorophyll b and lower arachidonic acid (20:4ω6) content in endolithic versus free-living Ostreobium. Photosynthetic carbon fixation and nitrate uptake, quantified via 8 h pulse-labeling with ¹³ C-bicarbonate and ¹⁵ N-nitrate, showed lower isotopic enrichment in endolithic compared to free-living filaments. Our results highlight the functional plasticity of Ostreobium phenotypes. The isotope tracer approach opens the way to further study the biogeochemical cycling and trophic ecology of these cryptic algae at coral holobiont and reef scales.

Host-symbiont coevolution, cryptic structure, and bleaching susceptibility, in a coral species complex (Scleractinia; Poritidae)

The 'species' is a key concept for conservation and evolutionary biology, yet the lines between population and species-level variation are often blurred, especially for corals. The 'Porites lobata species complex' consists of branching and mounding corals that form reefs across the Pacific. We used reduced representation meta-genomic sequencing to examine genetic relationships within this species complex and to identify candidate loci associated with colony morphology, cryptic genetic structure, and apparent bleaching susceptibility. We compared existing Porites data with bleached and unbleached colonies of the branching coral P. compressa collected in Kāne'ohe Bay Hawai'i during the 2015 coral bleaching event. Loci that mapped to coral, symbiont, and microbial references revealed genetic structure consistent with recent host-symbiont co-evolution. Cryptic genetic clades were resolved that previous work has associated with distance from shore, but no genetic structure was associated with bleaching. We identified many candidate loci associated with morphospecies, including candidate host and symbiont loci with fixed differences between branching and mounding corals. We also found many loci associated with cryptic genetic structure, yet relatively few loci associated with bleaching. Recent host-symbiont co-evolution and rapid diversification suggests that variation and therefore the capacity of these corals to adapt may be underappreciated.

A framework for in situ molecular characterization of coral holobionts using nanopore sequencing

Molecular characterization of the coral host and the microbial assemblages associated with it (referred to as the coral holobiont) is currently undertaken via marker gene sequencing. This requires bulky instruments and controlled laboratory conditions which are impractical for environmental experiments in remote areas. Recent advances in sequencing technologies now permit rapid sequencing in the field; however, development of specific protocols and pipelines for the effective processing of complex microbial systems are currently lacking. Here, we used a combination of 3 marker genes targeting the coral animal host, its symbiotic alga, and the associated bacterial microbiome to characterize 60 coral colonies collected and processed in situ, during the Tara Pacific expedition. We used Oxford Nanopore Technologies to sequence marker gene amplicons and developed bioinformatics pipelines to analyze nanopore reads on a laptop, obtaining results in less than 24 h. Reef scale network analysis of coral-associated bacteria reveals broadly distributed taxa, as well as host-specific associations. Protocols and tools used in this work may be applicable for rapid coral holobiont surveys, immediate adaptation of sampling strategy in the field, and to make informed and timely decisions in the context of the current challenges affecting coral reefs worldwide.

New approaches to species delimitation and population structure of anthozoans: Two case studies of octocorals using ultraconserved elements and exons

As coral populations decline worldwide in the face of ongoing environmental change, documenting their distribution, diversity and conservation status is now more imperative than ever. Accurate delimitation and identification of species is a critical first step. This task, however, is not trivial as morphological variation and slowly evolving molecular markers confound species identification. New approaches to species delimitation in corals are needed to overcome these challenges. Here, we test whether target enrichment of ultraconserved elements (UCEs) and exons can be used for delimiting species boundaries and population structure within species of corals by focusing on two octocoral genera, Alcyonium and Sinularia, as exemplary case studies. We designed an updated bait set (29,181 baits) to target-capture 3,023 UCE and exon loci, recovering a mean of 1,910 ± 168 SD per sample with a mean length of 1,055 ± 208 bp. Similar numbers of loci were recovered from Sinularia (1,946 ± 227 SD) and Alcyonium (1,863 ± 177 SD). Species-level phylogenies were highly supported for both genera. Clustering methods based on filtered single nucleotide polymorphisms delimited species and populations that are congruent with previous allozyme, DNA barcoding, reproductive and ecological data for Alcyonium, and offered further evidence of hybridization among species. For Sinularia, results were congruent with those obtained from a previous study using restriction site associated DNA sequencing. Both case studies demonstrate the utility of target-enrichment of UCEs and exons to address a wide range of evolutionary and taxonomic questions across deep to shallow timescales in corals.

Gene Expression and Photophysiological Changes in Pocillopora acuta Coral Holobiont Following Heat Stress and Recovery

The ability of corals to withstand changes in their surroundings is a critical survival mechanism for coping with environmental stress. While many studies have examined responses of the coral holobiont to stressful conditions, its capacity to reverse responses and recover when the stressor is removed is not well-understood. In this study, we investigated among-colony responses of Pocillopora acuta from two sites with differing distance to the mainland (Kusu (closer to the mainland) and Raffles Lighthouse (further from the mainland)) to heat stress through differential expression analysis of target genes and quantification of photophysiological metrics. We then examined how these attributes were regulated after the stressor was removed to assess the recovery potential of P. acuta. The fragments that were subjected to heat stress (2 °C above ambient levels) generally exhibited significant reduction in their endosymbiont densities, but the extent of recovery following stress removal varied depending on natal site and colony. There were minimal changes in chl a concentration and maximum quantum yield (Fv/Fm, the proportion of variable fluorescence (Fv) to maximum fluorescence (Fm)) in heat-stressed corals, suggesting that the algal endosymbionts’ Photosystem II was not severely compromised. Significant changes in gene expression levels of selected genes of interest (GOI) were observed following heat exposure and stress removal among sites and colonies, including Actin, calcium/calmodulin-dependent protein kinase type IV (Camk4), kinesin-like protein (KIF9), and small heat shock protein 16.1 (Hsp16.1). The most responsive GOIs were Actin, a major component of the cytoskeleton, and the adaptive immune-related Camk4 which both showed significant reduction following heat exposure and subsequent upregulation during the recovery phase. Our findings clearly demonstrate specific responses of P. acuta in both photophysiological attributes and gene expression levels, suggesting differential capacity of P. acuta corals to tolerate heat stress depending on the colony, so that certain colonies may be more resilient than others.

Lipopolysaccharide treatment stimulates Pocillopora coral genotype-specific immune responses but does not alter coral-associated bacteria communities

Corals are comprised of a coral host and associated microbes whose interactions are mediated by the coral innate immune system. The diversity of immune factors identified in the Pocillopora damicornis genome suggests that immunity is linked to maintaining microbial symbioses while also being able to detect pathogens. However, it is unclear which immune factors respond to specific microbe-associated molecular patterns and how these immune reactions simultaneously affect coral-associated bacteria. To investigate this, fragments of P. damicornis and P. acuta colonies from Taiwan were subjected to lipopolysaccharide (LPS) treatment to stimulate immune responses and measure bacteria community shifts. RNA-seq revealed genotype-specific immune responses to LPS involving the upregulation of immune receptors, transcription factors, and pore-forming toxins. Bacteria 16S sequencing revealed significantly different bacteria communities between coral genotypes but no differences in bacteria communities were caused by LPS. Our findings confirm that Pocillopora corals activate conserved immune factors in response to LPS and identify transcription factors coordinating Pocillopora corals' immune responses. Additionally, the strong effect of coral genotype on gene expression and bacteria communities highlights the importance of coral genotype in the investigation of coral host-microbe interactions.

The 3D Reconstruction of Pocillopora Colony Sheds Light on the Growth Pattern of This Reef-Building Coral

Coral reefs are formed by living polyps, and understanding the dynamic processes behind the reefs is crucial for marine ecosystem restoration. However, these processes are still unclear because the growth and budding patterns of living polyps are poorly known. Here, we investigate the growth pattern of a widely distributed reef-building coral Pocillopora damicornis from Xisha Islands using high-resolution computed tomography. We examine the corallites in a single corallum of the species in detail, to interpret the budding, growth, and distribution pattern of the polyps, to reconstruct the growth pattern of this important reef-building species. Our results reveal a three-stage growth pattern of P. damicornis, based on different growth bundles that are secreted by polyps along the dichotomous growth axes of the corallites. Our work on the three-dimensional reconstruction of calice and inter-septal space structure of P. damicornis sheds lights on its reef-building processes by reconstructing the budding patterns.

•

We use high-resolution computed tomography to investigate coral forming and polyp budding processes

•

The calice reconstruction shows coral growth patterns and budding information

•

Our work visualizes the growth pattern of Pocillopora damicornis

•

High-resolution computed tomography is a method for future reef-building coral studies

Epifaunal invertebrate assemblages associated with branching Pocilloporids in Moorea, French Polynesia

Reef-building corals can harbour high abundances of diverse invertebrate epifauna. Coral characteristics and environmental conditions are important drivers of community structure of coral-associated invertebrates; however, our current understanding of drivers of epifaunal distributions is still unclear. This study tests the relative importance of the physical environment (current flow speed) and host quality (e.g., colony height, surface area, distance between branches, penetration depth among branches, and background partial mortality) in structuring epifaunal communities living within branching Pocillopora colonies on a back reef in Moorea, French Polynesia. A total of 470 individuals belonging to four phyla, 16 families and 39 genera were extracted from 36 Pocillopora spp. colonies. Decapods were the most abundant epifaunal organisms (accounting for 84% of individuals) found living in Pocillopora spp. While coral host characteristics and flow regime are very important, these parameters were not correlated with epifaunal assemblages at the time of the study. Epifaunal assemblages associated with Pocillopora spp. were consistent and minimally affected by differences in host characteristics and flow regime. The consistency in abundance and taxon richness among colonies (regardless of habitat characteristics) highlighted the importance of total habitat availability. With escalating effects of climate change and other localized disturbances, it is critical to preserve branching corals to support epifaunal communities.

Transcriptome-based target-enrichment baits for stony corals (Cnidaria: Anthozoa: Scleractinia)

Despite the ecological and economic significance of stony corals (Scleractinia), a robust understanding of their phylogeny remains elusive due to patchy taxonomic and genetic sampling, as well as the limited availability of informative markers. To increase the number of genetic loci available for phylogenomic analyses in Scleractinia, we designed 15,919 DNA enrichment baits targeting 605 orthogroups (mean 565 ± SD 366 bp) over 1,139 exon regions. A further 236 and 62 barcoding baits were designed for COI and histone H3 genes respectively for quality and contamination checks. Hybrid capture using these baits was performed on 18 coral species spanning the presently understood scleractinian phylogeny, with two corallimorpharians as outgroup. On average, 74% of all loci targeted were successfully captured for each species. Barcoding baits were matched unambiguously to their respective samples and revealed low levels of cross-contamination in accordance with expectation. We put the data through a series of stringent filtering steps to ensure only scleractinian and phylogenetically informative loci were retained, and the final probe set comprised 13,479 baits, targeting 452 loci (mean 531 ± SD 307 bp) across 865 exon regions. Maximum likelihood, Bayesian and species tree analyses recovered maximally supported, topologically congruent trees consistent with previous phylogenomic reconstructions. The phylogenomic method presented here allows for consistent capture of orthologous loci among divergent coral taxa, facilitating the pooling of data from different studies and increasing the phylogenetic sampling of scleractinians in the future.

Living on the edge: Assessing the diversity of South African Pocillopora on the margins of the Southwestern Indian Ocean

Scleractinia of the Maputaland reef complex (MRC) in South Africa exist at the margins of the Western Indian Ocean (WIO) coral distribution and are the only substantial hermatypic coral communities in South Africa. Pocillopora species occupy a conspicuous component of the MRC, and previous investigations identified three species of Pocillopora utilizing conventional taxonomy. Thus, our aims were four-fold: to elucidate Pocillopora species diversity using genetic techniques, primarily using species delimitation methods based on the ORF gene; to test for the presence of hybridisation within the Pocillopora community on the South-West margin of distribution in the Indian Ocean using two nuclear and two mitochondrial markers; to test the presence of cryptic species, using 13 microsatellite markers, finally, to elucidate the degree of genetic diversity within each Pocillopora species found and compare this to communities in lower latitudes. We illustrate taxonomic inconsistencies between these inventories and our phylogenetic data. The MRC harbours unique populations of Pocillopora, consisting of three species hypothetically co-occurring throughout the south WIO, namely: P. meandrina/P. eydouxi, commonly misidentified as P. verrucosa, P. verrucosa, sometimes correctly identified, but also commonly misidentified as P. damicornis sensu lato, and P. villosa, almost always misidentified as P. eydouxi. The hypothesis that hybrid swarms of Pocillopora occur in marginal environments such as the MRC was not supported, with low levels of introgressive hybridization reported instead. Analyses illustrate low genetic diversity at the species and population resolutions, suggesting a small founder population for each species. Nevertheless, these populations are demographically unique, exhibiting high levels of ITS2 haplotype endemism compared to higher latitude populations and the rest of the WIO. Pocillopora diversity on the MRC represents a unique assemblage and warrants further protection.

Rare coral under the genomic microscope: timing and relationships among Hawaiian Montipora

Background

Evolutionary patterns of scleractinian (stony) corals are difficult to infer given the existence of few diagnostic characters and pervasive phenotypic plasticity. A previous study of Hawaiian Montipora (Scleractinia: Acroporidae) based on five partial mitochondrial and two nuclear genes revealed the existence of a species complex, grouping one of the rarest known species (M. dilatata, which is listed as Endangered by the International Union for Conservation of Nature - IUCN) with widespread corals of very different colony growth forms (M. flabellata and M. cf. turgescens). These previous results could result from a lack of resolution due to a limited number of markers, compositional heterogeneity or reflect biological processes such as incomplete lineage sorting (ILS) or introgression.

Results

All 13 mitochondrial protein-coding genes from 55 scleractinians (14 lineages from this study) were used to evaluate if a recent origin of the M. dilatata species complex or rate heterogeneity could be compromising phylogenetic inference. Rate heterogeneity detected in the mitochondrial data set seems to have no significant impacts on the phylogenies but clearly affects age estimates. Dating analyses show different estimations for the speciation of M. dilatata species complex depending on whether taking compositional heterogeneity into account (0.8 [0.05–2.6] Myr) or assuming rate homogeneity (0.4 [0.14–0.75] Myr). Genomic data also provided evidence of introgression among all analysed samples of the complex. RADseq data indicated that M. capitata colour morphs may have a genetic basis.

Conclusions

Despite the volume of data (over 60,000 SNPs), phylogenetic relationships within the M. dilatata species complex remain unresolved most likely due to a recent origin and ongoing introgression. Species delimitation with genomic data is not concordant with the current taxonomy, which does not reflect the true diversity of this group. Nominal species within the complex are either undergoing a speciation process or represent ecomorphs exhibiting phenotypic polymorphisms.

Electronic supplementary material

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

A next generation approach to species delimitation reveals the role of hybridization in a cryptic species complex of corals

BACKGROUND

Our ability to investigate processes shaping the evolutionary diversification of corals (Cnidaria: Anthozoa) is limited by a lack of understanding of species boundaries. Discerning species of corals has been challenging due to a multitude of factors, including homoplasious and plastic morphological characters and the use of molecular markers that are either not informative or have not completely sorted. Hybridization can also blur species boundaries by leading to incongruence between morphology and genetics. We used traditional DNA barcoding and restriction-site associated DNA sequencing combined with coalescence-based and allele-frequency methods to elucidate species boundaries and simultaneously examine the potential role of hybridization in a speciose genus of octocoral, Sinularia.

RESULTS

Species delimitations using two widely used DNA barcode markers, mtMutS and 28S rDNA, were incongruent with one another and with the morphospecies identifications. When mtMutS and 28S were concatenated, a 0.3% genetic distance threshold delimited the majority of morphospecies. In contrast, 12 of the 15 examined morphospecies formed well-supported monophyletic clades in both concatenated RAxML phylogenies and SNAPP species trees of > 6000 RADSeq loci. DAPC and Structure analyses also supported morphospecies assignments, but indicated the potential for two additional cryptic species. Three morphologically distinct species pairs could not, however, be distinguished genetically. ABBA-BABA tests demonstrated significant admixture between some of those species, suggesting that hybridization may confound species delimitation in Sinularia.

CONCLUSIONS

A genomic approach can help to guide species delimitation while simultaneously elucidating the processes generating coral diversity. Results support the hypothesis that hybridization is an important mechanism in the evolution of Anthozoa, including octocorals, and future research should examine the contribution of this mechanism in generating diversity across the coral tree of life.

Annotated checklist for stony corals of American Sāmoa with reference to mesophotic depth records

An annotated checklist of the stony corals (Scleractinia, Milleporidae, Stylasteridae, and Helioporidae) of American Sāmoa is presented. A total of 377 valid species has been reported from American Sāmoa with 342 species considered either present (251) or possibly present (91). Of these 342 species, 66 have a recorded geographical range extension and 90 have been reported from mesophotic depths (30-150 m). Additionally, four new species records (Acanthastreasubechinata Veron, 2000, Favitesparaflexuosus Veron, 2000, Echinophylliaechinoporoides Veron & Pichon, 1980, Turbinariairregularis Bernard, 1896) are presented. Coral species of concern include species listed under the US Endangered Species Act (ESA) and the International Union for Conservation of Nature's (IUCN) Red List of threatened species. Approximately 17.5% of the species present or possibly present are categorized as threatened by IUCN compared to 27% of the species globally. American Sāmoa has seven ESA-listed or ESA candidate species, including Acroporaglobiceps (Dana, 1846), Acroporajacquelineae Wallace, 1994, Acroporaretusa (Dana, 1846), Acroporaspeciosa (Quelch, 1886), Fimbriaphylliaparadivisa (Veron, 1990), Isoporacrateriformis (Gardiner, 1898), and Pocilloporameandrina Dana, 1846. There are two additional species possibly present, i.e., Pavonadiffluens (Lamarck, 1816) and Poritesnapopora Veron, 2000.

Computational Characterization of the mtORF of Pocilloporid Corals: Insights into Protein Structure and Function in Stylophora Lineages from Contrasting Environments

More than a decade ago, a new mitochondrial Open Reading Frame (mtORF) was discovered in corals of the family Pocilloporidae and has been used since then as an effective barcode for these corals. Recently, mtORF sequencing revealed the existence of two differentiated Stylophora lineages occurring in sympatry along the environmental gradient of the Red Sea (18.5°C to 33.9°C). In the endemic Red Sea lineage RS_LinB, the mtORF and the heat shock protein gene hsp70 uncovered similar phylogeographic patterns strongly correlated with environmental variations. This suggests that the mtORF too might be involved in thermal adaptation. Here, we used computational analyses to explore the features and putative function of this mtORF. In particular, we tested the likelihood that this gene encodes a functional protein and whether it may play a role in adaptation. Analyses of full mitogenomes showed that the mtORF originated in the common ancestor of Madracis and other pocilloporids, and that it encodes a transmembrane protein differing in length and domain architecture among genera. Homology-based annotation and the relative conservation of metal-binding sites revealed traces of an ancient hydrolase catalytic activity. Furthermore, signals of pervasive purifying selection, lack of stop codons in 1830 sequences analyzed, and a codon-usage bias similar to that of other mitochondrial genes indicate that the protein is functional, i.e., not a pseudogene. Other features, such as intrinsically disordered regions, tandem repeats, and signals of positive selection particularly in StylophoraRS_LinB populations, are consistent with a role of the mtORF in adaptive responses to environmental changes.

Comparative analysis of the Pocillopora damicornis genome highlights role of immune system in coral evolution

Comparative analysis of the expanding genomic resources for scleractinian corals may provide insights into the evolution of these organisms, with implications for their continued persistence under global climate change. Here, we sequenced and annotated the genome of Pocillopora damicornis, one of the most abundant and widespread corals in the world. We compared this genome, based on protein-coding gene orthology, with other publicly available coral genomes (Cnidaria, Anthozoa, Scleractinia), as well as genomes from other anthozoan groups (Actiniaria, Corallimorpharia), and two basal metazoan outgroup phlya (Porifera, Ctenophora). We found that 46.6% of P. damicornis genes had orthologs in all other scleractinians, defining a coral ‘core’ genome enriched in basic housekeeping functions. Of these core genes, 3.7% were unique to scleractinians and were enriched in immune functionality, suggesting an important role of the immune system in coral evolution. Genes occurring only in P. damicornis were enriched in cellular signaling and stress response pathways, and we found similar immune-related gene family expansions in each coral species, indicating that immune system diversification may be a prominent feature of scleractinian coral evolution at multiple taxonomic levels. Diversification of the immune gene repertoire may underlie scleractinian adaptations to symbiosis, pathogen interactions, and environmental stress.

Molecular assessment of Pocillopora verrucosa (Scleractinia; Pocilloporidae) distribution along a depth gradient in Ludao, Taiwan

It can be challenging to identify scleractinian corals from the genus Pocillopora Lamarck 1816 in the field because of their large range of inter- and intra-specific morphological variation that co-occur with changes in the physical environment. This task is made more arduous in the context of a depth gradient, where light and water current could greatly affect the morphology of the corallum. Pocillopora verrucosa (Ellis & Solander 1786) in Taiwan was previously reported exclusively from shallow waters (<10 m in depth), but a recent observation of this species in the mesophotic zone (>40 m in depth) questions this bathymetric distribution. We used the mitochondrial open reading frame and the histone 3 molecular markers to investigate the vertical and horizontal spatial distribution of P. verrucosa around Ludao (Green Island), Taiwan. We genotyped 101 P. verrucosa-like colonies collected from four depth zones, ranging from 7 to 45 m, at three sites around the island. Of the 101 colonies sampled, 85 were genotyped as P. verrucosa, 15 as P. meandrina, and one specimen as an undescribed Pocillopora species. P. verrucosa was found at all depths, while P. meandrina and the undescribed Pocillopora specimen were limited to 15 m depth. P. verrucosa has a large bathymetric distribution around Ludao and could benefit from the refuge that the mesophotic zone offers. This study illustrates the difficulty of identifying Pocillopora corals in the field and emphasizes the relevance of molecular taxonomy as an important and complementary tool to traditional taxonomy for clarifying vertical and horizontal species distribution. Our results also illustrate the need in conservation biology to target species genetic diversity rather than just species diversity.

Spatial variation in the morphological traits of Pocillopora verrucosa along a depth gradient in Taiwan

Pocillopora verrucosa is a widely distributed depth-generalist coral that presents plasticity in its skeletal macro- and microstructure in response to environmental gradients. Light and water movement, which covary with depth, are the main environmental drivers of morphological plasticity in this genus; however, assessing environmentally-induced plasticity may be confounded by the extent of interspecific variation in Pocillopora. We examine the morphology of 8 typed P. verrucosa specimens collected along a depth gradient ranging from 7 to 45 meters and comprising 3 sites throughout Ludao, Taiwan. We measured 36 morphological characters, 14 which are novel, in 3 regions on the corallum-the apex, branch and base-in order to quantify their relationship to site and depth. We found significant correlation between depth and 19 morphological characters, notably branch verruca area, branch verruca height, base verruca spacing, base spinule length, and branch corallite area. 60% of microstructural characters and 25% of macrostructural characters showed a correlative relation to depth, suggesting that depth acclimatization is manifested primarily at the microstructural level. Canonical discriminant analysis of all morphometric characters by depth supports clustering into 3 groups: an overlapping 7m and 15m group, a 23-30m group, and a 38-45m group. Canonical discriminant analysis by site supports clustering into low- and high-current sites, differentiated primarily by branch septa width, base septa width, pre-terminal branch width, terminal branch maximum length, and terminal branch minimum length. We conclude that distinctive patterns of morphological variation in mesophotic specimens of P. verrucosa could reflect the effects of abiotic parameters such as light and water flow. Elucidating the mechanisms behind the morphological changes that occur in response to environmental gradients can help clarify the role that physiological plasticity plays in the acclimatization of corals to the unique environmental settings of mesophotic coral ecosystems.

Complex signatures of genomic variation of two non-model marine species in a homogeneous environment

BACKGROUND

Genomic tools are increasingly being used on non-model organisms to provide insights into population structure and variability, including signals of selection. However, most studies are carried out in regions with distinct environmental gradients or across large geographical areas, in which local adaptation is expected to occur. Therefore, the focus of this study is to characterize genomic variation and selective signals over short geographic areas within a largely homogeneous region. To assess adaptive signals between microhabitats within the rocky shore, we compared genomic variation between the Cape urchin (Parechinus angulosus), which is a low to mid-shore species, and the Granular limpet (Scutellastra granularis), a high shore specialist.

RESULTS

Using pooled restriction site associated DNA (RAD) sequencing, we described patterns of genomic variation and identified outlier loci in both species. We found relatively low numbers of outlier SNPs within each species, and identified outlier genes associated with different selective pressures than those previously identified in studies conducted over larger environmental gradients. The number of population-specific outlier loci differed between species, likely owing to differential selective pressures within the intertidal environment. Interestingly, the outlier loci were highly differentiated within the two northernmost populations for both species, suggesting that unique evolutionary forces are acting on marine invertebrates within this region.

CONCLUSIONS

Our study provides a background for comparative genomic studies focused on non-model species, as well as a baseline for the adaptive potential of marine invertebrates along the South African west coast. We also discuss the caveats associated with Pool-seq and potential biases of sequencing coverage on downstream genomic metrics. The findings provide evidence of species-specific selective pressures within a homogeneous environment, and suggest that selective forces acting on small scales are just as crucial to acknowledge as those acting on larger scales. As a whole, our findings imply that future population genomic studies should expand from focusing on model organisms and/or studying heterogeneous regions to better understand the evolutionary processes shaping current and future biodiversity patterns, particularly when used in a comparative phylogeographic context.

Quantification of complex modular architecture in plants

Morphometrics, the assignment of quantities to biological shapes, is a powerful tool to address taxonomic, evolutionary, functional and developmental questions. We propose a novel method for shape quantification of complex modular architecture in thalloid plants, whose extremely reduced morphologies, combined with the lack of a formal framework for thallus description, have long rendered taxonomic and evolutionary studies extremely challenging. Using graph theory, thalli are described as hierarchical series of nodes and edges, allowing for accurate, homologous and repeatable measurements of widths, lengths and angles. The computer program MorphoSnake was developed to extract the skeleton and contours of a thallus and automatically acquire, at each level of organization, width, length, angle and sinuosity measurements. Through the quantification of leaf architecture in Hymenophyllum ferns (Polypodiopsida) and a fully worked example of integrative taxonomy in the taxonomically challenging thalloid liverwort genus Riccardia, we show that MorphoSnake is applicable to all ramified plants. This new possibility of acquiring large numbers of quantitative traits in plants with complex modular architectures opens new perspectives of applications, from the development of rapid species identification tools to evolutionary analyses of adaptive plasticity.

A simple molecular technique for distinguishing species reveals frequent misidentification of Hawaiian corals in the genus Pocillopora

Species within the scleractinian genus Pocillopora Lamarck 1816 exhibit extreme phenotypic plasticity, making identification based on morphology difficult. However, the mitochondrial open reading frame (mtORF) marker provides a useful genetic tool for identification of most species in this genus, with a notable exception of P. eydouxi and P. meandrina. Based on recent genomic work, we present a quick and simple, gel-based restriction fragment length polymorphism (RFLP) method for the identification of all six Pocillopora species occurring in Hawai‘i by amplifying either the mtORF region, a newly discovered histone region, or both, and then using the restriction enzymes targeting diagnostic sequences we unambiguously identify each species. Using this approach, we documented frequent misidentification of Pocillopora species based on colony morphology. We found that P. acuta colonies are frequently mistakenly identified as P. damicornis in Kāne‘ohe Bay, O‘ahu. We also found that P. meandrina likely has a northern range limit in the Northwest Hawaiian Islands, above which P. ligulata was regularly mistaken for P. meandrina.

From population connectivity to the art of striping Russian dolls: the lessons from Pocillopora corals

Here, we examined the genetic variability in the coral genus Pocillopora, in particular within the Primary Species Hypothesis PSH09, identified by Gélin, Postaire, Fauvelot and Magalon (2017) using species delimitation methods [also named Pocillopora eydouxi/meandrina complex sensu, Schmidt-Roach, Miller, Lundgren, & Andreakis (2014)] and which was found to split into three secondary species hypotheses (SSH09a, SSH09b, and SSH09c) according to assignment tests using multi-locus genotypes (13 microsatellites). From a large sampling (2,507 colonies) achieved in three marine provinces [Western Indian Ocean (WIO), Tropical Southwestern Pacific (TSP), and Southeast Polynesia (SEP)], genetic structuring analysis conducted with two clustering analyses (structure and DAPC) using 13 microsatellites revealed that SSH09a was restricted to the WIO while SSH09b and SSH09c were almost exclusively in the TSP and SEP. More surprisingly, each SSH split into two to three genetically differentiated clusters, found in sympatry at the reef scale, leading to a pattern of nested hierarchical levels (PSH > SSH > cluster), each level hiding highly differentiated genetic groups. Thus, rather than structured populations within a single species, these three SSHs, and even the eight clusters, likely represent distinct genetic lineages engaged in a speciation process or real species. The issue is now to understand which hierarchical level (SSH, cluster, or even below) corresponds to the species one. Several hypotheses are discussed on the processes leading to this pattern of mixed clusters in sympatry, evoking formation of reproductive barriers, either by allopatric speciation or habitat selection.

Integrative taxonomy of the ornamental 'peppermint' shrimp public market and population genetics of Lysmata boggessi, the most heavily traded species worldwide

The ornamental trade is a worldwide industry worth >15 billion USD with a problem of rampant product misidentification. Minimizing misidentification is critical in the face of overexploitation of species in the trade. We surveyed the peppermint shrimp ornamental marketplace in the southeastern USA, the most intense market for peppermint shrimps worldwide, to characterize the composition of species in the trade, reveal the extent of misidentification, and describe the population genetics of the true target species. Shrimps were bought from aquarium shops in FL, GA, SC, and NC. We demonstrated, contrary to popular belief (information from dealers), that the most heavily traded species in the market was Lysmata boggessi, an endemic species to the eastern Gulf of Mexico, and not Lysmata wurdemanni. Importantly, only when color pattern or genetic markers in conjunction with morphological traits were employed, was it was possible to unequivocally identify L. boggessi as the only species in the trade. The intensity of the market for peppermint shrimps in the USA has led to L. boggessi being the most traded species worldwide. Misidentification in the shrimp aquarium trade is accidental and involuntary, and is explained by remarkable similarity among congeneric species. Using sequences of the 16S-mt-DNA marker, we found no indication of population genetic structure in the endemic L. boggessi across  550 km of linear coast. Therefore, this species can be considered genetically homogeneous and a single fished stock. Still, we argue in favor of additional studies using more powerful markers (e.g., SNPs) capable of revealing genetic structure at a finer spatial-scale. Our results will help advance management and conservation policies in this lucrative yet understudied fishery. Future studies of other ornamental fisheries will benefit from using an integrative taxonomic approach, as we demonstrate here.