The Role of DNA Methylation in Genome Defense in Cnidaria and Other Invertebrates.

Considerable attention has recently been focused on the potential involvement of DNA methylation in regulating gene expression in cnidarians. Much of this work has been centered on corals, in the context of changes in methylation perhaps facilitating adaptation to higher seawater temperatures and other stressful conditions. Although first proposed more than 30 years ago, the possibility that DNA methylation systems function in protecting animal genomes against the harmful effects of transposon activity has largely been ignored since that time. Here, we show that transposons are specifically targeted by the DNA methylation system in cnidarians, and that the youngest transposons (i.e., those most likely to be active) are most highly methylated. Transposons in longer and highly active genes were preferentially methylated and, as transposons aged, methylation levels declined, reducing the potentially harmful side effects of CpG methylation. In Cnidaria and a range of other invertebrates, correlation between the overall extent of methylation and transposon content was strongly supported. Present transposon burden is the dominant factor in determining overall level of genomic methylation in a range of animals that diverged in or before the early Cambrian, suggesting that genome defense represents the ancestral role of CpG methylation.

Morphological stasis masks ecologically divergent coral species on tropical reefs.

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

Conservation and turnover of miRNAs and their highly complementary targets in early branching animals.

MicroRNAs (miRNAs) are crucial post-transcriptional regulators that have been extensively studied in Bilateria, a group comprising the majority of extant animals, where more than 30 conserved miRNA families have been identified. By contrast, bilaterian miRNA targets are largely not conserved. Cnidaria is the sister group to Bilateria and thus provides a unique opportunity for comparative studies. Strikingly, like their plant counterparts, cnidarian miRNAs have been shown to predominantly have highly complementary targets leading to transcript cleavage by Argonaute proteins. Here, we assess the conservation of miRNAs and their targets by small RNA sequencing followed by miRNA target prediction in eight species of Anthozoa (sea anemones and corals), the earliest-branching cnidarian class. We uncover dozens of novel miRNAs but only a few conserved ones. Further, given their high complementarity, we were able to computationally identify miRNA targets in each species. Besides evidence for conservation of specific miRNA target sites, which are maintained between sea anemones and stony corals across 500 Myr of evolution, we also find indications for convergent evolution of target regulation by different miRNAs. Our data indicate that cnidarians have only few conserved miRNAs and corresponding targets, despite their high complementarity, suggesting a high evolutionary turnover.

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.

A genomic view of the reef-building coral Porites lutea and its microbial symbionts.

Corals and the reef ecosystems that they support are in global decline due to increasing anthropogenic pressures such as climate change1. However, effective reef conservation strategies are hampered by a limited mechanistic understanding of coral biology and the functional roles of the diverse microbial communities that underpin coral health2,3. Here, we present an integrated genomic characterization of the coral species Porites lutea and its microbial partners. High-quality genomes were recovered from P. lutea, as well as a metagenome-assembled Cladocopium C15 (the dinoflagellate symbiont) and 52 bacterial and archaeal populations. Comparative genomic analysis revealed that many of the bacterial and archaeal genomes encode motifs that may be involved in maintaining association with the coral host and in supplying fixed carbon, B-vitamins and amino acids to their eukaryotic partners. Furthermore, mechanisms for ammonia, urea, nitrate, dimethylsulfoniopropionate and taurine transformation were identified that interlink members of the holobiont and may be important for nutrient acquisition and retention in oligotrophic waters. Our findings demonstrate the critical and diverse roles that microorganisms play within the coral holobiont and underscore the need to consider all of the components of the holobiont if we are to effectively inform reef conservation strategies.

Ocean acidification at a coastal CO2 vent induces expression of stress-related transcripts and transposable elements in the sea anemone Anemonia viridis.

Ocean acidification threatens to disrupt interactions between organisms throughout marine ecosystems. The diversity of reef-building organisms decreases as seawater CO2 increases along natural gradients, yet soft-bodied animals, such as sea anemones, are often resilient. We sequenced the polyA-enriched transcriptome of adult sea anemone Anemonia viridis and its dinoflagellate symbiont sampled along a natural CO2 gradient in Italy to assess stress levels in these organisms. We found that about 3.1% of the anemone transcripts, but <1% of the Symbiodinium sp. transcripts were differentially expressed. Processes enriched at high seawater CO2 were linked to cellular stress and inflammation, including significant up-regulation of protective cellular functions and down-regulation of metabolic pathways. Transposable elements were differentially expressed at high seawater CO2, with an extreme up-regulation (> 100-fold) of the BEL-family of long terminal repeat retrotransposons. Seawater acidified by CO2 generated a significant stress reaction in A. viridis, but no bleaching was observed and Symbiodinium sp. appeared to be less affected. These observed changes indicate the mechanisms by which A. viridis acclimate to survive chronic exposure to ocean acidification conditions. We conclude that many organisms that are common in acidified conditions may nevertheless incur costs due to hypercapnia and/or lowered carbonate saturation states.

Transcriptomic analyses highlight the likely metabolic consequences of colonization of a cnidarian host by native or non-native Symbiodinium species.

Reef-building corals and some other cnidarians form symbiotic relationships with members of the dinoflagellate family Symbiodinaceae. As Symbiodinaceae is a highly diverse taxon, the physiological interactions between its members and their hosts are assumed to differ between associations. The presence of different symbiont types is known to affect expression levels of specific host genes, but knowledge of the effects on the transcriptome more broadly remains limited. In the present study, transcriptome profiling was conducted on the tropical corallimorpharian, Ricordea yuma, following the establishment of symbiosis with either the 'homologous' symbiont Symbiodinium goreaui (also known as Cladocopium goreaui; ITS2 type C1) or 'heterologous' symbionts (predominantly S. trenchii, which is also known as Durusdinium trenchii; ITS2 type D1a) isolated from a different corallimorpharian host (Rhodactis indosinensis). Transcriptomic analyses showed that genes encoding host glycogen biosynthesis pathway components are more highly induced during colonization by the homologous symbiont than by the heterologous symbiont. Similar patterns were also observed for several other genes thought to facilitate symbiotic nutrient exchange, including those involved in lipid translocation/storage and metabolite transport. The gene expression results presented here imply that colonization by homologous or heterologous Symbiodinium types may have very different metabolic consequences for the Ricordea host, supporting the notion that even though some cnidarians may be able to form novel symbioses after bleaching, the metabolic performance of these may be compromised.This article has an associated First Person interview with the first author of the paper.

Transcriptomic analysis reveals protein homeostasis breakdown in the coral Acropora millepora during hypo-saline stress.

BACKGROUND: Coral reefs can experience salinity fluctuations due to rainfall and runoff; these events can have major impacts on the corals and lead to bleaching and mortality. On the Great Barrier Reef (GBR), low salinity events, which occur during summer seasons and can involve salinity dropping ~ 10 PSU correlate with declines in coral cover, and these events are predicted to increase in frequency and severity under future climate change scenarios. In other marine invertebrates, exposure to low salinity causes increased expression of genes involved in proteolysis, responses to oxidative stress, and membrane transport, but the effects that changes in salinity have on corals have so far received only limited attention. To better understand the coral response to hypo-osmotic stress, here we investigated the transcriptomic response of the coral Acropora millepora in both adult and juvenile life stages to acute (1 h) and more prolonged (24 h) exposure to low salinity. RESULTS: Differential gene expression analysis revealed the involvement of both common and specific response mechanisms in Acropora. The general response to environmental stressors included up-regulation of genes involved in the mitigation of macromolecular and oxidative damage, while up-regulation of genes involved in amino acid metabolism and transport represent specific responses to salinity stress. CONCLUSIONS: This study is the first comprehensive transcriptomic analysis of the coral response to low salinity stress and provides important insights into the likely consequences of heavy rainfall and runoff events on coral reefs.

Finding Nemo's Genes: A chromosome-scale reference assembly of the genome of the orange clownfish Amphiprion percula.

The iconic orange clownfish, Amphiprion percula, is a model organism for studying the ecology and evolution of reef fishes, including patterns of population connectivity, sex change, social organization, habitat selection and adaptation to climate change. Notably, the orange clownfish is the only reef fish for which a complete larval dispersal kernel has been established and was the first fish species for which it was demonstrated that antipredator responses of reef fishes could be impaired by ocean acidification. Despite its importance, molecular resources for this species remain scarce and until now it lacked a reference genome assembly. Here, we present a de novo chromosome-scale assembly of the genome of the orange clownfish Amphiprion percula. We utilized single-molecule real-time sequencing technology from Pacific Biosciences to produce an initial polished assembly comprised of 1,414 contigs, with a contig N50 length of 1.86 Mb. Using Hi-C-based chromatin contact maps, 98% of the genome assembly were placed into 24 chromosomes, resulting in a final assembly of 908.8 Mb in length with contig and scaffold N50s of 3.12 and 38.4 Mb, respectively. This makes it one of the most contiguous and complete fish genome assemblies currently available. The genome was annotated with 26,597 protein-coding genes and contains 96% of the core set of conserved actinopterygian orthologs. The availability of this reference genome assembly as a community resource will further strengthen the role of the orange clownfish as a model species for research on the ecology and evolution of reef fishes.

Expression of the neuropeptides RFamide and LWamide during development of the coral Acropora millepora in relation to settlement and metamorphosis.

Neuropeptides play critical roles in cnidarian development. However, although they are known to play key roles in settlement and metamorphosis, their temporal and spatial developmental expression has not previously been characterized in any coral. We here describe Acropora millepora LWamide and RFamide and their developmental expression from the time of their first appearance, using in situ hybridization and FMRFamide immunohistochemistry. AmRFamide transcripts first appear in the ectoderm toward the oral end of the planula larva following blastopore closure. This oral bias becomes less apparent as the planula develops. The cell bodies of AmRFamide-expressing cells are centrally located in the ectoderm, with narrow projections to the mesoglea and to the cell surface. As the planula approaches settlement, AmRFamide expression disappears and is undetectable in the newly settled polyp. Expressing cells then gradually reappear as the polyp develops, becoming particularly abundant on the tentacles. AmLWamide transcripts first appear in ectodermal cells of the developing planula, with minimal expression at its two ends. The cell bodies of expressing cells lie just above the mesoglea, in a position distinct from those of AmRFamide-expressing cells, and have a narrow projection extending across the ectoderm to its surface. AmLWamide-expressing cells persist for most of the planula stage, disappearing shortly before settlement, but later than AmRFamide-expressing cells. As is the case with AmRFamide, expressing cells are absent from the polyp immediately after settlement, reappearing later on its oral side. AmLWamide expression lags that of AmRFamide in both its disappearance and reappearance. Antibodies to FMRFamide stain cells in a pattern similar to that of the transcripts, but also cells in areas where there is no expression revealed by in situ hybridization, most notably at the aboral end of the planula and in the adult polyp. Adult polyps have numerous staining cells on the tentacles and oral discs, as well as an immunoreactive nerve ring around the mouth. There are scattered staining cells in the coenosarc between polyps and staining cells are abundant in the mesenterial filaments. The above results are discussed in the context of our knowledge of the behavior of coral planulae at the time of their settlement and metamorphosis. Corals are facing multiple environmental threats, and these results both highlight the need for, and bring us a step closer to, a mechanistic understanding of a process that is critical to their survival.

Comparative genomics reveals the distinct evolutionary trajectories of the robust and complex coral lineages.

BACKGROUND: Despite the biological and economic significance of scleractinian reef-building corals, the lack of large molecular datasets for a representative range of species limits understanding of many aspects of their biology. Within the Scleractinia, based on molecular evidence, it is generally recognised that there are two major clades, Complexa and Robusta, but the genomic bases of significant differences between them remain unclear. RESULTS: Draft genome assemblies and annotations were generated for three coral species: Galaxea fascicularis (Complexa), Fungia sp., and Goniastrea aspera (Robusta). Whilst phylogenetic analyses strongly support a deep split between Complexa and Robusta, synteny analyses reveal a high level of gene order conservation between all corals, but not between corals and sea anemones or between sea anemones. HOX-related gene clusters are, however, well preserved across all of these combinations. Differences between species are apparent in the distribution and numbers of protein domains and an apparent correlation between number of HSP20 proteins and stress tolerance. Uniquely amongst animals, a complete histidine biosynthesis pathway is present in robust corals but not in complex corals or sea anemones. This pathway appears to be ancestral, and its retention in the robust coral lineage has important implications for coral nutrition and symbiosis. CONCLUSIONS: The availability of three new coral genomes enabled recognition of a de novo histidine biosynthesis pathway in robust corals which is only the second identified biosynthetic difference between corals. These datasets provide a platform for understanding many aspects of coral biology, particularly the interactions of corals with their endosymbionts.

Symbiodinium genomes reveal adaptive evolution of functions related to coral-dinoflagellate symbiosis.

Symbiosis between dinoflagellates of the genus Symbiodinium and reef-building corals forms the trophic foundation of the world's coral reef ecosystems. Here we present the first draft genome of Symbiodinium goreaui (Clade C, type C1: 1.03 Gbp), one of the most ubiquitous endosymbionts associated with corals, and an improved draft genome of Symbiodinium kawagutii (Clade F, strain CS-156: 1.05 Gbp) to further elucidate genomic signatures of this symbiosis. Comparative analysis of four available Symbiodinium genomes against other dinoflagellate genomes led to the identification of 2460 nuclear gene families (containing 5% of Symbiodinium genes) that show evidence of positive selection, including genes involved in photosynthesis, transmembrane ion transport, synthesis and modification of amino acids and glycoproteins, and stress response. Further, we identify extensive sets of genes for meiosis and response to light stress. These draft genomes provide a foundational resource for advancing our understanding of Symbiodinium biology and the coral-algal symbiosis.

Erratum: Publisher Correction: Symbiodinium genomes reveal adaptive evolution of functions related to coral-dinoflagellate symbiosis.

[This corrects the article DOI: 10.1038/s42003-018-0098-3.].

Epigenome-associated phenotypic acclimatization to ocean acidification in a reef-building coral.

There are increasing concerns that the current rate of climate change might outpace the ability of reef-building corals to adapt to future conditions. Work on model systems has shown that environmentally induced alterations in DNA methylation can lead to phenotypic acclimatization. While DNA methylation has been reported in corals and is thought to associate with phenotypic plasticity, potential mechanisms linked to changes in whole-genome methylation have yet to be elucidated. We show that DNA methylation significantly reduces spurious transcription in the coral Stylophora pistillata. Furthermore, we find that DNA methylation also reduces transcriptional noise by fine-tuning the expression of highly expressed genes. Analysis of DNA methylation patterns of corals subjected to long-term pH stress showed widespread changes in pathways regulating cell cycle and body size. Correspondingly, we found significant increases in cell and polyp sizes that resulted in more porous skeletons, supporting the hypothesis that linear extension rates are maintained under conditions of reduced calcification. These findings suggest an epigenetic component in phenotypic acclimatization that provides corals with an additional mechanism to cope with environmental change.

Elucidating the Small Regulatory RNA Repertoire of the Sea Anemone Anemonia viridis Based on Whole Genome and Small RNA Sequencing.

Cnidarians harbor a variety of small regulatory RNAs that include microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs), but detailed information is limited. Here, we report the identification and expression of novel miRNAs and putative piRNAs, as well as their genomic loci, in the symbiotic sea anemone Anemonia viridis. We generated a draft assembly of the A. viridis genome with putative size of 313 Mb that appeared to be composed of about 36% repeats, including known transposable elements. We detected approximately equal fractions of DNA transposons and retrotransposons. Deep sequencing of small RNA libraries constructed from A. viridis adults sampled at a natural CO2 gradient off Vulcano Island, Italy, identified 70 distinct miRNAs. Eight were homologous to previously reported miRNAs in cnidarians, whereas 62 appeared novel. Nine miRNAs were recognized as differentially expressed along the natural seawater pH gradient. We found a highly abundant and diverse population of piRNAs, with a substantial fraction showing ping-pong signatures. We identified nearly 22% putative piRNAs potentially targeting transposable elements within the A. viridis genome. The A. viridis genome appeared similar in size to that of other hexacorals with a very high divergence of transposable elements resembling that of the sea anemone genus Exaiptasia. The genome encodes and expresses a high number of small regulatory RNAs, which include novel miRNAs and piRNAs. Differentially expressed small RNAs along the seawater pH gradient indicated regulatory gene responses to environmental stressors.

Universal target-enrichment baits for anthozoan (Cnidaria) phylogenomics: New approaches to long-standing problems.

Anthozoans (e.g., corals, anemones) are an ecologically important and diverse group of marine metazoans that occur from shallow to deep waters worldwide. However, our understanding of the evolutionary relationships among the ~7,500 species within this class is hindered by the lack of phylogenetically informative markers that can be reliably sequenced across a diversity of taxa. We designed and tested 16,306 RNA baits to capture 720 ultraconserved element loci and 1,071 exon loci. Library preparation and target enrichment were performed on 33 taxa from all orders within the class Anthozoa. Following Illumina sequencing and Trinity assembly, we recovered 1,774 of 1,791 targeted loci. The mean number of loci recovered from each species was 638 ± 222, with more loci recovered from octocorals (783 ± 138 loci) than hexacorals (475 ± 187 loci). Parsimony informative sites ranged from 26 to 49% for alignments at differing hierarchical taxonomic levels (e.g., Anthozoa, Octocorallia, Hexacorallia). The per cent of variable sites within each of three genera (Acropora, Alcyonium, and Sinularia) for which multiple species were sequenced ranged from 4.7% to 30%. Maximum-likelihood analyses recovered highly resolved trees with topologies matching those supported by other studies, including the monophyly of the order Scleractinia. Our results demonstrate the utility of this target-enrichment approach to resolve phylogenetic relationships from relatively old to recent divergences. Redesigning the baits with improved affinities to capture loci within each subclass will provide a valuable toolset to address systematic questions, further our understanding of the timing of diversifications and help resolve long-standing controversial relationships in the class Anthozoa.

Detecting rare asymmetrically methylated cytosines and decoding methylation patterns in the honeybee genome.

Context-dependent gene expression in eukaryotes is controlled by several mechanisms including cytosine methylation that primarily occurs in the CG dinucleotides (CpGs). However, less frequent non-CpG asymmetric methylation has been found in various cell types, such as mammalian neurons, and recent results suggest that these sites can repress transcription independently of CpG contexts. In addition, an emerging view is that CpG hemimethylation may arise not only from deregulation of cellular processes but also be a standard feature of the methylome. Here, we have applied a novel approach to examine whether asymmetric CpG methylation is present in a sparsely methylated genome of the honeybee, a social insect with a high level of epigenetically driven phenotypic plasticity. By combining strand-specific ultra-deep amplicon sequencing of illustrator genes with whole-genome methylomics and bioinformatics, we show that rare asymmetrically methylated CpGs can be unambiguously detected in the honeybee genome. Additionally, we confirm differential methylation between two phenotypically and reproductively distinct castes, queens and workers, and offer new insight into the heterogeneity of brain methylation patterns. In particular, we challenge the assumption that symmetrical methylation levels reflect symmetry in the underlying methylation patterns and conclude that hemimethylation may occur more frequently than indicated by methylation levels. Finally, we question the validity of a prior study in which most of cytosine methylation in this species was reported to be asymmetric.

Host modification of a bacterial quorum-sensing signal induces a phenotypic switch in bacterial symbionts.

Bacterial communities colonize epithelial surfaces of most animals. Several factors, including the innate immune system, mucus composition, and diet, have been identified as determinants of host-associated bacterial communities. Here we show that the early branching metazoan Hydra is able to modify bacterial quorum-sensing signals. We identified a eukaryotic mechanism that enables Hydra to specifically modify long-chain 3-oxo-homoserine lactones into their 3-hydroxy-HSL counterparts. Expression data revealed that Hydra's main bacterial colonizer, Curvibacter sp., responds differentially to N-(3-hydroxydodecanoyl)-l-homoserine lactone (3OHC12-HSL) and N-(3-oxododecanoyl)-l-homoserine lactone (3OC12-HSL). Investigating the impacts of the different N-acyl-HSLs on host colonization elucidated that 3OHC12-HSL allows and 3OC12-HSL represses host colonization of Curvibacter sp. These results show that an animal manipulates bacterial quorum-sensing signals and that this modification leads to a phenotypic switch in the bacterial colonizers. This mechanism may enable the host to manipulate the gene expression and thereby the behavior of its bacterial colonizers.

A comparative integrated gene-based linkage and locus ordering by linkage disequilibrium map for the Pacific white shrimp, Litopenaeus vannamei.

The Pacific whiteleg shrimp, Litopenaeus vannamei, is the most farmed aquaculture species worldwide with global production exceeding 3 million tonnes annually. Litopenaeus vannamei has been the focus of many selective breeding programs aiming to improve growth and disease resistance. However, these have been based primarily on phenotypic measurements and omit potential gains by integrating genetic selection into existing breeding programs. Such integration of genetic information has been hindered by the limited available genomic resources, background genetic parameters and knowledge on the genetic architecture of commercial traits for L. vannamei. This study describes the development of a comprehensive set of genomic gene-based resources including the identification and validation of 234,452 putative single nucleotide polymorphisms in-silico, of which 8,967 high value SNPs were incorporated into a commercially available Illumina Infinium ShrimpLD-24 v1.0 genotyping array. A framework genetic linkage map was constructed and combined with locus ordering by disequilibrium methodology to generate an integrated genetic map containing 4,817 SNPs, which spanned a total of 4552.5 cM and covered an estimated 98.12% of the genome. These gene-based genomic resources will not only be valuable for identifying regions underlying important L. vannamei traits, but also as a foundational resource in comparative and genome assembly activities.