Astrangia poculata is a temperate scleractinian coral that exists in facultative symbiosis with the dinoflagellate alga Breviolum psygmophilum across a range spanning the Gulf of Mexico to Cape Cod, Massachusetts. Our previous work on metabolic thermal performance of Virginia (VA) and Rhode Island (RI) populations of A. poculata revealed physiological signatures of cold (RI) and warm (VA) adaptation of these populations to their respective local thermal environments. Here, we used whole-transcriptome sequencing (mRNA-Seq) to evaluate genetic differences and identify potential loci involved in the adaptive signature of VA and RI populations. Sequencing data from 40 A. poculata individuals, including 10 colonies from each population and symbiotic state (VA-white, VA-brown, RI-white, and RI-brown), yielded a total of 1,808 host-associated and 59 algal symbiont-associated single nucleotide polymorphisms (SNPs) post filtration. Fst outlier analysis identified 66 putative high outlier SNPs in the coral host and 4 in the algal symbiont. Differentiation of VA and RI populations in the coral host was driven by putatively adaptive loci, not neutral divergence (Fst = 0.16, p = 0.001 and Fst = 0.002, p = 0.269 for outlier and neutral SNPs respectively). In contrast, we found evidence of neutral population differentiation in B. psygmophilum (Fst = 0.093, p = 0.001). Several putatively adaptive host loci occur on genes previously associated with the coral stress response. In the symbiont, three of four putatively adaptive loci are associated with photosystem proteins. The opposing pattern of neutral differentiation in B. psygmophilum, but not the A. poculata host, reflects the contrasting dynamics of coral host and algal symbiont population connectivity, dispersal, and gene by environment interactions.

Within the siphonous green algal order Bryopsidales, the size and gene arrangement of chloroplast genomes has been examined extensively, while mitochondrial genomes have been mostly overlooked. The recently published mitochondrial genome of Caulerpa lentillifera is large with expanded noncoding DNA, but it remains unclear if this is characteristic of the entire order. Our study aims to evaluate the evolutionary forces shaping organelle genome dynamics in the Bryopsidales based on the C. lentillifera and Ostreobium quekettii mitochondrial genomes. In this study, the mitochondrial genome of O. quekettii was characterised using a combination of long and short read sequencing, and bioinformatic tools for annotation and sequence analyses. We compared the mitochondrial and chloroplast genomes of O. quekettii and C. lentillifera to examine hypotheses related to genome evolution. The O. quekettii mitochondrial genome is the largest green algal mitochondrial genome sequenced (241,739 bp), considerably larger than its chloroplast genome. As with the mtDNA of C. lentillifera, most of this excess size is from the expansion of intergenic DNA and proliferation of introns. Inflated mitochondrial genomes in the Bryopsidales suggest effective population size, recombination and/or mutation rate, influenced by nuclear-encoded proteins, differ between the genomes of mitochondria and chloroplasts, reducing the strength of selection to influence evolution of their mitochondrial genomes.

The role of species' interactions in structuring biological communities remains unclear. Mutualistic symbioses, involving close positive interactions between two distinct organismal lineages, provide an excellent means to explore the roles of both evolutionary and ecological processes in determining how positive interactions affect community structure. In this study, we investigate patterns of co-diversification between fungi and algae for a range of New Zealand lichens at the community, genus, and species levels and explore explanations for possible patterns related to spatial scale and pattern, taxonomic diversity of the lichens considered, and the level sampling replication. We assembled six independent datasets to compare patterns in phylogenetic congruence with varied spatial extent of sampling, taxonomic diversity and level of specimen replication. For each dataset, we used the DNA sequences from the ITS regions of both the fungal and algal genomes from lichen specimens to produce genetic distance matrices. Phylogenetic congruence between fungi and algae was quantified using distance-based redundancy analysis and we used geographic distance matrices in Moran's eigenvector mapping and variance partitioning to evaluate the effects of spatial variation on the quantification of phylogenetic congruence. Phylogenetic congruence was highly significant for all datasets and a large proportion of variance in both algal and fungal genetic distances was explained by partner genetic variation. Spatial variables, primarily at large and intermediate scales, were also important for explaining genetic diversity patterns in all datasets. Interestingly, spatial structuring was stronger for fungal than algal genetic variation. As the spatial extent of the samples increased, so too did the proportion of explained variation that was shared between the spatial variables and the partners' genetic variation. Different lichen taxa showed some variation in their phylogenetic congruence and spatial genetic patterns and where greater sample replication was used, the amount of variation explained by partner genetic variation increased. Our results suggest that the phylogenetic congruence pattern, at least at small spatial scales, is likely due to reciprocal co-adaptation or co-dispersal. However, the detection of these patterns varies among different lichen taxa, across spatial scales and with different levels of sample replication. This work provides insight into the complexities faced in determining how evolutionary and ecological processes may interact to generate diversity in symbiotic association patterns at the population and community levels. Further, it highlights the critical importance of considering sample replication, taxonomic diversity and spatial scale in designing studies of co-diversification.

The genus Plocamium encompasses seaweeds that are widely distributed throughout the world's oceans, with Plocamium brasiliense found along the tropical and subtropical coasts of the Western Atlantic. This wide distribution can lead to structured populations due to environmental differences (e.g., light levels or temperature), restricted gene flow, and the presence of cryptic species. Abiotic variation can also affect gene expression, which consequently leads to differences in the seaweeds protein profile. This study aimed to analyze the genetic and proteomic profiles of P. brasiliense sampled in two geographically distinct sites on the coastline of Rio de Janeiro state, Brazil: Arraial do Cabo (P1) and Búzios (P2). The genetic profiles of macroalgal specimens from these two sites were indistinguishable as assessed by the markers UPA/23S, rbcL, and COI-5P; however, the protein profiles varied significantly between populations from the two sites. At both sites the ribulose-1,5-biphosphate carboxylase/oxygenase was the most abundant protein found in P. brasiliense specimens. The number of phycobiliproteins differed between both sites with the highest numbers being found at P1, possibly due to water depth. The differences in proteomic profiles of the two nearly identical populations of P. brasiliense suggest that environmental parameters such as light availability and desiccation might induce distinct protein expression, probably as a result of the phenotypic plasticity within this population of seaweed.

Microbial community structures of harmful algal bloom (HAB) caused by Heterosigma akashiwo in Geoje were analyzed using the MiSeq platform. To analyze phytoplankton communities without cross-reactivity with predominant bacteria, a new phytoplankton-specific 23S universal primer set was designed by modifying two previously used ones. The new universal primer set turned out to be a useful tool for the analysis of the phytoplankton community; it showed a high specificity for phytoplankton without cross-reactivity to bacterial sequences as well as the wide taxon coverage presenting from prokaryotic cyanobacteria to eukaryotic algae. Next Generation Sequencing (NGS) data generated by two universal primer sets (16S and 23S) provided useful information about the H. akashiwo bloom. According to the 23S universal primer set, proportions of H. akashiwo increased by more than 200-fold as the bloom occurred and its numbers were high enough to detect in control sites. Its operational taxonomic units (OTUs) were detected in the bloom sites at low proportions suggesting that the 16S universal primer set may not be as effective for monitoring harmful algal blooming (HAB) as the 23S universal primer set. In addition, several abundant OTUs in Chlorophyta were not presented by the 16S universal primer set in this study. However, the 16S primer set was useful for detecting decreases in Foraminifera as HAB occurred suggesting that genomic analyses using two universal primer sets would provide more reliable data for understanding microbial community changes by various environmental or ecological events, including HAB. Genomic analyses using two universal primer sets was also useful for determining a correlation between microbial components as HAB occurred. Heterosigma akashiwo was positively correlated with other bloom species, including Karenia mikimotoi, Teleaulax amphioxeia, and bacteria in Verrucomicrobia.

The dissolved organic matter (DOM) released from the cocoolithophores (Chrysotila dentata) was studied in laboratory experiments after co-culturing C. dentata with bacteria. Marinobacter hydrocarbonoclasticus (CA6)-γ-Proteobacteria and Bacillus firmus (CF2) were used to investigate the utilization and processing of the DOM derived from C. dentata, utilizing fluorescence excitation-emission matrix (EEM) combined with parallel factor analysis (EEM-PARAFAC), while measuring algal abundance and photosynthetic parameters. The experimental groups consisted of axenic C. dentata groups, filter cultured with bacteria (CA6 or CF2) groups, C. dentata co-cultured with bacteria (CA6 or CF2) groups and axenic bacteria (CA6 or CF2) groups. We then evaluated the processing of DOM by determining four fluorescence indices. The number of C. dentata cells and the photosynthetic capacity of microalgae were enhanced by CA6 and CF2. The main known fluorophores, including humic-like components and protein-like components, were present in all sample. The protein-like component of algal-bacterial co-cultures was effectively utilized by CA6 and CF2. The humic-like components increased at the end of the culture time for all cultures. Meanwhile, the average fluorescence intensity of protein-like in CA6 co-culture with algae was lower than that in CF2 co-culture with algae over time. On the other hand, the average fluorescence intensity of humic-like in CA6 was higher than CF2. However, the total change in fluorescence in humic-like and protein-like of axenic CF2 cultures was lower than that of CA6. Hence, the ability of CA6 to transform microalgal-derived DOM was superior to that of CF2, and CF2's ability to consume bacterial-derived DOM was superior to that of CA6.

Trichoplax adhaerens, the only known species of Placozoa is likely to be closely related to an early metazoan that preceded branching of Cnidaria and Bilateria. This animal species is surprisingly well adapted to free life in the World Ocean inhabiting tidal costal zones of oceans and seas with warm to moderate temperatures and shallow waters. The genome of T. adhaerens (sp. Grell) includes four nuclear receptors, namely orthologue of RXR (NR2B), HNF4 (NR2A), COUP-TF (NR2F) and ERR (NR3B) that show a high degree of similarity with human orthologues. In the case of RXR, the sequence identity to human RXR alpha reaches 81% in the DNA binding domain and 70% in the ligand binding domain. We show that T. adhaerens RXR (TaRXR) binds 9-cis retinoic acid (9-cis-RA) with high affinity, as well as high specificity and that exposure of T. adhaerens to 9-cis-RA regulates the expression of the putative T. adhaerens orthologue of vertebrate L-malate-NADP+ oxidoreductase (EC 1.1.1.40) which in vertebrates is regulated by a heterodimer of RXR and thyroid hormone receptor. Treatment by 9-cis-RA alters the relative expression profile of T. adhaerens nuclear receptors, suggesting the existence of natural ligands. Keeping with this, algal food composition has a profound effect on T. adhaerens growth and appearance. We show that nanomolar concentrations of 9-cis-RA interfere with T. adhaerens growth response to specific algal food and causes growth arrest. Our results uncover an endocrine-like network of nuclear receptors sensitive to 9-cis-RA in T. adhaerens and support the existence of a ligand-sensitive network of nuclear receptors at the base of metazoan evolution.

It is well-established that there is a hierarchy of susceptibilities amongst coral genera during heat-stress. However, molecular mechanisms governing these differences are still poorly understood. Here we explored if specific corals possessing different morphologies and different susceptibilities to heat stress may manifest varied gene expression patterns. We examined expression patterns of seven genes in the branching corals Stylophora pistillata and Acropora eurystoma and additionally in the massive robust coral, Porites sp. The tested genes are representatives of key cellular processes occurring during heat-stress in Cnidaria: oxidative stress, ER stress, energy metabolism, DNA repair and apoptosis. Varied response to the heat-stress, in terms of visual coral paling, algal maximum quantum yield and host gene expression was evident in the different growth forms. The two branching corals exhibited similar overall responses that differed from that of the massive coral. A. eurystoma that is considered as a susceptible species did not bleach in our experiment, but tissue sloughing was evident at 34 °C. Interestingly, in this species redox regulation genes were up-regulated at the very onset of the thermal challenge. In S. pistillata, bleaching was evident at 34 °C and most of the stress markers were already up-regulated at 32 °C, either remaining highly expressed or decreasing when temperatures reached 34 °C. The massive Porites species displayed severe bleaching at 32 °C but stress marker genes were only significantly elevated at 34 °C. We postulate that by expelling the algal symbionts from Porites tissues, oxidation damages are reduced and stress genes are activated only at a progressed stage. The differential gene expression responses exhibited here can be correlated with the literature well-documented hierarchy of susceptibilities amongst coral morphologies and genera in Eilat's coral reef.

In many freshwater habitats, green algae form intracellular symbioses with a variety of heterotrophic host taxa including several species of freshwater sponge. These sponges perform important ecological roles in their habitats, and the poriferan:green algae partnerships offers unique opportunities to study the evolutionary origins and ecological persistence of endosymbioses. We examined the association between Ephydatia muelleri and its chlorophyte partner to identify features of host cellular and genetic responses to the presence of intracellular algal partners. Chlorella-like green algal symbionts were isolated from field-collected adult E. muelleri tissue harboring algae. The sponge-derived algae were successfully cultured and subsequently used to reinfect aposymbiotic E. muelleri tissue. We used confocal microscopy to follow the fate of the sponge-derived algae after inoculating algae-free E. muelleri grown from gemmules to show temporal patterns of symbiont location within host tissue. We also infected aposymbiotic E. muelleri with sponge-derived algae, and performed RNASeq to study differential expression patterns in the host relative to symbiotic states. We compare and contrast our findings with work in other systems (e.g., endosymbiotic Hydra) to explore possible conserved evolutionary pathways that may lead to stable mutualistic endosymbioses. Our work demonstrates that freshwater sponges offer many tractable qualities to study features of intracellular occupancy and thus meet criteria desired for a model system.

Dinoflagellate cysts (i.e., dinocysts) are biologically and ecologically important as they can help dinoflagellate species survive harsh environments, facilitate their dispersal and serve as seeds for harmful algal blooms. In addition, dinocysts derived from some species can produce more toxins than vegetative forms, largely affecting species through their food webs and even human health. Consequently, accurate identification of dinocysts represents the first crucial step in many ecological studies. As dinocysts have limited or even no available taxonomic keys, molecular methods have become the first priority for dinocyst identification. However, molecular identification of dinocysts, particularly when using single cells, poses technical challenges. The most serious is the low success rate of PCR, especially for heterotrophic species.

In this study we compared genotypes of zoantharian host-associating algal symbionts among Palythoa species, which are among the dominant benthic reef organisms in the Ryukyu Archipelago, Japan, and evaluated Symbiodiniaceae diversities of closely related congeneric Palythoa species. We targeted a species complex of the zoantharian genus Palythoa (P. tuberculosa, P. sp. yoron, P. mutuki) living among different microhabitats in a narrow reef area of Tokunoshima Island. For phylogenetic analyses, we used two DNA marker regions; nuclear internal transcribed spacer (ITS) and plastid mini-circle non-coding region (psbAncr), both of which have previously been used to determine Symbiodiniaceae genotypes of zoantharian species. Our results showed that all Palythoa species hosted symbionts of the genus Cladocopium, with genotypic compositions of this genus showing some variations among the three different Palythoa species. Additionally, we found that the Cladocopium genotypic composition was statistically different among Palythoa species, and among P. tuberculosa specimens in different microhabitats. Our results suggest that ecological divergence among these three Palythoa species may be related to differing Symbiodiniaceae diversities that may in turn contribute to eco-physiological adaptation into different microhabitats on coral reefs.

Construction of hydropower stations has been an important approach to meet China's increasing power demand, but the impact of construction of hydropower stations on river microbiota is not fully understood. To evaluate this, the microbial composition from 18 sampling sites in the downstream of Jinsha River of China, upstream and downstream of two completed and two under-construction hydropower stations, were analyzed using high-throughput 16S rRNA gene sequencing. Three independent samples from each site were analyzed. A total of 18,683 OTUs from 1,350 genera were identified at 97% sequence similarity. Our results showed that the completion of hydropower stations would significantly increase the relative abundances of Acidobacteria, Chlorobi, Chloroflexi, Cyanobacteria, Nitrospirae, and Planctomycetes, especially the relative abundance of Synechococcus dOTUs and thus increase the risk of algal blooms. PCA based on all KEGG pathways and the significantly different KEGG pathways showed the predicted metabolic characteristics of the water microbiota by PICRUSt in the activated hydropower station group were significant difference to the other groups. Results from canonical correspondence analysis showed that water temperature and dissolved oxygen had significant effects on microbiota composition. These results are important for assessing the impact of hydropower stations on river microbiota and their potential environmental risks.

Current research posits that all multicellular organisms live in symbioses with associated microorganisms and form so-called metaorganisms or holobionts. Cnidarian metaorganisms are of specific interest given that stony corals provide the foundation of the globally threatened coral reef ecosystems. To gain first insight into viruses associated with the coral model system Aiptasia (sensu Exaiptasia pallida), we analyzed an existing RNA-Seq dataset of aposymbiotic, partially populated, and fully symbiotic Aiptasia CC7 anemones with Symbiodinium. Our approach included the selective removal of anemone host and algal endosymbiont sequences and subsequent microbial sequence annotation. Of a total of 297 million raw sequence reads, 8.6 million (∼3%) remained after host and endosymbiont sequence removal. Of these, 3,293 sequences could be assigned as of viral origin. Taxonomic annotation of these sequences suggests that Aiptasia is associated with a diverse viral community, comprising 116 viral taxa covering 40 families. The viral assemblage was dominated by viruses from the families Herpesviridae (12.00%), Partitiviridae (9.93%), and Picornaviridae (9.87%). Despite an overall stable viral assemblage, we found that some viral taxa exhibited significant changes in their relative abundance when Aiptasia engaged in a symbiotic relationship with Symbiodinium. Elucidation of viral taxa consistently present across all conditions revealed a core virome of 15 viral taxa from 11 viral families, encompassing many viruses previously reported as members of coral viromes. Despite the non-random selection of viral genetic material due to the nature of the sequencing data analyzed, our study provides a first insight into the viral community associated with Aiptasia. Similarities of the Aiptasia viral community with those of corals corroborate the application of Aiptasia as a model system to study coral holobionts. Further, the change in abundance of certain viral taxa across different symbiotic states suggests a role of viruses in the algal endosymbiosis, but the functional significance of this remains to be determined.

The ancestral dinoflagellate most likely established a peridinin-containing plastid, which have been inherited in the extant photosynthetic descendants. However, kareniacean dinoflagellates and Lepidodinium species were known to bear "non-canonical" plastids lacking peridinin, which were established through haptophyte and green algal endosymbioses, respectively. For plastid function and maintenance, the aforementioned dinoflagellates were known to use nucleus-encoded proteins vertically inherited from the ancestral dinoflagellates (vertically inherited- or VI-type), and those acquired from non-dinoflagellate organisms (including the endosymbiont). These observations indicated that the proteomes of the non-canonical plastids derived from a haptophyte and a green alga were modified by "exogenous" genes acquired from non-dinoflagellate organisms. However, there was no systematic evaluation addressing how "exogenous" genes reshaped individual metabolic pathways localized in a non-canonical plastid.

Seagrass beds provide a variety of ecosystem services, both within and outside the bounds of the habitat itself. Here we use environmental DNA (eDNA) amplicons to analyze a broad cross-section of taxa from ecological communities in and immediately surrounding eelgrass (Zostera marina). Sampling seawater along transects extending alongshore outward from eelgrass beds, we demonstrate that eDNA provides meter-scale resolution of communities in the field. We evaluate eDNA abundance indices for 13 major phylogenetic groups of marine and estuarine taxa along these transects, finding highly local changes linked with proximity to Z. marina for a diverse group of dinoflagellates, and for no other group of taxa. Eelgrass habitat is consistently associated with dramatic reductions in dinoflagellate abundance both within the contiguous beds and for at least 15 m outside, relative to nearby sites without eelgrass. These results are consistent with the hypothesis that eelgrass-associated communities have allelopathic effects on dinoflagellates, and that these effects can extend in a halo beyond the bounds of the contiguous beds. Because many dinoflagellates are capable of forming harmful algal blooms (HABs) toxic to humans and other animal species, the apparent salutary effect of eelgrass habitat on neighboring waters has important implications for public health as well as shellfish aquaculture and harvesting.

Algae with potential biotechnological applications in different industries are commonly isolated from the environment in order to obtain pure (axenic) stocks that can be safely stored for long periods of time. To obtain axenic cultures, antibiotics are frequently employed, and cryopreservation is applied to preserve standing stocks. However, many of these now standard methods were developed using strains derived from pristine to near-pristine environments and cold to temperate regions. The potential effect of the said methods on the life cycle and biochemical profile of algae isolates from hyper-eutrophic and constant high-temperature tropical regions is not well understood. These effects could potentially render them unsuitable for their intended biotechnological application. In this study, we conducted a genetic characterization (18S rRNA) and evaluated the effect of purification (the use of the antibiotic chloramphenicol, CAP) and cryopreservation (dimethyl sulfoxide; DMSO-sucrose mix and glycerol) on the growth rate and lipid content of three new tropical freshwater algal isolates: Chorella sp. M2, Chlorella sp. M6, and Scenedesmus sp. R3, obtained from the Ecuadorian coast. The genetic and morphological characterization revealed a clear discrimination between these strains. All strains cultured with CAP exhibited a lower growth rate. Subsequent to cryopreservation, Chorella sp. M2, Chlorella sp. M6, and Scenedesmus sp. R3 presented no significant difference in growth rate between the cryopreservants. Further, a significantly higher lipid content was observed in the biomass cryopreserved with glycerol in relation to the DMSO-sucrose, with Chorella sp. M2 and Chlorella sp. M6 having twice as much as they had in the first treatment. These results highlight the relevance of selecting an appropriate method for storage, as the materials used can affect the biological performance of different tropical species, although it is still to be determined if the effects observed in this study are long lasting in subsequent cultures of these algae.

Molecular stress responses associated with coral diseases represent an under-studied area of cnidarian transcriptome investigations. Caribbean Yellow Band Disease (CYBD) is considered a disease of Symbiodinium within the tissues of the coral host Orbicella faveolata. There is a paucity of diagnostic tools to assist in the early detection and characterization of coral diseases. The validity of a diagnostic test is determined by its ability to distinguish host organisms that have the disease from those that do not. The ability to detect and identify disease-affected tissue before visible signs of the disease are evident would then be a useful diagnostic tool for monitoring and managing disease outbreaks. Representational Difference Analysis (RDA) was utilized to isolate differentially expressed genes in O. faveolata exhibiting CYBD. Preliminary screening of RDA products identified a small number of genes of interest (GOI) which included an early growth response factor and ubiquitin ligase from the coral host as well as cytochrome oxidase from the algal symbiont. To further characterize the specificity of response, quantitative real-time PCR (qPCR) was utilized to compare the expression profiles of these GOIs within diseased tissues (visible lesions), tissues that precede visible lesions by 2-4 cm (transition area), and tissues from healthy-looking colonies with no signs of disease. Results show there are distinctive differences in the expression profiles of these three GOIs within each tissue examined. Collectively, this small suite of GOIs can provide a molecular "finger print" which is capable of differentiating between infected and uninfected colonies on reefs where CYBD is known to occur.