Gene expression response under thermal stress in two Hawaiian corals is dominated by ploidy and genotype.

Transcriptome data are frequently used to investigate coral bleaching; however, the factors controlling gene expression in natural populations of these species are poorly understood. We studied two corals, Montipora capitata and Pocillopora acuta, that inhabit the sheltered Kane'ohe Bay, Hawai'i. M. capitata colonies in the bay are outbreeding diploids, whereas P. acuta is a mixture of clonal diploids and triploids. Populations were sampled from six reefs and subjected to either control (no stress), thermal stress, pH stress, or combined pH and thermal stress treatments. RNA-seq data were generated to test two competing hypotheses: (1) gene expression is largely independent of genotype, reflecting a shared treatment-driven response (TDE) or, (2) genotype dominates gene expression, regardless of treatment (GDE). Our results strongly support the GDE model, even under severe stress. We suggest that post-transcriptional processes (e.g., control of translation, protein turnover) modify the signal from the transcriptome, and may underlie the observed differences in coral bleaching sensitivity via the downstream proteome and metabolome.

Evidence of a putative CO2 delivery system to the chromatophore in the photosynthetic amoeba Paulinella.

The photosynthetic amoeba, Paulinella provides a recent (ca. 120 Mya) example of primary plastid endosymbiosis. Given the extensive data demonstrating host lineage-driven endosymbiont integration, we analysed nuclear genome and transcriptome data to investigate mechanisms that may have evolved in Paulinella micropora KR01 (hereinafter, KR01) to maintain photosynthetic function in the novel organelle, the chromatophore. The chromatophore is of α-cyanobacterial provenance and has undergone massive gene loss due to Muller's ratchet, but still retains genes that encode the ancestral α-carboxysome and the shell carbonic anhydrase, two critical components of the biophysical CO2 concentrating mechanism (CCM) in cyanobacteria. We identified KR01 nuclear genes potentially involved in the CCM that arose via duplication and divergence and are upregulated in response to high light and downregulated under elevated CO2. We speculate that these genes may comprise a novel CO2 delivery system (i.e., a biochemical CCM) to promote the turnover of the RuBisCO carboxylation reaction and counteract photorespiration. We posit that KR01 has an inefficient photorespiratory system that cannot fully recycle the C2 product of RuBisCO oxygenation back to the Calvin-Benson cycle. Nonetheless, both these systems appear to be sufficient to allow Paulinella to persist in environments dominated by faster-growing phototrophs.

Diverse fates of ancient horizontal gene transfers in extremophilic red algae.

Horizontal genetic transfer (HGT) is a common phenomenon in eukaryotic genomes. However, the mechanisms by which HGT-derived genes persist and integrate into other pathways remain unclear. This topic is of significant interest because, over time, the stressors that initially favoured the fixation of HGT may diminish or disappear. Despite this, the foreign genes may continue to exist if they become part of a broader stress response or other pathways. The conventional model suggests that the acquisition of HGT equates to adaptation. However, this model may evolve into more complex interactions between gene products, a concept we refer to as the 'Integrated HGT Model' (IHM). To explore this concept further, we studied specialized HGT-derived genes that encode heavy metal detoxification functions. The recruitment of these genes into other pathways could provide clear examples of IHM. In our study, we exposed two anciently diverged species of polyextremophilic red algae from the Galdieria genus to arsenic and mercury stress in laboratory cultures. We then analysed the transcriptome data using differential and coexpression analysis. Our findings revealed that mercury detoxification follows a 'one gene-one function' model, resulting in an indivisible response. In contrast, the arsH gene in the arsenite response pathway demonstrated a complex pattern of duplication, divergence and potential neofunctionalization, consistent with the IHM. Our research sheds light on the fate and integration of ancient HGTs, providing a novel perspective on the ecology of extremophiles.

Hot springs viruses at Yellowstone National Park have ancient origins and are adapted to thermophilic hosts.

Geothermal springs house unicellular red algae in the class Cyanidiophyceae that dominate the microbial biomass at these sites. Little is known about host-virus interactions in these environments. We analyzed the virus community associated with red algal mats in three neighboring habitats (creek, endolithic, soil) at Lemonade Creek, Yellowstone National Park (YNP), USA. We find that despite proximity, each habitat houses a unique collection of viruses, with the giant viruses, Megaviricetes, dominant in all three. The early branching phylogenetic position of genes encoded on metagenome assembled virus genomes (vMAGs) suggests that the YNP lineages are of ancient origin and not due to multiple invasions from mesophilic habitats. The existence of genomic footprints of adaptation to thermophily in the vMAGs is consistent with this idea. The Cyanidiophyceae at geothermal sites originated ca. 1.5 Bya and are therefore relevant to understanding biotic interactions on the early Earth.

Facultative lifestyle drives diversity of coral algal symbionts.

The photosynthetic symbionts of corals sustain biodiverse reefs in nutrient-poor, tropical waters. Recent genomic data illuminate the evolution of coral symbionts under genome size constraints and suggest that retention of the facultative lifestyle, widespread among these algae, confers a selective advantage when compared with a strict symbiotic existence. We posit that the coral symbiosis is analogous to a 'bioreactor' that selects winner genotypes and allows them to rise to high numbers in a sheltered habitat prior to release by the coral host. Our observations lead to a novel hypothesis, the 'stepping-stone model', which predicts that local adaptation under both the symbiotic and free-living stages, in a stepwise fashion, accelerates coral alga diversity and the origin of endemic strains and species.

Peeling back the layers of coral holobiont multi-omics data.

The integration of multiple 'omics' datasets is a promising avenue for answering many important and challenging questions in biology, particularly those relating to complex ecological systems. Although multi-omics was developed using data from model organisms with significant prior knowledge and resources, its application to non-model organisms, such as coral holobionts, is less clear-cut. We explore, in the emerging rice coral model Montipora capitata, the intersection of holobiont transcriptomic, proteomic, metabolomic, and microbiome amplicon data and investigate how well they correlate under high temperature treatment. Using a typical thermal stress regime, we show that transcriptomic and proteomic data broadly capture the stress response of the coral, whereas the metabolome and microbiome datasets show patterns that likely reflect stochastic and homeostatic processes associated with each sample. These results provide a framework for interpreting multi-omics data generated from non-model systems, particularly those with complex biotic interactions among microbial partners.

Ploidy Variation and Its Implications for Reproduction and Population Dynamics in Two Sympatric Hawaiian Coral Species.

Standing genetic variation is a major driver of fitness and resilience and therefore of fundamental importance for threatened species such as stony corals. We analyzed RNA-seq data generated from 132 Montipora capitata and 119 Pocillopora acuta coral colonies collected from Kane'ohe Bay, O'ahu, Hawai'i. Our goals were to determine the extent of colony genetic variation and to study reproductive strategies in these two sympatric species. Surprisingly, we found that 63% of the P. acuta colonies were triploid, with putative independent origins of the different triploid clades. These corals have spread primarily via asexual reproduction and are descended from a small number of genotypes, whose diploid ancestor invaded the bay. In contrast, all M. capitata colonies are diploid and outbreeding, with almost all colonies genetically distinct. Only two cases of asexual reproduction, likely via fragmentation, were identified in this species. We report two distinct strategies in sympatric coral species that inhabit the largest sheltered body of water in the main Hawaiian Islands. These data highlight divergence in reproductive behavior and genome biology, both of which contribute to coral resilience and persistence.

A k-mer-Based Approach for Phylogenetic Classification of Taxa in Environmental Genomic Data.

In the age of genome sequencing, whole-genome data is readily and frequently generated, leading to a wealth of new information that can be used to advance various fields of research. New approaches, such as alignment-free phylogenetic methods that utilize k-mer-based distance scoring, are becoming increasingly popular given their ability to rapidly generate phylogenetic information from whole-genome data. However, these methods have not yet been tested using environmental data, which often tends to be highly fragmented and incomplete. Here, we compare the results of one alignment-free approach (which utilizes the D2 statistic) to traditional multi-gene maximum likelihood trees in 3 algal groups that have high-quality genome data available. In addition, we simulate lower-quality, fragmented genome data using these algae to test method robustness to genome quality and completeness. Finally, we apply the alignment-free approach to environmental metagenome assembled genome data of unclassified Saccharibacteria and Trebouxiophyte algae, and single-cell amplified data from uncultured marine stramenopiles to demonstrate its utility with real datasets. We find that in all instances, the alignment-free method produces phylogenies that are comparable, and often more informative, than those created using the traditional multi-gene approach. The k-mer-based method performs well even when there are significant missing data that include marker genes traditionally used for tree reconstruction. Our results demonstrate the value of alignment-free approaches for classifying novel, often cryptic or rare, species, that may not be culturable or are difficult to access using single-cell methods, but fill important gaps in the tree of life.

Building consensus around the assessment and interpretation of Symbiodiniaceae diversity.

Within microeukaryotes, genetic variation and functional variation sometimes accumulate more quickly than morphological differences. To understand the evolutionary history and ecology of such lineages, it is key to examine diversity at multiple levels of organization. In the dinoflagellate family Symbiodiniaceae, which can form endosymbioses with cnidarians (e.g., corals, octocorals, sea anemones, jellyfish), other marine invertebrates (e.g., sponges, molluscs, flatworms), and protists (e.g., foraminifera), molecular data have been used extensively over the past three decades to describe phenotypes and to make evolutionary and ecological inferences. Despite advances in Symbiodiniaceae genomics, a lack of consensus among researchers with respect to interpreting genetic data has slowed progress in the field and acted as a barrier to reconciling observations. Here, we identify key challenges regarding the assessment and interpretation of Symbiodiniaceae genetic diversity across three levels: species, populations, and communities. We summarize areas of agreement and highlight techniques and approaches that are broadly accepted. In areas where debate remains, we identify unresolved issues and discuss technologies and approaches that can help to fill knowledge gaps related to genetic and phenotypic diversity. We also discuss ways to stimulate progress, in particular by fostering a more inclusive and collaborative research community. We hope that this perspective will inspire and accelerate coral reef science by serving as a resource to those designing experiments, publishing research, and applying for funding related to Symbiodiniaceae and their symbiotic partnerships.

Algae obscura: The potential of rare species as model systems.

Model organism research has provided invaluable knowledge about foundational biological principles. However, most of these studies have focused on species that are in high abundance, easy to cultivate in the lab, and represent only a small fraction of extant biodiversity. Here, we present three examples of rare algae with unusual features that we refer to as "algae obscura." The Cyanidiophyceae (Rhodophyta), Glaucophyta, and Paulinella (rhizarian) lineages have all transitioned out of obscurity to become models for fundamental evolutionary research. Insights have been gained into the prevalence and importance of eukaryotic horizontal gene transfer, early Earth microbial community dynamics, primary plastid endosymbiosis, and the origin of Archaeplastida. By reviewing the research that has come from the exploration of these organisms, we demonstrate that underappreciated algae have the potential to help us formulate, refine, and substantiate core hypotheses and that such organisms should be considered when establishing future model systems.

Endosymbiotic ratchet accelerates divergence after organelle origin: The Paulinella model for plastid evolution: The Paulinella model for plastid evolution.

We hypothesize that as one of the most consequential events in evolution, primary endosymbiosis accelerates lineage divergence, a process we refer to as the endosymbiotic ratchet. Our proposal is supported by recent work on the photosynthetic amoeba, Paulinella, that underwent primary plastid endosymbiosis about 124 Mya. This amoeba model allows us to explore the early impacts of photosynthetic organelle (plastid) origin on the host lineage. The current data point to a central role for effective population size (Ne ) in accelerating divergence post-endosymbiosis due to limits to dispersal and reproductive isolation that reduce Ne , leading to local adaptation. We posit that isolated populations exploit different strategies and behaviors and assort themselves in non-overlapping niches to minimize competition during the early, rapid evolutionary phase of organelle integration. The endosymbiotic ratchet provides a general framework for interpreting post-endosymbiosis lineage evolution that is driven by disruptive selection and demographic and population shifts. Also see the video abstract here: https://youtu.be/gYXrFM6Zz6Q.

Multiple waves of viral invasions in Symbiodiniaceae algal genomes.

Dinoflagellates from the family Symbiodiniaceae are phototrophic marine protists that engage in symbiosis with diverse hosts. Their large and distinct genomes are characterized by pervasive gene duplication and large-scale retroposition events. However, little is known about the role and scale of horizontal gene transfer (HGT) in the evolution of this algal family. In other dinoflagellates, high levels of HGTs have been observed, linked to major genomic transitions, such as the appearance of a viral-acquired nucleoprotein that originated via HGT from a large DNA algal virus. Previous work showed that Symbiodiniaceae from different hosts are actively infected by viral groups, such as giant DNA viruses and ssRNA viruses, that may play an important role in coral health. Latent viral infections may also occur, whereby viruses could persist in the cytoplasm or integrate into the host genome as a provirus. This hypothesis received experimental support; however, the cellular localization of putative latent viruses and their taxonomic affiliation are still unknown. In addition, despite the finding of viral sequences in some genomes of Symbiodiniaceae, viral origin, taxonomic breadth, and metabolic potential have not been explored. To address these questions, we searched for putative viral-derived proteins in thirteen Symbiodiniaceae genomes. We found fifty-nine candidate viral-derived HGTs that gave rise to twelve phylogenies across ten genomes. We also describe the taxonomic affiliation of these virus-related sequences, their structure, and their genomic context. These results lead us to propose a model to explain the origin and fate of Symbiodiniaceae viral acquisitions.

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.

Loss of key endosymbiont genes may facilitate early host control of the chromatophore in Paulinella.

The primary plastid endosymbiosis (∼124 Mya) that occurred in the heterotrophic amoeba lineage, Paulinella, is at an earlier stage of evolution than in Archaeplastida, and provides an excellent model for studying organelle integration. Using genomic data from photosynthetic Paulinella, we identified a plausible mechanism for the evolution of host control of endosymbiont (termed the chromatophore) biosynthetic pathways and functions. Specifically, random gene loss from the chromatophore and compensation by nuclear-encoded gene copies enables host control of key pathways through a minimal number of evolutionary innovations. These gene losses impact critical enzymatic steps in nucleotide biosynthesis and the more peripheral components of multi-protein DNA replication complexes. Gene retention in the chromatophore likely reflects the need to maintain a specific stoichiometric balance of the encoded products (e.g., involved in DNA replication) rather than redox state, as in the highly reduced plastid genomes of algae and plants.

Development of a portable toolkit to diagnose coral thermal stress.

Coral bleaching, precipitated by the expulsion of the algal symbionts that provide colonies with fixed carbon is a global threat to reef survival. To protect corals from anthropogenic stress, portable tools are needed to detect and diagnose stress syndromes and assess population health prior to extensive bleaching. Here, medical grade Urinalysis strips, used to detect an array of disease markers in humans, were tested on the lab stressed Hawaiian coral species, Montipora capitata (stress resistant) and Pocillopora acuta (stress sensitive), as well as samples from nature that also included Porites compressa. Of the 10 diagnostic reagent tests on these strips, two appear most applicable to corals: ketone and leukocytes. The test strip results from M. capitata were explored using existing transcriptomic data from the same samples and provided evidence of the stress syndromes detected by the strips. We designed a 3D printed smartphone holder and image processing software for field analysis of test strips (TestStripDX) and devised a simple strategy to generate color scores for corals (reflecting extent of bleaching) using a smartphone camera (CoralDX). Our approaches provide field deployable methods, that can be improved in the future (e.g., coral-specific stress test strips) to assess reef health using inexpensive tools and freely available software.

Retrotransposition facilitated the establishment of a primary plastid in the thecate amoeba Paulinella.

The evolution of eukaryotic life was predicated on the development of organelles such as mitochondria and plastids. During this complex process of organellogenesis, the host cell and the engulfed prokaryote became genetically codependent, with the integration of genes from the endosymbiont into the host nuclear genome and subsequent gene loss from the endosymbiont. This process required that horizontally transferred genes become active and properly regulated despite inherent differences in genetic features between donor (endosymbiont) and recipient (host). Although this genetic reorganization is considered critical for early stages of organellogenesis, we have little knowledge about the mechanisms governing this process. The photosynthetic amoeba Paulinella micropora offers a unique opportunity to study early evolutionary events associated with organellogenesis and primary endosymbiosis. This amoeba harbors a "chromatophore," a nascent photosynthetic organelle derived from a relatively recent cyanobacterial association (∼120 million years ago) that is independent of the evolution of primary plastids in plants (initiated ∼1.5 billion years ago). Analysis of the genome and transcriptome of Paulinella revealed that retrotransposition of endosymbiont-derived nuclear genes was critical for their domestication in the host. These retrocopied genes involved in photoprotection in cyanobacteria became expanded gene families and were "rewired," acquiring light-responsive regulatory elements that function in the host. The establishment of host control of endosymbiont-derived genes likely enabled the cell to withstand photo-oxidative stress generated by oxygenic photosynthesis in the nascent organelle. These results provide insights into the genetic mechanisms and evolutionary pressures that facilitated the metabolic integration of the host­endosymbiont association and sustained the evolution of a photosynthetic organelle.

Genome-powered classification of microbial eukaryotes: focus on coral algal symbionts.

Modern microbial taxonomy generally relies on the use of single marker genes or sets of concatenated genes to generate a framework for the delineation and classification of organisms at different taxonomic levels. However, given that DNA is the 'blueprint of life', and hence the ultimate arbiter of taxonomy, classification systems should attempt to use as much of the blueprint as possible to capture a comprehensive phylogenetic signal. Recent analysis of whole-genome sequences from coral reef symbionts (dinoflagellates of the family Symbiodiniaceae) and other microalgal groups has uncovered extensive divergence not recognised by current algal taxonomic approaches. In the era of 'sequence everything', we argue that whole-genome data are pivotal to guide informed taxonomic inference, particularly for microbial eukaryotes.

Hypothesis: Trans-splicing Generates Evolutionary Novelty in the Photosynthetic Amoeba Paulinella.

Plastid primary endosymbiosis has occurred twice, once in the Archaeplastida ancestor and once in the Paulinella (Rhizaria) lineage. Both events precipitated massive evolutionary changes, including the recruitment and activation of genes that are horizontally acquired (HGT) and the redeployment of existing genes and pathways in novel contexts. Here we address the latter aspect in Paulinella micropora KR01 (hereafter, KR01) that has independently evolved spliced leader (SL) trans-splicing (SLTS) of nuclear-derived transcripts. We investigated the role of this process in gene regulation, novel gene origination, and endosymbiont integration. Our analysis shows that 20% of KR01 genes give rise to transcripts with at least one (but in some cases, multiple) sites of SL addition. This process, which often occurs at canonical cis-splicing acceptor sites (internal introns), results in shorter transcripts that may produce 5'-truncated proteins with novel functions. SL-truncated transcripts fall into four categories that may show: (i) altered protein localization, (ii) altered protein function, structure, or regulation, (iii) loss of valid alternative start codons, preventing translation, or (iv) multiple SL addition sites at the 5'-terminus. The SL RNA genes required for SLTS are putatively absent in the heterotrophic sister lineage of photosynthetic Paulinella species. Moreover, a high proportion of transcripts derived from genes of endosymbiotic gene transfer (EGT) and HGT origin contain SL sequences. We hypothesize that truncation of transcripts by SL addition may facilitate the generation and expression of novel gene variants and that SLTS may have enhanced the activation and fixation of foreign genes in the host genome of the photosynthetic lineages, playing a key role in primary endosymbiont integration.

Evidence for a robust photosystem II in the photosynthetic amoeba Paulinella.

Paulinella represents the only known case of an independent primary plastid endosymbiosis, outside Archaeplastida, that occurred c. 120 (million years ago) Ma. These photoautotrophs grow very slowly in replete culture medium with a doubling time of 6-7 d at optimal low light, and are highly sensitive to photodamage under moderate light levels. We used genomic and biophysical methods to investigate the extreme slow growth rate and light sensitivity of Paulinella, which are key to photosymbiont integration. All photosystem II (PSII) genes except psb28-2 and all cytochrome b6 f complex genes except petM and petL are present in Paulinella micropora KR01 (hereafter, KR01). Biophysical measurements of the water oxidation complex, variable chlorophyll fluorescence, and photosynthesis-irradiance curves show no obvious evidence of PSII impairment. Analysis of photoacclimation under high-light suggests that although KR01 can perform charge separation, it lacks photoprotection mechanisms present in cyanobacteria. We hypothesize that Paulinella species are restricted to low light environments because they are deficient in mitigating the formation of reactive oxygen species formed within the photosystems under peak solar intensities. The finding that many photoprotection genes have been lost or transferred to the host-genome during endosymbiont genome reduction, and may lack light-regulation, is consistent with this hypothesis.

Multi-omic characterization of the thermal stress phenome in the stony coral Montipora capitata.

BACKGROUND: Corals, which form the foundation of biodiverse reef ecosystems, are under threat from warming oceans. Reefs provide essential ecological services, including food, income from tourism, nutrient cycling, waste removal, and the absorption of wave energy to mitigate erosion. Here, we studied the coral thermal stress response using network methods to analyze transcriptomic and polar metabolomic data generated from the Hawaiian rice coral Montipora capitata. Coral nubbins were exposed to ambient or thermal stress conditions over a 5-week period, coinciding with a mass spawning event of this species. The major goal of our study was to expand the inventory of thermal stress-related genes and metabolites present in M. capitata and to study gene-metabolite interactions. These interactions provide the foundation for functional or genetic analysis of key coral genes as well as provide potentially diagnostic markers of pre-bleaching stress. A secondary goal of our study was to analyze the accumulation of sex hormones prior to and during mass spawning to understand how thermal stress may impact reproductive success in M. capitata. METHODS: M. capitata was exposed to thermal stress during its spawning cycle over the course of 5 weeks, during which time transcriptomic and polar metabolomic data were collected. We analyzed these data streams individually, and then integrated both data sets using MAGI (Metabolite Annotation and Gene Integration) to investigate molecular transitions and biochemical reactions. RESULTS: Our results reveal the complexity of the thermal stress phenome in M. capitata, which includes many genes involved in redox regulation, biomineralization, and reproduction. The size and number of modules in the gene co-expression networks expanded from the initial stress response to the onset of bleaching. The later stages involved the suppression of metabolite transport by the coral host, including a variety of sodium-coupled transporters and a putative ammonium transporter, possibly as a response to reduction in algal productivity. The gene-metabolite integration data suggest that thermal treatment results in the activation of animal redox stress pathways involved in quenching molecular oxygen to prevent an overabundance of reactive oxygen species. Lastly, evidence that thermal stress affects reproductive activity was provided by the downregulation of CYP-like genes and the irregular production of sex hormones during the mass spawning cycle. Overall, redox regulation and metabolite transport are key components of the coral animal thermal stress phenome. Mass spawning was highly attenuated under thermal stress, suggesting that global climate change may negatively impact reproductive behavior in this species.