GraphDNA: a Java program for graphical display of DNA composition analyses

Contrasting modes of mitochondrial genome evolution in sister taxa of wood-eating marine bivalves (Teredinidae and Xylophagaidae)

The bivalve families Teredinidae and Xylophagaidae include voracious consumers of wood in shallow and deep-water marine environments, respectively. The taxa are sister clades whose members consume wood as food with the aid of intracellular cellulolytic endosymbionts housed in their gills. This combination of adaptations is found in no other group of animals and was likely present in the common ancestor of both families. Despite these commonalities, the two families have followed dramatically different evolutionary paths with respect to anatomy, life history and distribution. Here we present 42 new mitochondrial genome sequences from Teredinidae and Xylophagaidae and show that distinct trajectories have also occurred in the evolution and organization of their mitochondrial genomes. Teredinidae display significantly greater rates of amino acid substitution but absolute conservation of protein-coding gene order, whereas Xylophagaidae display significantly less amino acid change but have undergone numerous and diverse changes in genome organization since their divergence from a common ancestor. As with many bivalves, these mitochondrial genomes encode two ribosomal RNAs, 12 protein coding genes, and 22 tRNAs; atp8 was not detected. We further show that their phylogeny, as inferred from amino acid sequences of 12 concatenated mitochondrial protein-coding genes, is largely congruent with those inferred from their nuclear genomes based on 18S and 28S ribosomal RNA sequences. Our results provide a robust phylogenetic framework to explore the tempo and mode of mitochondrial genome evolution and offer directions for future phylogenetic and taxonomic studies of wood-boring bivalves.

SkewDB, a comprehensive database of GC and 10 other skews for over 30,000 chromosomes and plasmids

GC skew denotes the relative excess of G nucleotides over C nucleotides on the leading versus the lagging replication strand of eubacteria. While the effect is small, typically around 2.5%, it is robust and pervasive. GC skew and the analogous TA skew are a localized deviation from Chargaff’s second parity rule, which states that G and C, and T and A occur with (mostly) equal frequency even within a strand. Different bacterial phyla show different kinds of skew, and differing relations between TA and GC skew. This article introduces an open access database (https://skewdb.org) of GC and 10 other skews for over 30,000 chromosomes and plasmids. Further details like codon bias, strand bias, strand lengths and taxonomic data are also included. The SkewDB can be used to generate or verify hypotheses. Since the origins of both the second parity rule and GC skew itself are not yet satisfactorily explained, such a database may enhance our understanding of prokaryotic DNA.

Measurement(s)	Imbalances in the use of DNA nucleotides
Technology Type(s)	Next Generation Sequencing
Factor Type(s)	Position within DNA sequence • Organism type
Sample Characteristic - Organism	bacterium • archaea
Sample Characteristic - Environment	Varying
Sample Characteristic - Location	World

Computational prediction of species-specific yeast DNA replication origin via iterative feature representation

Deoxyribonucleic acid replication is one of the most crucial tasks taking place in the cell, and it has to be precisely regulated. This process is initiated in the replication origins (ORIs), and thus it is essential to identify such sites for a deeper understanding of the cellular processes and functions related to the regulation of gene expression. Considering the important tasks performed by ORIs, several experimental and computational approaches have been developed in the prediction of such sites. However, existing computational predictors for ORIs have certain curbs, such as building only single-feature encoding models, limited systematic feature engineering efforts and failure to validate model robustness. Hence, we developed a novel species-specific yeast predictor called yORIpred that accurately identify ORIs in the yeast genomes. To develop yORIpred, we first constructed optimal 40 baseline models by exploring eight different sequence-based encodings and five different machine learning classifiers. Subsequently, the predicted probability of 40 models was considered as the novel feature vector and carried out iterative feature learning approach independently using five different classifiers. Our systematic analysis revealed that the feature representation learned by the support vector machine algorithm (yORIpred) could well discriminate the distribution characteristics between ORIs and non-ORIs when compared with the other four algorithms. Comprehensive benchmarking experiments showed that yORIpred achieved superior and stable performance when compared with the existing predictors on the same training datasets. Furthermore, independent evaluation showcased the best and accurate performance of yORIpred thus underscoring the significance of iterative feature representation. To facilitate the users in obtaining their desired results without undergoing any mathematical, statistical or computational hassles, we developed a web server for the yORIpred predictor, which is available at: http://thegleelab.org/yORIpred.

iRO-PsekGCC: Identify DNA Replication Origins Based on Pseudo k-Tuple GC Composition

Identification of replication origins is playing a key role in understanding the mechanism of DNA replication. This task is of great significance in DNA sequence analysis. Because of its importance, some computational approaches have been introduced. Among these predictors, the iRO-3wPseKNC predictor is the first discriminative method that is able to correctly identify the entire replication origins. For further improving its predictive performance, we proposed the Pseudo k-tuple GC Composition (PsekGCC) approach to capture the "GC asymmetry bias" of yeast species by considering both the GC skew and the sequence order effects of k-tuple GC Composition (k-GCC) in this study. Based on PseKGCC, we proposed a new predictor called iRO-PsekGCC to identify the DNA replication origins. Rigorous jackknife test on two yeast species benchmark datasets (Saccharomyces cerevisiae, Pichia pastoris) indicated that iRO-PsekGCC outperformed iRO-3wPseKNC. It can be anticipated that iRO-PsekGCC will be a useful tool for DNA replication origin identification. Availability and implementation: The web-server for the iRO-PsekGCC predictor was established, and it can be accessed at http://bliulab.net/iRO-PsekGCC/.

Novel diversity in mitochondrial genomes of deep-sea Pennatulacea (Cnidaria: Anthozoa: Octocorallia)

We present the first documented complete mitogenomes of deep-sea Pennatulacea, representing nine genera and eight families. These include one species each of the deep-sea genera Funiculina, Halipteris, Protoptilum and Distichoptilum, four species each of Umbellula and Pennatula, three species of Kophobelemnon and two species of Anthoptilum, as well as one species of the epi- and mesobenthic genus Virgularia. Seventeen circular genomes ranged from 18,513 bp (Halipteris cf. finmarchica) to 19,171 bp (Distichoptilum gracile) and contained all genes standard to octocoral mitochondrial genomes (14 protein-coding genes, two ribosomal RNA genes and one transfer RNA). We found at least three different gene orders in Pennatulacea: the ancestral gene order, the gene order found in bamboo corals (Family Isididae), and a novel gene order. The mitogenome of one species of Umbellula has a bipartite genome (∼13 kbp and ∼5 kbp), with good evidence that both parts are circular.

DNAvisualization.org: a serverless web tool for DNA sequence visualization

Raw DNA sequences contain an immense amount of meaningful biological information. However, these sequences are hard for humans to intuitively interpret. To solve this problem, a number of methods have been proposed to transform DNA sequences into two-dimensional visualizations. DNAvisualization.org implements several of these methods in a cost effective and performant manner via a novel, entirely serverless architecture. By taking advantage of recent developments in serverless parallel computing and selective data retrieval, the website is able to offer users the ability to visualize up to thirty 4.5 Mb DNA sequences simultaneously using one of five supported methods and to export these visualizations in a variety of publication-ready formats.

Mitochondrial Architecture Rearrangements Produce Asymmetrical Nonadaptive Mutational Pressures That Subvert the Phylogenetic Reconstruction in Isopoda

The phylogeny of Isopoda, a speciose order of crustaceans, remains unresolved, with different data sets (morphological, nuclear, mitochondrial) often producing starkly incongruent phylogenetic hypotheses. We hypothesized that extreme diversity in their life histories might be causing compositional heterogeneity/heterotachy in their mitochondrial genomes, and compromising the phylogenetic reconstruction. We tested the effects of different data sets (mitochondrial, nuclear, nucleotides, amino acids, concatenated genes, individual genes, gene orders), phylogenetic algorithms (assuming data homogeneity, heterogeneity, and heterotachy), and partitioning; and found that almost all of them produced unique topologies. As we also found that mitogenomes of Asellota and two Cymothoida families (Cymothoidae and Corallanidae) possess inversed base (GC) skew patterns in comparison to other isopods, we concluded that inverted skews cause long-branch attraction phylogenetic artifacts between these taxa. These asymmetrical skews are most likely driven by multiple independent inversions of origin of replication (i.e., nonadaptive mutational pressures). Although the PhyloBayes CAT-GTR algorithm managed to attenuate some of these artifacts (and outperform partitioning), mitochondrial data have limited applicability for reconstructing the phylogeny of Isopoda. Regardless of this, our analyses allowed us to propose solutions to some unresolved phylogenetic debates, and support Asellota are the most likely candidate for the basal isopod branch. As our findings show that architectural rearrangements might produce major compositional biases even on relatively short evolutionary timescales, the implications are that proving the suitability of data via composition skew analyses should be a prerequisite for every study that aims to use mitochondrial data for phylogenetic reconstruction, even among closely related taxa.

Linear Mitochondrial Genome in Anthozoa (Cnidaria): A Case Study in Ceriantharia

Sequences and structural attributes of mitochondrial genomes have played a critical role in the clarification of relationships among Cnidaria, a key phylum of early-diverging animals. Among the major lineages of Cnidaria, Ceriantharia (“tube anemones”) remains one of the most enigmatic in terms of its phylogenetic position. We sequenced the mitochondrial genomes of two ceriantharians to see whether the complete organellar genome would provide more support for the phylogenetic placement of Ceriantharia. For both Isarachnanthus nocturnus and Pachycerianthus magnus, the mitochondrial gene sequences could not be assembled into a single circular genome. Instead, our analyses suggest that both species have mitochondrial genomes consisting of multiple linear fragments. Linear mitogenomes are characteristic of members of Medusozoa, one of the major lineages of Cnidaria, but are unreported for Anthozoa, which includes the Ceriantharia. The inferred number of fragments and variation in gene order between species is much greater within Ceriantharia than among the lineages of Medusozoa. We identify origins of replication for each of the five putative chromosomes of the Isarachnanthus nocturnus mitogenome and for each of the eight putative chromosomes of the Pachycerianthus magnus mitogenome. At 80,923 bp, I. nocturnus now holds the record for the largest animal mitochondrial genome reported to date. The novelty of the mitogenomic structure in Ceriantharia highlights the distinctiveness of this lineage but, because it appears to be both unique to and diverse within Ceriantharia, it is uninformative about the phylogenetic position of Ceriantharia relative to other Anthozoa. The presence of tRNA^Met and tRNA^Trp in both ceriantharian mitogenomes supports a closer relationship between Ceriantharia and Hexacorallia than between Ceriantharia and any other cnidarian lineage, but phylogenetic analysis of the genes contained in the mitogenomes suggests that Ceriantharia is sister to a clade containing Octocorallia + Hexacorallia indicating a possible suppression of tRNA^Trp in Octocorallia.

Comparative genomic analysis for nucleotide, codon, and amino acid usage patterns of mycoplasmas

The evolutionary factors in influencing the genetic characteristics of nucleotide, synonymous codon, and amino acid usage of 18 mycoplasma species were analyzed. The nucleotide usage at the 1st and 2nd codon position which determines amino acid composition of proteins has a significant correlation with the total nucleotide composition of gene population of these mycoplasma species, however, the nucleotide usage at the 3rd codon position which affects synonymous codon usage patterns has a slight correlation with either the total nucleotide composition or the nucleotide usage at the 1st and 2nd codon position. Other evolutionary factors join in the evolutionary process of mycoplasma apart from mutation pressure caused by nucleotide usage constraint based on the relationships between effective number of codons/codon adaptation index and nucleotide usage at the 3rd codon position. Although nucleotide usage of gene population in mycoplasma dominates in forming the overall codon usage trends, the relative abundance of codon with nucleotide context and amino acid usage pattern show that translation selection involved in translation accuracy and efficiency play an important role in synonymous codon usage patterns. In addition, synonymous codon usage patterns of gene population have a bigger power to represent genetic diversity among different species than amino acid usage. These results suggest that although the mycoplasmas reduce its genome size during the evolutionary process and shape the form, which is opposite to their hosts, of AT usages at high levels, this kind organism still depends on nucleotide usage at the 1st and 2nd codon positions to control syntheses of the requested proteins for surviving in their hosts and nucleotide usage at the 3rd codon position to develop genetic diversity of different mycoplasma species. This systemic analysis with 18 mycoplasma species may provide useful clues for further in vivo genetic studies on the related species.

In vitro site-specific recombination mediated by the tyrosine recombinase XerA of Thermoplasma acidophilum

In organisms with circular chromosomes, such as bacteria and archaea, an odd number of homologous recombination events can generate a chromosome dimer. Such chromosome dimers cannot be segregated unless they are converted to monomers before cell division. In Escherichia coli, dimer-to-monomer conversion is mediated by the paralogous XerC and XerD recombinases at a specific dif site in the replication termination region. Dimer resolution requires the highly conserved cell division protein/chromosome translocase FtsK, and this site-specific chromosome resolution system is present or predicted in most bacteria. However, most archaea have only XerA, a homologue of the bacterial XerC/D proteins, but no homologues of FtsK. In addition, the molecular mechanism of XerA-mediated chromosome resolution in archaea has been less thoroughly elucidated than those of the corresponding bacterial systems. In this study, we identified two XerA-binding sites (dif1 and dif2) in the Thermoplasma acidophilum chromosome. In vitro site-specific recombination assays showed that dif2, but not dif1, serves as a target site for XerA-mediated chromosome resolution. Mutational analysis indicated that not only the core consensus sequence of dif2, but also its flanking regions play important roles in the recognition and recombination reactions mediated by XerA.

Phylogenetic relationships of Hemiptera inferred from mitochondrial and nuclear genes

Here, we reconstructed the Hemiptera phylogeny based on the expanded mitochondrial protein-coding genes and the nuclear 18S rRNA gene, separately. The differential rates of change across lineages may associate with long-branch attraction (LBA) effect and result in conflicting estimates of phylogeny from different types of data. To reduce the potential effects of systematic biases on inferences of topology, various data coding schemes, site removal method, and different algorithms were utilized in phylogenetic reconstruction. We show that the outgroups Phthiraptera, Thysanoptera, and the ingroup Sternorrhyncha share similar base composition, and exhibit "long branches" relative to other hemipterans. Thus, the long-branch attraction between these groups is suspected to cause the failure of recovering Hemiptera under the homogeneous model. In contrast, a monophyletic Hemiptera is supported when heterogeneous model is utilized in the analysis. Although higher level phylogenetic relationships within Hemiptera remain to be answered, consensus between analyses is beginning to converge on a stable phylogeny.

Expression of a non-coding RNA in ectromelia virus is required for normal plaque formation

Poxviruses are dsDNA viruses with large genomes. Many genes in the genome remain uncharacterized, and recent studies have demonstrated that the poxvirus transcriptome includes numerous so-called anomalous transcripts not associated with open reading frames. Here, we characterize the expression and role of an apparently non-coding RNA in orthopoxviruses, which we call viral hairpin RNA (vhRNA). Using a bioinformatics approach, we predicted expression of a transcript not associated with an open reading frame that is likely to form a stem-loop structure due to the presence of a 21 nt palindromic sequence. Expression of the transcript as early as 2 h post-infection was confirmed by northern blot and analysis of publicly available vaccinia virus infected cell transcriptomes. The transcription start site was determined by RACE PCE and transcriptome analysis, and early and late promoter sequences were identified. Finally, to test the function of the transcript we generated an ectromelia virus knockout, which failed to form plaques in cell culture. The important role of the transcript in viral replication was further demonstrated using siRNA. Although the function of the transcript remains unknown, our work contributes to evidence of an increasingly complex poxvirus transcriptome, suggesting that transcripts such as vhRNA not associated with an annotated open reading frame can play an important role in viral replication.

Visual ModuleOrganizer: a graphical interface for the detection and comparative analysis of repeat DNA modules

BACKGROUND

DNA repeats, such as transposable elements, minisatellites and palindromic sequences, are abundant in sequences and have been shown to have significant and functional roles in the evolution of the host genomes. In a previous study, we introduced the concept of a repeat DNA module, a flexible motif present in at least two occurences in the sequences. This concept was embedded into ModuleOrganizer, a tool allowing the detection of repeat modules in a set of sequences. However, its implementation remains difficult for larger sequences.

RESULTS

Here we present Visual ModuleOrganizer, a Java graphical interface that enables a new and optimized version of the ModuleOrganizer tool. To implement this version, it was recoded in C++ with compressed suffix tree data structures. This leads to less memory usage (at least 120-fold decrease in average) and decreases by at least four the computation time during the module detection process in large sequences. Visual ModuleOrganizer interface allows users to easily choose ModuleOrganizer parameters and to graphically display the results. Moreover, Visual ModuleOrganizer dynamically handles graphical results through four main parameters: gene annotations, overlapping modules with known annotations, location of the module in a minimal number of sequences, and the minimal length of the modules. As a case study, the analysis of FoldBack4 sequences clearly demonstrated that our tools can be extended to comparative and evolutionary analyses of any repeat sequence elements in a set of genomic sequences. With the increasing number of sequences available in public databases, it is now possible to perform comparative analyses of repeated DNA modules in a graphic and friendly manner within a reasonable time period.

AVAILABILITY

Visual ModuleOrganizer interface and the new version of the ModuleOrganizer tool are freely available at: http://lcb.cnrs-mrs.fr/spip.php?rubrique313.

Unusual large-scale chromosomal rearrangements in Mycobacterium tuberculosis Beijing B0/W148 cluster isolates

The Mycobacterium tuberculosis (MTB) Beijing family isolates are geographically widespread, and there are examples of Beijing isolates that are hypervirulent and associated with drug resistance. One-fourth of Beijing genotype isolates found in Russia belong to the B0/W148 group. The aim of the present study was to investigate features of these endemic strains on a genomic level. Four Russian clinical isolates of this group were sequenced, and the data obtained was compared with published sequences of various MTB strain genomes, including genome of strain W-148 of the same B0/W148 group. The comparison of the W-148 and H37Rv genomes revealed two independent inversions of large segments of the chromosome. The same inversions were found in one of the studied strains after deep sequencing using both the fragment and mate-paired libraries. Additionally, inversions were confirmed by RFLP hybridization analysis. The discovered rearrangements were verified by PCR in all four newly sequenced strains in the study and in four additional strains of the same Beijing B0/W148 group. The other 32 MTB strains from different phylogenetic lineages were tested and revealed no inversions. We suggest that the initial largest inversion changed the orientation of the three megabase (Mb) segment of the chromosome, and the second one occurred in the previously inverted region and partly restored the orientation of the 2.1 Mb inner segment of the region. This is another remarkable example of genomic rearrangements in the MTB in addition to the recently published of large-scale duplications. The described cases suggest that large-scale genomic rearrangements in the currently circulating MTB isolates may occur more frequently than previously considered, and we hope that further studies will help to determine the exact mechanism of such events.

Complete mitochondrial genomes of the Japanese pink coral (Corallium elatius) and the Mediterranean red coral (Corallium rubrum): a reevaluation of the phylogeny of the family Coralliidae based on molecular data

Precious corals are soft corals belonging to the family Coralliidae (Anthozoa: Octocorallia: Alcyonacea) and class Anthozoa, whose skeletal axes are used for jewelry. The family Coralliidae includes ca. 40 species and was originally thought to comprise of the single genus Corallium. In 2003, Corallium was split into two genera, Corallium and Paracorallium, and seven species were moved to this newly identified genus on the bases of morphological features. Previously, we determined the complete mitochondrial genome sequence of two precious corals Paracorallium japonicum and Corallium konojoi, in order to clarify their systematic positions. The two genomes showed high nucleotide sequence identity, but their gene order arrangements were not identical. Here, we determined three complete mitochondrial genome sequences from the one specimen of Mediterranean Corallium rubrum and two specimens of Corallium elatius coming from Kagoshima (South Japan). The circular mitochondrial genomes of C. rubrum and C. elatius are 18,915bp and 18,969-18,970bp in length, respectively, and encode 14 typical octocorallian protein-coding genes (nad1-6, nad4L, cox1-3, cob, atp6, atp8, and mtMutS, which is an octocoral-specific mismatch repair gene homologue), two ribosomal RNA genes (rns and rnl), and one transfer RNA (trnM). The overall nucleotide differences between C. konojoi and each C. elatius haplotype (T2007 and I2011) are only 10 and 11 nucleotides, respectively; this degree of similarity indicates that C. elatius and C. konojoi are very closely related species. Notably, the C. rubrum mitochondrial genome shows more nucleotide sequence identity to P. japonicum (99.5%) than to its congeneric species C. konojoi (95.3%) and C. elatius (95.3%). Moreover, the gene order arrangement of C. rubrum was the same as that of P. japonicum, while that of C. elatius was the same as C. konojoi. Phylogenetic analysis based on three mitochondrial genes from 24 scleraxonian species shows that the family Coralliidae is separated into two distinct groups, recovering Corallium as a paraphyletic genus. Our results indicate that the currently accepted generic classification of Coralliidae should be reconsidered.

The mitochondrial genome of Paraminabea aldersladei (Cnidaria: Anthozoa: Octocorallia) supports intramolecular recombination as the primary mechanism of gene rearrangement in octocoral mitochondrial genomes

Sequencing of the complete mitochondrial genome of the soft coral Paraminabea aldersladei (Alcyoniidae) revealed a unique gene order, the fifth mt gene arrangement now known within the cnidarian subclass Octocorallia. At 19,886 bp, the mt genome of P. aldersladei is the second largest known for octocorals; its gene content and nucleotide composition are, however, identical to most other octocorals, and the additional length is due to the presence of two large, noncoding intergenic regions. Relative to the presumed ancestral octocoral gene order, in P. aldersladei a block of three protein-coding genes (nad6–nad3–nad4l) has been translocated and inverted. Mapping the distribution of mt gene arrangements onto a taxonomically comprehensive phylogeny of Octocorallia suggests that all of the known octocoral gene orders have evolved by successive inversions of one or more evolutionarily conserved blocks of protein-coding genes. This mode of genome evolution is unique among Metazoa, and contrasts strongly with that observed in Hexacorallia, in which extreme gene shuffling has occurred among taxonomic orders. Two of the four conserved gene blocks found in Octocorallia are, however, also conserved in the linear mt genomes of Medusozoa and in one group of Demospongiae. We speculate that the rate and mechanism of gene rearrangement in octocorals may be influenced by the presence in their mt genomes of mtMutS, a putatively active DNA mismatch repair protein that may also play a role in mediating intramolecular recombination.

A unique tRNA gene family and a novel, highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS(-AGY) and trnQ genes, a trnM isomerism, but it lacks trnS(-CUN). We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP→trnS(-AGY)) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, "Oxidored_q1_N", which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.

Complete genome and proteome of Acholeplasma laidlawii

We present the complete genome sequence and proteogenomic map for Acholeplasma laidlawii PG-8A (class Mollicutes, order Acholeplasmatales, family Acholeplasmataceae). The genome of A. laidlawii is represented by a single 1,496,992-bp circular chromosome with an average G+C content of 31 mol%. This is the longest genome among the Mollicutes with a known nucleotide sequence. It contains genes of polymerase type I, SOS response, and signal transduction systems, as well as RNA regulatory elements, riboswitches, and T boxes. This demonstrates a significant capability for the regulation of gene expression and mutagenic response to stress. Acholeplasma laidlawii and phytoplasmas are the only Mollicutes known to use the universal genetic code, in which UGA is a stop codon. Within the Mollicutes group, only the sterol-nonrequiring Acholeplasma has the capacity to synthesize saturated fatty acids de novo. Proteomic data were used in the primary annotation of the genome, validating expression of many predicted proteins. We also detected posttranslational modifications of A. laidlawii proteins: phosphorylation and acylation. Seventy-four candidate phosphorylated proteins were found: 16 candidates are proteins unique to A. laidlawii, and 11 of them are surface-anchored or integral membrane proteins, which implies the presence of active signaling pathways. Among 20 acylated proteins, 14 contained palmitic chains, and six contained stearic chains. No residue of linoleic or oleic acid was observed. Acylated proteins were components of mainly sugar and inorganic ion transport systems and were surface-anchored proteins with unknown functions.

The complete mitochondrial genome of Calicogorgia granulosa (Anthozoa: Octocorallia): potential gene novelty in unidentified ORFs formed by repeat expansion and segmental duplication

Mitochondrial genomes of many nonbilaterian animals show high diversity of genome size and gene content, revealing many intergenic regions (IGRs), diverse repeats and additional genes. Here we present a new complete mitogenome of the cnidarian sea fan species, Calicogorgia granulosa (Anthozoa: Octocorallia) and its novel genomic features. The 20,246 bp of the complete mitogenome, which is the largest among the nine octocorals sequenced to date, contains 13 protein coding genes, 2 rRNAs and a tRNA within its circular form of mitochondrial DNA. We found an identical segmental duplication (S1 and S2, 913 bp) composed of an ORF (672 bp) coding for a hypothetical protein within which Direct Variant Repeat (DVR) expansions reside in-frame to the coding sequence. Additionally, the duplicated segmental DNA showed no variation in nucleotide sequences both between S1 and S2 and across multiple individual samples. Upon these observations, we discuss plausible causes for the intramitochondrial segmental duplication and the absence of sequence variation, and a need for further investigation of the novel ORF as well. In conclusion the present mitogenome of C. granulosa adds more information to our understanding of the diversity and evolution of mitogenomes of nonbilaterian animals.

Complete mitochondrial genomes of two Japanese precious corals, Paracorallium japonicum and Corallium konojoi (Cnidaria, Octocorallia, Coralliidae): notable differences in gene arrangement

Precious coral are taxonomically a group of corals that belong to the family Coralliidae within the order Alcyonacea, subclass Octocorallia, and class Anthozoa, whose skeletal axes are used for jewelry. They are distributed in the Mediterranean Sea and in waters adjacent to Japan, Taiwan, Midway Island and the Hawaiian Islands. The genus Corallium of the family Coralliidae was recently divided into two genera, Corallium and Paracorallium, based on morphological observations, but insufficient molecular evidence to support this classification has been presented to date. We determined for the first time the complete mitochondrial genome sequence of two precious corals P. japonicum and C. konojoi, in order to clarify their systematic positions. The circular mitochondrial genomes of P. japonicum and C. konojoi are 18,913bp and 18,969bp in length, respectively, and encode 13 typical energy pathway protein coding genes (nad1-6, nad4L, cox1-3, cob, atp6 and atp8), two ribosomal RNA genes (rns and rnl), a transfer RNA (trnM) and a mismatch repair gene homologue msh1. The two genomes have an overall nucleotide sequence identity of 97.5%, which is comparable to that between Acanella eburnea and Keratoisidinae sp. belonging to Octocorallia. Surprisingly, however, their gene arrangements were not identical. Phylogenetic analyses using seven complete mitochondrial genome sequences belonging to species in the subclass Octocorallia indicated that within the subclass, at least three gene order rearrangement events occurred during evolution. Our results support the validity of the morphological classification that separated the family Coralliidae into two genera, Corallium and Paracorallium.

ModuleOrganizer: detecting modules in families of transposable elements

Background

Most known eukaryotic genomes contain mobile copied elements called transposable elements. In some species, these elements account for the majority of the genome sequence. They have been subject to many mutations and other genomic events (copies, deletions, captures) during transposition. The identification of these transformations remains a difficult issue. The study of families of transposable elements is generally founded on a multiple alignment of their sequences, a critical step that is adapted to transposons containing mostly localized nucleotide mutations. Many transposons that have lost their protein-coding capacity have undergone more complex rearrangements, needing the development of more complex methods in order to characterize the architecture of sequence variations.

Results

In this study, we introduce the concept of a transposable element module, a flexible motif present in at least two sequences of a family of transposable elements and built on a succession of maximal repeats. The paper proposes an assembly method working on a set of exact maximal repeats of a set of sequences to create such modules. It results in a graphical view of sequences segmented into modules, a representation that allows a flexible analysis of the transformations that have occurred between them. We have chosen as a demonstration data set in depth analysis of the transposable element Foldback in Drosophila melanogaster. Comparison with multiple alignment methods shows that our method is more sensitive for highly variable sequences. The study of this family and the two other families AtREP21 and SIDER2 reveals new copies of very different sizes and various combinations of modules which show the potential of our method.

Conclusions

ModuleOrganizer is available on the Genouest bioinformatics center at http://moduleorganizer.genouest.org

Mechanisms and evolution of control logic in prokaryotic transcriptional regulation

A major part of organismal complexity and versatility of prokaryotes resides in their ability to fine-tune gene expression to adequately respond to internal and external stimuli. Evolution has been very innovative in creating intricate mechanisms by which different regulatory signals operate and interact at promoters to drive gene expression. The regulation of target gene expression by transcription factors (TFs) is governed by control logic brought about by the interaction of regulators with TF binding sites (TFBSs) in cis-regulatory regions. A factor that in large part determines the strength of the response of a target to a given TF is motif stringency, the extent to which the TFBS fits the optimal TFBS sequence for a given TF. Advances in high-throughput technologies and computational genomics allow reconstruction of transcriptional regulatory networks in silico. To optimize the prediction of transcriptional regulatory networks, i.e., to separate direct regulation from indirect regulation, a thorough understanding of the control logic underlying the regulation of gene expression is required. This review summarizes the state of the art of the elements that determine the functionality of TFBSs by focusing on the molecular biological mechanisms and evolutionary origins of cis-regulatory regions.

In silico prediction of the origin of replication among bacteria: a case study of Bacteroides thetaiotaomicron

The initiation of chromosomal replication occurs only once during the prokaryote cell cycle. Some origins of replication have been experimentally determined and have led to the development of in silico approaches to find the origin of replication among other prokaryotes. DNA base composition asymmetry is the basis of numerous in silico methods used to detect the origin and terminus of replication in prokaryotes. However, the composition asymmetry does not allow us to locate precisely the positions of the origin and terminus. Since DNA replication is a key step in the cell cycle it is important to determine properly the origin and terminus regions. Therefore, we have reviewed here the methods, tools, and databases for predicting the origins and terminuses of replication, and we have proposed some complementary analyses to reinforce these predictions. These analyses include finding the dnaA gene and its binding sites; making BLAST analyses of the intergenic sequences compared to related species; studying the gene order around the origin sequence; and studying the distribution of the genes encoded in the leading versus the lagging strand.

Genome sequence of Thermofilum pendens reveals an exceptional loss of biosynthetic pathways without genome reduction

We report the complete genome of Thermofilum pendens, a deeply branching, hyperthermophilic member of the order Thermoproteales in the archaeal kingdom Crenarchaeota. T. pendens is a sulfur-dependent, anaerobic heterotroph isolated from a solfatara in Iceland. It is an extracellular commensal, requiring an extract of Thermoproteus tenax for growth, and the genome sequence reveals that biosynthetic pathways for purines, most amino acids, and most cofactors are absent. In fact, T. pendens has fewer biosynthetic enzymes than obligate intracellular parasites, although it does not display other features that are common among obligate parasites and thus does not appear to be in the process of becoming a parasite. It appears that T. pendens has adapted to life in an environment rich in nutrients. T. pendens was known previously to utilize peptides as an energy source, but the genome revealed a substantial ability to grow on carbohydrates. T. pendens is the first crenarchaeote and only the second archaeon found to have a transporter of the phosphotransferase system. In addition to fermentation, T. pendens may obtain energy from sulfur reduction with hydrogen and formate as electron donors. It may also be capable of sulfur-independent growth on formate with formate hydrogen lyase. Additional novel features are the presence of a monomethylamine:corrinoid methyltransferase, the first time that this enzyme has been found outside the Methanosarcinales, and the presence of a presenilin-related protein. The predicted highly expressed proteins do not include proteins encoded by housekeeping genes and instead include ABC transporters for carbohydrates and peptides and clustered regularly interspaced short palindromic repeat-associated proteins.

From GC skews to wavelets: a gentle guide to the analysis of compositional asymmetries in genomic data

Compositional asymmetries are pervasive in DNA sequences. They are the result of the asymmetric interactions between DNA and cellular mechanisms such as replication and transcription. Here, we review many of the methods that have been proposed over the years to analyse compositional asymmetries in DNA sequences. Among these we list GC skews, oligonucleotide skews and wavelets, which among other uses have been extensively employed to delimitate origins and termini of replication in genomes. We also review the use of multivariate methods, such as factorial correspondence analysis, discriminant analysis and analysis of variance, which allow assigning compositional strand asymmetries to the different biological processes shaping sequence composition. Finally, we review methods that have been used to infer substitution matrices and allow understanding the mutational processes underlying strand asymmetry. We focus on replication asymmetries because they have been more thoroughly studied, but the methods may be adapted, and often are, to other problems. Although strand asymmetry has been studied more frequently through compositional skews of nucleotides or oligonucleotides, we recall that, depending on the goal of the analysis, other methods may be more appropriate to answer certain biological questions. We also refer to programs freely available to analyse strand asymmetry.