Replacement of the Arginine Biosynthesis Operon in Xanthomonadales by Lateral Gene Transfer

The role of lateral gene transfer (LGT) in prokaryotes has been shown to rapidly change the genome content, providing new gene tools for environmental adaptation. Features related to pathogenesis and resistance to strong selective conditions have been widely shown to be products of gene transfer between bacteria. The genomes of the γ-proteobacteria from the genus Xanthomonas, composed mainly of phytopathogens, have potential genomic islands that may represent imprints of such evolutionary processes. In this work, the evolution of genes involved in the pathway responsible for arginine biosynthesis in Xanthomonadales was investigated, and several lines of evidence point to the foreign origin of the arg genes clustered within a potential operon. Their presence inside a potential genomic island, bordered by a tRNA gene, the unusual ranking of sequence similarity, and the atypical phylogenies indicate that the metabolic pathway for arginine biosynthesis was acquired through LGT in the Xanthomonadales group. Moreover, although homologues were also found in Bacteroidetes (Flavobacteria group), for many of the genes analyzed close homologues are detected in different life domains (Eukarya and Archaea), indicating that the source of these arg genes may have been outside the Bacteria clade. The possibility of replacement of a complete primary metabolic pathway by LGT events supports the selfish operon hypothesis and may occur only under very special environmental conditions. Such rare events reveal part of the history of these interesting mosaic Xanthomonadales genomes, disclosing the importance of gene transfer modifying primary metabolism pathways and extending the scenario for bacterial genome evolution.     

The Mitochondrial Genome of Paramphistomum cervi (Digenea), the First Representative for the Family Paramphistomidae

We determined the complete mitochondrial DNA (mtDNA) sequence of a fluke, Paramphistomum cervi (Digenea: Paramphistomidae). This genome (14,014 bp) is slightly larger than that of Clonorchis sinensis (13,875 bp), but smaller than those of other digenean species. The mt genome of P. cervi contains 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and 2 non-coding regions (NCRs), a complement consistent with those of other digeneans. The arrangement of protein-coding and ribosomal RNA genes in the P. cervi mitochondrial genome is identical to that of other digeneans except for a group of Schistosoma species that exhibit a derived arrangement. The positions of some transfer RNA genes differ. Bayesian phylogenetic analyses, based on concatenated nucleotide sequences and amino-acid sequences of the 12 protein-coding genes, placed P. cervi within the Order Plagiorchiida, but relationships depicted within that order were not quite as expected from previous studies. The complete mtDNA sequence of P. cervi provides important genetic markers for diagnostics, ecological and evolutionary studies of digeneans.     

Construction of a Phylogenetic Tree of Photosynthetic Prokaryotes Based on Average Similarities of Whole Genome Sequences

Phylogenetic trees have been constructed for a wide range of organisms using gene sequence information, especially through the identification of orthologous genes that have been vertically inherited. The number of available complete genome sequences is rapidly increasing, and many tools for construction of genome trees based on whole genome sequences have been proposed. However, development of a reasonable method of using complete genome sequences for construction of phylogenetic trees has not been established. We have developed a method for construction of phylogenetic trees based on the average sequence similarities of whole genome sequences. We used this method to examine the phylogeny of 115 photosynthetic prokaryotes, i.e., cyanobacteria, Chlorobi, proteobacteria, Chloroflexi, Firmicutes and nonphotosynthetic organisms including Archaea. Although the bootstrap values for the branching order of phyla were low, probably due to lateral gene transfer and saturated mutation, the obtained tree was largely consistent with the previously reported phylogenetic trees, indicating that this method is a robust alternative to traditional phylogenetic methods.     

The Mitochondrial Genome of an Aquatic Plant,Spirodela polyrhiza

Background

Spirodela polyrhiza is a species of the order Alismatales, which represent the basal lineage of monocots with more ancestral features than the Poales. Its complete sequence of the mitochondrial (mt) genome could provide clues for the understanding of the evolution of mt genomes in plant.

Methods

Spirodela polyrhiza mt genome was sequenced from total genomic DNA without physical separation of chloroplast and nuclear DNA using the SOLiD platform. Using a genome copy number sensitive assembly algorithm, the mt genome was successfully assembled. Gap closure and accuracy was determined with PCR products sequenced with the dideoxy method.

Conclusions

This is the most compact monocot mitochondrial genome with 228,493 bp. A total of 57 genes encode 35 known proteins, 3 ribosomal RNAs, and 19 tRNAs that recognize 15 amino acids. There are about 600 RNA editing sites predicted and three lineage specific protein-coding-gene losses. The mitochondrial genes, pseudogenes, and other hypothetical genes (ORFs) cover 71,783 bp (31.0%) of the genome. Imported plastid DNA accounts for an additional 9,295 bp (4.1%) of the mitochondrial DNA. Absence of transposable element sequences suggests that very few nuclear sequences have migrated into Spirodela mtDNA. Phylogenetic analysis of conserved protein-coding genes suggests that Spirodela shares the common ancestor with other monocots, but there is no obvious synteny between Spirodela and rice mtDNAs. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly four-fifths of the Spirodela mitochondrial genome is of unknown origin and function. Although it contains a similar chloroplast DNA content and range of RNA editing as other monocots, it is void of nuclear insertions, active gene loss, and comprises large regions of sequences of unknown origin in non-coding regions. Moreover, the lack of synteny with known mitochondrial genomic sequences shed new light on the early evolution of monocot mitochondrial genomes.     

Characterization of Insertion Sites in Rainbow Papaya, the First Commercialized Transgenic Fruit Crop

Inserts and insert sites in transgenic, papaya ringspot virus (PRSV)-resistant commercial papaya Rainbow and SunUp, were characterized as part of a petition to Japan to allow import of fresh fruit of these cultivars from the U.S. and to provide data for a larger study aimed at understanding the global impact of DNA transformation on whole genome structure. The number and types of inserts were determined by Southern analysis using probes spanning the entire transformation plasmid and their sequences determined from corresponding clones or sequence reads from the whole-genome shotgun (WGS) sequence of SunUp papaya. All the functional transgenes, coding for the PRSV coat protein (CP), neophosphotransferase (nptII) and β-glucuronidase (uidA) were found in a single 9,789 basepair (bp) insert. Only two other inserts, one consisting of a 290 bp nonfunctional fragment of the nptII gene and a 1,533 bp plasmid-derived fragment containing a nonfunctional 222 bp segment of the tetA gene were detected in Rainbow and SunUp. Detection of the same three inserts in samples representing transgenic generations five to eight (R5 to R8) suggests that the three inserts are stably inherited. Five out of the six genomic DNA segments flanking the three inserts were nuclear plastid sequences (nupts). From the biosafety standpoint, no changes to endogenous gene function based on sequence structure of the transformation plasmid DNA insertion sites could be determined and no allergenic or toxic proteins were predicted from analysis of the insertion site and flanking genomic DNA.     

Usefulness of rpoB Gene Sequencing for Identification of Afipia and Bosea Species, Including a Strategy for Choosing Discriminative Partial Sequences

Bacteria belonging to the genera Afipia and Bosea are amoeba-resisting bacteria that have been recently reported to colonize hospital water supplies and are suspected of being responsible for intensive care unit-acquired pneumonia. Identification of these bacteria is now based on determination of the 16S ribosomal DNA sequence. However, the 16S rRNA gene is not polymorphic enough to ensure discrimination of species defined by DNA-DNA relatedness. The complete rpoB sequences of 20 strains were first determined by both PCR and genome walking methods. The percentage of homology between different species ranged from 83 to 97% and was in all cases lower than that observed with the 16S rRNA gene; this was true even for species that differed in only one position. The taxonomy of Bosea and Afipia is discussed in light of these results. For strain identification that does not require the complete rpoB sequence (4,113 to 4,137 bp), we propose a simple computerized method that allows determination of nucleotide positions of high variability in the sequence that are bordered by conserved sequences and that could be useful for design of universal primers. A fragment of 740 to 752 bp that contained the most highly variable area (positions 408 to 420) was amplified and sequenced with these universal primers for 47 strains. The variability of this sequence allowed identification of all strains and correlated well with results of DNA-DNA relatedness. In the future, this method could be also used for the determination of variability “hot spots” in sets of housekeeping genes, not only for identification purposes but also for increasing the discriminatory power of sequence typing techniques such as multilocus sequence typing.     

The complete mitochondrial genome of the intertidal copepod Tigriopus sp. (Copepoda, Harpactidae) from Korea and phylogenetic considerations

This paper reports on the basic characteristics of the mitochondrial genome of the intertidal copepod Tigriopus sp. from Korea, including its structural organization, base composition of rRNAs and protein-encoding genes, and the secondary structure of tRNAs. We amplified the complete mitochondrial DNA of the intertidal copepod Tigriopus sp. from Korea (sampling site Busan) by long-polymerase chain reaction (long-PCR) with conserved primers and sequenced this mitogenome by primer walking using flanking sequences as sequencing primers. The primer informations were obtained as expressed sequence tags (ESTs) from Tigriopus sp. The resultant Tigriopus sp. mitochondrial DNA sequence was 14,301 bp with a conserved structural organization, compared to that of T. japonicus from Japan with significant differences in several protein-coding regions including rRNAs, although the genomic organization of the mitochondrial genome was identical. In order to investigate biogeographic differences within the genus Tigriopus, we analyzed the CO1 gene by sequencing. This way, we compared several Tigriopus species from Korea, Japan, Hong Kong and Taiwan as well as other related species such as T. californicus, T. brevicornis and T. fulvus. The results further support the notion that the copepods display significantly different genomes within the same genus. These findings provide valuable genomic information for further studies on the population genetics and speciation processes within the genus Tigriopus.     

Genome-wide detection and analysis of homologous recombination among sequenced strains of Escherichia coli

Background  

Comparisons of complete bacterial genomes reveal evidence of lateral transfer of DNA across otherwise clonally diverging lineages. Some lateral transfer events result in acquisition of novel genomic segments and are easily detected through genome comparison. Other more subtle lateral transfers involve homologous recombination events that result in substitution of alleles within conserved genomic regions. This type of event is observed infrequently among distantly related organisms. It is reported to be more common within species, but the frequency has been difficult to quantify since the sequences under comparison tend to have relatively few polymorphic sites.     

Transgene constructs in coho salmon (Oncorhynchus kisutch) are repeated in a head-to-tail fashion and can be integrated adjacent to horizontally-transmitted parasite DNA

Currently, little information is available regarding the molecular organization of integrated transgenes in genetically-engineered fish. We performed a detailed structural analysis of an inserted transgene in one strain (M77) of transgenic coho salmon (Oncorhynchus kisutch) containing a salmon growth hormone gene construct (OnMTGH1). Microinjected DNA was found to have inserted into a single site in the coho salmon genome, and was organized with four complete internal copies and two partial terminal copies of the OnMTGH1 construct. All construct copies were organized in a direct-tandem (head-to-tail) repeat fashion in strain M77 and five additional strains (one also possessed a second recombinant junction fragment). For strain M77, the junctions between the transgene insert and the insertion point within the wild-type genome were cloned from strain-specific cosmid libraries and sequenced, revealing that the transgene insertion was accompanied by a deletion of 587 bp of wild-type DNA as well as a small insertion (19 bp) of unknown DNA upstream and a 14 bp direct- tandem duplication of sequence downstream. Upstream and downstream wild-type DNA sequence contained several repetitive sequence elements based on Southern blot analysis and homology to repetitive sequences in GenBank. In the downstream flank, a pseudogene sequence was also identified which has high homology to the CA membrane protein gene from Schistosoma japonicum, a parasite closely related to Sanguinicola sp. parasites which infect salmonids. Whether the presence of an inserted transgene and the presence of potentially horizontally-transmitted DNA are indicative of a genomic region with a predisposition for insertion of foreign DNA requires further study. The information derived from this transgene structure provides information useful for comparison to other transgenic organisms and for determination of the mechanism of transgene integration in lower vertebrates.     

Prediction of signal recognition particle RNA genes

We describe a method for prediction of genes that encode the RNA component of the signal recognition particle (SRP). A heuristic search for the strongly conserved helix 8 motif of SRP RNA is combined with covariance models that are based on previously known SRP RNA sequences. By screening available genomic sequences we have identified a large number of novel SRP RNA genes and we can account for at least one gene in every genome that has been completely sequenced. Novel bacterial RNAs include that of Thermotoga maritima, which, unlike all other non-gram-positive eubacteria, is predicted to have an Alu domain. We have also found the RNAs of Lactococcus lactis and Staphylococcus to have an unusual UGAC tetraloop in helix 8 instead of the normal GNRA sequence. An investigation of yeast RNAs reveals conserved sequence elements of the Alu domain that aid in the analysis of these RNAs. Analysis of the human genome reveals only two likely genes, both on chromosome 14. Our method for SRP RNA gene prediction is the first convenient tool for this task and should be useful in genome annotation.     

A genomic analysis of the archaeal system Ignicoccus hospitalis-Nanoarchaeum equitans

Background

The relationship between the hyperthermophiles Ignicoccus hospitalis and Nanoarchaeum equitans is the only known example of a specific association between two species of Archaea. Little is known about the mechanisms that enable this relationship.

Results

We sequenced the complete genome of I. hospitalis and found it to be the smallest among independent, free-living organisms. A comparative genomic reconstruction suggests that the I. hospitalis lineage has lost most of the genes associated with a heterotrophic metabolism that is characteristic of most of the Crenarchaeota. A streamlined genome is also suggested by a low frequency of paralogs and fragmentation of many operons. However, this process appears to be partially balanced by lateral gene transfer from archaeal and bacterial sources.

Conclusions

A combination of genomic and cellular features suggests highly efficient adaptation to the low energy yield of sulfur-hydrogen respiration and efficient inorganic carbon and nitrogen assimilation. Evidence of lateral gene exchange between N. equitans and I. hospitalis indicates that the relationship has impacted both genomes. This association is the simplest symbiotic system known to date and a unique model for studying mechanisms of interspecific relationships at the genomic and metabolic levels.     

Exploring the genome of the salt-marsh Spartina maritima (Poaceae, Chloridoideae) through BAC end sequence analysis

Spartina species play an important ecological role on salt marshes. Spartina maritima is an Old-World species distributed along the European and North-African Atlantic coasts. This hexaploid species (2n = 6x = 60, 2C = 3,700 Mb) hybridized with different Spartina species introduced from the American coasts, which resulted in the formation of new invasive hybrids and allopolyploids. Thus, S. maritima raises evolutionary and ecological interests. However, genomic information is dramatically lacking in this genus. In an effort to develop genomic resources, we analysed 40,641 high-quality bacterial artificial chromosome-end sequences (BESs), representing 26.7 Mb of the S. maritima genome. BESs were searched for sequence homology against known databases. A fraction of 16.91 % of the BESs represents known repeats including a majority of long terminal repeat (LTR) retrotransposons (13.67 %). Non-LTR retrotransposons represent 0.75 %, DNA transposons 0.99 %, whereas small RNA, simple repeats and low-complexity sequences account for 1.38 % of the analysed BESs. In addition, 4,285 simple sequence repeats were detected. Using the coding sequence database of Sorghum bicolor, 6,809 BESs found homology accounting for 17.1 % of all BESs. Comparative genomics with related genera reveals that the microsynteny is better conserved with S. bicolor compared to other sequenced Poaceae, where 37.6 % of the paired matching BESs are correctly orientated on the chromosomes. We did not observe large macrosyntenic rearrangements using the mapping strategy employed. However, some regions appeared to have experienced rearrangements when comparing Spartina to Sorghum and to Oryza. This work represents the first overview of S. maritima genome regarding the respective coding and repetitive components. The syntenic relationships with other grass genomes examined here help clarifying evolution in Poaceae, S. maritima being a part of the poorly-known Chloridoideae sub-family.     

Interkingdom gene fusions

Background  

Genome comparisons have revealed major lateral gene transfer between the three primary kingdoms of life - Bacteria, Archaea, and Eukarya. Another important evolutionary phenomenon involves the evolutionary mobility of protein domains that form versatile multidomain architectures. We were interested in investigating the possibility of a combination of these phenomena, with an invading gene merging with a pre-existing gene in the recipient genome.     

Horizontal gene transfer: sustaining pathogenicity and optimizing host–pathogen interactions

Successful host–pathogen interactions require the presence, maintenance and expression of gene cassettes called 'pathogenicity islands' (PAIs) and 'metabolic islands' (MAIs) in the respective pathogen. The products of these genes confer on the pathogen the means to recognize their host(s) and to efficiently evade host defences in order to colonize, propagate within the host and eventually disseminate from the host. Virulence effectors secreted by type III and type IV secretion systems, among others, play vital roles in sustaining pathogenicity and optimizing host–pathogen interactions. Complete genome sequences of plant pathogenic bacteria have revealed the presence of PAIs and MAIs. The genes of these islands possess mosaic structures with regions displaying differences in nucleotide composition and codon usage in relation to adjacent genome structures, features that are highly suggestive of their acquisition from a foreign donor. These donors can be other bacteria, as well as lower members of the Archaea and Eukarya. Genes that have moved from the domains Archaea and Eukarya to the domain Bacteria are true cases of horizontal gene transfer. They represent interdomain genetic transfer. Genetic exchange between distinct members of the domain Bacteria, however, represents lateral gene transfer, an intradomain event. Both horizontal and lateral gene transfer events have been used to facilitate survival fitness of the pathogen.     

A-tract clusters may facilitate DNA packaging in bacterial nucleoid

Molecular mechanisms of bacterial chromosome packaging are still unclear, as bacteria lack nucleosomes or other apparent basic elements of DNA compaction. Among the factors facilitating DNA condensation may be a propensity of the DNA molecule for folding due to its intrinsic curvature. As suggested previously, the sequence correlations in genome reflect such a propensity [Trifonov and Sussman (1980) Proc. Natl Acad. Sci. USA, 77, 3816–3820]. To further elaborate this concept, we analyzed positioning of A-tracts (the sequence motifs introducing the most pronounced DNA curvature) in the Escherichia coli genome. First, we observed that the A-tracts are over-represented and distributed ‘quasi-regularly’ throughout the genome, including both the coding and intergenic sequences. Second, there is a 10–12 bp periodicity in the A-tract positioning indicating that the A-tracts are phased with respect to the DNA helical repeat. Third, the phased A-tracts are organized in ~100 bp long clusters. The latter feature was revealed with the help of a novel approach based on the Fourier series expansion of the A-tract distance autocorrelation function. Since the A-tracts introduce local bends of the DNA duplex and these bends accumulate when properly phased, the observed clusters would facilitate DNA looping. Also, such clusters may serve as binding sites for the nucleoid-associated proteins that have affinities for curved DNA (such as HU, H-NS, Hfq and CbpA). Therefore, we suggest that the ~100 bp long clusters of the phased A-tracts constitute the ‘structural code’ for DNA compaction by providing the long-range intrinsic curvature and increasing stability of the DNA complexes with architectural proteins.     

Swine and poultry pathogens: the complete genome sequences of two strains of Mycoplasma hyopneumoniae and a strain of Mycoplasma synoviae

This work reports the results of analyses of three complete mycoplasma genomes, a pathogenic (7448) and a nonpathogenic (J) strain of the swine pathogen Mycoplasma hyopneumoniae and a strain of the avian pathogen Mycoplasma synoviae; the genome sizes of the three strains were 920,079 bp, 897,405 bp, and 799,476 bp, respectively. These genomes were compared with other sequenced mycoplasma genomes reported in the literature to examine several aspects of mycoplasma evolution. Strain-specific regions, including integrative and conjugal elements, and genome rearrangements and alterations in adhesin sequences were observed in the M. hyopneumoniae strains, and all of these were potentially related to pathogenicity. Genomic comparisons revealed that reduction in genome size implied loss of redundant metabolic pathways, with maintenance of alternative routes in different species. Horizontal gene transfer was consistently observed between M. synoviae and Mycoplasma gallisepticum. Our analyses indicated a likely transfer event of hemagglutinin-coding DNA sequences from M. gallisepticum to M. synoviae.