Genome sequencing, assembly, and annotation of the self-flocculating microalga Scenedesmus obliquus AS-6-11
BMC Genomics volume 21, Article number: 743 (2020)
Scenedesmus obliquus belongs to green microalgae and is widely used in aquaculture as feed, which is also explored for lipid production and bioremediation. However, genomic studies of this microalga have been very limited. Cell self-flocculation of microalgal cells can be used as a simple and economic method for harvesting biomass, and it is of great importance to perform genome-scale studies for the self-flocculating S. obliquus strains to promote their biotechnological applications.
We employed the Pacific Biosciences sequencing platform for sequencing the genome of the self-flocculating microalga S. obliquus AS-6-11, and used the MECAT software for de novo genome assembly. The estimated genome size of S. obliquus AS-6-11 is 172.3 Mbp with an N50 of 94,410 bp, and 31,964 protein-coding genes were identified. Gene Ontology (GO) and KEGG pathway analyses revealed 65 GO terms and 428 biosynthetic pathways. Comparing to the genome sequences of the well-studied green microalgae Chlamydomonas reinhardtii, Chlorella variabilis, Volvox carteri and Micractinium conductrix, the genome of S. obliquus AS-6-11 encodes more unique proteins, including one gene that encodes D-mannose binding lectin. Genes encoding the glycosylphosphatidylinositol (GPI)-anchored cell wall proteins, and proteins with fasciclin domains that are commonly found in cell wall proteins might be responsible for the self-flocculating phenotype, and were analyzed in detail. Four genes encoding both GPI-anchored cell wall proteins and fasciclin domain proteins are the most interesting targets for further studies.
The genome sequence of the self-flocculating microalgal S. obliquus AS-6-11 was annotated and analyzed. To our best knowledge, this is the first report on the in-depth annotation of the S. obliquus genome, and the results will facilitate functional genomic studies and metabolic engineering of this important microalga. The comparative genomic analysis here also provides new insights into the evolution of green microalgae. Furthermore, identification of the potential genes encoding self-flocculating proteins will benefit studies on the molecular mechanism underlying this phenotype for its better control and biotechnological applications as well.
Microalgae are widely studied for producing biofuels and mitigating greenhouse gas emissions . In addition, microalgae are also producers of various high-value biochemicals, such as lipids, proteins, polysaccharides, pigments, vitamins, and antioxidants . For economic bioproduction by microalgae, robust strains and optimized processes are both essential . Genome-scale studies of microalgae can provide in-depth information on intracellular metabolism from a global prospect [4,5,6], and benefit the development of robust microalgal strains and efficient processes. Therefore, studies on genome sequencing and annotation of microalgae have received increasing attention.
Due to small cell size, negative surface charge and low biomass concentration achieved during photosynthetic autotrophic culture, harvesting microalgae biomass from a large volume of culture medium is a great challenge for their biorefinery [7, 8]. Among various technologies developed for microalgal biomass recovery, gravity sedimentation facilitated by the flocculation of microalgal cells is more economically competitive . On the other hand, some microalgal strains can flocculate or aggregate spontaneously [9,10,11]. The self-flocculation of microalgal cells enables their harvest without the addition of exogenous flocculants, and biomass harvesting based on the self-flocculation of microalgal cells is more environmentally friendly compared to the flocculation of microalgal cells through physical and chemical methods, or induced by infochemicals from predators [8, 12].
Despite the progress in using cell flocculation of microalgae for biomass recovery, in-depth studies on mechanisms of microalgal cell flocculation are still lacking. So far, microalgal cell wall polysaccharide and other extracellular polymeric substances (EPS) containing sugars and proteins have been identified as flocculating agents [10, 11, 13], but genes involved in the process remain unexplored. Identification of genes involved in microalgal cell self-flocculation is of importance to develop robust strains with controlled flocculation phenotype for microalgal biorefinery, which can be enabled by the advancement of genetic engineering of microalgae .
Scenedesmus belongs to green microalgae and is useful for lipid and pigment production, wastewater treatment, heavy metal removal and CO2 fixation [14,15,16,17,18]. The self-flocculating microalga S. obliquus AS-6-1 showing great advantages in microalgal cell harvest and heavy metals adsorption was reported previously [10, 16]. However, genes encoding key protein(s) for the synthesis of flocculating agents in S. obliquus are still not clear. To date, the genome sequences of four Scenedesmus strains are available (NCBI BioProjects PRJNA498405, PRJNA394817, PRJNA394817 and PRJNA428298), but none of them have been annotated adequately. Three of these genomes lack annotation information, and the annotation of the remaining strain Scenedesmus sp. ARA is incomplete due to the large contig numbers of 4727 and a low N50 value of 37,561. Additionally, the three Scenedesmus genomes were sequenced by the second-generation sequencing (SGS) technology that might cause bias annotation due to shorter read lengths compared with Pacific Biosciences (Pacbio) sequencing technology .
In this study, we sequenced the self-flocculating S. obliquus AS-6-11 genome using the Pacbio technology and reported its genome assembly and annotation. We explored the metabolic potential of this microalga, and performed comparative genome analyses with the other four annotated microalgal genomes. We also comprehensively analyzed the cell wall proteins of S. obliquus AS-6-11 that might act as the flocculating agents for the self-flocculating phenotype. The knowledge obtained in this work can not only benefit understanding and control of the self-flocculation of microalgal cells, but will also provide insights for further genome-scale studies of S. obliquus and other related microalgae to explore their biotechnological potentials.
Morphological features and genome assembly
Cell self-flocculation of S. obliquus AS-6-11 was observed by SEM analysis. The microalgal cells are round and form aggregates through cell-cell contacts (Fig. 1), which is different from the other reported Scenedesmus strains that are in spindle shape .
The estimated genome size of S. obliquus AS-6-11 is 172.3 Mbp with 2772 contigs, and the N50 contig size is 94.4 kbp using MECAT for the genome assembly (Additional file 1: Table S1; NCBI BioProject ID: PRJNA593662). Results using the MECAT software showed a better assembly ability than that of SMRT Portal in S. obliquus AS-6-11, in which the contig numbers are 58.1% less, and the N50 value is 1.5-fold higher (Additional file 1: Table S1). The genome sizes of the released Scenedesmus strains [20,21,22,23,24] range from 23.4 to 208.0 Mbp (Table 1). Among the available results, the N50 contig sizes of S. obliquus AS-6-11 reported in this study and S. obliquus strain DOE0152z using Pacbio technology are significantly higher than the other Scenedesmus strains using SGS (Table 1). The N50 contig size of S. obliquus AS-6-11 is 1.2-fold and 10.7-fold higher than Scenedesmus sp. MC-1 and S. quadricauda LWG 002611, respectively. Besides, the GC content of Scenedesmus strains ranges from 52.0 to 63.2%, and S. obliquus AS-6-11 has the lowest GC content (Table 1). Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis showed that the assembly of S. obliquus AS-6-11 is 87.1% complete with 2168 BUSCO groups (Additional file 2).
A total of 31,964 protein-coding genes were predicted in the S. obliquus AS-6-11 genome (Table 2). The predicted gene number of S. obliquus AS-6-11 genome is dramatically higher than the other Scenedesmus strains (Table 1). According to the Non-redundant protein (NR), SWISS-PROT, and Pfam protein families databases, 19,847, 13,099, and 13,612 proteins were annotated, respectively (Table 2). The protein number annotated based on the NR database is the largest, which is 1.52-fold higher than that obtained based on the SWISS-PROT database. Besides, 65 GO terms and 428 pathways were predicted by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases in S. obliquus AS-6-11, respectively.
The top 20 GO terms and KEGG pathways enriched in gene function annotation of the S. obliquus AS-6-11 genome were illustrated in Fig. 2. The top 20 GO terms are mainly located in biological process (10) and cellular component (8), in which the cell, cell part, and organelle are the top three GO terms (Fig. 2a). The top 20 KEGG pathways are mainly related to genetic information processing (14), in which chromosome and associated proteins, membrane trafficking, and spliceosome are the top three KEGG pathways (Fig. 2b).
Comparative genomic analysis based on KEGG pathways
A total of 428 pathways were annotated in the S. obliquus AS-6-11 genome. In terms of lipid metabolism, the fewest genes (171) were annotated in S. obliquus AS-6-11, especially in glycerolipid metabolism, glycerophospholipid metabolism and arachidonic acid metabolism (Table 3). However, more genes related to fatty acid biosynthesis and elongation were identified in S. obliquus AS-6-11 than that in C. reinhardtii and V. carteri (Table 3). Moreover, genes in the carotenoid biosynthesis in S. obliquus AS-6-11 are the fewest.
Comparative genomic analysis of orthologous gene clusters
Comparing with the other four species, S. obliquus AS-6-11 has 15,879 gene clusters with 14,576 orthologous clusters and 1303 single-copy gene clusters (Fig. 3). There are 3357 overlapping orthologous gene clusters among the five microalgae. S. obliquus AS-6-11 has the most gene clusters and singletons (defined as the singleton genes for which no orthologs could be found in any of the other species ), and the number (8751) is 1.26-fold, 3.71-fold, 5.34-fold and 1.67-fold higher than that in C. reinhardtii, C. variabilis, M. conductrix and V. carteri, respectively (Fig. 3). Comparative orthologous gene cluster analysis also showed that the phylogenetic proximity of S. obliquus AS-6-11 is very similar to that of the other four microalgae (Additional file 3: Fig. S1).
Comparative genomic analysis based on gene families
A total of 3608 gene families were identified in S. obliquus AS-6-11, in which 136 unique gene families existed (Fig. 4). Both the total and unique gene families in S. obliquus AS-6-11 are more abundant than that in the other four microalgae (Fig. 4). The number of the unique gene families in S. obliquus AS-6-11 is 0.86, 1.19, 1.31 and 1.39-fold larger than C. reinhardtii, C. variabilis, M. conductrix and V. carteri, respectively (Fig. 4). In the S. obliquus AS-6-11 genome, the unique gene families include membrane protein (PF10160), red chlorophyll catabolite reductase (RCC reductase, PF06405), D-mannose binding lectin (PF01453), lipase maturation factor (PF06762), lipid-A-disaccharide synthetase (PF02684), thioesterase-like superfamily (PF13279) and so on. In addition, S. obliquus AS-6-11 and M. conductrix have the most common gene families (Fig. 4).
Analysis of the genome features related to cell self-flocculation
Cell self-flocculation of budding yeast Saccharomyces cerevisiae has been well-studied. The flocculation proteins, for example, Flo1p, Flo5p, Flo9p, and Flo10p, are cell wall proteins (CWPs) and also called lectin [26, 27]. GPI-anchor was reported as the common element in cell adhesion proteins and the GPI-anchored adhesins in yeast species of Candida albicans and S. cerevisiae are the well-known fungal adhesions . In S. obliquus AS-6-11, a total of 432 GPI-anchored CWPs are identified. Analysis of the top 10 GPI-anchored CWPs indicated that seven of them has the transmembrane region, and eight of them had the signal peptides (Table 4). The isoelectric point (pI) and molecular weight (Mw) of the GPI-anchored CWPs vary from 4.95 to 9.58 and 6.10 KDa to 78.84 KDa, respectively (Table 4).
Fasciclin (PF02469) is an extracellular domain (http://pfam.xfam.org/family/PF02469) that belongs to the ancient cell adhesion domain that is common to plants and animals. So far, fasciclin domain proteins have not been analyzed in microalgae. In the S. obliquus AS-6-11 genome, a total of 33 fasciclin domain proteins are identified, which are divided into three groups (Fig. 5a). Three main motifs are randomly distributed across the fasciclin domain proteins (Fig. 5b). The predicted pI values and Mw greatly differ among the fasciclin domain proteins (Additional file 4: Table S2). The subcellular localization prediction of fasciclin domain proteins indicated that most proteins have cytoplasmic (cyto) sites, and 15 of them have secreted (extr) sites (Additional file 4: Table S2). Further analysis of these 15 fasciclin domain proteins containing extr sites showed that six proteins are homologous to the reported fasciclin proteins of Monoraphidium neglectum (64.84%), Aquabacterium sp. (61.36%), Scenedesmus sp. Ki4 (48.09%), Pelomonas puraquae (46.94%) (Table 5). Additionally, two of the predicted proteins are annotated into the extracellular region part according to the GO database.
Combining analysis of GPI-anchored CWPs and fasciclin domain proteins, four fasciclin domain proteins were found to distribute in GPI-anchored CWPs (Fig. 6a; Additional file 5), in which one has two FAS1 domains (four repeated domains in the fasciclin I family of proteins), two have transmembrane regions, and one has signal peptide (Fig. 6a). Comparative genomic analysis of S. obliquus AS-6-11 and the other four microalgae species (C. reinhardtii, C. variabilis, M. conductrix and V. carteri) revealed no similar proteins to the four fasciclin domain proteins. We also performed comparative transcriptome analysis of S. obliquus AS-6-11 and the non-flocculating S. obliquus FSP-3, and the results showed that the four fasciclin domain protein-encoding genes (Fig. 6a) had transcription level in S. obliquus AS-6-11, but the transcription of these genes cannot be detected in S. obliquus FSP-3 (Additional file 6: Table S3).
The unique gene family D-mannose binding lectin was also analyzed (Additional file 5). One gene belongs to this unique gene family was identified, and the encoded protein has two conserved domains: CAP (cysteine-rich secretory proteins) domain and B_lectin (D-mannose binding lectin) domain. The putative D-mannose binding lectin of S. obliquus AS-6-11 is homologous to a secreted glycoprotein Pry1p of S. cerevisiae YJM693 (SGD ID: S000003615), and the identity is 58% (Fig. 6b). The similarity between Pry1p and D-mannose binding lectin attributes to the same CAP domain (Fig. 6b).
Genome feature of S. obliquus AS-6-11
We report here the genome sequence and annotation of S. obliquus AS-6-11, which is the first sequenced self-flocculating microalgal genome, and is also so far the most comprehensively annotated genome among the available S. obliquus genome information. The much larger gene numbers of S. obliquus AS-6-11 (Table 1) suggested its unique feature. Considering the relative completeness of S. obliquus AS-6-11 genome and the lack of S. obliquus genome annotation, S. obliquus AS-6-11 may serve as a model alga for supplying reference genome annotation and investigating the gene function, evolution, and biotechnology application of S. obliquus strains.
We obtained a larger contig N50 size of S. obliquus AS-6-11 genome than that of Scenedesmus sp. MC- and S. quadricauda, which may attribute to Pacbio sequencing technology and the assembly tool MECAT . The MECAT software showed higher computing efficiency with comparable or improved genome results than other current tools for de novo assembly of large genomes , and this study is the first application of this tool in the genome assembly of microalgae.
Comparative analysis of S. obliquus AS-6-11 genome with the other microalgae
Due to the lack of genome annotations of Scenedesmus strains, we compared the S. obliquus AS-6-11 genome with four represented green microalgae (Fig. 3; Fig. 4; Table 3). S. obliquus AS-6-11 has more singletons, unique gene families as well as additional KEGG pathways than the other four microalgae, revealing its special evolutionary status, genomic complexity, and metabolic characteristics. Comparative genomics analysis of orthologous clusters among multiple species is important for identifying the overlap among orthologous clusters that helps to elucidate the evolution and function of proteins . The similar phylogenetic proximity of S. obliquus AS-6-11 to the four green microalgae further supports that this species locates in the intermediate stage of green algae evolution .
Due to the potential applications of Scenedesmus strains in wastewater treatment and lipid production [31, 32], their lipid metabolism needs to be well-studied. S. obliquus AS-6-11 has fewer genes of lipid metabolism compared to S. quadricauda LWG002611 and the other four green microalgae, suggesting the difference of lipid biosynthesis or genome integrity among the strains. For example, the fewer genes of arachidonic acid metabolism in S. obliquus AS-6-11 may attribute to multiple genes encoding one enzyme in other microalgae (Table 3). Key enzymes involved in triacylglycerol (TAG) biosynthesis pathways and carbon fixation were identified based on the genome sequence of S. quadricauda LWG002611, providing targets for genetic and metabolic engineering to improve biofuel production and reconstruct the metabolic pathways of this strain . In the previous work, overexpression of type 2 diacylglycerol acyltransferase gene (DGTT1) of C. reinhardtii in S. obliquus CPC2 enhanced its lipid content by nearly two-fold . Overexpression of acetyl-CoA carboxylase gene in Scenedesmus sp. MC-1 increased its intracellular lipid by 28.6%, indicating its importance in the lipid metabolism of Scenedesmus . However, the endogenous genes involved in lipid biosynthesis in S. obliquus has not been investigated. The identified lipid biosynthetic genes in this work will facilitate improving lipid production in S. obliquus.
On the other hand, Scenedesmus can produce important pigments carotenoids as valuable products, and overexpression of synthetic phytoene synthase gene (PSY) in Scenedesmus sp. CPC2 increased β-carotene production to 30 mg g− 1-cell . So far no studies have been focused on the innate carotenoid biosynthetic genes in S. obliquus. Further investigation of the functional genome of Scenedesmus will benefit the metabolic engineering of this important microalga for carotenoids production.
Potential flocculation proteins in S. obliquus AS-6-11
Although bioflocculation has been widely accepted as a promising way to harvest microalgal biomass, studies on cell self-flocculation of microalgae are still very limited. In the previous study, cell wall polysaccharides were revealed to act as flocculating agents in S. obliquus AS-6-1 . In the recent study, we found that protease treatment led to the de-flocculation of S. obliquus AS-6-11 (data not shown), indicating their different mechanisms of cell self-flocculation. Therefore, we focused on the identification of flocculating proteins in S. obliquus AS-6-11.
Yeast cell flocculation has been widely used in the beer industry and biofuels production as a simple, economic and environmentally friendly way to harvest cells [35, 36]. The interaction of lectin-like receptors with adjacent mannose side chains in cell walls played a vital role in yeast flocculation . In addition to lectin-glycan interaction, glycan-glycan interactions also act a pivotal part in cell-cell adhesion, and the flocculation mechanism contributes to the self-interaction of Flo proteins in the coordination of Ca2+ . In our recent studies, mechanisms of cell flocculation in an ethanol-producing bacterium Zymomonas mobilis were also revealed, where cellulose was found to be important for the cell flocculation . Moreover, our previous studies [10, 11] found that cell wall polysaccharides play key roles in microalgal cell self-flocculation. Considering that mannose-specific lectin plays an important role in yeast flocculation , the D-mannose binding lectin family gene found in the S. obliquus AS-6-11 genome may be related to cell flocculation. However, we did not find the transcription of this gene in the transcriptome data of S. obliquus AS-6-11. In the water surface-floating microalga Chlorococcum sp. FFG039, one protein in the unique gene family (jacalin-like lectin domain, PF01419) was assumed to be related to biofilm formation with the help of the lectin domain . However, the gene encoding jacalin-like lectin domain protein was not found in the S. obliquus AS-6-11 genome, suggesting the different mechanisms of microalgal flocculation.
GPI-anchored proteins are widespread in eukaryotes for anchoring proteins to the extracellular surface of the plasma membrane, and are involved in multiple cellular functions . GPI-anchored CWP Flo1p in S. cerevisiae is the main flocculation protein, and the GPI-anchor is necessary for Flo1p to attach to the cell wall . However, the GPI-anchored CWPs in microalgae have not been analyzed. Therefore, the genome-scale analysis of GPI-anchored CWPs is of great importance. In the genome of the model plant Arabidopsis thaliana, 210 GPI-anchored proteins were identified, and most of these proteins were involved in the primary modification for targeting specific proteins to the cell surface for extracellular matrix remodeling and signaling . According to the GPI-anchor studies of S. cerevisiae and Arabidopsis, the identified cell wall proteins with GPI-anchor of S. obliquus AS-6-11 may be related to the attachment of proteins to the cell wall, and the potential flocculation proteins may be members among them.
Fasciclin 1 (FAS1) domain is an ancient motif in extracellular proteins widely exists in all kingdoms of life, and FAS1 proteins mediated the interactions between the cell surface and cell exterior . FAS1 domain protein of V. carteri that is homology to Drosophila fasciclin I was identified as the cell adhesion protein , and it also caused unicells of C. reinhardtii to flocculate . Therefore, the four proteins identified in this study (Fig. 6a), which are homologous to fasciclin domain-containing proteins of Synechocystis, D. melanogaster and Galdieria sulphuraria or hypothetical protein of M. neglectum and C. variabilis, are the most likely flocculation proteins. Although we have made great efforts to investigate the potential flocculating gene function by gene disruption, the genetic transformation method of S. obliquus AS-6-11 was not successful. Alternatively, the protein-encoding genes can also be tested in various microalgal species, which will be focused in future studies. It needs to point out that in addition to the function of cell adhesion and/or flocculation, the fasciclin domain proteins in microalgae also take parts in stress response and cell wall formation. The gene encoding one fasciclin domain protein (carotenoprotein) in microalga Ki-4 was overexpressed under salt, dehydration and high light stresses, showing its function in protecting cells against photooxidative stresses . In addition, the fasciclin-like arabinogalactan protein family in higher plant Eucalyptus grandis took parts in the growth and properties of the secondary cell wall . To the best of our knowledge, this is the first report on the analysis of fascilin domain proteins in microalgae. It will be interesting to further explore the functions of this category of proteins. The in-depth analysis of S. obliquus AS-6-11 genome can also provide a basis for functional genomic studies of other microalgae.
We present here the genome sequencing, annotation, and analysis of the self-flocculating microalga S. obliquus AS-6-11. Comparative genomic analysis between S. obliquus AS-6-11 and the other microalgae reveals its strain specificities, evolutionary status as well as metabolic characteristics. Through the analysis of the protein family, the GPI-anchored CWPs and fasciclin domain proteins were identified for the first time in microalgae, and four GPI-anchored CWPs with fasciclin domain are the most potential flocculation proteins for further studies. Taken together, the draft genome of S. obliquus AS-6-11 will provide a reliable reference for the microalgae genome studies, increasing the understanding of microalgal self-flocculation mechanisms for promoting the microalgal harvest, and benefit efficiency biorefinery using microalgae.
Microalgal strain and culture conditions
The self-flocculating S. obliquus AS-6-11 was isolated from the freshwater pond in the campus of National Cheng Kung University (22°99′74.29″N, 120°22′22.30″E) in southern Taiwan, and was preserved at National Cheng Kung University. Cells were cultured at 25 °C in BG11 medium with continuous illumination of 75 μmol m2 s− 1 and continuous air aeration.
The cell shape and surface morphology of S. obliquus AS-6-11 cells were observed by scanning electron microscope (SEM, Hitachi S-3400 N II). For SEM, cells during the exponential growth phase were harvested and fixed in 2.5% glutaraldehyde solution overnight. After that, the samples were washed two times using PBS buffer (10 mM, pH 7.4), dehydrated in different concentrations of ethanol solutions, and then referred to the method of Salim et al. .
Genome DNA preparation and quality assessment
Two hundred milliliter cells at the exponential growth phase (Day 6) were harvested. The genome DNA of S. obliquus AS-6-11 was extracted by EZ-10 Spin Column Plant Genomic DNA Purification Kit (NO. B518261, Sangon Biotech (Shanghai) Co., Ltd.). The DNA quality was monitored and controlled by Nanodrop (Thermo Scientific NanoDrop 2000) and DNA gel electrophoresis.
Genome sequencing and assembly
The obtained high-quality DNA was sequenced using the PacBio RS II System. The size of DNA libraries was 10–20 kb and 10 SMRT cells were sequenced. Filtered Pacbio subreads were de novo assembled by software MECAT  and SMRT Portal with recommending parameters, respectively. Genome assembly quality using MECAT was further assessed by BUSCO (3.0.2) . The following analyses were based on the MECAT assembly result.
Genome annotation and comparative transcriptome analysis
MAKER2 training and annotation pipeline  with AUGUST  and SNAP  was used for genome structural annotation of S. obliquus AS-6-11 according to the C. reinhardtii training set and S. obliquus AS-6-11 transcriptome data. The transcriptome sequencing was performed by Illumina Hiseq platform. Functional annotation was performed by BLASTp (2.7.1+)  according to NR (https://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/) and Swiss-Prot protein (https://ftp.ncbi.nlm.nih.gov/blast/db/) databases. GO (http://geneontology.org/docs/download-ontology/) and KEGG (http://www.genome.jp/kegg/ko.html) annotations were performed according to the analysis method of Tamanna Sharma and Rajinder Singh Chauhan . Gene family analysis was performed by InterProScan 5.36–75.0-64  using the amino acid sequences of S. obliquus AS-6-11, C. reinhardtii (NCBI accession number: ABCN00000000.2), C. variabilis (NCBI accession number: ADIC00000000.1), V. carteri (NCBI accession number: ACJH00000000.1) and M. conductrix (NCBI accession number: LHPF00000000.2). In addition, the orthologous gene clusters between these microalgal genomes were compared using OrthoVenn2 (https://orthovenn2.bioinfotoolkits.net/) .
For transcriptome analysis, S. obliquus AS-6-11 and S. obliquus FSP-3 cells grown for 48 h under the same condition were harvested at 6000 rpm for 5 min, and then washed three times with the sterilized water. The washed cells were stored at − 80 °C immediately and send to Novogene Co., Ltd. for initial sequencing and analysis.
Prediction of GPI-anchored cell wall proteins
GPI-anchored CWPs in S. obliquus AS-6-11 were predicted by the GPI-anchored protein predictor developed by the National Science and technology development agency of Thailand . The top 10 of predicted GPI-anchored CWPs with signal peptide were further analyzed using compute pI/Mw tools of ExPASy (https://web.expasy.org/compute_pi/) and SMART (http://smart.embl-heidelberg.de/). Protein subcellular localization sites were estimated by Protein Subcellular Localization Prediction (https://wolfpsort.hgc.jp/).
Analysis of fasciclin domain proteins
The fasciclin domain family proteins were extracted from protein family annotation files and further analyzed using compute pI/Mw tool of ExPASy, SMART, and Protein Subcellular Localization Prediction.
Phylogenetic analysis of fasciclin domain proteins
A phylogenetic tree of fasciclin domain proteins was constructed using MEGA 7 . The sequence alignment was performed using MUSCLE and the phylogenetic tree was constructed using the Maximum likelihood method with 1000 bootstrap replicates.
Analysis of the conserved motifs of fasciclin domain proteins
The motifs’ analysis of fasciclin domain proteins was performed by the online MEME website (http://meme-suite.org/tools/meme) with the default parameters.
Analysis of conserved domains and homologous protein
The analyses of conserved domains and homologous protein were performed using NCBI CD-Search (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) and NCBI-blastp choosing the UniProtKB/Swiss-Prot database (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastp&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome), respectively.
Availability of data and materials
The genome sequence information of S. obliquus AS-6-11 was submitted to NCBI with the accession number of PRJNA593662. The amino acid sequences of the potential flocculation proteins of S. obliquus AS-6-11 can be found in Additional file 5. The genome sequences of C. reinhardtii (NCBI accession number: ABCN00000000.2), C. variabilis (NCBI accession number: ADIC00000000.1), V. carteri (NCBI accession number: ACJH00000000.1) and M. conductrix (NCBI accession number: LHPF00000000.2) were downloaded from the related websites: https://www.ncbi.nlm.nih.gov/genome/?term=ABCN00000000.2%2C, https://www.ncbi.nlm.nih.gov/genome/?term=ADIC00000000.1, https://www.ncbi.nlm.nih.gov/genome/?term=ACJH00000000.1, and https://www.ncbi.nlm.nih.gov/genome/?term=LHPF00000000.2, respectively.
Benchmarking Universal Single-Copy Orthologs
cell wall proteins
- DGTT1 :
type 2 diacylglycerol acyltransferase gene
environmental information processing
extracellular polymeric substances
- FAS1 domain:
four repeated domains in the fasciclin I family of proteins
genetic information processing
- PSY :
synthetic phytoene synthase gene
signaling and cellular processes
Chen J, Li J, Dong W, Zhang X, Tyagi RD, Drogui P, Surampalli RY. The potential of microalgae in biodiesel production. Renew Sust Energ Rev. 2018;90:336–46.
Chew KW, Yap JY, Show PL, Suan NH, Juan JC, Ling TC, et al. Microalgae biorefinery: high value products perspectives. Bioresour Technol. 2017;229:53–62.
Georgianna DR, Mayfield SP. Exploiting diversity and synthetic biology for the production of algal biofuels. Nature. 2012;488:329–35.
Starkenburg SR, Kwon KJ, Jha RK, McKay C, Jacobs M, Chertkov O, et al. A pangenomic analysis of the Nannochloropsis organellar genomes reveals novel genetic variations in key metabolic genes. BMC Genomics. 2014;15(1):212.
Bagnato C, Prados MB, Franchini GR, Scaglia N, Miranda SE, Beligni MV. Analysis of triglyceride synthesis unveils a green algal soluble diacylglycerol acyltransferase and provides clues to potential enzymatic components of the chloroplast pathway. BMC Genomics. 2017;18(1):223.
Ogura A, Akizuki Y, Imoda H, Mineta K, Gojobori T, Nagai S. Comparative genome and transcriptome analysis of diatom, Skeletonema costatum, reveals evolution of genes for harmful algal bloom. BMC Genomics. 2018;19(1):1–12.
Chen CY, Yeh KL, Aisyah R, Lee DJ, Chang JS. Cultivation, photobioreactor design and harvesting of microalgae for biodiesel production: a critical review. Bioresour Technol. 2011;102:71–81.
Wan C, Alam MA, Zhao XQ, Zhang XY, Guo SL, Ho SH, et al. Current progress and future prospect of microalgal biomass harvest using various flocculation technologies. Bioresour Technol. 2015;184:251–7.
Salim S, Bosma R, Vermuë MH, Wijffels RH. Harvesting of microalgae by bio-flocculation. J Appl Phycol. 2011;23:849–55.
Guo SL, Zhao XQ, Wan C, Huang ZY, Yang YL, Alam AM, et al. Characterization of flocculating agent from the self-flocculating microalga Scenedesmus obliquus AS-6-1 for efficient biomass harvest. Bioresour Technol. 2013;145:285–9.
Alam MA, Wan C, Guo SL, Zhao XQ, Huang ZY, Yang YL, et al. Characterization of the flocculating agent from the spontaneously flocculating microalga Chlorella vulgaris JSC-7. J Biosci Bioeng. 2014;118(1):29–33.
Roccuzzo S, Couto N, Karunakaran E, Kapoore RV, Butler TO, Mukherjee J, et al. Metabolic insights into infochemicals induced colony formation and flocculation in Scenedesmus subspicatus unraveled by quantitative proteomics. Front Microbiol. 2020;11:792.
Aljuboori AHR, Uemura Y, Thanh NT. Flocculation and mechanism of self-flocculating lipid producer microalga Scenedesmus quadricauda for biomass harvesting. Biomass Bioenergy. 2016;93:38–42.
Fayyaz M, Chew KW, Show PL, Ling TC, Ng IS, Chang JS. Genetic engineering of microalgae for enhanced biorefinery capabilities. Biotechnol Adv. 2020;10:107554.
Ling Y, Sun LP, Wang SY, Lin CSK, Sun Z, Zhou ZG. Cultivation of oleaginous microalga Scenedesmus obliquus coupled with wastewater treatment for enhanced biomass and lipid production. Biochem Eng J. 2019;148:162–9.
Zhang XY, Zhao XQ, Wan C, Chen BL, Bai FW. Efficient biosorption of cadmium by the self-flocculating microalga Scenedesmus obliquus AS-6-1. Algal Res. 2016;16:427–33.
Sharma T, Chauhan RS. Comparative transcriptomics reveals molecular components associated with differential lipid accumulation between microalgal sp. Scenedesmus dimorphus and Scenedesmus quadricauda. Algal Res. 2016;19:109–22.
Ho SH, Li PJ, Liu CC, Chang JS. Bioprocess development on microalgae-based CO2 fixation and bioethanol production using Scenedesmus obliquus CNW-N. Bioresour Technol. 2013;145:142–9.
Rhoads A, Au KF. PacBio sequencing and its applications. Genom Proteom Bioinf. 2015;13(5):278–89.
Nelson DR, Chaiboonchoe A, Fu W, Hazzouri KM, Huang Z, Jaiswal A, et al. Potential for heightened sulfur-metabolic capacity in coastal subtropical microalgae. iScience. 2019;11:450–65.
Ma C, Ren H, Xing D, Xie G, Ren N, Liu B. Mechanistic understanding towards the effective lipid production of a microalgal mutant strain Scenedesmus sp. Z-4 by the whole genome bioinformation. J Hazard Mater. 2019;375:115–20.
Dasgupta CN, Nayaka S, Toppo K, Singh AK, Deshpande U, Mohapatra A. Draft genome sequence and detailed characterization of biofuel production by oleaginous microalga Scenedesmus quadricauda LWG002611. Biotechnol Biofuels. 2018;11(1):308.
Carreres BM, de Jaeger L, Springer J, Barbosa MJ, Breuer G, van den End EJ, et al. Draft genome sequence of the oleaginous green alga Tetradesmus obliquus UTEX 393. Genome Announc. 2017;5(3):e01449–16.
Starkenburg S, Polle J, Hovde B, Daligault H, Davenport K, Huang A, et al. Draft nuclear genome, complete chloroplast genome, and complete mitochondrial genome for the biofuel/bioproduct feedstock species Scenedesmus obliquus strain DOE0152z. Genome Announc. 2017;5(32):e00617.
Wang Y, Devin CD, Chen G, Gu YQ. OrthoVenn: a web server for genome wide comparison and annotation of orthologous clusters across multiple species. Nucleic Acids Res. 2015;43(W1):W78–84.
Bony M, Barre P, Blondin B. Distribution of the flocculation protein, flop, at the cell surface during yeast growth: the availability of flop determines the flocculation level. Yeast. 1998;14(1):25–35.
Bony M, Thines-Sempoux D, Barre P, Blondin B. Localization and cell surface anchoring of the Saccharomyces cerevisiae flocculation protein Flo1p. J Bacteriol. 1997;179(15):4929–36.
Lipke PN. What we do not know about fungal cell adhesion molecules. J Fungi. 2018;4(2):59.
Xiao CL, Chen Y, Xie SQ, Chen KN, Wang Y, Han Y, et al. MECAT: fast mapping, error correction, and de novo assembly for single-molecule sequencing reads. Nat Methods. 2017;14:1072–4.
Nedelcu AM, Lee RW, Lemieux C, Gray MW, Burger G. The complete mitochondrial DNA sequence of Scenedesmus obliquus reflects an intermediate stage in the evolution of the green algal mitochondrial genome. Genome Res. 2000;10(6):819–31.
Tripathi R, Gupta A, Thakur IS. An integrated approach for phycoremediation of wastewater and sustainable biodiesel production by green microalgae, Scenedesmus sp. ISTGA1. Renew Energ. 2019;135:617–25.
Ye S, Gao L, Zhao J, An M, Wu H, Li M. Simultaneous wastewater treatment and lipid production by Scenedesmus sp. HXY2. Bioresour Technol. 2020;302:122903.
Chen CY, Kao AL, Tsai ZC, Chow TJ, Chang HY, Zhao XQ, et al. Expression of type 2 diacylglycerol acyltransferase gene DGTT1 from Chlamydomonas reinhardtii enhances lipid production in Scenedesmus obliquus. Biotechnol J. 2016;11(3):336–44.
Chen CY, Kao AL, Tsai ZC, Shen YM, Kao PH, Ng IS, et al. Expression of synthetic phytoene synthase gene to enhance β-carotene production in Scenedesmus sp. CPC2. J Biotechnol. 2017;12(11):1700204.
Zhao XQ, Bai FW. Yeast flocculation: new story in fuel ethanol production. Biotechnol Adv. 2009;27(6):849–56.
Cilliers JJ, Harrison STL. Yeast flocculation aids the performance of yeast dewatering using mini-hydrocyclones. Sep Purif Techn. 2019;209:159–63.
Nayyar A, Walker G, Wardrop F, Adya AK. Flocculation in industrial strains of Saccharomyces cerevisiae: role of cell wall polysaccharides and lectin-like receptors. J I Brewing. 2017;123(2):211–8.
Goossens KVY, Ielasi FS, Nookaew I, Stals I, Alonso-Sarduy L, Daenen L, et al. Molecular mechanism of flocculation self-recognition in yeast and its role in mating and survival. Mbio. 2015;6(2):e00427–15.
Xia J, Liu CG, Zhao XQ, Yi X, Xia XX, Bai FW. Contribution of cellulose synthesis, formation of fibrils and their entanglement to the self-flocculation of Zymomonas mobilis. Biotechnol Bioeng. 2018;115(11):2714–25.
Javadekar VS, Sivaraman H, Sainkar SR, Khan MI. A mannose-binding protein from the cell surface of flocculent Saccharomyces cerevisiae (NCIM 3528): its role in flocculation. Yeast. 2000;16(2):99–110.
Maeda Y, Nojima D, Sakurai M, Nomaguchi T, Ichikawa M, Ishizuka Y, et al. Genome analysis and genetic transformation of a water surface-floating microalga Chlorococcum sp. FFG039. Sci Rep-UK. 2019;9(1):1–7.
Fankhauser N, Mäser P. Identification of GPI anchor attachment signals by a Kohonen self-organizing map. Bioinformatics. 2005;21(9):1846–52.
Borner GHH, Sherrier DJ, Stevens TJ, Arkin IT, Dupree P. Prediction of glycosylphosphatidylinositol-anchored proteins in Arabidopsis. A genomic analysis. Plant Physiol. 2002;129(2):486–99.
Georg S. Fascinating fasciclins: a surprisingly widespread family of proteins that mediate interactions between the cell exterior and the cell surface. Int J Mol Sci. 2018;19(6):1628.
Huber O, Sumper M. Algal-CAMs: isoforms of a cell adhesion molecule in embryos of the alga Volvox with homology to Drosophila fasciclin I. EMBO J. 1994;13(18):4212–22.
Lowder LG, Herbert SK. Heterologous expression of a Volvox cell adhesion molecule causes flocculation in Chlamydomonas reinhardtii. J Appl Phycol. 2015;27(2):721–31.
Kawasaki S, Mizuguchi K, Sato M, Kono T, Shimizu H. A novel astaxanthin-binding photooxidative stress-inducible aqueous carotenoprotein from a eukaryotic microalga isolated from asphalt in midsummer. Plant Cell Physiol. 2013;54(7):1027–40.
MacMillan CP, Taylor L, Bi Y, Southerton SG, Evans R, Spokevicius A. The fasciclin-like arabinogalactan protein family of Eucalyptus grandis contains members that impact wood biology and biomechanics. New Phytol. 2015;206(4):1314–27.
Salim S, Kosterink NR, Wacka NDT, Vermuë MH, Wijffels RH. Mechanism behind autoflocculation of unicellular green microalgae Ettlia texensis. J Biotechnol. 2014;174:34–8.
Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31(19):3210–2.
Holt C, Yandell M. MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinformatics. 2011;12(1):491.
Stanke M, Diekhans M, Baertsch R, Haussler D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008;24(5):637–44.
Korf I. Gene finding in novel genomes. BMC Bioinformatics. 2004;5(1):59.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
Philip J, David B, Hsin-Yu C, Matthew F, Weizhong L, Craig MA, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40.
Phienluphon A, Mhuantong W, Boonyapakron K, Deenarn P, Champreda V, Wichadakul D, et al. Identification and evaluation of novel anchoring proteins for cell surface display on Saccharomyces cerevisiae. Appl Microbiol Biot. 2019;103(7):3085–97.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
The authors are thankful to Dr. Zhongqu Duan at Shanghai Jiao Tong University for genome assembly and annotation. We also appreciate Dr. Chun Wan and Professor Fantao Kong at the Dalian University of Technology for helpful discussion.
Ethics approval and consent to participate
The self-flocculating S. obliquus AS-6-11 was isolated from the freshwater pond on the campus of National Cheng Kung University, Taiwan by researchers at Prof. Jo-Shu Chang’s lab, and no permission is required for sample collection to isolate microalgal strains.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Genomic features of S. obliquus AS-6-11 using MECAT and SMRT Portal. (DOCX 14 kb)
Assessing genome completeness with BUSCO. (TXT 880 bytes)
The pairwise heatmap of the overlapping cluster numbers between the pair-wise genomes. (DOCX 55 kb)
Analysis of the fasciclin domain proteins in S. obliquus AS-6-11. (DOCX 20 kb)
The amino acid sequences of the potential flocculation proteins of S. obliquus AS-6-11. (TXT 2 kb)
Cq value for the potential flocculating genes in S. obliquus AS-6-11 compared to S. obliquus FSP-3. (DOCX 14 kb)
About this article
Cite this article
Chen, BL., Mhuantong, W., Ho, SH. et al. Genome sequencing, assembly, and annotation of the self-flocculating microalga Scenedesmus obliquus AS-6-11. BMC Genomics 21, 743 (2020). https://doi.org/10.1186/s12864-020-07142-4