First Sequenced Mitochondrial Genome from the Phylum Acanthocephala(Leptorhynchoides thecatus) and Its Phylogenetic Position Within Metazoa

The complete sequence of the mitochondrial genome of Leptorhynchoides thecatus (Acanthocephala) was determined, and a phylogenetic analysis was carried out to determine its placement within Metazoa. The genome is circular, 13,888 bp, and contains at least 36 of the 37 genes typically found in animal mitochondrial genomes. The genes for the large and small ribosomal RNA subunits are shorter than those of most metazoans, and the structures of most of the tRNA genes are atypical. There are two significant noncoding regions (377 and 294 bp), which are the best candidates for a control region; however, these regions do not appear similar to any of the control regions of other animals studied to date. The amino acid and nucleotide sequences of the protein coding genes of L. thecatus and 25 other metazoan taxa were used in both maximum likelihood and maximum parsimony phylogenetic analyses. Results indicate that among taxa with available mitochondrial genome sequences, Platyhelminthes is the closest relative to L. thecatus, which together are the sister taxon of Nematoda; however, long branches and/or base composition bias could be responsible for this result. The monophyly of Ecdysozoa, molting organisms, was not supported by any of the analyses. This study represents the first mitochondrial genome of an acanthocephalan to be sequenced and will allow further studies of systematics, population genetics, and genome evolution.Reviewing Editor: Dr. Rafael Zardoya The entire genome sequence has been deposited with the GenBank Data Libraries under-accession number AY562383.     

Recent amplification of miniature inverted-repeat transposable elements in the vector mosquito Culex pipiens: characterization of the Mimo family

We describe a new family of repetitive elements, named Mimo, from the mosquito Culex pipiens. Structural characteristics of these elements fit well with those of miniature inverted-repeat transposable elements (MITEs), which are ubiquitous and highly abundant in plant genomes. The occurrence of Mimo in C. pipiens provides new evidence that MITEs are not restricted to plant genomes, but may be widespread in arthropods as well. The copy number of Mimo elements in C. pipiens (1000 copies in a 540 Mb genome) supports the hypothesis that there is a positive correlation between genome size and the magnitude of MITE proliferation. In contrast to most MITE families described so far, members of the Mimo family share a high sequence conservation, which may reflect a recent amplification history in this species. In addition, we found that Mimo elements are a frequent nest for other MITE-like elements, suggesting that multiple and successive MITE transposition events have occurred very recently in the C. pipiens genome. Despite evidence for recent mobility of these MITEs, no element has been found to encode a protein; therefore, we do not know how they have transposed and have spread in the genome. However, some sequence similarities in terminal inverted-repeats suggest a possible filiation of some of these mosquito MITEs with pogo-like DNA transposons.     

Complete sequence of three plasmids from Bacillus thuringiensis INTA-FR7-4 environmental isolate and comparison with related plasmids from the Bacillus cereus group

Bacillus thuringiensis is an insect pathogen used worldwide as a bioinsecticide. It belongs to the Bacillus cereus sensu lato group as well as Bacillus anthracis and B. cereus. Plasmids from this group of organisms have been implicated in pathogenicity as they carry the genes responsible for different types of diseases that affect mammals and insects. Some plasmids, like pAW63 and pBT9727, encode a functional conjugation machinery allowing them to be transferred to a recipient cell. They also share extensive homology with the non-functional conjugation apparatus of pXO2 from B. anthracis. In this study we report the complete sequence of three plasmids from an environmental B. thuringiensis isolate from Argentina, obtained by a shotgun sequencing method. We obtained the complete nucleotide sequence of plasmids pFR12 (12 095 bp), pFR12.5 (12 459 bp) and pFR55 (55 712 bp) from B. thuringiensis INTA-FR7-4. pFR12 and pFR12.5 were classified as cryptic as they do not code for any obvious functions besides replication and mobilization. Both small plasmids were classified as RCR plasmids due to similarities with the replicases they encode. Plasmid pFR55 showed a structural organization similar to that observed for plasmids pAW63, pBT9727 and pXO2. pFR55 also shares a tra region with these plasmids, containing genes related to T4SS and conjugation. A comparison between pFR55 and conjugative plasmids led to the postulation that pFR55 is a conjugative plasmid. Genes related to replication functions in pFR55 are different to those described for plasmids with known complete sequences. pFR55 is the first completely sequenced plasmid with a replication machinery related to that of ori44. The analysis of the complete sequence of plasmids from an environmental isolate of B. thuringiensis permitted the identification of a near complete conjugation apparatus in pFR55, resembling those of plasmids pAW63, pBT9727 and pXO2. The availability of this sequence is a step forward in the study of the molecular basis of the conjugative process in Gram positive bacteria, particularly due to the similarity with known conjugation systems. It is also a contribution to the expansion of the non-pathogenic B. cereus plasmid gene pool.     

The Mitochondrial Genome of Xiphinema americanum sensu stricto (Nematoda: Enoplea): Considerable Economization in the Length and Structural Features of Encoded Genes

The complete sequence of the mitochondrial genome of the plant parasitic nematode Xiphinema americanum sensu stricto has been determined. At 12626bp it is the smallest metazoan mitochondrial genome reported to date. Genes are transcribed from both strands. Genes coding for 12 proteins, 2 rRNAs and 17 putative tRNAs (with the tRNA-C, I, N, S1, S2 missing) are predicted from the sequence. The arrangement of genes within the X. americanum mitochondrial genome is unique and includes gene overlaps. Comparisons with the mtDNA of other nematodes show that the small size of the X. americanum mtDNA is due to a combination of factors. The two mitochondrial rRNA genes are considerably smaller than those of other nematodes, with most of the protein encoding and tRNA genes also slightly smaller. In addition, five tRNAs genes are absent, lengthy noncoding regions are not present in the mtDNA, and several gene overlaps are present. [Reviewing Editor: Dr. Yues van de Peer] F. Lamberti: Deceased, 2004     

Search for Drosophila genes encoding a conserved domain present in the achaete-scute complex and myc proteins

Summary Several genes of the achaete-scute complex (ASC) of Drosophila melanogaster encode a 60 amino acids long conserved domain which shares a significant homology with a region of the vertebrate myc proteins. Based on these results, the existence of a family of Drosophila genes that would share both this conserved domain and the neurogenic function of the AS-C has been postulated. To test this proposal, we have searched a D. melanogaster genomic library with a probe that encodes the conserved domain. Only under very low stringency hybridization conditions, clones not belonging to the AS-C cross-hybridized with the probe. Those that gave the strongest signals were characterized. Sequencing of the cross-hybridizing regions showed that they had no significant homology with the conserved domain, the sequence similarity extending at the most for 37 nucleotides. Although our results do not conclusively disprove the existence of a family of AS-C-like genes, they indicate that the conservation of the domain would be lower than that found for shared motifs in other families of Drosophila developmental genes.     

Thermostable RNase P RNAs lacking P18 identified in the Aquificales

The RNase P RNA (rnpB) and protein (rnpA) genes were identified in the two Aquificales Sulfurihydrogenibium azorense and Persephonella marina. In contrast, neither of the two genes has been found in the sequenced genome of their close relative, Aquifex aeolicus. As in most bacteria, the rnpA genes of S. azorense and P. marina are preceded by the rpmH gene coding for ribosomal protein L34. This genetic region, including several genes up- and downstream of rpmH, is uniquely conserved among all three Aquificales strains, except that rnpA is missing in A. aeolicus. The RNase P RNAs (P RNAs) of S. azorense and P. marina are active catalysts that can be activated by heterologous bacterial P proteins at low salt. Although the two P RNAs lack helix P18 and thus one of the three major interdomain tertiary contacts, they are more thermostable than Escherichia coli P RNA and require higher temperatures for proper folding. Related to their thermostability, both RNAs include a subset of structural idiosyncrasies in their S domains, which were recently demonstrated to determine the folding properties of the thermostable S domain of Thermus thermophilus P RNA. Unlike 16S rRNA phylogeny that has placed the Aquificales as the deepest lineage of the bacterial phylogenetic tree, RNase P RNA-based phylogeny groups S. azorense and P. marina with the green sulfur, cyanobacterial, and delta/epsilon proteobacterial branches.     

Functional constraint and divergence in the G protein family in <Emphasis Type="Italic">Caenorhabditis elegans</Emphasis> and <Emphasis Type="Italic">Caenorhabditis briggsae</Emphasis>

Part of the challenge of the post-genomic world is to identify functional elements within the wide array of information generated by genome sequencing. Although cross-species comparisons and investigation of rates of sequence divergence are an efficient approach, the relationship between sequence divergence and functional conservation is not clear. Here, we use a comparative approach to examine questions of evolutionary rates and conserved function within the guanine nucleotide-binding protein (G protein) gene family in nematodes of the genus Caenorhabditis. In particular, we show that, in cases where the Caenorhabditis elegans ortholog shows a loss-of-function phenotype, G protein genes of C. elegans and Caenorhabditis briggsae diverge on average three times more slowly than G protein genes that do not exhibit any phenotype when mutated in C. elegans, suggesting that genes with loss of function phenotypes are subject to stronger selective constraints in relation to their function in both species. Our results also indicate that selection is as strong on G proteins involved in environmental perception as it is on those controlling other important processes. Finally, using phylogenetic footprinting, we identify a conserved non-coding motif present in multiple copies in the genomes of four species of Caenorhabditis. The presence of this motif in the same intron in the gpa-1 genes of C. elegans, C. briggsae and Caenorhabditis remanei suggests that it plays a role in the regulation of gpa-1, as well as other loci.Electronic Supplementary Material Supplementary material is available for this article at     

Multicellular Ascomycetous Fungal Genomes Contain More Than 8000 Genes

Fungi comprise a large monophyletic group of uni- and multicellular eukaryotic organisms in which many species are of economic or medical importance. Fungal genomes are variable in size (13–42 Mb), and multicellular species support true spatial and temporal cell-type-specific regulation of gene expression. In a 38.8-kbAspergillus nidulanscontiguous genomic DNA region, a transposable element and 12 potential genes were identified, 7 similar to genes in other organisms. This observation is consistent with the prediction that multicellular ascomycetous fungi harbor 8000–9000 genes in a 36-Mb average genome. Thus, the genomic DNA sequence of filamentous fungi will provide substantial amounts of genetic and functional information that is not available in yeast, for the human and other metazoan minimal gene complement.     

Flagellin gene sequence variation in the genus Pseudomonas

Flagellin gene (fliC) sequences from 18 strains of Pseudomonas sensu stricto representing 8 different species, and 9 representative fliC sequences from other members of the gamma sub-division of proteobacteria, were compared. Analysis was performed on N-terminal, C-terminal and whole fliC sequences. The fliC analyses confirmed the inferred relationship between P. mendocina, P. oleovorans and P. aeruginosa based on 16S rRNA sequence comparisons. In addition, the analyses indicated that P. putida PRS2000 was closely related to P. fluorescens SBW25 and P. fluorescens NCIMB 9046T, but suggested that P. putida PaW8 and P. putida PRS2000 were more closely related to other Pseudomonas spp. than they were to each other. There were a number of inconsistencies in inferred evolutionary relationships between strains, depending on the analysis performed. In particular, whole flagellin gene comparisons often differed from those obtained using N- and C-terminal sequences. However, there were also inconsistencies between the terminal region analyses, suggesting that phylogenetic relationships inferred on the basis of fliC sequence should be treated with caution. Although the central domain of fliC is highly variable between Pseudomonas strains, there was evidence of sequence similarities between the central domains of different Pseudomonas fliC sequences. This indicates the possibility of recombination in the central domain of fliC genes within Pseudomonas species, and between these genes and those from other bacteria.     

Incongruence between primary sequence data and the distribution of a mitochondrial atp1 group II intron among ferns and horsetails

Using DNA sequence data from multiple genes (often from more than one genome compartment) to reconstruct phylogenetic relationships has become routine. Augmenting this approach with genomic structural characters (e.g., intron gain and loss, changes in gene order) as these data become available from comparative studies already has provided critical insight into some long-standing questions about the evolution of land plants. Here we report on the presence of a group II intron located in the mitochondrial atp1 gene of leptosporangiate and marattioid ferns. Primary sequence data for the atp1 gene are newly reported for 27 taxa, and results are presented from maximum likelihood-based phylogenetic analyses using Bayesian inference for 34 land plants in three data sets: (1) single-gene mitochondrial atp1 (exon+intron sequences); (2) five combined genes (mitochondrial atp1 [exon only]; plastid rbcL, atpB, rps4; nuclear SSU rDNA); and (3) same five combined genes plus morphology. All our phylogenetic analyses corroborate results from previous fern studies that used plastid and nuclear sequence data: the monophyly of euphyllophytes, as well as of monilophytes; whisk ferns (Psilotidae) sister to ophioglossoid ferns (Ophioglossidae); horsetails (Equisetopsida) sister to marattioid ferns (Marattiidae), which together are sister to the monophyletic leptosporangiate ferns. In contrast to the results from the primary sequence data, the genomic structural data (atp1 intron distribution pattern) would seem to suggest that leptosporangiate and marattioid ferns are monophyletic, and together they are the sister group to horsetails--a topology that is rarely reconstructed using primary sequence data.     

Opa-like repeats in the genome of the Medfly Ceratitis capitata

The sequence determination of several genomic clones isolated from the Mediterranean fruitfly Ceratitis capitata identified the existence of opa-like repeats, often more than one being clustered in small chromosomal segments. These repeats have previously been shown to consist of stretches of tandemly reiterated glutamine-encoding residues, and they are found in multiple genes of several organisms. Most of the repeats described here are flanked or interrupted by stop codons in all reading frames and, thus, could not possibly be part of protein-coding sequences. Furthermore, these repeats, of which there are several hundred in the genome of the Medfly, can be used effectively for the determination of sequence polymorphisms, providing a convenient approach to obtain additional landmarks for the construction of genomic maps of this economically important insect.This paper is dedicated to the memory of our colleague and friend Dr. Jim Flach who took part in the initial phase of this work and died during the course of the investigation.     

Extending the Bacillus cereus group genomics to putative food-borne pathogens of different toxicity

The Bacillus cereus group represents sporulating soil bacteria containing pathogenic strains which may cause diarrheic or emetic food poisoning outbreaks. Multiple locus sequence typing revealed a presence in natural samples of these bacteria of about 30 clonal complexes. Application of genomic methods to this group was however biased due to the major interest for representatives closely related to Bacillus anthracis. Albeit the most important food-borne pathogens were not yet defined, existing data indicate that they are scattered all over the phylogenetic tree. The preliminary analysis of the sequences of three genomes discussed in this paper narrows down the gaps in our knowledge of the B. cereus group. The strain NVH391-98 is a rare but particularly severe food-borne pathogen. Sequencing revealed that the strain should be a representative of a novel bacterial species, for which the name Bacillus cytotoxis or Bacillus cytotoxicus is proposed. This strain has a reduced genome size compared to other B. cereus group strains. Genome analysis revealed absence of sigma B factor and the presence of genes encoding diarrheic Nhe toxin, not detected earlier. The strain B. cereus F837/76 represents a clonal complex close to that of B. anthracis. Including F837/76, three such B. cereus strains had been sequenced. Alignment of genomes suggests that B. anthracis is their common ancestor. Since such strains often emerge from clinical cases, they merit a special attention. The third strain, KBAB4, is a typical facultative psychrophile generally found in soil. Phylogenic studies show that in nature it is the most active group in terms of gene exchange. Genomic sequence revealed high presence of extra-chromosomal genetic material (about 530kb) that may account for this phenomenon. Genes coding Nhe-like toxin were found on a big plasmid in this strain. This may indicate a potential mechanism of toxicity spread from the psychrophile strain community. The results of this genomic work and ecological compartments of different strains incite to consider a necessity of creating prophylactic vaccines against bacteria closely related to NVH391-98 and F837/76. Presumably developing of such vaccines can be based on the properties of non-pathogenic strains such as KBAB4 or ATCC14579 reported here or earlier. By comparing the protein coding genes of strains being sequenced in this project to others we estimate the shared proteome, or core genome, in the B. cereus group to be 3000+/-200 genes and the total proteome, or pan-genome, to be 20-25,000 genes.     

Malate dehydrogenase: a useful phylogenetic marker for the genus Aeromonas

The reconstruction of correct genealogies among biological entities, the estimation of the divergence time between organisms or the study of the different events that occur along evolutionary lineages are not always based on suitable genes. For reliable results, it is necessary to look at full-length sequences of genes under stabilizing selection (neutral or purifying) and behaving as good molecular clocks. In bacteria it has been proved that the malate dehydrogenase gene (mdh) can be used to determine the inter- and intraspecies divergence, and hence this gene constitutes a potential marker for phylogeny and bacterial population genetics. We have sequenced the full-length mdh gene in 36 type and reference strains of Aeromonas. The species grouping obtained in the phylogenetic tree derived from mdh sequences was in agreement with that currently accepted for the genus Aeromonas. The maximum likelihood models applied to our sequences indicated that the mdh gene is highly conserved among the Aeromonas species and the main evolutionary force acting on it is purifying selection. Only two sites under potential diversifying selection were identified (T 108 and S 193). In order to determine if these two residues could have an influence on the MDH structure, we mapped them in a three-dimensional model constructed from the sequence of A. hydrophila using the human mitochondrial MDH as a template. The presence of purifying selection together with the linear relationship between substitutions and gene divergence makes the mdh an excellent candidate gene for a phylogeny of Aeromonas and probably for other bacterial groups.     

Identification of distinct genes with restricted expression in the somitic mesoderm in Xenopus embryo

We have used whole-mount in situ hybridisation to identify genes expressed in the somitic mesoderm during Xenopus early development. We report here the analysis of eight genes whose expression pattern has not been described previously. They include the Xenopus homologues of eukaryotic initiation factor 2beta, methionine adenosyltransferase II, serine dehydratase, alpha-adducin, oxoglutarate dehydrogenase, fragile X mental retardation syndrome related protein 1, monocarboxylate transporter and voltage-dependent anion channel 1. Interestingly, these genes exhibit very dynamic expression pattern during early development. At early gastrula stages several genes do not show localised expression pattern, while other genes are expressed in the marginal mesoderm or in ectoderm. As development proceeds, the expression of these genes is gradually restricted to different compartments of somite. This study thus reveals an unexpected dynamic expression pattern for various genes with distinct function in vertebrates.     

Speciation in <Emphasis Type="Italic">Chlamydia</Emphasis>: Genomewide Phylogenetic Analyses Identified a Reliable Set of Acquired Genes

Horizontal gene transfer (HGT), a process through which genomes acquire sequences from distantly related organisms, is believed to be a major source of genetic diversity in bacteria. A central question concerning the impact of HGT on bacterial genome evolution is the proportion of horizontally transferred sequences within genomes. This issue, however, remains unresolved because the various methods developed to detect potential HGT events identify different sets of genes. The present-day consensus is that phylogenetic analysis of individual genes is still the most objective and accurate approach for determining the occurrence and directionality of HGT. Here we present a genome-scale phylogenetic analysis of protein-encoding genes from five closely related Chlamydia, identifying a reliable set of sequences that have arisen via HGT since the divergence of the Chlamydia lineage. According to our knowledge, this is the first systematic phylogenetic inference-based attempt to establish a reliable set of acquired genes in a bacterial genome. Although Chlamydia are obligate intracellular parasites of higher eukaryotes, and thus suspected to be isolated from HGT more than the free-living species, our results show that their diversification has involved the introduction of foreign sequences into their genome. Furthermore, we also identified a complete set of genes that have undergone deletion, duplication, or rearrangement during this evolutionary period leading to the radiation of Chlamydia species. Our analysis may provide a deeper insight into how these medically important pathogens emerged and evolved from a common ancestor.     

Egg production and fertility in Drosophila depend upon the number of yolk-protein gene copies

Summary The yolk proteins of Drosophila melanogaster comprise a family of three related yolk polypeptides each encoded by a single-copy gene. We show by genetic crosses that each gene makes an equivalent contribution to the fecundity and fertility of the female and they do not individually provide unique functions to the embryo. We show that the number of eggs laid by a female depends upon the number of genes encoding yolk polypeptides present in the genome and furthermore that the probability of an egg hatching into an adult also critically depends upon the number of yolk protein genes present in the mother. This suggests that the three yolk protein-encoding genes in Drosophila melanogaster may have arisen by duplication, then been maintained for quantitative reasons because they increased egg production and fertility, rather than each protein evolving a different function as is the case with most small gene families, such as tubulins and collagen genes.     

Diversity and evolution of 5S rRNA gene family organization in Pythium

The 5S rRNA gene family organization among 87 species and varieties of Pythium was investigated to assess evolutionary stability of the two patterns detected and to determine which pattern is likely the ancestral state in the genus. Species with filamentous sporangia (Groups A-C according to the ITS phylogenetic tree for Pythium) had 5S genes linked to the rDNA repeat that were predominantly coded for on the DNA strand opposite to the one with the other rRNA genes (‘inverted’ orientation). A small group of species with contiguous sporangia (Group D) is related to Groups A-C but had unlinked 5S genes. The main group of species with spherical zoosporangia (Groups E-J) generally had unlinked 5S genes in tandem arrays. The six species in Group K, although they also have spherical sporangia, had linked genes on the same strand as the other rRNA genes ‘non-inverted’ and most of them had pairs of tandem 5S genes. The evolutionary stability of 5S sequence organization was compared with the stability of morphological characters as interpreted from a phylogeny based on ITS sequence analysis. Features of 5S sequence organization were found to be just as consistent within groups as were the morphological characters. To determine the ancestral type of 5S family organization, a survey of Phytophthora strains was conducted to supply an outgroup reference. The most parsimonious interpretation of the data in this survey yielded the tentative conclusion that the linked condition of the 5S sequences was ancestral.     

The complete nucleotide sequence and gene organization of carp (Cyprinus carpio) mitochondrial genome

The complete sequence of the carp mitochondrial genome of 16,575 base pairs has been determined. The carp mitochondrial genome encodes the same set of genes (13 proteins, 2 rRNAs, and 22 tRNAs) as do other vertebrate mitochondrial DNAs. Comparison of this teleostean mitochondrial genome with those of other vertebrates reveals a similar gene order and compact genomic organization. The codon usage of proteins of carp mitochondrial genome is similar to that of other vertebrates. The phylogenetic relationship for mitochondrial protein genes is more apparent than that for the mitochondrial tRNA and rRNA genes.Correspondence to: F. Huang