Novel lineages of Prochlorococcus and Synechococcus in the global oceans

Picocyanobacteria represented by Prochlorococcus and Synechococcus have an important role in oceanic carbon fixation and nutrient cycling. In this study, we compared the community composition of picocyanobacteria from diverse marine ecosystems ranging from estuary to open oceans, tropical to polar oceans and surface to deep water, based on the sequences of 16S-23S rRNA internal transcribed spacer (ITS). A total of 1339 ITS sequences recovered from 20 samples unveiled diverse and several previously unknown clades of Prochlorococcus and Synechococcus. Six high-light (HL)-adapted Prochlorococcus clades were identified, among which clade HLVI had not been described previously. Prochlorococcus clades HLIII, HLIV and HLV, detected in the Equatorial Pacific samples, could be related to the HNLC clades recently found in the high-nutrient, low-chlorophyll (HNLC), iron-depleted tropical oceans. At least four novel Synechococcus clades (out of six clades in total) in subcluster 5.3 were found in subtropical open oceans and the South China Sea. A niche partitioning with depth was observed in the Synechococcus subcluster 5.3. Members of Synechococcus subcluster 5.2 were dominant in the high-latitude waters (northern Bering Sea and Chukchi Sea), suggesting a possible cold-adaptation of some marine Synechococcus in this subcluster. A distinct shift of the picocyanobacterial community was observed from the Bering Sea to the Chukchi Sea, which reflected the change of water temperature. Our study demonstrates that oceanic systems contain a large pool of diverse picocyanobacteria, and further suggest that new genotypes or ecotypes of picocyanobacteria will continue to emerge, as microbial consortia are explored with advanced sequencing technology.     

Acetate uptake by the unicellular cyanobacteria Synechococcus and Aphanocapsa

Acetate uptake by strains of Synechococcus and Aphanocapsa in short experiments required light, and was strongly inhibited by m-dichlorocarbonyl cyanide phenylhydrazone and dichlorophenyl dimethyl urea. Acetate carbon was distributed in amino acids and in the acyl portion of lipids in the same way as during growth experiments when CO₂ was available, but the reduced incorporation in the absence of CO₂ was primarily into the lipid fraction. An apparent K _m for uptake by Synechococcus and for Aphanocapsa 6308 of 20 and 180 M at pH 7.4 was obtained; corresponding V _max values were 6 and 11 nmol x min^-1 x mg protein^-1. Uptake with Synechococcus was affected by pH, with affinity decreased and maximal rate increase with rising pH. Acetate uptake was not affected by propionate or butyrate when both were added at the same time, but a light and concentration dependent inhibition developed if suspensions were preincubated with propionate. Acetate carbon moved rapidly into acid insoluble material, but after 10–15 s 75% or more of the recovered intracellular counts were in acetyl CoA. Counts in this compound were reduced by preincubation with propionate.Kinetic measurements of acetyl CoA synthetase in fractionated cell extracts gave values for K _m of about 50 M for acetate, 5 mM for propionate, 100 M for CoA and 0.38 mM for ATP. The internal pool of free CoA was measured to be about 20 M, and was reduced by preincubation with propionate. This suggests that the activity of CoA-mediated reactions may be regulated by the availability of this cofactor.Abbreviations Used CCCP m-Dichlorocarbonyl cyanide phenyl hydrazone - DCMU dichlorophenyl dimethyl urea - TCA trichloroacetic acid - Tris trishydroxymethyl amino methane - HEPES N-2-hydroxyethylpiperazine-N-2-ethane-sulfonic acid     

Codon usage patterns among genes for lepidopteran hemolymph proteins

Patterns in codon usage were examined for the coding regions of the 23 known lepidopteran hemolymph proteins. Coding triplets are GC rich at the third position and a significant linear relationship between GC content of silent and nonsilent (replacement) sites was demonstrated. Intron GC content was significantly lower than in coding regions and no relationship between intron GC content and the same at silent and nonsilent sites was found. Though hemolymph proteins are all produced by the same tissue—fat body—significantly less bias was observed when all moth sequences were pooled than when sequences of the two major species were analyzed separately, as predicted by the genome hypothesis. In cases where no statistically significant bias was observed, polar or acidic basic amino acids were almost exclusively involved. Calculation of codon adaptation indices (CAI) was of limited value in quantifying the degree of codon bias and probably reflects the complexity of multicellular-organism life cycles and the changing patterns of gene expression over different developmental stages. Correspondence to: D.R. Frohlich     

Codon usage patterns suggest independent evolution of two catabolic operons on toluene-degradative plasmid TOL pWW0 of Pseudomonas putida

TOL plasmid pWW0 of Pseudomonas putida encodes a set of enzymes responsible for the degradation of toluene. The structural genes for these catobolic enzymes are clustered into two operons—namely, the xylCMAB and xylXYZLTEGFJQKIH operons. We examined the codon usage patterns of these catabolic genes by measuring the codon-usage distances between pairs of these catabolic genes. The codon-usage distance, d, between gene 1 and gene 2 was defined as d = [(p _j – q _j )²]^1/2, where p _j > and q _j are the frequencies of the j-th codon in gene 1 and 2, respectively, j being any one of the 64 possible codons. We found that the genes in the same operon exhibit similar codon-usage patterns while genes in the different operons exhibit different codon bias. This observation suggests that genes in the same operon have coevolved, and that the ancestors of the xylCMAB and xylXYZLTEGFJQKIH operons evolved in different organisms. Correspondence to: S. Harayama     

Codon usage and base composition inRickettsia prowazekii

Codon usage and base composition in sequences from the A + T-rich genome ofRickettsia prowazekii, a member of the alpha Proteobacteria, have been investigated. Synonymous codon usage patterns are roughly similar among genes, even though the data set includes genes expected to be expressed at very different levels, indicating that translational selection has been ineffective in this species. However, multivariate statistical analysis differentiates genes according to their G + C contents at the first two codon positions. To study this variation, we have compared the amino acid composition patterns of 21R. prowazekii proteins with that of a homologous set of proteins fromEscherichia coli. The analysis shows that individual genes have been affected by biased mutation rates to very different extents: genes encoding proteins highly conserved among other species being the least affected. Overall, protein coding and intergenic spacer regions have G + C content values of 32.5% and 21.4%, respectively. Extrapolation from these values suggests thatR. prowazekii has around 800 genes and that 60–70% of the genome may be coding. Correspondence to: S.G.E. Andersson     

Evidence of selectively driven codon usage in rice: implications for GC content evolution of Gramineae genes

Two gene classes characterized by high and low GC content have been found in rice and other cereals, but not dicot genomes. We used paralogs with high and low GC contents in rice and found: (a) a greater increase in GC content at exonic fourfold-redundant sites than at flanking introns; (b) with reference to their orthologs in Arabidopsis, most substitution sites between the two kinds of paralogs are found at 2- and 4-degenerate sites with a T-->C mode, while A-->C and A-->G play major roles at 0-degenerate sites; and (c) high-GC genes have greater bias and codon usage is skewed toward codons that are preferred in highly expressed genes. We believe this is strong evidence for selectively driven codon usage in rice. Another cereal, maize, also showed the same trend as in rice. This represents a potential evolutionary process for the origin of genes with a high GC content in rice and other cereals.     

Analysis of codon usage pattern in the radioresistant bacterium Deinococcus radiodurans

The main factors shaping codon usage bias in the Deinococcus radiodurans genome were reported. Correspondence analysis (COA) was carried out to analyze synonymous codon usage bias. The results showed that the main trend was strongly correlated with gene expression level assessed by the "Codon Adaptation Index" (CAI) values, a result that was confirmed by the distribution of genes along the first axis. The results of correlation analysis, variance analysis and neutrality plot indicated that gene nucleotide composition was clearly contributed to codon bias. CDS length was also key factor in dictating codon usage variation. A general tendency of more biased codon usage of genes with longer CDS length to higher expression level was found. Further, the hydrophobicity of each protein also played a role in shaping codon usage in this organism, which could be confirmed by the significant correlation between the positions of genes placed on the first axis and the hydrophobicity values (r=-0.100, P<0.01). In summary, gene expression level played a crucial role, nucleotide mutational bias, CDS length and the hydrophobicity of each protein just in a minor way in shaping the codon usage pattern of D. radiodurans. Notably, 19 codons firstly defined as "optimal codons" may provide useful clues for molecular genetic engineering and evolutionary studying.     

Rapid evolutionary divergence of Photosystem I core subunits PsaA and PsaB in the marine prokaryote Prochlorococcus

The nucleotide sequences of the genes coding for the subunits of the Photosystem I (PS I) core, PsaA and PsaB were determined for the marine prokaryotic oxyphototrophs Prochlorococcus sp. MED4 (CCMP1378), P. marinus SS120 (CCMP1375) and Synechococcus sp. WH7803. Divergence of these sequences from those of both freshwater cyanobacteria and higher plants was remarkably high, given the conserved nature of PsaA and PsaB proteins. In particular, the PsaA of marine prokaryotes showed several specific insertions and deletions with regard to known PsaA sequences. Even in between the two Prochlorococcus strains, which correspond to two genetically different ecotypes with shifted growth irradiance optima, the sequence identity was only 80.2% for PsaA and 88.9% for PsaB. Possible causes and implications of the fast evolution rates of these two PS I core subunits are discussed. This revised version was published online in June 2006 with corrections to the Cover Date.     

Codon usage divergence of homologous vertebrate genes and codon usage clock

Summary This paper is concerned with the divergence of synonymous codon usage and its bias in three homologous genes within vertebrate species. Genetic distances among species are described in terms of synonymous codon usage divergence and the correlation is found between the genetic distances and taxonomic distances among species under study. A codon usage clock is reported in alphaglobin and beta-globin. A method is developed to define the synonymous codon preference bias and it is observed that the bias changes considerably among species.     

Potential for phosphite and phosphonate utilization by Prochlorococcus

Phosphonates (Pn) are diverse organic phosphorus (P) compounds containing C–P bonds and comprise up to 25% of the high-molecular weight dissolved organic P pool in the open ocean. Pn bioavailability was suggested to influence markedly bacterial primary production in low-P areas. Using metagenomic data from the Global Ocean Sampling expedition, we show that the main potential microbial contributor in Pn utilization in oceanic surface water is the globally important marine primary producer Prochlorococcus. Moreover, a number of Prochlorococcus strains contain two distinct putative Pn uptake operons coding for ABC-type Pn transporters. On the basis of microcalorimetric measurements, we find that each of the two different putative Pn-binding protein (PhnD) homologs transcribed from these operons possesses different Pn- as well as inorganic phosphite-binding specificities. Our results suggest that Prochlorococcus adapt to low-P environments by increasing the number of Pn transporters with different specificities towards phosphite and different Pns.     

Ecology of uncultured Prochlorococcus clades revealed through single-cell genomics and biogeographic analysis

Prochlorococcus is the numerically dominant photosynthetic organism throughout much of the world''s oceans, yet little is known about the ecology and genetic diversity of populations inhabiting tropical waters. To help close this gap, we examined natural Prochlorococcus communities in the tropical Pacific Ocean using a single-cell whole-genome amplification and sequencing. Analysis of the gene content of just 10 single cells from these waters added 394 new genes to the Prochlorococcus pan-genome—that is, genes never before seen in a Prochlorococcus cell. Analysis of marker genes, including the ribosomal internal transcribed sequence, from dozens of individual cells revealed several representatives from two uncultivated clades of Prochlorococcus previously identified as HNLC1 and HNLC2. While the HNLC clades can dominate Prochlorococcus communities under certain conditions, their overall geographic distribution was highly restricted compared with other clades of Prochlorococcus. In the Atlantic and Pacific oceans, these clades were only found in warm waters with low Fe and high inorganic P levels. Genomic analysis suggests that at least one of these clades thrives in low Fe environments by scavenging organic-bound Fe, a process previously unknown in Prochlorococcus. Furthermore, the capacity to utilize organic-bound Fe appears to have been acquired horizontally and may be exchanged among other clades of Prochlorococcus. Finally, one of the single Prochlorococcus cells sequenced contained a partial genome of what appears to be a prophage integrated into the genome.     

Rapid extraction of phycobiliproteins from cultured cyanobacteria samples

Cyanobacteria are a valuable and ubiquitous component of marine picophytoplankton that contribute significantly to total carbon biomass and primary productivity of the oceans. They contain water soluble, natively highly fluorescent proteins, phycobiliproteins, that can be considered ideal marker pigments for understanding the distribution and trophic dynamics of picoplankton populations. However, there is no standard protocol for extracting and quantitating these proteins from cyanobacterial cells. Ideally, the cells would be disrupted quickly and efficiently with complete extraction and recovery of the released proteins. For that purpose, we describe a method for extracting phycobiliproteins from a Synechococcus CCMP 833 cyanobacteria culture that utilizes 3% 3-[(3-cholamidopropyl)dimethyammonio]propanesulfonic acid (Chaps) 0.3% asolectin combined with nitrogen cavitation. Extraction efficiencies of greater than 85% were achieved by this method, which requires less than 3h. The analysis of the extracted samples was carried out by capillary electrophoresis with laser-induced fluorescence detection.     

Non-uniqueness of factors constraint on the codon usage in Bombyx mori

Background

The analysis of codon usage is a good way to understand the genetic and evolutionary characteristics of an organism. However, there are only a few reports related with the codon usage of the domesticated silkworm, Bombyx mori (B. mori). Hence, the codon usage of B. mori was analyzed here to reveal the constraint factors and it could be helpful to improve the bioreactor based on B. mori.

Results

A total of 1,097 annotated mRNA sequences from B. mori were analyzed, revealing there is only a weak codon bias. It also shows that the gene expression level is related to the GC content, and the amino acids with higher general average hydropathicity (GRAVY) and aromaticity (Aromo). And the genes on the primary axis are strongly positively correlated with the GC content, and GC3s. Meanwhile, the effective number of codons (ENc) is strongly correlated with codon adaptation index (CAI), gene length, and Aromo values. However, the ENc values are correlated with the second axis, which indicates that the codon usage in B. mori is affected by not only mutation pressure and natural selection, but also nucleotide composition and the gene expression level. It is also associated with Aromo values, and gene length. Additionally, B. mori has a greater relative discrepancy in codon preferences with Drosophila melanogaster (D. melanogaster) or Saccharomyces cerevisiae (S. cerevisiae) than with Arabidopsis thaliana (A. thaliana), Escherichia coli (E. coli), or Caenorhabditis elegans (C. elegans).

Conclusions

The codon usage bias in B. mori is relatively weak, and many influence factors are found here, such as nucleotide composition, mutation pressure, natural selection, and expression level. Additionally, it is also associated with Aromo values, and gene length. Among them, natural selection might play a major role. Moreover, the “optimal codons” of B. mori are all encoded by G and C, which provides useful information for enhancing the gene expression in B. mori through codon optimization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1596-z) contains supplementary material, which is available to authorized users.     

Codon bias and gene expression of mitochondrial ND2 gene in chordates

Background: Mitochondrial ND gene, which encodes NADH dehydrogenase, is the first enzyme of the mitochondrial electron transport chain. Leigh syndrome, a neurodegenerative disease caused by mutation in the ND2 gene (T4681C), is associated with bilateral symmetric lesions in basal ganglia and subcortical brain regions. Therefore, it is of interest to analyze mitochondrial DNA to glean information for evolutionary relationship. This study highlights on the analysis of compositional dynamics and selection pressure in shaping the codon usage patterns in the coding sequence of MT-ND2 gene across pisces, aves and mammals by using bioinformatics tools like effective number of codons (ENC), codon adaptation index (CAI), relative synonymous codon usage (RSCU) etc. Results: We observed a low codon usage bias as reflected by high ENC values in MT-ND2 gene among pisces, aves and mammals. The most frequently used codons were ending with A/C at the 3^rd position of codon and the gene was AT rich in all the three classes. The codons TCA, CTA, CGA and TGA were over represented in all three classes. The F1 correspondence showed significant positive correlation with G, T3 and CAI while the F2 axis showed significant negative correlation with A and T but significant positive correlation with G, C, G3, C3, ENC, GC, GC1, GC2 and GC3. Conclusions: The codon usage bias in MTND2 gene is not associated with expression level. Mutation pressure and natural selection affect the codon usage pattern in MT-ND 2 gene.     

A comparison of the major lipid classes and fatty acid composition of marine unicellular cyanobacteria with freshwater species

The fatty acid composition of two motile (strains WH 8113 and WH 8103) and one nonmotile (strain WH 7803) marine cyanobacteria has been determined and compared with two freshwater unicellular Synechocystis species (strain PCC 6308 and PCC 6803). The fatty acid composition of lipid extracts of isolated membranes from Synechocystis PCC 6803 was found to be identical to that of whole cells. All the marine strains contained myristic acid (14:0) as the major fatty acid, with only traces of polyunsaturated fatty acids. This composition is similar to Synechocystis PCC 6308. The major lipid classes of the nonmotile marine strain were identified as digalactosyl diacylglycerol, monogalactosyl diacylglycerol, phosphatidylglycerol, and sulfoquinovosyl diacylglycerol, identical to those found in other cyanobacteria.Abbreviations DGDG Digalactosyl diacylglycerol - MGDG Monogalactosyldiacylglycerol - PG Phosphatidylglycerol - SGDG sulfoquinovosyl diacylglycerol - gc gas chromatography - ms mass spectrometry     

Synonymous codon usage in Thermosynechococcus elongatus (cyanobacteria) identifies the factors shaping codon usage variation

Analysis of synonymous codon usage pattern in the genome of a thermophilic cyanobacterium, Thermosynechococcus elongatus BP-1 using multivariate statistical analysis revealed a single major explanatory axis accounting for codon usage variation in the organism. This axis is correlated with the GC content at third base of synonymous codons (GC3s) in correspondence analysis taking T. elongatus genes. A negative correlation was observed between effective number of codons i.e. Nc and GC3s. Results suggested a mutational bias as the major factor in shaping codon usage in this cyanobacterium. In comparison to the lowly expressed genes, highly expressed genes of this organism possess significantly higher proportion of pyrimidine-ending codons suggesting that besides, mutational bias, translational selection also influenced codon usage variation in T. elongatus. Correspondence analysis of relative synonymous codon usage (RSCU) with A, T, G, C at third positions (A3s, T3s, G3s, C3s, respectively) also supported this fact and expression levels of genes and gene length also influenced codon usage. A role of translational accuracy was identified in dictating the codon usage variation of this genome. Results indicated that although mutational bias is the major factor in shaping codon usage in T. elongatus, factors like translational selection, translational accuracy and gene expression level also influenced codon usage variation.