Novel transcriptome resources for three scleractinian coral species from the Indo-Pacific

Synergistic response to climate stressors in coral is associated with genotypic variation in baseline expression

As environments are rapidly reshaped due to climate change, phenotypic plasticity plays an important role in the ability of organisms to persist and is considered an especially important acclimatization mechanism for long-lived sessile organisms such as reef-building corals. Often, this ability of a single genotype to display multiple phenotypes depending on the environment is modulated by changes in gene expression, which can vary in response to environmental changes via two mechanisms: baseline expression and expression plasticity. We used transcriptome-wide expression profiling of eleven genotypes of common-gardened Acropora cervicornis to explore genotypic variation in the expression response to thermal and acidification stress, both individually and in combination. We show that the combination of these two stressors elicits a synergistic gene expression response, and that both baseline expression and expression plasticity in response to stress show genotypic variation. Additionally, we demonstrate that frontloading of a large module of coexpressed genes is associated with greater retention of algal symbionts under combined stress. These results illustrate that variation in the gene expression response of individuals to climate change stressors can persist even when individuals have shared environmental histories, affecting their performance under future climate change scenarios.

Meta-organism gene expression reveals that the impact of nitrate enrichment on coral larvae is mediated by their associated Symbiodiniaceae and prokaryotic assemblages

BACKGROUND

Coral meta-organisms consist of the coral, and its associated Symbiodiniaceae (dinoflagellate algae), bacteria, and other microbes. Corals can acquire photosynthates from Symbiodiniaceae, whilst Symbiodiniaceae uses metabolites from corals. Prokaryotic microbes provide Symbiodiniaceae with nutrients and support the resilience of corals as meta-organisms. Eutrophication is a major cause of coral reef degradation; however, its effects on the transcriptomic response of coral meta-organisms remain unclear, particularly for prokaryotic microbes associated with corals in the larval stage. To understand acclimation of the coral meta-organism to elevated nitrate conditions, we analyzed the physiological and transcriptomic responses of Pocillopora damicornis larvae, an ecologically important scleractinian coral, after 5 days of exposure to elevated nitrate levels (5, 10, 20, and 40 µM).

RESULTS

The major differentially expressed transcripts in coral, Symbiodiniaceae, and prokaryotic microbes included those related to development, stress response, and transport. The development of Symbiodiniaceae was not affected in the 5 and 20 µM groups but was downregulated in the 10 and 40 µM groups. In contrast, prokaryotic microbe development was upregulated in the 10 and 40 µM groups and downregulated in the 5 and 20 µM groups. Meanwhile, coral larval development was less downregulated in the 10 and 40 µM groups than in the 5 and 20 µM groups. In addition, multiple larval, Symbiodiniaceae, and prokaryotic transcripts were significantly correlated with each other. The core transcripts in correlation networks were related to development, nutrient metabolism, and transport. A generalized linear mixed model, using least absolute shrinkage and selection operator, demonstrated that the Symbiodiniaceae could both benefit and cost coral larval development. Furthermore, the most significantly correlated prokaryotic transcripts maintained negative correlations with the physiological functions of Symbiodiniaceae.

CONCLUSIONS

Results suggested that Symbiodiniaceae tended to retain more nutrients under elevated nitrate concentrations, thereby shifting the coral-algal association from mutualism towards parasitism. Prokaryotic microbes provided Symbiodiniaceae with essential nutrients and may control Symbiodiniaceae growth through competition, whereby prokaryotes can also restore coral larval development inhibited by Symbiodiniaceae overgrowth. Video Abstract.

Caryophylliids (Anthozoa, Scleractinia) and mitochondrial gene order: insights from mitochondrial and nuclear phylogenomics

Molecularly, the family Caryophylliidae is polyphyletic and different sets of genetic data converge towards a consensus that a taxonomic review of this family is necessary. Overall, the order of genes in the mitochondrial genome (mitogenome) together with DNA sequences have been used to successfully untangle evolutionary relationships in several groups of organisms. Published mitogenomes of two caryophylliid genera (Desmophyllum and Solenosmilia) present a transposition of the gene block containing cob, nad2, and nad6, which is located between nad5 5' exon and trnW, while that of Polycyathus chaishanensis presents the same gene order as the majority of scleractinian corals. In molecular-based evolutionary reconstructions, caryophylliids that have the mitochondrial gene rearrangement were recovered as a monophyletic lineage ("true" caryophylliids), while members of the genus Polycyathus were placed in a different position. In this study, additional mitogenomes of this family were assembled and included in evolutionary reconstructions of Scleractinia in order to improve our understanding on whether the mitogenome gene rearrangement is limited to and, therefore, could be a synapomorphy of the actual members of Caryophylliidae. Specimens of Caryophyllia scobinosa, Premocyathus sp., Heterocyathus sulcatus, and Trochocyathus caryophylloides, as well as Desmophyllum pertusum and Solenosmilia variabilis from the Southwest Atlantic were sequenced using Illumina platforms. Then, mitochondrial genomes were assembled and annotated, and nuclear datasets were recovered in-silico from assembled contigs using a previously published set of baits. Evolutionary reconstructions were performed using mitochondrial and nuclear datasets and based on Maximum Likelihood and Bayesian Inference. Obtained mitogenomes are circular and range between 15,816 and 18,225 bp in size and from 30.76% to 36.63% in GC content. The gene rearrangement is only seen in C. scobinosa, D. pertusum, Premocyathus sp., and S. variabilis, which were recovered as a monophyletic clade in both mitochondrial and nuclear phylogenies. On the other hand, the "caryophylliids" with the canonical mitogenome gene order were not recovered within this clade. Differences in features of the skeleton of "true" caryophylliids in comparison to traditional members of the family were observed and offer further support that the gene rearrangement might be seen as a synapomorphy of family Caryophylliidae.

Transcriptome-Wide Comparisons and Virulence Gene Polymorphisms of Host-Associated Genotypes of the Cnidarian Parasite Ceratonova shasta in Salmonids

Ceratonova shasta is an important myxozoan pathogen affecting the health of salmonid fishes in the Pacific Northwest of North America. Ceratonova shasta exists as a complex of host-specific genotypes, some with low to moderate virulence, and one that causes a profound, lethal infection in susceptible hosts. High throughput sequencing methods are powerful tools for discovering the genetic basis of these host/virulence differences, but deep sequencing of myxozoans has been challenging due to extremely fast molecular evolution of this group, yielding strongly divergent sequences that are difficult to identify, and unavoidable host contamination. We designed and optimized different bioinformatic pipelines to address these challenges. We obtained a unique set of comprehensive, host-free myxozoan RNA-seq data from C. shasta genotypes of varying virulence from different salmonid hosts. Analyses of transcriptome-wide genetic distances and maximum likelihood multigene phylogenies elucidated the evolutionary relationship between lineages and demonstrated the limited resolution of the established Internal Transcribed Spacer marker for C. shasta genotype identification, as this marker fails to differentiate between biologically distinct genotype II lineages from coho salmon and rainbow trout. We further analyzed the data sets based on polymorphisms in two gene groups related to virulence: cell migration and proteolytic enzymes including their inhibitors. The developed single-nucleotide polymorphism-calling pipeline identified polymorphisms between genotypes and demonstrated that variations in both motility and protease genes were associated with different levels of virulence of C. shasta in its salmonid hosts. The prospective use of proteolytic enzymes as promising candidates for targeted interventions against myxozoans in aquaculture is discussed. We developed host-free transcriptomes of a myxozoan model organism from strains that exhibited different degrees of virulence, as a unique source of data that will foster functional gene analyses and serve as a base for the development of potential therapeutics for efficient control of these parasites.

Maternal effects in gene expression of interspecific coral hybrids

Maternal effects have been well documented for offspring morphology and life history traits in plants and terrestrial animals, yet little is known about maternal effects in corals. Further, few studies have explored maternal effects in gene expression. In a previous study, F1 interspecific hybrid and purebred larvae of the coral species Acropora tenuis and Acropora loripes were settled and exposed to ambient or elevated temperature and pCO₂ conditions for 7 months. At this stage, the hybrid coral recruits from both ocean conditions exhibited strong maternal effects in several fitness traits. We conducted RNA-sequencing on these corals and showed that gene expression of the hybrid Acropora also exhibited clear maternal effects. Only 40 genes were differentially expressed between hybrids and their maternal progenitor. In contrast, ~2000 differentially expressed genes were observed between hybrids and their paternal progenitors, and between the reciprocal F1 hybrids. These results indicate that maternal effects in coral gene expression can be long-lasting. Unlike findings from most short-term stress experiments in corals, no genes were differentially expressed in the hybrid nor purebred offspring after seven months of exposure to elevated temperature and pCO₂ conditions.

Global gene expression patterns in Porites white patch syndrome: Disentangling symbiont loss from the thermal stress response in reef-building coral

The mechanisms resulting in the breakdown of the coral symbiosis once the process of bleaching has been initiated remain unclear. Distinguishing the process of symbiont loss from the thermal stress response may shed light on the cellular and molecular pathways involved in each process. This study examined physiological changes and global gene expression patterns associated with white patch syndrome (WPS) in Porites lobata, which manifests in localized bleaching independent of thermal stress. In addition, a meta-analysis of global gene expression studies in other corals and anemones was used to contrast differential regulation as a result of disease and thermal stress from patterns correlated with symbiotic state. Symbiont density, chlorophyll a content, holobiont productivity, instant calcification rate, and total host protein content were uniformly reduced in WPS relative to healthy tissue. While expression patterns associated with WPS were secondary to fixed effects of source colony, specific functional enrichments combined with a lack of immune regulation suggest that the viral infection putatively giving rise to this condition affects symbiont rather than host cells. Expression in response to WPS also clustered independently of patterns in white syndrome impacted A. hyacinthus, further supporting a distinct aetiology of this syndrome. Expression patterns in WPS-affected tissues were significantly correlated with prior studies that examined short-term thermal stress responses independent of symbiotic state, suggesting that the majority of expression changes reflect a nonspecific stress response. Across studies, the magnitude and direction of expression change among particular functional enrichments suggests unique responses to stressor duration and highlights distinct responses to bleaching in an anemone model.

Investigating the roles of transforming growth factor-beta in immune response of Orbicella faveolata, a scleractinian coral

Symbiotic relationships range from parasitic to mutualistic, yet all endosymbionts face similar challenges, including evasion of host immunity. Many symbiotic organisms have evolved similar mechanisms to face these challenges, including manipulation of the host's transforming growth factor-beta (TGFβ) pathway. Here we investigate the TGFβ pathway in scelaractinian corals which are dependent on symbioses with dinoflagellates from the family Symbiodiniaceae. Using the Caribbean coral, Orbicella faveolata, we explore the effects of enhancement and inhibition of the TGFβ pathway on host gene expression. Following transcriptomic analyses, we demonstrated limited effects of pathway manipulation in absence of immune stimulation. However, manipulation of the TGFβ pathway significantly affects the subsequent ability of host corals to mount an immune response. Enhancement of the TGFβ pathway eliminates transcriptomic signatures of host coral immune response, while inhibition of the pathway maintains the response. This is, to our knowledge, the first evidence of an immunomodulatory role for TGFβ in a scelaractinian coral. These findings suggest variation in TGFβ signaling may have implications in the face of increasing disease prevelance. Our results suggest that the TGFβ pathway can modulate tradeoffs between symbiosis and immunity. Further study of links between symbiosis, TGFβ, and immunity is needed to better understand the ecological implications of these findings.

The Transcriptome Landscape of Walnut Interspecies Hybrid (Juglans hindsii × Juglans regia) and Regulation of Cambial Activity in Relation to Grafting

Walnuts (Juglans, Juglandaceae) are known throughout the world as economically important trees that provide fat, protein, vitamins, and minerals as a food source, and produce high-quality timber. We have amended the purpose section to say "However," the omics resources are limited, which hampered the elucidation of molecular mechanisms resulting in their economically important traits (such as yield, fertility alternation, oil synthesis, and wood formation). To enrich the omics database of walnut, there is great need for analyses of its genomic and transcriptomic characteristics. In this study, we reported for the first time of the transcriptome landscape of six important organs or tissues in walnut interspecies hybrid using next-generation sequencing technology. Over 338 million clean reads were obtained. This yielded 74,072 unigenes with an average length of 782.71 bp. To develop an understanding of gene functions and regulatory pathways, 66,355 of the unigenes were identified as homologs of annotated genes and classified into three general categories with 61 functional subcategories. 2,288 out of 2,549 unmapped unigenes had at least one BLAST hit against the public databases. A total of 1,237 transcription factor-encoding genes (TFs) and 2,297 tissue-specific unigenes were identified. Interestingly, in the new shoot between an adult seedling and a grafted tree, the expression of 9,494 unigenes were significantly different, among which 4,388 were up-regulated and 5,106 were down-regulated. Of these, 195, 177, 232, 75, 114, and 68 unigenes were related to transcription factors, cell wall, defense response, transport, plant hormone biosynthesis, and other cambial activity-related functions, respectively. The obtained sequences and putative functional data constitute a resource for future functional analyses in walnut and other woody plants. These findings will be useful in further studies addressing the molecular mechanisms underlying grafting-related cambial activity.

Development of a transcriptomic database for 14 species of scleractinian corals

BACKGROUND

Scleractinian corals are important reef builders, but around the world they are under the threat of global climate change as well as local stressors. Molecular resources are critical for understanding a species' stress responses and resilience to the changing environment, but such resources are unavailable for most scleractinian corals, especially those distributed in the South China Sea. We therefore aimed to provide transcriptome resources for 14 common species, including a few structure forming species, in the South China Sea.

DESCRIPTION

We sequenced the transcriptome of 14 species of scleractinian corals using high-throughput RNA-seq and conducted de novo assembly. For each species, we produced 7.4 to 12.0 gigabases of reads, and assembled them into 271 to 762 thousand contigs with a N50 value of 629 to 1427 bp. These contigs included 66 to 114 thousand unigenes with a predicted open reading frame, and 74.3 to 80.5% of the unigenes were functionally annotated. In the azooxanthelate species Tubastraea coccinea, 41.5% of the unigenes had at least a best-hit sequence from corals. In the other thirteen species, 20.2 to 48.9% of the annotated unigenes had best-hit sequences from corals, and 28.3 to 51.6% from symbiotic algae belonging to the family Symbiodinaceae. With these resources, we developed a transcriptome database (CoralTBase) which features online BLAST and keyword search for unigenes/functional terms through a user friendly Internet interface.

SHORT CONCLUSION

We developed comprehensive transcriptome resources for 14 species of scleractinian corals and constructed a publicly accessible database ( www.comp.hkbu.edu.hk/~db/CoralTBase ). CoralTBase will facilitate not only functional studies using these corals to understand the molecular basis of stress responses and adaptation, but also comparative transcriptomic studies with other species of corals and more distantly related cnidarians.

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.

Molecular cloning and characterization of calmodulin-like protein CaLP from the Scleractinian coral Galaxea astreata

Optimal temperature and light are both necessary conditions for coral survival. Light enhances calcification, and thermal stress disrupts Ca2+ homeostasis. As calcium is involved in many important metabolic activities, in this study, we cloned the calmodulin-like protein (CaLP) gene of one of the scleractinian corals, Galaxea astreata. We also detected the relative mRNA expression levels of gaCaLP using the calcium channel blocker verapamil and CaCl2 treatment under conditions of light and dark, and compared expression levels under controlled temperature conditions. Full-length gaCaLP cDNA comprised 1290 nucleotides and contained 498 bp open reading frame that encoded a protein with 165 amino acids. With CaCl2, expression levels of gaCaLP only increased in the presence of light, suggesting that light may be a restrictive factor in CaLP expression when sufficient calcium is available in the environment. In addition, after verapami treatment, we noted that a down regulation of gaCaLP, suggesting that the expression of CaLP is closely related to extracellular Ca2+ influx. Under temperature stress at both high (30 °C) and low (20 °C) temperatures, expression levels of gaCaLP showed an initial increase, followed by a decreasing trend as treatment progressed. Expression levels reached their maximum value at 24 h. This result showed that CaLP participated in a temperature stress response, and Ca2+ homeostasis was disrupted during stress. The findings of the present study will help determine the function and regulatory mechanisms of gaCaLP.

De novo transcriptome assembly of the coral Agaricia lamarcki (Lamarck's sheet coral) from mesophotic depth in southwest Puerto Rico

The plating coral, Agaricia lamarcki is a widely distributed species inhabiting reefs across the Caribbean basin and Florida. This species is of interest since it is considered a depth-generalist, found from 10 to 70 m. Given the scope of contemporary studies on this coral's population dynamics and physiology, as well as, the potential of mesophotic reefs to be refuge habitats for deteriorated shallow water reefs, we present the first de novo transcriptome assembly of an important mesophotic coral. Using next-generation paired-end sequencing (Illumina Hiseq4000; 2 × 150 bp), we obtained a total of 82,506,058 raw reads. The novel transcriptome assembly strategy included the recently developed National Center for Genome Analysis Support de novo transcriptome assembly pipeline. Assembly produced a total of 101,322 biologically true, non-redundant transcripts with an average contig length of 959 and N50 of 1830. EvidentialGene and TransDecoder were used to identify open reading frames (ORFs) with homology insight provided by the UniProtKb and PFAM databases. ORF prediction resulted in 38,517 putative ORFs of which 12,107 ORFs were annotated as genes dealing with molecular function, 1266 with biological processes and 416 with cellular components.

Novel transcriptome resources for three scleractinian coral species from the Indo-Pacific

Transcriptomic resources for coral species can provide insight into coral evolutionary history and stress-response physiology. Goniopora columna, Galaxea astreata, and Galaxea acrhelia are scleractinian corals of the Indo-Pacific, representing a diversity of morphologies and life-history traits. G. columna and G. astreata are common and cosmopolitan, while G. acrhelia is largely restricted to the coral triangle and Great Barrier Reef. Reference transcriptomes for these species were assembled from replicate colony fragments exposed to elevated (31°C) and ambient (27°C) temperatures. Trinity was used to create de novo assemblies for each species from 92–102 million raw Illumina Hiseq 2 × 150 bp reads. Host-specific assemblies contained 65 460–72 405 contigs, representing 26 693–37 894 isogroups (∼genes) with an average N50 of 2254. Gene name and/or gene ontology annotations were possible for 58% of isogroups on average. Transcriptomes contained 93.1–94.3% of EuKaryotic Orthologous Groups comprising the core eukaryotic gene set, and 89.98–91.92% of the single-copy metazoan core gene set orthologs were complete, indicating fairly comprehensive assemblies. This work expands the complement of transcriptomic resources available for scleractinian coral species, including the first reference for a representative of Goniopora spp. as well as species with novel morphology.