Genetic and functional properties of uncultivated thermophilic crenarchaeotes from a subsurface gold mine as revealed by analysis of genome fragments

Within a phylum Crenarchaeota, only some members of the hyperthermophilic class Thermoprotei, have been cultivated and characterized. In this study, we have constructed a metagenomic library from a microbial mat formation in a subsurface hot water stream of the Hishikari gold mine, Japan, and sequenced genome fragments of two different phylogroups of uncultivated thermophilic Crenarchaeota: (i) hot water crenarchaeotic group (HWCG) I (41.2 kb), and (ii) HWCG III (49.3 kb). The genome fragment of HWCG I contained a 16S rRNA gene, two tRNA genes and 35 genes encoding proteins but no 23S rRNA gene. Among the genes encoding proteins, several genes for putative aerobic-type carbon monoxide dehydrogenase represented a potential clue with regard to the yet unknown metabolism of HWCG I Archaea. The genome fragment of HWCG III contained a 16S/23S rRNA operon and 44 genes encoding proteins. In the 23S rRNA gene, we detected a homing-endonuclease encoding a group I intron similar to those detected in hyperthermophilic Crenarchaeota and Bacteria, as well as eukaryotic organelles. The reconstructed phylogenetic tree based on the 23S rRNA gene sequence reinforced the intermediate phylogenetic affiliation of HWCG III bridging the hyperthermophilic and non-thermophilic uncultivated Crenarchaeota.     

Characterization of the rRNA genes of Ehrlichia chaffeensis and Anaplasma phagocytophila

The rRNA genes of Ehrlichia chaffeensis and Anaplasma phagocytophila have been analyzed. The 16S rRNA genes were previously characterized for both of these agents. Southern hybridization was used to show that there are single copies of both the 16S and 23S rRNA genes in the genomes of each organism, and that the 16S rRNA genes were upstream from the 23S rRNA genes by at least 16 and 11 Kb for E. chaffeensis and A. phagocytophila, respectively. PCR amplification and gene walking was used to sequence the 23S and 5S rRNA genes, and show that these genes are contiguous and are likely expressed as a single operon. The level of homology between the E. chaffeensis and A. phagocytophila 23S and 5S rRNA genes, and 23S-5S spacers, was 91.8, 81.5, and 40%, respectively. To confirm the hybridization data, genome walking was used to sequence downstream of the 16S rRNA genes, and although no tRNA genes were identified, open reading frames encoding homologues of the Escherichia coli succinate dehydrogenase, subunit C, were found in both E. chaffeensis and A. phagocytophila. Phylogenetic analysis using the 23S rRNA gene suggests that reorganization of the phylum Proteobacteria by division of the class Alphaproteobacteria into two separate subclasses, may be appropriate.     

Ribosomal RNA and ribosomal proteins in corynebacteria

Ribosomal RNAs (rRNAs) (16S, 23S, 5S) encoded by the rrn operons and ribosomal proteins play a very important role in the formation of ribosomes and in the control of translation. Five copies of the rrn operon were reported by hybridization studies in Brevibacterium (Corynebacterium) lactofermentum but the genome sequence of Corynebacterium glutamicum provided evidence for six rrn copies. All six copies of the C. glutamicum 16S rRNA have a size of 1523 bp and each of the six copies of the 5S contain 120 bp whereas size differences are found between the six copies of the 23S rRNA. The anti-Shine-Dalgarno sequence at the 3'-end of the 16S rRNA was 5'-CCUCCUUUC-3'. Each rrn operon is transcribed as a large precursor rRNA (pre-rRNA) that is processed by RNaseIII and other RNases at specific cleavage boxes that have been identified in the C. glutamicum pre-rRNA. A secondary structure of the C. glutamicum 16S rRNA is proposed. The 16S rRNA sequence has been used as a molecular evolution clock allowing the deduction of a phylogenetic tree of all Corynebacterium species. In C. glutamicum, there are 11 ribosomal protein gene clusters encoding 42 ribosomal proteins. The organization of some of the ribosomal protein gene cluster is identical to that of Escherichia coli whereas in other clusters the organization of the genes is rather different. Some specific ribosomal protein genes are located in a different cluster in C. glutamicum when compared with E. coli, indicating that the control of expression of these genes is different in E. coli and C. glutamicum.     

Comparative analysis of a genome fragment of an uncultivated mesopelagic crenarchaeote reveals multiple horizontal gene transfers

Marine planktonic crenarchaeota have escaped all cultivation attempts to date, all crenarchaeota growing in pure culture so far being hyperthermophiles. Here, we present a comparative genomic analysis of a 16S- plus 23S-rDNA-containing fragment of a crenarchaeote retrieved from an environmental genomic library constructed from picoplankton collected at 500 m depth in the Antarctic Polar Front. The clone DeepAnt-EC39 contained an insert of 33.3 kbp, which was completely sequenced. DeepAnt-EC39 appears to represent a lineage specific to deep-sea waters but widespread geographically, as revealed by the analysis of the 16S-23S-rDNA intergenic spacer region. A comparison with previously sequenced marine crenarchaeotal genomic clones also containing an rrn operon (74A4, 4B7 and Cenarchaeum symbiosum strains A and B) revealed a highly variable structure involving gene rearrangements and insertions/deletions. The surroundings of the rrn operon and the contiguous glutamate-1-semialdehyde aminotransferase gene appear hot spots for recombination. Phylogenetic analyses of all individual predicted proteins revealed the existence of several likely cases of horizontal gene transfer both, between the two archaeal kingdoms and between the two prokaryotic domains. The most frequent horizontal transfers appear to involve genes from mesophilic methanogenic euryarchaeota related to Methanosarcinales. We hypothesise that the acquisition of genes from mesophilic bacteria and euryarchaeota has played a major role in the adaptation of Group I crenarchaeota to life at lower temperatures.     

Nucleotide sequence of a 'truncated rRNA operon' of the Euglena gracilis chloroplast genome.

An extra 16S rRNA gene (s-16S rDNA) from the Euglena gracilis chloroplast genome and several hundred positions of its flanking regions have been sequenced. The structural part has 1486 positions and is to 98% homologous in its sequence with the 16S rRNA gene in functional chloroplast rRNA operons. Sequences of about 200 positions upstream and 15 positions downstream of the structural part of the s-16S rRNA gene region are highly homologous with corresponding parts in the functional operon. Neither tRNA genes (A1a, I1e) nor parts of the 23S and 5S rRNA genes are found within 557 positions after the 3' end of the s-16S rRNA gene, i.e., the 330 bp homology, observed in electron microscopic studies of heteroduplexes (4), between the s-16S rDNA downstream region and the 6.2 kb repeated segment containing the functional rRNA operon, must be due to a DNA stretch in the interoperon spacer. A structural model of the "truncated rRNA operon" is presented. Results from S-1 endonuclease analysis suggest that the s-16S rDNA region is probably not transcribed into stable s-16S rRNA.     

Molecular cloning and comparative sequence analyses of rRNA operons in Streptomyces nodosus ATCC 14899.

The genome of Streptomyces nodosus contains six ribosomal RNA (rRNA) operons. Four of the rRNA operons; rrnB, rrnD, rrnE and rrnF were cloned. We have completely sequenced all four operons, including a region 750 base pairs (bp) upstream of the 16S rRNA gene. The three rRNA genes present in each operon were closely linked in the order 16S-23S-5S. A sequence comparison of the four operons showed more than 99% sequence similarity between the corresponding 16S and 23S rRNA genes, and more than 97% similarity between 5S rRNA genes. The sequence differences observed between 23S rRNA genes appeared to be localized in two specific regions. Substantial sequence differences were found in the region upstream of the 16S rRNA gene as well as in the internal transcribed spacers. No tRNA gene was found in the 16S-23S spacer regions.     

Domain organization and phylogenetic analysis of the chitinase-like family of proteins in three species of insects

A bioinformatics-based investigation of three insect species with completed genome sequences has revealed that insect chitinase-like proteins (glycosylhydrolase family 18) are encoded by a rather large and diverse group of genes. We identified 16, 16 and 13 putative chitinase-like genes in the genomic databases of the red flour beetle, Tribolium castaneum, the fruit fly, Drosophila melanogaster, and the malaria mosquito, Anopheles gambiae, respectively. Chitinase-like proteins encoded by this gene family were classified into five groups based on phylogenetic analyses. Group I chitinases are secreted proteins that are the most abundant such enzymes in molting fluid and/or integument, and represent the prototype enzyme of the family, with a single copy each of the catalytic domain and chitin-binding domain (ChBD) connected by an S/T-rich linker polypeptide. Group II chitinases are unusually larger-sized secreted proteins that contain multiple catalytic domains and ChBDs. Group III chitinases contain two catalytic domains and are predicted to be membrane-anchored proteins. Group IV chitinases are the most divergent. They usually lack a ChBD and/or an S/T-rich linker domain, and are known or predicted to be secreted proteins found in gut or fat body. Group V proteins include the putative chitinase-like imaginal disc growth factors (IDGFs). In each of the three insect genomes, multiple genes encode group IV and group V chitinase-like proteins. In contrast, groups I-III are each represented by only a singe gene in each species.     

Restriction enzyme analysis of Bacillus subtilis ribosomal ribonucleic acid genes.

The organization of the ribosomal ribonucleic acid (rRNA) genes (rDNA) of Bacillus subtilis was examined by cleaving the genome with several restriction endonucleases. The rDNA sequences were assayed by hybridization with purified radioactive rRNA's. Our interpretation of the resulting electrophoretic patterns is strengthened by an analysis of a fragment of B. subtilis rDNA cloned in Escherichia coli. The results indicated that there are eight rRNA operons in B. subtilis. Each operon contains one copy of the sequences coding for 16S, 23S, and 5S rRNA. The sequences coding for 5S rRNA were shown to be more closely linked to the 23S rRNA genes than to the 16S rRNA genes.     

Acidobacteria form a coherent but highly diverse group within the bacterial domain: evidence from environmental genomics

Acidobacteria have been established as a novel phylum of Bacteria that is consistently detected in many different habitats around the globe by 16S rDNA-based molecular surveys. The phylogenetic diversity, ubiquity and abundance of this group, particularly in soil habitats, suggest an important ecological role and extensive metabolic versatility. However, the genetic and physiological information about Acidobacteria is scarce. In order to gain insight into genome structure, evolution and diversity of these microorganisms we have initiated an environmental genomic approach by constructing large insert libraries directly from DNA of a calcerous grassland soil. Genomic fragments of Acidobacteria were identified with specific 16S rDNA probes and sequence analyses of six independently identified clones were performed, representing in total more than 210,000 bp. The 16S rRNA genes of the genomic fragments differed between 2.3% and 19.9% and were placed into two different subgroups of Acidobacteria (groups III and V). Although partial co-linearity was found between genomic fragments, the gene content around the rRNA operons was generally not conserved. Phylogenetic reconstructions with orthologues that were encoded on two of the six genomic fragments (PurF, PurL, PurB and formamidopyrimidine-DNA glycosylase) confirmed the coherence of the acidobacterial phylum. One genomic fragment harboured a cluster of eight genes which was syntenic and highly homologous to genomic regions in Rhodopseudomonas palustris and Bradyrhizobium japonicum, including a conserved two-component system. Phylogenetic analysis of the putative response regulator confirmed that this similarity between Rhizobiales and Acidobacteria might be due to a horizontal gene transfer. In total, our data give first insight into the genome content and diversity of the ubiquitously distributed but poorly characterized phylum of Acidobacteria. Furthermore they support the phylogenetic inferences made from 16S rRNA gene libraries, suggesting that Acidobacteria form a broad group in the same sense and with a similar diversity as that of many well-studied bacterial phyla.     

Multiplex PCR assay for ureC and 16S rRNA genes clearly discriminates between both subspecies of Photobacterium damselae

A multiplex-PCR approach, employing 2 primer pairs directed to internal regions of the 16S rRNA and ureC genes, was utilized to analyze a collection of Photobacterium damselae strains, including 25 isolates of subspecies piscicida and 15 isolates of subspecies damselae. With this procedure, all the P. damselae subsp. damselae strains yielded 2 amplification products, one of 267 bp and the other of 448 bp, corresponding to internal fragments of the 16S rRNA and ureC genes, respectively. However, P. damselae subsp. piscicida isolates only showed the PCR product of 267 bp (16S rRNA fragment), indicating the absence of the urease gene in its genome. We have constructed a DNA probe directed to an internal region of the ureC gene, and corroborated by dot blot hybridization that the P. damselae subsp. piscicida lacks this gene, whereas it is present in the subspecies damselae. This constitutes the first successful discrimination between both subspecies using a PCR procedure, which could become a useful tool for diagnosis of pasteurellosis in the field. In addition, since these 2 subspecies have been shown to share nearly the same rrn operon sequence, our results provided evidence that one of the steps in the P. damselae speciation proccess included gain/loss events associated with the ure operon.     

Genes involved in the anaerobic degradation of ethylbenzene in a denitrifying bacterium,strain EbN1

Genes involved in anaerobic degradation of the petroleum hydrocarbon ethylbenzene in the denitrifying Azoarcus-like strain EbN1 were identified on a 56-kb DNA contig obtained from shotgun sequencing. Ethylbenzene is first oxidized via ethylbenzene dehydrogenase to (S)-1-phenylethanol; this is converted by (S)-1-phenylethanol dehydrogenase to acetophenone. Further degradation probably involves acetophenone carboxylase forming benzoylacetate, a ligase forming benzoylacetyl-CoA, and a thiolase forming acetyl-CoA and benzoyl-CoA. Genes of this pathway were identified via N-terminal sequences of proteins isolated from strain EbN1 and by sequence similarities to proteins from other bacteria. Ethylbenzene dehydrogenase is encoded by three genes (ebdABC), in accordance with the heterotrimeric enzyme structure. Binding domains for a molybdenum cofactor (in subunit EbdA) and iron/sulfur-clusters (in subunits EbdA and EbdB) were identified. The previously observed periplasmic location of the enzyme was corroborated by the presence of a twin-arginine leader peptide characteristic of the Tat system for protein export. A fourth gene (ebdD) was identified, the product of which may act as an enzyme-specific chaperone in the maturation of the molybdenum-containing subunit. A distinct gene (ped) coding for (S)-1-phenylethanol dehydrogenase apparently forms an operon with the ebdABCD genes. The ped gene product with its characteristic NAD(P)-binding motif in the N-terminal domain belongs to the short-chain dehydrogenase/reductase (SDR) superfamily. A further operon apparently contains five genes (apc1-5) suggested to code for subunits of acetophenone carboxylase. Four of the five gene products are similar to subunits of acetone carboxylase from Xanthobacter autotrophicus. Upstream of the apc genes, a single gene (bal) was identified which possibly codes for a benzoylacetate CoA-ligase and which is co-transcribed with the apc genes. In addition, an apparent operon containing almost all genes required for beta-oxidation of fatty acids was detected; one of the gene products may be involved in thiolytic cleavage of benzoylacetyl-CoA. The DNA fragment also included genes for regulatory systems; these were two sets of two-component systems, two LysR homologs, and a TetR homolog. Some of these proteins may be involved in ethylbenzene-dependent gene expression.     

Nucleotide sequence of the single ribosomal RNA operon of pea chloroplast DNA

The nucleotide sequence of an 8 kbp region of pea ( Pisum sativum L.) chloroplast DNA containing the rRNA operon and putative promoter sites has been determined and compared to the corresponding sequences from maize, tobacco and the liverwort Marchantia polymorpha . The chloroplast DNA species of all vascular plants investigated, with the exception of a few legumes including pea, and of Marchantia contain an inverted repeat with an rRNA operon. The pea rRNA operon is the first sequenced rRNA operon from a plant with only one copy of the rRNA genes per molecule of chloroplast DNA. The organization of the operon is the same as for maize, tobacco and Marchantia . i.e. tRNA-Val gene/16S rRNA gene/spacer with intron-containing genes for tRNA-Ile and tRNA-Ala/23S rRNA gene/4.5S rRNA gene/5S rRNA gene. Current evidence suggests that the tRNA-Val gene may not be contranscribed with the other genes. For pea 16S, 23S, 4.5S and 5S rRNA have 1488, 2813, 105 and 121 nucleotides, respectively. The homologies of the entire operon (the tRNA-Val gene - 5S rRNA region) to those from tobacco, maize and Marchantia are 88, 82 and 79%, respectively. The corresponding homologies for tobacco/maize, tobacco/ Marchantia and maize/ Marchantia have similar values. The 16S and 23S rRNA genes from pea are more than 90% homologous to those from the 3 other species. We conclude that the fact that pea only has one set of rRNA genes per molecule of chloroplast DNA is apparently not correlated with any significant difference between the pea operon and the rRNA operons from tobacco, maize and Marchantia .     

rrndb: the Ribosomal RNA Operon Copy Number Database

The Ribosomal RNA Operon Copy Number Database (rrndb) is an Internet-accessible database containing annotated information on rRNA operon copy number among prokaryotes. Gene redundancy is uncommon in prokaryotic genomes, yet the rRNA genes can vary from one to as many as 15 copies. Despite the widespread use of 16S rRNA gene sequences for identification of prokaryotes, information on the number and sequence of individual rRNA genes in a genome is not readily accessible. In an attempt to understand the evolutionary implications of rRNA operon redundancy, we have created a phylogenetically arranged report on rRNA gene copy number for a diverse collection of prokaryotic microorganisms. Each entry (organism) in the rrndb contains detailed information linked directly to external websites including the Ribosomal Database Project, GenBank, PubMed and several culture collections. Data contained in the rrndb will be valuable to researchers investigating microbial ecology and evolution using 16S rRNA gene sequences. The rrndb web site is directly accessible on the WWW at http://rrndb.cme. msu.edu.     

Toxicity phenotype does not correlate with phylogeny of Cylindrospermopsis raciborskii strains

Cylindrospermopsis raciborskii is a species of freshwater, bloom-forming cyanobacterium. C. raciborskii produces toxins, including cylindrospermopsin (hepatotoxin) and saxitoxin (neurotoxin), although non toxin-producing strains are also observed. In spite of differences in toxicity, C. raciborskii strains comprise a monophyletic group, based upon 16S rRNA gene sequence identities (greater than 99%). We performed phylogenetic analyses; 16S rRNA gene and 16S-23S rRNA gene internally transcribed spacer (ITS-1) sequence comparisons, and genomic DNA restriction fragment length polymorphism (RFLP), resolved by pulsed-field gel electrophoresis (PFGE), of strains of C. raciborskii, obtained mainly from the Australian phylogeographic cluster. Our results showed no correlation between toxic phenotype and phylogenetic association in the Australian strains. Analyses of the 16S rRNA gene and the respective ITS-1 sequences (long L, and short S) showed an independent evolution of each ribosomal operon. The genes putatively involved in the cylindrospermopsin biosynthetic pathway were present in one locus and only in the hepatotoxic strains, demonstrating a common genomic organization for these genes and the absence of mutated or inactivated biosynthetic genes in the non toxic strains. In summary, our results support the hypothesis that the genes involved in toxicity may have been transferred as an island by processes of gene lateral transfer, rather than convergent evolution.