A pangenomic analysis of the Nannochloropsisorganellar genomes reveals novel genetic variations in key metabolic genes
© Starkenburg et al.; licensee BioMed Central Ltd. 2014
Received: 24 April 2013
Accepted: 11 March 2014
Published: 19 March 2014
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.
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).
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.
KeywordsNannochloropsis Chloroplast Mitochondria Genome Stramenopiles Genome evolution Gene divergence
Stramenopiles encompass a broad array of golden brown algae that are morphologically diverse, ranging from unicells (e.g., diatoms) to large bladed species (e.g., kelps). These organisms acquired their chloroplast via secondary endosymbiosis, thus their evolutionary progression differs significantly from that of their green (Chorophyta) and red (Rhodophyta) primary endosymbiotic algal counterparts . Among the 17 classes of stramenopiles, the Eustigmatophyceae represent one of the smallest divisions. Members of this class, found in fresh, brackish, and marine waters, are minute in size, coccoid in shape, yellow-green in color, and essentially indistinguishable from one another given the lack of defining morphological characteristics .
Select unicellular photosynthetic microalgae have been targeted for commercial applications given their ability to efficiently accumulate biomass and/or lipids for conversion into renewable transportation fuels and other useful bioproducts. Algae within the Eustigmatophyceae, specifically within the genus Nannochloropsis are actively being evaluated for use in biofuel and aquaculture production systems due to their ability to convert a significant portion of their biomass (up to 60% dry weight) into lipids [3–5]. Although significant effort has been expended to characterize growth phenotypes and the fatty acid content within the genus Nannochloropsis [6–11], knowledge of the genetic and genomic basis that defines and controls their physiological behavior are still lacking; critical information required to support effective genetic engineering strategies. Recently, an analysis of the mitochondrial and chloroplast genomes of seven strains from six species of Nannochloropsis revealed that the genomic content was highly conserved between these species yet, evolutionarily divergent ‘hotspots’ were present, enabling an accurate phylotyping of these closely related species .
Here, we present the first analysis of the chloroplast and mitochondrial genomes from N. salina CCMP1776 and the resequencing and analysis of N. oculata CCMP525. To determine the unique features of these Nannochloropsis organelles, we compared these genomes to the complete organellar genomes of Nannochloropsis gaditana CCMP526 , an improved draft assembly and annotation of Nannochloropsis oceanica str. LAMB0001 , and to the six strains analyzed by Wei, et al. . Through these analyses, genomic variations and similarities were identified between Nannochloropsis and its stramenopile relatives. Striking similarity was observed between the organellar genomes of N. salina and N. gaditina. Additionally, novel modifications to key metabolic genes in the organelles of the genus Nannochloropsis were uncovered which further inform the physiological properties of this unique algal taxon.
Culturing and DNA purification
Nannochloropsis salina (CCMP1776) was grown at 30°C with a modified F/2-Si media with 10X nitrate and 7X phosphate  utilizing fluorescent plant grow lights at 1200 μEm-2 s-1 on a 16/ 8 hour light dark cycle. Dissolved O2 was maintained at 100% of base level from an 80% N2/ 20% O2 atmosphere through mass flow regulation of N2 or O2 gas input. Cell cultures were maintained at pH 8.2 utilizing pH controlled mass flow valves supplementing CO2 as needed into the continuous air supply. Optical density was continually monitored utilizing a Bugeye reading at 850 nm and cells were harvested during late log growth by centrifugation. Genomic DNA was isolated and purified utilizing the Qiagen DNeasy plant maxi kit. Cells were lysed by extraction in the Avestin Emulsiflex-B15 homogenizer at 30,000 psi prior to purification.
Nannochloropsis oculata (CCMP 525) was axenically maintained in 2.8 L wide-mouth Fernbach flasks that contained 1,000 ml F/2 medium . The flasks were plugged with cheesecloth-covered, hand rolled cotton stoppers and capped with #2 Kraft autoclave bags (Paper Mart, Orange, CA.). Cultures were maintained at 20°C on a 12 h light: 12 h dark photoperiod at 100 μEm-2 s-1 light intensity using full spectrum T12 fluorescent light bulbs (Pacific Lamp Supply Co., Seattle, WA.). Cell counts were accomplished using an Accuri C6 flow cytometer (BD Scientific, Ann Arbor, MI). Cultures were harvested at early stationary phase of growth and total high molecular weight DNA (greater than 500 kb in size) was extracted from N. oculata using the Qiagen Genomic-Tip 500G kit according to manufacturer’s directions (Qiagen, Valencia, CA, USA).
Sequencing and assembly
N. salina chloroplast and mitochondrial genomes were sequenced using a combination of Illumina  and 454 sequencing technologies . A 1 X 100 base pair shotgun library was prepared using standard TruSeq protocols and sequenced from bulk N. salina genomic DNA on an Illumina HiSeq2000 sequencer, generating approximately 100 million reads. Additional shotgun single-end and paired-end (11 kb insert) DNA libraries were prepared for sequencing on the 454 Titanium platform, generating 0.807 million and 3.23 million reads, respectively. The 454 single-end data and the 454 paired end data (insert size 4720 +/- 1180 bp) were assembled together using Newbler, version 2.3 (release 091027_1459). The Illumina-generated sequences were assembled separately with VELVET, version 1.0.13 . The resulting consensus sequences from both the VELVET and Newbler assemblies were computationally shredded into 10 kb fragments and were re-assembled with reads from the 454 paired end library using parallel phrap, version 1.080812 (High Performance Software, LLC). The chloroplast and mitochondrial replicons were identified in this final hybrid assembly based on: a) increased coverage in the 454 paired-end library (> 20 times higher than nuclear genome reads), b) the absence of paired end links to other contigs in assembly and, c) verification via homologous blast searches against the N. gaditana chloroplast and mitochondrial genomes. Sequence reads that belonged to each respective organelle were removed from the main project and re-assembled separately. Mis-assemblies in the contigs/scaffolds were corrected using gapResolution (Cliff Han, unpublished script) or Dupfinisher (Han, 2006) and repeat resolution was performed in Consed to generate the final circular consensus sequence. The final, fully assembled chloroplast genome was supported by > 500x average coverage from both sequencing platforms.
N. oculata chloroplast genome was sequenced by constructing large-insert fosmid clones from high molecular weight DNA as previously described in Raymond et al.  and as adapted in Cattolico et al. . Clones were plated using 12 μg/mL chloramphenicol selection, picked using the Q-pix automated colony picker (Genetix Ltd. UK) and inoculated into 384-well glycerol stock freezing plates.
Fosmid DNA was recovered using a standard alkaline-lysis protocol, and sequenced using standard dye-termination methods and capillary electrophoresis according to ABI manufacturer’s directions using a 3730xl Genome Analyzer. Vector sequences were removed and sequences were further trimmed to remove low quality bases. Sequences were compared to a custom database consisting of published chloroplast genomes using BLASTX. Fosmids in which both end sequences had high quality matches (E value < 10-4) to a chloroplast gene as judged by both BLAST analyses were identified as chloroplast-derived. All fosmid end sequences are available on our web site database (http://chloroplast.ocean.washington.edu). A total of fourteen 384-well plates were sequenced from three independent library preparations. Of those, 41 clones had end-sequences with chloroplast signatures, and these were subjected to Multiple Complete Digest (MCD) restriction analysis. Clones were analyzed by MCD analysis as previously described [19, 21]. Fosmid clones were digested using HindIII, BglII, NsiI, and EcoRI, subjected to electrophoretic separation on a 0.8% agarose gels, and visualized using a Typhoon 8600 Variable Mode Imager (Amersham Biosciences, Piscatawny, NJ). Automated band calling was performed using QGAP software (Quantitative Gel Analysis Program). Restriction data were analyzed using GenVal software  that compares DNA fingerprints and aligns end-sequence data for multiple clones, either against a reference genome or de novo. For this genome, three fosmids were initially sequenced that appeared to be spatially positioned to maximize genome coverage. Following sequencing and finishing (see below), two additional clones were selected for sequencing to extend the contig, but they did not complete the genome. Final finishing of the sequence was performed using experiments designed by Autofinish . Each fosmid clone was finished (mis-assemblies resolved, weak regions and gaps closed) separately and then assembled in Consed. Final validation was completed by expert finishers at the University of Washington using the MCD data from the fosmids. A final gap of approximately 15 kbp gap was not covered by fosmids and was closed by sequencing PCR products that were generated using primers designed using the N. oceanica genome .
The N. oculata mitochondrial genome was sequenced to ~50X coverage using the Illumina Hiseq 2000 according to manufacturer’s instructions. A paired end shotgun library was prepared from total genomic DNA using the Illumina Nextera DNA sample preparation kit (Catalog #FC-121-1030) using dual indexing . A total of 2.5 million 60mer reads were recovered following demultiplexing. The reads were assembled using Velvet version 1.2.03 . Assembly parameters were determined empirically using a custom script which explored velvet parameter space and compared resulting assemblies against the mitochondrial genome of N. oceanica. Once optimum Velvet parameters were determined, all contigs greater than 1000 bp were annotated using a custom auto-annotation pipeline and mitochondrial contigs were easily identified.
To enable an accurate comparative analysis, all chloroplast and mitochondrial genome sequences were syntenously aligned and reoriented upstream of the rrn23S and cox1 codons, respectively, and annotated using the same methods and employing standard genetic codes for translating open reading frames (ORFs). ORFs were initially predicted using Glimmer 3.0 , ribosomal RNA genes were identified with RNAmmer , and genes for tRNAs were identified using tRNASCAN-SE . Searches for tmRNAs and the signal recognition particle RNA employed ARAGORN  and SRPscan . Predicted gene functions were initially assigned using a BLASTP search of a custom chloroplast or mitochondrial genome database and refined manually with the aid of conserved protein motifs identified using the PFAM database . Tandem repeats were found with Tandem Repeat Finder  using default settings. Inverted repeats were found with E-inverted from the EMBOSS package  using the default settings and the additional constraint that repeats had to be more than 80% similar and the length of the loop shorter than the stem. Repeats were further examines using M-fold (http://mfold.rna.albany.edu/?q=mfold/DNA-Folding-Form) using default settings. Circular genome maps were created with OGDRAW . Manual corrections to the above automated structural and functional assignments were completed on an individual gene-by-gene basis as needed.
Sequences and annotations for the chloroplast and mitochondria genomes in Ectocarpus siliculosus , Aureococcus anophagefferens , Thalassiosira pseudonana [1, 35], Phaeodactylum tricornutum [1, 35], and Heterosigma akashiwo  were retrieved from Genbank [http://www.ncbi.nlm.nih.gov/pubmed].
Protein translations of all ORFs found on the Nannochloropsis organellar genomes were subjected to BLASTP searches against the NCBI non-redundant (NR) protein database (version 2012.10.19). Genes were considered ‘divergent’ based on the following criteria: a) e-value of the best hit was >1e-20, and b) the query and subject lengths varied by >20% or the aligned portion of the proteins had <40% similarity to the closest blast hit. Nucleotide alignments of all replicons were completed using Mauve 2.3.1 and the EMBOSS Stretcher  pairwise sequence alignment tool with default settings. Multi-protein sequence alignments were completed in MEGA  employing MUSCLE algorithms. Tertiary structure prediction for Nsk00142 (‘clpN’) was completed using I-Tasser  with default settings. Primary amino acid sequences alignments for AtpD were completed using clustalW and manually curated based on tertiary structure predictions (see Methods below). Phylogenetic trees of CbbX employed RAxML v 7.2.8 using 400 amino acid positions (excluding the C-terminal extension possessed only by Nannochloropsis) with rapid bootstrapping, a gamma model of rate heterogeneity and the RTREV substitution matrix.
Ab initio modeling  and comparative modeling  were completed using Rosetta to garner insight on structural changes encoded by the atpD, atpG and atpA-N terminus (first 20 amino acids of atpA sequence). Three and nine amino acid fragments were created from the protein database using the ROSETTA server . Secondary structure predictions were made for the sequences using psipred . For ab initio structure predictions, 16050 trajectories were run for AtpD sequence and 20400 trajectories for AtpG and the N terminus of AtpD. The models were clustered based on their RMSD and the top 20 clusters based on the total-score were visually evaluated.
For comparative modeling, the crystal structure of Ecoli delta-subunit (PDB code: 1abv)  and bovine OSCP (PDB code: 2bo5)  were used as templates. The sequence alignment of N. salina AtpD was completed using ClustalX  and gaps were removed manually based on secondary structure predictions for N. salina protein sequences using psipred and the secondary structure observed in E coli AtpD and bovine OSCP structures. The alignment was adjusted to place the gaps in the loop region of the template structures. 15300 trajectories were run against each template structure. The secondary structures of the extreme N and C terminal regions of all AtpD homologs not covered by crystal structures were predicted bases on consensus predictions from psipred  and Porter . Similarly, N. salina AtpG was modeled against multiple structures from PDB that had close sequence homology with N. salina atpG identified using HHpred server . A subset of structural hits were used as templates for comparative modeling (PDB codes: 3V6I, 1B9U, 1L2P, 2KHK, 2CLY, 2K88, 2KK7, 3VOU).
To determine if the N. salina atpA and its predicted protein structure would interact with AtpD and be in a similar orientation as the E. coli complex and bovine complexes, the N-terminus fragment (20 residues) from the NS-atpA sequence was used to estimate the structure. The top structural hit on HHpred server (PDB code: 3KKR) was used as a template for comparative modeling of NS-atpA N-terminus sequence. The lowest energy predicted structure was then docked in the expected pocket in the predicted AtpD structure using Rosetta docking protocol [48, 49]. A total of 10200 dock trajectories were run. During the dock, the predicted AtpD structure was truncated at the C-terminus beyond the structural overlap with the template PDB (1abv in this case). Using gnuplot , the total predicted full-atom energy  of each complex was then plotted aginst the RMS deviation of each complex from the best full-atom energy complex. Another random pocket on the predicted AtpD was chosen and 30600 trajectories of docking were completed, where the AtpA-N terminus structure was randomly positioned all over the AtpD predicted structure. The total predicted full-atom energy was recorded for each docked conformation (Additional file 1: Figure S1).
Transcript preparation, sequencing, and analysis
Nannochloropsis salina cells were grown as indicated above (see ‘Culturing and DNA purification’ methods). Samples (10 mL) were robotically removed on days 8, 9, and 13 during a N deprivation experiment, centrifuged at 3500 X g, flash frozen and stored at -80°C. Total RNA was extracted as follows: cells were lysed by addition of 3 mL ice cold Trizol with 1% w/v laurylsarcosine, passed three times through a cold Avestin pressure homogenizer at 36000 psi then vortexed after the addition of 750 μL of chloroform. The solution was allowed to equilibrate for 5 minutes, and then phase separated by centrifugation with phase lock gel at 13000 g for 10 min at 4°C. The aqueous phase was mixed with 100% ethanol to a final concentration of 70% and applied to an Invitrogen PureLink mRNA column (Life Technologies, Carlsbad, CA). Residual DNA was removed on column by treating with DNase.
Total RNA from each time point was separated into two aliquots. One aliquot was subjected to poly-A selection by hybridizing to poly-T coated beads using the Invitrogen Fastrack MAG mRNA Isolation Kit (Part number 45–7000; Life Technologies, Carlsbad, CA) according to the manufacturers instructions. Ribosomal RNA was removed from the second aliquot using both the Plant Leaf and Bacteria RiboZero rRNA Removal Kits (Part numbers MRZPL116 and MRZMB126; Epicentre, Madison, WI). Following these pre-treatments, both RNA aliquots were prepared for shotgun sequencing (2 X 100 base pairs) using the ScriptSeq v2 RNA-seq Library Preparation Kit (SSV21124, Epicentre, Madison, WI) and sequenced on the Illumina Hiseq 2000 platform , generating approximately 20 million reads per sample. Sequence reads were quality trimmed on both ends (Q > 10 sliding window), mapped to the N. salina chloroplast and mitochondrial genomes using Bowtie2, and RPKM values were calculated for each gene using Artemis . Resultant transcript expression profiles (.bam files) were visualized via Artemis and/or IGV [52, 53].
Results and discussion
Global characteristics and interspecies comparisons
General characteristics of the Nannochloropsis organellar genomes
N. gaditana #
N. oculata #
Nucleotide identity (%)†
Nucleotide identity (%)†
Seven and three novel ‘ORFans’ were annotated in N. salina and N. oculata, respectively. Two of these ORFans NskMp00219 and NskMp00232 are conserved in all four species examined and encode proteins of 323 and 231 amino acids, respectively. Based on BLASTP analysis, both of these genes do not have homologs (outside of the Nannochloropsis) in the NCBI non-redundant protein sequence database. Both genes appear to be transcribed as sequence reads from the transcriptome mapped to these regions. Unfortunately, tertiary structure analysis of the proteins encoded by either gene did not produce analogs with high structural similarity scores (data not shown).
Pangenomes of the Nannochloropsis organelles
Chloroplast* (138 genes)
Mitochondria* (48 genes)
psaA psaB psaC psaD psaE psaF psaI psaJ psaL psaM psb30 psbA psbB psbC psbD psbE psbF psbH psbI psbJ psbK psbL psbN psbT psbV psb28 psbX psbY psbZ
Cytochromes, chlorophyll, cofactor biosynthesis
chlI chlL chlN chlB acsF ycf54 ccsA css1 petA petB petD petF petG petJ petL petM petN thiG thiS
rbcL rbcS cbbX ilvB acpP
atpA atpB atpD atpE atpF atpG atpH atpI
atp1 apt6 apt8 atp9 cob cox1 (2) cox2 cox3 nad1 nad2 nad3 nad4 nad4L nad5 nad6 nad7 nad9 nad10
rpl1 rpl2 rpl3 rpl4 rpl5 rpl6 rpl11 rpl12 rpl13 rpl14 rpl16 rpl18 rpl19 rpl20 rpl21 rpl22 rpl23 rpl27 rpl29 rpl31 rpl32 rpl33 rpl34 rpl35 rpl36 rps2 rps3 rps4 rps5 rps6 rps7 rps8 rps9 rps10 rps11 rps12 rps13 rps14 rps16 rps17 rps18 rps19 rps20 tufA
rpl14 rpl16 rpl2 rpl5 rpl6 rps 2 rps3 rps4 rps7 rps8 rps10 rps11 rps12 rps13 rps14 rps19
ftsH dnaK groEL secA secY tatC clpC1 clpC2 (2) clpN sufB/ycf24 sufC/ycf16
rpoA rpoB rpoC1 rpoC2
ycf3 ycf4 ycf34 ycp36 ycf44 ycf46 ycf49 ycf66
Novel hypotheticals/ORFans +
Nsk00013 Nsk00014 Nsk00050b Nsk00053 Nsk00085 Nsk00087 Naoc00026 Naoc00069
Nsk00219 Nsk00223 Nsk00232
Nsk00202 Nsk00203 Nsk00204 Nsk00206 NaocMp0016 Naon00225 NaonMp0027 NaonMp0040
All Nannochloropsis strains encode small inverted repeats throughout their chloroplast genomes, almost exclusively within intergenic regions (Figures 2 & 3). Cruciform arrangements, formed by inverted repeats, represent alternative DNA structural elements that are known to impact a wide variety of cellular processes, including DNA replication, repair, protein association and gene expression. M-folding show these inverted repeats to have a very defined architecture wherein loop and stem sizes are highly conserved. Of the 66 inverted repeats examined, loop domains were found to be quite small. Seventy-four percent of the loop structures averaged 4.2+/- 0.8 bp in length; 23% were 7.7 +/- 1.7 bp in length while only 3% has a length of 11.5 +/- 0.7 bp. Stem size of the repeats appeared to fall into two categories. For example, those repeats servicing photosystem I genes (psaB, J and L) and energy conservation (atpH,G, E; petA,D, F), had an average stem length of 20.8 +/- 4.4 bp, while stem length of repeats servicing photosystem II genes (psbH, Y, N, I, T) had a longer length of 30.7 +/- 6.9 bp. Interestingly, psbH, N, and I also have the among the largest (~10 bp) loop domains. Several individual genes also have longer repeat stem structures. The rpoC2 and acfF (with stem lenths of 34 and 37 bp respectively) are good examples. Similar to bacterial gene regulation, we note that the small repeats may serve adjacent genes that are on opposite reading strands (e.g., petD-rpl12; petA-thiS; ccsA-rps6; psaJ- ThiG; petF- rps10). Such placement is often conserved for all four Nannochloropsis strains. Though these dual serving repeats are "shared" with nearest neighbors, we have found that specific genes, such as petD or rps10, appear to be targeted, regardless of taxon for repeat embellishment (data not shown; ). One may speculate that the proteins encoded by these genes are seminal players in photosynthesis or transcription and may be targets for regulation.
All four Nannochloropsis chloroplast genomes are divided into two approximately equal coding domains by the presence of a large inverted repeat (IR). The size of this repeat has been show to be strain dependent in Nannochloropsis . The N. salina and N. gaditana repeat encodes the 23, 16 and 5S ribosomal genes (4.9 kb). Confirming previous observations , N. oculata expands this repeat coding array to include three additional genes: clpC2, psbV and petJ (7.5 kb) (Figure 3 inset). Though an N. oceanica LAMB0001 repeat structure is evident and is likely similar to IR found in the other sequenced N. oceanica strains , the publicly available genome remains incompletely assembled in both of the repeat domains (Figure 1). As more genomes are completed, the new data suggests that stramenopile IR size may generally be taxon dependent with complete loss  or smaller IR’s (~ 6Kb) occurring within the eustigmatophytes, pelageophytes, pinguiophytes, and xanthophytes and larger repeats (10 – 22 kb) found in the raphidophytes and bacilliariophytes (http://chloroplast.ocean.washington.edu/home). Well-documented chloroplast genome IR size change has been extensively studied in the viridiplantae [58, 59]. Outside of the Nannochloropsis lineage, chloroplast genome strain comparisons have only been accomplished in the stramenopiles for Heterosigma akashiwo (strains CCMP 452 and NIES 293; ), and for species comparisons made between Thalassiosira oceanica (CCMP1005) and T. pseudonana (CCMP 1335; ). Unlike the observation in Nannochloropsis , no difference in large IR repeat size was observed either between strains or genera. Why chloroplast genomes maintain the IR domain remains undeciphered. Although the large IR structures promote the formation of molecular isomers within the chloroplast genome population [20, 60] via recombination, differential function for these isomorphic forms has not been determined. However, copy correction between IR domains may contribute to genome stability- an especially important fact when one considers that an algal cell may contain hundreds of chloroplast DNA molecules .
Due to the high level of protein similarity and synteny encoded in the organallar genomes of these Nannochloropsis species, we globally aligned each organellar genome and examined the relative nucleotide similarity of each species (Table 1, Figure 1). Although the N. oculata mitochondrial genome is closer in size to N. salina than N. oceanica, the nucleotide similarity of N. oculata is most similar to N. oceanica (N. oculata vs. N. salina; 76.2%, N. oculata vs. N. oceanica; 87.7%). Similarly, the entire N. salina and N. gaditana mitochondrial genomes share a 97% nucleotide identity and a 100% conservation of gene synteny. With respect to the chloroplast genomes, the N. salina and N. gaditana replicons only differ by 75 bp, are 98.4% identical at the nucleotide level, and contain an identical inventory of open reading frames. In contrast, the N. oculata chloroplast sequence is only 84.3% identical to N. salina but is 92.4% identical at the nucleotide level to N. oceanica. Taken together, these data indicate that the N. salina and N. gaditana replicons are more similar to each other than they are to the organellar genomes found in N. oculata and N. oceanica which is consistent with previously known evolutionary relationships  and a recent phylogenomic study of these organelles .
Because of the high degree of nucleotide similarity in the organellar genomes of N. salina and N. gaditana, a re-assessment of the phylogentic placement of N. gaditana is warranted. To our knowledge, no dogma has been established to phylogenetically classify single cell eukaryotes strictly based on the degree of nucleotide variation in highly conserved genes. As a general rule in bacteria, if two different bacterial isolates contain 16 s rDNA genes that are ≥ 97% similar, they are classified as the same species. The chloroplast ribosomal RNAs in N. salina and N. gaditana only differ by 7 nucleotides (99.76% identical). As stated above, we observed ≥ 97% nucleotide similarity across the entire mitochondrial and chloroplast replicons of N. salina CCMP1776 and N. gaditana CCMP526. If NADH dehydrogenase subunit 5 (nad5) is used as a strain discriminator, a higher resolution among organisms can be achieved (Black and Cattolico, unpublished). Little difference in nucleotide sequence diversity in nad5 is observed when either N. salina and N. gaditana (1.6% difference) or N. oculata and N. oceanica (5.8% difference) are compared. In contrast, comparisons between N. salina and N. oceanica or N. gaditana and N. oculata nad5 indicate 14.6% and 15.5% sequence variation, respectively. These data provide further support that N. salina and N. gaditana are closely related. In conclusion, the identical gene synteny and high degree of nucleotide identities suggest that N. gaditana could be reclassified as a strain of N. salina (i.e. "Nannochloropsis salina strain gaditana"). The availability and consequent comparative analysis of the nuclear genomes from both isolates will undoubtedly provide clarifying evidence to support this proposition.
To gain further insight into the unique features conserved within the genus Nannochloropsis, we compared the gene content of the Nannochloropsis organellar pangenomes (Table 2) to other representative sequences found in the same phyla (Thalassiosira pseudonana (Coscinodiscophyceae), Phaeodactylum tricornutum (Bacilliariophyceae), Ectocarpus siliculosus (Phaeophyceae), Aureococcus anophagefferens (Pelagophyceae), and Heterosigma akashiwo (Raphidophyceae) . Within this set of stramenopiles, the global gene inventory of Nannochloropsis is most similar to H. akashiwo (data not shown). Unlike A. anophagefferens, P. tricornutum, T. thalassiosira, and H. akashiwo, the Nannochloropsis and Ectocarpus chloroplast genomes both contain chlB, chlN, chlL (light independent protochlorophyllide reduction), the acsF/chl27 (Mg-protoporphyrin IX monomethyl ester cyclase) gene as well as ycf54 (demonstrated to play a critical role in AcsF synthesis/maturation or in the process of cyclase assembly ). This gene assemblage suggests that these stramenopile genera (Nannochloropsis, Ectocarpus, and others with the same gene complement) may share similar mechanisms of chlorophyll biosynthesis. The Nannochloropsis chloroplasts have also maintained single copies of petJ, ycf49, ycf36, genes more typically conserved in cyanobacteria, rhodophytes and some stramenopiles (e.g., xanthophytes and raphidophytes) but are usually found to be transferred to the nucleus in the bacilliariophytes. Additionally, all four Nannochloropsis mitochondria encode atp1, a subunit of the F1F0 ATP synthase. This gene is absent in all other stramenopile mitochondrial genomes sequenced to date.
In all four Nannochloropsis mitochondria, the gene which encodes for subunit ‘G’ of the NADH dehydrogensase, nad11, is shorter than what is canonically known, containing only the molybdopterin cofactor binding domain but lacking the NADH iron-sulfur (Fe-S) binding region. A gene which encodes a very similar Fe-S binding domain was located in the drafted N. salina nuclear genome, indicating that this portion of the protein is now encoded by the nuclear genome. In P. littoralis, the opposite transfer occurred as only the Fe-S domain is present in the mitochondrial genome and the molybdopterin binding domain is encoded in the nucleus . The fact that the P. tricornutum nad11 is split into two parts corresponding to these two domains in N. salina and P. littoralis, but that the domains still reside on the mitochondrial genome [1, 35] suggests that this protein is a vulnerable target for nuclear transfer.
As in all other chloroplasts , many structural subunits of Photosystems (PS) I and II are conserved in the Nannochloropsis chloroplast genomes. Nevertheless, the PS subunits that have been lost from the chloroplast (through migration or deletion) follow previous deletion patterns observed in several stramenopile and rhodophytic representatives. Similar to what has been observed in stramenopiles , the PSI subunit genes psaG, psaH, psaK, psaN, psaO, psaP, psaX and the PSII subunits psbM, psbP, psbQ, psbR, psbS have been removed from the chloroplast genomes of Nannochloropsis. As seen in rhodophytic algae , the genes encoding PsbO, PsbU, and Psb27 are also absent in the Nannochloropsis genome.
Carbon dioxide fixation in Nannochloropsis is mediated by a ‘red-type’ Form 1 ribulose-1,5-bisphosphate carboxylase-oxygenase (RuBisCO) , shown to have a high affinity for CO2 yet a low specificity factor due to poor discrimination between O2 and CO2 . Single copies of rbcL and rbcS are found on each of the chloroplast genomes but a gene which encodes for the transcriptional regulator, rbcR was not identified in N. salina and is consistent with previous observations in Nannochloropsis [12, 13] and in some of the other stramenopiles (; unpublished results). In viridiplanta and some algae, RuBisCO is post-translationally regulated via nitrosylation of conserved cysteine residues in RbcL (Cys 460 and Cys181 in G. suphuraria), resulting in inactivation of the enzyme at the active site [68–70]. Interestingly, the Nannochloropsis RbcL does contain a cysteine at position 460 but does not encode a cysteine near the active site at position 181, which suggests that this type of post-translational control may not be functioning in Nannochloropsis.
The Nannochloropsis chloroplast pangenome contains an ortholog of the large subunit of an acetohydroxyacid synthase (i.e. ilvB, Nsk0066), which is the only known enzyme to catalyze the first step in biosynthesis of branched chain amino acids; valine, leucine and isoleucine. Surprisingly, the accompanying ‘small subunit’ regulator, ilvH/N, required for negative feedback regulation and optimum activity [71–74], appears to have been uniquely lost from this genera as an ortholog of ilvH was not identified in any of the sequenced Nannochloropsis genomes (nuclear, mitochondria, or plastids). With respect to all publicly available stramenopile choloroplast genomes, either a.) ilvB and ilvH have both been maintained (i.e.,. H. akashiwo, E. siliculosus, A. anophagefferens) or b.) both subunits have been transferred to the nuclear genome (i.e., T. psuedonana and P. tricornutum). Searching broadly across photosynthetic organisms in other eukaryotic phyla, we could not identify another instance where ilvH or ilvB had been lost from any chloroplast genome independent of its partner gene, which is a striking occurance considering a recent review indicated all known acetohydroxyacid synthases contain both subunits . Therefore, the absence of ilvH suggests that Nannochloropsis has either lost its ability to negatively regulate IlvB or has evolved a novel regulator.
Highly divergent genes on the Nannochloropsis organellar genomes
rps6; 30S ribosomal protein S6
30S ribosomal protein S6 [Thalassiosira pseudonana];|YP_874616.1|
atpD; Atp synthase delta subunit
Hypothetical protein MldDRAFT_4321 [delta proteobacterium MLMS-1]; |ZP_01290127.1|
atpF; ATP synthase b subunit
CF0 subunit I of ATP synthase [Oltmannsiellopsis viridis]; |YP_635887.1|
atpG; ATP synthase b’ subunit
ATP synthase CF0 subunit II [Vaucheria litorea];|YP_002327468.1|
Hypothetical; putative peroxidase
Hypothetical protein tlr1577 [Thermosynechococcus elongatus BP-1];|NP_682367.1|
psb28; photosystem II protein (ycf79)
Photosystem II protein W [Guillardia theta]; |NP_050669.1|
ycf4; photosystem I assembly protein
Photosystem I assembly protein Ycf4 [Coccomyxa subellipsoidea C-169]; |YP_004222004.1|
Unknown DUF2499 [Picea sitchensis]; |ABK25760.1|
Hypothetical protein SPPN_02855 [Streptococcus pseudopneumoniae IS7493
rpoA; RNA polymerase alpha chain
RNA polymerase alpha subunit [Cryptomonas paramecium]; |YP_003359271.1|
Chloroplast protein Ycf34 [Gloeobacter violaceus PCC 7421]; |NP_927340.1|
ATP-dependent Clp protease ATP-binding subunit ClpA [Desulfobulbus propionicus DSM 2032; |YP_004196194.1|
Predicted protein with ABC transporter signatures [Fibroporia radiculosa]; |CCM01526.1|
Hyp. periplasmic binding protein MARHY3762 [Marinobacter hydrocarbonoclasticus ATCC 49840]; |YP_005431639.1|
Coiled-coil domain-containing protein 141 [Nomascus leucogenys]; |XP_003253834.1|
rps10; 30S ribosomal protein S10
30S ribosomal protein S10 [Spirochaeta smaragdinae DSM 11293]; |YP_003802682.1|
rps11; 30S ribosomal protein S11
30S ribosomal protein S11, partial [uncultured bacterium]; |EKD46317.1|
rps2; 30S ribosomal protein S2
Hypothetical protein [Batrachochytrium dendrobatidis JAM81];|EGF78568.1|
rps4: 30S ribosomal protein S4
Ribosomal protein S4 [Synedra acus]; |YP_003359457.1|
Hypothetical protein [Trichomonas vaginalis G3]; |XP_001579587.1|
rpl5; 50S ribosomal protein L5
Ribosomal protein L5 [Thalassiosira pseudonana]; |YP_316605.1|
atp8; ATP synthase F0 subunit 8
ATP synthase F0 subunit 8 [Fucus vesiculosus; |YP_448633.1|
fmhA protein [Staphylococcus saprophyticus ATCC 15305]; |YP_300577.1|
rps13; 30S ribosomal protein S13
NADH dehydrogenase s9- S13 fusion protein [endosymbiont of Durinskia baltica] |gb|AEP20701.1|
Nsk00013, Nsk0014, Nsk00150, Nsk00085, Nsk00203, Nsk00223
ATP synthase gene cluster
As indicated previously, many of the Nannochloropsis ATP synthase genes have diverged significantly from all other F1F0 type ATP synthases across the tree of life. Like most other algae, the Nannochloropsis chloroplast genomes encode an F1F0 type ATP synthase, a multimeric complex that catalyzes the synthesis of ATP from energy conserved through photosynthesis [76, 77]. The F1 complex (stator), which houses the catalytic site, is encoded by the alpha, beta, and delta subunits (AtpA, AtpB, AtpD) and interacts structurally with the central (gamma subunit) and peripheral (b/b’ subunits; AtpF/G) stalks to connect and stabilize the F1 to the membrane bound F0 complex (rotor). All previously published annotations of the Nannochloropsis chloroplast genomes genomes indicated that atpD was not present on the replicon [12, 13, 78], and analysis of the Nannochloropsis oceanica CCMP 1779 genome indicated that atpD was possibly located in the nuclear genome . Although atpD is often tandemly transferred to the nuclear genome with atpG in other algae , a close homolog of atpD could not be identified in the drafted nuclear genomes of N. gaditana, N. salina, or either N. oceanica genome. Because AtpD has been shown to be essential for function of the ATP synthase complex in yeast  and bacteria [81, 82] and similarly, loss of atpD expression in Arabidopsis disabled photoautotrophic growth , we hypothesized that a functional replacement (or a highly diverged ortholog) must be present on the chloroplast or nuclear genome.
Given the extreme level of divergence in the atpD nucleotide and amino acid translation, we also investigated changes in the main ATP synthase subunits known to interact with the delta subunit: AtpA and AtpG. Overall, the amino acid sequence alignments of the Nannochloropsis AtpA display a high level of conservation with other AtpA proteins (Additional file 3: Figure S3) yet, the N-terminal amino acids, which have been shown to interact with AtpD [43, 84, 85] have diverged. Similarly, the N-terminus of the b’ subunit (atpG), which anchors the protein in the chloroplast membrane, is conserved although the C-terminal end, which interacts with AtpD, aligns poorly with canonical AtpG sequences (Additional file 4: Figure S4).
Structure prediction and comparative modeling of the ATP synthase subunits
To gain insight into the tertiary structures of Ns-AtpD, Ns-AtpG and Ns-AtpA-N terminus, molecular docking and comparative modeling using known structures in the PDB database were conducted. Because the Ns-AtpD amino acid sequence was very divergent from any known structures, only low scoring homologies were observed on the HHpred server. Therefore, known structures of homologs from E. coli (PDB code: 1abv) and bovine (PDB code: 2bo5) AtpDs were used for comparative modeling. The predicted models for Ns-AtpD consistently acquired similar folds as those observed in the E. coli and bovine homologs (Figure 6A-C). With respect to Ns-AtpG, comparative modeling identified 43 PDB structures with some degree of sequence homology. The selected model for Ns-AtpG (Additional file 7: Figure S6F, Additional file 8) is based on the template structure 2K88 (Additional file 6: Figure S7B) and is similar to canonical AtpG structures with a long helix with breaks only towards the ends. This model strongly suggests that NS-AtpG sequence is an ortholog of the b’ subunit of the ATP synthase. The Ns-AtpA-N terminus model was arbitrarily placed in the proximity of Ns-AtpD model between the helices which correspond to the helices that interact with AtpA sequence in E. coli and bovine complexes (Figure 7B & C, Additional file 9). Rigorous random local docking accompanied by complete randomization of the Ns-AtpA N terminus generated a top scoring conformation similar to known AtpA-AtpD interactions. Another set of 30600 trajectories of local docking but from a different starting point failed to produce a Ns-AtpD/Ns-AtpA-N terminal conformation with a better total energy of the complex.
Approaching this analysis critically, ab initio modeling minimizes the structural energy by producing maximum interactions, which results in compacted AtpD, Atp-A, and AtpG structures (Additional file 7: Figure S6A-C) that deviate from the structures of the known homologs. In an ATP synthase structure, there are multiple subunits and each interact with one another to provide a stable complex . Thus, in the absence of intermolecular domain-domain interactions, the predicted top scoring ab initio tertiary structures are likely artificial. Nevertheless, given that the ab initio secondary structure predictions and the the comparative modeling of tertiary structure were remarkably similar to known homologs, and the fact that the modified Nannochloropsis atpD is transcribed and present in the same canonical location, strongly suggests that Nsk00027 encodes a functional AtpD.
Clp protease complex
ATP dependent chaperone-protease complexes (Clp) play a critical role in protein homeostasis in both photosynthetic and non-photosynthetic bacteria and eukaryotes. All extant Clp complexes contain two functional elements: a chaperone protein and a proteolytic core. The bacterial chaperones (or ‘unfoldases’) ClpA, ClpC, and ClpX are members of the Clp/Hsp100 family of AAA + proteins, which function to recognize, unfold, and deliver polypeptides to the ClpP protease for degradation. Functional ClpCP complexes require an adaptor, MecA, to recruit specific protein substrates to ClpC . Similarly, the related but distinct ClpAP complex utilizes an adaptor, ClpS, to recruit N-end rule substrates to ClpA [87, 88].
Intriguingly, a homolog for a MecA adaptor gene (as indicated above, MecA interacts with ClpC chaperones) could not be identified on any N. salina replicon although two genes containing the conserved domains for ClpS (which normally interact with ClpA type chaperones) were found in the drafted nuclear genome assembly. With respect to ClpP, the Nannochloropsis chloroplast and mitochondria pangenomes are likewise devoid of genes which encode the ClpP protease, yet the drafted nuclear genome was found to contain five separate ORFs with putative ClpP protease domains (Additional file 10: Table S1).
To our knowledge, this is the first observation of a complete disassembly of individual ClpC domains into separate reading frames in any organism. Although it is currently unknown if these new chloroplast encoded ‘subunits’ still function collectively with the other nuclear encoded Clp components to create an active protease complex, it is still interesting to speculate on how the Clp homologs present in Nannochloropsis may interact and/or how these modifications change the function of the proteins. If we first envisage a traditionally functioning ClpCP, Nannochloropsis could have adapted to utilize ClpS as an adaptor for ClpC due to the absence of MecA (and other known bacterial) orthologs. This suggestion is not without precedent since interactions between MecA and the N-domain of ClpC were shown to resemble those of ClpS and the N-domain of ClpA , and in the cyanobacteria Synechococcus elongatus (which is also devoid of a MecA), ClpC was shown to interact directly with ClpS in vitro . Furthermore, because the N-domain is thought to partially mask the pore, separation of the ‘clp-N’ domain could increase the degradation efficiency by other mechanisms. For example, in the ClpAP system, SsrA-tagged substrates compete with ClpS recognized proteins for delivery to the unfoldase . If an SsrA-dependent system was present in Nannochloropsis, physical separation of the N domain would enable unhindered access to the active site and freely enable ClpS-independent proteins to be degraded.
Because the Nannochloropsis ClpC2 does not contain an obvious ClpP binding loop and the M-domain in D1 is either modified (or missing), we must also consider the alternative that a canonical ClpCP complex does not function in Nannochloropsis and that the single domain ClpC proteins have developed specialized functions and may act independent of adaptors. The Nannochloropsis ClpC2 protein has a conserved D2 loop, and if hexamerized, may continue to function as an unfoldase and promiscuously deliver substrates to the ClpP protease. Furthermore, as has been shown for several bacterial ClpC orthologs [91, 93, 94], clpC1 could compliment the activity of clpC2 (acting dependent or independent of adaptors) by stabilizing and preventing aggregation of newly synthesized, unfolded proteins; a function that is essential to effectively assemble (or dispose of) large multimeric complexes in the chloroplast. Clearly, further experimentation is required to determine if these or other scenarios explain the functional role of these novel Clp orthologs.
A pangenomic comparison of the Nannochloropsis with other stramenopiles revealed an extreme divergence in several key metabolic genes/systems: amino acid synthesis, carbon fixation, energy conservation, and protein homeostasis. These observations and further discovery of (as yet) currently unidentified genetic and structural modifications to critical cellular components will explain the unique physiological properties found in the genus Nannochloropsis. It is worthy to note that the high degree of divergence in the amino acid sequences of many Nannochloropsis proteins led to false annotations. Thus, implementation of tertiary structure prediction during annotation will be crucial to improve de novo gene calls in all newly sequenced organisms. Finally, the extraordinary similarity of the N. salina and N. gaditana organellar genomes suggests that these two isolates should be reclassified as different strains of the same species.
Availability of supporting data
The protein models (.pdb files) for AtpG and AtpD supporting the results of this article are available as Additional files 8 and 9, respectively. The genome assemblies and annotation data sets for each organelle are available in the GenBank repository; N. salina CCMP1776 organelles; accession numbers KJ410685 and KJ410689; N. oculata CCMP525 organelles; accession numbers KJ410684 and KJ410688; N. gaditana CCMP526 organelles; accession numbers KJ410682 and KJ410686; N. oceanica LAMB0001 draft organelles; accession numbers KJ410683 and KJ410687 (http://www.ncbi.nlm.nih.gov/genbank). The transcript mapping data shown in Figure 6 is available in the NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/bioproject/242770).
Protein data bank
Oligomycin-sensitivity conferring protein
Open reading frame
ATPases associated with various cellular activities
Root mean square deviation.
We would like to thank the Sequencing Technologies Team at Los Alamos National Laboratory for generating the raw sequence reads for this project and particularly Momo Vuyisch, Kim McMurry, and Cheryl Gleasner for technical asssistance in preparing and sequencing the DNA and RNA-seq libraries. Special thanks are also extended to Chloe Deodato at the U. of Washington for support in the N. oculata sequence study as well as Matt Posewitz and Robert Jinkerson for providing the sequence and annotation files for the N. gaditana chloroplast and mitochondrial genomes. We are also thankful to the High Performance Computing facility at Los Alamos National Laboratory for computational resources. This work was funded by the US Department of Energy under Contract DE-EE0003046 awarded to the National Alliance for Advanced Biofuels and Bioproducts and the US National Science Foundation grants 0523756 and 0629521 to RAC and GR, respectively. Work completed by RJ was fund by the US Defense Threat Reduction Agency (DTRA) grant CBCALL12-LS6-1-0622, and was enabled by Los Alamos National Laboratory internal resource grant WSYN_BIO. RJ thanks Charlie E. M. Strauss for supervision and guidance.
- Green BR: Chloroplast genomes of photosynthetic eukaryotes. Plant J. 2011, 66: 34-44. 10.1111/j.1365-313X.2011.04541.x.PubMedView ArticleGoogle Scholar
- Andersen RA, Brett RW, Potter D, Sexton JP: Phylogeny of the Eustigmatophyceae based upon 18S rDNA, with emphasis on Nannochloropsis. Protist. 1998, 149: 61-74. 10.1016/S1434-4610(98)70010-0.PubMedView ArticleGoogle Scholar
- Rodolfi L, Chini Zittelli G, Bassi N, Padovani G, Biondi N, Bonini G, Tredici MR: Microalgae for oil: strain selection, induction of lipid synthesis and outdoor mass cultivation in a low-cost photobioreactor. Biotechnol Bioeng. 2009, 102: 100-112. 10.1002/bit.22033.PubMedView ArticleGoogle Scholar
- Chiu SY, Kao CY, Tsai MT, Ong SC, Chen CH, Lin CS: Lipid accumulation and CO2 utilization of Nannochloropsis oculata in response to CO2 aeration. Bioresour Technol. 2009, 100: 833-838. 10.1016/j.biortech.2008.06.061.PubMedView ArticleGoogle Scholar
- Van Vooren G, Le Grand F, Legrand J, Cuine S, Peltier G, Pruvost J: Investigation of fatty acids accumulation in Nannochloropsis oculata for biodiesel application. Bioresour Technol. 2012, 124: 421-432.PubMedView ArticleGoogle Scholar
- Bucy H, Marchese AJ: Oxidative stability of algae derived methyl esters. J Eng Gas Turbines Power Trans. 2012, 134:Google Scholar
- Simionato D, Block MA, La Rocca N, Jouhet J, Marechal E, Finazzi G, Morosinotto T: Response of nannochloropsis gaditana to nitrogen starvation includes a De novo biosynthesis of triacylglycerols, a decrease of chloroplast galactolipids and a reorganization of the photosynthetic apparatus. Eukaryot Cell. 2013, 12 (5): 665-76. 10.1128/EC.00363-12.PubMed CentralPubMedView ArticleGoogle Scholar
- Mohammady NG: Characterization of the fatty acid composition of Nannochloropsis salina as a determinant of biodiesel properties. J Biosci. 2011, 66: 328-332.Google Scholar
- Quinn JC, Yates T, Douglas N, Weyer K, Butler J, Bradley TH, Lammers PJ: Nannochloropsis production metrics in a scalable outdoor photobioreactor for commercial applications. Bioresour Technol. 2012, 117: 164-171.PubMedView ArticleGoogle Scholar
- Wahidin S, Idris A, Shaleh SR: The influence of light intensity and photoperiod on the growth and lipid content of microalgae Nannochloropsis sp. Bioresour Technol. 2013, 129: 7-11.PubMedView ArticleGoogle Scholar
- Hoffmann M, Marxen K, Schulz R, Vanselow KH: TFA and EPA productivities of Nannochloropsis salina influenced by temperature and nitrate stimuli in turbidostatic controlled experiments. Mar Drugs. 2010, 8: 2526-2545. 10.3390/md8092526.PubMedView ArticleGoogle Scholar
- Wei L, Xin Y, Wang D, Jing X, Zhou Q, Su X, Jia J, Ning K, Chen F, Hu Q, Xu J: Nannochloropsis plastid and mitochondrial phylogenomes reveal organelle diversification mechanism and intragenus phylotyping strategy in microalgae. BMC Genomics. 2013, 14: 534-10.1186/1471-2164-14-534.PubMed CentralPubMedView ArticleGoogle Scholar
- Radakovits R, Jinkerson RE, Fuerstenberg SI, Tae H, Settlage RE, Boore JL, Posewitz MC: Draft genome sequence and genetic transformation of the oleaginous alga Nannochloropis gaditana. Nat Commun. 2012, 3: 686-PubMed CentralPubMedView ArticleGoogle Scholar
- Pan K, Qin JJ, Li S, Dai WK, Zhu BH, Jin YC, Yu WG, Yang GP, Li DF: Nuclear monoploidy and asexual propagation of Nannochloropsis oceanica (Eustigmatophyceae) as revealed by its genome sequence. J Psychol. 2011, 47: 1425-1432.Google Scholar
- Guillard RR, Ryther JH: Studies of marine planktonic diatoms. I. Cyclotella nana Hustedt, and Detonula confervacea (cleve) Gran. Can J Microbiol. 1962, 8: 229-239. 10.1139/m62-029.PubMedView ArticleGoogle Scholar
- Bennett S: Solexa Ltd. Pharmacogenomics. 2004, 5: 433-438. 10.1517/146224126.96.36.1993.PubMedView ArticleGoogle Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer ML, Jarvie TP, Jirage KB, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, et al: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437: 376-380.PubMed CentralPubMedGoogle Scholar
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829. 10.1101/gr.074492.107.PubMed CentralPubMedView ArticleGoogle Scholar
- Raymond CK, Subramanian S, Paddock M, Qiu RL, Deodato C, Palmieri A, Chang J, Radke T, Haugen E, Kas A, Waring D, Bovee D, Stacy R, Kaul R, Olson MV: Targeted, haplotype-resolved resequencing of long segments of the human genome. Genomics. 2005, 86: 759-766. 10.1016/j.ygeno.2005.08.013.PubMedView ArticleGoogle Scholar
- Cattolico RA, Jacobs MA, Zhou Y, Chang J, Duplessis M, Lybrand T, McKay J, Ong HC, Sims E, Rocap G: Chloroplast genome sequencing analysis of Heterosigma akashiwo CCMP452 (West Atlantic) and NIES293 (West Pacific) strains. BMC Genomics. 2008, 9: 211-10.1186/1471-2164-9-211.PubMed CentralPubMedView ArticleGoogle Scholar
- Hayden HS, Lim R, Brittnacher MJ, Sims EH, Ramage ER, Fong C, Wu Z, Crist E, Chang J, Zhou Y, Radey M, Rohmer L, Haugen E, Gillett W, Wuthiekanun V, Peacock SJ, Kaul R, Miller SI, Manoil C, Jacobs MA: Large-insert genome analysis technology detects structural variation in Pseudomonas aeruginosa clinical strains from cystic fibrosis patients. Genomics. 2008, 91: 530-537. 10.1016/j.ygeno.2008.02.005.PubMed CentralPubMedView ArticleGoogle Scholar
- Gordon D, Desmarais C, Green P: Automated finishing with autofinish. Genome Res. 2001, 11: 614-625. 10.1101/gr.171401.PubMed CentralPubMedView ArticleGoogle Scholar
- Hayden HS, Gillett W, Saenphimmachak C, Lim R, Zhou Y, Jacobs MA, Chang J, Rohmer L, D'Argenio DA, Palmieri A, Levy R, Haugen E, Wong GK, Brittnacher MJ, Burns JL, Miller SI, Olson MV, Kaul R: Evolution of Burkholderia pseudomallei in recurrent melioidosis. PLoS One. 2012, 7: e36507-10.1371/journal.pone.0036507.PubMed CentralPubMedView ArticleGoogle Scholar
- Delcher AL, Bratke KA, Powers EC, Salzberg SL: Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics. 2007, 23: 673-679. 10.1093/bioinformatics/btm009.PubMed CentralPubMedView ArticleGoogle Scholar
- Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW: RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35: 3100-3108. 10.1093/nar/gkm160.PubMed CentralPubMedView ArticleGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964. 10.1093/nar/25.5.0955.PubMed CentralPubMedView ArticleGoogle Scholar
- Laslett D, Canback B: ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 2004, 32: 11-16. 10.1093/nar/gkh152.PubMed CentralPubMedView ArticleGoogle Scholar
- Regalia M, Rosenblad MA, Samuelsson T: Prediction of signal recognition particle RNA genes. Nucleic Acids Res. 2002, 30: 3368-3377. 10.1093/nar/gkf468.PubMed CentralPubMedView ArticleGoogle Scholar
- Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, Holm L, Sonnhammer EL, Eddy SR, Bateman A: The Pfam protein families database. Nucleic Acids Res. 2010, 38: D211-D222. 10.1093/nar/gkp985.PubMed CentralPubMedView ArticleGoogle Scholar
- Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999, 27: 573-580. 10.1093/nar/27.2.573.PubMed CentralPubMedView ArticleGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.PubMedView ArticleGoogle Scholar
- Lohse M, Drechsel O, Bock R: OrganellarGenomeDRAW (OGDRAW): a tool for the easy generation of high-quality custom graphical maps of plastid and mitochondrial genomes. Curr Genet. 2007, 52: 267-274. 10.1007/s00294-007-0161-y.PubMedView ArticleGoogle Scholar
- Le Corguille G, Pearson G, Valente M, Viegas C, Gschloessl B, Corre E, Bailly X, Peters AF, Jubin C, Vacherie B, Cock JM, Leblanc C: Plastid genomes of two brown algae. Ectocarpus siliculosus and Fucus vesiculosus: further insights on the evolution of red-algal derived plastids. BMC Evol Biol. 2009, 9: 253-10.1186/1471-2148-9-253.PubMed CentralPubMedView ArticleGoogle Scholar
- Ong HC, Wilhelm SW, Gobler CJ, Bullerjahn G, Jacobs MA, McKay J, Sims EH, Gillett WG, Zhou Y, Haugen E, Rocap G, Cattolico RA: Analyses of the complete chloroplast genome sequences of two members of the Pelagophyceae: Aureococcus anophagefferens Ccmp 1984 and Aureoumbra lagunensis Ccmp1507. J Phycol. 2010, 46: 602-615. 10.1111/j.1529-8817.2010.00841.x.View ArticleGoogle Scholar
- Oudot-Le Secq MP, Grimwood J, Shapiro H, Armbrust EV, Bowler C, Green BR: Chloroplast genomes of the diatoms Phaeodactylum tricornutum and Thalassiosira pseudonana: comparison with other plastid genomes of the red lineage. Mol Genet Genomics. 2007, 277: 427-439. 10.1007/s00438-006-0199-4.PubMedView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinforma. 2008, 9: 40-10.1186/1471-2105-9-40.View ArticleGoogle Scholar
- Bradley P, Misura KM, Baker D: Toward high-resolution de novo structure prediction for small proteins. Science. 2005, 309: 1868-1871. 10.1126/science.1113801.PubMedView ArticleGoogle Scholar
- Raman S, Vernon R, Thompson J, Tyka M, Sadreyev R, Pei J, Kim D, Kellogg E, DiMaio F, Lange O, Kinch L, Sheffler W, Kim BH, Das R, Grishin NV, Baker D: Structure prediction for CASP8 with all-atom refinement using Rosetta. Proteins. 2009, 77 (Suppl 9): 89-99.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim DE, Chivian D, Baker D: Protein structure prediction and analysis using the Robetta server. Nucleic Acids Res. 2004, 32: W526-W531. 10.1093/nar/gkh468.PubMed CentralPubMedView ArticleGoogle Scholar
- Jones DT: Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol. 1999, 292: 195-202. 10.1006/jmbi.1999.3091.PubMedView ArticleGoogle Scholar
- Wilkens S, Dunn SD, Chandler J, Dahlquist FW, Capaldi RA: Solution structure of the N-terminal domain of the delta subunit of the E. coli ATP synthase. Nat Struct Biol. 1997, 4: 198-201. 10.1038/nsb0397-198.PubMedView ArticleGoogle Scholar
- Carbajo RJ, Kellas FA, Runswick MJ, Montgomery MG, Walker JE, Neuhaus D: Structure of the F1-binding domain of the stator of bovine F1Fo-ATPase and how it binds an alpha-subunit. J Mol Biol. 2005, 351: 824-838. 10.1016/j.jmb.2005.06.012.PubMedView ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
- Buchan DW, Ward SM, Lobley AE, Nugent TC, Bryson K, Jones DT: Protein annotation and modelling servers at University College London. Nucleic Acids Res. 2010, 38: W563-W568. 10.1093/nar/gkq427.PubMed CentralPubMedView ArticleGoogle Scholar
- Pollastri G, McLysaght A: Porter: a new, accurate server for protein secondary structure prediction. Bioinformatics. 2005, 21: 1719-1720. 10.1093/bioinformatics/bti203.PubMedView ArticleGoogle Scholar
- Soding J, Biegert A, Lupas AN: The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 2005, 33: W244-W248. 10.1093/nar/gki408.PubMed CentralPubMedView ArticleGoogle Scholar
- Chaudhury S, Berrondo M, Weitzner BD, Muthu P, Bergman H, Gray JJ: Benchmarking and analysis of protein docking performance in Rosetta v3.2. PLoS One. 2011, 6: e22477-10.1371/journal.pone.0022477.PubMed CentralPubMedView ArticleGoogle Scholar
- Gray JJ, Moughon S, Wang C, Schueler-Furman O, Kuhlman B, Rohl CA, Baker D: Protein-protein docking with simultaneous optimization of rigid-body displacement and side-chain conformations. J Mol Biol. 2003, 331: 281-299. 10.1016/S0022-2836(03)00670-3.PubMedView ArticleGoogle Scholar
- Racine J: gnuplot 4.0: A portable interactive plotting utility. J Appl Econ. 2006, 21: 133-141. 10.1002/jae.885.View ArticleGoogle Scholar
- Rohl CA, Strauss CEM, Misura KMS, Baker D: Protein structure prediction using rosetta. Num Comput Methods. 2004, 383: 66-+-View ArticleGoogle Scholar
- Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J: ACT: the artemis comparison tool. Bioinformatics. 2005, 21: 3422-3423. 10.1093/bioinformatics/bti553.PubMedView ArticleGoogle Scholar
- Thorvaldsdottir H, Robinson JT, Mesirov JP: Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2012Google Scholar
- Lambowitz AM, Zimmerly S: Group II introns: mobile ribozymes that invade DNA. Cold Spring Harb Perspect Biol. 2011, 3: a003616-PubMed CentralPubMedView ArticleGoogle Scholar
- Oudot-Le Secq MP, Fontaine JM, Rousvoal S, Kloareg B, Loiseaux-De Goer S: The complete sequence of a brown algal mitochondrial genome, the ectocarpale Pylaiella littoralis (L.) Kjellm. J Mol Evol. 2001, 53 (2): 80-88. 10.1007/s002390010196.PubMedView ArticleGoogle Scholar
- Ikuta K, Kawai H, Muller DG, Ohama T: Recurrent invasion of mitochondrial group II introns in specimens of Pylaiella littoralis (brown alga), collected worldwide. Curr Genet. 2008, 53: 207-216. 10.1007/s00294-008-0178-x.PubMedView ArticleGoogle Scholar
- Ravin NV, Galachyants YP, Mardanov AV, Beletsky AV, Petrova DP, Sherbakova TA, Zakharova YR, Likhoshway YV, Skryabin KG, Grachev MA: Complete sequence of the mitochondrial genome of a diatom alga Synedra acus and comparative analysis of diatom mitochondrial genomes. Curr Genet. 2010, 56: 215-223. 10.1007/s00294-010-0293-3.PubMedView ArticleGoogle Scholar
- Wang RJ, Cheng CL, Chang CC, Wu CL, Su TM, Chaw SM: Dynamics and evolution of the inverted repeat-large single copy junctions in the chloroplast genomes of monocots. BMC Evol Biol. 2008, 8:Google Scholar
- Wicke S, Schneeweiss GM, dePamphilis CW, Muller KF, Quandt D: The evolution of the plastid chromosome in land plants: gene content, gene order, gene function. Plant Mol Biol. 2011, 76: 273-297. 10.1007/s11103-011-9762-4.PubMed CentralPubMedView ArticleGoogle Scholar
- Bourne CM, Palmer JD, Stoermer EF: Organization of the chloroplast genome of the fresh-water centric diatom Cyclotella meneghiniana. J Psychol. 1992, 28: 347-355.Google Scholar
- Udy DB, Belcher S, Williams-Carrier R, Gualberto JM, Barkan A: Effects of reduced chloroplast gene copy number on chloroplast gene expression in Maize. Plant Physiol. 2012, 160: 1420-1431. 10.1104/pp.112.204198.PubMed CentralPubMedView ArticleGoogle Scholar
- Hollingshead S, Kopecna J, Jackson PJ, Canniffe DP, Davison PA, Dickman MJ, Sobotka R, Hunter CN: Conserved chloroplast open-reading frame ycf54 is required for activity of the magnesium protoporphyrin monomethylester oxidative cyclase in Synechocystis PCC 6803. J Biol Chem. 2012, 287: 27823-27833. 10.1074/jbc.M112.352526.PubMed CentralPubMedView ArticleGoogle Scholar
- Oudot MP, Kloareg B, Loiseaux-de Goer S: The mitochondrial Pylaiella littoralis nad11 gene contains only the N-terminal FeS-binding domain. Gene. 1999, 235: 131-137. 10.1016/S0378-1119(99)00194-8.PubMedView ArticleGoogle Scholar
- Vanselow C, Weber AP, Krause K, Fromme P: Genetic analysis of the Photosystem I subunits from the red alga, Galdieria sulphuraria. Biochim Biophys Acta. 2009, 1787: 46-59. 10.1016/j.bbabio.2008.10.004.PubMedView ArticleGoogle Scholar
- Mueller-Cajar O, Stotz M, Wendler P, Hartl FU, Bracher A, Hayer-Hartl M: Structure and function of the AAA + protein CbbX, a red-type Rubisco activase. Nature. 2011, 479: 194-199. 10.1038/nature10568.PubMedView ArticleGoogle Scholar
- Tchernov D, Livne A, Kaplan A, Sukenik A: The kinetic properties of ribulose-1,5-bisphosphate carboxylase/oxygenase may explain the high apparent photosynthetic affinity of Nannochloropsis sp to ambient inorganic carbon. Isr J Plant Sci. 2008, 56: 37-44. 10.1560/IJPS.56.1-2.37.View ArticleGoogle Scholar
- Rumpho ME, Worful JM, Lee J, Kannan K, Tyler MS, Bhattacharya D, Moustafa A, Manhart JR: Horizontal gene transfer of the algal nuclear gene psbO to the photosynthetic sea slug Elysia chlorotica. Proc Natl Acad Sci USA. 2008, 105: 17867-17871. 10.1073/pnas.0804968105.PubMed CentralPubMedView ArticleGoogle Scholar
- Stec B: Structural mechanism of RuBisCO activation by carbamylation of the active site lysine. Proc Natl Acad Sci USA. 2012, 109: 18785-18790. 10.1073/pnas.1210754109.PubMed CentralPubMedView ArticleGoogle Scholar
- Moreno J, Garcia-Murria MJ, Marin-Navarro J: Redox modulation of Rubisco conformation and activity through its cysteine residues. J Exp Bot. 2008, 59: 1605-1614.PubMedView ArticleGoogle Scholar
- Sakihama Y, Nakamura S, Yamasaki H: Nitric oxide production mediated by nitrate reductase in the green alga Chlamydomonas reinhardtii: an alternative NO production pathway in photosynthetic organisms. Plant Cell Physiol. 2002, 43: 290-297. 10.1093/pcp/pcf034.PubMedView ArticleGoogle Scholar
- Barak Z, Chipman DM: Allosteric regulation in Acetohydroxyacid Synthases (AHASs)–different structures and kinetic behavior in isozymes in the same organisms. Arch Biochem Biophys. 2012, 519: 167-174. 10.1016/j.abb.2011.11.025.PubMedView ArticleGoogle Scholar
- Belenky I, Steinmetz A, Vyazmensky M, Barak Z, Tittmann K, Chipman DM: Many of the functional differences between acetohydroxyacid synthase (AHAS) isozyme I and other AHASs are a result of the rapid formation and breakdown of the covalent acetolactate-thiamin diphosphate adduct in AHAS I. FEBS J. 2012, 279: 1967-1979. 10.1111/j.1742-4658.2012.08577.x.PubMedView ArticleGoogle Scholar
- Lee YT, Duggleby RG: Regulatory interactions in Arabidopsis thaliana acetohydroxyacid synthase. FEBS Lett. 2002, 512: 180-184. 10.1016/S0014-5793(02)02253-6.PubMedView ArticleGoogle Scholar
- Duggleby RG, McCourt JA, Guddat LW: Structure and mechanism of inhibition of plant acetohydroxyacid synthase. Plant Physiol Biochem. 2008, 46: 309-324. 10.1016/j.plaphy.2007.12.004.PubMedView ArticleGoogle Scholar
- Boczar BA, Delaney TP, Cattolico RA: Gene for the ribulose-1,5-bisphosphate carboxylase small subunit protein of the marine chromophyte Olisthodiscus luteus is similar to that of a chemoautotrophic bacterium. Proc Natl Acad Sci USA. 1989, 86: 4996-4999. 10.1073/pnas.86.13.4996.PubMed CentralPubMedView ArticleGoogle Scholar
- von Ballmoos C, Wiedenmann A, Dimroth P: Essentials for ATP synthesis by F1F0 ATP synthases. Annu Rev Biochem. 2009, 78: 649-672. 10.1146/annurev.biochem.78.081307.104803.PubMedView ArticleGoogle Scholar
- von Ballmoos C, Cook GM, Dimroth P: Unique rotary ATP synthase and its biological diversity. Annu Rev Biophys. 2008, 37: 43-64. 10.1146/annurev.biophys.37.032807.130018.PubMedView ArticleGoogle Scholar
- Vieler A, Wu G, Tsai CH, Bullard B, Cornish AJ, Harvey C, Reca IB, Thornburg C, Achawanantakun R, Buehl CJ, Campbell MS, Cavalier D, Childs KL, Clark TJ, Deshpande R, Erickson E, Armenia Ferguson A, Handee W, Kong Q, Li X, Liu B, Lundback S, Peng C, Roston RL, Sanjaya , Simpson JP, Terbush A, Warakanont J, Zauner S, Farre EM: Genome, functional gene annotation, and nuclear transformation of the Heterokont Oleaginous Alga Nannochloropsis oceanica CCMP1779. PLoS Genet. 2012, 8: e1003064-10.1371/journal.pgen.1003064.PubMed CentralPubMedView ArticleGoogle Scholar
- Stoebe B, Kowallik KV: Gene-cluster analysis in chloroplast genomics. Trends Genet. 1999, 15: 344-347. 10.1016/S0168-9525(99)01815-6.PubMedView ArticleGoogle Scholar
- Prescott M, Bush NC, Nagley P, Devenish RJ: Properties of yeast cells depleted of the OSCP subunit of mitochondrial ATP synthase by regulated expression of the ATP5 gene. Biochem Mol Biol Int. 1994, 34: 789-799.PubMedGoogle Scholar
- Hazard AL, Senior AE: Defective energy coupling in delta-subunit mutants of Escherichia coli F1F0-ATP synthase. J Biol Chem. 1994, 269: 427-432.PubMedGoogle Scholar
- Stack AE, Cain BD: Mutations in the delta subunit influence the assembly of F1F0 ATP synthase in Escherichia coli. J Bacteriol. 1994, 176: 540-542.PubMed CentralPubMedGoogle Scholar
- Maiwald D, Dietzmann A, Jahns P, Pesaresi P, Joliot P, Joliot A, Levin JZ, Salamini F, Leister D: Knock-out of the genes coding for the Rieske protein and the ATP-synthase delta-subunit of Arabidopsis. Effects on photosynthesis, thylakoid protein composition, and nuclear chloroplast gene expression. Plant Physiol. 2003, 133: 191-202. 10.1104/pp.103.024190.PubMed CentralPubMedView ArticleGoogle Scholar
- Rees DM, Leslie AG, Walker JE: The structure of the membrane extrinsic region of bovine ATP synthase. Proc Natl Acad Sci USA. 2009, 106: 21597-21601. 10.1073/pnas.0910365106.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilkens S, Borchardt D, Weber J, Senior AE: Structural characterization of the interaction of the delta and alpha subunits of the Escherichia coli F1F0-ATP synthase by NMR spectroscopy. Biochemistry. 2005, 44: 11786-11794. 10.1021/bi0510678.PubMedView ArticleGoogle Scholar
- Wang F, Mei Z, Qi Y, Yan C, Hu Q, Wang J, Shi Y: Structure and mechanism of the hexameric MecA-ClpC molecular machine. Nature. 2011, 471: 331-335. 10.1038/nature09780.PubMedView ArticleGoogle Scholar
- Zeth K, Ravelli RB, Paal K, Cusack S, Bukau B, Dougan DA: Structural analysis of the adaptor protein ClpS in complex with the N-terminal domain of ClpA. Nat Struct Biol. 2002, 9: 906-911. 10.1038/nsb869.PubMedView ArticleGoogle Scholar
- Schuenemann VJ, Kralik SM, Albrecht R, Spall SK, Truscott KN, Dougan DA, Zeth K: Structural basis of N-end rule substrate recognition in Escherichia coli by the ClpAP adaptor protein ClpS. Embo Rep. 2009, 10: 508-514. 10.1038/embor.2009.62.PubMed CentralPubMedView ArticleGoogle Scholar
- Kirstein J, Moliere N, Dougan DA, Turgay K: Adapting the machine: adaptor proteins for Hsp100/Clp and AAA + proteases. Nat Rev Microbiol. 2009, 7: 589-599. 10.1038/nrmicro2185.PubMedView ArticleGoogle Scholar
- Tryggvesson A, Stahlberg FM, Mogk A, Zeth K, Clarke AK: Interaction specificity between the chaperone and proteolytic components of the cyanobacterial Clp protease. Biochem J. 2012, 446: 311-320. 10.1042/BJ20120649.PubMedView ArticleGoogle Scholar
- Andersson FI, Blakytny R, Kirstein J, Turgay K, Bukau B, Mogk A, Clarke AK: Cyanobacterial ClpC/HSP100 protein displays intrinsic chaperone activity. J Biol Chem. 2006, 281: 5468-5475.PubMedView ArticleGoogle Scholar
- Kress W, Maglica Z, Weber-Ban E: Clp chaperone-proteases: structure and function. Res Microbiol. 2009, 160: 618-628. 10.1016/j.resmic.2009.08.006.PubMedView ArticleGoogle Scholar
- Wickner S, Gottesman S, Skowyra D, Hoskins J, McKenney K, Maurizi MR: A molecular chaperone, ClpA, functions like DnaK and DnaJ. Proc Natl Acad Sci USA. 1994, 91: 12218-12222. 10.1073/pnas.91.25.12218.PubMed CentralPubMedView ArticleGoogle Scholar
- Weibezahn J, Schlieker C, Bukau B, Mogk A: Characterization of a trap mutant of the AAA + chaperone ClpB. J Biol Chem. 2003, 278: 32608-32617. 10.1074/jbc.M303653200.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.