De novo metatranscriptome assembly and coral gene expression profile of Montipora capitata with growth anomaly.

BACKGROUND: Scleractinian corals are a vital component of coral reef ecosystems, and of significant cultural and economic value worldwide. As anthropogenic and natural stressors are contributing to a global decline of coral reefs, understanding coral health is critical to help preserve these ecosystems. Growth anomaly (GA) is a coral disease that has significant negative impacts on coral biology, yet our understanding of its etiology and pathology is lacking. In this study we used RNA-seq along with de novo metatranscriptome assembly and homology assignment to identify coral genes that are expressed in three distinct coral tissue types: tissue from healthy corals ("healthy"), GA lesion tissue from diseased corals ("GA-affected") and apparently healthy tissue from diseased corals ("GA-unaffected"). We conducted pairwise comparisons of gene expression among these three tissue types to identify genes and pathways that help us to unravel the molecular pathology of this coral disease. RESULTS: The quality-filtered de novo-assembled metatranscriptome contained 76,063 genes, of which 13,643 were identified as putative coral genes. Overall gene expression profiles of coral genes revealed high similarity between healthy tissue samples, in contrast to high variance among diseased samples. This indicates GA has a variety of genetic effects at the colony level, including on seemingly healthy (GA-unaffected) tissue. A total of 105 unique coral genes were found differentially expressed among tissue types. Pairwise comparisons revealed the greatest number of differentially expressed genes between healthy and GA-affected tissue (93 genes), followed by healthy and GA-unaffected tissue (33 genes), and GA-affected and -unaffected tissue (7 genes). The putative function of these genes suggests GA is associated with changes in the activity of genes involved in developmental processes and activation of the immune system. CONCLUSION: This is one of the first transcriptome-level studies to investigate coral GA, and the first metatranscriptome assembly for the M. capitata holobiont. The gene expression data, metatranscriptome assembly and methodology developed through this study represent a significant addition to the molecular information available to further our understanding of this coral disease.

Symbiont type and environmental factors affect transcriptome-wide gene expression in the coral Montipora capitata.

Reef-building corals may harbor genetically distinct lineages of endosymbiotic dinoflagellates in the genus Symbiodinium, which have been shown to affect important colony properties, including growth rates and resilience against environmental stress. However, the molecular processes underlying these differences are not well understood. In this study, we used whole transcriptome sequencing (RNA-seq) to assess gene expression differences between 27 samples of the coral Montipora capitata predominantly hosting two different Symbiodinium types in clades C and D. The samples were further characterized by their origin from two field sites on Hawai'i Island with contrasting environmental conditions. We found that transcriptome-wide gene expression profiles clearly separated by field site first, and symbiont clade second. With 273 differentially expressed genes (DEGs, 1.3% of all host transcripts), symbiont clade had a measurable effect on host gene expression, but the effect of field site proved almost an order of magnitude higher (1,957 DEGs, 9.6%). According to SNP analysis, we found moderate evidence for host genetic differentiation between field sites (F ST = 0.046) and among corals harboring alternative symbiont clades (F ST = 0.036), suggesting that site-related gene expression differences are likely due to a combination of local adaptation and acclimatization to environmental factors. The correlation between host gene expression and symbiont clade may be due to several factors, including host genotype or microhabitat selecting for alternative clades, host physiology responding to different symbionts, or direct modulation of host gene expression by Symbiodinium. However, the magnitude of these effects at the level of transcription was unexpectedly small considering the contribution of symbiont type to holobiont phenotype.