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.     

The unusual cell biology of the hyperthermophilic Crenarchaeon Ignicoccus hospitalis

The Crenarchaeon Ignicoccus hospitalis is an anaerobic, obligate chemolithoautotrophic hyperthermophile, growing by reduction of elemental sulfur using molecular hydrogen as electron donor. Together with Nanoarchaeum equitans it forms a unique, archaeal biocoenosis, in which I. hospitalis serves as host for N. equitans. Both organisms can be cultivated in a stable coculture which is mandatory for N. equitans but not for I. hospitalis. This strong dependence is affirmed by the fact that N. equitans obtains its lipids and amino acids from the host. I. hospitalis cells exhibit several unique features: they can adhere to surfaces by extracellular appendages (‘fibers’) which are not used for motility; they use a novel CO₂ fixation pathway, the dicarboxylate/4-hydroxybutyrate pathway; and they exhibit a unique cell envelope for Archaea consisting of two membranes but lacking an S-layer. These membranes form two cell compartments, a tightly packed cytoplasm surrounded by a weakly staining intermembrane compartment (IMC) with a variable width from 20 to 1,000 nm. In this IMC, many round or elongated vesicles are found which may function as carriers of lipids or proteins out of the cytoplasm. Based on immuno-EM analyses and immuno-fluorescence experiments it was demonstrated recently that the A₁A_O ATP synthase, the H₂:sulfur oxidoreductase complex and the acetyl-CoA synthetase (ACS) of I. hospitalis are located in its outermost membrane. Therefore, this membrane is energized and is here renamed as “outer cellular membrane” (OCM). Among all prokaryotes possessing two membranes in their cell envelope, I. hospitalis is the first organism with an energized outermost membrane and ATP synthesis outside the cytoplasm. Since DNA and ribosomes are localized in the cytoplasm, energy conservation is separated from information processing and protein biosynthesis in I. hospitalis. This raises questions concerning the function and characterization of the two membranes, the two cell compartments and of a possible ATP transfer to N. equitans.     

Rätselhafte Lebensgemeinschaft im Reich der Archaea. Die Feuerkugel und ihr Urzwerg

An enigmatical association of two Archaea The obligate anaerobic hyperthermophilic Nanoarchaeum equitans and Ignicoccus hospitalis represent a unique, purely archaeal biocoenosis which is mandatory for N. equitans. Its strong dependence on I. hospitalis is affirmed by the fact that its lipids and amino acids are obtained exclusively from the host. The Crenarchaeon I. hospitalis is characterized by energy production via reduction of elemental sulfur with molecular hydrogen and a novel CO₂‐fixation pathway. It possesses a unique cell envelope for Archaea with an inner and an outer membrane, forming two cell compartments, the cytoplasm and a huge intermembrane compartment. By immuno‐analyses we demonstrated that the ATP synthase and H₂:sulfur oxidoreductase complexes of I. hospitalis are located in the outer membrane. Thus I. hospitalis is the first Prokaryote with an energized outer membrane, ATP synthesis outside the cytoplasm, and spatial separation of energy conservation from information processing and protein biosynthesis. This raises many questions on the function and characterization of the two membranes, the two cell compartments, and a possible ATP transfer to N. equitans.     

Characterization of two key enzymes for aromatic amino acid biosynthesis in symbiotic archaea

Biosynthesis of L-tyrosine (L-Tyr) and L-phenylalanine (L-Phe) is directed by the interplay of three enzymes. Chorismate mutase (CM) catalyzes the rearrangement of chorismate to prephenate, which can be either converted to hydroxyphenylpyruvate by prephenate dehydrogenase (PD) or to phenylpyruvate by prephenate dehydratase (PDT). This work reports the first characterization of a trifunctional PD-CM-PDT from the smallest hyperthermophilic archaeon Nanoarchaeum equitans and a bifunctional CM-PD from its host, the crenarchaeon Ignicoccus hospitalis. Hexa-histidine tagged proteins were expressed in Escherichia coli and purified by affinity chromatography. Specific activities determined for the trifunctional enzyme were 21, 80, and 30 U/mg for CM, PD, and PDT, respectively, and 47 and 21 U/mg for bifunctional CM and PD, respectively. Unlike most PDs, these two archaeal enzymes were insensitive to regulation by L-Tyr and preferred NADP⁺ to NAD⁺ as a cofactor. Both the enzymes were highly thermally stable and exhibited maximal activity at 90 °C. N. equitans PDT was feedback inhibited by L-Phe (K_i = 0.8 µM) in a non-competitive fashion consistent with L-Phe’s combination at a site separate from that of prephenate. Our results suggest that PD from the unique symbiotic archaeal pair encompass a distinct subfamily of prephenate dehydrogenases with regard to their regulation and co-substrate specificity.     

Composition of the lipids of <Emphasis Type="Italic">Nanoarchaeum equitans</Emphasis> and their origin from its host <Emphasis Type="Italic">Ignicoccus</Emphasis> sp. strain KIN4/I

The contents and nature of the membrane lipids of Nanoarchaeum equitans and Ignicoccus sp. strain KIN4/I, grown at 90°C, and Ignicoccus sp. strain KIN4/I, cultivated at its lowest and highest growth temperatures (75°C and 95°C) were analyzed. Both organisms contained very simple and qualitatively identical assemblages of glycerol ether lipids, showing only differences in the amounts of certain components. LC–MS analyses of the total lipid extracts revealed that archaeol and caldarchaeol were the main core lipids. The predominant polar headgroups consisted of one or more sugar residues attached either directly to the core lipid or via a phosphate group. GC–MS analyses of hydrolyzed total lipid extracts revealed that the co-culture of N. equitans and Ignicoccus sp. strain KIN4/I, as well as Ignicoccus sp. strain KIN4/I grown at 90°C, contained phytane and biphytane in a ratio of approximately 4:1. Purified N. equitans cells and Ignicoccus sp. strain KIN4/I cultivated at 75°C and 95°C had a phytane to biphytane ratio of 10:1. Sugar residues were mainly mannose and small amounts of glucose. Consistent ¹³C fractionation patterns of isoprenoid chains of N. equitans and its host indicated that the N. equitans lipids were synthesized in the host cells.     

NotI clones in the analysis of the human genome

NotI linking clones contain sequences flanking NotI recognition sites and were previously shown to be tightly associated with CpG islands and genes. To directly assess the value of NotI clones in genome research, high density grids with 50 000 NotI linking clones originating from six representative NotI linking libraries were constructed. Altogether, these libraries contained nearly 100 times the total number of NotI sites in the human genome. A total of 3437 sequences flanking NotI sites were generated. Analysis of 3265 unique sequences demonstrated that 51% of the clones displayed significant protein similarity to SWISSPROT and TREMBL database proteins based on MSPcrunch filtering with stringent parameters. Of the 3265 sequences, 1868 (57.2%) were new sequences, not present in the EMBL and EST databases (similarity ≤ 90%). Among these new sequences, 795 (24.3%) showed similarity to known proteins and 712 (21.8%) displayed an identity of >75% at the nucleotide level to sequences from EMBL or EST databases. The remaining 361 (11.1%) sequences were completely new, i.e. <75% identical. The work also showed tight, specific association of NotI sites with the first exon and suggest that the so-called 3′ ESTs can actually be generated from 5′-ends of genes that contain NotI sites in their first exon.     

Characterization of a Single-Stranded DNA-Binding-Like Protein from Nanoarchaeum equitans—A Nucleic Acid Binding Protein with Broad Substrate Specificity

BackgroundSSB (single-stranded DNA-binding) proteins play an essential role in all living cells and viruses, as they are involved in processes connected with ssDNA metabolism. There has recently been an increasing interest in SSBs, since they can be applied in molecular biology techniques and analytical methods. Nanoarchaeum equitans, the only known representative of Archaea phylum Nanoarchaeota, is a hyperthermophilic, nanosized, obligatory parasite/symbiont of Ignicoccus hospitalis.ResultsThis paper reports on the ssb-like gene cloning, gene expression and characterization of a novel nucleic acid binding protein from Nanoarchaeum equitans archaeon (NeqSSB-like protein). This protein consists of 243 amino acid residues and one OB fold per monomer. It is biologically active as a monomer like as SSBs from some viruses. The NeqSSB-like protein displays a low sequence similarity to the Escherichia coli SSB, namely 10% identity and 29% similarity, and is the most similar to the Sulfolobus solfataricus SSB (14% identity and 32% similarity). The NeqSSB-like protein binds to ssDNA, although it can also bind mRNA and, surprisingly, various dsDNA forms, with no structure-dependent preferences as evidenced by gel mobility shift assays. The size of the ssDNA binding site, which was estimated using fluorescence spectroscopy, is 7±1 nt. No salt-dependent binding mode transition was observed. NeqSSB-like protein probably utilizes a different model for ssDNA binding than the SSB proteins studied so far. This protein is highly thermostable; the half-life of the ssDNA binding activity is 5 min at 100°C and melting temperature (T_m) is 100.2°C as shown by differential scanning calorimetry (DSC) analysis.ConclusionNeqSSB-like protein is a novel highly thermostable protein which possesses a unique broad substrate specificity and is able to bind all types of nucleic acids.     

Genome Sequence of the Chemolithoautotrophic Nitrite-Oxidizing Bacterium Nitrobacter winogradskyi Nb-255

The alphaproteobacterium Nitrobacter winogradskyi (ATCC 25391) is a gram-negative facultative chemolithoautotroph capable of extracting energy from the oxidation of nitrite to nitrate. Sequencing and analysis of its genome revealed a single circular chromosome of 3,402,093 bp encoding 3,143 predicted proteins. There were extensive similarities to genes in two alphaproteobacteria, Bradyrhizobium japonicum USDA110 (1,300 genes) and Rhodopseudomonas palustris CGA009 CG (815 genes). Genes encoding pathways for known modes of chemolithotrophic and chemoorganotrophic growth were identified. Genes encoding multiple enzymes involved in anapleurotic reactions centered on C₂ to C₄ metabolism, including a glyoxylate bypass, were annotated. The inability of N. winogradskyi to grow on C₆ molecules is consistent with the genome sequence, which lacks genes for complete Embden-Meyerhof and Entner-Doudoroff pathways, and active uptake of sugars. Two gene copies of the nitrite oxidoreductase, type I ribulose-1,5-bisphosphate carboxylase/oxygenase, cytochrome c oxidase, and gene homologs encoding an aerobic-type carbon monoxide dehydrogenase were present. Similarity of nitrite oxidoreductases to respiratory nitrate reductases was confirmed. Approximately 10% of the N. winogradskyi genome codes for genes involved in transport and secretion, including the presence of transporters for various organic-nitrogen molecules. The N. winogradskyi genome provides new insight into the phylogenetic identity and physiological capabilities of nitrite-oxidizing bacteria. The genome will serve as a model to study the cellular and molecular processes that control nitrite oxidation and its interaction with other nitrogen-cycling processes.     

Complete Genome Sequence of the Plant Pathogen Erwinia amylovora Strain ATCC 49946

Erwinia amylovora causes the economically important disease fire blight that affects rosaceous plants, especially pear and apple. Here we report the complete genome sequence and annotation of strain ATCC 49946. The analysis of the sequence and its comparison with sequenced genomes of closely related enterobacteria revealed signs of pathoadaptation to rosaceous hosts.Erwinia amylovora, a plant-associated member of the Enterobacteriaceae, causes fire blight, a devastating disease of rosaceous plants, especially pear and apple (6). The complete genome of Ea273 (ATCC 49946), a virulent strain isolated from an infected apple tree in New York State, was sequenced. Total DNA was extracted and prepared in pMAQ1 shotgun libraries. The complete shotgun sequence was obtained by using dye terminator chemistry in ABI 3730 automated sequencers and contains 88,457 reads (11.12-fold coverage), yielding a theoretical coverage of the genome of 99.99%. The sequence was assembled, finished, and annotated as described previously (1, 5), using Artemis (4) to collate data and facilitate annotation.The genome of E. amylovora consists of a circular chromosome of 3,805,874 bp and two plasmids, AMYP1 (28,243 bp) and AMYP2 (71,487 bp). Coding regions in the chromosome account for 85.1% of the total sequence, with 3,483 identified coding sequences (CDS). Two hundred fifty-four (7%) of the CDSs do not have any matches in current NCBI databases; 114 (3.3%) correspond to conserved hypothetical proteins. Forty-nine CDSs (1.4%) are similar to genes from mobile elements such as integrases, transposases, and bacteriophages, and 110 CDSs (3.2%) were classified as pseudogenes due to interruptions or truncations of the CDSs. The remaining 2,956 annotated CDSs include among other categories genes involved in biosynthesis of the cellular envelope and modifications of surface proteins (299 CDSs [11%]) and genes involved in signal transduction and regulation (228 CDSs [8%]). Seven rRNA operons and 78 tRNA sequences were identified in the chromosome; two new clusters were identified (AMY1550-1575 and AMY2648-2676) that resemble the T3SS-encoding SSR-1 island of Sodalis glossinidius (2), and four clusters that contain genes for biosynthesis of flagella, which based on their location might be regulated independently.The smaller plasmid, AMYP1, had been reported as pEA29 (3); its sequence is nearly identical to the one reported here. The larger plasmid, AMYP2, renamed pEA72 for consistency in nomenclature, contains 87 predicted CDSs, with two predicted mobile-element-related CDSs and one pseudogene. Among the CDSs with annotated functions are a cluster of genes (AMYP2_49 to AMYP2_62) that encode a putative type IV fimbrial system (pil genes).The genome of E. amylovora is only 3.8 Mb long, whereas most free-living enterobacteria, including plant pathogens, have genomes of 4.5 Mb to 5.5 Mb. Comparison of the genome of Ea273 with the sequenced genomes of 15 closely related enterobacteria identified 21 lineage-specific regions, which might be considered genomic islands. E. amylovora has many more predicted pseudogenes, relative to other enterobacteria with similar lifestyles. Given its size and the preponderance of pseudogenes, genome reduction may have occurred via mutational inactivation and subsequent deletion with the following consequences: E. amylovora has fewer genes involved in anaerobic respiration and fermentation than are found in typical related enterobacteria; this likely result in a reduced capacity to live in anaerobic environments.The genome sequence of E. amylovora has revealed clear signs of pathoadaptation to the rosaceous plant environment. For example, T3SS-related proteins are present that are more similar to proteins of other plant pathogens than to proteins of closely related enterobacteria. These include type III effectors, homologous to those of plant-pathogenic pseudomonads, which confer virulence to E. amylovora in plants, and a sorbitol-metabolizing cluster that may confer a competitive advantage for survival in rosaceous plants. The reduced genome size and erosion or loss of genes involved in anaerobic respiration and nitrate assimilation are remarkable, relative to other plant- and animal-pathogenic members of the Enterobacteriaceae.     

Comparative proteomic analysis of a sea urchin (Heliocidaris erythrogramma) antibacterial response revealed the involvement of apextrin and calreticulin

Echinoderms evolved early in the deuterostome lineage, and as such constitute model organisms for comparative physiology and immunology. The sea urchin genome sequence (Strongylocentrotus purpuratus) revealed a complex repertoire of genes with similarities to the immune response genes of other species. To complement these genomic data, we investigated the responses of sea urchins to the injection of bacteria using a comparative proteomics approach on a closely related species. In the sea urchin, Heliocidaris erythrogramma, the relative abundance of many proteins was altered in response to the injection of both bacteria and saline, suggesting their involvement in wounding responses, while others were differentially altered in response to bacteria only. The identities of 15 proteins that differed in relative abundance were determined by mass spectrometry. These proteins revealed a significant modification in energy metabolism in coelomocytes towards the consumption of glutamate and the production of NADPH after injection, as well as an increased concentration of cell signalling molecules, such as heterotrimeric guanine nucleotide-binding protein. The injection of bacteria specifically increased the abundance of apextrin and calreticulin, suggesting that these two proteins are involved in the sequestration or inactivation of bacteria.     

Development and Evaluation of Functional Gene Arrays for Detection of Selected Genes in the Environment

To determine the potential of DNA array technology for assessing functional gene diversity and distribution, a prototype microarray was constructed with genes involved in nitrogen cycling: nitrite reductase (nirS and nirK) genes, ammonia mono-oxygenase (amoA) genes, and methane mono-oxygenase (pmoA) genes from pure cultures and those cloned from marine sediments. In experiments using glass slide microarrays, genes possessing less than 80 to 85% sequence identity were differentiated under hybridization conditions of high stringency (65°C). The detection limit for nirS genes was approximately 1 ng of pure genomic DNA and 25 ng of soil community DNA using our optimized protocol. A linear quantitative relationship (r² = 0.89 to 0.94) was observed between signal intensity and target DNA concentration over a range of 1 to 100 ng for genomic DNA (or genomic DNA equivalent) from both pure cultures and mixed communities. However, the quantitative capacity of microarrays for measuring the relative abundance of targeted genes in complex environmental samples is less clear due to divergent target sequences. Sequence divergence and probe length affected hybridization signal intensity within a certain range of sequence identity and size, respectively. This prototype functional gene array did reveal differences in the apparent distribution of nir and amoA and pmoA gene families in sediment and soil samples. Our results indicate that glass-based microarray hybridization has potential as a tool for revealing functional gene composition in natural microbial communities; however, more work is needed to improve sensitivity and quantitation and to understand the associated issue of specificity.     

Analysis of complete genome sequence of Neorickettsia risticii: causative agent of Potomac horse fever

Neorickettsia risticii is an obligate intracellular bacterium of the trematodes and mammals. Horses develop Potomac horse fever (PHF) when they ingest aquatic insects containing encysted N. risticii-infected trematodes. The complete genome sequence of N. risticii Illinois consists of a single circular chromosome of 879 977 bp and encodes 38 RNA species and 898 proteins. Although N. risticii has limited ability to synthesize amino acids and lacks many metabolic pathways, it is capable of making major vitamins, cofactors and nucleotides. Comparison with its closely related human pathogen N. sennetsu showed that 758 (88.2%) of protein-coding genes are conserved between N. risticii and N. sennetsu. Four-way comparison of genes among N. risticii and other Anaplasmataceae showed that most genes are either shared among Anaplasmataceae (525 orthologs that generally associated with housekeeping functions), or specific to each genome (>200 genes that are mostly hypothetical proteins). Genes potentially involved in the pathogenesis of N. risticii were identified, including those encoding putative outer membrane proteins, two-component systems and a type IV secretion system (T4SS). The bipolar localization of T4SS pilus protein VirB2 on the bacterial surface was demonstrated for the first time in obligate intracellular bacteria. These data provide insights toward genomic potential of N. risticii and intracellular parasitism, and facilitate our understanding of PHF pathogenesis.     

GST-PRIME: a genome-wide primer design software for the generation of gene sequence tags

The availability of sequenced genomes has generated a need for experimental approaches that allow the simultaneous analysis of large, or even complete, sets of genes. To facilitate such analyses, we have developed GST-PRIME, a software package for retrieving and assembling gene sequences, even from complex genomes, using the NCBI public database, and then designing sets of primer pairs for use in gene amplification. Primers were designed by the program for the direct amplification of gene sequence tags (GSTs) from either genomic DNA or cDNA. Test runs of GST-PRIME on 2000 randomly selected Arabidopsis and Drosophila genes demonstrate that 93 and 88% of resulting GSTs, respectively, fulfilled imposed length criteria. GST-PRIME primer pairs were tested on a set of 1900 Arabidopsis genes coding for chloroplast-targeted proteins: 95% of the primer pairs used in PCRs with genomic DNA generated the correct amplicons. GST-PRIME can thus be reliably used for large-scale or specific amplification of intron-containing genes of multicellular eukaryotes.     

Expression profiling of major histocompatibility and natural killer complex genes reveals candidates for controlling risk of graft versus host disease

Background

The major histocompatibility complex (MHC) is the most important genomic region that contributes to the risk of graft versus host disease (GVHD) after haematopoietic stem cell transplantation. Matching of MHC class I and II genes is essential for the success of transplantation. However, the MHC contains additional genes that also contribute to the risk of developing acute GVHD. It is difficult to identify these genes by genetic association studies alone due to linkage disequilibrium in this region. Therefore, we aimed to identify MHC genes and other genes involved in the pathophysiology of GVHD by mRNA expression profiling.

Methodology/Principal Findings

To reduce the complexity of the task, we used genetically well-defined rat inbred strains and a rat skin explant assay, an in-vitro-model of the graft versus host reaction (GVHR), to analyze the expression of MHC, natural killer complex (NKC), and other genes in cutaneous GVHR. We observed a statistically significant and strong up or down regulation of 11 MHC, 6 NKC, and 168 genes encoded in other genomic regions, i.e. 4.9%, 14.0%, and 2.6% of the tested genes respectively. The regulation of 7 selected MHC and 3 NKC genes was confirmed by quantitative real-time PCR and in independent skin explant assays. In addition, similar regulations of most of the selected genes were observed in GVHD-affected skin lesions of transplanted rats and in human skin explant assays.

Conclusions/Significance

We identified rat and human MHC and NKC genes that are regulated during GVHR in skin explant assays and could therefore serve as biomarkers for GVHD. Several of the respective human genes, including HLA-DMB, C2, AIF1, SPR1, UBD, and OLR1, are polymorphic. These candidates may therefore contribute to the genetic risk of GVHD in patients.