Nannochloropsis plastid and mitochondrial phylogenomes reveal organelle diversification mechanism and intragenus phylotyping strategy in microalgae

Conjunctive BSA-Seq and BSR-Seq to Map the Genes of Yellow Leaf Mutations in Hot Peppers (Capsicum annuum L.)

Yellow leaf mutations have been widely used to study the chloroplast structures, the pigment synthesis, the photosynthesis mechanisms and the chlorophyll biosynthesis pathways across various species. For this study, a spontaneous mutant with the yellow leaf color named 96-140YBM was employed to explore the primary genetic elements that lead to the variations in the leaf color of hot peppers. To identify the pathways and genes associated with yellow leaf phenotypes, we applied sequencing-based Bulked Segregant Analysis (BSA-Seq) combined with BSR-Seq. We identified 4167 differentially expressed genes (DEGs) in the mutant pool compared with the wild-type pool. The results indicated that DEGs were involved in zeatin biosynthesis, plant hormone signal transduction, signal transduction mechanisms, post-translational modification and protein turnover. A total of 437 candidates were identified by the BSA-Seq, while the BSR-Seq pinpointed four candidate regions in chromosomes 8 and 9, containing 222 candidate genes. Additionally, the combination of BSA-Seq and BSR-Seq showed that there were 113 overlapping candidate genes between the two methods, among which 8 common candidates have been previously reported to be related to the development of chloroplasts, the photomorphogenesis and chlorophyll formation of plant chloroplasts and chlorophyll biogenesis. qRT-PCR analysis of the 8 common candidates showed higher expression levels in the mutant pool compared with the wild-type pool. Among the overlapping candidates, the DEG analysis showed that the CaKAS2 and CaMPH2 genes were down-regulated in the mutant pool compared to the wild type, suggesting that these genes may be key contributors to the yellow leaf phenotype of 96-140YBM. This research will deepen our understanding of the genetic basis of leaf color formation and provide valuable information for the breeding of hot peppers with diverse leaf colors.

Comparative structure and evolution of the organellar genomes of Padina usoehtunii (Dictyotales) with the brown algal crown radiation clade

BACKGROUND

Organellar genomes have become increasingly essential for studying genetic diversity, phylogenetics, and evolutionary histories of seaweeds. The order Dictyotales (Dictyotophycidae), a highly diverse lineage within the Phaeophyceae, is long-term characterized by a scarcity of organellar genome datasets compared to orders of the brown algal crown radiation (Fucophycidae).

RESULTS

We sequenced the organellar genomes of Padina usoehtunii, a representative of the order Dictyotales, to investigate the structural and evolutionary differences by comparing to five other major brown algal orders. Our results confirmed previously reported findings that the rate of structural rearrangements in chloroplast genomes is higher than that in mitochondria, whereas mitochondrial sequences exhibited a higher substitution rate compared to chloroplasts. Such evolutionary patterns contrast with land plants and green algae. The expansion and contraction of the inverted repeat (IR) region in the chloroplast correlated with the changes in the number of boundary genes. Specifically, the size of the IR region influenced the position of the boundary gene rpl21, with complete rpl21 genes found within the IR region in Dictyotales, Sphacelariales and Ectocarpales, while the rpl21 genes in Desmarestiales, Fucales, and Laminariales span both the IR and short single copy (SSC) regions. The absence of the rbcR gene in the Dictyotales may indicate an endosymbiotic transfer from the chloroplast to the nuclear genome. Inversion of the SSC region occurred at least twice in brown algae. Once in a lineage only represented by the Ectocarpales in the present study and once in a lineage only represented by the Fucales. Photosystem genes in the chloroplasts experienced the strongest signature of purifying selection, while ribosomal protein genes in both chloroplasts and mitochondria underwent a potential weak purifying selection.

CONCLUSIONS

Variations in chloroplast genome structure among different brown algal orders are evolutionarily linked to their phylogenetic positions in the Phaeophyceae tree. Chloroplast genomes harbor more structural rearrangements than the mitochondria, despite mitochondrial genes exhibiting faster mutation rates. The position and the change in the number of boundary genes likely shaped the IR regions in the chloroplast, and the produced structural variability is important mechanistically to create gene diversity in brown algal chloroplast.

Complete mitochondrial genome of an oleaginous microalga Vischeria punctata (Eustigmatophyceae: Chlorobotryaceae) and phylogenetic analysis

Vischeria punctata, as first described by Vischer in 1945, is a member of the family Chlorobotryaceae, within the order Eustigmatales. This species is recognized for its potential as a source of biofuels and other high-value products. In the present investigation, the whole genome of V. punctata was sequenced utilizing the Illumina HiSeq 4000 platform, enabling the assembly and annotation of its complete mitochondrial genome. The resulting circular genome spans 41,528 base pairs (bp) with a guanine-cytosine (GC) content of 27.3%. This genome encompasses 36 protein-coding genes, alongside 28 transfer RNA (tRNA), and three ribosomal RNA (rRNA) genes. The evolutionary trajectory of V. punctata was further explored by constructing a phylogenetic tree derived from the mitochondrial 33 gene dataset of 16 Ochrophyta species. Comparative analysis reveals that V. punctata bears closer ties to Vischeria sp. CAUP Q202 than to Vischeria stellata strain SAG 33.83, suggesting shared evolutionary pathways and phenotypic traits. This investigation constitutes the inaugural study into the mitochondrial evolution and phylogenetic patterning of the mitogenome in V. punctata. The outcomes from this research bolster our understanding of the genetic diversity and evolutionary processes within the class Eustigmatophyceae. In particular, the mitochondrial genome of V. punctata serves as a valuable resource in elucidating these aspects.

A k-mer-Based Approach for Phylogenetic Classification of Taxa in Environmental Genomic Data

In the age of genome sequencing, whole-genome data is readily and frequently generated, leading to a wealth of new information that can be used to advance various fields of research. New approaches, such as alignment-free phylogenetic methods that utilize k-mer-based distance scoring, are becoming increasingly popular given their ability to rapidly generate phylogenetic information from whole-genome data. However, these methods have not yet been tested using environmental data, which often tends to be highly fragmented and incomplete. Here, we compare the results of one alignment-free approach (which utilizes the D2 statistic) to traditional multi-gene maximum likelihood trees in 3 algal groups that have high-quality genome data available. In addition, we simulate lower-quality, fragmented genome data using these algae to test method robustness to genome quality and completeness. Finally, we apply the alignment-free approach to environmental metagenome assembled genome data of unclassified Saccharibacteria and Trebouxiophyte algae, and single-cell amplified data from uncultured marine stramenopiles to demonstrate its utility with real datasets. We find that in all instances, the alignment-free method produces phylogenies that are comparable, and often more informative, than those created using the traditional multi-gene approach. The k-mer-based method performs well even when there are significant missing data that include marker genes traditionally used for tree reconstruction. Our results demonstrate the value of alignment-free approaches for classifying novel, often cryptic or rare, species, that may not be culturable or are difficult to access using single-cell methods, but fill important gaps in the tree of life.

Identification and characterization of CsSRP43, a major gene controlling leaf yellowing in cucumber

Mutants are crucial to extending our understanding of genes and their functions in higher plants. In this study a spontaneous cucumber mutant, yf, showed yellow color leaves, had significant decreases in related physiological indexes of photosynthesis characteristics, and had more abnormal chloroplasts and thylakoids. Inheritance analysis indicated that the yellow color of the leaf was controlled by a recessive nuclear locus, yf. A candidate gene, CsSRP43, encoding a chloroplast signal recognition particle 43 protein, was identified through map-based cloning and whole-genome sequence analysis. Alignment of the CsSRP43 gene homologs between both parental lines revealed a 7-kb deletion mutation including the promoter region and the coding sequence in the yf mutant. In order to determine if the CsSRP43 gene was involved in the formation of leaf color, the CRISPR/Cas9-mediate system was used to modify CsSRP43 in the 9930 background; two independent transgenic lines, srp43-1 and srp43-2, were generated, and they showed yellow leaves with abnormal chloroplasts and thylakoids. Transcriptomic analysis revealed that differentially expressed genes associated with the photosynthesis-related pathway were highly enriched between srp43-1 and wild type, most of which were significantly downregulated in line srp43-1. Furthermore, yeast two-hybrid and biomolecular fluorescence complementation assays were used to confirm that CsSRP43 directly interacted with LHCP and cpSRP54 proteins. A model was established to explain the molecular mechanisms by which CsSRP43 participates in the leaf color and photosynthesis pathway, and it provides a valuable basis for understanding the molecular and genetic mechanisms of leaf color in cucumber.

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

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

Large Differences in the Haptophyte Phaeocystis globosa Mitochondrial Genomes Driven by Repeat Amplifications

The haptophyte Phaeocystis globosa is a well-known species for its pivotal role in global carbon and sulfur cycles and for its capability of forming harmful algal blooms (HABs) with serious ecological consequences. Its mitochondrial genome (mtDNA) sequence has been reported in 2014 but it remains incomplete due to its long repeat sequences. In this study, we constructed the first full-length mtDNA of P. globosa, which was a circular genome with a size of 43,585 bp by applying the PacBio single molecular sequencing method. The mtDNA of this P. globosa strain (CNS00066), which was isolated from the Beibu Gulf, China, encoded 19 protein-coding genes (PCGs), 25 tRNA genes, and two rRNA genes. It contained two large repeat regions of 6.7 kb and ∼14.0 kb in length, respectively. The combined length of these two repeat regions, which were missing from the previous mtDNA assembly, accounted for almost half of the entire mtDNA and represented the longest repeat region among all sequenced haptophyte mtDNAs. In this study, we tested the hypothesis that repeat unit amplification is a driving force for different mtDNA sizes. Comparative analysis of mtDNAs of five additional P. globosa strains (four strains obtained in this study, and one strain previously published) revealed that all six mtDNAs shared identical numbers of genes but with dramatically different repeat regions. A homologous repeat unit was identified but with hugely different numbers of copies in all P. globosa strains. Thus, repeat amplification may represent an important driving force of mtDNA evolution in P. globosa.

Genome engineering of Nannochloropsis with hundred-kilobase fragment deletions by Cas9 cleavages

Industrial microalgae are promising photosynthetic cell factories, yet tools for large-scale targeted genome engineering are limited. Here for the model industrial oleaginous microalga Nannochloropsis oceanica, we established a method to precisely and serially delete large genome fragments of ~100 kb from its 30.01 Mb nuclear genome. We started by identifying the 'non-essential' chromosomal regions (i.e. low expression region or LER) based on minimal gene expression under N-replete and N-depleted conditions. The largest such LER (LER1) is ~98 kb in size, located near the telomere of the 502.09-kb-long Chromosome 30 (Chr 30). We deleted 81 kb and further distal and proximal deletions of up to 110 kb (21.9% of Chr 30) in LER1 by dual targeting the boundaries with the episome-based CRISPR/Cas9 system. The telomere-deletion mutants showed normal telomeres consisting of CCCTAA repeats, revealing telomere regeneration capability after losing the distal part of Chr 30. Interestingly, the deletions caused no significant alteration in growth, lipid production or photosynthesis (transcript-abundance change for < 3% genes under N depletion). We also achieved double-deletion of both LER1 and LER2 (from Chr 9) that total ~214 kb at maximum, which can result in slightly higher growth rate and biomass productivity than the wild-type. Therefore, loss of the large, yet 'non-essential' regions does not necessarily sacrifice important traits. Such serial targeted deletions of large genomic regions had not been previously reported in microalgae, and will accelerate crafting minimal genomes as chassis for photosynthetic production.

Occurrence, Evolution and Specificities of Iron-Sulfur Proteins and Maturation Factors in Chloroplasts from Algae

Iron-containing proteins, including iron-sulfur (Fe-S) proteins, are essential for numerous electron transfer and metabolic reactions. They are present in most subcellular compartments. In plastids, in addition to sustaining the linear and cyclic photosynthetic electron transfer chains, Fe-S proteins participate in carbon, nitrogen, and sulfur assimilation, tetrapyrrole and isoprenoid metabolism, and lipoic acid and thiamine synthesis. The synthesis of Fe-S clusters, their trafficking, and their insertion into chloroplastic proteins necessitate the so-called sulfur mobilization (SUF) protein machinery. In the first part, we describe the molecular mechanisms that allow Fe-S cluster synthesis and insertion into acceptor proteins by the SUF machinery and analyze the occurrence of the SUF components in microalgae, focusing in particular on the green alga Chlamydomonas reinhardtii. In the second part, we describe chloroplastic Fe-S protein-dependent pathways that are specific to Chlamydomonas or for which Chlamydomonas presents specificities compared to terrestrial plants, putting notable emphasis on the contribution of Fe-S proteins to chlorophyll synthesis in the dark and to the fermentative metabolism. The occurrence and evolutionary conservation of these enzymes and pathways have been analyzed in all supergroups of microalgae performing oxygenic photosynthesis.

Comparative Plastid Genomics of Cryptomonas Species Reveals Fine-Scale Genomic Responses to Loss of Photosynthesis

Loss of photosynthesis is a recurring theme in eukaryotic evolution. In organisms that have lost the ability to photosynthesize, nonphotosynthetic plastids are retained because they play essential roles in processes other than photosynthesis. The unicellular algal genus Cryptomonas contains both photosynthetic and nonphotosynthetic members, the latter having lost the ability to photosynthesize on at least three separate occasions. To elucidate the evolutionary processes underlying the loss of photosynthesis, we sequenced the plastid genomes of two nonphotosynthetic strains, Cryptomonas sp. CCAC1634B and SAG977-2f, as well as the genome of the phototroph Cryptomonas curvata CCAP979/52. These three genome sequences were compared with the previously sequenced plastid genome of the nonphotosynthetic species Cryptomonas paramecium CCAP977/2a as well as photosynthetic members of the Cryptomonadales, including C. curvata FBCC300012D. Intraspecies comparison between the two C. curvata strains showed that although their genome structures are stable, the substitution rates of their genes are relatively high. Although most photosynthesis-related genes, such as the psa and psb gene families, were found to have disappeared from the nonphotosynthetic strains, at least ten pseudogenes are retained in SAG977-2f. Although gene order is roughly shared among the plastid genomes of photosynthetic Cryptomonadales, genome rearrangements are seen more frequently in the smaller genomes of the nonphotosynthetic strains. Intriguingly, the light-independent protochlorophyllide reductase comprising chlB, L, and N is retained in nonphotosynthetic SAG977-2f and CCAC1634B. On the other hand, whereas CCAP977/2a retains ribulose-1,5-bisphosphate carboxylase/oxygenase-related genes, including rbcL, rbcS, and cbbX, the plastid genomes of the other two nonphotosynthetic strains have lost the ribulose-1,5-bisphosphate carboxylase/oxygenase protein-coding genes.

The NanDeSyn database for Nannochloropsis systems and synthetic biology

Nannochloropsis species, unicellular industrial oleaginous microalgae, are model organisms for microalgal systems and synthetic biology. To facilitate community-based annotation and mining of the rapidly accumulating functional genomics resources, we have initiated an international consortium and present a comprehensive multi-omics resource database named Nannochloropsis Design and Synthesis (NanDeSyn; http://nandesyn.single-cell.cn). Via the Tripal toolkit, it features user-friendly interfaces hosting genomic resources with gene annotations and transcriptomic and proteomic data for six Nannochloropsis species, including two updated genomes of Nannochloropsis oceanica IMET1 and Nannochloropsis salina CCMP1776. Toolboxes for search, Blast, synteny view, enrichment analysis, metabolic pathway analysis, a genome browser, etc. are also included. In addition, functional validation of genes is indicated based on phenotypes of mutants and relevant bibliography. Furthermore, epigenomic resources are also incorporated, especially for sequencing of small RNAs including microRNAs and circular RNAs. Such comprehensive and integrated landscapes of Nannochloropsis genomics and epigenomics will promote and accelerate community efforts in systems and synthetic biology of these industrially important microalgae.

Optimal Growth Temperature and Intergenic Distances in Bacteria, Archaea, and Plastids of Rhodophytic Branch

The lengths of intergenic regions between neighboring genes that are convergent, divergent, or unidirectional were calculated for plastids of the rhodophytic branch and complete archaeal and bacterial genomes. Statistically significant linear relationships between any pair of the medians of these three length types have been revealed in each genomic group. Exponential relationships between the optimal growth temperature and each of the three medians have been revealed as well. The leading coefficients of the regression equations relating all pairs of the medians as well as temperature and any of the medians have the same sign and order of magnitude. The results obtained for plastids, archaea, and bacteria are also similar at the qualitative level. For instance, the medians are always low at high temperatures. At low temperatures, the medians tend to statistically significant greater values and scattering. The original model was used to test our hypothesis that the intergenic distances are optimized in particular to decrease the competition of RNA polymerases within the locus that results in transcribing shortened RNAs. Overall, this points to an effect of temperature for both remote and close genomes.

Plastid Genomes and Proteins Illuminate the Evolution of Eustigmatophyte Algae and Their Bacterial Endosymbionts

Eustigmatophytes, a class of stramenopile algae (ochrophytes), include not only the extensively studied biotechnologically important genus Nannochloropsis but also a rapidly expanding diversity of lineages with much less well characterized biology. Recent discoveries have led to exciting additions to our knowledge about eustigmatophytes. Some proved to harbor bacterial endosymbionts representing a novel genus, Candidatus Phycorickettsia, and an operon of unclear function (ebo) obtained by horizontal gene transfer from the endosymbiont lineage was found in the plastid genomes of still other eustigmatophytes. To shed more light on the latter event, as well as to generally improve our understanding of the eustigmatophyte evolutionary history, we sequenced plastid genomes of seven phylogenetically diverse representatives (including new isolates representing undescribed taxa). A phylogenomic analysis of plastid genome-encoded proteins resolved the phylogenetic relationships among the main eustigmatophyte lineages and provided a framework for the interpretation of plastid gene gains and losses in the group. The ebo operon gain was inferred to have probably occurred within the order Eustigmatales, after the divergence of the two basalmost lineages (a newly discovered hitherto undescribed strain and the Pseudellipsoidion group). When looking for nuclear genes potentially compensating for plastid gene losses, we noticed a gene for a plastid-targeted acyl carrier protein that was apparently acquired by horizontal gene transfer from Phycorickettsia. The presence of this gene in all eustigmatophytes studied, including representatives of both principal clades (Eustigmatales and Goniochloridales), is a genetic footprint indicating that the eustigmatophyte-Phycorickettsia partnership started no later than in the last eustigmatophyte common ancestor.

Advanced genetic tools enable synthetic biology in the oleaginous microalgae Nannochloropsis sp

Nannochloropsis is a genus of fast-growing microalgae that are regularly used for biotechnology applications. Nannochloropsis species have a high triacylglycerol content and their polar lipids are rich in the omega-3 long-chain polyunsaturated fatty acid, eicosapentaenoic acid. Placed in the heterokont lineage, the Nannochloropsis genus has a complex evolutionary history. Genome sequences are available for several species, and a number of transcriptomic datasets have been produced, making this genus a facile model for comparative genomics. There is a growing interest in Nannochloropsis species as models for the study of microalga lipid metabolism and as a chassis for synthetic biology. Recently, techniques for gene stacking, and targeted gene disruption and repression in the Nannochloropsis genus have been developed. These tools enable gene-specific, mechanistic studies and have already allowed the engineering of improved Nannochloropsis strains with superior growth, or greater bioproduction.

Optimized methods of chromatin immunoprecipitation for profiling histone modifications in industrial microalgae Nannochloropsis spp

Epigenetic factors such as histone modifications play integral roles in plant development and stress response, yet their implications in algae remain poorly understood. In the industrial oleaginous microalgae Nannochloropsis spp., the lack of an efficient methodology for chromatin immunoprecipitation (ChIP), which determines the specific genomic location of various histone modifications, has hindered probing the epigenetic basis of their photosynthetic carbon conversion and storage as oil. Here, a detailed ChIP protocol was developed for Nannochloropsis oceanica, which represents a reliable approach for the analysis of histone modifications, chromatin state, and transcription factor-binding sites at the epigenetic level. Using ChIP-qPCR, genes related to photosynthetic carbon fixation in this microalga were systematically assessed. Furthermore, a ChIP-Seq protocol was established and optimized, which generated a genome-wide profile of histone modification events, using histone mark H3K9Ac as an example. These results are the first step for appreciation of the chromatin landscape in industrial oleaginous microalgae and for epigenetics-based microalgal feedstock development.

Recurrent Loss, Horizontal Transfer, and the Obscure Origins of Mitochondrial Introns in Diatoms (Bacillariophyta)

We sequenced mitochondrial genomes from five diverse diatoms (Toxarium undulatum, Psammoneis japonica, Eunotia naegelii, Cylindrotheca closterium, and Nitzschia sp.), chosen to fill important phylogenetic gaps and help us characterize broadscale patterns of mitochondrial genome evolution in diatoms. Although gene content was strongly conserved, intron content varied widely across species. The vast majority of introns were of group II type and were located in the cox1 or rnl genes. Although recurrent intron loss appears to be the principal underlying cause of the sporadic distributions of mitochondrial introns across diatoms, phylogenetic analyses showed that intron distributions superficially consistent with a recurrent-loss model were sometimes more complicated, implicating horizontal transfer as a likely mechanism of intron acquisition as well. It was not clear, however, whether diatoms were the donors or recipients of horizontally transferred introns, highlighting a general challenge in resolving the evolutionary histories of many diatom mitochondrial introns. Although some of these histories may become clearer as more genomes are sampled, high rates of intron loss suggest that the origins of many diatom mitochondrial introns are likely to remain unclear.

Complete Genome Sequence of the Model Oleaginous Alga Nannochloropsis gaditana CCMP1894

The model oleaginous alga Nannochloropsis gaditana was completely sequenced using a combination of optical mapping and next-generation sequencing technologies to generate one of the most complete eukaryotic genomes published to date. The assembled genome is 30.7 Mb long.

Engineering the Chloroplast Genome of Oleaginous Marine Microalga Nannochloropsis oceanica

Plastid engineering offers an important tool to fill the gap between the technical and the enormous potential of microalgal photosynthetic cell factory. However, to date, few reports on plastid engineering in industrial microalgae have been documented. This is largely due to the small cell sizes and complex cell-wall structures which make these species intractable to current plastid transformation methods (i.e., biolistic transformation and polyethylene glycol-mediated transformation). Here, employing the industrial oleaginous microalga Nannochloropsis oceanica as a model, an electroporation-mediated chloroplast transformation approach was established. Fluorescent microscopy and laser confocal scanning microscopy confirmed the expression of the green fluorescence protein, driven by the endogenous plastid promoter and terminator. Zeocin-resistance selection led to an acquisition of homoplasmic strains of which a stable and site-specific recombination within the chloroplast genome was revealed by sequencing and DNA gel blotting. This demonstration of electroporation-mediated chloroplast transformation opens many doors for plastid genome editing in industrial microalgae, particularly species of which the chloroplasts are recalcitrant to chemical and microparticle bombardment transformation.

The plastid genome of some eustigmatophyte algae harbours a bacteria-derived six-gene cluster for biosynthesis of a novel secondary metabolite

Acquisition of genes by plastid genomes (plastomes) via horizontal gene transfer (HGT) seems to be a rare phenomenon. Here, we report an interesting case of HGT revealed by sequencing the plastomes of the eustigmatophyte algae Monodopsis sp. MarTras21 and Vischeria sp. CAUP Q 202. These plastomes proved to harbour a unique cluster of six genes, most probably acquired from a bacterium of the phylum Bacteroidetes, with homologues in various bacteria, typically organized in a conserved uncharacterized putative operon. Sequence analyses of the six proteins encoded by the operon yielded the following annotation for them: (i) a novel family without discernible homologues; (ii) a new family within the superfamily of metallo-dependent hydrolases; (iii) a novel subgroup of the UbiA superfamily of prenyl transferases; (iv) a new clade within the sugar phosphate cyclase superfamily; (v) a new family within the xylose isomerase-like superfamily; and (vi) a hydrolase for a phosphate moiety-containing substrate. We suggest that the operon encodes enzymes of a pathway synthesizing an isoprenoid–cyclitol-derived compound, possibly an antimicrobial or other protective substance. To the best of our knowledge, this is the first report of an expansion of the metabolic capacity of a plastid mediated by HGT into the plastid genome.

RNAi-based targeted gene knockdown in the model oleaginous microalgae Nannochloropsis oceanica

Microalgae are promising feedstock for renewable fuels such as biodiesel, yet development of industrial oleaginous strains has been hindered by the paucity and inefficiency of reverse genetics tools. Here we established an efficient RNAi-based targeted gene-knockdown method for Nannochloropsis spp., which are emerging model organisms for industrial microalgal oil production. The method achieved a 40-80% success rate in Nannochloropsis oceanica strain IMET1. When transcript level of one carbonic anhydrase (CA) was inhibited by 62-83% via RNAi, mutant cells exhibited photosynthetic oxygen evolution (POE) rates that were 68-100% higher than wild-type (WT) at pH 6.0, equivalent to WT at pH 8.2, yet 39-45% lower than WT at pH 9.0. Moreover, the mutant POE rates were negatively correlated with the increase of culture pH, an exact opposite of WT. Thus, a dynamic carbon concentration mechanism (CCM) that is highly sensitive to pH homeostasis was revealed, where the CA inhibition likely partially abrogated the mechanism that normally deactivates CCM under a high level of dissolved CO₂ . Extension of the method to another sequenced N. oceanica strain of CCMP 1779 demonstrated comparable performance. Finally, McrBC-PCR followed by bisulfite sequencing revealed that the gene knockdown is mediated by the CG, CHG and CHH types of DNA methylation at the coding region of the targeted gene. The efficiency, robustness and general applicability of this reverse genetics approach suggested the possibility of large-scale RNAi-based gene function screening in industrial microalgae.

Membrane Proteomic Insights into the Physiology and Taxonomy of an Oleaginous Green Microalga

Ettlia oleoabundans is a nonsequenced oleaginous green microalga. Despite the significant biotechnological interest in producing value-added compounds from the acyl lipids of this microalga, a basic understanding of the physiology and biochemistry of oleaginous microalgae is lacking, especially under nitrogen deprivation conditions known to trigger lipid accumulation. Using an RNA sequencing-based proteomics approach together with manual annotation, we are able to provide, to our knowledge, the first membrane proteome of an oleaginous microalga. This approach allowed the identification of novel proteins in E. oleoabundans, including two photoprotection-related proteins, Photosystem II Subunit S and Maintenance of Photosystem II under High Light1, which were considered exclusive to higher photosynthetic organisms, as well as Retinitis Pigmentosa Type 2-Clathrin Light Chain, a membrane protein with a novel domain architecture. Free-flow zonal electrophoresis of microalgal membranes coupled to liquid chromatography-tandem mass spectrometry proved to be a useful technique for determining the intracellular location of proteins of interest. Carbon-flow compartmentalization in E. oleoabundans was modeled using this information. Molecular phylogenetic analyses of protein markers and 18S ribosomal DNA support the reclassification of E. oleoabundans within the trebouxiophycean microalgae, rather than with the Chlorophyceae class, in which it is currently classified, indicating that it may not be closely related to the model green alga Chlamydomonas reinhardtii A detailed survey of biological processes taking place in the membranes of nitrogen-deprived E. oleoabundans, including lipid metabolism, provides insights into the basic biology of this nonmodel organism.

Genome editing of model oleaginous microalgae Nannochloropsis spp. by CRISPR/Cas9

Microalgae are promising feedstock for biofuels yet mechanistic probing of their cellular network and industrial strain development have been hindered by lack of genome-editing tools. Nannochloropsis spp. are emerging model microalgae for scalable oil production and carbon sequestration. Here we established a CRISPR/Cas9-based precise genome-editing approach for the industrial oleaginous microalga Nannochloropsis oceanica, using nitrate reductase (NR; g7988) as example. A new screening procedure that compares between restriction enzyme-digested nested PCR (nPCR) products derived from enzyme-digested and not-digested genomic DNA of transformant pools was developed to quickly, yet reliably, detect genome-engineered mutants. Deep sequencing of nPCR products directly amplified from pooled genomic DNA revealed over an 1% proportion of 5-bp deletion mutants and a lower frequency of 12-bp deletion mutants, with both types of editing precisely located at the targeted site. The isolated mutants, in which precise deletion of five bases caused a frameshift in NR translation, grow normally under NH₄ Cl but fail to grow under NaNO₃ , and thus represent a valuable chassis strain for transgenic-strain development. This demonstration of CRISPR/Cas9-based genome editing in industrial microalgae opens many doors for microalgae-based biotechnological applications.

Comparative genomics of mitochondria in chlorarachniophyte algae: endosymbiotic gene transfer and organellar genome dynamics

Chlorarachniophyte algae possess four DNA-containing compartments per cell, the nucleus, mitochondrion, plastid and nucleomorph, the latter being a relic nucleus derived from a secondary endosymbiont. While the evolutionary dynamics of plastid and nucleomorph genomes have been investigated, a comparative investigation of mitochondrial genomes (mtDNAs) has not been carried out. We have sequenced the complete mtDNA of Lotharella oceanica and compared it to that of another chlorarachniophyte, Bigelowiella natans. The linear mtDNA of L. oceanica is 36.7 kbp in size and contains 35 protein genes, three rRNAs and 24 tRNAs. The codons GUG and UUG appear to be capable of acting as initiation codons in the chlorarachniophyte mtDNAs, in addition to AUG. Rpl16, rps4 and atp8 genes are missing in L.oceanica mtDNA, despite being present in B. natans mtDNA. We searched for, and found, mitochondrial rpl16 and rps4 genes with spliceosomal introns in the L. oceanica nuclear genome, indicating that mitochondrion-to-host-nucleus gene transfer occurred after the divergence of these two genera. Despite being of similar size and coding capacity, the level of synteny between L. oceanica and B. natans mtDNA is low, suggesting frequent rearrangements. Overall, our results suggest that chlorarachniophyte mtDNAs are more evolutionarily dynamic than their plastid counterparts.

Stress-induced neutral lipid biosynthesis in microalgae - Molecular, cellular and physiological insights

Photosynthetic microalgae have promise as biofuel feedstock. Under certain conditions, they produce substantial amounts of neutral lipids, mainly in the form of triacylglycerols (TAGs), which can be converted to fuels. Much of our current knowledge on the genetic and molecular basis of algal neutral lipid metabolism derives mainly from studies of plants, i.e. seed tissues, and to a lesser extent from direct studies of algal lipid metabolism. Thus, the knowledge of TAG synthesis and the cellular trafficking of TAG precursors in algal cells is to a large extent based on genome predictions, and most aspects of TAG metabolism have yet to be experimentally verified. The biofuel prospects of microalgae have raised the interest in mechanistic studies of algal TAG biosynthesis in recent years and resulted in an increasing number of publications on lipid metabolism in microalgae. In this review we summarize the current findings on genetic, molecular and physiological studies of TAG accumulation in microalgae. Special emphasis is on the functional analysis of key genes involved in TAG synthesis, molecular mechanisms of regulation of TAG biosynthesis, as well as on possible mechanisms of lipid droplet formation in microalgal cells. This article is part of a Special Issue entitled: Plant Lipid Biology edited by Kent D. Chapman and Ivo Feussner.

A Comparative Analysis of Mitochondrial Genomes in Eustigmatophyte Algae

Eustigmatophyceae (Ochrophyta, Stramenopiles) is a small algal group with species of the genus Nannochloropsis being its best studied representatives. Nuclear and organellar genomes have been recently sequenced for several Nannochloropsis spp., but phylogenetically wider genomic studies are missing for eustigmatophytes. We sequenced mitochondrial genomes (mitogenomes) of three species representing most major eustigmatophyte lineages, Monodopsis sp. MarTras21, Vischeria sp. CAUP Q 202 and Trachydiscus minutus, and carried out their comparative analysis in the context of available data from Nannochloropsis and other stramenopiles, revealing a number of noticeable findings. First, mitogenomes of most eustigmatophytes are highly collinear and similar in the gene content, but extensive rearrangements and loss of three otherwise ubiquitous genes happened in the Vischeria lineage; this correlates with an accelerated evolution of mitochondrial gene sequences in this lineage. Second, eustigmatophytes appear to be the only ochrophyte group with the Atp1 protein encoded by the mitogenome. Third, eustigmatophyte mitogenomes uniquely share a truncated nad11 gene encoding only the C-terminal part of the Nad11 protein, while the N-terminal part is encoded by a separate gene in the nuclear genome. Fourth, UGA as a termination codon and the cognate release factor mRF2 were lost from mitochondria independently by the Nannochloropsis and T. minutus lineages. Finally, the rps3 gene in the mitogenome of Vischeria sp. is interrupted by the UAG codon, but the genome includes a gene for an unusual tRNA with an extended anticodon loop that we speculate may serve as a suppressor tRNA to properly decode the rps3 gene.

Algorithms for reconstruction of chromosomal structures

BACKGROUND

One of the main aims of phylogenomics is the reconstruction of objects defined in the leaves along the whole phylogenetic tree to minimize the specified functional, which may also include the phylogenetic tree generation. Such objects can include nucleotide and amino acid sequences, chromosomal structures, etc. The structures can have any set of linear and circular chromosomes, variable gene composition and include any number of paralogs, as well as any weights of individual evolutionary operations to transform a chromosome structure. Many heuristic algorithms were proposed for this purpose, but there are just a few exact algorithms with low (linear, cubic or similar) polynomial computational complexity among them to our knowledge. The algorithms naturally start from the calculation of both the distance between two structures and the shortest sequence of operations transforming one structure into another. Such calculation per se is an NP-hard problem.

RESULTS

A general model of chromosomal structure rearrangements is considered. Exact algorithms with almost linear or cubic polynomial complexities have been developed to solve the problems for the case of any chromosomal structure but with certain limitations on operation weights. The computer programs are tested on biological data for the problem of mitochondrial or plastid chromosomal structure reconstruction. To our knowledge, no computer programs are available for this model.

CONCLUSIONS

Exactness of the proposed algorithms and such low polynomial complexities were proved. The reconstructed evolutionary trees of mitochondrial and plastid chromosomal structures as well as the ancestral states of the structures appear to be reasonable.

Transcriptional coordination of physiological responses in Nannochloropsis oceanica CCMP1779 under light/dark cycles

Nannochloropsis oceanica CCMP1779 is a marine unicellular stramenopile and an emerging reference species for basic research on oleogenic microalgae with biotechnological relevance. We investigated its physiology and transcriptome under light/dark cycles. We observed oscillations in lipid content and a predominance of cell division in the first half of the dark phase. Globally, more than 60% of the genes cycled in N. oceanica CCMP1779, with gene expression peaking at different times of the day. Interestingly, the phase of expression of genes involved in certain biological processes was conserved across photosynthetic lineages. Furthermore, in agreement with our physiological studies we found the processes of lipid metabolism and cell division enriched in cycling genes. For example, there was tight coordination of genes involved in the lower part of glycolysis, fatty acid synthesis and lipid production at dawn preceding lipid accumulation during the day. Our results suggest that diel lipid storage plays a key role for N. oceanica CCMP1779 growth under natural conditions making this alga a promising model to gain a basic mechanistic understanding of triacylglycerol production in photosynthetic cells. Our data will help the formulation of new hypotheses on the role of cyclic gene expression in cell growth and metabolism in Nannochloropsis.

Rapid Characterization of Microalgae and Microalgae Mixtures Using Matrix-Assisted Laser Desorption Ionization Time-Of-Flight Mass Spectrometry (MALDI-TOF MS)

Current molecular methods to characterize microalgae are time-intensive and expensive. Matrix Assisted Laser Desorption/Ionization Time-of-Flight Mass Spectrometry (MALDI-TOF MS) may represent a rapid and economical alternative approach. The objectives of this study were to determine whether MALDI-TOF MS can be used to: 1) differentiate microalgae at the species and strain levels and 2) characterize simple microalgal mixtures. A common protein extraction sample preparation method was used to facilitate rapid mass spectrometry-based analysis of 31 microalgae. Each yielded spectra containing between 6 and 56 peaks in the m/z 2,000 to 20,000 range. The taxonomic resolution of this approach appeared higher than that of 18S rDNA sequence analysis. For example, two strains of Scenedesmus acutus differed only by two 18S rDNA nucleotides, but yielded distinct MALDI-TOF mass spectra. Mixtures of two and three microalgae yielded relatively complex spectra that contained peaks associated with members of each mixture. Interestingly, though, mixture-specific peaks were observed at m/z 11,048 and 11,230. Our results suggest that MALDI-TOF MS affords rapid characterization of individual microalgae and simple microalgal mixtures.

Guanchochroma wildpretii gen. et spec. nov. (Ochrophyta) Provides New Insights into the Diversification and Evolution of the Algal Class Synchromophyceae

A new relative of the chrysophyte genus Chrysopodocystis was found in Tenerife and termed Guanchochroma wildpretii. This unicellular alga was most noticeably discernible from Chrysopodocystis socialis (the only species of this genus) by the presence of a cyst-like stage with a multilayered lorica, which also functions as a dispersal unit and shows secondary wall growth. Secondary expansion of loricae (cell casings not involved in cell division, usually with a more or less pronounced opening) has never been observed previously and marks a unique feature of the new taxon. Plastids are non-randomly distributed within cells of G. wildpretii. 18S rRNA gene analyses identified the two species as sister lineages and placed them in a monophyletic group with the Synchromophyceae, a heterokont algal (Ochrophyta) class characterized by the presence of chloroplast complexes. Yet, neither Chrysopodocystis nor Guanchochroma showed this feature in ultrastructure analyses. Additionally, their 18S rRNA genes possessed distinct inserts, the highest GC-content known for Ochrophyta and exceptionally long branches on the Ochrophyta 18S rDNA phylogenetic tree, suggesting substantially increased substitution rates along their branch compared to Synchromophyceae. Plastid marker data (rbcL) recovered a monophyletic clade of Chrysopodocystis, Guanchochroma and Synchromophyceae as well, yet with lower supports for internal split order due to limited resolution of the marker. Evidence for the sequence of events leading to the formation of the plastid complex of Synchromophyceae still remains ambiguous because of the apparently short timeframe in which they occurred.

Antagonistic roles of abscisic acid and cytokinin during response to nitrogen depletion in oleaginous microalga Nannochloropsis oceanica expand the evolutionary breadth of phytohormone function

The origin of phytohormones is poorly understood, and their physiological roles in microalgae remain elusive. Genome comparison of photosynthetic autotrophic eukaryotes has revealed that the biosynthetic pathways of abscisic acid (ABA) and cytokinins (CKs) emerged in unicellular algae. While ABA and CK degradation mechanisms emerged broadly in algal lineages, complete vascular plant-type conjugation pathways emerged prior to the rise of Streptophyta. In microalgae, a complete set of proteins from the canonical ABA and CK sensing and signaling pathways is not essential, but individual components are present, suggesting stepwise recruitment of phytohormone signaling components. In the oleaginous eustigmatophyte Nannochloropsis oceanica IMET1, UHPLC-MS/MS detected a wide array of plant hormones, despite a phytohormone profile that is very distinct from that of flowering plants. Time-series transcriptional analysis during nitrogen depletion revealed activation of the ABA biosynthetic pathway and antagonistic transcription of CK biosynthetic genes. Correspondingly, the ABA level increases while the dominant bioactive CK forms decrease. Moreover, exogenous CKs stimulate cell-cycle progression while exogenous ABA acts as both an algal growth repressor and a positive regulator in response to stresses. The presence of such functional flowering plant-like phytohormone signaling systems in Nannochloropsis sp. suggests a much earlier origin of phytohormone biosynthesis and degradation than previously believed, and supports the presence in microalgae of as yet unknown conjugation and sensing/signaling systems that may be exploited for microalgal feedstock development.

Genome-wide identification of transcription factors and transcription-factor binding sites in oleaginous microalgae Nannochloropsis

Nannochloropsis spp. are a group of oleaginous microalgae that harbor an expanded array of lipid-synthesis related genes, yet how they are transcriptionally regulated remains unknown. Here a phylogenomic approach was employed to identify and functionally annotate the transcriptional factors (TFs) and TF binding-sites (TFBSs) in N. oceanica IMET1. Among 36 microalgae and higher plants genomes, a two-fold reduction in the number of TF families plus a seven-fold decrease of average family-size in Nannochloropsis, Rhodophyta and Chlorophyta were observed. The degree of similarity in TF-family profiles is indicative of the phylogenetic relationship among the species, suggesting co-evolution of TF-family profiles and species. Furthermore, comparative analysis of six Nannochloropsis genomes revealed 68 “most-conserved” TFBS motifs, with 11 of which predicted to be related to lipid accumulation or photosynthesis. Mapping the IMET1 TFs and TFBS motifs to the reference plant TF-“TFBS motif” relationships in TRANSFAC enabled the prediction of 78 TF-“TFBS motif” interaction pairs, which consisted of 34 TFs (with 11 TFs potentially involved in the TAG biosynthesis pathway), 30 TFBS motifs and 2,368 regulatory connections between TFs and target genes. Our results form the basis of further experiments to validate and engineer the regulatory network of Nannochloropsis spp. for enhanced biofuel production.

Choreography of Transcriptomes and Lipidomes of Nannochloropsis Reveals the Mechanisms of Oil Synthesis in Microalgae

To reveal the molecular mechanisms of oleaginousness in microalgae, transcriptomic and lipidomic dynamics of the oleaginous microalga Nannochloropsis oceanica IMET1 under nitrogen-replete (N+) and N-depleted (N-) conditions were simultaneously tracked. At the transcript level, enhanced triacylglycerol (TAG) synthesis under N- conditions primarily involved upregulation of seven putative diacylglycerol acyltransferase (DGAT) genes and downregulation of six other DGAT genes, with a simultaneous elevation of the other Kennedy pathway genes. Under N- conditions, despite downregulation of most de novo fatty acid synthesis genes, the pathways that shunt carbon precursors from protein and carbohydrate metabolic pathways into glycerolipid synthesis were stimulated at the transcript level. In particular, the genes involved in supplying carbon precursors and energy for de novo fatty acid synthesis, including those encoding components of the pyruvate dehydrogenase complex (PDHC), glycolysis, and PDHC bypass, and suites of specific transporters, were substantially upregulated under N- conditions, resulting in increased overall TAG production. Moreover, genes involved in the citric acid cycle and β-oxidation in mitochondria were greatly enhanced to utilize the carbon skeletons derived from membrane lipids and proteins to produce additional TAG or its precursors. This temporal and spatial regulation model of oil accumulation in microalgae provides a basis for improving our understanding of TAG synthesis in microalgae and will also enable more rational genetic engineering of TAG production.

A pangenomic analysis of the Nannochloropsis organellar genomes reveals novel genetic variations in key metabolic genes

BACKGROUND

Microalgae in the genus Nannochloropsis are photosynthetic marine Eustigmatophytes of significant interest to the bioenergy and aquaculture sectors due to their ability to efficiently accumulate biomass and lipids for utilization in renewable transportation fuels, aquaculture feed, and other useful bioproducts. To better understand the genetic complement that drives the metabolic processes of these organisms, we present the assembly and comparative pangenomic analysis of the chloroplast and mitochondrial genomes from Nannochloropsis salina CCMP1776.

RESULTS

The chloroplast and mitochondrial genomes of N. salina are 98.4% and 97% identical to their counterparts in Nannochloropsis gaditana. Comparison of the Nannochloropsis pangenome to other algae within and outside of the same phyla revealed regions of significant genetic divergence in key genes that encode proteins needed for regulation of branched chain amino synthesis (acetohydroxyacid synthase), carbon fixation (RuBisCO activase), energy conservation (ATP synthase), protein synthesis and homeostasis (Clp protease, ribosome).

CONCLUSIONS

Many organellar gene modifications in Nannochloropsis are unique and deviate from conserved orthologs found across the tree of life. Implementation of secondary and tertiary structure prediction was crucial to functionally characterize many proteins and therefore should be implemented in automated annotation pipelines. The exceptional similarity of the N. salina and N. gaditana organellar genomes suggests that N. gaditana be reclassified as a strain of N. salina.

Regulation of the cholesterol biosynthetic pathway and its integration with fatty acid biosynthesis in the oleaginous microalga Nannochloropsis oceanica

BACKGROUND

Sterols are vital structural and regulatory components in eukaryotic cells; however, their biosynthetic pathways and functional roles in microalgae remain poorly understood.

RESULTS

In the oleaginous microalga Nannochloropsis oceanica, the sterol biosynthetic pathway produces phytosterols as minor products and cholesterol as the major product. The evidence together with their deduced biosynthetic pathways suggests that N. oceanica exhibits features of both higher plants and mammals. Temporal tracking of sterol profiles and sterol-biosynthetic transcripts in response to changes in light intensity and nitrogen supply reveal that sterols play roles in cell proliferation, chloroplast differentiation, and photosynthesis. Furthermore, the dynamics of fatty acid (FA) and FA-biosynthetic transcripts upon chemical inhibitor-induced sterol depletion reveal possible co-regulation of sterol production and FA synthesis, in that the squalene epoxidase inhibitor terbinafine reduces sterol content yet significantly elevates free FA production. Thus, a feedback regulation of sterol and FA homeostasis is proposed, with the 1-deoxy-D-xylulose 5-phosphate synthase (DXS, the committed enzyme in isoprenoid and sterol biosynthesis) gene potentially subject to feedback regulation by sterols.

CONCLUSION

These findings reveal features of sterol function and biosynthesis in microalgae and suggest new genetic engineering or chemical biology approaches for enhanced oil production in microalgae.