Characterizing the Complete Mitochondrial Genome of Arma custos and Picromerus lewisi (Hemiptera: Pentatomidae: Asopinae) and Conducting Phylogenetic Analysis.

We characterized the mitochondrial genome (mitogenome) and conducted phylogenetic analyses of 48 Hemiptera species by sequencing and analyzing the mitogenome of Arma custos (Fabricius) and Picromerus lewisi (Scott). The complete mitogenomes of the two predators were 16,024 bp and 19,587 bp in length, respectively, and it contained 37 classical genes, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs), and a control region. Most PCGs in these predators use ATN as the start codon. This research revealed that the genes of the two natural enemy species have an A + T content of 75.40% and all tRNAs have a typical cloverleaf structure, with the exception of trnS1, which lacks a dihydrouridine arm. This is the first study to compare the mitochondrial genetic structure of two predatory insects; the mitochondrial genetic structure of individual predatory insects has been sequenced in previous studies. Here, phylogenetic analysis on the basis of amino acid and nucleotide sequences of 13 mitochondrial PCGs using Bayesian inference and maximum likelihood methods were conducted to generate similar tree topologies, which suggested that the two predators with close genetic relationships belong to Asopinae subfamily. Furthermore, the monophyly of the Pentatomoidea superfamily is well accepted despite limited taxon and species sampling. Finally, their complete mitogenome provided data to establish a predator-prey food web, which is the foundation of effective pest management. Our results also enhanced the database of natural enemy insects.

Sequencing and annotation of the endangered wild buffalo (Bubalus arnee) mitogenome for taxonomic assessment.

The wild water buffalo (Bubalus arnee) is one of the most endangered and least studied large bovid in the Indian subcontinent. India retains 90% of the estimated global population of >4000 individuals as two fragmented populations in Assam and Chhattisgarh, both threatened by habitat loss and degradation, hunting, disease from livestock, and hybridization with the domestic buffalos. Small, fragmented population size and potential hybridisation pressures from co-occurring domestic buffalos are the major conservation challenges. For the first time, we sequenced the 16,357 bp long mitogenome of three opportunistically collected wild water buffalo samples from Assam (n = 1) and Chhattishgarh (n = 2). The annotated sequence has a base composition of 26.4% T, 26.6% C, 33.1% A and 13.9% G depicting an AT-rich mitogenome composition, including 13 protein-coding genes (11,361 bp), 22 transfer RNA (tRNA) (1514 bp), two ribosomal genes (2525 bp), and a non-coding control region (928 bp). The gene order is conserved with other bovid species. Comparative mitogenome analyses showed both populations are genetically similar but significantly different from domestic buffalo. We also identified structural differences in seven tRNA secondary structures between both species. The genetic distance between wild buffalo and other bovids varied between 0.103 and 0.122. Multiple Bayesian phylogenetic trees showed that both wild and domestic water buffalo formed sister clades which were paraphyletic to other potentially sympatric species of genus Bos. This study provides baseline information on wild buffalo mitogenome for further research on phylogeny, phylogeography and hybrid assessment and help conserving this endangered species.

Hierarchical natural move Monte Carlo refines flexible RNA structures into cryo-EM densities.

Ribonucleic acids (RNAs) play essential roles in living cells. Many of them fold into defined three-dimensional (3D) structures to perform functions. Recent advances in single-particle cryo-electron microscopy (cryo-EM) have enabled structure determinations of RNA to atomic resolutions. However, most RNA molecules are structurally flexible, limiting the resolution of their structures solved by cryo-EM. In modeling these molecules, several computational methods are limited by the requirement of massive computational resources and/or the low efficiency in exploring large-scale structural variations. Here we use hierarchical natural move Monte Carlo (HNMMC), which takes advantage of collective motions for groups of nucleic acid residues, to refine RNA structures into their cryo-EM maps, preserving atomic details in the models. After validating the method on a simulated density map of tRNA, we applied it to objectively obtain the model of the folding intermediate for the specificity domain of ribonuclease P from Bacillus subtilis and refine a flexible ribosomal RNA (rRNA) expansion segment from the Mycobacterium tuberculosis (Mtb) ribosome in different conformational states. Finally, we used HNMMC to model atomic details and flexibility for two distinct conformations of the complete genomic RNA (gRNA) inside MS2, a single-stranded RNA virus, revealing multiple pathways for its capsid assembly.

tRNA Fragments Show Intertwining with mRNAs of Specific Repeat Content and Have Links to Disparities.

tRNA-derived fragments (tRF) are a class of potent regulatory RNAs. We mined the datasets from The Cancer Genome Atlas (TCGA) representing 32 cancer types with a deterministic and exhaustive pipeline for tRNA fragments. We found that mitochondrial tRNAs contribute disproportionally more tRFs than nuclear tRNAs. Through integrative analyses, we uncovered a multitude of statistically significant and context-dependent associations between the identified tRFs and mRNAs. In many of the 32 cancer types, these associations involve mRNAs from developmental processes, receptor tyrosine kinase signaling, the proteasome, and metabolic pathways that include glycolysis, oxidative phosphorylation, and ATP synthesis. Even though the pathways are common to multiple cancers, the association of specific mRNAs with tRFs depends on and differs from cancer to cancer. The associations between tRFs and mRNAs extend to genomic properties as well; specifically, tRFs are positively correlated with shorter genes that have a higher density in repeats, such as ALUs, MIRs, and ERVLs. Conversely, tRFs are negatively correlated with longer genes that have a lower repeat density, suggesting a possible dichotomy between cell proliferation and differentiation. Analyses of bladder, lung, and kidney cancer data indicate that the tRF-mRNA wiring can also depend on a patient's sex. Sex-dependent associations involve cyclin-dependent kinases in bladder cancer, the MAPK signaling pathway in lung cancer, and purine metabolism in kidney cancer. Taken together, these findings suggest diverse and wide-ranging roles for tRFs and highlight the extensive interconnections of tRFs with key cellular processes and human genomic architecture. SIGNIFICANCE: Across 32 TCGA cancer contexts, nuclear and mitochondrial tRNA fragments exhibit associations with mRNAs that belong to concrete pathways, encode proteins with particular destinations, have a biased repeat content, and are sex dependent.

tiRNAs as a novel biomarker for cell damage assessment in in vitro ischemia-reperfusion model in rat neuronal PC12 cells.

The disruption of appropriate cellular stress responses is implicated in the pathogenesis of different neurological disorders including ischemic injury. Early diagnosis and treatment are often associated with better prognosis in ischemic stroke patients. Thus, there is an urgent need to improve the speed and accuracy of stroke diagnosis by developing highly sensitive stroke biomarkers. We recently reported that transfer RNA (tRNA) was involved in cell stress response pathways. Under cell stress conditions, mature tRNA is cleaved by a specific ribonuclease, angiogenin, generating tRNA-derived stress-induced RNA (tiRNA). To study tiRNA generation in an in vitro model of ischemic-reperfusion injury, we used the rat neuronal cell line, PC12, in combination with analysis of SYBR staining and immuno-northern blotting using anti-1-methyladenosine antibody, which detects 1-methyladenosine (m1A) modification of tRNA. We demonstrated that oxygen-glucose deprivation induced tRNA cleavage and tiRNA generation. Time course analysis showed a dramatic up-regulation of tiRNA generation by oxygen-glucose deprivation (OGD) which started a few minutes after reperfusion. Minocycline, a neuroprotective antibiotic, treatment protected PC12 cells against OGD-reperfusion cell damage resulting in a marked down-regulation of the generated tiRNA. Our findings show that cleavage of tRNA and tiRNA generation in rat neuronal PC12 cells occurs with reperfusion injury and the detection of tiRNA could be used as a potential cell damage marker and treatment effect indicator for this type of injury.

Complete mitochondrial genome of northern Indian red muntjac (Muntiacus vaginalis) and its phylogenetic analysis.

We report complete mitochondrial genome of Northern Indian red muntjac, Muntiacus vaginalis, and its phylogenetic inferences. Mitogenome composition was 16,352 bp in length and its overall base composition in the circular genome was A = 33.2%, T = 29.0%, C = 24.50% and G = 13.30%. It exhibited a typical mitogenome structure, including 22 transfer RNA genes, 13 protein-coding genes, two ribosomal RNA genes and a major non-coding control region (D-loop region). All the genes except ND6 and eight tRNA's were encoded on the heavy strand. Phylogenetic analyses showed that M. vaginalis is closely related to M. muntjak and formed a sister relationship with Elaphodus cephalophus. In view of the unclear distribution range and escalating habitat loss, it is important to identify its population genetic status. The complete mitogenome described in this study can be used in further phylogenetics, identification of extant maternal lineage, evolutionary significance unit and its genetic conservation.

Selection-driven cost-efficiency optimization of transcripts modulates gene evolutionary rate in bacteria.

BACKGROUND: Most amino acids are encoded by multiple synonymous codons. However, synonymous codons are not used equally, and this biased codon use varies between different organisms. It has previously been shown that both selection acting to increase codon translational efficiency and selection acting to decrease codon biosynthetic cost contribute to differences in codon bias. However, it is unknown how these two factors interact or how they affect molecular sequence evolution. RESULTS: Through analysis of 1320 bacterial genomes, we show that bacterial genes are subject to multi-objective selection-driven optimization of codon use. Here, selection acts to simultaneously decrease transcript biosynthetic cost and increase transcript translational efficiency, with highly expressed genes under the greatest selection. This optimization is not simply a consequence of the more translationally efficient codons being less expensive to synthesize. Instead, we show that transfer RNA gene copy number alters the cost-efficiency trade-off of synonymous codons such that, for many species, selection acting on transcript biosynthetic cost and translational efficiency act in opposition. Finally, we show that genes highly optimized to reduce cost and increase efficiency show reduced rates of synonymous and non-synonymous mutation. CONCLUSIONS: This analysis provides a simple mechanistic explanation for variation in evolutionary rate between genes that depends on selection-driven cost-efficiency optimization of the transcript. These findings reveal how optimization of resource allocation to messenger RNA synthesis is a critical factor that determines both the evolution and composition of genes.

Race Disparities in the Contribution of miRNA Isoforms and tRNA-Derived Fragments to Triple-Negative Breast Cancer.

Triple-negative breast cancer (TNBC) is a breast cancer subtype characterized by marked differences between White and Black/African-American women. We performed a systems-level analysis on datasets from The Cancer Genome Atlas to elucidate how the expression patterns of mRNAs are shaped by regulatory noncoding RNAs (ncRNA). Specifically, we studied isomiRs, that is, isoforms of miRNAs, and tRNA-derived fragments (tRF). In normal breast tissue, we observed a marked cohesiveness in both the ncRNA and mRNA layers and the associations between them. This cohesiveness was widely disrupted in TNBC. Many mRNAs become either differentially expressed or differentially wired between normal breast and TNBC in tandem with isomiR or tRF dysregulation. The affected pathways included energy metabolism, cell signaling, and immune responses. Within TNBC, the wiring of the affected pathways with isomiRs and tRFs differed in each race. Multiple isomiRs and tRFs arising from specific miRNA loci (e.g., miR-200c, miR-21, the miR-17/92 cluster, the miR-183/96/182 cluster) and from specific tRNA loci (e.g., the nuclear tRNAGly and tRNALeu, the mitochondrial tRNAVal and tRNAPro) were strongly associated with the observed race disparities in TNBC. We highlight the race-specific aspects of transcriptome wiring by discussing in detail the metastasis-related MAPK and the Wnt/ß-catenin signaling pathways, two of the many key pathways that were found differentially wired. In conclusion, by employing a data- and knowledge-driven approach, we comprehensively analyzed the normal and cancer transcriptomes to uncover novel key contributors to the race-based disparities of TNBC.Significance: This big data-driven study comparing normal and cancer transcriptomes uncovers RNA expression differences between Caucasian and African-American patients with triple-negative breast cancer that might help explain disparities in incidence and aggressive character. Cancer Res; 78(5); 1140-54. ©2017 AACR.

Assessment of mitochondrial DNA mutations in Chinese family with essential hypertension.

Mitochondrial DNA (MtDNA) mutations played an important role in the development of essential hypertension. Mitochondrial tRNA point mutations, caused the failure in tRNA metabolism, responsible for the pathogenesis of this complex disease. In this study, we evaluated the possible role of the 4329C >G mutation in the clinical expression of hypertension in a Chinese family. Analysis of the complete mtDNA sequence variants showed that other mutations may play synergic roles in the phenotypic manifestation of hypertension. In addition, other potential pitfalls were also discussed in this context.

Noncanonical secondary structure stabilizes mitochondrial tRNA(Ser(UCN)) by reducing the entropic cost of tertiary folding.

Mammalian mitochondrial tRNA(Ser(UCN)) (mt-tRNA(Ser)) and pyrrolysine tRNA (tRNA(Pyl)) fold to near-canonical three-dimensional structures despite having noncanonical secondary structures with shortened interhelical loops that disrupt the conserved tRNA tertiary interaction network. How these noncanonical tRNAs compensate for their loss of tertiary interactions remains unclear. Furthermore, in human mt-tRNA(Ser), lengthening the variable loop by the 7472insC mutation reduces mt-tRNA(Ser) concentration in vivo through poorly understood mechanisms and is strongly associated with diseases such as deafness and epilepsy. Using simulations of the TOPRNA coarse-grained model, we show that increased topological constraints encoded by the unique secondary structure of wild-type mt-tRNA(Ser) decrease the entropic cost of folding by ∼2.5 kcal/mol compared to canonical tRNA, offsetting its loss of tertiary interactions. Further simulations show that the pathogenic 7472insC mutation disrupts topological constraints and hence destabilizes the mutant mt-tRNA(Ser) by ∼0.6 kcal/mol relative to wild-type. UV melting experiments confirm that insertion mutations lower mt-tRNA(Ser) melting temperature by 6-9 °C and increase the folding free energy by 0.8-1.7 kcal/mol in a largely sequence- and salt-independent manner, in quantitative agreement with our simulation predictions. Our results show that topological constraints provide a quantitative framework for describing key aspects of RNA folding behavior and also provide the first evidence of a pathogenic mutation that is due to disruption of topological constraints.

Kinetic Monte Carlo approach to RNA folding dynamics using structure-based models.

RNA molecules form three-dimensional structures via base pairing that determine the function and biochemical activity of the molecule. Here we introduce a structure-based method for studying the folding dynamics of RNA secondary structures. The approach focuses on native contacts that are parametrized with standard empirical free energies. Kinetic Monte Carlo simulations for free folding of simple hairpins and complex structures such as a tRNA as well as for folding in the presence of an external force show good agreement with experimental data. A systematic comparison of simulated and experimental folding rates for various structures shows a strong correlation, indicating that the approach can predict folding rates within about an order of magnitude.

Decoding properties of tRNA leave a detectable signal in codon usage bias.

MOTIVATION: The standard genetic code translates 61 codons into 20 amino acids using fewer than 61 transfer RNAs (tRNAs). This is possible because of the tRNA's ability to 'wobble' at the third base to decode more than one codon. Although the anticodon-codon mapping of tRNA to mRNA is a prerequisite for certain codon usage indices and can contribute to the understanding of the evolution of alternative genetic codes, it is usually not determined experimentally because such assays are prohibitively expensive and elaborate. Instead, the codon reading is approximated from theoretical inferences of nucleotide binding, the wobble rules. Unfortunately, these rules fail to capture all of the nuances of codon reading. This study addresses the codon reading properties of tRNAs and their evolutionary impact on codon usage bias. RESULTS: Using three different computational methods, the signal of tRNA decoding in codon usage bias is identified. The predictions by the methods generally agree with each other and compare well with experimental evidence of codon reading. This analysis suggests a revised codon reading for cytosolic tRNA in the yeast genome (Saccharomyces cerevisiae) that is more accurate than the common assignment by wobble rules. The results confirm the earlier observation that the wobble rules are not sufficient for a complete description of codon reading, because they depend on genome-specific factors. The computational methods presented here are applicable to any fully sequenced genome. AVAILABILITY: By request from the author. CONTACT: alexander.roth@isb-sib.ch.

Phylogeny and historical biogeography of ancient assassin spiders (Araneae: Archaeidae) in the Australian mesic zone: evidence for Miocene speciation within Tertiary refugia.

The rainforests, wet sclerophyll forests and temperate heathlands of the Australian mesic zone are home to a diverse and highly endemic biota, including numerous old endemic lineages restricted to refugial, mesic biomes. A growing number of phylogeographic studies have attempted to explain the origins and diversification of the Australian mesic zone biota, in order to test and better understand the mode and tempo of historical speciation within Australia. Assassin spiders (family Archaeidae) are a lineage of iconic araneomorph spiders, characterised by their antiquity, remarkable morphology and relictual biogeography on the southern continents. The Australian assassin spider fauna is characterised by a high diversity of allopatric species, many of which are restricted to individual mountains or montane systems, and all of which are closely tied to mesic and/or refugial habitats in the east and extreme south-west of mainland Australia. We tested the phylogeny and vicariant biogeography of the Australian Archaeidae (genus Austrarchaea Forster & Platnick), using a multi-locus molecular approach. Fragments from six mitochondrial genes (COI, COII, tRNA-K, tRNA-D, ATP8, ATP6) and one nuclear protein-coding gene (Histone H3) were used to infer phylogenetic relationships and to explore the phylogeographic origins of the diverse Australian fauna. Bayesian analyses of the complete molecular dataset, along with differentially-partitioned Bayesian and parsimony analyses of a smaller concatenated dataset, revealed the presence of three major Australian lineages, each with non-overlapping distributions in north-eastern Queensland, mid-eastern Australia and southern Australia, respectively. Divergence date estimation using mitochondrial data and a rate-calibrated relaxed molecular clock revealed that major lineages diverged in the early Tertiary period, prior to the final rifting of Australia from East Antarctica. Subsequent speciation occurred during the Miocene (23-5.3 million years ago), with tropical and subtropical taxa diverging in the early-mid Miocene, prior to southern and temperate taxa in the mid-late Miocene. Area cladograms reconciled with Bayesian chronograms for all known Archaeidae in southern and south-eastern Australia revealed seven potentially vicariant biogeographic barriers in eastern Queensland, New South Wales and southern Australia, each proposed and discussed in relation to other mesic zone taxa. Five of these barriers were inferred as being of early Miocene age, and implicated in the initial vicariant separation of endemic regional clades. Phylogeographic results for Australian Archaeidae are congruent with a model of sequential allopatric speciation in Tertiary refugia, as driven by the contraction and fragmentation of Australia's mesic biomes during the Miocene. Assassin spiders clearly offer great potential for further testing historical biogeographic processes in temperate and eastern Australia, and are a useful group for better understanding the biology and biogeography of the Australian mesic zone.

Relative roles of Neogene vicariance and Quaternary climate change on the historical diversification of bunchgrass lizards (Sceloporus scalaris group) in Mexico.

Neogene vicariance during the Miocene and Pliocene and Quaternary climate change have synergistically driven diversification in Mexican highland taxa. We investigated the impacts of these processes on genetic diversification in the widely distributed bunchgrass lizards in the Sceloporus scalaris group. We searched for correlations between timing in diversification and timing of (1) a period of marked volcanism across the Trans-Mexican Volcanic Belt in central Mexico 3-7.5million years ago (Ma) and (2) a transition to larger glacial-interglacial cycles during the mid-Pleistocene. From our phylogenetic analyses of mitochondrial DNA we identified two major clades that contained 13 strongly supported lineages. One clade contained lineages from the two northern sierras of Mexico, and the other clade included lineages associated with the Trans-Mexican Volcanic Belt and Central Mexican Plateau. Results provided support for Neogene divergences within the S. scalaris group in response to uplift of the Trans-Mexican Volcanic Belt, a pattern observed in several co-distributed taxa, and suggested that Quaternary climate change likely had little effect on diversification between lineages. Uplift of the Trans-Mexican Volcanic Belt during specific time periods appears to have strongly impacted diversification in Mexican highland taxa.

The dynamics of supply and demand in mRNA translation.

We study the elongation stage of mRNA translation in eukaryotes and find that, in contrast to the assumptions of previous models, both the supply and the demand for tRNA resources are important for determining elongation rates. We find that increasing the initiation rate of translation can lead to the depletion of some species of aa-tRNA, which in turn can lead to slow codons and queueing. Particularly striking "competition" effects are observed in simulations of multiple species of mRNA which are reliant on the same pool of tRNA resources. These simulations are based on a recent model of elongation which we use to study the translation of mRNA sequences from the Saccharomyces cerevisiae genome. This model includes the dynamics of the use and recharging of amino acid tRNA complexes, and we show via Monte Carlo simulation that this has a dramatic effect on the protein production behaviour of the system.

Interspecific and intragenic differences in codon usage bias among vertebrate myosin heavy-chain genes.

Synonymous codon usage bias is a broadly observed phenomenon in bacteria, plants, and invertebrates and may result from selection. However, the role of selective pressures in shaping codon bias is still controversial in vertebrates, particularly for mammals. The myosin heavy-chain (MyHC) gene family comprises multiple isoforms of the major force-producing contractile protein in cardiac and skeletal muscles. Slow and fast genes are tandemly arrayed on separate chromosomes, and have distinct patterns of functionality and expression in muscle. We analyze both full-length MyHC genes (~5400 bp) and a larger collection of partial sequences at the 3' end (~500 bp). The MyHC isoforms are an interesting system in which to study codon usage bias because of their length, expression, and critical importance to organismal mobility. Codon bias and GC content differs among MyHC genes with regards to functional type, isoform, and position within the gene. Codon bias even varies by isoform within a species. We find evidence in favor of both chromosomal influences on nucleotide composition and selection against nonsense errors (SANE) acting on codon usage in MyHC genes. Intragenic variation in codon bias and elongation rate is significant, with a strong trend for increasing codon bias and elongation rate towards the 3' end of the gene, although the trend is dependent upon the degeneracy class of the codons. Therefore, patterns of codon usage in MyHC genes are consistent with models supporting SANE as a major force shaping codon usage.

Evolutionary systems biology of amino acid biosynthetic cost in yeast.

Every protein has a biosynthetic cost to the cell based on the synthesis of its constituent amino acids. In order to optimise growth and reproduction, natural selection is expected, where possible, to favour the use of proteins whose constituents are cheaper to produce, as reduced biosynthetic cost may confer a fitness advantage to the organism. Quantifying the cost of amino acid biosynthesis presents challenges, since energetic requirements may change across different cellular and environmental conditions. We developed a systems biology approach to estimate the cost of amino acid synthesis based on genome-scale metabolic models and investigated the effects of the cost of amino acid synthesis on Saccharomyces cerevisiae gene expression and protein evolution. First, we used our two new and six previously reported measures of amino acid cost in conjunction with codon usage bias, tRNA gene number and atomic composition to identify which of these factors best predict transcript and protein levels. Second, we compared amino acid cost with rates of amino acid substitution across four species in the genus Saccharomyces. Regardless of which cost measure is used, amino acid biosynthetic cost is weakly associated with transcript and protein levels. In contrast, we find that biosynthetic cost and amino acid substitution rates show a negative correlation, but for only a subset of cost measures. In the economy of the yeast cell, we find that the cost of amino acid synthesis plays a limited role in shaping transcript and protein expression levels compared to that of translational optimisation. Biosynthetic cost does, however, appear to affect rates of amino acid evolution in Saccharomyces, suggesting that expensive amino acids may only be used when they have specific structural or functional roles in protein sequences. However, as there appears to be no single currency to compute the cost of amino acid synthesis across all cellular and environmental conditions, we conclude that a systems approach is necessary to unravel the full effects of amino acid biosynthetic cost in complex biological systems.

Limited resources in a driven diffusion process.

The advance of particles in many driven diffusion systems depends on the availability of resources in the surrounding environment. In the balance between supply and demand of such resources we are confronted with a regime in which, under limited resource availability, the flow is markedly reduced. In the context of mRNA translation this represents the finite availability of amino acid-tRNA molecules. In this limited resources regime a severe depletion of amino acid tRNAs is also observed. These dramatic effects are vital to our understanding of translation, and are likely to also be important for the many other applications of driven diffusion models.

Unique cost dynamics elucidate the role of frameshifting errors in promoting translational robustness.

There is now considerable evidence supporting the view that codon usage is frequently under selection for translational accuracy. There are, however, multiple forms of inaccuracy (missense, premature termination, and frameshifting errors) and pinpointing a particular error process behind apparently adaptive mRNA anatomy is rarely straightforward. Understanding differences in the fitness costs associated with different types of translational error can help us devise critical tests that can implicate one error process to the exclusion of others. To this end, we present a model that captures distinct features of frameshifting cost and apply this to 641 prokaryotic genomes. We demonstrate that, although it is commonly assumed that the ribosome encounters an off-frame stop codon soon after the frameshift and costs of mis-elongation are therefore limited, genomes with high GC content typically incur much larger per-error costs. We go on to derive the prediction, unique to frameshifting errors, that differences in translational robustness between the 5' and 3' ends of genes should be less pronounced in genomes with higher GC content. This prediction we show to be correct. Surprisingly, this does not mean that GC-rich organisms necessarily carry a greater fitness burden as a consequence of accidental frameshifting. Indeed, increased per-error costs are often more than counterbalanced by lower predicted error rates owing to more diverse anticodon repertoires in GC-rich genomes. We therefore propose that selection on tRNA repertoires may operate to reduce frameshifting errors.

Genomic structure of an economically important cyanobacterium, Arthrospira (Spirulina) platensis NIES-39.

A filamentous non-N(2)-fixing cyanobacterium, Arthrospira (Spirulina) platensis, is an important organism for industrial applications and as a food supply. Almost the complete genome of A. platensis NIES-39 was determined in this study. The genome structure of A. platensis is estimated to be a single, circular chromosome of 6.8 Mb, based on optical mapping. Annotation of this 6.7 Mb sequence yielded 6630 protein-coding genes as well as two sets of rRNA genes and 40 tRNA genes. Of the protein-coding genes, 78% are similar to those of other organisms; the remaining 22% are currently unknown. A total 612 kb of the genome comprise group II introns, insertion sequences and some repetitive elements. Group I introns are located in a protein-coding region. Abundant restriction-modification systems were determined. Unique features in the gene composition were noted, particularly in a large number of genes for adenylate cyclase and haemolysin-like Ca(2+)-binding proteins and in chemotaxis proteins. Filament-specific genes were highlighted by comparative genomic analysis.