Evidence for bacterial origin of heat shock RNA-1
- 1Department of Life Science, College of Natural Science, Chung-Ang University, Seoul 156-756, South Korea
- 2Research Center for Biomolecules and Biosystems, Chung-Ang University, Seoul 156-756, South Korea
- 3Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 305-701, South Korea
Abstract
The heat shock RNA-1 (HSR1) is a noncoding RNA (ncRNA) reported to be involved in mammalian heat shock response. HSR1 was shown to significantly stimulate the heat-shock factor 1 (HSF1) trimerization and DNA binding. The hamster HSR1 sequence was reported to consist of 604 nucleotides (nt) plus a poly(A) tail and to have only a 4-nt difference with the human HSR1. In this study, we present highly convincing evidence for bacterial origin of the HSR1. No HSR1 sequence was found by exhaustive sequence similarity searches of the publicly available eukaryotic nucleotide sequence databases at the NCBI, including the expressed sequence tags, genome survey sequences, and high-throughput genomic sequences divisions of GenBank, as well as the Trace Archive database of whole genome shotgun sequences, and genome assemblies. Instead, a putative open reading frame (ORF) of HSR1 revealed strong similarity to the amino-terminal region of bacterial chloride channel proteins. Furthermore, the 5′ flanking region of the putative HSR1 ORF showed similarity to the 5′ upstream regions of the bacterial protein genes. We propose that the HSR1 was derived from a bacterial genome fragment either by horizontal gene transfer or by bacterial infection of the cells. The most probable source organism of the HSR1 is a species belonging to the order Burkholderiales.
Keywords
INTRODUCTION
A large number of mRNA-like long noncoding RNAs (ncRNAs) have been reported to be involved in crucial biological processes in mammalian cells including transcriptional regulation and epigenetic gene regulation (Prasanth and Spector 2007; St. Laurent and Wahlestedt 2007; Ponting et al. 2009; Wilusz et al. 2009). For example, ncRNA XIST is associated exclusively with the inactive X chromosome (Brown et al. 1992) and ncRNA H19 is expressed only from the maternal allele of the imprinted IGF2/H19 locus (Brannan et al. 1990; Webber et al. 1998). Many long ncRNAs have been reported to be directly involved in regulation of activities of the associated proteins. The ncRNA Evf-2, for example, forms a stable complex with the Dlx-2 protein to increase the transcriptional activity of the Dlx-5/6 enhancer (Feng et al. 2006). The maternally expressed imprinted ncRNA MEG3 activates p53 and functions as a tumor suppressor (Zhou et al. 2007). Recently it has been shown that large intergenic ncRNAs such as XIST and HOTAIR guide chromatin-modifying complexes to specific genomic loci and act as epigenetic regulators of gene expression (Khalil et al. 2009).
Many of ncRNAs show clear evolutionary conservation, indicating that they are subject to strong purifying selection (Guttman et al. 2009; Ponjavic et al. 2009). Some of these highly conserved ncRNA genes such as HAR1F show significant evolutionary acceleration in the human genome, possibly resulting in emergence of human-specific traits (Pollard et al. 2006). Interestingly, ncRNA XIST evolved from a protein-coding gene and a set of transposable elements in placental mammals after the divergence of placentals and marsupials (Duret et al. 2006; Elisaphenko et al. 2008).
It has been reported that a novel mammalian ncRNA called heat shock RNA-1 (HSR1) activates the heat shock transcription factor 1 (HSF1), which is essential for the induction of expression of heat shock proteins (HSPs) (Shamovsky et al. 2006). The HSR1 was initially isolated from the hamster kidney cell line BHK-21. The hamster HSR1 was detected as an RNA species of ∼2 kb, and it is composed of 604 nucleotides (nt) and a poly(A) tail. The human HSR1 was reported to be almost identical to the hamster HSR1 with a 4-nt difference (Shamovsky et al. 2006). Subsequently, it has been reported by the same group of researchers that noncoding RNAs homologous to mammalian HSR1 are present in other eukaryotic organisms including Xenopus, Drosophila, and Caenorhabditis elegans, and that the homologs are functionally interchangeable (abstract number 2C_03_S at the World Conference of Stress in 2007; available at http://www.stress07.com/051abs2.htm). However, the HSR1 homolog sequences of the latter three species are not publicly available.
The importance of HSR1 was immediately recognized by the ncRNA research field (Kugel and Goodrich 2006; Costa 2007). Involvement of HSR1 in the trimerization of HSF1 was referred as favorable support for existence of noncoding-RNA regulators of other important processes such as RNA polymerase II transcription (Goodrich and Kugel 2006).
In this study, we present that the HSR1 sequence shows strong sequence similarity to the 5′ flanking region and part of the coding region of chloride channel proteins of bacterial genomes. We also discuss possible mechanism for detection of the HSR1 molecules in the mammalian cells.
RESULTS AND DISCUSSION
As an initial attempt to perform functional analysis of the human HSR1 gene, we set to find its sequence information from the publicly available nucleotide sequence databases. Initially, we searched the National Center for Biotechnology Information (NCBI) nonredundant nucleotide sequence database with the BLASTN program using the hamster HSR1 sequence as a query. We could not obtain any match from human or other mammalian genomic clones or mRNA sequences. Instead, we obtained many partial-length hits to bacterial genomes with HSR1 regions 140–210 and 250–310 with the E-values ranging from 10−9 to 10−4. Next, we searched genome assemblies of mammalian species including human, chimpanzee, orangutan, rhesus macaque, marmoset, mouse, rat, guinea pig, cat, dog, horse, cow, opossum, and platypus with the BLAT program available at the University of California Santa Cruz (UCSC) Genome Browser (http://genome.ucsc.edu/). No match was found in these genome assemblies either.
Then we speculated that the HSR1 gene could be embedded in a heterochromatic region which has been known to be refractory to cloning and sequencing (International Human Genome Sequencing Consortium 2004). Heterochromatic regions are reported to contain many transcriptionally active genes (Lyle et al. 2007). Some constitutively heterochromatic regions become transcriptionally active in response to stress stimuli such as heat shock (Rizzi et al. 2004). The HSR1 RNA was reported to be constitutively expressed in the hamster BHK-21 cells and the human HeLa cells (Shamovsky et al. 2006). The HSR1 RNA molecules were polyadenylated and detected in size of ∼2 kb. Therefore, we assumed that cDNA sequences for HSR1 RNA could have been isolated as expressed sequence tags (ESTs) and searched the database for ESTs (dbEST). We also searched other high-throughput genome sequence databases such as genomic survey sequences (GSS), high-throughput genomic sequences (HTGS), and whole-genome shotgun reads (WGS) at the NCBI BLAST website. To find matches to unassembled, raw sequence data for ESTs and genomes, furthermore, we searched the NCBI Trace Archive databases using the Discontiguous MegaBLAST program. The searched trace data included human, chimpanzee, bonobo, rhesus macaque, marmoset, mouse, rat, cow, dog, cat, and horse. However, we were not able to find any significant nucleotide matches from all the high-throughput sequence data, suggesting that no sequence record for the HSR1 has been deposited in the current sequence databases.
Interestingly, however, a BLASTX search of the NCBI nonredundant protein database yielded about 100 meaningful hits to bacterial proteins. Inspection of titles of these proteins revealed that almost all the matches were chloride channel proteins; the rest were labeled as hypothetical proteins. The matches were between a putative open reading frame (ORF) of the 3′ part of the HSR1 sequence (nucleotide positions 390–604) and the amino-terminal regions of the bacterial chloride channel proteins. The predicted HSR1 ORF would encode a total of 71 amino acid residues and runs off at the end of the sequence which was reported to be polyadenylated in the mammalian cells.
The matched bacterial protein sequences were downloaded and sorted according to the genus name. The protein sequences were from 33 bacterial genera (Table 1). From each genus, only one sequence was selected based on FASTA alignment score between the bacterial protein and the HSR1 ORF sequence. A multiple sequence alignment of the HSR1 ORF and the selected 33 bacterial chloride channel proteins with MUSCLE software revealed many conserved residues (Fig. 1). The HSR1 ORF shares 33 (with Burkholderia, Lutiella, and Pseudomonas proteins) to 18 identical amino acid residues with the bacterial proteins out of its 71 residues. We also found the similar positioning of the translation start codon and the conserved amino acid sequence motifs, which clearly indicates that the HSR1 ORF would encode a truncated version of the chloride channel protein. It is unlikely that the HSR1 sequence shows this high degree of sequence similarity to the bacterial proteins simply by chance. Rather, it suggests that the HSR1 is actually derived from a bacterial protein coding gene.
List of organisms and proteins studied
Multiple sequence alignment of predicted HSR ORF and bacterial chloride channel proteins. The protein sequence deduced from the HSR1 ORF (top) and bacterial channel proteins with sequence similarity from 33 bacterial genera were aligned multiply by using the MUSCLE program. A nucleotide was removed to produce the HSR1 ORF (see Fig. 2). Detailed organism names and NCBI accession numbers for the proteins can be found in Table 1. AA, number of amino acid residues aligned; ID, number of identical residues to the HSR1 ORF.
A phylogenetic tree was constructed to infer the possible phylogenetic position of the HSR1 ORF with amino acid sequence alignments. Bootstrap consensus trees were constructed using maximum parsimony (MP) and neighbor-joining (NJ) methods. The result was phylogenetically uninformative. Almost all the branches were not supported by bootstrap values, and the position of HSR1 ORF was ambiguous. This is because the ORF sequences are mutationally saturated and highly diverged except for the functionally important amino acid positions, or because the size of the sequence data is not sufficient for reliable results.
To examine whether the sequence similarity can be found at the nucleotide level, we retrieved genome sequences of the bacteria showing high protein sequence identity to the HSR1 ORF. When the HSR1 nucleotide sequence and the genomic sequences of species of Burkholderia, Lutiella, Ralstonia, and Delftia genera were multiply aligned, the 3′ half of the HSR1 sequence (from 240 to 604) showed a high level of identity with these bacterial genomes (Fig. 2): 57.0% with the Burkholderia species, and 56.7% with the Lutiella and the Ralstonia species, and 55.3% with the Delftia species. The highest pairwise identity among these five sequences was 69.0% obtained between the Burkholderia and the Ralstonia species. The nucleotide sequence conservation extends to the 5′ upstream region (from 240 to 389) of the HSR1 ORF. In fact, sequence conservation is higher in the 5′ upstream region than in the ORF itself, indicating that the region may harbor transcriptional and/or translational regulatory elements. Interestingly, the HSR1 regions 140–210 and 250–310 yielded meaningful hits to many bacterial genomes with significant E-values (up to 10−9) by BLASTN search. Moreover, the region 140–210 showed multiple matches in each genome. Many of them were located within or adjacent to an integrase/transposase of an insertion sequence (IS) element, implying that the 5′ half of the HSR1 may be part of a bacterial IS element. The high nucleotide sequence identity between the HSR1 sequence and the bacterial genomes further supports that the HSR1 is derived from a genomic fragment of currently uncharacterized bacterium. All the four bacterial species above belong to Betaproteobacteria: Burkholderia, Ralstonia, and Delftia species are in the order Burkholderiales and Lutiella in the order Neisseriales.
Multiple alignment of bacterial genome sequences and a part of HSR1. A part of HSR1 sequence from 240 to 604 was aligned with four representative bacterial genomes. The start codons for the chloride channel protein ORFs are indicated by three consecutive asterisks. An additional nucleotide (arrowhead) in the HSR1 ORF, which would be present due to an insertion or perhaps a sequencing error, was removed to deduce the protein sequence. Residues identical in four or five sequences are highlighted on black background. GenBank accession numbers and positions for the nucleotide sequences are: Burkholderia, NZ_ABBH01000052.1:1764-2131; Lutiella, NZ_ACIS01000003.1:194423-194791; Ralstonia, NZ_AAKL01000043.1:37423-37060; Delftia, NC_010002.1: 1926518-1926158.
It is questionable how the bacterial genome sequence-like HSR1 molecule was discovered in the hamster and human cell lines and other eukaryotic cells. One possible mechanism is the horizontal gene transfer. As demonstrated in Wolbachia-insect cases, integration of bacterial genome fragments to eukaryotic nuclear genome is rather common (Nikoh et al. 2008). Some of the integrated Wolbachia genes were even detected by RT-PCR, indicating that they were transcriptionally active in the host insect genome. Similarly, the HSR1 gene might have been horizontally transferred from a bacterium that had infected an early metazoan ancestor. The integrated bacterial genomic fragment could have been exapted as a regulator of the heat shock process. As in the case of XIST, which evolved from a protein-coding gene (Duret et al. 2006; Elisaphenko et al. 2008), the HSR1 sequence might have lost its protein-coding potential.
Another possible mechanism is a simple bacterial infection of the cells examined. Infection of eukaryotic cell lines with microbial organisms such as Mycoplasma species is rather common and problematic during cell culture research (Schmitt and Pawlita 2009). The gene sequences derived from the pathogenic viruses, bacteria, fungi, and protozoa can be detected in the human cells and tissues (Weber et al. 2002). If this is the case, the HSR1 RNA molecules would have been directly derived from some bacterial cells or transcribed from the nuclear genome after the integration. The absence of HSR1 sequence in mammalian genomes and transcripts despite their extensive sequence data deposited in the current sequence databases may justify the simple bacterial infection. Furthermore, the fact that the most similar sequence to HSR1 was found in the intracellular pathogens Ralstonia and Burkholderia (Valvano et al. 2005) could be suggestive of this possibility.
In conclusion, we have found that the mammalian HSR1 sequence shows strong similarity to bacterial genomes, implying that the HSR1 molecules have bacterial origin. Judging from all the data, we propose that the most probable source organism of the HSR1 is a species belonging to the order Burkholderiales. We advise that the origin of the HSR1 should be further scrutinized for better understanding of the evolution and the mechanism of mammalian heat shock process.
MATERIALS AND METHODS
Database search and analysis of the HSR1 sequence
BLAST programs were used to search the NCBI sequence databases (Altschul et al. 1997). All the NCBI nucleotide databases including the nonredundant database, EST database (dbEST), genomic survey sequences (GSS), high-throughput genomic sequences (HTGS), whole-genome shotgun (WGS) reads, and environmental samples (ENV_NT) were searched with the BLASTN program (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Discontiguous MegaBLAST searches of the WGS trace databases were carried out at the NCBI Trace Archive (http://www.ncbi.nlm.nih.gov/Traces/). BLASTX search of the NCBI nonredundant protein database was performed at the NCBI website with default parameters. All the sequence similarity searches were conducted on September 24, 2009.
BLAT searches of the mammalian genome assemblies such as human (hg18), chimpanzee (panTro2), orangutan (ponAbe2), rhesus macaque (rheMac2), marmoset (calJac1), mouse (mm9), rat (rn4), guinea pig (cavPor3), cat (felCat3), dog (canFam2), horse (equCab2), cow (bosTau4), opossum (monDom5), and platypus (ornAna1) were performed at the UCSC Genome Browser server (http://genome.ucsc.edu/) (Kent 2002; Kent et al. 2002).
The FASTA program was used to obtain pairwise alignment (Pearson and Lipman 1988). Multiple sequence alignment was generated using MUSCLE (Edgar 2004) and visualized using ClustalX2 (Larkin et al. 2007) or BOXSHADE (http://www.ch.embnet.org/software/BOX_form.html). Sequence Logo representation (Schneider and Stephens 1990) was generated with the WebLogo server (Crooks et al. 2004).
Molecular phylogenetic analysis
Molecular phylogenetic analysis was performed by using MEGA program version 4 (http://www.megasoftware.net/) (Kumar et al. 2008). We used MP and NJ methods to infer phylogenetic trees. Positions with gaps were deleted from the alignment. Bootstrap consensus trees were constructed with 1000 replicates.
We also used PAUP* software (version 4.0b10 for Linux platform) (http://paup.scs.fsu.edu/) (Swofford 2002). We constructed bootstrap 50% majority-rule consensus trees using MP and NJ methods with 100 replicates. The gaps were treated as missing data.
ACKNOWLEDGMENTS
This work was supported in part by the Korea Research Foundation Grant funded by the Korean Government (MOEHRD, Basic Research Promotion Fund) (KRF-2008-331-C00263), and in part by Grant (R01-2008-000-11660-0) from the Korea Science and Engineering Foundation.
Footnotes
-
Reprint requests to: Yoonsoo Hahn, Department of Life Science, College of Natural Science, Chung-Ang University, 47 Heukseok-ro, Dongjak-gu, Seoul 156-756, South Korea; e-mail: hahny{at}cau.ac.kr; fax: 82-2-826-5812.
-
Article published online ahead of print. Article and publication date are at http://www.rnajournal.org/cgi/doi/10.1261/rna.1879610.
-
- Received August 13, 2009.
- Accepted October 27, 2009.
- Copyright © 2010 RNA Society
