| Home | E-Submission | Sitemap | Contact Us |  
top_img
Algae > Volume 39(4); 2024 > Article
Kim, Park, Nam, Lee, and Jin: Differential expression of silicon transporter genes in frustule formation of the marine diatom Thalassiosira eccentrica LIMS-PS-3165

ABSTRACT

Diatoms are photosynthetic unicellular organisms surrounded by silica-based cell walls known as frustules. Given their widespread distribution in aquatic systems and their prolific nature, diatoms play a pivotal role in global silica cycling. Silicon transporter (SIT) proteins play a critical role in this process by facilitating the import and intracellular transport of dissolved silica. Therefore, it is important to understand the mechanisms controlling the uptake and utilization of silica below the concentration range of Si, where free diffusion of Si is insufficient for proliferation. In our investigation, we observed varied frustule synthesis responses of Thalassiosira eccentrica LIMS-PS-3165 under different nutrient limitations through fluorescence labeling of newly synthesized frustules. Additionally, we examined the expression patterns of SIT genes following silicon replenishment and analyzed the sequences of two SITs, TeSIT1/2 and TeSIT3, from T. eccentrica. Upon silicon resupply after a period of silicon deprivation, TeSIT1/2 transcription was downregulated, whereas TeSIT3 expression remained unaffected. Analysis of the predicted amino acid sequences revealed that TeSIT1/2 and TeSIT3 share 10 transmembrane domains and two conserved GXQ motifs, indicating structural similarities and differences in their silicon transport mechanisms. These findings enhance our understanding of silicon transport responses in T. eccentrica under nutrient-limited conditions and contribute to broader knowledge of silica biogeochemistry in diatoms.

INTRODUCTION

As photosynthesizing unicellular microorganisms, diatoms play an important role as primary producers in marine ecosystems and they are responsible for 20–25% of carbon recycling (Smetacek 1999, Scala and Bowler 2001). Another element that is recycled by diatoms on a global scale, affecting biogeochemistry, is silicon (Si), which is estimated to be recycled at 240 ± 40 Tmol per year (Nelson et al. 1995). Silicon constitutes the frustule which is silicified cell walls of diatoms. The intricately patterned nanostructures of this tiny shell have inspired many studies on the biomineralization process of Si, as well as applications in nanotechnology, such as bio-manufacturing, fabrication of nanostructures (Parkinson and Gordon 1999, Huang and Daboussi 2017, Roychoudhury et al. 2022), and biomedical applications (Sardo et al. 2021).
Observations using electron microscopy have greatly contributed to our understanding of the morphology and structure of frustules, and this knowledge has been applied to the classification of diatoms. Additionally, monitoring the exocytosis of the silica frustules has provided deeper insights into biosilicification and frustule formation. Fluorescent probes are highly effective tools for this investigation, because they can be used to study the incorporation of Si into newly synthesized frustules. LysoTracker HCK-123, introduced by Desclés et al. (2008), is an excellent tracer for newly synthesized biosilica shells. This probe offers several advantages over other fluorescence markers, including no need for ionophores, compatibility with common filter sets for visualization, an enhanced signal-to-noise ratio, and stability in acidic cellular compartments.
The absorption of Si into cells is a fundamental step in frustule formation. In aquatic environments, Si exists as dissolved silicates (DSi), such as silicic acid and Si(OH)4. Above certain concentrations, silicic acid can freely diffuse into the cells (Raven 1983, Thamatrakoln and Hildebrand 2008). However, the Si concentration in the ocean is less than 30 μM, and therefore, active silicon transporters are required (Tréguer et al. 1995, Thamatrakoln and Hildebrand 2008). In diatoms, silicon transporter proteins, known as SITs are reported to be responsible for silicon uptake. These integral membrane proteins are believed to function as symporters of Na+, demonstrating a relatively high affinity for Na+ ions (Knight et al. 2016).
In this study, we identified and analyzed SIT sequences and expression patterns of the centric diatom Thalassiosira eccentrica LIMS-PS-3165, isolated from the southern coastal region of South Korea. We also traced frustule synthesis via cell-cycle synchronization using a fluorescent probe under nutrient-limited conditions. The species in this study exhibited a unique shape and arrangement of areolae compared to other T. eccentrica reported earlier, which may be related to the frustule formation process. Although this study only examined some characteristics of the SIT proteins involved in Si absorption and transport, it provides important molecular biological insights into the local strain of T. eccentrica, a species that has been less studied. In particular, when compared with the well-studied model species T. pseudonana, this research could be valuable for future comparative studies.

MATERIALS AND METHODS

Algal strain and culture conditions

The T. eccentrica LIMS-PS-3165 strain was obtained from the Korea Institute of Ocean Science and Technology (KIOST) and was isolated from the coastal regions of Jeju, South Korea. Stock cells were batch cultured in L1 medium (Guillard and Hargraves 1993) at 20°C under a 12 : 12 light : dark regime at 55 μmol photons m−2 s−1 at 50% humidity. For pre-experimental culture, cells were harvested by centrifugation (1,500 rpm, 20°C, 5 min), and washed with L1 medium twice, and then inoculated into fresh medium at a density of 0.1 × 104 cells mL−1. When the cells reached to exponential growth phase (approximately 0.5 × 104 cells mL−1), cells were harvested by centrifugation, washed twice with fresh experimental medium, and inoculated into new media for experiments. The media used for the nutrient deprivation experiment lacked one of the three nutrients from the L1 medium: nitrogen (N), phosphorus (P), or silicon (Si). For the Si availability experiment, the washed cells were arrested by incubating in Si-lacking medium for 1 d, followed by transfer to either Si+ or Si− media.
Cell density was measured by counting under light microscopy using hemocytometer (Marienfeld, Gütersloh, Germany), eight times as technical replicates, and the counts were averaged after excluding the highest and lowest values. All growth experiments included a minimum of three biological replicates.

Spectrophotometric measurement of mono silicic acid via molybdate blue method

Monosilicic acid, Si(OH)4 was determined based on the molybdate blue method using a scaled-down protocol provided with Hach silica testing kit (Hach Gmbh, Düsseldorf, Germany). One milliliter of the cell culture medium was centrifuged at 3,200 rpm to remove the cells. The remaining media samples were diluted by a dilution factor of 5. Forty microliters of molybdate 3 reagent was added and vortexed. The samples were incubated for 4 min. Forty microliters of citric acid reagent was added, vortexed, and incubated for 1 min. One packet of amino acid F reagent was dissolved in 1 mL dH2O and 40 μL of amino acid F reagent was added and vortexed. The sample was measured using an Ultrospec 2100 Pro (Amersham, GE Healthcare, Piscataway, NJ, USA) at OD803 using distilled H2O as a blank. A reagent blank sample was also prepared and measured but was omitted from the calculations owing to its negligible values. The optical density was measured for distilled H2O spiked with silicic acid at the range of 0–10 mg L−1. A standard curve was fitted to the standard samples to estimate silicic acid content.

Fluorescence microscopy observation

For fluorescence labeling of newly synthesized frustules, LysoTracker Yellow HCK-123 (Thermo Fisher Scientific, Waltham, MA, USA) was added to the experimental culture. HCK-123 stock solution (1 mM) was prepared in dimethyl sulfoxide, and added to the culture to final 1 μM. A 1 μL aliquot of cells was taken at specified time points for fluorescence microscopy using an Eclipse Ni-U microscope (Nikon, Tokyo, Japan) with a PhotoFluor LM-75 illumination system (89 North, Williston, VT, USA). Images were captured using a DS-Fi3 camera (Nikon) using NIS-Elements version 4.60, build 1171 (Nikon), as described in Won et al. (2023). Images of live cells, chlorophyll, and newly synthesized frustules were taken for each dataset. Live cells and chlorophyll were quantified to confirm cell viability. A minimum of 100 cells were analyzed in each experiment. To observe chlorophyll, a TxRed filter set was used, and a fluorescein isothiocyanate (FITC) filter was used to observe newly formed frustules stained with HCK-123.

Bioinformatic analyses of SIT proteins

Two SIT homologs, TeSIT1/2 and TeSIT3, were identified in the de novo assembled transcripts of T. eccentrica LIMS-PS-3165 (unpublished data). The coding sequences of the two SITs were deposited in GenBank with ID numbers PP526761 and PP526762, respectively (https://www.ncbi.nlm.nih.gov/gene). The predicted amino acid sequences of SITs listed in Table 1 were aligned and analyzed using the ClustalOmega web service from EMBL-EBI (Madeira et al. 2022). A phylogenetic tree of the SIT proteins was constructed using MEGA11 v11.0.10 (Tamura et al. 2021). The distances were computed using the Poisson correction method (Zuckerkandl and Pauling 1965), and the tree was evaluated through 500 bootstrap replicates (Kim et al. 2022). The phylogenetic tree was rooted with the SIT protein from Stephanoeca diplocostata, serving as a non-diatom SIT outgroup.
The structure of each SIT protein was predicted using Colabfold, which combines MMseqs2 and Alphafold2 to offer faster prediction of protein structures or complexes (Mirdita et al. 2022), and then visualized using UCSF ChimeraX (Pettersen et al. 2021). DeepLoc-2.0 web server (Thumuluri et al. 2022), Predotar v.1, and HECTAR (Gschloessl et al. 2008) were used to predict subcellular localization. Four web servers (Small et al. 2004) were used. Transmembrane regions were predicted using TMHMM 2.0, a web server (Krogh et al. 2001).

Real-time reverse transcriptase-polymerase chain reaction

Total RNA was isolated using a Hybrid-R RNA purification kit (Geneall Biotechnology, Seoul, Korea) and used as a template for cDNA synthesis using a High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Waltham, MA, USA), following the manufacturer’s instructions. For the transcription abundance, real-time PCR was performed using a PCR premix containing TB Green (TB Green Premix Ex Taq II; Takara Bio, Kusatsu, Japan). The thermal cycling profile comprised of 30 s pre-denaturing step at 95°C followed by 45 cycles of denaturing for 5 s at 95°C and annealing / extension for 30 s at 60°C. TATA-box binding protein was used as an internal control to normalize the data for relative quantification (Kang et al. 2012). Thermal cycling and quantification were performed using a Thermal Cycler Dice Real-Time System TP800 (Takara Bio), and each experiment was performed in triplicate. The primers used in these experiments are listed in Table 2.

RESULTS AND DISCUSSION

Thalassiosira eccentrica LIMS-PS-3165, henceforth referred to as T. eccentrica, exhibited characteristics typical of centric diatoms, including circular valves with radial symmetry. Colony formation was not observed under normal or stress conditions. The cells generally displayed a golden-brown color (Fig. 1A). The diameter of the valve was approximately 10 μm, while the height of a single cell ranged from approximately 5 to 10 μm, suggesting various stages of cell growth. Scanning electron micrographic images of the valve revealed irregular patterning of the areolae (Fig. 1B). T. eccentrica was reported to be found globally except for Arctic and Antarctic areas, and also frequently found in Korean coastal area (Park et al. 2016). It was reported that this species has hexagonal areolae, distributed regularly over the valve face with one row of marginal spines (Park et al. 2009, 2016). However, the cells observed in this study had areolae in irregular shape and size which were randomly distributed over the valve face, and no spines was observed (Fig. 1B). In diatoms, cell maturity can generally influence the shape and arrangement of areolae, which tend to become more regular as cells mature (Sato 2010). However, based on repeated observations, these characteristics in this study do not appear to be due to cell maturity.

Nutrient deprivation and cell cycle arrest in diatoms observed in newly synthesized silica frustules

After inoculating cultures with HCK-123, we stained newly synthesized frustules and observed them using a fluorescence microscope equipped with a FITC filter set. The frustules were categorized based on three observable states of fluorescence staining, as illustrated in Fig. 2B. Middle staining indicates the initiation of the G2 + M phase, during which valve fragments begin to be synthesized at the center of the mother cell. Tip-stained cells were observed when the newly synthesized frustule fragments fully expanded across the cell, forming the hypotheca of the upper daughter cells and the epitheca of the lower daughter cells. This tip fluorescence indicated that the cells had completed one division cycle and were at the end of G1 phase. Fully stained cells suggested that the cells had undergone another cycle of division, as only one of the two valves was synthesized per division cycle (Fattorini and Maier 2021, Bulankova et al. 2022).
Cell cycle progression, as indicated by HCK-123 staining, varies under different nutrient limitations. Under Si-deprived conditions, frustule staining ceased at the tip stage. However, under phosphorus and nitrogen deprivation, full staining of frustules occurred on days 3 and 4, respectively, implying the resumption of the cell division cycle (Fig. 2, Supplementary Fig. S1).
At the 24-h mark, the most representative staining stage for all three treatments was middle-stained cells. However, at 72 and 96 h, tip-stained frustules were most common. Notably, in phosphorus-depleted P(−) media, a small portion of cells showing fully stained frustules was observed starting on day 3 (Fig. 2, Supplementary Fig. S1). When cells are supplied with Si after 1 d of arrest, around 90% of cells were labeled with fluorescence in one of the theca within 24 h, implying that 90% of cells had gone through one cycle of cell division (Fig. 3A). In Si-deplete medium, staining of one theca was observed in about 10% of cells within 24 h, while approximately 80% remained with middle staining. This suggests that, although the remaining Si in the medium was absorbed into the cells, cell division did not occur during the first 24 h of Si depletion. Nevertheless, as shown in Fig. 2B, the gradual increase in tip staining over 4 d in Si-deplete medium suggests that cell division is proceeding slowly using the absorbed Si.
Silicon (Si) is a critical element in diatoms and constitutes the silica cell wall. Its absence significantly affects cell division. Phosphorus (P) and nitrogen (N) are essential for cell growth. Nitrogen, a key component of proteins and nucleic acids, is vital for maintaining cell structure. Its deficiency not only reduces photosynthesis and chlorophyll synthesis, but also leads to the remodeling of carbon metabolism, affecting overall energy processes (Alipanah et al. 2015). Phosphorus is another major component of nucleic acids and phospholipids. Phosphorus limitation also affects photosynthesis, because the photosynthetic machinery relies on the phospholipid bilayer in the thylakoid membrane. Phosphate deprivation also affects N uptake and utilization, leading to growth retardation (Alipanah et al. 2018).
We found that deprivation of these three elements affected cell division. Especially, N and Si deprivation resulted in complete arrest of cell division (Supplementary Fig. S1), and fluorescent signal did not increase after 24 h (data not shown). However, cell density under P deprivation showed slight increase after the third day of deprivation. This modest increase is likely due to the recycling of these essential elements from reserves and metabolic remodeling in response to nutrient shortages (Kim et al. 2017). Unlike nitrogen and phosphorus, Si, which is available only as a soluble pool, cannot be accessed from storage, making it a more effective inhibitor of cell division (Parslow et al. 1984, Harrison et al. 1990).

Phylogeny and conserved domains of SIT proteins: predicting functional differences

In our transcriptome analysis of T. eccentrica (unpublished laboratory data), we identified two SIT-coding sequences. One sequence was similar to SIT1 and SIT2 found in other diatoms, whereas the other was closer to SIT3; therefore, these two SITs were designated as TeSIT1/2 and TeSIT3, respectively. Sequence similarity was analyzed, and the protein structure was predicted based on the deduced amino acid sequences of the two SITs.
A phylogenetic tree constructed from deduced amino acid sequences from various diatoms including that of T. eccentrica is shown in Fig. 4. Fourteen SIT sequences were chosen from databases based on previous reports (Thamatrakoln and Hildebrand 2008, Durkin et al. 2016), and compared with TeSIT1/2 and TeSIT3.
The SIT1 and SIT2 sequences from T. pseudonana clustered with TeSIT1/2, with 100% bootstrap support, confirming a close relationship. The separation between the SIT1/2 and SIT3 sequences within Thalassiosira was well-supported with a 96% bootstrap value. The split between centric diatom SITs and most pennate diatom SITs was supported by 52% of the bootstrap values, with the SIT1 sequences of Fragilariopsis cylindrus and Phaeodactylum tricornutum diverging earlier than the separation between pennate and centric diatom SITs. Notably, SIT1 and SIT2 grouped separately from SIT3, indicating distinct evolutionary paths.
Previous studies on SIT sequences across various diatom species have shown that SITs from the Thalassiosira genus are clearly grouped within clade E (Hendry et al. 2018). Our findings suggest that within Thalassiosira, there are two distinct SIT lineages with different evolutionary histories. Evolutionary research on SITs indicates that SIT3 likely arose from a genome duplication event of SIT1-2, following a unique evolutionary path, particularly within Thalassiosirales (Bryłka et al. 2023).
In the alignments of the deduced amino acid sequences of TeSIT1/2 with those of SIT1 and SIT2 of the closely related species T. pseudonana, the pairwise identities were 83% and 79%, respectively. A comparison between the SIT3 of both species revealed 74% amino acid identity. In contrast, TeSIT1/2 and TeSIT3 showed only 50% similarity.
Despite relatively low similarity, 10 conserved transmembrane domains were detected in both TeSITs (Fig. 5). The 10 transmembrane domains predicted based on the hydropathy analysis by Thamatrakoln et al. (2006) was identified. However, when we attempted prediction using TMHMM, one of the 10 TM domains TeSIT1/2 had a low score, leading to an uncertain 3D model of the protein (Supplementary Fig. S3). Knight et al. (2023) found eight authentic transmembrane helices using the AlphaFold model, and the other two helices were partially submerged parallel to the membrane plane. Despite the small uncertainty due to differences in the algorithms used in the prediction, the combination of bioinformatics data and low-resolution experimental results collectively indicated the presence of 10 transmembrane domains. There were slight variations between the two TeSITs, which may contribute to functional differences between them.
Four GXQ motifs were conserved in all the aligned sequences, along with the positions of the motifs relative to the transmembrane domains. The GXQ motifs are reported to be located in pairs at the ends of helices 2 and 3 and helices 7 and 8, making them pairs that facilitate translocation if Si crosses the membrane (Thamatrakoln 2006, Knight et al. 2016), and are also found at the same position in T. eccentrica. It is presumed that the C-terminus of SIT1 interacts with other proteins and is less conserved (Hildebrand et al. 1997, 1998, Thamatrakoln et al. 2006, Durkin et al. 2016), which is also less conserved in SIT3 (Fig. 5).

Functional differences presumed by the transcriptional pattern and cellular localization of TeSITs

Although no significant differences in functional motifs were found between the two TeSITs, the transcription patterns in response to silicon availability suggest that they play different roles in silicon-related metabolic processes.
The relative mRNA expression levels of two silicon transporter genes, TeSIT1/2 and TeSIT3, were measured in conjunction with cell growth and Si concentration in the culture media at various time points.
The cell density of the experimental cultures remained relatively stable up to 9 h in both Si(+) and Si(−) conditions (Fig. 6A & B). The cell densities did not increase in the Si(−) condition as the Si concentration in the culture media slightly decreased. Conversely, under Si(+) conditions, the cell density showed a sharp increase between 9 and 12 h, followed by a steady increase until 48 h, albeit at a slower rate. The Si concentration in the Si(+) cultures decreased sharply in the first 3 h after Si was replenished and then continued to decrease at a stable rate until 48 h.
The transcription patterns of the two TeSIT genes, TeSIT1/2 and TeSIT3, differed. The transcription of the TeSIT1/2 decreased during the first 12 h following Si resupply after 1 d of Si starvation (Supplementary Fig. S2A). This is in line with a previous study by Shrestha and Hildebrand (2015), which reported a decrease in the transcription of TpSIT1 and TpSIT2 during the first 9 h following Si resupply. However, the transcription of TeSIT1/2 gradually increased over the next 36 h. This is considered a response to the exhaustion of Si in the medium (Fig. 6A). In contrast, the transcription ratio of TeSIT3 remained relatively constant, and was not affected by the presence of Si (Fig. 6C, Supplementary Fig. S2B). This is in line with a previous study by Shrestha and Hildebrand (Shrestha and Hildebrand 2015), who reported differences between TpSIT1, TpSIT2, and TpSIT3. In their study, the sequences of TpSIT1 and TpSIT2 proteins showed higher similarity (88%) than those of TpSIT3 (46 and 74% with TpSIT1 and TpSIT2, respectively), and the transcription levels of TpSIT1 and TpSIT2 changed according to the cell cycle, in contrast to that of TpSIT3 which showed a relatively low and constant transcription level.
The transcription patterns in response to the absence or presence of Si also differed in previous studies. The transcription of TpSTI1 and TpSIT2 is upregulated under Si depletion and is downregulated upon Si replenishment (Shrestha and Hildebrand 2015, Hildebrand et al. 2018). In contrast, SIT3 expression is unaffected by Si levels (Hildebrand et al. 2018). Since silicic acid can easily diffuse into the cell, it is presumed that some SIT proteins do not actively contribute to the uptake of Si (Thamatrakoln et al. 2006), but rather sense internal Si levels and regulate intracellular Si metabolism. In addition, most diatoms possess multiple SIT genes, the expression levels of which change in response to Si availability; however, the abundance of transcripts and proteins does not always correspond, suggesting multiple levels of regulation (Thamatrakoln and Hildebrand 2008, Shrestha and Hildebrand 2015). In this study, although the differences were slight, the transcription patterns of TeSIT1/2 and TeSIT3 were nearly opposite (Fig. 6C, Supplementary Fig. S2), suggesting a difference in the functionality of the two proteins. This divergence may also reflect slight variations in the predicted protein structure and localization (Supplementary Fig. S3). Notably, TeSIT3 was predicted to localize to the endoplasmic reticulum (ER) (Supplementary Table S1). When predicted using HECTAR with the heterokont range, SIT3 is predicted to localize to the mitochondria and other locations. Although the predictions for subcellular localization vary depending on the methods, both DeepLoc and Predotar indicate a high probability of ER localization for SIT3. Previous studies have reported that ER localization predictions can vary depending on the analytical method used (Anand et al. 2018, Darabi and Seddigh 2018). Kumar et al. (2017) reported that only 50% of the proteins in an ER protein dataset had an ER signal, suggesting that these signals are not essential for their presence in the ER. Based on the predicted localization and transcription patterns, TeSIT3, unlike TeSIT1/2, is presumed to play an important role in concentrating silicon in silicon-deposited vesicles and in frustule formation, rather than in silicon uptake.
Further investigation is necessary to fully elucidate the functions of TeSIT1/2 and TeSIT3. However, based on their distinctive transcription patterns, one may be involved in the internal regulation of Si transport, whereas the other may regulate the uptake of Si. This is critical because the timing of Si uptake and incorporation must be closely coordinated, particularly because silica deposition in the cell wall is tightly linked to Si transport and influenced by intracellular components (Thamatrakoln 2006).

CONCLUSION

Thalassiosira eccentrica is prevalent worldwide and is especially common in the coastal areas of South Korea. However, it has been relatively understudied in biomolecular studies. Our experiments, which involved labeling newly synthesized frustules with a fluorescent probe under nutrient-deprived conditions, revealed that Si was the most potent element for arresting the cell cycle. Furthermore, the two TeSITs appeared to have different functions in Si uptake and utilization, as suggested by their sequence similarities and markedly different transcription levels. Employing genetic manipulation techniques to target these SITs will be the key to defining the specific functions of each TeSIT protein.

ACKNOWLEDGEMENTS

This research was supported by a grant of the Korea Health Technology R & D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (grant number: HP23C0011).

CONFLICTS OF INTEREST

The authors declare that they have no potential conflicts of interest.

SUPPLEMENTARY MATERIALS

Supplementary Table S1
Probability for localization of silicon transporter (SIT) proteins (https://www.e-algae.org)
algae-2024-39-11-22-Supplementary-Table-S1.pdf
Supplementary Fig. S1
Growth rate and change of silicon (Si) concentration in the culture of Thalassiosira eccentrica LIMS-PS-3165 at different nutrient limitations (https://www.e-algae.org)
algae-2024-39-11-22-Supplementary-Fig-S1.pdf
Supplementary Fig. S2
Relative transcription level of TeSIT1/2 and TeSIT3 (https://www.e-algae.org)
algae-2024-39-11-22-Supplementary-Fig-S2,3.pdf

Fig. 1
Microscopic images of Thalassiosira eccentrica. (A) Light microscopy image capturing a single T. eccentrica cell from various angles. (B) Scanning electron microscopy image of T. eccentrica on the valve face. Scale bars represent: A, 5 μm; B, 10 μm. [Colour figure can be viewed at http://www.e-algae.org].
algae-2024-39-11-22f1.jpg
Fig. 2
Quantification of the proportion of cells exhibiting newly synthesized frustules when cultured under various nutrient limitations for 4 d. (A) Newly formed parts of the frustules, observed under a fluorescence microscope following staining with HCK-123. (B) The newly synthesized frustules were expressed as a percentage of fluorescence-labeled cells by counting 100 cells per each replication, and the results were averaged from three replicates. The white bars represent percentage of cells emitting green fluorescence in the middle, grey bars represent fluorescence in one of the two thecae (tip), and black bars represent fluorescence in the entire frustule (full). Scale bars represent: A, 5 μm. [Colour figure can be viewed at http://www.e-algae.org].
algae-2024-39-11-22f2.jpg
Fig. 3
Staining progression of newly synthesized frustules in cells arrested for 1 d, followed by treatment with silicon (Si) addition (A) or no Si addition (B). Cells were quantified according to the categories depicted in Fig. 2, allowing for a comparative analysis of frustule development in the presence and absence of Si. The white bars represent percentage of cells emitting green fluorescence in the middle, grey bars in one of the two valves. Cell counting and percentage calculations were conducted as described in Fig. 2.
algae-2024-39-11-22f3.jpg
Fig. 4
The phylogenetic relationships of silicon transporter (SIT) proteins inferred using the Neighbor-Joining method (Saitou and Nei 1987) based on amino acid sequences listed in Table 1. The two SIT proteins of Thalassiosira eccentrica in this study are marked with asterisk (*). Bootstrap values of 500 replicates are shown below the branches. A and E refer to clade A and clade E, respectively (adapted from Durkin et al. 2016, CC BY-NC-ND 4.0). The designation of the SIT proteins follows the names provided in the database where these sequences are listed.
algae-2024-39-11-22f4.jpg
Fig. 5
Alignment of degenerate amino acid sequences of silicon transporter (SIT) proteins. (A) Multiple alignments of SIT1 and 2 from Thalassiosira pseudonana and SIT1/2 of T. eccentrica. (B) Pairwise alignment of SIT3 from T. pseudonana and T. eccentrica. In both alignments, 10 transmembrane domains are highlighted with shading, and conserved GXQ domains are marked in red. [Colour figure can be viewed at http://www.e-algae.org].
algae-2024-39-11-22f5.jpg
Fig. 6
Growth and transcriptional responses to silicon (Si) availability. Growth of Thalassiosira eccentrica and Si concentration remaining in the growth media under Si-replete (A) and Si-depleted conditions (B). The grey areas represent the remaining Si concentration in the growth medium, expressed as a percentage, with 100% corresponding to the Si content in the fresh culture medium. The ratios of the transcript abundance of TeSIT1/2 and TeSIT3 under Si-replete and Si-depleted conditions (C). The white bars in (A) and (B) represent Si concentration (%) in the culture medium.
algae-2024-39-11-22f6.jpg
Table 1
SIT1, SIT2, and SIT3 genes with confirmed full-length sequences
Type Species ID SIT#
Pennate diatom Cylindrotheca fusiformis NCBI: AF052454 SIT1
NCBI: AF052455 SIT2
NCBI: AF052456 SIT3
Fragilariopsis cylindrus JGI: 157255 SIT1
JGI: 212017 SIT3
Nitzschia alba NCBI: DQ256058 SIT1
Phaeodactylum tricornutum JGI: 23423 SIT1
Centric diatom Skeletonema costatum NCBI: DQ256065 SIT2
Thalassiosira pseudonana JGI: 268895 SIT1
JGI: 41392 SIT2
JGI: 35133 SIT3
T. eccentrica LIMS-PS-3165 This study SIT1/2
This study SIT3
Non-diatom Stephanoeca diplocostata NCBI: HE981737 SITc
NCBI: HE981736 SITb
NCBI: HE981735 SITa

SIT, silicon transporter.

Table 2
Sequences of primers used in real-time PCR
Protein Gene Primer sequence 5′-3′ Amplicon size (bp)
ATA-box binding protein TBP FW CGCTGTGATTATGCGTCTTCGC 112
RV TTGTGCGTGCTTTTCGTTCCAG
SIT1/2 C53894_g1_i1_SIT1/2 FW CCCTTCACCCTCTTTCGCTACG 129
RV GATCAGGGCGACAATCGGACTC
SIT3 C49972_g1_i1_SIT3 FW GATTGGATCGGCATTTGCAGCG 101
RV TTCTCGGGTGTGCCAACGTAAG

SIT, silicon transporter.

REFERENCES

Alipanah, L., Rohloff, J., Winge, P., Bones, A. M. & Brembu, T. 2015. Whole-cell response to nitrogen deprivation in the diatom Phaeodactylum tricornutum. J. Exp. Bot. 66:6281–6296. doi.org/10.1093/jxb/erv340
crossref pmid pmc
Alipanah, L., Winge, P., Rohloff, J., Najafi, J., Brembu, T. & Bones, A. M. 2018. Molecular adaptations to phosphorus deprivation and comparison with nitrogen deprivation responses in the diatom Phaeodactylum tricornutum. PLoS ONE. 13:e0193335. doi.org/10.1371/journal.pone.0193335
crossref pmid pdf
Anand, G., Yadav, S., Tanveer, A., et al. 2018. Genome-wide assessment of polygalacturonases-like (PGL) genes of Medicago truncatula, Sorghum bicolor, Vitis vinifera and Oryza sativa using comparative genomics approach. Interdiscip. Sci. 10:704–721. doi.org/10.1007/s12539-017-0230-y
pmid pmc
Bryłka, K., Pinseel, E., Roberts, W. R., Ruck, E. C., Conley, D. J. & Alverson, A. J. 2023. Gene duplication, shifting selection, and dosage balance of silicon transporter proteins in marine and freshwater diatoms. Genome Biol. Evol. 15:evad212. doi.org/10.1093/gbe/evad212
crossref
Bulankova, P., Bilcke, G., Vyverman, W. & De Veylder, L. 2022. Cellular hallmarks and regulation of the diatom cell cycle. In : Falciatore A., Mock T., editors The Molecular Life of Diatoms. Springer, Cham, 229–263.
crossref
Darabi, M. & Seddigh, S. 2018. Structural, functional, and phylogenetic characterization of phosphoenolpyruvate carboxylase (PEPC) in C4 and CAM plants. Caryologia. 71:272–288. doi.org/10.1080/00087114.2018.1465762
crossref pmid pdf
Desclés, J., Vartanian, M., El Harrak, A., et al. 2008. New tools for labeling silica in living diatoms. New Phytol. 177:822–829. doi.org/10.1111/j.1469-8137.2007.02303.x
crossref pmid pmc pdf
Durkin, C. A., Koester, J. A., Bender, S. J. & Armbrust, E. V. 2016. The evolution of silicon transporters in diatoms. J. Phycol. 52:716–731. doi.org/10.1111/jpy.12441
crossref pdf
Fattorini, N. & Maier, U. G. 2021. Targeting of proteins to the cell wall of the diatom Thalassiosira pseudonana. Discover Materials. 1:5. doi.org/10.1007/s43939-021-00005-z
crossref pmid pmc pdf
Gschloessl, B., Guermeur, Y. & Cock, J. M. 2008. HECTAR: a method to predict subcellular targeting in heterokonts. BMC Bioinformatics. 9:393. doi.org/10.1186/1471-2105-9-393
crossref
Guillard, R. R. L. & Hargraves, P. E. 1993. Stichochrysis immobilis is a diatom, not a chrysophyte. Phycologia. 32:234–236. doi.org/10.2216/i0031-8884-32-3-234.1
crossref pdf
Harrison, P. J., Thompson, P. A. & Calderwood, G. S. 1990. Effects of nutrient and light limitation on the biochemical composition of phytoplankton. J. Appl. Phycol. 2:45–56. doi.org/10.1007/BF02179768
crossref
Hendry, K. R., Marron, A. O., Vincent, F., et al. 2018. Competition between silicifiers and non-silicifiers in the past and present ocean and its evolutionary impacts. Front. Mar. Sci. 5:22. doi.org/10.3389/fmars.2018.00022
crossref pmid pdf
Hildebrand, M., Dahlin, K. & Volcani, B. E. 1998. Characterization of a silicon transporter gene family in Cylindrotheca fusiformis: sequences, expression analysis, and identification of homologs in other diatoms. Mol. Gen. Genet. 260:480–486. doi.org/10.1007/s004380050920
crossref pmid pdf
Hildebrand, M., Lerch, S. J. L. & Shrestha, R. P. 2018. Understanding diatom cell wall silicification: moving forward. Front. Mar. Sci. 5:125. doi.org/10.3389/fmars.2018.00125
crossref
Hildebrand, M., Volcani, B. E., Gassmann, W. & Schroeder, J. I. 1997. A gene family of silicon transporters. Nature. 385:688–689. doi.org/10.1038/385688b0
crossref pdf
Huang, W. & Daboussi, F. 2017. Genetic and metabolic engineering in diatoms. Philos. Trans. R. Soc. B. 372:20160411. doi.org/10.1098/rstb.2016.0411
crossref pmid pmc pdf
Kang, L.-K., Tsui, F.-H. & Chang, J. 2012. Quantification of diatom gene expression in the sea by selecting uniformly transcribed mRNA as the basis for normalization. Appl. Environ. Microbiol. 78:6051–6058. doi.org/10.1128/AEM.00935-12
crossref pmid pmc
Kim, J., Brown, C. M., Kim, M. K., et al. 2017. Effect of cell cycle arrest on intermediate metabolism in the marine diatom Phaeodactylum tricornutum. Proc. Natl. Acad. Sci. U. S. A. 114:E8007–E8016. doi.org/10.1073/pnas.1711642114
crossref pdf
Kim, M., Oh, H. J., Nguyen, K. & Jin, E. 2022. Identification and characterization of Dunaliella salina OH214 strain newly isolated from a saltpan in Korea. Algae. 37:317–329. doi.org/10.4490/algae.2022.37.9.13
crossref pmid
Knight, M. J., Hardy, B. J., Wheeler, G. L. & Curnow, P. 2023. Computational modelling of diatom silicic acid transporters predicts a conserved fold with implications for their function and evolution. Biochim. Biophys. Acta Biomembr. 1865:184056. doi.org/10.1016/j.bbamem.2022.184056
crossref pmid pmc pdf
Knight, M. J., Senior, L., Nancolas, B., Ratcliffe, S. & Curnow, P. 2016. Direct evidence of the molecular basis for biological silicon transport. Nat. Commun. 7:11926. doi.org/10.1038/ncomms11926
pmid
Krogh, A., Larsson, B., von Heijne, B. & Sonnhammer, E. L. L. 2001. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J. Mol. Biol. 305:567–580. doi.org/10.1006/jmbi.2000.4315
crossref pmid pmc pdf
Kumar, R., Kumari, B. & Kumar, M. 2017. Prediction of endoplasmic reticulum resident proteins using fragmented amino acid composition and support vector machine. PeerJ. 5:e3561. doi.org/10.7717/peerj.3561
crossref pmid pmc pdf
Madeira, F., Pearce, M., Tivey, A. R. N., et al. 2022. Search and sequence analysis tools services from EMBL-EBI in 2022. Nucleic Acids Res. 50:W276–W279. doi.org/10.1093/nar/gkac240
crossref pmid pmc pdf
Mirdita, M., Schütze, K., Moriwaki, Y., Heo, L., Ovchinnikov, S. & Steinegger, M. 2022. ColabFold: making protein folding accessible to all. Nat. Methods. 19:679–682. doi.org/10.1038/s41592-022-01488-1

Nelson, D. M., Treguer, P., Brezezinski, M. A., Leynaert, A. & Quéguiner, B. 1995. Production and dissolution of biogenic silica in the ocean: revised global estimates, compairson with regional data and relationship to biogenic sedimentation. Global Biogeochem. Cycles. 9:359–372. doi.org/10.1029/95GB01070
crossref
Park, J. S., Jung, S. W. & Lee, J. H. 2009. A study on the fine structure of the marine diatoms of Korean Coastal waters - genus Thalassiosira 4. Algae. 24:67–77. doi.org/10.4490/algae.2009.24.2.067
crossref
Park, J. S., Jung, S. W., Lee, S. D., Yun, S. M. & Lee, J. H. 2016. Species diversity of the genus Thalassiosira (Thalassiosirales, Bacillariophyta) in South Korea and its biogeographical distribution in the world. Phycologia. 55:403–423. doi.org/10.2216/15-66.1
crossref pmid
Parkinson, J. & Gordon, R. 1999. Beyond micromachining: the potential of diatoms. Trends Biotechnol. 17:190–196. doi.org/10.1016/S0167-7799(99)01321-9
crossref pdf
Parslow, J. S., Harrison, P. J. & Thompson, P. A. 1984. Saturated uptake kinetics: transient response of the marine diatom Thalassiosira pseudonana to ammonium, nitrate, silicate or phosphate starvation. Mar. Biol. 83:51–59. doi.org/10.1007/BF00393085
pmid
Pettersen, E. F., Goddard, T. D., Huang, C. C., et al. 2021. UCSF ChimeraX: structure visualization for researchers, educators, and developers. Protein Sci. 30:70–82. doi.org/10.1002/pro.3943
crossref
Raven, J. A. 1983. The transport and function of silicon in plants. Biol. Rev. 58:179–207. doi.org/10.1111/j.1469-185X.1983.tb00385.x
crossref pmid pmc
Roychoudhury, P., Bose, R., Dąbek, P. & Witkowski, A. 2022. Photonic nano-/microstructured diatom based biosilica in metal modification and removal: a review. Materials. 15:6597. doi.org/10.3390/ma15196597
pmid
Saitou, N. & Nei, M. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406–425. doi.org/10.1093/oxfordjournals.molbev.a040454
crossref
Sardo, A., Orefice, I., Balzano, S., Barra, L. & Romano, G. 2021. Mini-review: potential of diatom-derived silica for biomedical applications. Appl. Sci. 11:4533. doi.org/10.3390/app11104533
crossref
Sato, S. 2010. Valve and girdle band morphogenesis in a bipolar centric diatom Plagiogrammopsis vanheurckii (Cymatosiraceae, Bacillariophyta). Eur. J. Phycol. 45:167–176. doi.org/10.1080/09670260903568483
crossref pmid pmc pdf
Scala, S. & Bowler, C. 2001. Molecular insights into the novel aspects of diatom biology. Cell. Mol. Life Sci. 58:1666–1673. doi.org/10.1007/PL00000804
crossref pmid pdf
Shrestha, R. P. & Hildebrand, M. 2015. Evidence for a regulatory role of diatom silicon transporters in cellular silicon responses. Eukaryot. Cell. 14:29–40. doi.org/10.1128/EC.00209-14
crossref pmid
Small, I., Peeters, N., Legeai, F. & Lurin, C. 2004. Predotar: a tool for rapidly screening proteomes for N-terminal targeting sequences. Proteomics. 4:1581–1590. doi.org/10.1002/pmic.200300776
crossref pmid
Smetacek, V. 1999. Diatoms and the ocean carbon cycle. Protist. 150:25–32. doi.org/10.1016/S1434-4610(99)70006-4
crossref pmid pmc pdf
Tamura, K., Stecher, G. & Kumar, S. 2021. MEGA11: molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 38:3022–3027. doi.org/10.1093/molbev/msab120

Thamatrakoln, K. 2006. Molecular insights into the function and regulation of diatom silicon transporters. Ph.D. dissertation. University of Calofornia, San Diego, CA, USA, 11 pp.

Thamatrakoln, K., Alverson, A. J. & Hildebrand, M. 2006. Comparative sequence analysis of diatom silicon transporters: toward a mechanistic model of silicon transport. J. Phycol. 42:822–834. doi.org/10.1111/j.1529-8817.2006.00233.x
crossref pmid pmc
Thamatrakoln, K. & Hildebrand, M. 2008. Silicon uptake in diatoms revisited: a model for saturable and nonsaturable uptake kinetics and the role of silicon transporters. Plant Physiol. 146:1397–1407. doi.org/10.1104/pp.107.107094
crossref pmid pmc pdf
Thumuluri, V., Almagro Armenteros, J. J., Johansen, A. R., Nielsen, H. & Winther, O. 2022. DeepLoc 2.0: multi-label subcellular localization prediction using protein language models. Nucleic Acids Res. 50:W228–W234. doi.org/10.1093/nar/gkac278
crossref pmid
Tréguer, P., Nelson, D. M., Van Bennekom, A. J., Demaster, D. J., Leynaert, A. & Quéguiner, B. 1995. The silica balance in the world ocean: a reestimate. Science. 268:375–379. doi.org/10.1126/science.268.5209.375
crossref pmid
Won, H., Ro, E., Seo, S., Kim, B.-H. & Jin, E. 2023. Isolation and cultivation of freshwater diatom Nitzschia palea HY1 for increasing biomass and fucoxanthin production. Algae. 38:191–202. doi.org/10.4490/algae.2023.38.9.3
crossref pdf
Zuckerkandl, E. & Pauling, L. 1965. Evolutionary divergence and convergence in proteins. In : Bryson V., Vogel H. J., editors Evolving Genes and Proteins. Elsevier, New York, 97–166. doi.org/10.1016/B978-1-4832-2734-4.50017-6
crossref
TOOLS
PDF Links  PDF Links
PubReader  PubReader
ePub Link  ePub Link
Full text via DOI  Full text via DOI
Download Citation  Download Citation
Supplement Supplement1
Supplement Supplement2
Supplement Supplement3
Supplement Supplement4
  Print
Share:      
METRICS
0
Crossref
0
Scopus
428
View
63
Download
Related article
Editorial Office
[14348] A-1716, Gwangmyeong Trade Center, 72 Iljik-ro Gwangmyeong-si. Gyeonggi-do, Korea
Tel: +82-2-899-5980  Fax: +82-2-899-5922    E-mail: editalgae@gmail.com
About |  Browse Articles |  Current Issue |  For Authors and Reviewers
Copyright © The Korean Society of Phycology.                 Developed in M2PI