The discovery of microbial rhodopsins in marine proteobacteria changed the dogma that photosynthesis is the only pathway to use the solar energy for biological utilization in the marine environment. Although homologs of these rhodopsins have been identified in dinoflagellates, the diversity of the encoding genes and their physiological roles remain unexplored. As an initial step toward addressing the gap, we conducted high-throughput transcriptome sequencing on Oxyrrhis marina to retrieve rhodopsin transcripts, rapid amplification of cDNA ends to isolate full-length cDNAs of dominant representatives, and quantitative reverse-transcription PCR to investigate their expression under varying conditions. Our phylogenetic analyses showed that O. marina contained both the proton-pumping type (PR) and sensory type (SR) rhodopsins, and the transcriptome data showed that the PR type dominated over the SR type. We compared rhodopsin gene expression for cultures kept under light: dark cycle and continuous darkness in a time course of 24 days without feeding. Although both types of rhodopsin were expressed under the two conditions, the expression levels of PR were much higher than SR, consistent with the transcriptomic data. Furthermore, relative to cultures kept in the dark, rhodopsin expression levels and cell survival rate were both higher in cultures grown in the light. This is the first report of light-dependent promotion of starvation survival and concomitant promotion of PR expression in a eukaryote. While direct evidence needs to come from functional test on rhodopsins in vitro or gene knockout/knockdown experiments, our results suggest that the proton-pumping rhodopsin might be responsible for the light-enhanced survival of O. marina, as previously demonstrated in bacteria.

The genus Perkinsus occupies a precarious phylogenetic position. To gain a better understanding of the relationship between perkinsids, dinoflagellates and other alveolates, we analyzed the nuclear-encoded spliced-leader (SL) RNA and mitochondrial genes, intron prevalence, and multi-protein phylogenies. In contrast to the canonical 22-nt SL found in dinoflagellates (DinoSL), P. marinus has a shorter (21-nt) and a longer (22-nt) SL with slightly different sequences than DinoSL. The major SL RNA transcripts range in size between 80–83 nt in P. marinus, and ∼83 nt in P. chesapeaki, significantly larger than the typical ≤56-nt dinoflagellate SL RNA. In most of the phylogenetic trees based on 41 predicted protein sequences, P. marinus branched at the base of the dinoflagellate clade that included the ancient taxa Oxyrrhis and Amoebophrya, sister to the clade of apicomplexans, and in some cases clustered with apicomplexans as a sister to the dinoflagellate clade. Of 104 Perkinsus spp. genes examined 69.2% had introns, a higher intron prevalence than in dinoflagellates. Examination of Perkinsus spp. mitochondrial cytochrome B and cytochrome C oxidase subunit I genes and their cDNAs revealed no mRNA editing, but these transcripts can only be translated when frameshifts are introduced at every AGG and CCC codon as if AGGY codes for glycine and CCCCU for proline. These results, along with the presence of the numerous uncharacterized ‘marine alveolate group I' and Perkinsus-like lineages separating perkinsids from core dinoflagellates, expand support for the affiliation of the genus Perkinsus with an independent lineage (Perkinsozoa) positioned between the phyla of Apicomplexa and Dinoflagellata.

Coral reef-associated Symbiodiniaceae live in tropical and oligotrophic environments and are prone to heat and nutrient stress. How their metabolic pathways respond to pulses of warming and phosphorus (P) depletion is underexplored. Here, we conducted RNA-seq analysis to investigate transcriptomic responses to thermal stress, phosphate deprivation, and organic phosphorus (OP) replacement in Fugacium kawagutii. Using dual-algorithm (edgeR and NOIseq) to remedy the problem of no replicates, we conservatively found 357 differentially expressed genes (DEGs) under heat stress, potentially regulating cell wall modulation and the transport of iron, oxygen, and major nutrients. About 396 DEGs were detected under P deprivation and 671 under OP utilization, both mostly up-regulated and potentially involved in photosystem and defensome, despite different KEGG pathway enrichments. Additionally, we identified 221 genes that showed relatively stable expression levels across all conditions (likely core genes), mostly catalytic and binding proteins. This study reveals a wide range of, and in many cases previously unrecognized, molecular mechanisms in F. kawagutii to cope with heat stress and phosphorus-deficiency stress. Their quantitative expression dynamics, however, requires further verification with triplicated experiments, and the data reported here only provide clues for generating testable hypotheses about molecular mechanisms underpinning responses and adaptation in F. kawagutii to temperature and nutrient stresses.

Toward understanding the genetic diversity and distribution of copepod-associated symbiotic ciliates and the evolutionary relationships with their hosts in the marine environment, we developed a small subunit ribosomal RNA gene (18S rDNA)-based molecular method and investigated the genetic diversity and genotype distribution of the symbiotic ciliates on copepods. Of the 10 copepod species representing six families collected from six locations of Pacific and Atlantic Oceans, 9 were found to harbor ciliate symbionts. Phylogenetic analysis of the 391 ciliate 18S rDNA sequences obtained revealed seven groups (ribogroups), six (containing 99% of all the sequences) belonging to subclass Apostomatida, the other clustered with peritrich ciliate Vorticella gracilis. Among the Apostomatida groups, Group III were essentially identical to Vampyrophrya pelagica, and the other five groups represented the undocumented ciliates that were close to Vampyrophrya/Gymnodinioides/Hyalophysa. Group VI ciliates were found in all copepod species but one (Calanus sinicus), and were most abundant among all ciliate sequences obtained, indicating that they are the dominant symbiotic ciliates universally associated with copepods. In contrast, some ciliate sequences were found only in some of the copepods examined, suggesting the host selectivity and geographic differentiation of ciliates, which requires further verification by more extensive sampling. Our results reveal the wide occurrence and high genetic diversity of symbiotic ciliates on marine copepods and highlight the need to systematically investigate the host- and geography-based genetic differentiation and ecological roles of these ciliates globally.