An expanding universe of noncoding RNAs

GenRate: a generative model that reveals novel transcripts in genome-tiling microarray data

Genome-wide microarray designs containing millions to hundreds of millions of probes are available for a variety of mammals, including mouse and human. These genome tiling arrays can potentially lead to significant advances in science and medicine, e.g., by indicating new genes and alternative primary and secondary transcripts. While bottom-up pattern matching techniques (e.g., hierarchical clustering) can be used to find gene structures in microarray data, we believe the many interacting hidden variables and complex noise patterns more naturally lead to an analysis based on generative models. We describe a generative model of tiling data and show how the sum-product algorithm can be used to infer hybridization noise, probe sensitivity, new transcripts, and alternative transcripts. The method, called GenRate, maximizes a global scoring function that enables multiple transcripts to compete for ownership of putative probes. We apply GenRate to a new exon tiling dataset from mouse chromosome 4 and show that it makes significantly more predictions than a previously described hierarchical clustering method at the same false positive rate. GenRate correctly predicts many known genes and also predicts new gene structures. As new problems arise, additional hidden variables can be incorporated into the model in a principled fashion, so we believe that GenRate will prove to be a useful tool in the new era of genome-wide tiling microarray analysis.

RNA thermometers

Temperature is an important parameter that free-living cells monitor constantly. The expression of heat-shock, cold-shock and some virulence genes is coordinated in response to temperature changes. Apart from protein-mediated transcriptional control mechanisms, translational control by RNA thermometers is a widely used regulatory strategy. RNA thermometers are complex RNA structures that change their conformation in response to temperature. Most, but not all, RNA thermometers are located in the 5'-untranslated region and mask ribosome-binding sites by base pairing at low temperatures. Melting of the structure at increasing temperature permits ribosome access and translation initiation. Different cis-acting RNA thermometers and a trans-acting thermometer will be presented.

Molecular characterization of a novel translocation t(5;14)(q21;q32) in a patient with congenital abnormalities

Chromosomal translocations are frequently found to be associated with various malignant disorders as well as congenital abnormalities. We report the characterization of a novel reciprocal translocation t(5;14)(q21;q32) in a patient with congenital abnormalities manifested by severe mental retardation, athetotic tetraplegia, microcephaly, peculiar facies (upward slanting of palpebral fissures), clinodactyly of the fifth fingers, and overlapping toes. Using a JHGP24 lymphoblast cell line derived from this patient, metaphase fluorescence in situ hybridization with bacterial artificial chromosome and cosmid probes and subsequent molecular analysis mapped the translocation breakpoint to the nucleotide level. Sequence analysis of the breakpoint junctions revealed the presence of a homologous sequence, GTGGC, along with a single nucleotide substitution and an insertion in der(14), and a single nucleotide deletion in the der(5) chromosome. We also attempted to identify and characterize the transcripts near the breakpoint by 5' and 3' rapid amplification of cDNA ends. Although we found several transcripts near the breakpoint of chromosome 14, the lack of significant ORFs within these transcripts suggests they are likely to be non-coding RNAs. These transcripts may have an important role in the neurogenesis or differentiation.

Experimental approaches to identify non-coding RNAs

Cellular RNAs that do not function as messenger RNAs (mRNAs), transfer RNAs (tRNAs) or ribosomal RNAs (rRNAs) comprise a diverse class of molecules that are commonly referred to as non-protein-coding RNAs (ncRNAs). These molecules have been known for quite a while, but their importance was not fully appreciated until recent genome-wide searches discovered thousands of these molecules and their genes in a variety of model organisms. Some of these screens were based on biocomputational prediction of ncRNA candidates within entire genomes of model organisms. Alternatively, direct biochemical isolation of expressed ncRNAs from cells, tissues or entire organisms has been shown to be a powerful approach to identify ncRNAs both at the level of individual molecules and at a global scale. In this review, we will survey several such wet-lab strategies, i.e. direct sequencing of ncRNAs, shotgun cloning of small-sized ncRNAs (cDNA libraries), microarray analysis and genomic SELEX to identify novel ncRNAs, and discuss the advantages and limits of these approaches.

Features of Arabidopsis genes and genome discovered using full-length cDNAs

Arabidopsis is currently the reference genome for higher plants. A new, more detailed statistical analysis of Arabidopsis gene structure is presented including intron and exon lengths, intergenic distances, features of promoters, and variant 5'-ends of mRNAs transcribed from the same transcription unit. We also provide a statistical characterization of Arabidopsis transcripts in terms of their size, UTR lengths, 3'-end cleavage sites, splicing variants, and coding potential. These analyses were facilitated by scrutiny of our collection of sequenced full-length cDNAs and much larger collection of 5'-ESTs, together with another set of full-length cDNAs from Salk/Stanford/Plant Gene Expression Center/RIKEN. Examples of alternative splicing are observed for transcripts from 7% of the genes and many of these genes display multiple spliced isoforms. Most splicing variants lie in non-coding regions of the transcripts. Non-canonical splice sites constitute less than 1% of all splice sites. Genes with fewer than four introns display reduced average mRNA levels. Putative alternative transcription start sites were observed in 30% of highly expressed genes and in more than 50% of the genes with low expression. Transcription start sites correlate remarkably well with a CG skew peak in the DNA sequences. The intergenic distances vary considerably, those where genes are transcribed towards one another being significantly shorter. New transcripts, missing in the current TIGR genome annotation and ESTs that are non-coding, including those antisense to known genes, are derived and cataloged in the Supplementary Material. They identify 148 new loci in the Arabidopsis genome. The conclusions drawn provide a better understanding of the Arabidopsis genome and how the gene transcripts are processed. The results also allow better predictions to be made for, as yet, poorly defined genes and provide a reference for comparisons with other plant genomes whose complete sequences are currently being determined. Some comparisons with rice are included in this paper.

Chronic lymphocytic leukemia: molecular genetics and animal models

Chronic lymphocytic leukemia accounts for almost 30% of all adult leukemia cases in the United States and Western Europe. Although several common genomic abnormalities in CLL have been identified, mutational and functional analysis of corresponding genes so far have not proved their involvement in CLL. Our latest studies demonstrated functional involvement of Tcl1 oncoprotein and microRNA genes in the pathogenesis of CLL. Deregulated expression of Tcl1 in transgenic mice resulted in CLL. These CLL tumors showed abnormalities in expression of murine microRNA genes mmu-mir-15a and mmu-mir-16-1. Interestingly, human homologs of these genes, mir-15a and mir-16-1, located at the chromosome 13q14 are also deleted in human CLL samples. In this review we summarize and discuss these new developments. These recently emerged insights into the molecular mechanisms of CLL will allow for the development of new approaches to treat this disease.

Three-dimensional structures of fibrillar Sm proteins: Hfq and other Sm-like proteins

Hfq is a nucleic acid-binding protein that functions as a global regulator of gene expression by virtue of its interactions with several small, non-coding RNA species. Originally identified as an Escherichia coli host factor required for RNA phage Qbeta replication, Hfq is now known to post-transcriptionally regulate bacterial gene expression by modulating both mRNA stability and translational activity. Recently shown to be a member of the diverse Sm protein family, Hfq adopts the OB-like fold typical of other Sm and Sm-like (Lsm) proteins, and also assembles into toroidal homo-oligomers that bind single-stranded RNA. Similarities between the structures, functions, and evolution of Sm/Lsm proteins and Hfq are continually being discovered, and we now report an additional, unexpected biophysical property that is shared by Hfq and other Sm proteins: E.coli Hfq polymerizes into well-ordered fibres whose morphologies closely resemble those found for Sm-like archaeal proteins (SmAPs). However, the hierarchical assembly of these fibres is dissimilar: whereas SmAPs polymerize into polar tubes (and striated bundles of such tubes) by head-to-tail stacking of individual homo-heptamers, helical Hfq fibres are formed by cylindrical slab-like layers that consist of 36 subunits arranged as a hexamer of Hfq homo-hexamers (i.e. protofilaments in a 6 x 6 arrangement). The different fibrillar ultrastructures formed by Hfq and SmAP are presented and examined herein, with the overall goal of elucidating another similarity amongst the diverse members of the Sm protein family.

Sequence-based heuristics for faster annotation of non-coding RNA families

MOTIVATION

Non-coding RNAs (ncRNAs) are functional RNA molecules that do not code for proteins. Covariance Models (CMs) are a useful statistical tool to find new members of an ncRNA gene family in a large genome database, using both sequence and, importantly, RNA secondary structure information. Unfortunately, CM searches are extremely slow. Previously, we created rigorous filters, which provably sacrifice none of a CM's accuracy, while making searches significantly faster for virtually all ncRNA families. However, these rigorous filters make searches slower than heuristics could be.

RESULTS

In this paper we introduce profile HMM-based heuristic filters. We show that their accuracy is usually superior to heuristics based on BLAST. Moreover, we compared our heuristics with those used in tRNAscan-SE, whose heuristics incorporate a significant amount of work specific to tRNAs, where our heuristics are generic to any ncRNA. Performance was roughly comparable, so we expect that our heuristics provide a high-quality solution that--unlike family-specific solutions--can scale to hundreds of ncRNA families.

AVAILABILITY

The source code is available under GNU Public License at the supplementary web site.

Predicting candidate genomic sequences that correspond to synthetic functional RNA motifs

Riboswitches and RNA interference are important emerging mechanisms found in many organisms to control gene expression. To enhance our understanding of such RNA roles, finding small regulatory motifs in genomes presents a challenge on a wide scale. Many simple functional RNA motifs have been found by in vitro selection experiments, which produce synthetic target-binding aptamers as well as catalytic RNAs, including the hammerhead ribozyme. Motivated by the prediction of Piganeau and Schroeder [(2003) Chem. Biol., 10, 103–104] that synthetic RNAs may have natural counterparts, we develop and apply an efficient computational protocol for identifying aptamer-like motifs in genomes. We define motifs from the sequence and structural information of synthetic aptamers, search for sequences in genomes that will produce motif matches, and then evaluate the structural stability and statistical significance of the potential hits. Our application to aptamers for streptomycin, chloramphenicol, neomycin B and ATP identifies 37 candidate sequences (in coding and non-coding regions) that fold to the target aptamer structures in bacterial and archaeal genomes. Further energetic screening reveals that several candidates exhibit energetic properties and sequence conservation patterns that are characteristic of functional motifs. Besides providing candidates for experimental testing, our computational protocol offers an avenue for expanding natural RNA's functional repertoire.

Evidence for a preferential targeting of 3'-UTRs by cis-encoded natural antisense transcripts

Although both the 5'- and 3'-untranslated regions (5'- and 3'-UTRs) of eukaryotic mRNAs may play a crucial role in posttranscriptional gene regulation, we observe that cis-encoded natural antisense RNAs have a striking preferential complementarity to the 3'-UTRs of their target genes in mammalian (human and mouse) genomes. A null neutral model, evoking differences in the rate of 3'-UTR and 5'-UTR extension, could potentially explain high rates of 3'-to-3' overlap compared with 5'-to-5' overlap. However, employing a simulation model we show that this null model probably cannot explain the finding that 3'-to-3' overlapping pairs have a much higher probability (>5 times) of conservation in both mouse and human genomes with the same overlapping pattern than do 5'-to-5' overlaps. Furthermore, it certainly cannot explain the finding that overlapping pairs seen in both genomes have a significantly higher probability of having co-expression and inverse expression (i.e. characteristic of sense-antisense regulation) than do overlapping pairs seen in only one of the two species. We infer that the function of many 3'-to-3' overlaps is indeed antisense regulation. These findings underscore the preference for, and conservation of, 3'-UTR-targeted antisense regulation, and the importance of 3'-UTRs in gene regulation.

Uncertainty principle of genetic information in a living cell

Background

Formal description of a cell's genetic information should provide the number of DNA molecules in that cell and their complete nucleotide sequences. We pose the formal problem: can the genome sequence forming the genotype of a given living cell be known with absolute certainty so that the cell's behaviour (phenotype) can be correlated to that genetic information? To answer this question, we propose a series of thought experiments.

Results

We show that the genome sequence of any actual living cell cannot physically be known with absolute certainty, independently of the method used. There is an associated uncertainty, in terms of base pairs, equal to or greater than μs (where μ is the mutation rate of the cell type and s is the cell's genome size).

Conclusion

This finding establishes an "uncertainty principle" in genetics for the first time, and its analogy with the Heisenberg uncertainty principle in physics is discussed. The genetic information that makes living cells work is thus better represented by a probabilistic model rather than as a completely defined object.

RNase P: role of distinct protein cofactors in tRNA substrate recognition and RNA-based catalysis

The Escherichia coli ribonuclease P (RNase P) has a protein component, termed C5, which acts as a cofactor for the catalytic M1 RNA subunit that processes the 5′ leader sequence of precursor tRNA. Rpp29, a conserved protein subunit of human RNase P, can substitute for C5 protein in reconstitution assays of M1 RNA activity. To better understand the role of the former protein, we compare the mode of action of Rpp29 to that of the C5 protein in activation of M1 RNA. Enzyme kinetic analyses reveal that complexes of M1 RNA–Rpp29 and M1 RNA–C5 exhibit comparable binding affinities to precursor tRNA but different catalytic efficiencies. High concentrations of substrate impede the activity of the former complex. Rpp29 itself exhibits high affinity in substrate binding, which seems to reduce the catalytic efficiency of the reconstituted ribonucleoprotein. Rpp29 has a conserved C-terminal domain with an Sm-like fold that mediates interaction with M1 RNA and precursor tRNA and can activate M1 RNA. The results suggest that distinct protein folds in two unrelated protein cofactors can facilitate transition from RNA- to ribonucleoprotein-based catalysis by RNase P.

sRNAPredict: an integrative computational approach to identify sRNAs in bacterial genomes

Small non-coding bacterial RNAs (sRNAs) play important regulatory roles in a variety of cellular processes. Nearly all known sRNAs have been identified in Escherichia coli and most of these are not conserved in the majority of other bacterial species. Many of the E.coli sRNAs were initially predicted through bioinformatic approaches based on their common features, namely that they are encoded between annotated open reading frames and are flanked by predictable transcription signals. Because promoter consensus sequences are undetermined for most species, the successful use of bioinformatics to identify sRNAs in bacteria other than E.coli has been limited. We have created a program, sRNAPredict, which uses coordinate-based algorithms to integrate the respective positions of individual predictive features of sRNAs and rapidly identify putative intergenic sRNAs. Relying only on sequence conservation and predicted Rho-independent terminators, sRNAPredict was used to search for sRNAs in Vibrio cholerae. This search identified 9 of the 10 known or putative V.cholerae sRNAs and 32 candidates for novel sRNAs. Small transcripts for 6 out of 9 candidate sRNAs were observed by Northern analysis. Our findings suggest that sRNAPredict can be used to efficiently identify novel sRNAs even in bacteria for which promoter consensus sequences are not available.

Silence of the transcripts: RNA interference in medicine

Silencing of gene expression by ribonucleic acid (RNA), known as RNA interference (RNAi), is now recognized as a major means of gene regulation in biology. In this mechanism, small noncoding double-stranded RNA molecules knock down gene expression through a variety of mechanisms that include messenger RNA (mRNA) degradation, inhibition of mRNA translation, or chromatin remodeling. The posttranscriptional mechanism of RNAi has been embraced by researchers as a powerful tool for generating deficient phenotypes without mutating the gene. In parallel, exciting recent results have promised its application in disease therapy. This review aims to summarize the current knowledge in this area and provide a roadmap that may eventually launch RNAi from the research bench to the medicine chest.

The role of small RNAs in human diseases: potential troublemaker and therapeutic tools

Small RNAs, including short interfering RNAs (siRNAs) and microRNAs (miRNAs), are ubiquitous, versatile repressors of gene expression in plants, animals, and many fungi. They can trigger destruction of homologous mRNA or inhibition of cognate mRNA translation and play an important role in maintaining the stable state of chromosome structure and regulating the expression of protein-coding genes. Furthermore, the recent research showed that there exists close relationship between small RNAs and human diseases. Several human diseases have surfaced in which miRNAs or their machinery might be implicated, such as some neurological diseases and cancers. The specificity and potency of small RNAs suggest that they might be promising as therapeutic agents. This article will review the role of small RNAs in some human diseases etiology and investigations of taking siRNAs as therapeutic tools for treating viral infection, cancer, and other diseases. We also discuss the potential of miRNAs in gene therapy.

Meiotic messenger RNA and noncoding RNA targets of the RNA-binding protein Translin (TSN) in mouse testis

In postmeiotic male germ cells, TSN, formerly known as testis brain-RNA binding protein, is found in the cytoplasm and functions as a posttranscriptional regulator of a group of genes transcribed by the transcription factor CREM-tau. In contrast, in pachytene spermatocytes, TSN is found predominantly in nuclei. Tsn-null males show a reduced sperm count and high levels of apoptosis in meiotic cells, suggesting a critical function for TSN during meiosis. To identify meiotic target RNAs that associate in vivo with TSN, we reversibly cross-linked TSN to RNA in testis extracts from 17-day-old and adult mice and immunoprecipitated the complexes with an affinity-purified TSN antibody. Extracts from Tsn-null mice were used as controls. Cloning and sequencing the immunoprecipitated RNAs, we identified four new TSN target mRNAs, encoding diazepam-binding inhibitor-like 5, arylsulfatase A, a tetratricopeptide repeat structure-containing protein, and ring finger protein 139. In contrast to the population of postmeiotic translationally delayed mRNAs that bind TSN, these four mRNAs are initially expressed in pachytene spermatocytes. In addition, anti-TSN also precipitated a nonprotein-coding RNA (ncRNA), which is abundant in nuclei of pachytene spermatocytes and has a putative polyadenylation signal, but no open reading frame. A second similar ncRNA is adjacent to a GGA repeat, a motif frequently associated with recombination hot spots. RNA gel-shift assays confirm that the four new target mRNAs and the ncRNA specifically bind to TSN in testis extracts. These studies have, for the first time, identified both mRNAs and a ncRNA as TSN targets expressed during meiosis.

Evidence for abundant transcription of non-coding regions in the Saccharomyces cerevisiae genome

Background

Recent studies in a growing number of organisms have yielded accumulating evidence that a significant portion of the non-coding region in the genome is transcribed. We address this issue in the yeast Saccharomyces cerevisiae.

Results

Taking into account the absence of a significantly large yeast EST database, we use microarray expression data collected for genomic regions erroneously believed to be coding to study the expression pattern of non-coding regions in the Saccharomyces cerevisiae genome. We find that at least 164 out of 589 (28%) such regions are expressed under specific biological conditions. In particular, looking at the probes that are located opposing other known genes at the same genomic locus, we find that 88 out of 341 (26%) of these genes support antisense transcription. The expression patterns of these antisense genes are positively correlated. We validate these results using RT-PCR on a sample of 6 non-coding transcripts.

Conclusion

1. The yeast genome is transcribed on a scale larger than previously assumed. 2. Correlated transcription of antisense genes is abundant in the yeast genome. 3. Antisense genes in yeast are non-coding.

Bacterial transcription elongation factors: new insights into molecular mechanism of action

Like transcription initiation, the elongation and termination stages of transcription cycle serve as important targets for regulatory factors in prokaryotic cells. In this review, we discuss the recent progress in structural and biochemical studies of three evolutionarily conserved elongation factors, GreA, NusA and Mfd. These factors affect RNA polymerase (RNAP) processivity by modulating transcription pausing, arrest, termination or anti-termination. With structural information now available for RNAP and models of ternary elongation complexes, the interaction between these factors and RNAP can be modelled, and possible molecular mechanisms of their action can be inferred. The models suggest that these factors interact with RNAP at or near its three major, nucleic acid-binding channels: Mfd near the upstream opening of the primary (DNA-binding) channel, NusA in the vicinity of both the primary channel and the RNA exit channel, and GreA within the secondary (backtracked RNA-binding) channel, and support the view that these channels are involved in the maintenance of RNAP processivity.

Predicting non-coding RNA genes in Escherichia coli with boosted genetic programming

Several methods exist for predicting non-coding RNA (ncRNA) genes in Escherichia coli (E.coli). In addition to about sixty known ncRNA genes excluding tRNAs and rRNAs, various methods have predicted more than thousand ncRNA genes, but only 95 of these candidates were confirmed by more than one study. Here, we introduce a new method that uses automatic discovery of sequence patterns to predict ncRNA genes. The method predicts 135 novel candidates. In addition, the method predicts 152 genes that overlap with predictions in the literature. We test sixteen predictions experimentally, and show that twelve of these are actual ncRNA transcripts. Six of the twelve verified candidates were novel predictions. The relatively high confirmation rate indicates that many of the untested novel predictions are also ncRNAs, and we therefore speculate that E.coli contains more ncRNA genes than previously estimated.

Deconstructing PTI-1: PTI-1 is a truncated, but not mutated, form of translation elongatin factor 1A1, eEF1A1

The prostate tumor-inducing gene 1 (PTI-1) transcript is detected in various human carcinoma cells. PTI-1 is reported to consist of a 5' untranslated region (5' UTR) homologous to mycoplasma 23S rRNA and a coding region corresponding to a truncated and mutated form of the translation elongation factor 1A, eEF1A. We have found that the PTI-1 transcript may encode a truncated, but not mutated, form of the human isoform eEF1A1. Additionally, the 5' UTR sequence of PTI-1 from genomic DNA of different cell lines and blood samples varies from the original sequence. This 5' -UTR region of PTI-1 presents a fusion of E. coli and Mycoplasma hyorhinis 23S rRNA. We have overexpressed the potential PTI-1 protein in E. coli and various human cell lines. The resulting protein could be detected by western blotting using anti-eEF1A antibodies. However, we were unable to detect the PTI-1 protein in LNCaP cell extracts. The potential roles of the PTI-1 protein in carcinogenesis and the origin of the PTI-1 gene in the human genome are discussed.

Blockade of murine erythroleukemia cell differentiation by hypomethylating agents causes accumulation of discrete small poly(A)- RNAs hybridized to 3'-end flanking sequences of beta(major) globin gene

Induction of murine erythroleukemia (MEL) cell differentiation is accompanied by transcriptional activation of globin genes and biosynthesis of hemoglobin. In this study, we observed cytoplasmic accumulation of relatively small RNAs of different size (150-600 nt) hybridized to alpha1 and beta(major) globin DNA probes in MEL cells blocked to differentiate by hypomethylating agents (neplanocin A, 3-deazaneplanocin A and cycloleucine). These RNAs lack poly(A) tail and appear to be quite stable. Search within the 3'-end flanking sequences of beta(major) globin gene revealed the presence of a B1 repeat element, several ATG initiation codons, a GATA-1 consensus sequence and sequences recognized by AP-1/NF-E2 and erythroid Krüppel-like factor (EKLF) transcription factors. These data taken together indicate that exposure of MEL cells to hypomethylating agents promotes accumulation of relatively small discrete RNA transcripts lacking poly(A) tail regardless of the presence or absence of inducer dimethylsulfoxide (DMSO). However, the relative steady-state level of small RNAs was comparatively higher in cells co-exposed to inducer and each one of the hypomethylating agents. Although the orientation of these RNAs has not been established as yet, the possibility these small poly(A)- RNAs which are induced by hypomethylating agents may be involved in the blockade of MEL cell differentiation program is discussed.

NONCODE: an integrated knowledge database of non-coding RNAs

NONCODE is an integrated knowledge database dedicated to non-coding RNAs (ncRNAs), that is to say, RNAs that function without being translated into proteins. All ncRNAs in NONCODE were filtered automatically from literature and GenBank, and were later manually curated. The distinctive features of NONCODE are as follows: (i) the ncRNAs in NONCODE include almost all the types of ncRNAs, except transfer RNAs and ribosomal RNAs. (ii) All ncRNA sequences and their related information (e.g. function, cellular role, cellular location, chromosomal information, etc.) in NONCODE have been confirmed manually by consulting relevant literature: more than 80% of the entries are based on experimental data. (iii) Based on the cellular process and function, which a given ncRNA is involved in, we introduced a novel classification system, labeled process function class, to integrate existing classification systems. (iv) In addition, some 1100 ncRNAs have been grouped into nine other classes according to whether they are specific to gender or tissue or associated with tumors and diseases, etc. (v) NONCODE provides a user-friendly interface, a visualization platform and a convenient search option, allowing efficient recovery of sequence, regulatory elements in the flanking sequences, secondary structure, related publications and other information. The first release of NONCODE (v1.0) contains 5339 non-redundant sequences from 861 organisms, including eukaryotes, eubacteria, archaebacteria, virus and viroids. Access is free for all users through a web interface at http://noncode.bioinfo.org.cn.

Seven fundamental, unsolved questions in molecular biology. Cooperative storage and bi-directional transfer of biological information by nucleic acids and proteins: an alternative to "central dogma"

The Human Genome Mapping Project provided us a large amount of sequence data. However our understanding of these data did not grow proportionally, because old dogmas still set the limits of our thinking. The gene-centric, reductionistical side of molecular biology is reviewed and seven problems are formulated, each indicating the insufficiency of the "central dogma". The following is concluded and suggested: 1. Genes are located and expressed on both DNA strands; 2. Introns are the source of important biological regulation and diversity; 3. Repeats are the frame of the chromatin structure and participate in the chromatin regulation; 4. The molecular accessibility of the canonical dsDNA structure is poor; 5. The genetic code is co-evolved with the amino acids and there is a stereochemical matching between the codes andamino acids; 6. The flow of information between nucleic acids and proteins is bi-directional and reverse translation might exist; 7. Complex genetic information is always carried and stored by nucleic acids and proteins together.

Renal carcinomas with the t(6;11)(p21;q12): clinicopathologic features and demonstration of the specific alpha-TFEB gene fusion by immunohistochemistry, RT-PCR, and DNA PCR

A highly distinctive subset of renal neoplasms of children and young adults contains a t(6;11)(p21;q12), a translocation recently been shown to result in fusion of Alpha, a gene on 11q12, with the transcription factor gene TFEB on 6p21. To define the clinicopathologic spectrum of this nascent entity and to establish immunohistochemical (IHC) and molecular methods for the detection of the specific Alpha-TFEB fusion, we studied 7 renal neoplasms that showed the t(6;11) by cytogenetic or molecular analysis (patient age: range, 9-33 years; mean, 17 years). While all tumors were confined to the kidney, 3 tumors demonstrated vascular invasion. In limited follow-up, none has metastasized. We postulated that the Alpha-TFEB gene fusion may result in deregulated expression of TFEB protein that would be detectable by IHC. Using a polyclonal antibody to TFEB on formalin-fixed, paraffin-embedded tissue sections, we found that all 7 renal neoplasms with the t(6;11) demonstrated moderate (2 cases) or strong (5 cases) nuclear TFEB immunoreactivity. In contrast, none of 1089 other tumors (of 74 histologic types from 16 sites) labeled significantly for TFEB. Nuclear immunoreactivity for TFEB in normal tissues was extremely rare, limited to weak labeling of scattered benign lymphocytes. We also show that the Alpha-TFEB fusion RNAs are highly variable in size and structure, making detection by reverse-transcriptase polymerase chain reaction (RT-PCR) less reliable than for other gene fusions. Because Alpha is an intronless gene and therefore lacks splice signals, we hypothesized that DNA PCR and RT-PCR products would be identical, allowing for the use of more robust molecular assays based on genomic DNA. Indeed, in 2 cases with available frozen tissue, we showed the genomic Alpha-TFEB junction detected by DNA PCR to be identical to the Alpha-TFEB fusion mRNA detected by RT-PCR. In summary, renal neoplasms with the t(6;11) are a distinctive neoplastic entity with many similarities to the Xp11 translocation carcinomas, and together with the latter form a growing "MiTF/TFE family" of translocation carcinomas. Nuclear immunoreactivity for TFEB protein is a highly sensitive and specific diagnostic marker for these renal neoplasms. Finally, the special molecular features of the Alpha-TFEB gene fusion allow its molecular detection by DNA PCR as a robust alternative to RT-PCR in clinical tumor samples.

Modular RNA architecture revealed by computational analysis of existing pseudoknots and ribosomal RNAs

Modular architecture is a hallmark of RNA structures, implying structural, and possibly functional, similarity among existing RNAs. To systematically delineate the existence of smaller topologies within larger structures, we develop and apply an efficient RNA secondary structure comparison algorithm using a newly developed two-dimensional RNA graphical representation. Our survey of similarity among 14 pseudoknots and subtopologies within ribosomal RNAs (rRNAs) uncovers eight pairs of structurally related pseudoknots with non-random sequence matches and reveals modular units in rRNAs. Significantly, three structurally related pseudoknot pairs have functional similarities not previously known: one pair involves the 3′ end of brome mosaic virus genomic RNA (PKB134) and the alternative hammerhead ribozyme pseudoknot (PKB173), both of which are replicase templates for viral RNA replication; the second pair involves structural elements for translation initiation and ribosome recruitment found in the viral internal ribosome entry site (PKB223) and the V4 domain of 18S rRNA (PKB205); the third pair involves 18S rRNA (PKB205) and viral tRNA-like pseudoknot (PKB134), which probably recruits ribosomes via structural mimicry and base complementarity. Additionally, we quantify the modularity of 16S and 23S rRNAs by showing that RNA motifs can be constructed from at least 210 building blocks. Interestingly, we find that the 5S rRNA and two tree modules within 16S and 23S rRNAs have similar topologies and tertiary shapes. These modules can be applied to design novel RNA motifs via build-up-like procedures for constructing sequences and folds.

Fast and reliable prediction of noncoding RNAs

We report an efficient method for detecting functional RNAs. The approach, which combines comparative sequence analysis and structure prediction, already has yielded excellent results for a small number of aligned sequences and is suitable for large-scale genomic screens. It consists of two basic components: (i) a measure for RNA secondary structure conservation based on computing a consensus secondary structure, and (ii) a measure for thermodynamic stability, which, in the spirit of a z score, is normalized with respect to both sequence length and base composition but can be calculated without sampling from shuffled sequences. Functional RNA secondary structures can be identified in multiple sequence alignments with high sensitivity and high specificity. We demonstrate that this approach is not only much more accurate than previous methods but also significantly faster. The method is implemented in the program rnaz, which can be downloaded from www.tbi.univie.ac.at/~wash/RNAz. We screened all alignments of length n > or = 50 in the Comparative Regulatory Genomics database, which compiles conserved noncoding elements in upstream regions of orthologous genes from human, mouse, rat, Fugu, and zebrafish. We recovered all of the known noncoding RNAs and cis-acting elements with high significance and found compelling evidence for many other conserved RNA secondary structures not described so far to our knowledge.

Engineering a ligand-dependent RNA transcriptional activator

RNA has recently been shown to play diverse roles in gene regulation, including the small molecule-dependent inhibition of translation in prokaryotes. To create an artificial genetic switch that acts at the level of transcription, we fused a small molecule binding aptamer to a previously evolved RNA that activates transcription when localized to a promoter. We designed a conformational shift in which a helical element required for transcriptional activation was stabilized upon ligand binding. Selection and screening in S. cerevisiae optimized the linker region, generating an RNA that is 10-fold more active in the presence of tetramethylrosamine (TMR). TMR increases the activity of this evolved RNA in a graded, dose-dependent manner. Our results exemplify a strategy for controlling the activity of laboratory-evolved RNAs in living cells.

Cycling of the Sm-like protein Hfq on the DsrA small regulatory RNA

Small RNAs (sRNAs) regulate bacterial genes involved in environmental adaptation. This RNA regulation requires Hfq, a bacterial Sm-like protein that stabilizes sRNAs and enhances RNA-RNA interactions. To understand the mechanism of target recognition by sRNAs, we investigated the interactions between Hfq, the sRNA DsrA, and its regulatory target rpoS mRNA, which encodes the stress response sigma factor. Nuclease footprinting revealed that Hfq recognized multiple sites in rpoS mRNA without significantly perturbing secondary structure in the 5' leader that inhibits translation initiation. Base-pairing with DsrA, however, made the rpoS ribosome binding site fully accessible, as predicted by genetic data. Hfq bound DsrA four times more tightly than the DsrA.rpoS RNA complex in gel mobility-shift assays. Consequently, Hfq is displaced rapidly from its high-affinity binding site on DsrA by conformational changes in DsrA, when DsrA base-pairs with rpoS mRNA. Hfq accelerated DsrA.rpoS RNA association and stabilized the RNA complex up to twofold. Hybridization of DsrA and rpoS mRNA was optimal when Hfq occupied its primary binding site on free DsrA, but was inhibited when Hfq associated with the DsrA.rpoS RNA complex. We conclude that recognition of rpoS mRNA is stimulated by binding of Hfq to free DsrA sRNA, followed by release of Hfq from the sRNA.mRNA complex.

Profiling and Searching for RNA Pseudoknot Structures in Genomes

We developed a new method that can profile and efficiently search for pseudoknot structures in noncoding RNA genes. It profiles interleaving stems in pseudoknot structures with independent Covariance Model (CM) components. The statistical alignment score for searching is obtained by combining the alignment scores from all CM components. Our experiments show that the model can achieve excellent accuracy on both random and biological data. The efficiency achieved by the method makes it possible to search for structures that contain pseudoknot in genomes of a variety of organisms.

The small RNA regulators of Escherichia coli: roles and mechanisms*

Small noncoding RNAs have been found in all organisms, primarily as regulators of translation and message stability. The most exhaustive searches have taken place in E. coli, resulting in identification of more than 50 small RNAs, or 1%-2% of the number of protein-coding genes. One large class of these small RNAs uses the RNA chaperone Hfq; members of this class act by pairing to target messenger RNAs. Among the members of this class are DsrA and RprA, which positively regulate rpoS translation, OxyS, which negatively regulates rpoS translation and fhlA translation, RyhB, which reapportions iron use in the cell by downregulating translation of many genes that encode Fe-containing proteins, and Spot 42, which changes the polarity of translation in the gal operon. The promoters of these small RNAs are tightly regulated, frequently as part of well-understood regulons. Lessons learned from the study of small RNAs in E. coli can be applied to finding these important regulators in other organisms.

Porcine MYF6 gene: sequence, homology analysis, and variation in the promoter region

MYF6 gene codes for the bHLH transcription factor belonging to MyoD family. Its expression accompanies the processes of differentiation and maturation of myotubes during embriogenesis and continues on a relatively high level after birth, affecting the muscle phenotype. The porcine MYF6 gene was amplified and sequenced and compared with MYF6 gene sequences of other species. The amino acid sequence was deduced and an interspecies homology analysis was performed. Myf-6 protein shows a high conservation among species of 99 and 97% identity when comparing pig with cow and human, respectively, and of 93% when comparing pig with mouse and rat. The single nucleotide polymorphism (SNP) was revealed within the promoter region, which appeared to be T --> C transition recognized by a MspI restriction enzyme.

Strategies for RNA folding and assembly

RNA is structurally very flexible, which provides the basis for its functional diversity. An RNA molecule can often adopt different conformations, which enables the regulation of its function through folding. Proteins help RNAs reach their functionally active conformation by increasing their structural stability or by chaperoning the folding process. Large, dynamic RNA-protein complexes, such as the ribosome or the spliceosome, require numerous proteins that coordinate conformational switches of the RNA components during assembly and during their respective activities.

LongSAGE analysis revealed the presence of a large number of novel antisense genes in the mouse genome

MOTIVATION

Despite the increasing notions of the functional importance of antisense transcripts in gene regulation, the genome-wide overview on the ontology of antisense genes has not been obtained. Therefore, we tried to find novel antisense genes genome-wide by using our LongSAGE dataset of 202 015 tags (consisting of 41 718 unique tags), experimentally generated from mouse embryonic tail libraries.

RESULTS

We identified 1260 potential antisense genes, of which 1001 are not annotated in EnsEMBL, thereby being regarded as novel. Interestingly their sense counterparts were co-expressed in the majority of the cases.

CONCLUSIONS

The use of LongSAGE transcriptome data is extremely powerful in the identification of thus-far unknown antisense transcripts, even in the case of well-characterized organisms like the mouse.

CONTACT

imai@gsf.de.

Human cytomegalovirus 5-kilobase immediate-early RNA is a stable intron

Immediate-early viral gene products of human cytomegalovirus (HCMV) are derived from several genomic loci and largely serve to establish a cellular environment conducive to viral replication. We have further examined an unusual immediate-early transcript known as the 5-kb RNA, concluding that it is a stable intron encoded by HCMV. The 5-kb RNA is highly AT rich in sequence and lacks open reading frames likely to be translated into protein. We confirmed the absence of polyadenylation of the transcript and showed that it is primarily nuclear localized during viral infection. We mapped the 5' end of the 5-kb RNA to a consensus splice donor site and localized the 3' end in the vicinity of a splice acceptor site. In transfection studies, we showed that the 5-kb RNA can be spliced from a heterologous primary transcript. Using bacterial artificial chromosome technology, we constructed a viral recombinant containing a mutation in the 5' splice donor site that defines the 5' end of the RNA and found that this mutation eliminates expression of the 5-kb RNA during viral infection. This mutant grows in human fibroblasts without complementation. Taken together, these data support the conclusion that the 5-kb RNA is a stable intron expressed by HCMV.

Finishing the euchromatic sequence of the human genome

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers approximately 99% of the euchromatic genome and is accurate to an error rate of approximately 1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human genome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead.

Escherichia coli Hfq has distinct interaction surfaces for DsrA, rpoS and poly(A) RNAs

The bacterial Sm-like protein Hfq facilitates RNA-RNA interactions involved in post-transcriptional regulation of the stress response. Specifically, Hfq helps pair noncoding RNAs (ncRNAs) with complementary regions of target mRNAs. To probe the mechanism of this pairing, we generated a series of Hfq mutants and measured their affinity for RNAs like those with which Hfq must associate in vivo. We tested the mutants' DsrA-dependent activation of rpoS, and their ability to stabilize DsrA ncRNA against degradation in vivo. Our results suggest that Hfq has two independent RNA-binding surfaces. In addition to a well-known site around the core of the Hfq hexamer, we observe interactions with the distal face of Hfq, a new locus with which mRNAs and poly(A) sequences associate. Our model explains how Hfq can simultaneously bind a ncRNA and its mRNA target to facilitate the strand displacement reaction required for Hfq-dependent translational regulation.

Revealing the world of RNA interference

The recent discoveries of RNA interference and related RNA silencing pathways have revolutionized our understanding of gene regulation. RNA interference has been used as a research tool to control the expression of specific genes in numerous experimental organisms and has potential as a therapeutic strategy to reduce the expression of problem genes. At the heart of RNA interference lies a remarkable RNA processing mechanism that is now known to underlie many distinct biological phenomena.

Candidates for novel RNA topologies

Because the functional repertiore of RNA molecules, like proteins, is closely linked to the diversity of their shapes, uncovering RNA's structural repertoire is vital for identifying novel RNAs, especially in genomic sequences. To help expand the limited number of known RNA families, we use graphical representation and clustering analysis of RNA secondary structures to predict novel RNA topologies and their abundance as a function of size. Representing the essential topological properties of RNA secondary structures as graphs enables enumeration, generation, and prediction of novel RNA motifs. We apply a probabilistic graph-growing method to construct the RNA structure space encompassing the topologies of existing and hypothetical RNAs and cluster all RNA topologies into two groups using topological descriptors and a standard clustering algorithm. Significantly, we find that nearly all existing RNAs fall into one group, which we refer to as "RNA-like"; we consider the other group "non-RNA-like". Our method predicts many candidates for novel RNA secondary topologies, some of which are remarkably similar to existing structures; interestingly, the centroid of the RNA-like group is the tmRNA fold, a pseudoknot having both tRNA-like and mRNA-like functions. Additionally, our approach allows estimation of the relative abundance of pseudoknot and other (e.g. tree) motifs using the "edge-cut" property of RNA graphs. This analysis suggests that pseudoknots dominate the RNA structure universe, representing more than 90% when the sequence length exceeds 120 nt; the predicted trend for <100 nt agrees with data for existing RNAs. Together with our predictions for novel "RNA-like" topologies, our analysis can help direct the design of functional RNAs and identification of novel RNA folds in genomes through an efficient topology-directed search, which grows much more slowly in complexity with RNA size compared to the traditional sequence-based search.

A comparative method for finding and folding RNA secondary structures within protein-coding regions

Existing computational methods for RNA secondary-structure prediction tacitly assume RNA to only encode functional RNA structures. However, experimental studies have revealed that some RNA sequences, e.g. compact viral genomes, can simultaneously encode functional RNA structures as well as proteins, and evidence is accumulating that this phenomenon may also be found in Eukaryotes. We here present the first comparative method, called RNA-DECODER, which explicitly takes the known protein-coding context of an RNA-sequence alignment into account in order to predict evolutionarily conserved secondary-structure elements, which may span both coding and non-coding regions. RNA-DECODER employs a stochastic context-free grammar together with a set of carefully devised phylogenetic substitution-models, which can disentangle and evaluate the different kinds of overlapping evolutionary constraints which arise. We show that RNA-DECODER's parameters can be automatically trained to successfully fold known secondary structures within the HCV genome. We scan the genomes of HCV and polio virus for conserved secondary-structure elements, and analyze performance as a function of available evolutionary information. On known secondary structures, RNA-DECODER shows a sensitivity similar to the programs MFOLD, PFOLD and RNAALIFOLD. When scanning the entire genomes of HCV and polio virus for structure elements, RNA-DECODER's results indicate a markedly higher specificity than MFOLD, PFOLD and RNAALIFOLD.

A comprehensive transcript index of the human genome generated using microarrays and computational approaches

BACKGROUND

Computational and microarray-based experimental approaches were used to generate a comprehensive transcript index for the human genome. Oligonucleotide probes designed from approximately 50,000 known and predicted transcript sequences from the human genome were used to survey transcription from a diverse set of 60 tissues and cell lines using ink-jet microarrays. Further, expression activity over at least six conditions was more generally assessed using genomic tiling arrays consisting of probes tiled through a repeat-masked version of the genomic sequence making up chromosomes 20 and 22.

RESULTS

The combination of microarray data with extensive genome annotations resulted in a set of 28,456 experimentally supported transcripts. This set of high-confidence transcripts represents the first experimentally driven annotation of the human genome. In addition, the results from genomic tiling suggest that a large amount of transcription exists outside of annotated regions of the genome and serves as an example of how this activity could be measured on a genome-wide scale.

CONCLUSIONS

These data represent one of the most comprehensive assessments of transcriptional activity in the human genome and provide an atlas of human gene expression over a unique set of gene predictions. Before the annotation of the human genome is considered complete, however, the previously unannotated transcriptional activity throughout the genome must be fully characterized.

Circular box C/D RNAs in Pyrococcus furiosus

Box C/D RNAs are small, noncoding RNAs that function in RNA modification in eukaryotes and archaea. Here, we report that box C/D RNAs exist in the rare biological form of RNA circles in the hyperthermophilic archaeon Pyrococcus furiosus. Northern analysis of box C/D RNAs reveals two prominent RNA species of different electrophoretic mobilities in total P. furiosus RNA preparations. Together, the results of Northern, ribozyme, RT-PCR, and lariat debranching analyses indicate that the two species are circular and linear RNAs of similar length and abundance. It seems that most, if not all, species of box C/D RNAs exist as circles in P. furiosus. In addition, the circular RNAs are found in complexes with proteins required for box C/D RNA function. Our finding places box C/D RNAs among the extremely few circular RNAs known to exist in nature. Moreover, the unexpected discovery of circular box C/D RNAs points to the existence of a previously unrecognized biogenesis pathway for box C/D RNAs in archaea.

Mouse-centric comparative transcriptomics of protein coding and non-coding RNAs

The largest transcriptome reported so far comprises 60,770 mouse full-length cDNA clones, and is an effective reference data set for comparative transcriptomics. The number of mouse cDNAs identified greatly exceeds the number of genes predicted from the sequenced human and mouse genomes. This is largely because of extensive alternative splicing and the presence of many non-coding RNAs (ncRNAs), which are difficult to predict from genomic sequences. Notably, ncRNAs are a major component of the transcriptomes of higher organisms, and many sense-antisense pairs have been identified. The ncRNAs function in a range of regulatory mechanisms for gene expression and other biological processes. They might also have contributed to the increased functional diversification of genomes during evolution. In this review, we discuss aspects of the transcriptome of various organisms in relation to the mouse data, in order to shed light on the regulatory mechanisms and physiological significance of these abundant RNAs.

6S RNA function enhances long-term cell survival

6S RNA was identified in Escherichia coli >30 years ago, but the physiological role of this RNA has remained elusive. Here, we demonstrate that 6S RNA-deficient cells are at a disadvantage for survival in stationary phase, a time when 6S RNA regulates transcription. Growth defects were most apparent as a decrease in the competitive fitness of cells lacking 6S RNA. To decipher the molecular mechanisms underlying the growth defects, we have expanded studies of 6S RNA effects on transcription. 6S RNA inhibition of sigma(70)-dependent transcription was not ubiquitous, in spite of the fact that the vast majority of sigma(70)-RNA polymerase is bound by 6S RNA during stationary phase. The sigma(70)-dependent promoters inhibited by 6S RNA contain an extended -10 promoter element, suggesting that this feature may define a class of 6S RNA-regulated genes. We also discovered a secondary effect of 6S RNA in the activation of sigma(S)-dependent transcription at several promoters. We conclude that 6S RNA regulation of both sigma(70) and sigma(S) activities contributes to increased cell persistence during nutrient deprivation.

Isolation of eight novel Caenorhabditis elegans small RNAs

Eight novel small RNAs that were encoded in the regions corresponding to the introns of protein-coding genes were isolated from Caenorhabditis elegans. Seven of them showed a typical snoRNA secondary structure: one C/D snoRNA and six H/ACA snoRNAs. The remaining one RNA did not show any homology to other RNAs in a database. Four of the seven isolated snoRNAs could form base pairings with parts of rRNAs, suggesting that they are potential pseudouridilation sites and methylation sites. The results of our study suggest that there are more as-yet-unidentified small ncRNAs of which genes are located in the intron regions of protein-coding genes in C. elegans.

MSARI: multiple sequence alignments for statistical detection of RNA secondary structure

We present a highly accurate method for identifying genes with conserved RNA secondary structure by searching multiple sequence alignments of a large set of candidate orthologs for correlated arrangements of reverse-complementary regions. This approach is growing increasingly feasible as the genomes of ever more organisms are sequenced. A program called msari implements this method and is significantly more accurate than existing methods in the context of automatically generated alignments, making it particularly applicable to high-throughput scans. In our tests, it discerned clustalw-generated multiple sequence alignments of signal recognition particle or RNaseP orthologs from controls with 89.1% sensitivity at 97.5% specificity and with 74.4% sensitivity with no false positives in 494 controls. We used msari to conduct a comprehensive scan for secondary structure in mRNAs of coding genes, and we found many genes with known mRNA secondary structure and compelling evidence for secondary structure in other genes. msari uses a method for coping with sequence redundancy that is likely to have applications in a large set of other comparison-based search methods. The program is available for download from http://theory.csail.mit.edu/MSARi.

Inhibition of HIV-1 fusion with small interfering RNAs targeting the chemokine coreceptor CXCR4

RNA interference (RNAi) is an evolutionarily conserved process by which plants and animals protect their genomes utilizing small, double-stranded RNAs to degrade target RNAs in a sequence-specific manner. Post-transcriptional gene silencing by these moieties can lead to degradation of both cellular and viral RNAs. It has recently been shown that double-stranded, small interfering RNAs (siRNAs) of 21-25 nucleotides can be transfected into relevant cells to target specific RNAs. This approach was utilized to inhibit human immunodeficiency virus type I (HIV-1) infection in human cells. siRNAs with homology to a motif in the mRNA that encodes for the HIV-1 chemokine coreceptor CXCR4 was utilized. Complementary studies via immunofluorescence microscopy and fluorescence-activated cell sorting demonstrated downregulation of CXCR4 from the surface of cells transfected with the specific siRNAs. As well, siRNAs without sequence homology to CXCR4 were used as controls and demonstrated no downregulation of CXCR4. siRNAs targeted to another chemokine coreceptor, APJ, showed specificity for downregulation of APJ but had no effects on CXCR4. Transfections with siRNAs targeting CXCR4 mRNA were shown to inhibit HIV-1 envelope fusion, which is relatively resistant to most viral inhibitors targeting chemokine coreceptors. The specificity of this effect was demonstrated by the inhibition of fusion by CXCR4-tropic and dual-tropic (CXCR4 and CCR5) envelope glycoproteins from HIV-1 on CXCR4+ indicator cells, but the lack of effects by siRNAs targeting CXCR4 mRNA on dual-tropic HIV-1 envelopes in CCR5+ indicator cells utilizing these fusion assays. Interestingly, siRNAs targeting CXCR4 selectively inhibited CXCR4-tropic cell-free virus infection of human cells but at only modest levels as compared to cell:cell fusion. siRNA may be a potential molecular therapeutic approach to alter a cellular cofactor critical for infection of human cells by relevant strains of HIV-1. The targeting of a cellular cofactor, rather than the HIV-1-specific mRNAs or genomic RNA, holds promise as the rapid mutational ability of the HIV-1 genome may obviate the potential clinical use of RNAi directly against this virus.