A Potential Role for Epigenetic Processes in the Acclimation Response to Elevated pCO2 in the Model Diatom Phaeodactylum tricornutum

Transcriptional responses to salinity-induced changes in cell wall morphology of the euryhaline diatom Pleurosira laevis

Diatoms are unicellular algae with morphologically diverse silica cell walls, which are called frustules. The mechanism of frustule morphogenesis has attracted attention in biology and nanomaterials engineering. However, the genetic regulation of the morphology remains unclear. We therefore used transcriptome sequencing to search for genes involved in frustule morphology in the centric diatom Pleurosira laevis, which exhibits morphological plasticity between flat and domed valve faces in salinity 2 and 7, respectively. We observed differential expression of transposable elements (TEs) and transporters, likely due to osmotic response. Up-regulation of mechanosensitive ion channels and down-regulation of Ca2+-ATPases in cells with flat valves suggested that cytosolic Ca2+ levels were changed between the morphologies. Calcium signaling could be a mechanism for detecting osmotic pressure changes and triggering morphological shifts. We also observed an up-regulation of ARPC1 and annexin, involved in the regulation of actin filament dynamics known to affect frustule morphology, as well as the up-regulation of genes encoding frustule-related proteins such as BacSETs and frustulin. Taken together, we propose a model in which salinity-induced morphogenetic changes are driven by upstream responses, such as the regulation of cytosolic Ca2+ levels, and downstream responses, such as Ca2+-dependent regulation of actin dynamics and frustule-related proteins. This study highlights the sensitivity of euryhaline diatoms to environmental salinity and the role of active cellular processes in controlling gross valve morphology under different osmotic pressures.

Complex dynamics of coral gene expression responses to low pH across species

Coral capacity to tolerate low pH affects coral community composition and, ultimately, reef ecosystem function. Low pH submarine discharges ('Ojo'; Yucatán, México) represent a natural laboratory to study plasticity and acclimatization to low pH in relation to ocean acidification. A previous >2-year coral transplant experiment to ambient and low pH common garden sites revealed differential survivorship across species and sites, providing a framework to compare mechanistic responses to differential pH exposures. Here, we examined gene expression responses of transplants of three species of reef-building corals (Porites astreoides, Porites porites and Siderastrea siderea) and their algal endosymbiont communities (Symbiodiniaceae) originating from low pH (Ojo) and ambient pH native origins (Lagoon or Reef). Transplant pH environment had the greatest effect on gene expression of Porites astreoides hosts and symbionts and P. porites hosts. Host P. astreoides Ojo natives transplanted to ambient pH showed a similar gene expression profile to Lagoon natives remaining in ambient pH, providing evidence of plasticity in response to ambient pH conditions. Although origin had a larger effect on host S. siderea gene expression due to differences in symbiont genera within Reef and Lagoon/Ojo natives, subtle effects of low pH on all origins demonstrated acclimatization potential. All corals responded to low pH by differentially expressing genes related to pH regulation, ion transport, calcification, cell adhesion and stress/immune response. This study demonstrates that the magnitude of coral gene expression responses to pH varies considerably among populations, species and holobionts, which could differentially affect acclimatization to and impacts of ocean acidification.

Lessons from the deep: mechanisms behind diversification of eukaryotic protein complexes

Genetic variation is the major mechanism behind adaptation and evolutionary change. As most proteins operate through interactions with other proteins, changes in protein complex composition and subunit sequence provide potentially new functions. Comparative genomics can reveal expansions, losses and sequence divergence within protein-coding genes, but in silico analysis cannot detect subunit substitutions or replacements of entire protein complexes. Insights into these fundamental evolutionary processes require broad and extensive comparative analyses, from both in silico and experimental evidence. Here, we combine data from both approaches and consider the gamut of possible protein complex compositional changes that arise during evolution, citing examples of complete conservation to partial and total replacement by functional analogues. We focus in part on complexes in trypanosomes as they represent one of the better studied non-animal/non-fungal lineages, but extend insights across the eukaryotes by extensive comparative genomic analysis. We argue that gene loss plays an important role in diversification of protein complexes and hence enhancement of eukaryotic diversity.

Epigenetic Regulation in Response to CO2 Fluctuation in Marine Microalga Nannochloropsis oceanica

Microalgae often undergo different CO2 experiment in their habitat. To adapt to low CO2, carbon concentrating mechanism (CCM) could be launched in majority of microalgae and CCM are regulated at RNA level are well known. However, epigenetic modifications and their potential regulation of the transcription of masked genes at the genome level in response to CO2 fluctuation remain unclear. Here epigenetic regulation in response to CO2 fluctuation and epigenome-association with phenotypic plasticity of CCM are firstly uncovered in marine microalga Nannochloropsis oceanica IMET1. The result showed that lysine butyrylation (Kbu) and histone H3K9m2 modifications were present in N. oceanica IMET1. Moreover, Kbu modification positively regulated gene expression. In response to CO2 fluctuation, there were 5,438 and 1,106 genes regulated by Kbu and H3K9m2 in Nannochloropsis, respectively. Gained or lost histone methylations were closely associated with activating or repressing gene expressions. Differential modifications were mainly enriched in carbon fixation, photorespiration, photosynthesis, and lipid metabolism etc. Massive genome-wide epigenetic reprogramming was observed after N. oceanica cells shifted from high CO2 to low CO2. Particularly, we firstly noted that the transcription of the key low CO2 responsive carbonic anhydrase (CA5), a key component involved in CCM stress signaling, was potentially regulated by bivalent Kbu-H3K9m2 modifications in microalgae. This study provides novel insights into the relationship between gene transcription and epigenetic modification in Nannochloropsis, which will lay foundation on genetic improvement of CCM at epigenetic level.

Light-dependent signal transduction in the marine diatom Phaeodactylum tricornutum

Unlike most higher plants, unicellular algae can acclimate to changes in irradiance on time scales of hours to a few days. The process involves an enigmatic signaling pathway originating in the plastid that leads to coordinated changes in plastid and nuclear gene expression. To deepen our understanding of this process, we conducted functional studies to examine how the model diatom, Phaeodactylum tricornutum, acclimates to low light and sought to identify the molecules responsible for the phenomenon. We show that two transformants with altered expression of two putative signal transduction molecules, a light-specific soluble kinase and a plastid transmembrane protein, that appears to be regulated by a long noncoding natural antisense transcript, arising from the opposite strand, are physiologically incapable of photoacclimation. Based on these results, we propose a working model of the retrograde feedback in the signaling and regulation of photoacclimation in a marine diatom.

Differential expression patterns of long noncoding RNAs in a pleiomorphic diatom and relation to hyposalinity

Long non-coding (lnc)RNAs have been shown to have central roles in stress responses, cell identity and developmental processes in multicellular organisms as well as in unicellular fungi. Previous works have shown the occurrence of lncRNAs in diatoms, namely in Phaeodactylum tricornutum, many of which being expressed under specific stress conditions. Interestingly, P. tricornutum is the only known diatom that has a demonstrated morphological plasticity, occurring in three distinct morphotypes: fusiform, triradiate and oval. Although the morphotypes are interchangeable, the fusiform is the dominant one while both the triradiate and the oval forms are less common, the latter often being associated with stress conditions such as low salinity and solid culture media, amongst others. Nonetheless, the molecular basis underpinning morphotype identity in P. tricornutum remains elusive. Using twelve previously published transcriptomic datasets originating from the three morphotypes of P. tricornutum, we sought to investigate the expression patterns of lncRNAs (lincRNAs and NATs) in these distinct morphotypes, using pairwise comparisons, in order to explore the putative involvement of these noncoding molecules in morphotype identity. We found that differentially expressed lncRNAs cluster according to morphotype, indicating that lncRNAs are not randomly expressed, but rather seem to provide a specific (noncoding) transcriptomic signature of the morphotype. We also present evidence to suggest that the major differences in DE genes (both noncoding and coding) between the stress related oval morphotype and the most common fusiform morphotype could be due, to a large extent, to the hyposaline culture conditions rather than to the morphotype itself. However, several lncRNAs associated to each one of the three morphotypes were identified, which could have a potential role in morphotype (or cell) identity in P. tricornutum, similar to what has been found in both animals and plant development.

Environmental Adaptation of Genetically Uniform Organisms with the Help of Epigenetic Mechanisms-An Insightful Perspective on Ecoepigenetics

Organisms adapt to different environments by selection of the most suitable phenotypes from the standing genetic variation or by phenotypic plasticity, the ability of single genotypes to produce different phenotypes in different environments. Because of near genetic identity, asexually reproducing populations are particularly suitable for the investigation of the potential and molecular underpinning of the latter alternative in depth. Recent analyses on the whole-genome scale of differently adapted clonal animals and plants demonstrated that epigenetic mechanisms such as DNA methylation, histone modifications and non-coding RNAs are among the molecular pathways supporting phenotypic plasticity and that epigenetic variation is used to stably adapt to different environments. Case studies revealed habitat-specific epigenetic fingerprints that were maintained over subsequent years pointing at the existence of epigenetic ecotypes. Environmentally induced epimutations and corresponding gene expression changes provide an ideal means for fast and directional adaptation to changing or new conditions, because they can synchronously alter phenotypes in many population members. Because microorganisms inclusive of human pathogens also exploit epigenetically mediated phenotypic variation for environmental adaptation, this phenomenon is considered a universal biological principle. The production of different phenotypes from the same DNA sequence in response to environmental cues by epigenetic mechanisms also provides a mechanistic explanation for the "general-purpose genotype hypothesis" and the "genetic paradox of invasions".

Transcriptomic analysis reveals distinct mechanisms of adaptation of a polar picophytoplankter under ocean acidification conditions

Human emissions of carbon dioxide are causing irreversible changes in our oceans and impacting marine phytoplankton, including a group of small green algae known as picochlorophytes. Picochlorophytes grown in natural phytoplankton communities under future predicted levels of carbon dioxide have been demonstrated to thrive, along with redistribution of the cellular metabolome that enhances growth rate and photosynthesis. Here, using next-generation sequencing technology, we measured levels of transcripts in a picochlorophyte Chlorella, isolated from the sub-Antarctic and acclimated under high and current ambient CO2 levels, to better understand the molecular mechanisms involved with its ability to acclimate to elevated CO2. Compared to other phytoplankton taxa that induce broad transcriptomic responses involving multiple parts of their cellular metabolism, the changes observed in Chlorella focused on activating gene regulation involved in different sets of pathways such as light harvesting complex binding proteins, amino acid synthesis and RNA modification, while carbon metabolism was largely unaffected. Triggering a specific set of genes could be a unique strategy of small green phytoplankton under high CO2 in polar oceans.

Identification and expression analysis of genome-wide long noncoding RNA responsive CO2 fluctuated environment in marine microalga Nannochloropsis oceanica

Long non-coding RNAs (lncRNAs) have been demonstrated to participate in plant growth and development as well as response to different biotic and abiotic stresses. However, the knowledge of lncRNA was limited in microalgae. In this study, by RNA deep sequencing, 134 lncRNAs were identified in marine Nannochloropsis oceanica in response to carbon dioxide fluctuation. Among them, there were 51 lncRNAs displayed differentially expressed between low and high CO₂ treatments, including 33 upregulation and 18 downregulation lncRNAs. Cellulose metabolic process, glucan metabolic process, polysaccharide metabolic process, and transmembrane transporter activity were functionally enriched. Multiple potential target genes of lncRNA and lncRNA-mRNA co-located gene network were analyzed. Subsequent analysis had demonstrated that lncRNAs would participate in many biological molecular processes, including gene expression, transcriptional regulation, protein expression and epigenetic regulation. In addition, alternative splicing events were firstly analyzed in response to CO₂ fluctuation. There were 2051 alternative splicing (AS events) identified, which might be associated with lncRNA. These observations will provide a novel insight into lncRNA function in Nannochloropsis and provide a series of targets for lncRNA-based gene editing in future.

Diurnal transcript profiling of the diatom Seminavis robusta reveals adaptations to a benthic lifestyle

Coastal regions contribute an estimated 20% of annual gross primary production in the oceans, despite occupying only 0.03% of their surface area. Diatoms frequently dominate coastal sediments, where they experience large variations in light regime resulting from the interplay of diurnal and tidal cycles. Here, we report on an extensive diurnal transcript profiling experiment of the motile benthic diatom Seminavis robusta. Nearly 90% (23 328) of expressed protein-coding genes and 66.9% (1124) of expressed long intergenic non-coding RNAs showed significant expression oscillations and are predominantly phasing at night with a periodicity of 24 h. Phylostratigraphic analysis found that rhythmic genes are enriched in highly conserved genes, while diatom-specific genes are predominantly associated with midnight expression. Integration of genetic and physiological cell cycle markers with silica depletion data revealed potential new silica cell wall-associated gene families specific to diatoms. Additionally, we observed 1752 genes with a remarkable semidiurnal (12-h) periodicity, while the expansion of putative circadian transcription factors may reflect adaptations to cope with highly unpredictable external conditions. Taken together, our results provide new insights into the adaptations of diatoms to the benthic environment and serve as a valuable resource for the study of diurnal regulation in photosynthetic eukaryotes.

Ultraviolet Radiation Stimulates Activity of CO2 Concentrating Mechanisms in a Bloom-Forming Diatom Under Reduced CO2 Availability

The diatom Skeletonema costatum is cosmopolitan and forms algal blooms in coastal waters, being exposed to varying levels of solar UV radiation (UVR) and reduced levels of carbon dioxide (CO₂). While reduced CO₂ availability is known to enhance CO₂ concentrating mechanisms (CCMs) in this diatom and others, little is known on the effects of UV on microalgal CCMs, especially when CO₂ levels fluctuate in coastal waters. Here, we show that S. costatum upregulated its CCMs in response to UVR (295–395 nm), especially to UVA (320–395 nm) in the presence and absence of photosynthetically active radiation (PAR). The intensity rise of UVA and/or UVR alone resulted in an increase of the activity of extracellular carbonic anhydrase (CAe); and the addition of UVA enhanced the activity of CCMs-related CAe by 23–27% when PAR levels were low. Such UV-stimulated CCMs activity was only significant at the reduced CO₂ level (3.4 μmol L⁻¹). In addition, UVA alone drove active HCO₃⁻ uptake although it was not as obvious as CAe activity, another evidence for its role in enhancing CCMs activity. In parallel, the addition of UVA enhanced photosynthetic carbon fixation only at the lower CO₂ level compared to PAR alone. In the absence of PAR, carbon fixation increased linearly with increased intensities of UVA or UVR regardless of the CO₂ levels. These findings imply that during S. costatum blooming period when CO₂ and PAR availability becomes lower, solar UVR (mainly UVA) helps to upregulate its CCMs and thus carbon fixation, enabling its success of frequent algal blooms.

Long-Term m5C Methylome Dynamics Parallel Phenotypic Adaptation in the Cyanobacterium Trichodesmium

A major challenge in modern biology is understanding how the effects of short-term biological responses influence long-term evolutionary adaptation, defined as a genetically determined increase in fitness to novel environments. This is particularly important in globally important microbes experiencing rapid global change, due to their influence on food webs, biogeochemical cycles, and climate. Epigenetic modifications like methylation have been demonstrated to influence short-term plastic responses, which ultimately impact long-term adaptive responses to environmental change. However, there remains a paucity of empirical research examining long-term methylation dynamics during environmental adaptation in nonmodel, ecologically important microbes. Here, we show the first empirical evidence in a marine prokaryote for long-term m5C methylome modifications correlated with phenotypic adaptation to CO2, using a 7-year evolution experiment (1,000+ generations) with the biogeochemically important marine cyanobacterium Trichodesmium. We identify m5C methylated sites that rapidly changed in response to high (750 µatm) CO2 exposure and were maintained for at least 4.5 years of CO2 selection. After 7 years of CO2 selection, however, m5C methylation levels that initially responded to high-CO2 returned to ancestral, ambient CO2 levels. Concurrently, high-CO2 adapted growth and N2 fixation rates remained significantly higher than those of ambient CO2 adapted cell lines irrespective of CO2 concentration, a trend consistent with genetic assimilation theory. These data demonstrate the maintenance of CO2-responsive m5C methylation for 4.5 years alongside phenotypic adaptation before returning to ancestral methylation levels. These observations in a globally distributed marine prokaryote provide critical evolutionary insights into biogeochemically important traits under global change.

The Aureochrome Photoreceptor PtAUREO1a Is a Highly Effective Blue Light Switch in Diatoms

Aureochromes represent a unique type of blue light photoreceptors that possess a blue light sensing flavin-binding LOV-domain and a DNA-binding bZIP domain, thus being light-driven transcription factors. The diatom Phaeodactylum tricornutum, a member of the essential marine primary producers, possesses four aureochromes (PtAUREO1a, 1b, 1c, 2). Here we show a dramatic change in the global gene expression pattern of P. tricornutum wild-type cells after a shift from red to blue light. About 75% of the genes show significantly changed transcript levels already after 10 and 60 min of blue light exposure, which includes genes of major transcription factors as well as other photoreceptors. Very surprisingly, this light-induced regulation of gene expression is almost completely inhibited in independent PtAureo1a knockout lines. Such a massive and fast transcriptional change depending on one single photoreceptor is so far unprecedented. We conclude that PtAUREO1a plays a key role in diatoms upon blue light exposure.

•

Blue light induces a very fast transcriptional response in the diatom P. tricornutum

•

This strong response is almost completely inhibited when Aureochrome 1a is absent

•

The results imply a key role of PtAureo1a in blue light-induced responses in diatoms

Genome-wide identification and transcriptional modulation of histone variants and modification related genes in the low pH-exposed marine rotifer Brachionus koreanus

Histone modification is considered to be a major epigenetic control mechanism. These modifications (e.g. acetylation, phosphorylation, and methylation) may affect the interaction of histones with DNA and/or regulate DNA-based processes (e.g., recombination, repair, replication, and transcription) and chromatin remodeling complexes. Despite their significance in metazoan life and evolution, few studies have been conducted to identify genes undergoing epigenetic control modification in aquatic invertebrates. In this study, we identified whole core histones (70 total genes) and post-translational modification (PTM) histone genes (63 total genes) in the marine rotifer Brachionus koreanus through whole-genome analysis, and annotated them according to the human nomenclature. Notably, upon comparative analysis of cis-regulatory motif sequences, we found that B. koreanus core histone protein structures were similar to those of mammals. Furthermore, to examine the effect of parental low pH stress on the offspring's epigenetic regulation, we investigated the expression of PTM genes in two generations of B. koreanus exposed to low pH conditions. Given that the B. koreanus genome does not possess DNA methyltransferase 1 and 3 genes, we concluded that histone genes could be involved as an important epigenetic mechanism in B. koreanus. Therefore, the histone-associated genes identified in this study could be useful for ecotoxicological studies and facilitate the application of chromatin immunoprecipitation sequencing using high-throughput DNA sequencing based on the genome-wide identification of transcription factor binding sites in rotifers.

Discovering the indigenous microbial communities associated with the natural fermentation of sap from the cider gum Eucalyptus gunnii

Over the course of human history and in most societies, fermented beverages have had a unique economic and cultural importance. Before the arrival of the first Europeans in Australia, Aboriginal people reportedly produced several fermented drinks including mangaitch from flowering cones of Banksia and way-a-linah from Eucalyptus tree sap. In the case of more familiar fermented beverages, numerous microorganisms, including fungi, yeast and bacteria, present on the surface of fruits and grains are responsible for the conversion of the sugars in these materials into ethanol. Here we describe native microbial communities associated with the spontaneous fermentation of sap from the cider gum Eucalyptus gunnii, a Eucalyptus tree native to the remote Central Plateau of Tasmania. Amplicon-based phylotyping showed numerous microbial species in cider gum samples, with fungal species differing greatly to those associated with winemaking. Phylotyping also revealed several fungal sequences which do not match known fungal genomes suggesting novel yeast species. These findings highlight the vast microbial diversity associated with the Australian Eucalyptus gunnii and the native alcoholic beverage way-a-linah.

Global identification of a marine diatom long noncoding natural antisense transcripts (NATs) and their response to phosphate fluctuations

Often ignored and regarded as mere transcriptional noise, long noncoding RNAs (lncRNAs) are starting to be considered key regulators of gene expression across the Eukarya domain of life. In the model diatom Phaeodactylum tricornutum, we have previously reported the occurrence of 1,510 intergenic lncRNAs (lincRNAs), many of which displaying specific patterns of expression under phosphate fluctuation (Pi). Using strand-specific RNA-sequencing data we now expand the repertoire of P. tricornutum lncRNAs by identifying 2,628 novel natural antisense transcripts (NATs) that cover 21.5% of the annotated genomic loci. We found that NAT expression is tightly regulated by phosphate depletion and other naturally occurring environmental stresses. Furthermore, we identified 121 phosphate stress responsive NAT-mRNA pairs, the great majority of which showing a positive correlation (concordant pairs) and a small fraction with negative correlation (discordant pairs). Taken together our results show that NATs are highly abundant transcripts in P. tricornutum and that their expression is under tight regulation by nutrient and environmental stresses. Furthermore, our results suggest that in P. tricornutum Pi stress response NAT pairs predominantly regulate positively the expression of their cognate sense genes, the latter being involved in several biological processes underlying the control of cellular homeostasis under stress.

Toxic effect and mechanism of tris (1,3-dichloro-2-propyl)phosphate (TDCPP) on the marine alga Phaeodactylum tricornutum

Tris (1,3-dichloro-2-propyl)phosphate (TDCPP) is an organophosphate-based plasticizer and flame retardant with a high production volume. The ubiquitous distribution and persistence of TDCPP in aquatic environment have led to concerns over its possible toxic effects on aquatic organism. However, data regarding the toxicity of TDCPP on algae are limited, and the molecular mechanism remains largely unknown. Therefore, we determined the growth characteristics, physiological changes and transcriptome profiles of Phaeodactylum tricornutum in response to 4 mg L^-1 TDCPP for 24 h. TDCPP caused morphological damage and growth inhibition with an EC₅₀ value of 3.71 mg L^-1 at 96 h. A decline in pigments and photosynthetic activity was observed, indicating the occurrence of photosynthesis inhibition. Although the activities of both glutathione peroxidase and glutathione reductase were stimulated, oxidative stress was not relieved in the algal cells, as evidenced by the elevated levels of reactive oxygen species and lipid peroxidation. Transcriptomic analyses revealed 3312 differentially expressed genes (DEGs), and photosynthesis was a key target, as genes related to this process were greatly altered under TDCPP stress. Moreover, some DEGs were also enriched in amino acid metabolism, nitrogen metabolism, nucleotide metabolism and lipid metabolism, implying that TDCPP-induced damage towards algae by various pathways. Additionally, several TFs related to stress signaling were differentially expressed, suggesting roles in the TDCPP stress response. The results will provide critical data to understand the ecological risks and toxic mechanism of OPFRs entering into marine habitat.

Transgenerational acclimation to changes in ocean acidification in marine invertebrates

The rapid pace of increasing oceanic acidity poses a major threat to the fitness of the marine ecosystem, as well as the buffering capacity of the oceans. Disruption in chemical equilibrium in the ocean leads to decreased carbonate ion precipitation, resulting in calcium carbonate saturation. If these trends continue, calcifying invertebrates will experience difficultly maintaining their calcium carbonate exoskeleton and shells. Because malfunction of exoskeleton formation by calcifiers in response to ocean acidification (OA) will have non-canonical biological cascading results in the marine ecosystem, many studies have investigated the direct and indirect consequences of OA on ecosystem- and physiology-related traits of marine invertebrates. Considering that evolutionary adaptation to OA depends on the duration of OA effects, long-term exposure to OA stress over multi-generations may result in adaptive mechanisms that increase the potential fitness of marine invertebrates in response to OA. Transgenerational studies have the potential to elucidate the roles of acclimation, carryover effects, and evolutionary adaptation within and over generations in response to OA. In particular, understanding mechanisms of transgenerational responses (e.g., antioxidant responses, metabolic changes, epigenetic reprogramming) to changes in OA will enhance our understanding of marine invertebrate in response to rapid climate change.

Non-Coding RNAs Participate in the Regulation of CRY-DASH in the Growth and Early Development of Saccharina japonica (Laminariales, Phaeophyceae)

CRY-DASH, a new cryptochrome blue light receptor, can repair damaged DNA and regulate secondary metabolism and development of fungus. However, its role in regulation during the growth of Saccharina japonica is still unclear. After cloning the full-length of CRY-DASH from S. japonica (sjCRY-DASH), we deduced that its open reading frame was 1779 bp long and encoded 592 amino acids. sjCRY-DASH transcription was rapidly upregulated within 30 min in response to blue light and exhibited 24 h periodicity with different photoperiods. Moreover, sjCRY-DASH maintained the same periodicity in suitable growth temperature, suggesting a close relationship between this periodicity and circadian rhythm regulation. Novel-m3234-5p, which was targeted to sjCRY-DASH, decreased with increasing sjCRY-DASH transcription, acting as a negative modulator of sjCRY-DASH. Six long non-coding RNAs classified as long intergenic non-coding RNAs (lincRNAs) exhibited co-expression with sjCRY-DASH. A miRNA sjCRY DASH lincRNA network was consequently identified. By predicting the endogenous competing mRNAs of novel-m3234-5p, we found that sjCRY-DASH indirectly participated in the regulation of DNA damage repair, protein synthesis and processing, and actin transport. In conclusion, our results revealed that non-coding RNAs participate in the regulation of sjCRY-DASH, which played vital roles in the growth and early development of S. japonica.

Evolution of Codon Usage Bias in Diatoms

Codon usage bias (CUB)-preferential use of one of the synonymous codons, has been described in a wide range of organisms from bacteria to mammals, but it has not yet been studied in marine phytoplankton. CUB is thought to be caused by weak selection for translational accuracy and efficiency. Weak selection can overpower genetic drift only in species with large effective population sizes, such as Drosophila that has relatively strong CUB, while organisms with smaller population sizes (e.g., mammals) have weak CUB. Marine plankton species tend to have extremely large populations, suggesting that CUB should be very strong. Here we test this prediction and describe the patterns of codon usage in a wide range of diatom species belonging to 35 genera from 4 classes. We report that most of the diatom species studied have surprisingly modest CUB (mean Effective Number of Codons, ENC = 56), with some exceptions showing stronger codon bias (ENC = 44). Modest codon bias in most studied diatom species may reflect extreme disparity between astronomically large census and modest effective population size (N_e), with fluctuations in population size and linked selection limiting long-term N_e and rendering selection for optimal codons less efficient. For example, genetic diversity (pi ~0.02 at silent sites) in Skeletonema marinoi corresponds to N_e of about 10 million individuals, which is likely many orders of magnitude lower than its census size. Still, N_e ~10⁷ should be large enough to make selection for optimal codons efficient. Thus, we propose that an alternative process-frequent changes of preferred codons, may be a more plausible reason for low CUB despite highly efficient selection for preferred codons in diatom populations. The shifts in the set of optimal codons should result in the changes of the direction of selection for codon usage, so the actual codon usage never catches up with the moving target of the optimal set of codons and the species never develop strong CUB. Indeed, we detected strong shifts in preferential codon usage within some diatom genera, with switches between preferentially GC-rich and AT-rich 3^rd codon positions (GC3). For example, GC3 ranges from 0.6 to 1 in most Chaetoceros species, while for Chaetoceros dichaeta GC3 = 0.1. Both variation in selection intensity and mutation spectrum may drive such shifts in codon usage and limit the observed CUB. Our study represents the first genome-wide analysis of CUB in diatoms and the first such analysis for a major phytoplankton group.

A Potential Role for Epigenetic Processes in the Acclimation Response to Elevated pCO2 in the Model Diatom Phaeodactylum tricornutum

Understanding of the molecular responses underpinning diatom responses to ocean acidification is fundamental for predicting how important primary producers will be shaped by the continuous rise in atmospheric CO2. In this study, we have analyzed global transcriptomic changes of the model diatom Phaeodactylum tricornutum following growth for 15 generations in elevated pCO2 by strand-specific RNA sequencing (ssRNA-seq). Our results indicate that no significant effects of elevated pCO2 and associated carbonate chemistry changes on the physiological performance of the cells were observed after 15 generations whereas the expression of genes encoding histones and other genes involved in chromatin structure were significantly down-regulated, while the expression of transposable elements (TEs) and genes encoding histone acetylation enzymes were significantly up-regulated. Furthermore, we identified a series of long non-protein coding RNAs (lncRNAs) specifically responsive to elevated pCO2, suggesting putative regulatory roles for these largely uncharacterized genome components. Taken together, our integrative analyses reveal that epigenetic elements such as TEs, histone modifications and lncRNAs may have important roles in the acclimation of diatoms to elevated pCO2 over short time scales and thus may influence longer term adaptive processes in response to progressive ocean acidification.