Morphological and molecular phylogenetic analysis of the high triglyceride‐producing marine diatom,Fistulifera solaris sp. nov. (Bacillariophyceae)

Fistulifera sp. strain JPCC DA0580, a marine pennate diatom, contains extremely high levels of intracellular triglyceride and has been suggested as a promising source of feedstock for biodiesel fuels. JPCC DA0580 was isolated from a mangrove swamp located in Sumiyo Bay below the mouths of the Sumiyo and Yakugachi Rivers in Amami‐Ohshima, Kagoshima, Japan. In this study, the morphology and the 18S rDNA sequence of JPCC DA0580 were compared with those of other Fistulifera strains. JPCC DA0580 possesses morphological characters of the genus Fistulifera, namely lightly silicified frustules, a distinct median costa (raphe sternum), and a wart‐like central pore (fistula). Morphometric analysis revealed that JPCC DA0580 differs from other Fistulifera species by the presence of a valve with coarser striation and coarser areolation. On the basis of 18S rDNA phylogeny, JPCC DA0580 formed a well‐supported clade with other members of the Fistulifera species complex, although the number of nucleotide substitutions was highest in JPCC DA0580. Our results led us to propose the taxonomic name Fistulifera solaris sp. nov. for JPCC DA0580.     

Establishment of a Genetic Transformation System for the Marine Pennate Diatom Fistulifera sp. Strain JPCC DA0580—A High Triglyceride Producer

A genetic transformation system for the marine pennate diatom, Fistulifera sp. JPCC DA0580, was established using microparticle bombardment methods. Strain JPCC DA0580 has been recently identified as the highest triglyceride (60 % w/w) producer from a culture collection of 1,393 strains of marine microalgae, and it is expected to be a feasible source of biodiesel fuel. The transformation conditions for strain JPCC DA0580 were optimised using the green fluorescent protein gene (gfp) and the gene encoding neomycin phosphotransferase II (nptII). The most efficient rate of transformation was attained when tungsten particles (0.6 μm in diameter) were used for microparticle bombardment. The effect of endogenous and exogenous promoters on the expression of nptII was examined. Endogenous promoters were more efficient for obtaining transformants compared with exogenous promoters. Southern hybridisation analysis suggested that nptII integrated into the nuclear genome. This genetic manipulation technique should allow us to understand the mechanisms of high triglyceride accumulation in this strain, thereby contributing to improving BDF production.     

Profiling of fatty acid methyl esters from the oleaginous diatom Fistulifera sp. strain JPCC DA0580 under nutrition-sufficient and -deficient conditions

The marine oleaginous diatom, Fistulifera sp. strain JPCC DA0580, is a promising candidate for biodiesel production due to its high lipid content. In order to truly evaluate the potential of this strain as biodiesel feedstock as well as the impact of nutrition-deficiency to this strain, the proportion of the lipid fractions and fatty acid methyl ester (FAME) derived from Fistulifera sp. cultured under nutrition-sufficient or -deficient conditions were analyzed. The nutrition deficiency led to the increase of the total lipid content in the form of neutral lipids (NLs) accumulation and the decline of polar lipids compared with nutrition-sufficiency. Meanwhile, the total lipid productivity was not significantly changed under two nutrition conditions while the NL productivity under nutrition-deficient condition was much higher than nutrition-sufficient condition. The major FAME components, C14:0, C16:0, C16:1, and C20:5, contribute to over 90 % of total FAMEs under both nutrition conditions. A lower polyunsaturated FAME level were observed in the nutrition-deficient condition (9.9?±?0.2 %) compared with the nutrition-sufficient condition (19.8?±?1.2 %), suggesting the availability of the nutrition stress on the strain JPCC DA0580 for improvement of fuel quality as well as productivity. The lipid quality estimation based on the FAME profile revealed that the nutrition-deficiency could further improve the lipid quality of both total lipids and NL fraction. In addition, direct infusion ESI-Q-TRAP-MS/MS was carried out for the fractionated NL in order to estimate triacylglycerol (TAG) composition, suggesting a crucial role of the chloroplast in TAG synthesis.     

High-throughput pyrosequencing of the chloroplast genome of a highly neutral-lipid-producing marine pennate diatom, Fistulifera sp. strain JPCC DA0580

The chloroplast genome of the highly neutral-lipid-producing marine pennate diatom Fistulifera sp. strain JPCC DA0580 was fully sequenced using high-throughput pyrosequencing. The general features and gene content were compared with three other complete diatom chloroplast genomes. The chloroplast genome is 134,918 bp with an inverted repeat of 13,330 bp and is slightly larger than the other diatom chloroplast genomes due to several low gene-density regions lacking similarity to the other diatom chloroplast genomes. Protein-coding genes were nearly identical to those from Phaeodactylum tricornutum. On the other hand, we found unique sequence variations in genes of photosystem II which differ from the consensus in other diatom chloroplasts. Furthermore, five functional unknown ORFs and a putative serine recombinase gene, serC2, are located in the low gene-density regions. SerC2 was also identified in the plasmids of another pennate diatom, Cylindrotheca fusiformis, and in the plastid genome of the diatom endosymbiont of Kryptoperidinium foliaceum. Exogenous plasmids might have been incorporated into the chloroplast genome of Fistulifera sp. by lateral gene transfer. Chloroplast genome sequencing analysis of this novel diatom provides many important insights into diatom evolution.     

Reconstruction and In Silico Analysis of Metabolic Network for an Oleaginous Yeast,Yarrowia lipolytica

With the emergence of energy scarcity, the use of renewable energy sources such as biodiesel is becoming increasingly necessary. Recently, many researchers have focused their minds on Yarrowia lipolytica, a model oleaginous yeast, which can be employed to accumulate large amounts of lipids that could be further converted to biodiesel. In order to understand the metabolic characteristics of Y. lipolytica at a systems level and to examine the potential for enhanced lipid production, a genome-scale compartmentalized metabolic network was reconstructed based on a combination of genome annotation and the detailed biochemical knowledge from multiple databases such as KEGG, ENZYME and BIGG. The information about protein and reaction associations of all the organisms in KEGG and Expasy-ENZYME database was arranged into an EXCEL file that can then be regarded as a new useful database to generate other reconstructions. The generated model iYL619_PCP accounts for 619 genes, 843 metabolites and 1,142 reactions including 236 transport reactions, 125 exchange reactions and 13 spontaneous reactions. The in silico model successfully predicted the minimal media and the growing abilities on different substrates. With flux balance analysis, single gene knockouts were also simulated to predict the essential genes and partially essential genes. In addition, flux variability analysis was applied to design new mutant strains that will redirect fluxes through the network and may enhance the production of lipid. This genome-scale metabolic model of Y. lipolytica can facilitate system-level metabolic analysis as well as strain development for improving the production of biodiesels and other valuable products by Y. lipolytica and other closely related oleaginous yeasts.     

Gene Set Enrichment Analysis (GSEA) of Toxoplasma gondii expression datasets links cell cycle progression and the bradyzoite developmental program

Background

Large amounts of microarray expression data have been generated for the Apicomplexan parasite Toxoplasma gondii in an effort to identify genes critical for virulence or developmental transitions. However, researchers’ ability to analyze this data is limited by the large number of unannotated genes, including many that appear to be conserved hypothetical proteins restricted to Apicomplexa. Further, differential expression of individual genes is not always informative and often relies on investigators to draw big-picture inferences without the benefit of context. We hypothesized that customization of gene set enrichment analysis (GSEA) to T. gondii would enable us to rigorously test whether groups of genes serving a common biological function are co-regulated during the developmental transition to the latent bradyzoite form.

Results

Using publicly available T. gondii expression microarray data, we created Toxoplasma gene sets related to bradyzoite differentiation, oocyst sporulation, and the cell cycle. We supplemented these with lists of genes derived from community annotation efforts that identified contents of the parasite-specific organelles, rhoptries, micronemes, dense granules, and the apicoplast. Finally, we created gene sets based on metabolic pathways annotated in the KEGG database and Gene Ontology terms associated with gene annotations available at http://www.toxodb.org. These gene sets were used to perform GSEA analysis using two sets of published T. gondii expression data that characterized T. gondii stress response and differentiation to the latent bradyzoite form.

Conclusions

GSEA provides evidence that cell cycle regulation and bradyzoite differentiation are coupled. Δgcn5A mutants unable to induce bradyzoite-associated genes in response to alkaline stress have different patterns of cell cycle and bradyzoite gene expression from stressed wild-type parasites. Extracellular tachyzoites resemble a transitional state that differs in gene expression from both replicating intracellular tachyzoites and in vitro bradyzoites by expressing genes that are enriched in bradyzoites as well as genes that are associated with the G1 phase of the cell cycle. The gene sets we have created are readily modified to reflect ongoing research and will aid researchers’ ability to use a knowledge-based approach to data analysis facilitating the development of new insights into the intricate biology of Toxoplasma gondii.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-515) contains supplementary material, which is available to authorized users.     

Soluble Methane Monooxygenase Gene Clusters from Trichloroethylene-Degrading Methylomonas sp. Strains and Detection of Methanotrophs during In Situ Bioremediation

The soluble MMO (sMMO) gene clusters from group I methanotrophs were characterized. An 8.1-kb KpnI fragment from Methylomonas sp. strain KSWIII and a 7.5-kb SalI fragment from Methylomonas sp. strain KSPIII which contained the sMMO gene clusters were cloned and sequenced. The sequences of these two fragments were almost identical. The sMMO gene clusters in the fragment consisted of six open reading frames which were 52 to 79% similar to the corresponding genes of previously described sMMO gene clusters of the group II and group X methanotrophs. The phylogenetic analysis of the predicted amino acid sequences of sMMO demonstrated that the sMMOs from these strains were closer to that from M. capsulatus Bath in the group X methanotrophs than to those from Methylosinus trichosporium OB3b and Methylocystis sp. strain M in the group II methanotrophs. Based on the sequence data of sMMO genes of our strains and other methanotrophs, we designed a new PCR primer to amplify sMMO gene fragments of all the known methanotrophs harboring the mmoX gene. The primer set was successfully used for detecting methanotrophs in the groundwater of trichloroethylene-contaminated sites during in situ-biostimulation treatments.     

PAGE: Parametric Analysis of Gene Set Enrichment

Background  

Gene set enrichment analysis (GSEA) is a microarray data analysis method that uses predefined gene sets and ranks of genes to identify significant biological changes in microarray data sets. GSEA is especially useful when gene expression changes in a given microarray data set is minimal or moderate.     

Construction of expression vectors for metabolic engineering of the vanillin-producing actinomycete Amycolatopsis sp. ATCC 39116

Amycolatopsis sp. ATCC 39116 is able to synthesize the important flavoring agent vanillin from cheap natural substrates. The bacterium is therefore of great interest for the industry and used for the fermentative production of vanillin. In order to improve the production of natural vanillin with Amycolatopsis sp. ATCC 39116, the strain has been genetically engineered to optimize the metabolic flux towards the desired product. Extensive metabolic engineering was hitherto hampered, due to the lack of genetic tools like functional promoters and expression vectors. In this study, we report the establishment of a plasmid-based gene expression system for Amycolatopsis sp. ATCC 39116 that allows a further manipulation of the genotype. Four new Escherichia coli–Amycolatopsis shuttle vectors harboring different promoter elements were constructed, and the functionality of these regulatory elements was proven by the expression of the reporter gene gusA, encoding a β-glucuronidase. Glucuronidase activity was detected in all plasmid-harboring strains, and remarkable differences in the expression strength of the reporter gene depending on the used promoter were observed. The new expression vectors will promote the further genetic engineering of Amycolatopsis sp. ATCC 39116 to get insight into the metabolic network and to improve the strain for a more efficient industrial use.     

Draft Genome Assembly and Annotation for Cutaneotrichosporon dermatis NICC30027, an Oleaginous Yeast Capable of Simultaneous Glucose and Xylose Assimilation

The identification of oleaginous yeast species capable of simultaneously utilizing xylose and glucose as substrates to generate value-added biological products is an area of key economic interest. We have previously demonstrated that the Cutaneotrichosporon dermatis NICC30027 yeast strain is capable of simultaneously assimilating both xylose and glucose, resulting in considerable lipid accumulation. However, as no high-quality genome sequencing data or associated annotations for this strain are available at present, it remains challenging to study the metabolic mechanisms underlying this phenotype. Herein, we report a 39,305,439 bp draft genome assembly for C. dermatis NICC30027 comprised of 37 scaffolds, with 60.15% GC content. Within this genome, we identified 524 tRNAs, 142 sRNAs, 53 miRNAs, 28 snRNAs, and eight rRNA clusters. Moreover, repeat sequences totaling 1,032,129 bp in length were identified (2.63% of the genome), as were 14,238 unigenes that were 1,789.35 bp in length on average (64.82% of the genome). The NCBI non-redundant protein sequences (NR) database was employed to successfully annotate 11,795 of these unigenes, while 3,621 and 11,902 were annotated with the Swiss-Prot and TrEMBL databases, respectively. Unigenes were additionally subjected to pathway enrichment analyses using the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Cluster of Orthologous Groups of proteins (COG), Clusters of orthologous groups for eukaryotic complete genomes (KOG), and Non-supervised Orthologous Groups (eggNOG) databases. Together, these results provide a foundation for future studies aimed at clarifying the mechanistic basis for the ability of C. dermatis NICC30027 to simultaneously utilize glucose and xylose to synthesize lipids.     

Computation of significance scores of unweighted Gene Set Enrichment Analyses

Background  

Gene Set Enrichment Analysis (GSEA) is a computational method for the statistical evaluation of sorted lists of genes or proteins. Originally GSEA was developed for interpreting microarray gene expression data, but it can be applied to any sorted list of genes. Given the gene list and an arbitrary biological category, GSEA evaluates whether the genes of the considered category are randomly distributed or accumulated on top or bottom of the list. Usually, significance scores (p-values) of GSEA are computed by nonparametric permutation tests, a time consuming procedure that yields only estimates of the p-values.     

Genes Involved in Degradation of para-Nitrophenol Are Differentially Arranged in Form of Non-Contiguous Gene Clusters in Burkholderia sp. strain SJ98

Biodegradation of para-Nitrophenol (PNP) proceeds via two distinct pathways, having 1,2,3-benzenetriol (BT) and hydroquinone (HQ) as their respective terminal aromatic intermediates. Genes involved in these pathways have already been studied in different PNP degrading bacteria. Burkholderia sp. strain SJ98 degrades PNP via both the pathways. Earlier, we have sequenced and analyzed a ~41 kb fragment from the genomic library of strain SJ98. This DNA fragment was found to harbor all the lower pathway genes; however, genes responsible for the initial transformation of PNP could not be identified within this fragment. Now, we have sequenced and annotated the whole genome of strain SJ98 and found two ORFs (viz., pnpA and pnpB) showing maximum identity at amino acid level with p-nitrophenol 4-monooxygenase (PnpM) and p-benzoquinone reductase (BqR). Unlike the other PNP gene clusters reported earlier in different bacteria, these two ORFs in SJ98 genome are physically separated from the other genes of PNP degradation pathway. In order to ascertain the identity of ORFs pnpA and pnpB, we have performed in-vitro assays using recombinant proteins heterologously expressed and purified to homogeneity. Purified PnpA was found to be a functional PnpM and transformed PNP into benzoquinone (BQ), while PnpB was found to be a functional BqR which catalyzed the transformation of BQ into hydroquinone (HQ). Noticeably, PnpM from strain SJ98 could also transform a number of PNP analogues. Based on the above observations, we propose that the genes for PNP degradation in strain SJ98 are arranged differentially in form of non-contiguous gene clusters. This is the first report for such arrangement for gene clusters involved in PNP degradation. Therefore, we propose that PNP degradation in strain SJ98 could be an important model system for further studies on differential evolution of PNP degradation functions.     

Engineering Aspergillus oryzae A-4 through the Chromosomal Insertion of Foreign Cellulase Expression Cassette to Improve Conversion of Cellulosic Biomass into Lipids

A genetic modification scheme was designed for Aspergillus oryzae A-4, a natural cellulosic lipids producer, to enhance its lipid production from biomass by putting the spotlight on improving cellulase secretion. Four cellulase genes were separately expressed in A-4 under the control of hlyA promoter, with the help of the successful development of a chromosomal genetic manipulation system. Comparison of cellulase activities of PCR-positive transformants showed that these transformants integrated with celA gene and with celC gene had significantly (p<0.05) higher average FPAase activities than those strains integrated with celB gene and with celD gene. Through the assessment of cellulosic lipids accumulating abilities, celA transformant A2-2 and celC transformant D1-B1 were isolated as promising candidates, which could yield 101%–133% and 35.22%–59.57% higher amount of lipids than the reference strain A-4 (WT) under submerged (SmF) conditions and solid-state (SSF) conditions, respectively. Variability in metabolism associated to the introduction of cellulase gene in A2-2 and D1-B1 was subsequently investigated. It was noted that cellulase expression repressed biomass formation but enhanced lipid accumulation; whereas the inhibitory effect on cell growth would be shielded during cellulosic lipids production owing to the essential role of cellulase in substrate utilization. Different metabolic profiles also existed between A2-2 and D1-B1, which could be attributed to not only different transgene but also biological impacts of different integration. Overall, both simultaneous saccharification and lipid accumulation were enhanced in A2-2 and D1-B1, resulting in efficient conversion of cellulose into lipids. A regulation of cellulase secretion in natural cellulosic lipids producers could be a possible strategy to enhance its lipid production from lignocellulosic biomass.     

Regulation of Gene Expression in Diverse Cyanobacterial Species by Using Theophylline-Responsive Riboswitches

Cyanobacteria are photosynthetic bacteria that are currently being developed as biological production platforms. They derive energy from light and carbon from atmospheric carbon dioxide, and some species can fix atmospheric nitrogen. One advantage of developing cyanobacteria for renewable production of biofuels and other biological products is that they are amenable to genetic manipulation, facilitating bioengineering and synthetic biology. To expand the currently available genetic toolkit, we have demonstrated the utility of synthetic theophylline-responsive riboswitches for effective regulation of gene expression in four diverse species of cyanobacteria, including two recent isolates. We evaluated a set of six riboswitches driving the expression of a yellow fluorescent protein reporter in Synechococcus elongatus PCC 7942, Leptolyngbya sp. strain BL0902, Anabaena sp. strain PCC 7120, and Synechocystis sp. strain WHSyn. We demonstrated that riboswitches can offer regulation of gene expression superior to that of the commonly used isopropyl-β-d-thiogalactopyranoside induction of a lacI^q-P_trc promoter system. We also showed that expression of the toxic protein SacB can be effectively regulated, demonstrating utility for riboswitch regulation of proteins that are detrimental to biomass accumulation. Taken together, the results of this work demonstrate the utility and ease of use of riboswitches in the context of genetic engineering and synthetic biology in diverse cyanobacteria, which will facilitate the development of algal biotechnology.