Next Article in Journal
Improving the Traits of Perilla frutescens (L.) Britt Using Gene Editing Technology
Next Article in Special Issue
Photosynthetic Electron Flows and Networks of Metabolite Trafficking to Sustain Metabolism in Photosynthetic Systems
Previous Article in Journal
Seed-Coat Pigmentation Plays a Crucial Role in Partner Selection and N2 Fixation in Legume-Root–Microbe Associations in African Soils
Previous Article in Special Issue
Channels of Evolution: Unveiling Evolutionary Patterns in Diatom Ca2+ Signalling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Knockout of the Photoreceptor PtAureo1a Results in Altered Diel Expression of Diatom Clock Components

by
Shvaita Madhuri
,
Bernard Lepetit
,
Alexander Helmut Fürst
and
Peter G. Kroth
*
Fachbereich Biologie, Universität Konstanz, 78457 Konstanz, Germany
*
Author to whom correspondence should be addressed.
Current address: Institut für Biowissenschaften, Mathematisch-Naturwissenschaftliche Fakultät, Universität Rostock, Albert-Einstein-Straße 3, 18059 Rostock, Germany.
Plants 2024, 13(11), 1465; https://doi.org/10.3390/plants13111465
Submission received: 10 January 2024 / Revised: 10 May 2024 / Accepted: 17 May 2024 / Published: 25 May 2024

Abstract

:
Plants and algae use light not only for driving photosynthesis but also to sense environmental cues and to adjust their circadian clocks via photoreceptors. Aureochromes are blue-light-dependent photoreceptors that also function as transcription factors, possessing both a LOV and a bZIP domain. Aureochromes so far have only been detected in Stramenopile algae, which include the diatoms. Four paralogues of aureochromes have been identified in the pennate model diatom Phaeodactylum tricornutum: PtAureo1a, 1b, 1c, and 2. While it was shown recently that diatoms have a diel rhythm, the molecular mechanisms and components regulating it are still largely unknown. Diel gene expression analyses of wild-type P. tricornutum, a PtAureo1a knockout strain, and the respective PtAureo1 complemented line revealed that all four aureochromes have a different diel regulation and that PtAureo1a has a strong co-regulatory influence on its own transcription, as well as on that of other genes encoding different blue-light photoreceptors (CPF1, 2 and 4), proteins involved in photoprotection (Lhcx1), and specific bHLH transcription factors (RITMO1). Some of these genes completely lost their circadian expression in the PtAureo1a KO mutant. Our results suggest a major involvement of aureochromes in the molecular clock of diatoms.

1. Introduction

Photoautotrophic algae are of major ecological relevance [1,2], representing the base of the aquatic food web. Diatoms are unicellular algae belonging to the Stramenopiles, which are found ubiquitously in both fresh water and in marine systems. These algae represent genetically chimeric eukaryotic organisms. While red algae, green algae, and land plants evolved by uptake of a cyanobacterium that was taken up and converted to an organelle, diatoms evolved by secondary endosymbiosis, which is the uptake of a eukaryotic red alga into a eukaryotic host cell [3,4]. They also display a remarkable capacity to acclimate to changes in environmental conditions [5,6,7,8,9,10].
Light is an essential environmental factor for plants and algae, either as an energy source for photosynthesis or by providing the organisms with information about the surrounding environment [11,12,13]. The effects of fluctuating or oscillating light on algal photophysiology [14,15,16,17] and on competition dynamics are well studied [18,19,20]. Cellular motility can also be affected by specific light conditions, like, for instance, migration of diatoms in sediment mats [21,22] or intracellular plastid migration [23]. Light quality and intensity in these organisms can be sensed by photoreceptors, which usually are pigment-binding proteins that absorb specific wavelengths and which may prompt several cellular responses via signalling cascades or direct regulatory effects. Photoreceptors also have been shown to entrain the circadian clock during diurnal cycles in plants and algae [24,25,26,27].
Blue light is reaching deeper within the water column than other wavelengths of the visible spectrum because of an increased absorption of red light by water. While the ratio of red and far red light can have a strong impact on plant development in land plants [28], the contribution and regulation of red and blue light in diatoms is poorly understood [29]. A peculiar type of algal blue-light receptor was identified in 2007 in the xanthophyte alga Vaucheria frigida; it was named Aureochrome [30]. This protein, which also is found in diatoms [31], contains a blue-light responsive flavin-binding LOV domain (Light, Oxygen, Voltage) together with a bZIP (basic leucin zipper) domain, which is a typical domain of bZIP transcription factors [31,32]. Aureochromes have subsequently been detected in other Stramenopiles, but not in other organisms. Interestingly, they are missing in the oomycetes [31,33,34], a non-photosynthetic subgroup of the Stramenopiles [35]. Four paralogues of aureochromes (PtAureo1a, b, c, and 2) were identified in the model diatom P. tricornutum [34], with PtAureo1a being the currently best-studied aureochrome both in vivo [34,36,37,38] and in vitro [39,40,41,42]. A transcriptomic study on PtAureo1a knockout mutants revealed that this photoreceptor plays an essential role for a cellular reprogramming after a shift from red light to blue light [43]. The diel expression pattern of individual aureochromes apparently is differently regulated [39]. Therefore, aureochromes are highly interesting candidates for investigating a possible involvement in the regulation of diel biological rhythms of P. tricornutum. It was shown that PtAureo1a affects the expression of diatom-specific cyclin (dsCYC2), which is implicated to be important for the G1 to S phase cell cycle transition after dark arrest, and that PtAureo1a is involved in photoacclimation [34,37,44,45].
The diel light−dark cycle is a significant environmental parameter for most organisms. In phytoplankton, it can affect or regulate metabolic processes via the photoperiod, maximum irradiance, spectral composition, or oscillation patterns [17,46]. Microalgae can also show robust diel rhythms in regard to various processes such as growth, gene expression, pigment synthesis, etc., which has been shown in a variety of phytoplankton organisms [27,47,48,49,50]. In the centric diatom Thalassiosira pseudonana, about 40% of the transcriptome may change during day and night transitions [51]. In the pennate diatom Phaeodactylum tricornutum, it was shown that the cell cycle is following an endogenous circadian pattern, but its synchronization depends on the light−dark cycle [44,52]; however, the molecular mechanism of endogenous circadian clocks in the algae is still unknown. Recently, the RITMO1 protein (PtbHLH1a), a member of the bHLH-PAS transcription factor family, was found to be involved in the regulation of the circadian rhythms in P. tricornutum [53]. Overexpression of this protein in this diatom results in a shift of circadian fluorescence patterns. In the well-studied land plant Arabidopsis thaliana, various blue-light photoreceptors, such as ZTL (ZEITLUPE) and FKF1 (FLAVIN BINDING, KELCH REPEAT, F-BOX 1), interact with GIGANTEA (GI), a clock-associated protein [54,55,56], controlling processes like photoperiodism. In diatoms, the blue-light receptor PtCPF1 is the only photoreceptor that so far has been shown to be involved in circadian rhythms. It can repress CLOCK:BMAL1-regulated expression of a reporter gene in a heterologous system [57].
Here, we compared expression patterns of a P. tricornutum WT strain with the respective PtAureo1a knockout mutant KO8 [37] and the PtAureo1a complemented mutant Co.48 [38], which had been cultivated in photobioreactors under turbidostat conditions. We demonstrated that a circadian expression of specific genes is directly dependent on the presence of PtAureo1a, indicating that aureochrome proteins may play an important role as input signal for the diatom clock.

2. Results

2.1. Growth Characteristics of WT, KO8 and Co.48

In this study, we compare the diel expression of specific genes of P. tricornutum wild-type strain Pt4, the Aureochrome 1a knockout strain PtAureo1aKO8 [37], and a PtAureo1a complementation mutant, PtAureo1aCo.48. PtAureo1aCo.48 is one of a number of mutant lines which have previously been created by complementing a PtAureo1aKO8 strain with the PtAureo1a gene driven by its endogenous promoter [38]. PtAureo1aCo.48 expresses the PtAureo1a protein to a very similar extent as the wild-type strain [38] and shows the same phenotype.
As a prerequisite for comparing diel and potential circadian gene expression in the selected strains, we studied their growth characteristics. The growth rates during the exponential phase were similar between all three strains (one-way ANOVA, p = 0.156, Figure 1A). We furthermore determined the maximum photosynthetic efficiency of photosystem II (PSII) by measuring Fv/Fm of the three diatom strains in order to detect any differences in photosynthetic and metabolic activities. Figure 1B demonstrates that there was no significant difference in the maximum quantum yield of PSII between all three strains (one-way ANOVA, p = 0.083).

2.2. Diel Regulation of Diatom Genes by PtAureo1a

Two major factors may limit data quality when studying rhythmic or periodic gene expression patterns in cell cultures: (i) a variable cell density in the samples, and (ii) the degree of culture synchrony that influences the congruity of gene expression patterns related to specific processes, such as cell division. In order to eliminate these factors as much as possible, we have used a turbidostat system for continuous cultivation of the diatoms. Based on optical density measurements, cell densities corresponding to cultures during exponential growth were stabilized via controlled dilution of the culture. Such a dilution with fresh media prevents changes in nutrient concentrations, which is known to affect general gene expression [58]. After 5 days of acclimation to standard conditions (16:8 h light/dark, L/D), we took samples at specific time points for a period of 52 h, each under L/D or continuous darkness (D/D) conditions, and monitored the expression of ten different genes encoding photoreceptors and regulators (PtAureo1a, 1b, 1c, 2, PtCPF1, 2, 4, PtbHLH1a, 1b, and PtLhcx1, listed with accession numbers in Supplementary Table S2) (Figure 2).
PtAureo1a transcript abundance (here only the WT and complemented line are studied) peaks 4 h after light onset at Day1 T12/16, and this expression pattern is recurring the next day both in WT and Co.48 (Day2 T12/16), albeit with a higher transcript level in the complemented line. Interestingly, under D/D conditions, the expression of PtAureo1a peaks at Day4 T12/16, which is equivalent to Day1 T12/16 in L/D conditions, and this is recurring the next day (Day5 T12/16), implying that the oscillation of PtAureo1a expression continues in the dark with no considerable phase shift. Western blot analyses show that the amount of PtAureo1a protein peaks at time points 6 and 24 (Figure S1). In the dark, this fits to a delayed protein synthesis, as there is an increase in transcripts detectable in the middle-of-the-night phase, while the peak of transcription during the day does not result in a direct protein increase, which may be due to light-dependent proteasomal degradation of PtAureo1a [59].
For the other aureochrome genes, the diel mRNA oscillations are weaker but still detectable. The transcript profile of PtAureo1b shows a diel regulation during the day, while, in the dark, nearly no induction of transcription is observable, indicating that light is required for induction of transcription. In the PtAureo1a knockout strain, the PtAureo1b transcripts show weaker changes, which is indicating that PtAureo1a might be required for the diel expression of PtAureo1b. PtAureo1c mRNA oscillation in WT shows a similar pattern as PtAureo1a, peaking at Day1 T12, four hours after light onset. Expression of PtAureo1c in the PtAureoKO8 mutant line shows a weaker oscillation, which can partially be rescued in PtAureo1aCo48. Transcripts of both PtAureo1b and 1c show a weaker oscillation in the dark phase. Also, PtAureo2 mRNA expression in WT and PtAureo8Co.48 shows a changed pattern in the knockout mutant, which is in contrast to previous reports showing mostly constant diel expression of PtAureo2 [39].
We also studied the expression profile of further blue-light photoreceptors belonging to the cryptochrome/photolyase family which supposedly are involved in regulating the clock in diatoms [57]. PtCPF1 transcription is induced at the beginning of the day (Day1 T8 and Day2 T8), with an mRNA peak one hour after dark to light transition in WT, which is consistent with results reported by [57,60]. In D/D, the transcript expression increased at the same time points as at the beginning of the light phase in L/D, which is abolished in KO8 in the second D/D cycle. The swift increase in expression at the end of the dark period and the endurance of the expression level without light input suggests the ability to anticipate light onset, which is typical for circadian-regulated genes and circadian clock components [57]. A more intriguing observation was a disoriented diel expression of PtCPF1 in KO8, which was recovered in Co.48 under L/D, indicating a prominent direct or indirect involvement of PtAureo1a in its regulation. PtCPF2 mRNA is generally only lowly expressed and shows no major differential expression pattern in all strains under L/D. Under constant darkness, no transcripts could be detected, probably due to mRNA abundance below the detection limit. PtCPF4 mRNA expression shows stable expression throughout the day, with a peak at the end of the light period (Day1 T20 and Day2 T20) in WT and Co.48, while it is feeble in KO8. In constant darkness, there are modest mRNA oscillations in WT which are absent in KO8.
The PtLhcx1 protein provides energy dependent quenching capacity in cells grown under low light conditions [61]. The mRNA oscillation of PtLhcx1 shows a strong light-dependent diurnal regulation, with an expression peak at Day1 T8 and Day2 T8 both in WT and Co.48. In the KO8 line though, this strong fluctuation is completely abolished, suggesting a regulatory involvement of PtAureo1a in PtLhcx1 transcription. Under D/D conditions, the oscillation of PtLhcx1 gene is strongly dampened, in line with previous results [62], implying the importance of light entrainment for regulation of this gene.
The PtbHLH1a and 1b genes show a pronounced and regular diurnal fluctuation of transcript levels in WT. The rhythm is obscured in KO8 but rescued in Co.48. On the first day of constant darkness (Day 4), we observed an increase in PtbHLH1a and 1b mRNA abundance without light induction. This effect is lost the recurring day, while the level of expression remains constantly high in WT. In KO8, this pattern completely disappeared, advocating a direct/indirect involvement of PtAureo1a in regulation of the diel biological rhythm even under dark conditions.
Although, for most genes described here, the wild-type phenotype can mostly be rescued in the complemented mutant, the transcript abundances are not completely identical, which is probably due to a very sensitive response to already minor differences of the Aureochrome 1a protein. Additionally, a weaker rhythm or a reduced oscillation amplitude under continuous darkness could result from the lack of photosynthesis under constant darkness in P. tricornutum [53]. To avoid misleading conclusions from the diel expression of the complemented mutant under constant darkness, we did not interpolate the complemented (Co.48) mutant in D/D rhythm analysis.
In order to reveal statistically significant patterns, we analysed the transcript abundances using the BioDare2 software package (https://biodare2.ed.ac.uk/) while calculating Kendal rank correlation coefficients (τ), which describes the ordinal association between two measured quantities. Under L/D conditions, transcript abundances of PtAureo1a, 1c, PtCPF1 and 4, PtbHLH1a and 1b, and PtLhcx1 genes were found to be significantly rhythmic throughout the cycle (shown in Table 1). Here, the Kendall rank coefficients for PtCPF1, PtLhcx1 and both PtbHLH1 genes exhibited a very high coefficient of >0.8. For D/D, we did not calculate the coefficient, because the chosen sampling intervals were not regular enough.
A lack of PtAureo1a protein in the knockout mutant obviously alters the mRNA oscillation of PtCPF1, PtLhcx1, and PtbHLH1a and 1b significantly (p < 0.01), while rhythmicity is again restored in the complemented mutant (p < 0.05).
To further investigate the impact of PtAureo1a only on the significantly diurnally regulated genes with very high Kendall rank coefficients >0.8 (PtCPF1, PtLhcx1, PtbHLH1a and PtbHLH1b), we used the arithmetic mean of gene expression to determine the rhythmic parameter period, phase and amplitude by two different methods. The Fast Fourier Transform Non-Linear Least Squares algorithm (FFT-NLLS) is a classical algorithm that tries to fit the sum of cosine functions on a potential regulatory pattern [63]. It has the disadvantage of not recognizing a pattern when the rhythm does not follow a sinusoidal waveform [64]. For the latter, the Maximum Entropy Spectral Analysis (MESA) is better suited, as it does not assume any underlying waveform of the rhythm but is based on stochastic spectrum modelling [64]. Both modelling approaches are suitable also for a relatively short and noisy data series, such as in our case [64].
FFT-NLLS analysis allowed us to fit the data for PtCPF1, PtbHLH1a and PtbHLH1b in WT and Co.48 with a period roughly around 24 h (Table 2). In KO8, no fitting was possible for PtbHLH1a. This indicated that the sine pattern of PtbHLH1a transcript abundance observed and fitted in WT and Co.49 was lost in the KO8 mutant. Also, in this line the period for PtCPF1 and PtbHLH1b rhythmicity was longer, with >30 h clearly deviating from a circadian period (=24 h). Interestingly, the completely different MESA approach provided similar period estimates as FFT-NLLS analyses for all three genes in all three lines, also in terms of phase (i.e., time of the peak) and amplitude. The latter was strongly dampened in KO8 compared to both WT and Co.48, indicating that PtAureo1a does not only influence the timing of diurnal circulation but also its amplitude.
FFT-NLLS did not fit the rhythm of PtLhcx1 gene regulation to the sum of cosine functions. Here, we could only rely on MESA analysis. The period estimates were quite distant from 24 h (WT 29.2 h, KO8 27.7 h, Co.48 29.5 h). However, the amplitude was again largely dampened in KO8 compared to WT and the complemented mutant, and a clear phase shift was observed (4.9 h and 6.4 h in WT and Co.49, respectively, versus 10.8 h in KO8). At present, these data indicated an involvement of PtAureo1a in the oscillating daily PtLhcx1 mRNA pattern, but further and especially longer time series will be needed to provide more robust statistical support. Overall, the FFT_NLLS and MESA analyses supported the results of the Kendal rank analysis: PtAureo1a had a direct or indirect effect on the diurnal mRNA pattern of PtCPF1, PtbHLH1a, PtbHLH1b and PtLhcx1.

2.3. Correlation Analysis of Diel Expression in P. tricornutum

In addition to the rhythmicity analysis, we performed a correlation analysis with all studied genes in WT cells to elaborate on potential network connections under light/dark and in constant darkness. Figure 3 demonstrates that, under L/D conditions, the positive correlation between diel expression of aureochrome isoforms is significant, while the expression of the PtCPF1 gene correlates only to the PtLhcx1 gene. The expression of PtCPF2 is unrelated to any other tested genes, while PtCPF4 transcripts significantly correlate to those of PtbHLH1b, which again is strongly correlated to PtAureo1a and 1c. PtbHLH1a expression is positively correlated with PtAureo1b, 1c, PtbHLH1b, and PtCPF4. All the described correlations are not observed in the PtAureo1a knock out mutant (KO8) but are restored in the PtAureo1a complemented mutant (Co.48), indicating that PtAureo1a is an important regulatory factor. Furthermore, the transcripts of the PtAureo2 protein correlate either positively or negatively with PtbHLH1a and 1b expression, respectively, in the KO8 mutant, which is not seen in WT cells, indicating that this regulation is suppressed in the WT cells.
In D/D conditions, PtAureo1b and 1c expression is positively correlated in WT and KO8, independent of presence/absence of PtAureo1a (Figure 4). PtAureo1c transcript abundance is strongly correlated with that of PtLhcx1 in constant darkness in WT, but not in KO8. PtbHLH1a and 1b are found to be not correlated with each other in WT but positively correlated in KO8 under constant darkness.

3. Discussion

Biological rhythms, including diurnal regulation of gene expression, have been described for different algal taxa and for different environments [50,65], while the mechanisms of algal clocks are largely unknown. Auto-regulatory transcription−translation regulatory feedback loops (TTFLs) that result in an oscillation of mRNA and proteins are key components of biological clocks [66]. One potential regulator of TTFLS can be light, which is perceived by photoreceptors, that may regulate or adjust the biological clock in plants [67] and green algae [68]. Even non-photosynthetic organisms like fungi may have such TTFLs, involving the White Collar Complex (WCC) photoreceptors of the fungus Neurospora crassa [69]. However, the genes encoding TTFL components lack a genetic conservation across the kingdoms [46], making it difficult to identify them. TTFLs in fungi, animals, and plants often comprise proteins that post-translationally modify other proteins, including kinases and methyltransferases—examples are casein kinases, glycogen synthase kinase 3, and protein arginine-N-methyltransferase [70,71,72]. For diatoms, no orthologues of bacterial, animal, or plant circadian clock components have yet been found in the diatom genomes, except for CRYs and casein kinases [73], although the uptake of silicic acid and the general photophysiology [48] of diatoms depend on the diel cycle. The cell division of P. tricornutum is regulated by blue light, which is hypothesised to be the mechanism for synchronized cell division under light−dark cycles [44]. The recent reports on proteins such as bHLH-PAS (basic Helix-Loop-Helix-Per/Arnt/Sim) proteins, like RITMO1 (PtbHLH1a), that possess domain structures that are also found in clock components of animals [53,74], now open up the possibility to identify potential interaction partners. RITMO1 displays a robust diel gene transcript pattern, which is adjusted in a photoperiod-dependent manner, and the rhythmic expression persists in cells exposed to continuous light. The gene for the related PtbHLH1b shows a similar expression pattern as RITMO1 in a light−dark cycle, and thus might have a similar function [53]. Our data show that the genes of both PtbHLH1a and 1b proteins show a rhythmic expression in the light that persists in the dark at least for one day. PtbHLH1a and 1b genes lose their diel expression patterns in the PtAureo1a knockout mutant, which can be rescued in the PtAureo1a complemented mutant in LD conditions (Figure 2). Together with the finding that a knockout of PtAureo1a has a strong effect on the two PtbHLH proteins during a shift from red to blue light [43], it suggests that PtAureo1a is involved in the blue light-dependent regulation of RITMO1 and, thus, the diatom circadian clock.
It has been proposed that blue light photoreceptors may play a key role for aquatic photosynthetic organisms [49], because blue light penetrates deeper into the water column than light of longer wavelengths. In a metagenomic approach, Coesel et al. studied the diel transcription of genes in marine algae and found that, among Stramenopiles, the genes of aureochromes and other photoreceptors may respond to and anticipate fluctuating light conditions [50]. Therefore, we studied the diel gene expression of other blue light photoreceptors, i.e., cryptochrome/photolyase family 1 (CPF1) and Cry-DASH-like CPF (2 & 4) [60], and the role of PtAureo1a in regard to their expression. Consistent with Oliveri et al. [60], we found that all CPF transcripts peak during the day, but at different times and with different intensities. Our finding that the clear diel expression pattern of CPF1 is lost in the PtAureo1a knockout mutant and rescued in PtAureo1a complemented lines is indicative for a direct or indirect effect of PtAureo1a on CPF gene expression.
The demonstration that, in PtAureo1a KO mutants, the expression of a very large number of genes is affected as a response to both a shift from red to blue light [43] and also the diel expression (this study) may, in part, be attributed to the fact that aureochromes, in contrast to other photoreceptors, are both photoreceptors and transcription factors. This allows them to transmit their light responses directly to gene expression, avoiding intermediate phosphorylation steps of further messenger proteins. Accordingly, a strong effect on gene expression in response to blue light can already be observed after less than 10 min [43]. While a knockout of PtAureo1a may have a positive effect on the expression of PtCryP and a larger set of other photoreceptors [43], silencing of the PtCryP gene itself only changes the gene expression of diatom phytochrome PtDPH and PtCPF1 [75], again indicating a prominent role of PtAureo1a in the regulation of diel expression. We observed a significant correlation between the diel expression of the aureochrome isoforms with that of the PtbHLH1a/1b TFs. Specifically, the matrix analyses revealed a compelling correlation of PtAureo1a and 1b with PtbHLH1a/1b genes, as well as PtAureo1c under L/D conditions in WT, which is only detected in the presence of PtAureo1a. Taken together, our data indicate that PtAureo1a may control or adjust the diatom clock by controlling the diel expression of PtbHLH1a/1b. Furthermore, the expression of the gene of the flavin-lacking PtAureo2 protein correlates under L/D conditions with the PtbHLH1 genes only in the KO8 line, which indicates that PtAureo1a may function as a repressor of PtAureo2 diel expression.
As the most prominent known phenotype of PtAureo1a mutants is a change in photoacclimation [34,37,38], we studied the impact of PtAureo1a on the expression of the gene of the photoprotective PtLhcx1 protein [61]. Indeed, we could observe that the strongly light-dependent diel expression pattern of PtLhcx1 observed in WT and Co.48 lines is mostly absent in the PtAureo1a KO mutant. At the same time, we do not have any indication yet that Lhcx1 itself might be involved in clock regulation, but rather may be light-induced to facilitate dynamic tracking of light fluctuations in turbulent waters [76]. In line, beyond its control via PtAureo1a, PtLhcx1 expression responds to light intensity changes transduced via the redox state of the plastidic plastoquinone pool [62].
As a conclusion, our results indicate that PtAureo1a may be an important light-driven regulator of rhythmic gene expression in diatoms (Figure 5). In land plants, the time point of the transcript peaks can be very important for the regulation of clocks, as it reflects the synthesis of the respective proteins and its regulatory activity shortly thereafter. Coesel et al. [50] demonstrate that, among diatoms and related algal groups, individual aureochrome genes may be transcribed at night, as well as in the morning and up to noon, indicating that the aureochromes are ready as soon as light is available, stressing the role of aureochromes as “early” regulators. In P. tricornutum, we found the different aureochrome transcripts peak during the day, with PtAureo1a being the first of them being expressed. The investigations on PtAureo1a may represent an entrée into understanding the diatom clock and revealing unknown components of the clock in the environmentally important group of algae that evolved by secondary endosymbiosis. Future work will therefore have to focus on the subsequent processes initiated by PtAureo1a.

4. Materials and Methods

4.1. Experimental Set-Up and Growth Conditions

P. tricornutum (Bohlin) wild-type strain UTEX646 (Pt4) was initially obtained from the culture collection of algae of the University of Texas (UTEX, Austin, TX, USA, https://utex.org/). The TALEN mediated knockout strain PtAureo1aKO8 had been produced by [37], while its complemented strain PtAureo1aCo.48 is described in [38]. The cultivation for the day/night cycle experiments was carried out in photobioreactors (PBR102-S, Phenometrics, Inc., East Lansing, MI, USA) containing up to 400 mL of F/2 medium without added silica and at 50% sea water concentration medium. Cultures were grown at 20 °C with bubbling by using a mini diaphragm pump (VP86, VWR International Ltd., Lutterworth, UK), as well as by magnetic stirring at 300 rpm. The pH of the medium was adjusted to 7.4 before the start of the experiment using 1N HCl/NaOH. The white LED light (resembling full sunlight spectrum) was programmed to provide a continuous light of 100 μmol photons m−2 s−1 in a 16:8 h L/D cycle (light onset at 7 a.m. and light offset at 11 p.m.). For the continuous dark condition (D/D), the lights were programmed to switch off after acclimation to a 16:8 h L/D cycle. The cultures were maintained in their logarithmic phase at a Chl a concentration of approximately 1 μg/mL by regular automatic measurement of the optical density and accordingly variable dilution with fresh F/2 medium by the Phenometrics turbidostat pumps. Cells were acclimated to these conditions for at least 5–7 generations (5–7 day−night cycles) before sample collection for the experiment. Samples were harvested after dawn, starting at Day1 T8 (8 a.m.) in 4 h intervals until the 3rd day before dawn, i.e., at Day3 T4. Notably, 30 mL samples were taken under sterile conditions, frozen in liquid nitrogen and stored at −80 °C for further analyses. After the regular light−dark incubation, another 16:8 h light−dark cycle was added to allow cultures to reach sufficient cell numbers after sampling before switching off the light during the subjective day photoperiod. During the dark phase, samples were harvested similarly.
Along with the wild-type, the TALEN mediated knockout mutant strains PtAureo1aKO8 of P. tricornutum [37] and the complemented strain PtAureo1aCo.48 [38] were grown axenically in liquid F/2 medium without added silica and at 50% sea water concentration. Cells in liquid F/2 medium were cultivated in a 16 h/8 h light/dark cycle in Erlenmeyer flasks under continuous shaking at 20 °C and with an illumination of 45 μmol photons m−2 s−1 (Osram Lumilux L58W/840, Munich, Germany).

4.2. Determination of Growth Rates

Cell numbers were counted using a Multisizer 4 Coulter Counter (Beckman, Brea, CA, USA) and the growth curve experiment was set up with an equal number of cells (100,000 cells/mL). Cell number was measured over a period of 10 days in the bioreactor and without diluting with fresh medium to observe growth until the stationary phase.
The growth rate (µ) was calculated during the exponential phase using the following formula:
μ = ln N t ln N 0 Δ t
with:
  • ln Nt = Natural Logarithm of cell no. at day 5;
  • ln N0 = Natural logarithm of cell no. at day 3.
Δt = t1 − t0

4.3. Chlorophyll Fluorescence Measurements

Using an AquaPen-C AP100 (Photon Systems Instruments, Drásov, Czech Republic), mid-exponential phase diatom cultures were analysed to determine the maximum photosynthetic efficiency of PSII via the Fv/Fm ratio [Fv/Fm = (Fm − F0)/Fm], where Fm = maximum fluorescence measured in very low light adapted (30 min) cells induced by a saturating blue light pulse of 2100 μmol photons m−2 s−1 and F0 = minimum fluorescence of low light adapted cells.

4.4. RNA Isolation and cDNA Synthesis

RNA was isolated using Peqgold RNAPure and Peqgold Total RNA Kit S-Line with an on-column DNA digest using the Peqgold DNase I Digest kit (VWR, Darmstadt, Germany). RNA was extracted from samples normalized to 1 µg/mL Chlorophyll a, and concentration was measured spectrophotometrically using a Nanodrop 2000 UV/VIS Spectrophotometer (Thermo Fisher, Schwerte, Germany). Additionally, cDNA synthesis was performed using the Primescript kit (Takara Bio Europe, Saint-Germain-en-Laye, France).

4.5. Quantitative Real-Time PCR (qPCR)

We analysed the transcript abundance of 10 genes: PtAureo1a, PtAureo1b, PtAureo1c, PtAureo2, PtCPF1, PtCPF2, PtCPF4, PtbHLH1a, PtbHLH1b, and PtLhcx1 using HPLC grade primers (specified in Supplementary Tables S1 and S2) by qPCR. Additionally, 18S, TBP (TATA Binding Protein) and RPS (ribosomal protein S1) were used as reference genes for qPCR analysis. qPCR was run on a 7500 Fast RT-PCR system (Applied Biosystems, Waltham, MA, USA) and, for technical reasons, a few samples were run on CFX ConnectTM Real-Time system (Bio-Rad, Hercules, CA, USA). There were no significant differences between the obtained data from the two different devices. Samples were measured in biological triplicates as well as technical duplicates for each gene. Ct values (cycle threshold) and gene amplification efficiencies were calculated using PCR Miner 4.0 [77] from the obtained raw data. Relative mRNA transcript levels were calculated according to [78].

4.6. Analysis of Circadian Rhythm and Its Parameters

Circadian rhythm analysis was performed on time-series data using the online tool BioDare2 (https://biodare2.ed.ac.uk/) [64], based on the empirical JTK_CYCLE method [79]. The Kendal rank correlation coefficient (τ) was calculated in the range of [−1,1], where, τ = +1 implies a perfectly correlated series score, τ = −1 implies an anti-correlated series score, while τ = 0 implies an uncorrelated series score. The scores obtained were further optimized by applying a Benjamini−Hochberg correction (BH-factor) to reduce the false positive discovery rate (FDR). Further, a representative circadian rhythm curve is depicted in a particular measure that varies according to time. This difference between the crest and trough values is the amplitude of the rhythm. The timing of a reference point in the cycle (e.g., the peak) relative to a fixed event, such as the beginning of the night phase, is the phase parameter of the rhythm. In our case time point 0 h at midnight was chosen as the fixed event. The time difference between two peaks is called the period. Thereby, the parameters of circadian rhythm, i.e., amplitude, phase, and period, were calculated using two different methods available in the BioDare2 online tool [64]: 1. The Fast Fourier Transform Non-Linear Least Squares algorithm (FFT-NLLS), an algorithm that fits the sum of cosine functions on a sinusoidal time pattern; and 2. MESA (Maximum Entropy Spectral Analysis), an algorithm that applies a stochastic spectrum modelling approach.

4.7. Correlation Analysis

We used statistical tool “R” and the packages “corrplot” and “Hmisc” to perform the Pearson product moment correlation coefficient analysis (http://www.R-project.org/) of tested genes in order to hypothesize their diel expression network (script in the Supplementary Figure S2).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants13111465/s1, Table S1: Primers used in the study for RT-qPCR describing respective target genes and sequences, Table S2: JGI Protein IDs and Ensembl IDs of genes investigated in this study; Table S3: Pearson’s correlation coefficient (r) analysis of diurnal rhythm of test genes used in this study between 3 groups as: WT and PtAureo1a knock out mutant (KO8); KO8 and PtAureo1a complemented mutant (Co.48); WT and Co.48 under L/D condition; Figure S1: Western blot showing the diurnal abundance of PtAureo1 protein in P. tricornutum WT cells using anti-PtAureo1a antiserum and, as a loading control, D1 antiserum. Figure S2: R Script for Correlation analysis. R script used in this study to perform Correlation analysis on 10 test genes mRNA expression data under D/D condition for WT.

Author Contributions

Conceptualization, S.M. and P.G.K.; Methodology, S.M.; Software, S.M.; Validation, S.M., B.L. and P.G.K.; Formal Analysis, S.M., A.H.F., B.L. and P.G.K.; Investigation, S.M. and A.H.F.; Resources, P.G.K.; Data Curation, S.M. and B.L.; Writing, S.M., B.L. and P.G.K.; Visualization, S.M. and B.L.; Supervision, P.G.K.; Project Administration, P.G.K.; Funding Acquisition, S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Deutscher Akademischer Austauschdienst (DAAD, to S.M.), the Deutsche Forschungsgemeinschaft (DFG, projects KR 1661-18, KR 1661-20 to P.G.K., and LE 3358/3 and Heisenberg programme to B.L.).

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Acknowledgments

We would like to acknowledge the financial support by the University of Konstanz. The work shown here is a modified part of the dissertation by Shvaita Madhuri “Elucidation of functional roles of blue light receptors: Aureochromes in Phaeodactylum tricornutum” (Universität Konstanz, 2020, https://kops.uni-konstanz.de/entities/publication/ab5f4f58-333c-4f39-bdca-74fcac0b090c) [80].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Falkowski, P.G.; Katz, M.E.; Knoll, A.H.; Quigg, A.; Raven, J.A.; Schofield, O.; Taylor, F.J. The evolution of modern eukaryotic phytoplankton. Science 2004, 305, 354–360. [Google Scholar] [CrossRef] [PubMed]
  2. Karlusich JJ, P.; Ibarbalz, F.M.; Bowler, C. Phytoplankton in the Tara Ocean. Annu. Rev. Mar. Sci. 2020, 12, 233–265. [Google Scholar] [CrossRef]
  3. Archibald, J.M. Endosymbiosis and Eukaryotic Cell Evolution. Curr. Biol. 2015, 25, R911–R921. [Google Scholar] [CrossRef] [PubMed]
  4. Burki, F.; Roger, A.J.; Brown, M.W.; Simpson, A.G.B. The New Tree of Eukaryotes. Trends Ecol. Evol. 2020, 35, 43–55. [Google Scholar] [CrossRef] [PubMed]
  5. Moustafa, A.; Beszteri, B.; Maier, U.G.; Bowler, C.; Valentin, K.; Bhattacharya, D. Genomic Footprints of a Cryptic Plastid Endosymbiosis in Diatoms. Science 2009, 324, 1724–1726. [Google Scholar] [CrossRef] [PubMed]
  6. Bowler, C.; Allen, A.E.; Badger, J.H.; Grimwood, J.; Jabbari, K.; Kuo, A.; Maheswari, U.; Martens, C.; Maumus, F.; Otillar, R.P.; et al. The Phaeodactylum genome reveals the evolutionary history of diatom genomes. Nature 2008, 456, 239–244. [Google Scholar] [CrossRef] [PubMed]
  7. Flori, S.; Jouneau, P.-H.; Bailleul, B.; Gallet, B.; Estrozi, L.F.; Moriscot, C.; Bastien, O.; Eicke, S.; Schober, A.; Bártulos, C.R.; et al. Plastid thylakoid architecture optimizes photosynthesis in diatoms. Nat. Commun. 2017, 8, 15885. [Google Scholar] [CrossRef] [PubMed]
  8. Allen, A.E.; Dupont, C.L.; Obornik, M.; Horak, A.; Nunes-Nesi, A.; McCrow, J.P.; Zheng, H.; Johnson, D.A.; Hu, H.; Fernie, A.R.; et al. Evolution and metabolic significance of the urea cycle in photosynthetic diatoms. Nature 2011, 473, 203–207. [Google Scholar] [CrossRef] [PubMed]
  9. Depauw, F.A.; Rogato, A.; Riberad, A.M.; Falciatore, A. Exploring the molecular basis of responses to light in marine diatoms. J. Exp. Bot. 2012, 63, 1575–1591. [Google Scholar] [CrossRef]
  10. Bailleul, B.; Berne, N.; Murik, O.; Petroutsos, D.; Prihoda, J.; Tanaka, A.; Villanova, V.; Bligny, R.; Flori, S.; Falconet, D.; et al. Energetic coupling between plastids and mitochondria drives CO2 assimilation in diatoms. Nature 2015, 524, 366–369. [Google Scholar] [CrossRef]
  11. Duanmu, D.; Rockwell, N.C.; Clark Lagarias, J. Algal light sensing and photoacclimation in aquatic environments. Plant Cell Environ. 2017, 40, 2558–2570. [Google Scholar] [CrossRef]
  12. Oakenfull, R.J.; Davis, S.J. Shining a light on the Arabidopsis circadian clock. Plant Cell Environ. 2017, 40, 2571–2585. [Google Scholar] [CrossRef] [PubMed]
  13. Kianianmomeni, A.; Hallmann, A. Algal photoreceptors: In vivo functions and potential applications. Planta 2014, 239, 1–26. [Google Scholar] [CrossRef] [PubMed]
  14. Litchman, E. Growth rates of phytoplankton under fluctuating light. Freshw. Biol. 2000, 44, 223–235. [Google Scholar] [CrossRef]
  15. Wagner, H.; Jakob, T.; Wilhelm, C. Balancing the energy flow from captured light to biomass under fluctuating light conditions. New Phytol. 2006, 169, 95–108. [Google Scholar] [CrossRef] [PubMed]
  16. Lepetit, B.; Gélin, G.; Lepetit, M.; Sturm, S.; Vugrinec, S.; Rogato, A.; Kroth, P.G.; Falciatore, A.; Lavaud, J. The diatom Phaeodactylum tricornutum adjusts nonphotochemical fluorescence quenching capacity in response to dynamic light via fine-tuned Lhcx and xanthophyll cycle pigment synthesis. New Phytol. 2017, 214, 205–218. [Google Scholar] [CrossRef] [PubMed]
  17. Orefice, I.; Chandrasekaran, R.; Smerilli, A.; Corato, F.; Caruso, T.; Casillo, A.; Corsaro, M.M.; Piaz, F.D.; Ruban, A.V.; Brunet, C. Light-induced changes in the photosynthetic physiology and biochemistry in the diatom Skeletonema marinoi. Algal Res. 2016, 17, 1–13. [Google Scholar] [CrossRef]
  18. Huisman, J.; Jonker, R.R.; Zonneveld, C.; Weissing, F.J. Competition for light between phtoplankton species: Experimental tests of mechanistic theory. Ecology 1999, 80, 211–222. [Google Scholar] [CrossRef]
  19. Litchman, E.; Klausmeier, C.A. Competition of phytoplankton under fluctuating light. Am. Nat. 2001, 157, 170–187. [Google Scholar] [CrossRef]
  20. Litchman, E.; Klausmeier, C.A.; Bossard, P. Phytoplankton nutrient competition under dynamic light regimes. Limnol. Oceanogr. 2004, 49, 1457–1462. [Google Scholar] [CrossRef]
  21. Villareal, T.A.; Pilskaln, C.; Brzezinski, M.; Lipschultz, F.; Dennett, M.; Gardner, G.B. Upward transport of oceanic nitrate by migrating diatom mats. Nature 1999, 397, 423–425. [Google Scholar] [CrossRef]
  22. McLachlan, D.H.; Brownlee, C.; Taylor, A.R.; Geider, R.J.; Underwood, G.J.C. Light-induced motile responses of the estuarine benthic diatoms Navicula perminuta and Cylindrotheca closterium (Bacillariophyceae). J. Phycol. 2009, 45, 592–599. [Google Scholar] [CrossRef] [PubMed]
  23. Shihira-Ishikawa, I.; Nakamura, T.; Higashi, S.; Watanabe, M. Distinct responses of chloroplasts to blue and green laser microbeam irradiations in the centric diatom Pleurosira laevis. Photochem. Photobiol. 2007, 83, 1101–1109. [Google Scholar] [CrossRef] [PubMed]
  24. Fankhauser, C.; Staiger, D. Photoreceptors in Arabidopsis thaliana: Light perception, signal transduction and entrainment of the endogenous clock. Planta 2002, 216, 1–16. [Google Scholar] [CrossRef] [PubMed]
  25. Schmalen, I.; Reischl, S.; Wallach, T.; Klemz, R.; Grudziecki, A.; Prabu, J.R.; Benda, C.; Kramer, A.; Wolf, E. Interaction of Circadian Clock Proteins CRY1 and PER2 Is Modulated by Zinc Binding and Disulfide Bond Formation. Cell 2014, 157, 1203–1215. [Google Scholar] [CrossRef] [PubMed]
  26. Kottke, T.; Oldemeyer, S.; Wenzel, S.; Zou, Y.; Mittag, M. Cryptochrome photoreceptors in green algae: Unexpected versatility of mechanisms and functions. J. Plant Physiol. 2017, 217, 4–14. [Google Scholar] [CrossRef] [PubMed]
  27. Noordally, Z.B.; Millar, A.J. Clocks in Algae. Biochemistry 2015, 54, 171–183. [Google Scholar] [CrossRef] [PubMed]
  28. Chen, M.; Chory, J. Phytochrome signaling mechanisms and the control of plant development. Trends Cell Biol. 2011, 21, 664–671. [Google Scholar] [CrossRef] [PubMed]
  29. Jaubert, M.; Duchêne, C.; Kroth, P.G.; Rogato, A.; Bouly, J.-P.; Falciatore, A. Sensing and Signalling in Diatom Responses to Abiotic Cues. In The Molecular Life of Diatoms; Falciatore, A., Mock, T., Eds.; Springer International Publishing: Cham, Switzerland, 2022; pp. 607–639. [Google Scholar]
  30. Takahashi, F.; Yamagata, D.; Ishikawa, M.; Fukamatsu, Y.; Ogura, Y.; Kasahara, M.; Kiyosue, T.; Kikuyama, M.; Wada, M.; Kataoka, H. Aureochrome, a photoreceptor required for photomorphogenesis in stramenopiles. Proc. Natl. Acad. Sci. USA 2007, 104, 19625–19630. [Google Scholar] [CrossRef]
  31. Kroth, P.G.; Wilhelm, C.; Kottke, T. An update on aureochromes: Phylogeny–mechanism–function. J. Plant Physiol. 2017, 217, 20–26. [Google Scholar] [CrossRef]
  32. Rodríguez-Martínez, J.A.; Reinke, A.W.; Bhimsaria, D.; Keating, A.E.; Ansari, A.Z. Combinatorial bZIP dimers display complex DNA-binding specificity landscapes. eLife 2017, 6, e19272. [Google Scholar] [CrossRef] [PubMed]
  33. Ishikawa, M.; Takahashi, F.; Nozaki, H.; Nagasato, C.; Motomura, T.; Kataoka, H. Distribution and phylogeny of the blue light receptors aureochromes in eukaryotes. Planta 2009, 230, 543–552. [Google Scholar] [CrossRef] [PubMed]
  34. Schellenberger Costa, B.; Sachse, M.; Jungandreas, A.; Río Bártulos, C.; Gruber, A.; Jakob, T.; Kroth, P.G.; Wilhelm, C. Aureochrome 1a is involved in the photoacclimation of the diatom Phaeodactylum tricornutum. PLoS ONE 2013, 8, e74451. [Google Scholar] [CrossRef] [PubMed]
  35. Grattepanche, J.-D.; Walker, L.M.; Ott, B.M.; Paim Pinto, D.L.; Delwiche, C.F.; Lane, C.E.; Katz, L.A. Microbial Diversity in the Eukaryotic SAR Clade: Illuminating the Darkness Between Morphology and Molecular Data. BioEssays 2018, 40, 1700198. [Google Scholar] [CrossRef] [PubMed]
  36. Schellenberger, C.B.; Jungandreas, A.; Jakob, T.; Weisheit, W.; Mittag, M.; Wilhelm, C. Blue light is essential for high light acclimation and photoprotection in the diatom Phaeodactylum tricornutum. J. Exp. Bot. 2012, 64, 483–493. [Google Scholar] [CrossRef] [PubMed]
  37. Serif, M.; Lepetit, B.; Weißert, K.; Kroth, P.G.; Río Bártulos, C. A fast and reliable strategy to generate TALEN-mediated gene knockouts in the diatom Phaeodactylum tricornutum. Algal Res. 2017, 23, 186–195. [Google Scholar] [CrossRef]
  38. Madhuri, S.; Río Bártulos, C.; Serif, M.; Lepetit, B.; Kroth, P.G. A strategy to complement PtAUREO1a in TALEN knockout strains of Phaeodactylum tricornutum. Algal Res. 2019, 39, 101469. [Google Scholar] [CrossRef]
  39. Banerjee, A.; Herman, E.; Serif, M.; Maestre-Reyna, M.; Hepp, S.; Pokorny, R.; Kroth, P.G.; Essen, L.-O.; Kottke, T. Allosteric communication between DNA-binding and light-responsive domains of diatom class I aureochromes. Nucleic Acids Res. 2016, 44, 5957–5970. [Google Scholar] [CrossRef] [PubMed]
  40. Goett-Zink, L.; Klocke, J.L.; Bögeholz LA, K.; Kottke, T. In-cell infrared difference spectroscopy of LOV photoreceptors reveals structural responses to light altered in living cells. J. Biol. Chem. 2020, 295, 11729–11741. [Google Scholar] [CrossRef]
  41. Herman, E.; Kottke, T. Allosterically Regulated Unfolding of the A′α Helix Exposes the Dimerization Site of the Blue-Light-Sensing Aureochrome-LOV Domain. Biochemistry 2015, 54, 1484–1492. [Google Scholar] [CrossRef]
  42. Tian, H.; Trozzi, F.; Zoltowski, B.D.; Tao, P. Deciphering the Allosteric Process of the Phaeodactylum tricornutum Aureochrome 1a LOV Domain. J. Phys. Chem. B 2020, 41, 8960–8972. [Google Scholar] [CrossRef] [PubMed]
  43. Mann, M.; Serif, M.; Wrobel, T.; Eisenhut, M.; Madhuri, S.; Flachbart, S.; Weber AP, M.; Lepetit, B.; Wilhelm, C.; Kroth, P.G. The Aureochrome Photoreceptor PtAUREO1a Is a Highly Effective Blue Light Switch in Diatoms. iScience 2020, 23, 101730. [Google Scholar] [CrossRef] [PubMed]
  44. Huysman, M.J.; Fortunato, A.E.; Matthijs, M.; Costa, B.S.; Vanderhaeghen, R.; Van den Daele, H.; Sachse, M.; Inze, D.; Bowler, C.; Kroth, P.G.; et al. AUREOCHROME1a-mediated induction of the diatom-specific cyclin dsCYC2 controls the onset of cell division in diatoms (Phaeodactylum tricornutum). Plant Cell 2013, 25, 215–228. [Google Scholar] [CrossRef] [PubMed]
  45. Im, S.H.; Lepetit, B.; Mosesso, N.; Shrestha, S.; Weiss, L.; Nymark, M.; Roellig, R.; Wilhelm, C.; Isono, E.; Kroth, P.G. Identification of promoter targets by Aureochrome 1a in the diatom Phaeodactylum tricornutum. J. Exp. Bot. 2023, 75, 1834–1851. [Google Scholar] [CrossRef] [PubMed]
  46. Farré, E.M. The brown clock: Circadian rhythms in stramenopiles. Physiol. Plant. 2020, 169, 430–441. [Google Scholar] [CrossRef] [PubMed]
  47. Moulager, M.; Monnier, A.; Jesson, B.; Bouvet, R.; Mosser, J.; Schwartz, C.; Garnier, L.; Corellou, F.; Bouget, F.O.-Y. Light-Dependent Regulation of Cell Division in Ostreococcus: Evidence for a Major Transcriptional Input. Plant Physiol. 2007, 144, 1360–1369. [Google Scholar] [CrossRef] [PubMed]
  48. Ragni, M.; d’Alcala, M.R. Circadian variability in the photobiology of Phaeodactylum tricornutum: Pigment content. J. Plankt Res. 2007, 29, 141–156. [Google Scholar] [CrossRef]
  49. Poliner, E.; Cummings, C.; Newton, L.; Farré, E.M. Identification of circadian rhythms in Nannochloropsis species using bioluminescence reporter lines. Plant J. 2019, 99, 112–127. [Google Scholar] [CrossRef]
  50. Coesel, S.N.; Durham, B.P.; Groussman, R.D.; Hu, S.K.; Caron, D.A.; Morales, R.L.; Ribalet, F.; Armbrust, E.V. Diel transcriptional oscillations of light-sensitive regulatory elements in open-ocean eukaryotic plankton communities. Proc. Natl. Acad. Sci. USA 2021, 118, e2011038118. [Google Scholar] [CrossRef]
  51. Ashworth, J.; Coesel, S.; Lee, A.; Armbrust, E.V.; Orellana, M.V.; Baliga, N.S. Genome-wide diel growth state transitions in the diatom Thalassiosira pseudonana. Proc. Natl. Acad. Sci. USA 2013, 110, 7518–7523. [Google Scholar] [CrossRef]
  52. Chauton, M.S.; Winge, P.; Brembu, T.; Vadstein, O.; Bones, A.M. Gene regulation of carbon fixation, storage and utilization in the diatom Phaeodactylum tricornutum acclimated to light/dark cycles. Plant Physiol. 2012, 161, 1034–1048. [Google Scholar] [CrossRef] [PubMed]
  53. Annunziata, R.; Ritter, A.; Fortunato, A.E.; Manzotti, A.; Cheminant-Navarro, S.; Agier, N.; Huysman MJ, J.; Winge, P.; Bones, A.M.; Bouget, F.-Y.; et al. bHLH-PAS protein RITMO1 regulates diel biological rhythms in the marine diatom Phaeodactylum tricornutum. Proc. Natl. Acad. Sci. USA 2019, 116, 13137–13142. [Google Scholar] [CrossRef] [PubMed]
  54. Fowler, S.; Lee, K.; Onouchi, H.; Samach, A.; Richardson, K.; Morris, B.; Coupland, G.; Putterill, J. GIGANTEA: A circadian clock-controlled gene that regulates photoperiodic flowering in Arabidopsis and encodes a protein with several possible membrane-spanning domains. EMBO J. 1999, 18, 4679–4688. [Google Scholar] [CrossRef] [PubMed]
  55. Kim, S.-G.; Kim, S.-Y.; Park, C.-M. A membrane-associated NAC transcription factor regulates salt-responsive flowering via FLOWERING LOCUS T in Arabidopsis. Planta 2007, 226, 647–654. [Google Scholar] [CrossRef] [PubMed]
  56. Sawa, M.; Nusinow, D.A.; Kay, S.A.; Imaizumi, T. FKF1 and GIGANTEA Complex Formation Is Required for Day-Length Measurement in Arabidopsis. Science 2007, 318, 261–265. [Google Scholar] [CrossRef] [PubMed]
  57. Coesel, S.; Mangogna, M.; Ishikawa, T.; Heijde, M.; Rogato, A.; Finazzi, G.; Todo, T.; Bowler, C.; Falciatore, A. Diatom PtCPF1 is a new cryptochrome/photolyase family member with DNA repair and transcription regulation activity. EMBO Rep. 2009, 10, 655–661. [Google Scholar] [CrossRef]
  58. Matthijs, M.; Fabris, M.; Broos, S.; Vyverman, W.; Goossens, A. Profiling of the early nitrogen stress response in the diatom Phaeodactylum tricornutum reveals a novel family of RING-domain transcription factors. Plant Physiol. 2016, 170, 489–498. [Google Scholar] [CrossRef] [PubMed]
  59. Im, S.H.; Madhuri, S.; Lepetit, B.; Kroth, P.G. Functional demonstration of Aureochrome 1a proteasomal degradation after blue light incubation in the diatom Phaeodactylum tricornutum. J. Plant Physiol. 2024, 292, 154148. [Google Scholar] [CrossRef] [PubMed]
  60. Oliveri, P.; Fortunato, A.E.; Petrone, L.; Ishikawa-Fujiwara, T.; Kobayashi, Y.; Todo, T.; Antonova, O.; Arboleda, E.; Zantke, J.; Tessmar-Raible, K.; et al. The Cryptochrome/Photolyase family in aquatic organisms. Mar. Genom. 2014, 14, 23–37. [Google Scholar] [CrossRef]
  61. Buck, J.M.; Sherman, J.; Río Bártulos, C.; Serif, M.; Halder, M.; Henkel, J.; Falciatore, A.; Lavaud, J.; Gorbunov, M.Y.; Kroth, P.G.; et al. Lhcx proteins provide photoprotection via thermal dissipation of absorbed light in the diatom Phaeodactylum tricornutum. Nat. Commun. 2019, 10, 4167. [Google Scholar] [CrossRef]
  62. Lepetit, B.; Sturm, S.; Rogato, A.; Gruber, A.; Sachse, M.; Falciatore, A.; Kroth, P.G.; Lavaud, J. High light acclimation in the secondary plastids containing diatom Phaeodactylum tricornutum is triggered by the redox state of the plastoquinone pool. Plant Physiol. 2013, 161, 853–865. [Google Scholar] [CrossRef] [PubMed]
  63. Straume, M.; Frasier-Cadoret, S.G.; Johnson, M.L. Least-Squares Analysis of Fluorescence Data. In Topics in Fluorescence Spectroscopy: Principles; Lakowicz, J.R., Ed.; Springer: Boston, MA, USA, 2002; pp. 177–240. [Google Scholar]
  64. Zielinski, T.; Moore, A.M.; Troup, E.; Halliday, K.J.; Millar, A.J. Strengths and Limitations of Period Estimation Methods for Circadian Data. PLoS ONE 2014, 9, e96462. [Google Scholar] [CrossRef]
  65. Kolody, B.C.; McCrow, J.P.; Allen, L.Z.; Aylward, F.O.; Fontanez, K.M.; Moustafa, A.; Moniruzzaman, M.; Chavez, F.P.; Scholin, C.A.; Allen, E.E.; et al. Diel transcriptional response of a California Current plankton microbiome to light, low iron, and enduring viral infection. ISME J. 2019, 13, 2817–2833. [Google Scholar] [CrossRef]
  66. Saini, R.; Jaskolski, M.; Davis, S.J. Circadian oscillator proteins across the kingdoms of life: Structural aspects. BMC Biol. 2019, 17, 13. [Google Scholar] [CrossRef]
  67. Creux, N.; Harmer, S. Circadian Rhythms in Plants. Cold Spring Harb. Perspect. Biol. 2019, 11, a034611. [Google Scholar] [CrossRef]
  68. Forbes-Stovall, J.; Howton, J.; Young, M.; Davis, G.; Chandler, T.; Kessler, B.; Rinehart, C.A.; Jacobshagen, S. Chlamydomonas reinhardtii strain CC-124 is highly sensitive to blue light in addition to green and red light in resetting its circadian clock, with the blue-light photoreceptor plant cryptochrome likely acting as negative modulator. Plant Physiol. Biochem. 2014, 75, 14–23. [Google Scholar] [CrossRef]
  69. Brunner, M.; Káldi, K. Interlocked feedback loops of the circadian clock of Neurospora crassa. Mol. Microbiol. 2008, 68, 255–262. [Google Scholar] [CrossRef] [PubMed]
  70. Hong, S.; Song, H.-R.; Lutz, K.; Kerstetter, R.A.; Michael, T.P.; McClung, C.R. Type II protein arginine methyltransferase 5 (PRMT5) is required for circadian period determination in Arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 2010, 107, 21211–21216. [Google Scholar] [CrossRef] [PubMed]
  71. Lu, S.X.; Knowles, S.M.; Webb, C.J.; Celaya, R.B.; Cha, C.; Siu, J.P.; Tobin, E.M. The Jumonji C Domain-Containing Protein JMJ30 Regulates Period Length in the Arabidopsis Circadian Clock. Plant Physiol. 2010, 155, 906–915. [Google Scholar] [CrossRef] [PubMed]
  72. Ode, K.L.; Ueda, H.R. Design Principles of Phosphorylation-Dependent Timekeeping in Eukaryotic Circadian Clocks. Cold Spring Harb. Perspect. Biol. 2018, 10, a028357. [Google Scholar] [CrossRef]
  73. Chisholm, S.W.; Azam, F.; Eppley, R.W. Silicic acid incorporation in marine diatoms on light: Dark cycles: Use as an assay for phased cell division. Limnol. Oceanogr. 1978, 23, 518–529. [Google Scholar] [CrossRef]
  74. Crane, B.R.; Young, M.W. Interactive Features of Proteins Composing Eukaryotic Circadian Clocks. Annu. Rev. Biochem. 2014, 83, 191–219. [Google Scholar] [CrossRef] [PubMed]
  75. König, S.; Eisenhut, M.; Bräutigam, A.; Kurz, S.; Weber, A.P.M.; Büchel, C. The Influence of a Cryptochrome on the Gene Expression Profile in the Diatom Phaeodactylum tricornutum under Blue Light and in Darkness. Plant Cell Physiol. 2017, 58, 1914–1923. [Google Scholar] [CrossRef] [PubMed]
  76. Giovagnetti, V.; Jaubert, M.; Shukla, M.K.; Ungerer, P.; Bouly, J.-P.; Falciatore, A.; Ruban, A.V. Biochemical and molecular properties of LHCX1, the essential regulator of dynamic photoprotection in diatoms. Plant Physiol. 2021, 188, 509–525. [Google Scholar] [CrossRef] [PubMed]
  77. Zhao, S.; Fernald, R.D. Comprehensive Algorithm for Quantitative Real-Time Polymerase Chain Reaction. J. Comput. Biol. 2005, 12, 1047–1064. [Google Scholar] [CrossRef] [PubMed]
  78. Pfaffl, M.W. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29, e45. [Google Scholar] [CrossRef]
  79. Hutchison, A.L.; Maienschein-Cline, M.; Chiang, A.H.; Tabei SM, A.; Gudjonson, H.; Bahroos, N.; Allada, R.; Dinner, A.R. Improved Statistical Methods Enable Greater Sensitivity in Rhythm Detection for Genome-Wide Data. PLoS Comput. Biol. 2015, 11, e1004094. [Google Scholar] [CrossRef]
  80. Madhuri, S. Elucidation of Functional Roles of Blue Light Receptors: Aureochromes in Phaeodactylum tricornutum. Ph.D. Thesis, Universität Konstanz, Konstanz, Germany, 2020. Available online: https://kops.uni-konstanz.de/entities/publication/ab5f4f58-333c-4f39-bdca-74fcac0b090c (accessed on 26 October 2020).
Figure 1. (A) Comparison of growth characteristics and (B) maximum photosystem II yield (Fv/Fm) of wild-type (WT), PtAureo1a knockout (KO8) and PtAureo1a8−complemented mutant (Co.48) during their exponential phase. Indicated values represent the arithmetic means (n = 3 for growth rate, n = 5–11 for Fv/Fm), and the error bars represent the standard error. A one-way ANOVA was performed, showing no statistical difference neither for growth (p = 0.154) nor for maximum PSII yield (p = 0.083) between the three strains.
Figure 1. (A) Comparison of growth characteristics and (B) maximum photosystem II yield (Fv/Fm) of wild-type (WT), PtAureo1a knockout (KO8) and PtAureo1a8−complemented mutant (Co.48) during their exponential phase. Indicated values represent the arithmetic means (n = 3 for growth rate, n = 5–11 for Fv/Fm), and the error bars represent the standard error. A one-way ANOVA was performed, showing no statistical difference neither for growth (p = 0.154) nor for maximum PSII yield (p = 0.083) between the three strains.
Plants 13 01465 g001
Figure 2. Diel expression of aureochrome genes (PtAureo1a, 1b, 1c, 2), cryptochrome photolyase family (PtCPF1, 2, 4), Lhcx1, and PtbHLH1a, 1b transcription factors in P. tricornutum WT (black), PtAureo1aKO8 knock out mutant (KO8, red), and PtAureo1a8Co.48−complemented mutant (Co.48, dotted blue line) under 16:8 h L/D and D/D conditions. The relative gene expression is represented as mean Log2 value with SEM (Standard error mean), n = 3 (WT & KO8); n = 2 (Co.48), normalized on 18S, TBP and RPS. The white and grey regions represent light and dark periods, respectively. The numbers represent the time points of sampling with respect to each day. Note that the y-scale is equidistant with respect to the relative log2 values for all genes.
Figure 2. Diel expression of aureochrome genes (PtAureo1a, 1b, 1c, 2), cryptochrome photolyase family (PtCPF1, 2, 4), Lhcx1, and PtbHLH1a, 1b transcription factors in P. tricornutum WT (black), PtAureo1aKO8 knock out mutant (KO8, red), and PtAureo1a8Co.48−complemented mutant (Co.48, dotted blue line) under 16:8 h L/D and D/D conditions. The relative gene expression is represented as mean Log2 value with SEM (Standard error mean), n = 3 (WT & KO8); n = 2 (Co.48), normalized on 18S, TBP and RPS. The white and grey regions represent light and dark periods, respectively. The numbers represent the time points of sampling with respect to each day. Note that the y-scale is equidistant with respect to the relative log2 values for all genes.
Plants 13 01465 g002
Figure 3. Pearson’s correlation matrix analysis of diel expression under L/D conditions of genes of aureochromes (PtAureo1a, b, c, 2), cryptochrome/photolyase family protein (PtCPF1), Cry−DASH−like family protein (PtCPF2 & 4), PtbHLH−PAS TFs (PtbHLH1a and 1b) and photoprotection protein PtLhcx1 of wild−type (WT), PtAureo1a knockout (KO8) and PtAureo1a−complemented mutant (Co.48): The lower correlograms show correlations with a significance level of an alpha value <0.05 where insignificant correlations are shown as blank boxes in the correlogram. The circle size indicates strength of correlation coefficient, r value (small circle ≤ ±0.2; medium circle ≤ ±0.5; normal circle ≤ ±1); for detailed values see Table S3. The matrix analysis was performed using R, packages corrplot and Hmisc (see Section 4).
Figure 3. Pearson’s correlation matrix analysis of diel expression under L/D conditions of genes of aureochromes (PtAureo1a, b, c, 2), cryptochrome/photolyase family protein (PtCPF1), Cry−DASH−like family protein (PtCPF2 & 4), PtbHLH−PAS TFs (PtbHLH1a and 1b) and photoprotection protein PtLhcx1 of wild−type (WT), PtAureo1a knockout (KO8) and PtAureo1a−complemented mutant (Co.48): The lower correlograms show correlations with a significance level of an alpha value <0.05 where insignificant correlations are shown as blank boxes in the correlogram. The circle size indicates strength of correlation coefficient, r value (small circle ≤ ±0.2; medium circle ≤ ±0.5; normal circle ≤ ±1); for detailed values see Table S3. The matrix analysis was performed using R, packages corrplot and Hmisc (see Section 4).
Plants 13 01465 g003
Figure 4. Pearson’s correlation matrix analysis of diel expression of genes of aureochromes (PtAureo1a, b, c, 2) under D/D conditions; cryptochrome/photolyase family protein (PtCPF1), Cry−DASH-like family protein (PtCPF2 & 4), bHLH−PAS TFs (PtbHLH1a and 1b) and photoprotection protein PtLhcx1 of wild−type (WT) P. tricornutum and PtAureo1a knock out mutant (KO8) (Upper). Correlogram showing correlation with significance level of an alpha value < 0.05 where insignificant correlations are left with blank boxes in the correlogram (lower figures). The circle size indicates strength of correlation coefficient, r value (small circle ≤ ±0.2; medium circle ≤ ±0.5; normal circle ≤ ±1); for detailed values see Table S3. The matrix analysis was performed using R, packages corrplot and Hmisc (see Section 4).
Figure 4. Pearson’s correlation matrix analysis of diel expression of genes of aureochromes (PtAureo1a, b, c, 2) under D/D conditions; cryptochrome/photolyase family protein (PtCPF1), Cry−DASH-like family protein (PtCPF2 & 4), bHLH−PAS TFs (PtbHLH1a and 1b) and photoprotection protein PtLhcx1 of wild−type (WT) P. tricornutum and PtAureo1a knock out mutant (KO8) (Upper). Correlogram showing correlation with significance level of an alpha value < 0.05 where insignificant correlations are left with blank boxes in the correlogram (lower figures). The circle size indicates strength of correlation coefficient, r value (small circle ≤ ±0.2; medium circle ≤ ±0.5; normal circle ≤ ±1); for detailed values see Table S3. The matrix analysis was performed using R, packages corrplot and Hmisc (see Section 4).
Plants 13 01465 g004
Figure 5. Scheme depicting and summarizing the obtained results on correlated expression of genes and the dependence in PtAureo1a, including the other aureochrome isoforms (PtAureo1b, c and PtAureo2), cryptochromes (PtCPF1 and 4) and PtbHLH1a/b). PtAureo1a might directly/indirectly influence the expression of the diatom clock component RITMO1 as well as the photoreceptors PtCPF1 and PtCPF4. Blue lines represent results from biological experiment (qPCR analysis), while grey dotted line represent results from Pearson’s correlation analysis (statistical).
Figure 5. Scheme depicting and summarizing the obtained results on correlated expression of genes and the dependence in PtAureo1a, including the other aureochrome isoforms (PtAureo1b, c and PtAureo2), cryptochromes (PtCPF1 and 4) and PtbHLH1a/b). PtAureo1a might directly/indirectly influence the expression of the diatom clock component RITMO1 as well as the photoreceptors PtCPF1 and PtCPF4. Blue lines represent results from biological experiment (qPCR analysis), while grey dotted line represent results from Pearson’s correlation analysis (statistical).
Plants 13 01465 g005
Table 1. Rhythm analysis of transcript abundances determined via qRT-PCR by amplifying transcripts of PtAureochromes (1a, b, c and 2), PtCPF blue-light photoreceptors (1, 2 and 4), bHLH-PAS proteins (PtbHLH1a and b), and the PtLhcx1 gene from WT, KO8, and Co.48 cultivated in photobioreactors under L/D conditions. Kendal rank correlation coefficients (τ) were calculated using empirical JTK_CYCLE method from the online tool BioDare2. τ = 1 would state a perfectly correlated time-series. Statistical significance (p-value) is represented as * < 0.1; ** < 0.05; *** < 0.01; **** < 0.005 and ***** < 0.001 with Benjamini−Hochberg correction.
Table 1. Rhythm analysis of transcript abundances determined via qRT-PCR by amplifying transcripts of PtAureochromes (1a, b, c and 2), PtCPF blue-light photoreceptors (1, 2 and 4), bHLH-PAS proteins (PtbHLH1a and b), and the PtLhcx1 gene from WT, KO8, and Co.48 cultivated in photobioreactors under L/D conditions. Kendal rank correlation coefficients (τ) were calculated using empirical JTK_CYCLE method from the online tool BioDare2. τ = 1 would state a perfectly correlated time-series. Statistical significance (p-value) is represented as * < 0.1; ** < 0.05; *** < 0.01; **** < 0.005 and ***** < 0.001 with Benjamini−Hochberg correction.
Proteins and RoleGene NameKendall Rank Correlation Coefficient (τ)
L/D
WTKO8Co.48
Aureochrome IsoformsPtAureo1a0.64 * 0.66 *
PtAureo1b0.590.360.71 **
PtAureo1c0.73 **0.520.78 **
PtAureo 20.540.40.64 **
Cryptochrome/photolyase family protein1PtCPF10.81 ****0.40.74 **
Cry-DASH-like protein familyPtCPF20.460.590.72 **
PtCPF40.64 *0.560.83 ****
Circadian clock-relatedPtbHLH1a (Ritmo1) 0.88 ****0.420.92 *****
PtbHLH1b0.88 ****0.460.89 ****
PhotoprotectionPtLhcx10.87 ****0.380.9 ****
Table 2. Period (h), phase (h) and amplitude analysis of PtCPF1, PtLhcx1, PtbHLH1a and 1b in WT, KO8 and Co.48 under 16:8 h (L/D). The BioDare2 online tool was used to apply FFT-NLLS and MESA analyses to calculate period and phase estimates.
Table 2. Period (h), phase (h) and amplitude analysis of PtCPF1, PtLhcx1, PtbHLH1a and 1b in WT, KO8 and Co.48 under 16:8 h (L/D). The BioDare2 online tool was used to apply FFT-NLLS and MESA analyses to calculate period and phase estimates.
PeriodPhaseAmplitude
FFT-NLLSMESAFFT-NLLSMESAFFT-NLLSMESA
CPF1WT26.326.15.45.71.51.5
KO831.434.210.68.70.60.7
Co.4824.325.17.66.71.11.1
bHLH1aWT2524.816.616.82.32.2
KO8no fit26.7no fit18.3no fit0.6
Co.4823.423.819.919.32.42.3
bHLH1bWT24.424.915.414.82.22.1
KO831.525.65.211.20.50.4
Co.4825.325.315.1152.42.3
Lhcx1WTno fit29.2no fit4.9no fit1.3
KO8no fit27.7no fit10.6no fit0.5
Co.48no fit29.5no fit6.4no fit2.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Madhuri, S.; Lepetit, B.; Fürst, A.H.; Kroth, P.G. A Knockout of the Photoreceptor PtAureo1a Results in Altered Diel Expression of Diatom Clock Components. Plants 2024, 13, 1465. https://doi.org/10.3390/plants13111465

AMA Style

Madhuri S, Lepetit B, Fürst AH, Kroth PG. A Knockout of the Photoreceptor PtAureo1a Results in Altered Diel Expression of Diatom Clock Components. Plants. 2024; 13(11):1465. https://doi.org/10.3390/plants13111465

Chicago/Turabian Style

Madhuri, Shvaita, Bernard Lepetit, Alexander Helmut Fürst, and Peter G. Kroth. 2024. "A Knockout of the Photoreceptor PtAureo1a Results in Altered Diel Expression of Diatom Clock Components" Plants 13, no. 11: 1465. https://doi.org/10.3390/plants13111465

APA Style

Madhuri, S., Lepetit, B., Fürst, A. H., & Kroth, P. G. (2024). A Knockout of the Photoreceptor PtAureo1a Results in Altered Diel Expression of Diatom Clock Components. Plants, 13(11), 1465. https://doi.org/10.3390/plants13111465

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop