Nidovirus subgenomic mRNAs contain a leader sequence derived from the 5′ end of the genome fused to different sequences (‘bodies’) derived from the 3′ end. Their generation involves a unique mechanism of discontinuous subgenomic RNA synthesis that resembles copy-choice RNA recombination. During this process, the nascent RNA strand is transferred from one site in the template to another, during either plus or minus strand synthesis, to yield subgenomic RNA molecules. Central to this process are transcription-regulating sequences (TRSs), which are present at both template sites and ensure the fidelity of strand transfer. Here we present results of a comprehensive co-variation mutagenesis study of equine arteritis virus TRSs, demonstrating that discontinuous RNA synthesis depends not only on base pairing between sense leader TRS and antisense body TRS, but also on the primary sequence of the body TRS. While the leader TRS merely plays a targeting role for strand transfer, the body TRS fulfils multiple functions. The sequences of mRNA leader–body junctions of TRS mutants strongly suggested that the discontinuous step occurs during minus strand synthesis.

Background Many RNA viruses do not have a single, representative genome but instead form a set of related variants that has been called a quasispecies. The sequence variability of such viruses presents a significant bioinformatics challenge. In order for the sequence information to be understood, the complete mutational spectrum needs to be distilled to a biologically relevant and analyzable representation. Results Here, we develop a "selection mapping" algorithm--QUASI--that identifies the positively selected variants of viral proteins. The key to the selection mapping algorithm is the identification of particular replacement mutations that are overabundant relative to silent mutations at each codon (e.g., threonine at hemagglutinin position 262). Selection mapping identifies such replacement mutations as positively selected. Conversely, selection mapping recognizes negatively selected variants as mutational "noise" (e.g., serine at hemagglutinin position 262). Conclusion Selection mapping is a fundamental improvement over earlier methods (e.g., dN/dS) that identify positive selection at codons but do not identify which amino acids at these codons confer selective advantage. Using QUASI's selection maps, we characterize the selected mutational landscapes of influenza A H3 hemagglutinin, HIV-1 reverse transcriptase, and HIV-1 gp120.

Background Post-transcriptional gene silencing (PTGS) by short interfering RNA has opened up new directions in the phenotypic mutation of cellular genes. However, its efficacy on non-nuclear genes and its effect on the interferon pathway remain unexplored. Since directed mutation of RNA genomes is not possible through conventional mutagenesis, we have tested sequence-specific 21-nucleotide long double-stranded RNAs (dsRNAs) for their ability to silence cytoplasmic RNA genomes. Results Short dsRNAs were generated against specific mRNAs of respiratory syncytial virus, a nonsegmented negative-stranded RNA virus with a cytoplasmic life cycle. At nanomolar concentrations, the dsRNAs specifically abrogated expression of the corresponding viral proteins, and produced the expected mutant phenotype ex vivo. The dsRNAs did not induce an interferon response, and did not inhibit cellular gene expression. The ablation of the viral proteins correlated with the loss of the specific mRNAs. In contrast, viral genomic and antigenomic RNA, which are encapsidated, were not directly affected. Conclusions Synthetic inhibitory dsRNAs are effective in specific silencing of RNA genomes that are exclusively cytoplasmic and transcribed by RNA-dependent RNA polymerases. RNA-directed RNA gene silencing does not require cloning, expression, and mutagenesis of viral cDNA, and thus, will allow the generation of phenotypic null mutants of specific RNA viral genes under normal infection conditions and at any point in the infection cycle. This will, for the first time, permit functional genomic studies, attenuated infections, reverse genetic analysis, and studies of host-virus signaling pathways using a wild type RNA virus, unencumbered by any superinfecting virus.

The ability of transfected synthetic small interfering (si) RNAs to suppress the expression of specific transcripts has proved a useful technique to probe gene function in mammalian cells. However, high production costs limit this technology's utility for many laboratories and experimental situations. Recently, several DNA-based plasmid vectors have been developed that direct transcription of small hairpin RNAs, which are processed into functional siRNAs by cellular enzymes. Although these vectors provide certain advantages over chemically synthesized siRNAs, numerous disadvantages remain including merely transient siRNA expression and low and variable transfection efficiency.

Background The genome of invertebrates is rich in retroelements which are structurally reminiscent of the retroviruses of vertebrates. Those containing three open reading frames (ORFs), including an env-like gene, may well be considered as endogenous retroviruses. Further support to this similarity has been provided by the ability of the env-like gene of DmeGypV (the Gypsy endogenous retrovirus of Drosophila melanogaster) to promote infection of Drosophila cells by a pseudotyped vertebrate retrovirus vector. Results To gain insights into their evolutionary story, a sample of thirteen insect endogenous retroviruses, which represents the largest sample analysed until now, was studied by computer-assisted comparison of the translated products of their gag, pol and env genes, as well as their LTR structural features. We found that the three phylogenetic trees based respectively on Gag, Pol and Env common motifs are congruent, which suggest a monophyletic origin for these elements. Conclusions We showed that most of the insect endogenous retroviruses belong to a major clade group which can be further divided into two main subgroups which also differ by the sequence of their primer binding sites (PBS). We propose to name IERV-K and IERV-S these two major subgroups of Insect Endogenous Retro Viruses (or Insect ERrantiVirus, according to the ICTV nomenclature) which respectively use Lys and Ser tRNAs to prime reverse transcription.

Overlapping but oppositely oriented transcripts have the potential to form sense-antisense perfect double-stranded (ds) RNA duplexes. A bioinformatics approach has identified over 217 candidate overlapping transcriptional units, bringing the total number of predicted and validated examples of overlapping but oppositely oriented transcripts to over 300.

Background Noncoding RNA genes produce transcripts that exert their function without ever producing proteins. Noncoding RNA gene sequences do not have strong statistical signals, unlike protein coding genes. A reliable general purpose computational genefinder for noncoding RNA genes has been elusive. Results We describe a comparative sequence analysis algorithm for detecting novel structural RNA genes. The key idea is to test the pattern of substitutions observed in a pairwise alignment of two homologous sequences. A conserved coding region tends to show a pattern of synonymous substitutions, whereas a conserved structural RNA tends to show a pattern of compensatory mutations consistent with some base-paired secondary structure. We formalize this intuition using three probabilistic "pair-grammars": a pair stochastic context free grammar modeling alignments constrained by structural RNA evolution, a pair hidden Markov model modeling alignments constrained by coding sequence evolution, and a pair hidden Markov model modeling a null hypothesis of position-independent evolution. Given an input pairwise sequence alignment (e.g. from a BLASTN comparison of two related genomes) we classify the alignment into the coding, RNA, or null class according to the posterior probability of each class. Conclusions We have implemented this approach as a program, QRNA, which we consider to be a prototype structural noncoding RNA genefinder. Tests suggest that this approach detects noncoding RNA genes with a fair degree of reliability.

Background Ribosomal proteins are encoded in all genomes of cellular life forms and are, generally, well conserved during evolution. In prokaryotes, the genes for most ribosomal proteins are clustered in several highly conserved operons, which ensures efficient co-regulation of their expression. Duplications of ribosomal-protein genes are infrequent, and given their coordinated expression and functioning, it is generally assumed that ribosomal-protein genes are unlikely to undergo horizontal transfer. However, with the accumulation of numerous complete genome sequences of prokaryotes, several paralogous pairs of ribosomal protein genes have been identified. Here we analyze all such cases and attempt to reconstruct the evolutionary history of these ribosomal proteins. Results Complete bacterial genomes were searched for duplications of ribosomal proteins. Ribosomal proteins L36, L33, L31, S14 are each duplicated in several bacterial genomes and ribosomal proteins L11, L28, L7/L12, S1, S15, S18 are so far duplicated in only one genome each. Sequence analysis of the four ribosomal proteins, for which paralogs were detected in several genomes, two of the ribosomal proteins duplicated in one genome (L28 and S18), and the ribosomal protein L32 showed that each of them comes in two distinct versions. One form contains a predicted metal-binding Zn-ribbon that consists of four conserved cysteines (in some cases replaced by histidines), whereas, in the second form, these metal-chelating residues are completely or partially replaced. Typically, genomes containing paralogous genes for these ribosomal proteins encode both versions, designated C+ and C-, respectively. Analysis of phylogenetic trees for these seven ribosomal proteins, combined with comparison of genomic contexts for the respective genes, indicates that in most, if not all cases, their evolution involved a duplication of the ancestral C+ form early in bacterial evolution, with subsequent alternative loss of the C+ and C- forms in different lineages. Additionally, evidence was obtained for a role of horizontal gene transfer in the evolution of these ribosomal proteins, with multiple cases of gene displacement 'in situ', that is, without a change of the gene order in the recipient genome. Conclusions A more complex picture of evolution of bacterial ribosomal proteins than previously suspected is emerging from these results, with major contributions of lineage-specific gene loss and horizontal gene transfer. The recurrent theme of emergence and disruption of Zn-ribbons in bacterial ribosomal proteins awaits a functional interpretation.

Background Alternative DNA conformations are of particular interest as potential signals to mark important sites on the genome. The structural variability of CA microsatellites is particularly pronounced; these are repetitive poly(CA) · poly(TG) DNA sequences spread in all eukaryotic genomes as tracts of up to 60 base pairs long. Many in vitro studies have shown that the structure of poly(CA) · poly(TG) can vary markedly from the classical right handed DNA double helix and adopt diverse alternative conformations. Here we have studied the mechanism of formation and the structure of an alternative DNA structure, named Form X, which was observed previously by polyacrylamide gel electrophoresis of DNA fragments containing a tract of the CA microsatellite poly(CA) · poly(TG) but had not yet been characterized. Results Formation of Form X was found to occur upon reassociation of the strands of a DNA fragment containing a tract of poly(CA) · poly(TG), in a process strongly stimulated by the nuclear proteins HMG1 and HMG2. By inserting Form X into DNA minicircles, we show that the DNA strands do not run fully side by side but instead form a DNA knot. When present in a closed DNA molecule, Form X becomes resistant to heating to 100°C and to alkaline pH. Conclusions Our data strongly support a model of Form X consisting in a DNA loop at the base of which the two DNA duplexes cross, with one of the strands of one duplex passing between the strands of the other duplex, and reciprocally, to form a semicatenated DNA junction also called a DNA hemicatenane.

Recombinant adenoviruses are straightforward to produce at high titres, have a promiscuous host-range, and, because of their ability to infect nondividing cells, lend themselves to in vivo gene delivery. Such advantages have led to their widespread and successful use in preclinical studies of arthritis gene therapy. While adenoviral vectors are well suited to 'proof of principle' experiments in laboratory animals, there are several barriers to their use in human studies at this time. Transient transgene expression limits their application to strategies, such as synovial ablation, which do not require extended periods of gene expression. Moreover, there are strong immunological barriers to repeat dosing. In addition, safety concerns predicate local, rather than systemic, delivery of the virus. Continued engineering of the adenoviral genome is producing vectors with improved properties, which may eventually overcome these issues. Promising avenues include the development of 'gutted' vectors encoding no endogenous viral genes and of adenovirus–AAV chimeras. Whether these will offer advantages over existing vectors, which may already provide safe, long-term gene expression following in vivo delivery, remains to be seen.