Next Article in Journal
FOXM1, MEK, and CDK4/6: New Targets for Malignant Peripheral Nerve Sheath Tumor Therapy
Next Article in Special Issue
Ion fluxes Involved in the Adaptation of the Estuarine Diatom Coscinodiscus centralis Ehrenberg to Salinity Stress
Previous Article in Journal
β-Catenin and SOX2 Interaction Regulate Visual Experience-Dependent Cell Homeostasis in the Developing Xenopus Thalamus
Previous Article in Special Issue
Analysis of Predicted Amino Acid Sequences of Diatom Microtubule Center Components
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phenotypic Analysis and Molecular Characterization of Enlarged Cell Size Mutant in Nannochloropsis oceanica

1
State Key Laboratory of Marine Environmental Science, Xiamen University, Xiamen 361000, China
2
College of Ocean and Earth Sciences, Xiamen University, Xiamen 361000, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(17), 13595; https://doi.org/10.3390/ijms241713595
Submission received: 3 August 2023 / Revised: 22 August 2023 / Accepted: 25 August 2023 / Published: 2 September 2023
(This article belongs to the Special Issue Advances in Research of Algae, Cyanobacteria, and Phytoplankton)

Abstract

:
The cell cycle is the fundamental cellular process of eukaryotes. Although cell-cycle-related genes have been identified in microalgae, their cell cycle progression differs from species to species. Cell enlargement in microalgae is an essential biological trait. At the same time, there are various causes of cell enlargement, such as environmental factors, especially gene mutations. In this study, we first determined the phenotypic and biochemical characteristics of a previously obtained enlarged-cell-size mutant of Nannochloropsis oceanica, which was designated ECS. Whole-genome sequencing analysis of the insertion sites of ECS indicated that the insertion fragment is integrated inside the 5′-UTR of U/P-type cyclin CYCU;1 and significantly decreases the gene expression of this cyclin. In addition, the transcriptome showed that CYCU;1 is a highly expressed cyclin. Furthermore, cell cycle analysis and RT-qPCR of cell-cycle-related genes showed that ECS maintains a high proportion of 4C cells and a low proportion of 1C cells, and the expression level of CYCU;1 in wild-type (WT) cells is significantly increased at the end of the light phase and the beginning of the dark phase. This means that CYCU;1 is involved in cell division in the dark phase. Our results explain the reason for the larger ECS size. Mutation of CYCU;1 leads to the failure of ECS to fully complete cell division in the dark phase, resulting in an enlargement of the cell size and a decrease in cell density, which is helpful to understand the function of CYCU;1 in the Nannochloropsis cell cycle.

1. Introduction

The cell cycle is a fundamental cellular process in eukaryotes that includes cell growth, DNA replication, and cell division [1]. The cell cycle consists of pre-DNA synthesis (G1), DNA synthesis (S), post-DNA synthesis (G2), and mitosis (M phase). There are two checkpoints (eukaryotes) between the G1/S and G2/M phases: commitment (Chlamydomonas reinhardtii; S/M phase) and start (yeast; G1–S phase). These checkpoints are present to ensure cell size and DNA replication readiness [2,3,4]. The cell cycle is vital for microalgae. When microalgae are in a suitable environment, cells accumulate sufficient energy and the population rapidly increases in number after one cell cycle [5].
With the rapid development of omics and molecular biology, we can identify cell-cycle-related genes in different eukaryotes via homology and verify their functions. Specific cyclin and cyclin-dependent kinase (CDK) complexes are involved in specific cell cycle phases [4]. Phosphorylation-activated cyclin and CDK complexes can phosphorylate retinoblastoma (Rb) protein and release E2F/DP transcriptional activator to stimulate the expression of downstream S-phase-related genes, which promotes cell entry into the S phase [6]. Although the cell-cycle-related genes in Chlamydomonas reinhardtii are conserved [2,7], with more research, researchers have found new regulatory mechanisms. CDKA1 regulates the cell size commitment and may be required to set a higher-cell-size threshold controlling the transition between nondividing and dividing states. The MAT3/E2F/Dp1 pathway regulates the division number as well as the commitment size, but CDKA1 and the MAT3-E2F/DP1 pathway influence cell size and commitment in distinct ways [8]. CDKG1, a sizer protein that enables mother cells to execute the correct number of mitotic cell divisions according to mother cell size, binds to D-type cyclins and phosphorylates the retinoblastoma (RB) tumor suppressor pathway homolog MAT3 to regulate mitotic counting [9]. Atkins and Cross proposed a negative feedback loop model based on cyclins A-CDKA and B-CDKB, and CDKA activates cyclin B-CDKB; in turn, cyclin B-CDKB promotes mitotic entry and inactivates cyclin A-CDKA [10]. Similarly, new cell regulatory mechanisms have been found in Phaeodactylum tricornutum. CDKA2 was originally assigned to A-type CDKs, but displays some typical characteristics of plant-specific B-type CDKs and plays an essential role in the G2/M phase [11]. Diatom-specific cyclins respond to nutrient availability [12]. AUREOCHROME1a and the transcription factor bZIP10 induce dsCYC2 transcriptional expression under blue light, light-dependent dsCYC2 interacts with CDKA1 to control the G1/S cell size checkpoint, and dsCYC2 silencing decreases the rate of cell division in light–dark cycles [13].
The genus Nannochloropsis is assigned to the class Eustigmatophyceae, Heterokonta. Nannochloropsis species have a high triacylglycerol content, and their polar lipids are rich in eicosapentaenoic acid (EPA), an omega-3 long-chain polyunsaturated fatty acid. With the use of updated genetic engineering toolkits, Nannochloropsis has also been considered as a model for the study of microalgal lipid metabolism and as a chassis for synthetic biology [14]. Nannochloropsis loses its flagella and mobility because the genome has few flagellum-associated and movement-associated genes. Meanwhile, the absence of meiosis-specific and mating-associated genes leads to the loss of their ability to sexually reproduce [15]. Cell-cycle-related genes of Nannochloropsis have been reported [16]; however, the functions of these genes have not been described. The unique evolutionary characteristics of Nannochloropsis indicate that the cell cycle may differ from what is known. The phenomenon of cell enlargement is universal in algae research, and dry biomass, lipid, and fatty acid methyl ester (FAME) contents change in terms of cell size [17,18,19]. Meanwhile, the cell size is affected by environmental changes, including phosphorus stress [20,21], nitrogen stress [21], temperature [22], and especially light quality [22,23,24]. Blue light stimulates cell enlargement by delaying DNA replication and mitosis, while red light has the opposite effect [22,25]. In addition, gene mutations lead to cell enlargement, including UV mutagenesis [17], laboratory adaptation evolution [18], natural modification under laboratory conditions [19], and artificial selection [26]. Such phenotypic traits are caused by multigene mutations according to genome sequencing analysis [19]. Thus, the cell size regulators for Nannochloropsis are still unclear, and whether having a large cell size affects the photosynthetic characteristics, nutrient storage, and potential molecular mechanism of the large-cell-size phenotype has yet to be explored.
In this study, we used previously obtained ECS, which is a randomly inserted mutant whose genome contains the exogenous resistance gene (Ble). ECS was selected according to the FSC value (indicated cell size) via flow cytometry, and it has a large-cell-size phenotype. Thus, we first studied the phenotypic and biochemical characteristics of ECS. Furthermore, to analyze the phenomenon of cell enlargement, we integrated whole-genome sequencing analysis of insertion sites, an analysis of gene families and publicly available transcriptional data, flow cytometry analysis of the cell cycle, and real-time reverse transcription quantitative (RT-qPCR) analysis of cell-cycle-related genes. It was found that mutation of CYCU;1 caused ECS’s failure to fully complete cell division in the dark phase, resulting in a larger size and altering the expression level of cell-cycle-related genes. Therefore, our results explain the reason for ECS mutant cell enlargement and reveal that CYCU;1 regulates cell division, and thus controls cell size and cell density, which provides a target gene for genetic engineering to regulate microalgae cell size and density.

2. Results

2.1. Cell Density and Cell Size in ECS

Based on the growth curve, we observed that the growth of ECS was much slower than that of its wild type after culture to the middle logarithmic phase (Figure 1A). On day 8, the WT cell density was 4.3 × 107 cells mL−1, while that of ECS was 3.6 × 107 cells mL−1, and the latter was significantly lower than the former after day 8, decreasing by 17.4%, 14.6%, 14.2%, 13.1%, and 18.6%. In addition, we also noted changes in cell size, with ECS being significantly larger than WT (Figure 1B,C). ECS had a significantly increased cell size compared to WT at different time points based on an analysis of the forward scatter area (FSC-A) signal via flow cytometry, with an increased range of 19.3% to 45.9% (p < 0.01) and an average increase of approximately 29%. It was also found that ECS was larger than WT at each time point through fluorescence microscopy. The cell area of ECS increased significantly compared to WT, increasing from 5.9% to 30.6% (p < 0.01; Supplementary Table S2 and Figure S1), and the average increase was approximately 16.1%.

2.2. Chlorophyll Fluorescence and Chlorophyll Content in ECS

We found that the Fv/Fm was significantly higher for ECS than WT (Figure 2A), but Y (II) was significantly lower on day 6 (Figure 2B). The chlorophyll a content per milliliter was significantly higher for ECS than WT, and the cell level was significantly higher on day 6 (Figure 2C).

2.3. Dry Biomass, Protein, Total Carbohydrate, and Lipid Contents in ECS

In addition, we analyzed the biomass accumulation of ECS and WT and found that ECS increased by approximately 6.5% per unit volume compared with WT and the cell level increased by approximately 45% (p < 0.0001) (Figure 3A), and the protein and lipid contents were also higher, increasing by approximately 9.5% (p > 0.05) and 5.8% (p < 0.05), respectively (Figure 3B). It is worth noting that the protein, carbohydrate, and lipid contents per cell were significantly higher in ECS than in WT, especially the protein and lipid contents, which increased by 58.8% (p < 0.0001) and 53.5% (p < 0.0001), respectively (Figure 3C).

2.4. Insertion Is Identified to Be in Locus NO14G02130, Encoding a Cyclin Protein

Analysis of whole-genome sequencing of ECS showed that the insertion fragment was integrated inside the 5′-UTR of NO14G02130 at 643,397 of chromosome 14, which is −871 bp from the ATG start codon (Figure 4A,B). According to conserved domain prediction analysis using SMART and CDD prediction, NO14G02130 encodes a cyclin protein (Figure 4C) and has high amino acid sequence similarity with the Nannochloropsis gaditana cyclin-dependent protein EWM24594.1 (E-value = 3 × 10116, 52.58% identity on 74% sequence coverage using BLASTP). The fragment was successfully integrated into the genome (Figure 4D), and CYCU;1 was significantly downregulated (p < 0.0001, Figure 4E).

2.5. Identification of Cyclin Gene Family in Nannochloropsis oceanica

To identify the cyclin genes in N. oceanica IMET1 proteins, the N. oceanica IMET1 protein sequences were subjected to hmmsearch and BLASTP analysis using the Pfam database and the five species cyclin genes. A total of 19 cyclin genes were identified in N. oceanica IMET1 (Table 1), all including at least one cyclin domain (PF00134, PF02984, or PF08613). There are 2, 3, 2, 1, 4, 1, 4, and 2 members in the A-, B-, D-, H-, L-, T-, L-, P-, and U-type cyclins, respectively. Furthermore, we found that 6 genes contained both N- and C-domains, 10 genes contained N-domains, and 3 genes contained cyclin domains. The cyclin genes were located on 15 of 30 chromosomes; in particular, chromosomes 2, 3, and 12 had more than 2 cyclin genes, while the others had only 1 cyclin gene. The number of introns in cyclin genes ranged from 0 to 10. The length of cyclin proteins ranged from 310 to 843 amino acids (aa), and proteins ranged in size from 34 to 89 kDa. All analysis results of cyclin genes are presented in Table 1.

2.6. Motif Analysis and Light/Dark Transcriptomic Analysis of Cyclin Genes

The N. oceanica IMET1 phylogenetic tree was constructed using the maximum-likelihood approach (Figure 5A). The 19 cyclin genes were divided into three clusters: cluster 1 contained A-, B-, and D-type cyclins; cluster 2 contained P- and U-type cyclins; and cluster 3 mainly contained H-, L-, and T-type cyclins. Based on the results of the MEME algorithm, all cyclin genes contained at least one conserved cyclin motif (Figure 5B, Supplementary Table S3). The CYCU;1 gene had 1 motif 1, 4 motif 2, and 1 motif 4. To understand the expression of cyclins under light/dark conditions, we downloaded the light/dark transcriptome from NCBI. According to our analysis results, there was an oscillatory expression of cyclins; some cyclins were highly expressed in the light phase, such as CYCA;1 (ZT6), CYCB;1 (ZT12), and CYCB;3 (ZT9); and some were highly expressed in the dark phase, such as CYCL;1 (ZT18), and CYCU;2 (ZT21); however, CYCU;1 and CYCP;4 were highly expressed at different times of the day (Figure 5C, Supplementary Table S4), which suggests that these two cyclins may play significant roles in the cell cycle.
A phylogenetic study was performed using the maximum-likelihood tree approach to understand the evolutionary relationship of N. oceanica IMET1 cyclins with several Nannochloropsis species: N. oceanica CCMP1779 (No1779), Ochrophyta P. tricornutum (Phatr2), Chlorophyta C. reinhardtii (Chlre5), monocot O. sativa (Orysa), and eudicot A. thaliana (Arath). The tree was built using 19, 10, 24, 14, 34, and 36 cyclin genes from N. oceanica IMET1, N. oceanica CCMP1779, P. tricornutum, C. reinhardtii, O. sativa, and A. thaliana, respectively. The tree was divided into five clusters (Figure 6), including A/B-type, D-type, L/H/T/C-type, and U/P-type cyclins and cyclins. The U/P-type cyclin cluster contained O. sativa, P. tricornutum, and N. oceanica. On the one hand, CYCP5 and CYCP6 of P. tricornutum were considered to be markers of the G1 phase [27]. On the other hand, six P-type cyclins of P. tricornutum are thought to play a role in phosphate signaling because they are clustered with the PHO80-like proteins [12]; in addition, OsCYCP1;1, OsCYCP4;1, OsCYCP4;2, and OsCYCP4;3 of O. sativa respond to Pi starvation stress, and it was found that overexpression of Pi-starvation-induced OsCYCP4 interferes with the interaction of OsCYC and OsCDK [28,29], indicating that U/P-type cyclins of N. oceanica may be involved in phosphate signaling.

2.7. Failure of ECS to Fully Complete Cell Division in the Dark Phase

Flow cytometry was used to observe the signal intensity of DAPI and FSC-A; WT and ECS cell growth and DNA replication were observed in the light phase, and cell division was observed in the dark phase (Figure 7). At ZT0, the main peak of WT was located between the 1C and 2C peaks, and the main peak of ECS was located at the 2C peak. Partial cells of ECS completed the second round of DNA synthesis, and the 4C peak was obvious at ZT12, while the 4C peak of WT had not yet formed, and it was lower in height than that of ECS. ECS and WT still retained the 2C and 4C peaks when entering the dark phase, and with time, the 4C peaks of both ECS and WT gradually decreased. It is noteworthy that when entering the dark phase, the 2C peak of WT shifted to the 1C peak; however, ECS remained at the 2C peak and was higher than that of WT (Figure 7A). These results indicate that ECS may have failed to fully complete cell division in the dark phase, causing the main peak of ECS to be close to 2C at ZT0; thus, ECS could complete the second round of DNA synthesis to form the 4C peak earlier.
The cells were divided into different states based on DNA content. We observed that the proportion of 4C cells in ECS and WT gradually increased in the light phase, reaching the highest at ZT12, at 16.57% and 9.07%, respectively, and the proportion of 4C cells gradually decreased in the dark phase, and were 8.81% and 5.91%, respectively, at ZT21. The proportion of 1C cells in ECS and WT showed a gradually increasing trend, at 18.43% and 31.83% at ZT21, respectively. The proportion of S2-4C cells gradually decreased, which may have converted into 4C cells or undergone direct cell division. The proportion of S1-2C cells was relatively stable in WT, varying from 42.17% to 46.13%. S1-2C cells in ECS varied between 44.07% and 49.7%, except at ZT0 and ZT21, when they were 34.43% and 59.93%, respectively (Figure 7B). These results show that ECS may be unable to fully complete cell division in the dark phase (Figure S2), and as a result, it always maintains a high proportion of 4C cells and a low proportion of 1C cells.
The cell size changes of ECS and WT during light/dark cycles were analyzed using the mean FSC-A signal value. We found that the cell size of ECS and WT gradually increased in the light phase, reached a maximum at ZT12, and gradually decreased in the dark phase. The mean FSC-A of ECS at ZT15 was 40.05 × 104, while the mean FSC-A of WT was 34.34 × 104, and ECS increased by 16.62% compared to WT (Figure 7C). With an additional 24 h of culturing in the dark after sampling at the same time point, the mean FSC-A value of ECS was 35 × 104 at ZT15, while that of WT was 30.89 × 104 (Figure 7D).

2.8. RT-qPCR Analysis of Cell-Cycle-Related Genes

According to the light/dark cycle transcriptomic data analysis, we selected the low expression levels of CYCP;3 and CYCB;3; the medium expression levels of CYCD;1, CYCL;2, CYCH;1, and CYCP;2; and the high expression levels of CYCP;4, CYCU;1, and CYCP;3 as controls. We observed higher expression levels of CYCP;4 and CYCU;1 than the other cyclins at each time point (Figure 8A), consistent with the results of transcriptomic data analysis (Figure 5C). The expression levels of CYCU;1 and CYCP;4 were 4.2- to 153.7-fold and 11.1- to 128.2-fold higher, respectively, than that of CYCP;3.
Furthermore, we found that, compared to ZT0, the expression levels at the last time point of the light phase (ZT12) and the first two time points of the dark phase (ZT15 and ZT18) were significantly increased by 177%, 228%, and 285%, respectively. The expression levels at the first three time points of the light phase (ZT3, ZT6, and ZT9) were significantly decreased by 47.6%, 75.3%, and 69.1%, respectively, compared with ZT0 (Figure 8B).
To analyze the expression of cell-cycle-related genes of ECS, we drew a heatmap of the gene expression levels based on the average 2−ΔΔCT method. ECS shows an increased expression of cycle-related genes at time points ZT3 and ZT15 and a decreased expression at time points ZT18 and ZT21 compared with WT (Figure 8C). CYCP;3 at ZT0 and CYCP;4 and CDKC;1 at ZT3 were significantly increased, and CYCB;3, CYCH;1, and CDKA;1 at ZT21 were significantly decreased.
These results indicate that CYCU;1 and CYCP;4 were the most highly expressed cyclins at different times of the day. The expression levels of CYCU;1 at ZT12 in the light phase and ZT15 and ZT18 in the dark phase were significantly increased, indicating that CYCU;1 is involved in cell division in the dark cycle. Mutation of CYCU;1 induced CYCP;3 and CYCP;4 to significantly upregulate their expression to compensate for its missing function, and the expression level of CDKA;1 was significantly decreased, which affected cell cycle progression.

3. Discussion

3.1. Cell Size Affects Photosynthetic Characteristics and Nutrient Storage

Internal and external factors, such as genome, nutrients, and growth factors, determine cell size, which affects cellular signaling pathways and metabolic activity [30]. In our study, ECS, which was enlarged as a result of the mutation of CYCU;1, had increased chlorophyll a and Fv/Fm contents and decreased Y (II) after cell enlargement, which was probably caused by the package effect [31]. Thus, photosynthetic characteristics were correlated with cell size. Malerba et al. used Dunaliellla teriolecta as experimental material and obtained algal cells with different sizes through an artificial selection approach. They found that larger sizes led to reduced oxygen production per chlorophyll molecule; increased oxygen production, Fv/Fm, and photosynthetic pigment; and smaller light-harvesting antennae [26]. They also found that larger sizes led to an upregulation of CO2-concentrating mechanisms (CCMs), which improved the DIC uptake and led to faster growth and higher maximum biovolume density [32]. In experiments with different light qualities, it was also found that photosynthetic characteristics changed after cell enlargement, including an increase in chlorophyll content, especially the expression of photosynthesis-related genes [24]. In addition, larger cells achieved a balance between growth and photoprotection by sacrificing the growth rate when exposed to strong light [33]. Stephen et al. proposed that the downregulation of cell size maximizes light absorption under limited N to enable the large-scale production of algal oil in continuous output [34]. This means that cell size impacts photosynthetic characteristics.
Nutrient storage was correlated with cell size. In our study, the total biomass of ECS was slightly increased; however, the cell dry weight, protein content, and lipid content per cell were significantly increased (Figure 3C). Similarly, enlarged cell size in Parachlorella sp. through adaptive laboratory evolution under salt stress increased fatty acid (FA) and FAME contents and FA productivity and decreased biomass productivity [18]. The combined conditions of blue light and a temperature of 24 or 28 °C induced C. reinhardtii to have a large cell size, resulting in the highest protein content; on the contrary, the combined conditions of red-orange light and a temperature of 24 °C promoted carbohydrate content [22]. However, the two smaller-size mutants of Chlorella vulgaris were isolated through UV-C irradiation. They showed a significant decrease in biomass and a significant increase in cell concentration and lipid and triacylglycerol (TAG) contents [17]. This means that changes in cell size lead to changes in nutrient storage. Whether cell size can be regulated to obtain specific bioproducts is worth further research.

3.2. Cyclins Respond to Environmental Signals

Cyclins not only activate their CDK partners but also respond to environmental signals: Cyclins respond to light. dsCYC2, CYCP5, and CYCP6 respond to the dark–light transition in P. tricornutum [12,35] and the PHO80-like cyclin PC3933011 of Porphyridium cruentum [36]. Cyclins respond to silica availability, and dsCYC9 transcript levels were higher in cultures grown in the presence of silica than in those grown without silica [37]. Cyclins respond to phosphates. In our experiments, P/U-type cyclins of Nannochloropsis were divided into the same cluster as O. sativa and P. tricornutum based on the phylogenetic tree (Figure 6). Meanwhile, P/U-type cyclins are homologous with yeast PHO80 [38]; thus, these cyclins are thought to be involved in phosphate signaling [12,28,29]. OsCYCP1;1 could partially restore the phosphate signaling pathway in the yeast pho80 mutant [29], and SiPHO80 of Serendipita indica was upregulated in high-phosphate conditions and also restored Pi homeostasis of the yeast pho80 mutant [39]. We supposed that CYCU;1 was responsive to phosphate; hence, we also analyzed the mRNA levels of cyclins based on the published phosphate deprivation transcriptome and found that CYCU;1 and CYCP;4 still had the highest cyclin expression under phosphate stress, while the expression of CYCU;1 was significantly lower in PD48h than in PD0h (log2FC = −0.49, p < 0.01). Moreover, it was also reported that the expression level of U/P-type cyclins was stable mainly due to the rich inorganic phosphorus in f/2 medium [40]. Interestingly, by analyzing the mRNA levels of CYCU;1 under different experimental conditions, we found that CYCU;1 exhibited high cyclin expression (Supplementary Table S4), indicating that it is an important cyclin for Nannochloropsis. Whether CYCU;1 is involved in the phosphate signaling pathway and its interacting proteins remains to be validated.

3.3. Cyclins Play a Role in Cell Division

The cell cycle is a regulatory network composed of a series of cell-cycle-related regulatory genes, in which cyclins are involved in the whole process, especially cell division. In this study, although all cyclins were found to have at least one cyclin box, only CYCU;1 had the CD20558 structural domain through CDD prediction (Table 1), and it contained Y(L/A)(E/A)RI(F/A)(R/K)(Y/F) and (N/S)VHRLL(V/I)T motifs [29,41]. In Arabidopsis, CYCPs may be involved in cell division, cell differentiation, and the nutritional status of the cell by interacting with CDKA1 [42]. CYCP2;1 is transcriptionally activated by carbohydrate signals, and it interacts with three of the five mitotic CDKs to promote the G2-to-M transition of cell cycle progression [43]. In rice, overexpressed OsCYCP4 induced by phosphate induction could compete with the other cyclins for binding with CDKs, and suppress growth by reducing cell proliferation [28]. OsCYCP3;1 is specifically expressed in the root meristem epidermis and lateral root cap; it regulates meristem cell division by associating with and activating CKDB2;1 [41]. In addition, different types of cyclins have been reported to be involved in cell division, including A-type, B-type, and D-type cyclins.
In microalgae, Chlamydomonas CYCA regulates the timing of cell division [10], CYCB1 is required for spindle formation, and CYCB1 is synthesized before each division in the multiple fission cycle and then is rapidly degraded before division occurs. CYCB1/CDKB1 and APC modulate microtubule function and assembly while regulating mitotic progression [44]. CDKG1 binds to D-type cyclins and phosphorylates MAT3 to regulate mitotic counting [9]. CYCB1, CYCB2, dsCYC3, and dsCYC4 of P. tricornutum are expressed at the G2/M phase, and CYCB1 is a mitotic biomarker cyclin [12,27]. In plants, D-type cyclins are crucial for growth and development and regulate the cell division process through CDK-CYCD. Overexpression of the Populus D-type cyclin PsnCYCD1;1 gene in Arabidopsis can promote cell division and lead to small cell generation [45], and PsnCYCD1;1 overexpression in plants could accelerate cell division, causing the generation of small cells and severe morphological changes in the vascular bundles [46]. Overexpression of PtoCYCD3;3 increases the thickness of secondary xylem and phloem by increasing cambium cell activity, and it may interact with 12 PtoCDK proteins to regulate cell cycle programming [47]. Anantha et al. used overexpressing lines of PpCYCD1, PpCDKA2, PpCYCD2, and PpCDKA1, and the phenotypic data confirmed their controlled G1 to S and G2 to M transitions; it was also found that overexpression of PpCYCD1 or PpCDKA2 led to larger gametophytes [48]. Meanwhile, the Arabidopsis CYCD4;2 gene has a promotive function in cell division by binding and activating CDKA;1 even if it lacks the Rb-binding motif and the PEST sequence, and overexpressing plants showed faster callus formation in a medium containing lower concentrations of auxin [49].
In this study, we observed that the main peak of ECS was maintained at 2C throughout an additional 24 h of dark culture (Figure S2), which suggests that ECS cell division was affected. Furthermore, RT-qPCR demonstrated that the expression level of CYCU;1 significantly increased at the end of the light phase and the beginning of the dark phase (Figure 8B), indicating that the mutation of CYCU;1 caused incomplete cell division of ECS in the dark phase, resulting in lower cell density and larger cells. In addition, in our study, ECS downregulated the expression level of CDKA;1. Whether or not CDKA;1 partners with CYCU;1 to regulate cell division is worth further research.
Cell size plays an essential role in cell division, especially in microalgae with multiple divisions, and the threshold cell size determines the number of successive cell divisions [50,51]. The current understanding is that the controlling cell division model includes timers, sizers, and adders [52]. Our study shows that ECS is larger during cell division in the dark phase (Figure 7C,D), but whether this indicates a high threshold for cell division requires further verification. In addition, Nannochloropsis is a chassis organism of synthetic biology, and how to increase cell growth is an important research topic. CYCU;1 provides a target gene for genetic engineering to regulate microalgae cell size and cell density.

4. Materials and Methods

4.1. Microalgal Strain and Culture Conditions

The microalgal strain used in this study, N. oceanica IMET1, was a kind gift from Danxiang Han (Institute of Hydrobiology, Chinese Academy of Sciences). ECS, a randomly inserted mutant whose genome contains the exogenous resistance gene (Ble), was obtained through electroporation, and was selected according to the FSC value (indicating cell size) via flow cytometry. This means ECS possesses a large-cell-size phenotype. In addition, ECS is resistant to Zeocin, and the screening concentration was 1 μg mL−1.
Microalgae were inoculated into f/2 medium of artificial seawater (S9983, Sigma-Aldrich, St. Louis, MO, USA) supplemented with nutrient enrichment (G0154, Sigma-Aldrich, USA). The cells were cultivated in an artificial climate incubator provided with LED light of 50 μmol photons m−2 s−1 under a 12/12 h light/dark photoperiod at 25 °C.

4.2. Measurement of Growth and Cell Size

Measurement of cell density and cell size: The growth curve of microalgae was determined using a CytoFLEX S flow cytometer (Beckman, Brea, CA, USA) with an initial concentration of 106, and cells were collected at 2-day intervals. The cells were collected and analyzed through flow cytometry to measure cell size (channel FSC-A, forward scatter area). Cell density and mean FSC-A data were analyzed through two-way ANOVA (GraphPad Prism). The cells were measured using fluorescence microscopy (DM6 B, Leica, Wetzlar, Germany); 100 cell areas were calculated using Python scripts, and cell area data were analyzed using t-test (GraphPad Prism).

4.3. Determination of Chlorophyll Content and Chlorophyll Fluorescence

Determination of chlorophyll content: Chlorophyll a in fresh cells (3 mL culture) was extracted overnight with 3 mL of pure methanol at 4 °C. Before measurement, samples were centrifuged at 4000 rpm for 10 min, and the absorbance of the supernatant was measured using a scanning spectrophotometer (TU-1810, PERSEE, Beijing, China). The Chl a concentration was calculated using the following equation: Chl a (μg mL−1) = 16.29 × (A665 − A750) − 8.54 × (A652 − A750) [53].
Chlorophyll fluorescence measurement: The maximum quantum yield of PSII (Fv/Fm) and the effective quantum yield of PSII (Y(II)) were assessed using a pulse-amplitude-modulated (PAM) chlorophyll fluorometer (MULTI-COLOR-PAM, Walz, Germany) [54]. After dark acclimation for 15 min at the culture temperature, measurements were taken at a light intensity of 50 µmol photons m−2 s−1, which was similar to the growth light level. Chlorophyll content and chlorophyll fluorescence data were analyzed using two-way ANOVA (GraphPad Prism).

4.4. Dry Biomass Measurement and Analysis of Protein, Total Carbohydrate, and Lipid Content

Dry biomass measurement: Samples were harvested after 14 days of cultivation. The wet biomass was centrifuged at 4000 rpm at 4 °C for 10 min and then washed twice with Milli-Q water to remove the salt. The cell pellets were lyophilized with vacuum freeze-drying equipment (FreeZone-18, Labconco, Kansas City, MO, USA) for 24 h. After drying, the cell pellets were weighed and stored at −80 °C until proteins, carbohydrates, and lipids were extracted and analyzed.
For protein content, total protein was extracted from the lyophilized algal biomass and determined using the Bradford method (P0006C, Beyotime, Shanghai, China). For total carbohydrate content, total carbohydrates were determined by using a commercially available Total Carbohydrate Assay Kit (No. BC2715, Beijing Solarbio Science & Technology Co., Ltd., Beijing, China). For total lipid content, total lipids were extracted from the lyophilized algal biomass in chloroform–methanol (2:1, v/v). After centrifugation, the supernatant was collected into a new glass bottle, dried under a continuous stream of nitrogen gas, and then weighed. The above data were analyzed using t-test (GraphPad Prism).

4.5. Whole-Genome Sequencing and Analysis of Insertion Sites

Genomic DNA was extracted using the improved CTAB method, and total DNA was qualified and quantified using a NanoDrop and Qubit 2.0 fluorometer. Genomic DNA was sequenced using next-generation sequencing (NGS) by Novogene Co., Ltd. (Beijing, China); a library with an average 350 bp insertion size was constructed, and 2 × 150 bp paired-end sequencing was implemented on a NovaSeq 6000 system. Identification of insertion sites: Briefly, clean reads of filter sequencing data from ECS alignments were performed using the Burrows–Wheeler aligner (BWA) [55] using default parameters for paired-end reads and mapped to wild-type N. oceanica IMET1 genome version 2 (https://nandesyn.single-cell.cn/; accessed on 29 July 2022). The unmapped and mate unmapped reads were extracted and assembled into contigs using Spades [56]. Contigs with coverage depth below 20 were filtered (ref_unmapped), and plasmids containing Zeocin gene sequences were searched against ref_unmapped by BLAST [57]. The sites of the inserted fragment were determined using the E-value and the matching length. The insertion sites were manually inspected and visualized using Integrated Genomics Viewer (IGV) [58].

4.6. Confirmation of Insertion Sites

According to the insertion site results, specific primers were designed to detect insertion sites in the ECS genome. The primer sequences are presented in Supplementary Table S1. PCR was carried out by using 2× Accurate Taq Master Mix (dye plus) (Accurate Biotechnology Co., Ltd., Changsha, China) following the manufacturer’s instructions. A total of 50 µL of PCR amplification reaction volume contained 1 μL of DNA, 25 μL of 2× Tag mix, 22 μL of ddH2O, and 1 μL of 10 μM forward and reverse primers. The PCR program was initiated at 94 °C for 90 s, followed by 35 cycles of 98 °C for 10 s, 55 °C for 30 s, and 72 °C for 2 min.
Total RNA was extracted from the samples using the SteadyPure Universal RNA Extraction Kit (Accurate Biotechnology Co., Ltd., Changsha, China) according to the manufacturer’s instructions. Total RNA was qualified and quantified using a NanoDrop. cDNA was synthesized using the Evo M-MLV RT Kit with gDNA Clean for RT–qPCR (Accurate Biotechnology Co., Ltd., Changsha, China) following the manufacturer’s instructions. RT-qPCR was carried out by using an SYBR® Green Premix Pro Taq HS qPCR Kit (Accurate Biotechnology Co., Ltd., Changsha, China) following the manufacturer’s instructions. A total reaction volume of 20 μL contained 1 μL of cDNA, 10 μL of SYBR Green (ROX mixed), 8.2 μL of ddH2O, and 0.4 μL of 10 μM forward and reverse primers. The PCR program was initiated at 95 °C for 20 s, followed by 40 cycles of 95 °C for 5 s and 60 °C for 30 s. mRNA levels were quantified using the 2−ΔΔCT method [59]. Ubiquitin conjugating enzyme (UBCE) gene was used as an internal control. RT-qPCR data were analyzed using t-test (GraphPad Prism).

4.7. Identification and Characterization of Cyclin Genes in Nannochloropsis oceanica

Genomes, proteins, and gene structure annotation files of N. oceanica IMET1 were available at NanDeSyn (http://nandesyn.single-cell.cn/; accessed on 29 July 2022), and it was found to contain 10,333 genes. The cyclin protein sequence of N. oceanica IMET1 was identified through the following steps: First, hmmsearch [60] software was used with the hidden Markov model (HMM) for PF00134 (cyclin, N-terminal domain), PF02984 (cyclin, C-terminal domain), and PF08613 (cyclin, Pfam database, http://pfam.xfam.org/; accessed on 29 July 2022) to search the cyclin domain of the proteins of N. oceanica IMET1 proteomes, and the E-value threshold was set at 103. Second, published cyclin proteins from Arabidopsis thaliana, Oryza sativa, Nannochloropsis oceanica CCMP1779, Chlamydomonas reinhardtii, and Phaeodactylum tricornutum were downloaded and used as query sequences to search against the N. oceanica IMET1 proteomes, and the E-value threshold was set at 105. Third, the N. oceanica IMET1 annotation files were compared. Finally, the screening results of the three steps were merged, followed by manual removal of redundant or repetitive sequences. Prediction of the cyclin structure domain was carried out within the amino acid sequence of the candidate cyclin protein family member in N. oceanica IMET1 with the help of the Conserved Domain Database (CDD, http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi; accessed on 29 July 2022) of the National Center for Biotechnology Information (NCBI, http://www.ncbi.nlm.nih.gov/; accessed on 29 July 2022), the Simple Modular Architecture Research Tool software (SMART, http://smart.embl-heidelberg.de; accessed on 29 July 2022), and the database of protein domains, families, and functional sites (Prosite, https://prosite.expasy.org/; accessed on 29 July 2022). Candidate genes without the cyclin structure domain were removed. The physicochemical properties of the cyclin protein of N. oceanica IMET1, including hydropathicity, molecular mass, instability index, and so forth, were predicted via ProtParam software from the Expasy database (http://web.expasy.org/protparam/; accessed on 29 July 2022).

4.8. Motif Analysis, Multiple Alignment, and Phylogenetic Analysis

The conserved protein motifs of the N. oceanica IMET1 cyclin protein sequences were identified using Multiple Em for Motif Elicitation (MEME) [60]. The analysis was performed using the following parameters: any number of repetitions, a maximum of 10 motifs, and optimum motif width from 10 to 60 amino acid residues. The conserved motifs were annotated using InterProScan (http://www.ebi.ac.uk/interpro/; accessed on 29 July 2022) [61]. Multiple alignment of the N. oceanica IMET1cyclin was performed using MAFFT [62], and the phylogenetic tree was inferred under maximum likelihood (RaxML) [63] with 1000 bootstrap replicates. The tree was visualized using EvolView (https://www.evolgenius.info/evolview/; accessed on 29 July 2022) [64].
To comprehensively understand the evolutionary relationships of cyclin members, a phylogenetic tree was generated based on Chlamydomonas reinhardtii cyclin (Chlre5) proteins, Nannochloropsis oceanica CCMP1779 cyclin (No1779) proteins, and Phaeodactylum tricornutum cyclin (Phatr2) proteins, especially Oryza sativa cyclin (Orysa) and Arabidopsis thaliana cyclin (Arath) proteins, as references to categorize the N. oceanica IMET1 cyclin proteins.

4.9. Transcriptomic Data Analysis

To gain insight into the expression profile of N. oceanica IMET1 cyclin genes, transcriptomic data were analyzed during the light/dark cycles. We downloaded transcriptomic data of PRJNA285666 (light/dark cycles) from the NCBI database. The expression levels of N. oceanica IMET1 cyclin genes were quantified based on their fragments per kilobase of transcript per million mapped reads (FPKM). Briefly, using the N. oceanica IMET1 genome (IMET1v2) as the reference genome, the example transcriptomic data were aligned to the reference genome with TopHat [65], and gene expression was measured as the number of reads aligned to annotated genes through Cufflinks [66] and normalized to FPKM. The FPKM values were log2-transformed and used to construct a heatmap.

4.10. Cell Cycle Analysis Using Flow Cytometry

The cells were cultured to the middle of the logarithmic phase, and Zeitgeber time 0 (ZT0) was used as the starting point. Two samples were collected every 3 h, and one sample was moved to a dark environment and cultured for 24 h. The ZT0 samples were collected in the dark, and the ZT12 samples were collected in the light. The cells were centrifuged at 4000× g for 5 min at 4 °C, and the cell pellet was fixed overnight in 70% ethanol at 4 °C, washed twice in phosphate-buffered saline (PBS) at pH 7.5, and stained in PBS 0.1% Triton-X and 1 μg mL−1 of 4′,6′-diamidino-2-phenylindole (DAPI, MBD0015, Sigma-Aldrich, USA) for 20 min. Flow CytoFLEX S flow cytometry analysis of the cell cycle was performed using a 375 nm ultraviolet (UV) light laser and a 450/45 bandpass filter, and 20,000 cells were analyzed per sample. FlowJo (v10.8.1) was used for data analysis.

4.11. RT–qPCR Analysis of N. oceanica IMET1 Cell-Cycle-Related Genes

To measure the expression of cell-cycle-related genes in cultured cells across a 24 h day–night cycle, the cells were cultured to the middle of the logarithmic phase and, at ZT0 as the starting point, were harvested every 3 h by centrifugation, washed twice in Milli-Q water, flash-frozen in liquid nitrogen, and stored at −80 °C. For each time point, three biological repeats were performed. Total RNA was extracted, cDNA was synthesized, and gene expression analysis was performed using the RT-qPCR referenced in Materials and Methods, Section 2.4. The primer sequences are presented in Supplementary Table S1. All experiments were performed with three biological repeats, with two technical repeats generated for each one.

5. Conclusions

In the present study, we observed that ECS cell enlargement significantly decreased cell density, and increased Fv/Fm, chlorophyll a content, cell dry weight, and protein and lipid content per cell. Furthermore, we found and demonstrated that the insertion fragment was integrated inside the 5′-UTR of the U/P-type cyclin CYCU;1, which is a highly expressed cyclin. Furthermore, ECS could not fully complete cell division in the dark phase, and the expression level of CYCU;1 increased significantly at the end of the light phase and the beginning of the dark phase, indicating that this cyclin regulates cell division in the dark phase. Our results explain the reason for ECS cell enlargement and reveal that CYCU;1 controls cell size and cell density by regulating cell division, providing a target gene for genetic engineering to regulate microalgae cell density.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms241713595/s1.

Author Contributions

The manuscript was written by W.X.; physiological experiments were designed by Y.L. (Yihua Lin) and Y.W.; W.X. assisted with the analysis of mutant genome data and cyclin gene family; W.X. assisted with the RT-qPCR of mutant and WT under different light/dark conditions; W.X. and Y.L. (Yanyan Li) assisted with cell cycle experiments; H.Z. (Hongmei Zhu) assisted with downloading the transcriptome data; H.Z. (Hantao Zhou) assisted with the experimental design and manuscript editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the National Natural Science Foundation of China (No. 31371988).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The whole-genome-sequencing clean reads of ECS and WT were deposited in the GenBank SRA database under accession number PRJNA1001626.

Acknowledgments

We wish to thank all members of Zhou Lab for providing valuable suggestions and help during the research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Umen, J.G. Sizing up the cell cycle: Systems and quantitative approaches in Chlamydomonas. Curr. Opin. Plant Biol. 2018, 46, 96–103. [Google Scholar] [CrossRef] [PubMed]
  2. Cross, F.R.; Umen, J.G. The Chlamydomonas cell cycle. Plant J. 2015, 82, 370–392. [Google Scholar] [CrossRef] [PubMed]
  3. Huysman, M.J.; Vyverman, W.; De Veylder, L. Molecular regulation of the diatom cell cycle. J. Exp. Bot. 2014, 65, 2573–2584. [Google Scholar] [CrossRef] [PubMed]
  4. Harashima, H.; Dissmeyer, N.; Schnittger, A. Cell cycle control across the eukaryotic kingdom. Trends Cell Biol. 2013, 23, 345–356. [Google Scholar] [CrossRef] [PubMed]
  5. Zhao, P.Z.; Ouyang, L.L.; Shen, A.L.; Wang, Y.L. The cell cycle of phytoplankton: A review. J. World Aquac. Soc. 2022, 53, 799–815. [Google Scholar] [CrossRef]
  6. He, C.; Liang, J.; Wu, Z.; Zhuge, X.; Xu, N.; Yang, H. Study on the interaction preference between CYCD subclass and CDK family members at the poplar genome level. Sci. Rep. 2022, 12, 16805. [Google Scholar] [CrossRef]
  7. Bisova, K.; Krylov, D.M.; Umen, J.G. Genome-wide annotation and expression profiling of cell cycle regulatory genes in Chlamydomonas reinhardtii. Plant Physiol. 2005, 137, 475–491. [Google Scholar] [CrossRef]
  8. Cross, F.R. Regulation of Multiple Fission and Cell-Cycle-Dependent Gene Expression by CDKA1 and the Rb-E2F Pathway in Chlamydomonas. Curr. Biol. 2020, 30, 1855–1865.e4. [Google Scholar] [CrossRef]
  9. Li, Y.; Liu, D.; Lopez-Paz, C.; Olson, B.J.; Umen, J.G. A new class of cyclin dependent kinase in Chlamydomonas is required for coupling cell size to cell division. eLife 2016, 5, e10767. [Google Scholar] [CrossRef]
  10. Atkins, K.C.; Cross, F.R. Interregulation of CDKA/CDK1 and the Plant-Specific Cyclin-Dependent Kinase CDKB in Control of the Chlamydomonas Cell Cycle. Plant Cell 2018, 30, 429–446. [Google Scholar] [CrossRef]
  11. Huysman, M.J.; Tanaka, A.; Bowler, C.; Vyverman, W.; De Veylder, L. Functional characterization of the diatom cyclin-dependent kinase A2 as a mitotic regulator reveals plant-like properties in a non-green lineage. BMC Plant Biol. 2015, 15, 86. [Google Scholar] [CrossRef] [PubMed]
  12. Huysman, M.J.; Martens, C.; Vandepoele, K.; Gillard, J.; Rayko, E.; Heijde, M.; Bowler, C.; Inze, D.; Van de Peer, Y.; De Veylder, L.; et al. Genome-wide analysis of the diatom cell cycle unveils a novel type of cyclins involved in environmental signaling. Genome Biol. 2010, 11, R17. [Google Scholar] [CrossRef] [PubMed]
  13. 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]
  14. Poliner, E.; Farre, E.M.; Benning, C. Advanced genetic tools enable synthetic biology in the oleaginous microalgae Nannochloropsis sp. Plant Cell Rep. 2018, 37, 1383–1399. [Google Scholar] [CrossRef] [PubMed]
  15. Guo, L.; Liang, S.; Zhang, Z.; Liu, H.; Wang, S.; Pan, K.; Xu, J.; Ren, X.; Pei, S.; Yang, G. Genome assembly of Nannochloropsis oceanica provides evidence of host nucleus overthrow by the symbiont nucleus during speciation. Commun. Biol. 2019, 2, 249. [Google Scholar] [CrossRef] [PubMed]
  16. Poliner, E.; Panchy, N.; Newton, L.; Wu, G.; Lapinsky, A.; Bullard, B.; Zienkiewicz, A.; Benning, C.; Shiu, S.H.; Farre, E.M. Transcriptional coordination of physiological responses in Nannochloropsis oceanica CCMP1779 under light/dark cycles. Plant J. 2015, 83, 1097–1113. [Google Scholar] [CrossRef]
  17. Carino, J.D.; Vital, P.G. Characterization of isolated UV-C-irradiated mutants of microalga Chlorella vulgaris for future biofuel application. Environ. Dev. Sustain. 2023, 25, 1258–1275. [Google Scholar] [CrossRef]
  18. Kim, Z.H.; Kim, K.; Park, H.; Lee, C.S.; Nam, S.W.; Yim, K.J.; Jung, J.Y.; Hong, S.-J.; Lee, C.-G. Enhanced Fatty Acid Productivity by Parachlorella sp. a Freshwater Microalga, via Adaptive Laboratory Evolution Under Salt Stress. Biotechnol. Bioprocess Eng. 2021, 26, 223–231. [Google Scholar] [CrossRef]
  19. Shin, S.E.; Koh, H.G.; Kang, N.K.; Suh, W.I.; Jeong, B.R.; Lee, B.; Chang, Y.K. Isolation, phenotypic characterization and genome wide analysis of a Chlamydomonas reinhardtii strain naturally modified under laboratory conditions: Towards enhanced microalgal biomass and lipid production for biofuels. Biotechnol. Biofuels 2017, 10, 308. [Google Scholar] [CrossRef] [PubMed]
  20. Li, Q.; Fu, L.; Wang, Y.; Zhou, D.; Rittmann, B.E. Excessive phosphorus caused inhibition and cell damage during heterotrophic growth of Chlorella regularis. Bioresour. Technol. 2018, 268, 266–270. [Google Scholar] [CrossRef]
  21. Arora, N.; Patel, A.; Pruthi, P.A.; Pruthi, V. Synergistic dynamics of nitrogen and phosphorous influences lipid productivity in Chlorella minutissima for biodiesel production. Bioresour. Technol. 2016, 213, 79–87. [Google Scholar] [CrossRef] [PubMed]
  22. Li, X.; Manuel, J.; Slavens, S.; Crunkleton, D.W.; Johannes, T.W. Interactive effects of light quality and culturing temperature on algal cell size, biomass doubling time, protein content, and carbohydrate content. Appl. Microbiol. Biotechnol. 2021, 105, 587–597. [Google Scholar] [CrossRef]
  23. Wagner, I.; Steinweg, C.; Posten, C. Mono- and dichromatic LED illumination leads to enhanced growth and energy conversion for high-efficiency cultivation of microalgae for application in space. Biotechnol. J. 2016, 11, 1060–1071. [Google Scholar] [CrossRef] [PubMed]
  24. Kim, D.G.; Lee, C.; Park, S.M.; Choi, Y.E. Manipulation of light wavelength at appropriate growth stage to enhance biomass productivity and fatty acid methyl ester yield using Chlorella vulgaris. Bioresour. Technol. 2014, 159, 240–248. [Google Scholar] [CrossRef]
  25. Oldenhof, H.; Bisova, K.; van den Ende, H.; Zachleder, V. Effect of red and blue light on the timing of cyclin-dependent kinase activity and the timing of cell division in Chlamydomonas reinhardtii. Plant Physiol. Biochem. 2004, 42, 341–348. [Google Scholar] [CrossRef] [PubMed]
  26. Malerba, M.E.; Palacios, M.M.; Palacios Delgado, Y.M.; Beardall, J.; Marshall, D.J. Cell size, photosynthesis and the package effect: An artificial selection approach. New Phytol. 2018, 219, 449–461. [Google Scholar] [CrossRef]
  27. Li, J.; Zhang, K.; Li, L.; Wang, Y.; Lin, S. Unsuspected functions of alkaline phosphatase PhoD in the diatom Phaeodactylum tricornutum. Algal Res. 2022, 68, 102873. [Google Scholar] [CrossRef]
  28. Xu, L.; Wang, F.; Li, R.; Deng, M.; Fu, M.; Teng, H.; Yi, K. OsCYCP4s coordinate phosphate starvation signaling with cell cycle progression in rice. J. Integr. Plant Biol. 2020, 62, 1017–1033. [Google Scholar] [CrossRef]
  29. Deng, M.; Hu, B.; Xu, L.; Liu, Y.; Wang, F.; Zhao, H.; Wei, X.; Wang, J.; Yi, K. OsCYCP1;1, a PHO80 homologous protein, negatively regulates phosphate starvation signaling in the roots of rice (Oryza sativa L.). Plant Mol. Biol. 2014, 86, 655–669. [Google Scholar] [CrossRef]
  30. Bjorklund, M. Cell size homeostasis: Metabolic control of growth and cell division. Biochim. Biophys. Acta Mol. Cell Res. 2019, 1866, 409–417. [Google Scholar] [CrossRef]
  31. Wilhelm, C.; Jungandreas, A.; Jakob, T.; Goss, R. Light acclimation in diatoms: From phenomenology to mechanisms. Mar. Genom. 2014, 16, 5–15. [Google Scholar] [CrossRef] [PubMed]
  32. Malerba, M.E.; Marshall, D.J.; Palacios, M.M.; Raven, J.A.; Beardall, J. Cell size influences inorganic carbon acquisition in artificially selected phytoplankton. New Phytol. 2021, 229, 2647–2659. [Google Scholar] [CrossRef]
  33. Yan, D.; Beardall, J.; Gao, K. Variation in cell size of the diatom Coscinodiscus granii influences photosynthetic performance and growth. Photosynth. Res. 2018, 137, 41–52. [Google Scholar] [CrossRef] [PubMed]
  34. Slocombe, S.P.; Huete-Ortega, M.; Kapoore, R.V.; Okurowska, K.; Mair, A.; Day, J.G.; Stanley, M.S.; Vaidyanathan, S. Enabling large-scale production of algal oil in continuous output mode. iScience 2021, 24, 102743. [Google Scholar] [CrossRef] [PubMed]
  35. Smith, S.R.; Gillard, J.T.; Kustka, A.B.; McCrow, J.P.; Badger, J.H.; Zheng, H.; New, A.M.; Dupont, C.L.; Obata, T.; Fernie, A.R.; et al. Transcriptional Orchestration of the Global Cellular Response of a Model Pennate Diatom to Diel Light Cycling under Iron Limitation. PLoS Genet. 2016, 12, e1006490. [Google Scholar] [CrossRef]
  36. Obando-Montoya, E.J.; Zapata-Ocampo, P.A.; Cuesta-Astroz, Y.; Atehortua, L. Impact of light wavelength on the transcriptome of Porphyridium cruentum and culture yield. Algal Res. 2022, 67, 102856. [Google Scholar] [CrossRef]
  37. Montsant, A.; Allen, A.E.; Coesel, S.; Martino, A.D.; Falciatore, A.; Mangogna, M.; Siaut, M.; Heijde, M.; Jabbari, K.; Maheswari, U.; et al. Identification and Comparative Genomic Analysis of Signaling and Regulatory Components in the Diatomthalassiosira pseudonana. J. Phycol. 2007, 43, 585–604. [Google Scholar] [CrossRef]
  38. Persson, B.L.; Lagerstedt, J.O.; Pratt, J.R.; Pattison-Granberg, J.; Lundh, K.; Shokrollahzadeh, S.; Lundh, F. Regulation of phosphate acquisition in Saccharomyces cerevisiae. Curr. Genet. 2003, 43, 225–244. [Google Scholar] [CrossRef]
  39. Loha, A.; Kashyap, A.K.; Sharma, P. A putative cyclin, SiPHO80 from root endophytic fungus Serendipita indica regulates phosphate homeostasis, salinity and heavy metal toxicity tolerance. Biochem. Biophys. Res. Commun. 2018, 507, 414–419. [Google Scholar] [CrossRef]
  40. Cato, M.L.; Jester, H.D.; Lavertu, A.; Lyman, A.; Tallent, L.M.; Mitchell, G.C. Genome-Wide Analysis of Cell Cycle-Regulating Genes in the Symbiotic Dinoflagellate Breviolum minutum. G3 2019, 9, 3843–3853. [Google Scholar] [CrossRef]
  41. Chen, Y.; Sun, S.; Wang, X. The epidermis-specific cyclin CYCP3;1 is involved in the excess brassinosteroid signaling-inhibited root meristem cell division. J. Integr. Plant Biol. 2020, 62, 1674–1687. [Google Scholar] [CrossRef] [PubMed]
  42. Torres Acosta, J.A.; de Almeida Engler, J.; Raes, J.; Magyar, Z.; De Groodt, R.; Inze, D.; De Veylder, L. Molecular characterization of Arabidopsis PHO80-like proteins, a novel class of CDKA;1-interacting cyclins. Cell. Mol. Life Sci. 2004, 61, 1485–1497. [Google Scholar] [CrossRef] [PubMed]
  43. Peng, L.; Skylar, A.; Chang, P.L.; Bisova, K.; Wu, X. CYCP2;1 integrates genetic and nutritional information to promote meristem cell division in Arabidopsis. Dev. Biol. 2014, 393, 160–170. [Google Scholar] [CrossRef]
  44. Pecani, K.; Lieberman, K.; Tajima-Shirasaki, N.; Onishi, M.; Cross, F.R. Control of division in Chlamydomonas by cyclin B/CDKB1 and the anaphase-promoting complex. PLoS Genet. 2022, 18, e1009997. [Google Scholar] [CrossRef] [PubMed]
  45. Zheng, T.; Dai, L.; Liu, Y.; Li, S.; Zheng, M.; Zhao, Z.; Qu, G.Z. Overexpression Populus d-Type Cyclin Gene PsnCYCD1;1 Influences Cell Division and Produces Curved Leaf in Arabidopsis thaliana. Int. J. Mol. Sci. 2021, 22, 5837. [Google Scholar] [CrossRef]
  46. Zheng, T.; Dai, L.; Li, S.; Liu, Y.; Zhao, Z.; Yang, C.; Qu, G. Populus D-type cyclin gene PsnCYCD1;1 accelerates cell division and participates in secondary growth of vascular bundles. J. Exp. Bot. 2023, 74, 4077–4092. [Google Scholar] [CrossRef]
  47. Guan, C.; Xue, Y.; Jiang, P.; He, C.; Zhuge, X.; Lan, T.; Yang, H. Overexpression of PtoCYCD3;3 Promotes Growth and Causes Leaf Wrinkle and Branch Appearance in Populus. Int. J. Mol. Sci. 2021, 22, 1288. [Google Scholar] [CrossRef]
  48. Peramuna, A.; Lopez, C.Q.; Rios, F.J.A.; Bae, H.; Fangel, J.U.; Batth, R.; Harholt, J.; Simonsen, H.T. Overexpression of Physcomitrium patens cell cycle regulators leads to larger gametophytes. Sci. Rep. 2023, 13, 4301. [Google Scholar] [CrossRef]
  49. Kono, A.; Ohno, R.; Umeda-Hara, C.; Uchimiya, H.; Umeda, M. A distinct type of cyclin D, CYCD4;2, involved in the activation of cell division in Arabidopsis. Plant Cell Rep. 2006, 25, 540–545. [Google Scholar] [CrossRef]
  50. Jong, L.W.; Fujiwara, T.; Nozaki, H.; Miyagishima, S.Y. Cell size for commitment to cell division and number of successive cell divisions in multicellular volvocine green algae Tetrabaena socialis and Gonium pectorale. Proc. Jpn. Acad. Ser. B Phys. Biol. Sci. 2017, 93, 832–840. [Google Scholar] [CrossRef]
  51. Jong, L.W.; Fujiwara, T.; Hirooka, S.; Miyagishima, S.Y. Cell size for commitment to cell division and number of successive cell divisions in cyanidialean red algae. Protoplasma 2021, 258, 1103–1118. [Google Scholar] [CrossRef] [PubMed]
  52. Jones, A.R.; Band, L.R.; Murray, J.A.H. Double or Nothing? Cell Division and Cell Size Control. Trends Plant Sci. 2019, 24, 1083–1093. [Google Scholar] [CrossRef] [PubMed]
  53. Wellburn, A.R. The Spectral Determination of Chlorophylls a and b, as well as Total Carotenoids, Using Various Solvents with Spectrophotometers of Different Resolution. J. Plant Physiol. 1994, 144, 307–313. [Google Scholar] [CrossRef]
  54. Schreiber, U.; Klughammer, C.; Kolbowski, J. High-end chlorophyll fluorescence analysis with the MULTI-COLOR-PAM. I. Various light qualities and their applications. PAM Appl. Notes 2011, 1, 1–21. [Google Scholar]
  55. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef]
  56. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef]
  57. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  58. Thorvaldsdottir, H.; Robinson, J.T.; Mesirov, J.P. Integrative Genomics Viewer (IGV): High-performance genomics data visualization and exploration. Brief Bioinform. 2013, 14, 178–192. [Google Scholar] [CrossRef]
  59. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  60. Finn, R.D.; Clements, J.; Eddy, S.R. HMMER web server: Interactive sequence similarity searching. Nucleic Acids Res. 2011, 39, W29–W37. [Google Scholar] [CrossRef]
  61. Jones, P.; Binns, D.; Chang, H.Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [PubMed]
  62. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed]
  63. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [PubMed]
  64. Subramanian, B.; Gao, S.; Lercher, M.J.; Hu, S.; Chen, W.H. Evolview v3: A webserver for visualization, annotation, and management of phylogenetic trees. Nucleic Acids Res. 2019, 47, W270–W275. [Google Scholar] [CrossRef]
  65. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef]
  66. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Phenotypes of WT and ECS. (A) Growth curve of WT and ECS after inoculation. (B) Flow cytometry signal of mean forward scatter per cell. (C) Statistical analysis of cell area through fluorescence microscopy. ns indicates p > 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. Error bars indicate standard deviations. Cell density and mean FSC-A values were analyzed using two-way ANOVA (GraphPad Prism), while area was analyzed through t-test (GraphPad Prism).
Figure 1. Phenotypes of WT and ECS. (A) Growth curve of WT and ECS after inoculation. (B) Flow cytometry signal of mean forward scatter per cell. (C) Statistical analysis of cell area through fluorescence microscopy. ns indicates p > 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. Error bars indicate standard deviations. Cell density and mean FSC-A values were analyzed using two-way ANOVA (GraphPad Prism), while area was analyzed through t-test (GraphPad Prism).
Ijms 24 13595 g001
Figure 2. Chlorophyll fluorescence and chlorophyll content of WT and ECS. (A) Maximum PSII quantum yield (Fv/Fm). (B) Effective quantum yield of PSII (YII). (C) Chlorophyll a content per unit volume and per cell. ns indicates p > 0.05, * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. Error bars represent standard deviations calculated from three independent biological replicates. All data were analyzed through two-way ANOVA (GraphPad Prism).
Figure 2. Chlorophyll fluorescence and chlorophyll content of WT and ECS. (A) Maximum PSII quantum yield (Fv/Fm). (B) Effective quantum yield of PSII (YII). (C) Chlorophyll a content per unit volume and per cell. ns indicates p > 0.05, * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. Error bars represent standard deviations calculated from three independent biological replicates. All data were analyzed through two-way ANOVA (GraphPad Prism).
Ijms 24 13595 g002
Figure 3. Dry biomass, protein, carbohydrate, and lipid contents of WT and ECS. (A) Dry weight and dry weight per cell on day 14. (B) Protein, carbohydrate, and lipid contents on day 14. (C) Protein, carbohydrate, and lipid contents per cell on day 14. ns indicates p > 0.05, * p < 0.05, *** p < 0.001, **** p < 0.0001. Error bars represent standard deviations calculated from three independent biological replicates. All data were analyzed through t-test (GraphPad Prism).
Figure 3. Dry biomass, protein, carbohydrate, and lipid contents of WT and ECS. (A) Dry weight and dry weight per cell on day 14. (B) Protein, carbohydrate, and lipid contents on day 14. (C) Protein, carbohydrate, and lipid contents per cell on day 14. ns indicates p > 0.05, * p < 0.05, *** p < 0.001, **** p < 0.0001. Error bars represent standard deviations calculated from three independent biological replicates. All data were analyzed through t-test (GraphPad Prism).
Ijms 24 13595 g003
Figure 4. Identification and confirmation of insertion sites of WT and ECS. (A) Visualization of whole-genome sequencing analysis of insertion site by IGV. (B) Schematic diagram of insertion site. (C) Gene annotation of NO14G02130 using SMART and CDD. (D) Identification of inserted fragment through PCR amplification. (E) Gene expression levels of CYCU;1 between WT and ECS determined through RT-qPCR. Arrows indicate insertion sites. M indicates 5000 marker, **** p < 0.0001. UBCE gene as internal control and RT-qPCR were analyzed through t-test (GraphPad Prism).
Figure 4. Identification and confirmation of insertion sites of WT and ECS. (A) Visualization of whole-genome sequencing analysis of insertion site by IGV. (B) Schematic diagram of insertion site. (C) Gene annotation of NO14G02130 using SMART and CDD. (D) Identification of inserted fragment through PCR amplification. (E) Gene expression levels of CYCU;1 between WT and ECS determined through RT-qPCR. Arrows indicate insertion sites. M indicates 5000 marker, **** p < 0.0001. UBCE gene as internal control and RT-qPCR were analyzed through t-test (GraphPad Prism).
Ijms 24 13595 g004
Figure 5. Cyclin gene families of N. oceanica IMET1. (A) Phylogenetic tree was constructed based on full-length sequences of cyclin proteins. (B) Distribution of conserved protein motifs. (C) Heatmap of cyclin expression under light/dark cycles based on FPKM value. ZT0, ZT15, ZT18, and ZT21 correspond to dark phase; ZT3, ZT6, ZT9, and ZT12 correspond to light phase.
Figure 5. Cyclin gene families of N. oceanica IMET1. (A) Phylogenetic tree was constructed based on full-length sequences of cyclin proteins. (B) Distribution of conserved protein motifs. (C) Heatmap of cyclin expression under light/dark cycles based on FPKM value. ZT0, ZT15, ZT18, and ZT21 correspond to dark phase; ZT3, ZT6, ZT9, and ZT12 correspond to light phase.
Ijms 24 13595 g005
Figure 6. Phylogenetic tree analysis of cyclins in A. thaliana, O. sativa, N. oceanica IMET1, N. oceanica CCMP1779, C. reinhardtii, and P. tricornutum. Tree was divided into five clusters indicated by different colors.
Figure 6. Phylogenetic tree analysis of cyclins in A. thaliana, O. sativa, N. oceanica IMET1, N. oceanica CCMP1779, C. reinhardtii, and P. tricornutum. Tree was divided into five clusters indicated by different colors.
Ijms 24 13595 g006
Figure 7. Flow cytometry analysis of cell cycle. (A) DNA content per cell measured as DAPI fluorescence using flow cytometry at different times of day (ZT, hours after lights on). (B) Estimated percentage of cells in 1C (unreplicated haploid genome), S1 (first round of DNA synthesis), 2C (replicated haploid genome), S2 (second round of DNA synthesis), and 4C (cells with four haploid genomes). (C) Flow cytometry signal of mean forward scatter per cell at different times of day. (D) Flow cytometry signal of mean forward scatter per cell with an additional 24 h of culturing in the dark after sampling at the same time point. ns indicates p > 0.05, *** p < 0.001, **** p < 0.0001. All data were analyzed using two-way ANOVA (GraphPad Prism); error bars indicate standard deviations calculated from three independent biological replicates.
Figure 7. Flow cytometry analysis of cell cycle. (A) DNA content per cell measured as DAPI fluorescence using flow cytometry at different times of day (ZT, hours after lights on). (B) Estimated percentage of cells in 1C (unreplicated haploid genome), S1 (first round of DNA synthesis), 2C (replicated haploid genome), S2 (second round of DNA synthesis), and 4C (cells with four haploid genomes). (C) Flow cytometry signal of mean forward scatter per cell at different times of day. (D) Flow cytometry signal of mean forward scatter per cell with an additional 24 h of culturing in the dark after sampling at the same time point. ns indicates p > 0.05, *** p < 0.001, **** p < 0.0001. All data were analyzed using two-way ANOVA (GraphPad Prism); error bars indicate standard deviations calculated from three independent biological replicates.
Ijms 24 13595 g007
Figure 8. Expression levels of cell-cycle-related genes between WT and ECS. (A) Expression levels of cyclin genes of WT at different times of day compared with CYCP;3. (B) Gene expression levels of CYCU;1 of WT at different times of day compared with ZT0. (C) Heatmap of expression of cell-cycle-related genes at different times of day based on average using 2−ΔΔCT method. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. RT-qPCR data are from three independent biological replicates, and error bars indicate standard deviations. UBCE gene as an internal control and RT-qPCR were analyzed through t-test (GraphPad Prism).
Figure 8. Expression levels of cell-cycle-related genes between WT and ECS. (A) Expression levels of cyclin genes of WT at different times of day compared with CYCP;3. (B) Gene expression levels of CYCU;1 of WT at different times of day compared with ZT0. (C) Heatmap of expression of cell-cycle-related genes at different times of day based on average using 2−ΔΔCT method. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. RT-qPCR data are from three independent biological replicates, and error bars indicate standard deviations. UBCE gene as an internal control and RT-qPCR were analyzed through t-test (GraphPad Prism).
Ijms 24 13595 g008
Table 1. Characteristics of cyclin genes identified in Nannochloropsis oceanica.
Table 1. Characteristics of cyclin genes identified in Nannochloropsis oceanica.
Gene Name Gene ID in N.o IMET1Locus in
Genome
Gene FeaturesProtein FeaturesPfam (Position)SMART (Position)CDD (Position)Prosite
(Position)
CDS LengthNo. of
Introns
Protein LengthpIMWInstability
Index
GRAVY
CYCL;1NO01G01210chr1:315,224-316,537131404379.4848,434.2760.81−0.516 PF00134 (135-204)SM00385 (138-246)cd00043 (137-192)
CYCH;1NO02G04130chr2:1,126,722-1,128,36711851394643,296.1651.85−0.378PF00134 (142-214)
PF16899 (218-308)
SM00385 (115-208)cd00043 (141-207)
CYCT;1NO02G05390chr2:1,437,666-1,439,100105013497.0339,336.5654.86−0.445PF00134 (26-106)SM00385 (29-181, 194-336)cd00043 (31-110)
CYCP;1NO03G03610chr3:1,070,086-1,072,470219907325.1882,325.0574.32−0.839PF00134 (401-446) cd00043 (398-445)
CYCL;2NO03G03880chr3:1,139,762-1,141,883177905926.0364,464.1468.61−0.65PF00134 (71-206)SM00385 (80-198, 217-299)cd00043 (74-169, 214-298)
CYCP;2NO03G05460chr3:1,562,402-1,566,971234047795.9486,517.3369.46−1.086PF08613 (371-677)SM00385 (576-672)cd20540 (572-675)
CYCD;2NO05G01050chr5:313,215-314,34796013199.2234,130.2357.4−0.086PF00134 (52-145)
PF02984 (154-231)
SM00385 (59-145)
SM01332 (98-248)
cd00043 (54-143)PS00292 (54-85)
CYCD;1NO07G04620chr7:1,302,238-1,305,696248728289.2586,631.0266.54−0.358PF00134 (474-592)
PF02984 (603-739)
SM00385 (505-590, 607-705)
SM01332 (603-740)
cd00043 (499-589, 603-662)
CYCP;3NO10G00130chr10:40,340-43,3081539105127.7558,056.1374.76−0.77PF00134 (328-433)SM00385 (341-426)cd00043 (339-424)
CYCB;1NO12G01680chr12:502,593-504,248141624715.6151,922.1362.33−0.29PF00134 (227-350)
PF02984 (354-465)
SM00385 (260-344, 357-438)
SM01332 (353-468)
cd00043 (254-343, 351-436)PS00292 (255-286)
CYCP;4NO12G03480chr12:951,333-954,950195316506.371,251.3572.85−0.641PF00134 (333-439)SM00385 (349-434)cd00043 (347-433)
CYCA;1NO13G01690chr13:499,826-502,404165015497.5760,194.3567.26−0.412PF00134 (274-402)
PF02984 (405-534)
PF08613 (371-677)
SM00385 (312-396, 409-506)
SM01332 (405-537)
cd00043 (306-395, 406-504)
CYCU;1NO14G02130chr14:640,825-642,525170105668.2361,891.9767.84−0.483PF08613 (47-159)SM00385 (69-153)cd20558 (63-162)
CYCL;3NO15G01200chr15:342,170-343,33993923129.3134,799.3246.11−0.265PF00134 (16-139)SM00385 (33-132)cd00043 (28-105)
CYCA;2NO19G00890chr19:301,869-304,149221117364.7282,335.8160.24−0.593PF00134 (145-268)SM00385 (175-261)cd00043 (169-260)
CYCL;4NO22G00250chr22:89,639-94,300167795588.7460,996.2473.95−0.673PF00134 (150-269)SM00385 (163-263)
CYCB;2NO23G01320chr23:476,595-480,938253268435.5989,728.790.25−0.666PF00134 (524-630)SM00385 (537-625, 652-743)
SM01332 (634-778)
cd00043 (531-563PS00292 (532-563)
CYCB;3NO25G01730chr25:508,155-511,392166265535.7859,436.7556.39−0.45PF00134 (301-422)
PF02984 (425-541)
SM00385 (332-416, 429-510)
SM01332 (425-545)
cd00043 (326-415, 425-508)PS00292 (327-358)
CYCU;2NO28G01340chr28:421,745-424,03893333109.3134,082.2972.84−0.352PF08613 (173-282)SM00385 (190-276) PS00292 (185-216)
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

Xu, W.; Lin, Y.; Wang, Y.; Li, Y.; Zhu, H.; Zhou, H. Phenotypic Analysis and Molecular Characterization of Enlarged Cell Size Mutant in Nannochloropsis oceanica. Int. J. Mol. Sci. 2023, 24, 13595. https://doi.org/10.3390/ijms241713595

AMA Style

Xu W, Lin Y, Wang Y, Li Y, Zhu H, Zhou H. Phenotypic Analysis and Molecular Characterization of Enlarged Cell Size Mutant in Nannochloropsis oceanica. International Journal of Molecular Sciences. 2023; 24(17):13595. https://doi.org/10.3390/ijms241713595

Chicago/Turabian Style

Xu, Weinan, Yihua Lin, Yu Wang, Yanyan Li, Hongmei Zhu, and Hantao Zhou. 2023. "Phenotypic Analysis and Molecular Characterization of Enlarged Cell Size Mutant in Nannochloropsis oceanica" International Journal of Molecular Sciences 24, no. 17: 13595. https://doi.org/10.3390/ijms241713595

APA Style

Xu, W., Lin, Y., Wang, Y., Li, Y., Zhu, H., & Zhou, H. (2023). Phenotypic Analysis and Molecular Characterization of Enlarged Cell Size Mutant in Nannochloropsis oceanica. International Journal of Molecular Sciences, 24(17), 13595. https://doi.org/10.3390/ijms241713595

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