Genome-wide identification of transcriptional start sites in the haloarchaeon Haloferax volcanii based on differential RNA-Seq (dRNA-Seq)

BACKGROUND

Differential RNA-Seq (dRNA-Seq) is a recently developed method of performing primary transcriptome analyses that allows for the genome-wide mapping of transcriptional start sites (TSSs) and the identification of novel transcripts. Although the transcriptomes of diverse bacterial species have been characterized by dRNA-Seq, the transcriptome analysis of archaeal species is still rather limited. Therefore, we used dRNA-Seq to characterize the primary transcriptome of the model archaeon Haloferax volcanii.

RESULTS

Three independent cultures of Hfx. volcanii grown under optimal conditions to the mid-exponential growth phase were used to determine the primary transcriptome and map the 5'-ends of the transcripts. In total, 4749 potential TSSs were detected. A position weight matrix (PWM) was derived for the promoter predictions, and the results showed that 64 % of the TSSs were preceded by stringent or relaxed basal promoters. Of the identified TSSs, 1851 belonged to protein-coding genes. Thus, fewer than half (46 %) of the 4040 protein-coding genes were expressed under optimal growth conditions. Seventy-two percent of all protein-coding transcripts were leaderless, which emphasized that this pathway is the major pathway for translation initiation in haloarchaea. A total of 2898 of the TSSs belonged to potential non-coding RNAs, which accounted for an unexpectedly high fraction (61 %) of all transcripts. Most of the non-coding TSSs had not been previously described (2792) and represented novel sequences (59 % of all TSSs). A large fraction of the potential novel non-coding transcripts were cis-antisense RNAs (1244 aTSSs). A strong negative correlation between the levels of antisense transcripts and cognate sense mRNAs was found, which suggested that the negative regulation of gene expression via antisense RNAs may play an important role in haloarchaea. The other types of novel non-coding transcripts corresponded to internal transcripts overlapping with mRNAs (1153 iTSSs) and intergenic small RNA (sRNA) candidates (395 TSSs).

CONCLUSION

This study provides a comprehensive map of the primary transcriptome of Hfx. volcanii grown under optimal conditions. Fewer than half of all protein-coding genes have been transcribed under these conditions. Unexpectedly, more than half of the detected TSSs belonged to several classes of non-coding RNAs. Thus, RNA-based regulation appears to play a more important role in haloarchaea than previously anticipated.

Quantitative proteomic and microarray analysis of the archaeon Methanosarcina acetivorans grown with acetate versus methanol

Methanosarcina acetivorans strain C2A is an acetate- and methanol-utilizing methane-producing organism for which the genome, the largest yet sequenced among the Archaea, reveals extensive physiological diversity. LC linear ion trap-FTICR mass spectrometry was employed to analyze acetate- vs methanol-grown cells metabolically labeled with 14N vs 15N, respectively, to obtain quantitative protein abundance ratios. DNA microarray analyses of acetate- vs methanol-grown cells was also performed to determine gene expression ratios. The combined approaches were highly complementary, extending the physiological understanding of growth and methanogenesis. Of the 1081 proteins detected, 255 were > or =3-fold differentially abundant. DNA microarray analysis revealed 410 genes that were > or =2.5-fold differentially expressed of 1972 genes with detected expression. The ratios of differentially abundant proteins were in good agreement with expression ratios of the encoding genes. Taken together, the results suggest several novel roles for electron transport components specific to acetate-grown cells, including two flavodoxins each specific for growth on acetate or methanol. Protein abundance ratios indicated that duplicate CO dehydrogenase/acetyl-CoA complexes function in the conversion of acetate to methane. Surprisingly, the protein abundance and gene expression ratios indicated a general stress response in acetate- vs methanol-grown cells that included enzymes specific for polyphosphate accumulation and oxidative stress. The microarray analysis identified transcripts of several genes encoding regulatory proteins with identity to the PhoU, MarR, GlnK, and TetR families commonly found in the Bacteria domain. An analysis of neighboring genes suggested roles in controlling phosphate metabolism (PhoU), ammonia assimilation (GlnK), and molybdopterin cofactor biosynthesis (TetR). Finally, the proteomic and microarray results suggested roles for two-component regulatory systems specific for each growth substrate.

An endangered oasis of aquatic microbial biodiversity in the Chihuahuan desert

The Cuatro Cienegas basin in the Chihuahuan desert is a system of springs, streams, and pools. These ecosystems support >70 endemic species and abundant living stromatolites and other microbial communities, representing a desert oasis of high biodiversity. Here, we combine data from molecular microbiology and geology to document the microbial biodiversity of this unique environment. Ten water samples from locations within the Cuatro Cienegas basin and two neighboring valleys as well as three samples of wet sediments were analyzed. The phylogeny of prokaryotic populations in the samples was determined by characterizing cultured organisms and by PCR amplification and sequencing of 16S rRNA genes from total community DNA. The composition of microbial communities was also assessed by determining profiles of terminal restriction site polymorphisms of 16S rRNA genes in total community DNA. There were 250 different phylotypes among the 350 cultivated strains. Ninety-eight partial 16S rRNA gene sequences were obtained and classified. The clones represented 38 unique phylotypes from ten major lineages of Bacteria and one of Archaea. Unexpectedly, 50% of the phylotypes were most closely related to marine taxa, even though these environments have not been in contact with the ocean for tens of millions of years. Furthermore, terminal restriction site polymorphism profiles and geological data suggest that the aquatic ecosystems of Cuatro Cienegas are hydrologically interconnected with adjacent valleys recently targeted for agricultural intensification. The findings underscore the conservation value of desert aquatic ecosystems and the urgent need for study and preservation of freshwater microbial communities.

Heterotrophic Archaea dominate sedimentary subsurface ecosystems off Peru

Studies of deeply buried, sedimentary microbial communities and associated biogeochemical processes during Ocean Drilling Program Leg 201 showed elevated prokaryotic cell numbers in sediment layers where methane is consumed anaerobically at the expense of sulfate. Here, we show that extractable archaeal rRNA, selecting only for active community members in these ecosystems, is dominated by sequences of uncultivated Archaea affiliated with the Marine Benthic Group B and the Miscellaneous Crenarchaeotal Group, whereas known methanotrophic Archaea are not detectable. Carbon flow reconstructions based on stable isotopic compositions of whole archaeal cells, intact archaeal membrane lipids, and other sedimentary carbon pools indicate that these Archaea assimilate sedimentary organic compounds other than methane even though methanotrophy accounts for a major fraction of carbon cycled in these ecosystems. Oxidation of methane by members of Marine Benthic Group B and the Miscellaneous Crenarchaeotal Group without assimilation of methane-carbon provides a plausible explanation. Maintenance energies of these subsurface communities appear to be orders of magnitude lower than minimum values known from laboratory observations, and ecosystem-level carbon budgets suggest that community turnover times are on the order of 100-2,000 years. Our study provides clues about the metabolic functionality of two cosmopolitan groups of uncultured Archaea.

Highly divergent genes for methanopterin-linked C1 transfer reactions in Lake Washington, assessed via metagenomic analysis and mRNA detection

The origins and the evolutionary history of tetrahydromethanopterin-linked C1 transfer reactions that are part of two environmentally important biotransformations, methylotrophy and methanogenesis, are still not well understood. In previous studies, we have expanded the known phylogenetic diversity of these reactions by identifying genes highly diverging from the ones associated with cultivated Proteobacteria, Planctomycetes, or Archaea (M. G. Kalyuzhnaya, M. E. Lidstrom, and L. Chistoserdova, Microb. Ecol. 48:463-472, 2004; M. G. Kalyuzhnaya, O. Nercessian, M. E. Lidstrom, and L. Chistoserdova, Environ. Microbiol. 7:1269-1274, 2005). Here we used a metagenomic approach to demonstrate that these divergent genes are present with high abundance in the microbial community inhabiting Lake Washington sediment. We also gained preliminary insights into the genomic composition of the organisms possessing these genes by sequencing genomic fragments from three uncultured microbes possessing the genes of interest. Phylogenetic analyses suggested that, although distantly related to each other, these organisms deeply diverge from known Bacteria and Archaea, with more relation to the former, suggesting their affiliation with a new bacterial phylum. We also demonstrate, via specific mRNA detection, that these divergent genes are expressed in the environment, pointing toward their potential role in local carbon cycling.

Diversity of microorganisms within rock varnish in the Whipple Mountains, California

Rock varnish from Arizona's Whipple Mountains harbors a microbial community containing about 10(8) microorganisms g(-1) of varnish. Analyses of varnish phospholipid fatty acids and rRNA gene libraries reveal a community comprised of mostly Proteobacteria but also including Actinobacteria, eukaryota, and a few members of the Archaea. Rock varnish represents a significant niche for microbial colonization.

Metabolically active Crenarchaeota in Altamira Cave

Altamira Cave contains valuable paleolithic paintings dating back to 15,000 years. The conservation of these unique paintings is attracting increasing interest, and so, understanding microbial proliferation in Altamira Cave represents a prioritary objective. Here, we show for the first time that members of the Crenarchaeota were metabolically active components of developing microbial communities. RNA was extracted directly from the studied environment, and a number of 16S rRNA gene sequences belonging to the low-temperature Crenarchaeota were detected. Although low-temperature Crenarchaeota detected in a variety of ecosystems by using molecular techniques remain uncultured, this RNA-based study confirms an active participation of the Crenarchaeota in cave biogeochemical cycles.

Characterization of microbial community structure in Gulf of Mexico gas hydrates: comparative analysis of DNA- and RNA-derived clone libraries

The characterization of microbial assemblages within solid gas hydrate, especially those that may be physiologically active under in situ hydrate conditions, is essential to gain a better understanding of the effects and contributions of microbial activities in Gulf of Mexico (GoM) hydrate ecosystems. In this study, the composition of the Bacteria and Archaea communities was determined by 16S rRNA phylogenetic analyses of clone libraries derived from RNA and DNA extracted from sediment-entrained hydrate (SEH) and interior hydrate (IH). The hydrate was recovered from an exposed mound located in the northern GoM continental slope with a hydrate chipper designed for use on the manned-submersible Johnson Sea Link (water depth, 550 m). Previous geochemical analyses indicated that there was increased metabolic activity in the SEH compared to the IH layer (B. N. Orcutt, A. Boetius, S. K. Lugo, I. R. Macdonald, V. A. Samarkin, and S. Joye, Chem. Geol. 205:239-251). Phylogenetic analysis of RNA- and DNA-derived clones indicated that there was greater diversity in the SEH libraries than in the IH libraries. A majority of the clones obtained from the metabolically active fraction of the microbial community were most closely related to putative sulfate-reducing bacteria and anaerobic methane-oxidizing archaea. Several novel bacterial and archaeal phylotypes for which there were no previously identified closely related cultured isolates were detected in the RNA- and DNA-derived clone libraries. This study was the first phylogenetic analysis of the metabolically active fraction of the microbial community extant in the distinct SEH and IH layers of GoM gas hydrate.

The expanding world of small RNAs in the hyperthermophilic archaeon Sulfolobus solfataricus

Archaeal L7Ae is a multifunctional protein that binds to a distinctive K-turn motif in RNA and is found as a component in the large subunit of the ribosome, and in ribose methylation and pseudouridylation guide RNP particles. A collection of L7Ae-associated small RNAs were isolated from Sulfolobus solfataricus cell extracts and used to construct a cDNA library; 45 distinct cDNA sequences were characterized and divided into six groups. Group 1 contained six RNAs that exhibited the features characteristic of the canonical C/D box archaeal sRNAs, two RNAs that were atypical C/D box sRNAs and one RNA representative of archaeal H/ACA sRNA family. Group 2 contained 13 sense strand RNA sequences that were encoded either within, or overlapping annotated open reading frames (ORFs). Group 3 contained three sequences form intergenic regions. Group 4 contained antisense sequences from within or overlapping sense strand ORFs or antisense sequences to C/D box sRNAs. More than two-thirds of these sequences possessed K-turn motifs. Group 5 contained two sequences corresponding to internal regions of 7S RNA. Group 6 consisted of 11 sequences that were fragments from the 5' or 3' ends of 16S and 23S ribosomal RNA and from seven different tRNAs. Our data suggest that S. solfataricus contains a plethora of small RNAs. Most of these are bound directly by the L7Ae protein; the others may well be part of larger, transiently stable RNP complexes that contain the L7Ae protein as core component.

Poly(A) polymerase modification and reverse transcriptase PCR amplification of environmental RNA

We describe a combination of two established techniques for a novel application for constructing full-length cDNA clone libraries from environmental RNA. The cDNA was cloned without the use of prescribed primers that target specific genes, and the procedure did not involve random priming. Purified RNA was first modified by addition of a poly(A) tail and then was amplified by using a commercially available reverse transcriptase PCR (RT-PCR) cDNA synthesis kit. To demonstrate the feasibility of this approach, a cDNA clone library was constructed from size-fractionated RNA (targeting 16S rRNA) purified from a geothermally heated soil in Yellowstone National Park in Wyoming. The resulting cDNA library contained clones representing Bacteria and Eukarya taxa and several mRNAs. There was no exact clone match between this library and a separate cDNA library generated from an RT-PCR performed with unmodified rRNA and Bacteria-specific forward and universal reverse primers that were designed from cultivated organisms; however, both libraries contained representatives of the Firmicutes and the alpha-Proteobacteria. Unexpectedly, there were no Archaea clones in the library generated from poly(A)-modified RNA. Additional RT-PCRs performed with universal and Archaea-biased primers and unmodified RNA demonstrated the presence of novel Archaea in the soil. Experiments with pure cultures of Sulfolobus solfataricus and Halobacterium halobium revealed that some Archaea rRNA may not be a suitable substrate for the poly(A) tail modification step. The protocol described here demonstrates the feasibility of directly accessing prokaryote RNA (rRNA and/or mRNA) in environmental samples, but the results also illustrate potentially important problems.

Nanoarchaeum equitans creates functional tRNAs from separate genes for their 5'- and 3'-halves

Analysis of the genome sequence of the small hyperthermophilic archaeal parasite Nanoarchaeum equitans has not revealed genes encoding the glutamate, histidine, tryptophan and initiator methionine transfer RNA species. Here we develop a computational approach to genome analysis that searches for widely separated genes encoding tRNA halves that, on the basis of structural prediction, could form intact tRNA molecules. A search of the N. equitans genome reveals nine genes that encode tRNA halves; together they account for the missing tRNA genes. The tRNA sequences are split after the anticodon-adjacent position 37, the normal location of tRNA introns. The terminal sequences can be accommodated in an intervening sequence that includes a 12-14-nucleotide GC-rich RNA duplex between the end of the 5' tRNA half and the beginning of the 3' tRNA half. Reverse transcriptase polymerase chain reaction and aminoacylation experiments of N. equitans tRNA demonstrated maturation to full-size tRNA and acceptor activity of the tRNA(His) and tRNA(Glu) species predicted in silico. As the joining mechanism possibly involves tRNA trans-splicing, the presence of an intron might have been required for early tRNA synthesis.

Stable-isotope probing of microorganisms thriving at thermodynamic limits: syntrophic propionate oxidation in flooded soil

Propionate is an important intermediate of the degradation of organic matter in many anoxic environments. In methanogenic environments, due to thermodynamic constraints, the oxidation of propionate requires syntrophic cooperation of propionate-fermenting proton-reducing bacteria and H(2)-consuming methanogens. We have identified here microorganisms that were active in syntrophic propionate oxidation in anoxic paddy soil by rRNA-based stable-isotope probing (SIP). After 7 weeks of incubation with [(13)C]propionate (<10 mM) and the oxidation of approximately 30 micromol of (13)C-labeled substrate per g dry weight of soil, we found that archaeal nucleic acids were (13)C labeled to a larger extent than those of the bacterial partners. Nevertheless, both terminal restriction fragment length polymorphism and cloning analyses revealed Syntrophobacter spp., Smithella spp., and the novel Pelotomaculum spp. to predominate in "heavy" (13)C-labeled bacterial rRNA, clearly showing that these were active in situ in syntrophic propionate oxidation. Among the Archaea, mostly Methanobacterium and Methanosarcina spp. and also members of the yet-uncultured "rice cluster I" lineage had incorporated substantial amounts of (13)C label, suggesting that these methanogens were directly involved in syntrophic associations and/or thriving on the [(13)C]acetate released by the syntrophs. With this first application of SIP in an anoxic soil environment, we were able to clearly demonstrate that even guilds of microorganisms growing under thermodynamic constraints, as well as phylogenetically diverse syntrophic associations, can be identified by using SIP. This approach holds great promise for determining the structure and function relationships of further syntrophic or other nutritional associations in natural environments and for defining metabolic functions of yet-uncultivated microorganisms.

Differential response of archaeal and bacterial communities to nitrogen inputs and pH changes in upland pasture rhizosphere soil

Grassland management regimens influence the structure of archaeal communities in upland pasture soils, which appear to be dominated by as yet uncultivated non-thermophilic Crenarchaeota. In an attempt to determine which grassland management factors select for particular crenarchaeal community structures, soil microcosm experiments were performed examining the effect of increased pH, application of inorganic fertilizer (ammonium nitrate) and sheep urine deposition on both archaeal and bacterial communities in unmanaged grassland soil. As grassland management typically increases pH, a further experiment examined the effect of a reduction in pH, to that typical of unimproved grassland soils, on archaeal and bacterial communities. The RT-PCR amplification of 16S rRNA followed by denaturing gradient gel electrophoresis analysis demonstrated a distinct and reproducible effect on bacterial communities after incubation for 28 or 30 days. In contrast, none of the treatments had a significant effect on the structure of the crenarchaeal community, indicating that these factors are not major drivers of crenarchaeal community structures in grassland soils.

The 492-bp RM07 DNA fragment from the halophilic Archaea confers promoter activity in all three domains of life

A 492-bp DNA fragment, designated RM07, was isolated from the chromosomal DNA of the halophilic Archaea, Halobacterium halobium, and was shown to confer promoter activity in Escherichia coli. Sequence analysis revealed that RM07 contained three consensus sequences of the archaeal distal promoter element as well as the typical -35 and -10 box sequences of bacterial promoters. Promoter probe analysis confirmed that RM07 conferred promoter activity in all three domains of life: Archaea (Haloferax volcanii), Eukarya (Saccharomyces cerevisiae) and Bacteria (Escherichia coli). Deletion analysis and site-directed mutagenesis further identified the functional regions within RM07 required for promoter activity. This is the first report of a DNA fragment from Archaea that confers promoter activity in all three domains of life, suggesting that the promoter structure and activity may be viewed as a bridge narrowing the gaps among the different domains of life.

Enhanced sensitivity of DNA- and rRNA-based stable isotope probing by fractionation and quantitative analysis of isopycnic centrifugation gradients

Stable isotope probing (SIP) of nucleic acids allows the detection and identification of active members of natural microbial populations that are involved in the assimilation of an isotopically labelled compound into nucleic acids. SIP is based on the separation of isotopically labelled DNA or rRNA by isopycnic density gradient centrifugation. We have developed a highly sensitive protocol for the detection of 'light' and 'heavy' nucleic acids in fractions of centrifugation gradients. It involves the fluorometric quantification of total DNA or rRNA, and the quantification of either 16S rRNA genes or 16S rRNA in gradient fractions by real-time PCR with domain-specific primers. Using this approach, we found that fully 13C-labelled DNA or rRNA of Methylobacterium extorquens was quantitatively resolved from unlabelled DNA or rRNA of Methanosarcina barkeri by cesium chloride or cesium trifluoroacetate density gradient centrifugation respectively. However, a constant low background of unspecific nucleic acids was detected in all DNA or rRNA gradient fractions, which is important for the interpretation of environmental SIP results. Consequently, quantitative analysis of gradient fractions provides a higher precision and finer resolution for retrieval of isotopically enriched nucleic acids than possible using ethidium bromide or gradient fractionation combined with fingerprinting analyses. This is a prerequisite for the fine-scale tracing of microbial populations metabolizing 13C-labelled compounds in natural ecosystems.

In vivo analysis of an essential archaeal signal recognition particle in its native host

The evolutionarily conserved signal recognition particle (SRP) plays an integral role in Sec-mediated cotranslational protein translocation and membrane protein insertion, as it has been shown to target nascent secretory and membrane proteins to the bacterial and eukaryotic translocation pores. However, little is known about its function in archaea, since characterization of the SRP in this domain of life has thus far been limited to in vitro reconstitution studies of heterologously expressed archaeal SRP components identified by sequence comparisons. In the present study, the genes encoding the SRP54, SRP19, and 7S RNA homologs (hv54h, hv19h, and hv7Sh, respectively) of the genetically and biochemically tractable archaeon Haloferax volcanii were cloned, providing the tools to analyze the SRP in its native host. As part of this analysis, an hv54h knockout strain was created. In vivo characterization of this strain revealed that the archaeal SRP is required for viability, suggesting that cotranslational protein translocation is an essential process in archaea. Furthermore, a method for the purification of this SRP employing nickel chromatography was developed in H. volcanii, allowing the successful copurification of (i) Hv7Sh with a histidine-tagged Hv54h, as well as (ii) Hv54h and Hv7Sh with a histidine-tagged Hv19h. These results provide the first in vivo evidence that these components interact in archaea. Such copurification studies will provide insight into the significance of the similarities and differences of the protein-targeting systems of the three domains of life, thereby increasing knowledge about the recognition of translocated proteins in general.

Widespread distribution of a 7S RNA in archaebacteria

An RNA of nonribosomal origin was found in the extreme halophilic bacteria. This novel small RNA was found to be a homogeneous species by RNA fingerprinting. Analysis of the ribonuclease T1 oligonucleotides gave no evidence of the presence of posttranscriptional modifications. Comparisons of electrophoretic mobility with other RNAs of known size suggest that this is a 7S RNA containing 325-375 nucleotides. An RNA of similar mobility was found in all major divisions of the archaebacteria. Insufficient sequence information is available to determine whether these RNAs are homologs of any other known small RNA.

Complexes with truncated RNAs from the large domain of Archaeoglobus fulgidus signal recognition particle

Protein SRP19 is an important component of the signal recognition particle (SRP) as it promotes assembly of protein SRP54 with SRP RNA and recognizes a tetranucleotide loop. Structural features and RNA binding activities of SRP19 of the hyperthermophilic archaeon Archaeoglobus fulgidus were investigated. An updated alignment of SRP19 sequences predicted three conserved regions and two alpha-helices. With Af-SRP RNA the Af-SRP54 protein assembled into an A. fulgidus SRP which remained intact for many hours. Stable complexes were formed between Af-SRP19 and truncated SRP RNAs, including a 36-residue fragment representing helix 6 of A. fulgidus SRP RNA.

Isolation of a minD-like gene in the hyperthermophilic archaeon Pyrococcus AL585, and phylogenetic characterization of related proteins in the three domains of life

The region upstream of the dinF-like gene of the hyperthermophilic archaeon Pyrococcus strain AL585 has been cloned and sequenced. This region contains an open reading frame (ORF) that encodes a polypeptide with a high similarity to MinD proteins and their Mrp paralogues. Transcripts of the dinF-like and the minD-like genes were detected by RT-PCR, indicating that they are both expressed in vivo. The MinD and MinD-like proteins belong to a broad family of ATPases involved in chromosome and plasmid partitioning. MinD-like proteins can be defined by specific amino-acid sequence signatures. A systematic search for proteins sharing these signatures in current databases and newly sequenced genomes show that MinD-like proteins are present in all archaeal genomes sequenced so far, often in several copies. Phylogenetic analysis identifies two groups of MinD-like proteins which are also characterized by more conserved amino-acid motifs. A first group, which includes the Escherichia coli MinD and the Pyrococcus AL585 MinDL protein, contains only procaryotic proteins. This group can be further divided into a subgroup of archaeal proteins and two subgroups of bacterial proteins. A second group includes proteins more related to the E. coli Mrp protein and contains representants of the three domains of life. The conservation of MinD-like proteins in the three domains of life suggests that these proteins play a central role in cellular metabolism.