Completion of the genome sequence of Lettuce necrotic yellows virus, type species of the genus Cytorhabdovirus

We completed the genome sequence of Lettuce necrotic yellows virus (LNYV) by determining the nucleotide sequences of the 4a (putative phosphoprotein), 4b, M (matrix protein), G (glycoprotein) and L (polymerase) genes. The genome consists of 12,807 nucleotides and encodes six genes in the order 3' leader-N-4a(P)-4b-M-G-L-5' trailer. Sequences were derived from clones of a cDNA library from LNYV genomic RNA and from fragments amplified using reverse transcription-polymerase chain reaction. The 4a protein has a low isoelectric point characteristic for rhabdovirus phosphoproteins. The 4b protein has significant sequence similarities with the movement proteins of capillo- and trichoviruses and may be involved in cell-to-cell movement. The putative G protein sequence contains a predicted 25 amino acids signal peptide and endopeptidase cleavage site, three predicted glycosylation sites and a putative transmembrane domain. The deduced L protein sequence shows similarities with the L proteins of other plant rhabdoviruses and contains polymerase module motifs characteristic for RNA-dependent RNA polymerases of negative-strand RNA viruses. Phylogenetic analysis of this motif among rhabdoviruses placed LNYV in a group with other sequenced cytorhabdoviruses, most closely related to Strawberry crinkle virus. (c) 2005 Elsevier B.V. All rights reserved.

Fate of hairpin transcript components during RNA silencing and its suppression in transgenic virus-resistant tobacco

Transgenic tobacco plants, carrying a Potato virus Y (PVY)-NIa hairpin sequence separated by a unique unrelated spacer sequence were specifically silenced and highly resistant to PVY infection. In such plants neither PVY-NIa nor spacer transgene transcripts were detectable by specific quantitative real time reverse transcriptase PCR (RT-qPCR) assays of similar relative efficiencies developed for direct comparative analysis. However, small interfering RNAs (siRNAs) specific for the PVY sequence of the transgene and none specific for the LNYV spacer sequence were detected. Following infection with Cucumber mosaic virus (CMV), which suppresses dsRNA-induced RNA silencing, transcript levels of PVY-NIa as well as spacer sequence increased manifold with the same time course. The cellular abundance of the single-stranded (ss) spacer sequence was consistently higher than that of PVY dsRNA in all cases. The results show that during RNA silencing and its suppression of a hairpin transcript in transgenic tobacco, the ssRNA spacer sequence is affected differently than the dsRNA. In PVY-silenced plants. the spacer is efficiently degraded by a mechanism not involving the accumulation of siRNAs, while following suppression of RNA silencing by CMV, the spacer appears protected from degradation. Crown Copyright (c) 2006 Published by Elsevier B.V. All rights reserved.

Mutations in the MYB intron I regulatory sequence increase transcription in colon cancers

Although MYB overexpression in colorectal cancer (CRC) is known to be a prognostic indicator for poor survival, the basis for this overexpression is unclear. Among multiple levels of MYB regulation, the most dynamic is the control of transcriptional elongation by sequences within intron I. The authors have proposed that this regulatory sequence is transcribed into an RNA stem-loop and 19-residue polyuridine tract, and is subject to mutation in CRC. When this region was examined in colorectal and breast carcinoma cell lines and tissues, the authors found frequent mutations only in CRC. It was determined that these mutations allowed increased transcription compared with the wild type sequence. These data suggest that this MYB regulatory region within intron I is subject to mutations in CRC but not breast cancer, perhaps consistent with the mutagenic insult that occurs within the colon and not mammary tissue. In CRC, these mutations may contribute to MYB overexpression, highlighting the importance of noncoding sequences in the regulation of key cancer genes. (c) 2006 Wiley-Liss, Inc.

Pseudo-messenger RNA: Phantoms of the transcriptome

The mammalian transcriptome harbours shadowy entities that resist classification and analysis. In analogy with pseudogenes, we define pseudo-messenger RNA to be RNA molecules that resemble protein- coding mRNA, but cannot encode full-length proteins owing to disruptions of the reading frame. Using a rigorous computational pipeline, which rules out sequencing errors, we identify 10,679 pseudo - messenger RNAs ( approximately half of which are transposonassociated) among the 102,801 FANTOM3 mouse cDNAs: just over 10% of the FANTOM3 transcriptome. These comprise not only transcribed pseudogenes, but also disrupted splice variants of otherwise protein- coding genes. Some may encode truncated proteins, only a minority of which appear subject to nonsense- mediated decay. The presence of an excess of transcripts whose only disruptions are opal stop codons suggests that there are more selenoproteins than currently estimated. We also describe compensatory frameshifts, where a segment of the gene has changed frame but remains translatable. In summary, we survey a large class of non- standard but potentially functional transcripts that are likely to encode genetic information and effect biological processes in novel ways. Many of these transcripts do not correspond cleanly to any identifiable object in the genome, implying fundamental limits to the goal of annotating all functional elements at the genome sequence level.

Sequence variation in the Newcastle disease virus genome

Full-length genome sequences of five virulent and five avirulent strains of Newcastle disease virus isolated between 1998 and 2002 in Victoria and New South Wales, Australia were determined. Comparisons between these strains revealed that coding sequence variability in the haemagglutinin-neuraminidase (HN), matrix (M) and phosphoprotein (P) gene sequences appeared to be more variable than in the fusion (F), nucleocapsid (N) and RNA dependent-RNA replicase (L) genes. Sequence analysis of a number of other isolates made during the recent virulent NDV outbreaks, also identified the presence of a number of variants with altered F gene cleavage sites, which resulted in altered biological properties of those viruses. Quasispecies analysis of a number of field isolates indicated the presence of virulent virus in one particular isolate. Gene sequence analysis of the progenitor virus isolated in 1998 showed very little sequence variation when compared to that of a progenitor-like virus isolated in 2001 demonstrating that in the field. viral genome sequence variation appears to be biologically restricted to that of a consensus sequence. (c) 2005 Elsevier B.V. All rights reserved.

Binding of an RNA trafficking response element to heterogeneous nuclear ribonucleoproteins A1 and A2

Heterogeneous nuclear ribonucleoprotein (hnRNP) A2 binds a 21-nucleotide myelin basic protein mRNA response element, the A2RE, and A2RE-like sequences in other localized mRNAs, and is a trans-acting factor in oligodendrocyte cytoplasmic RNA trafficking. Recombinant human hnRNPs A1 and A2 were used in a biosensor to explore interactions with A2RE and the cognate oligodeoxyribonucleotide. Both proteins have a single site that bound oligonucleotides with markedly different sequences but did not bind in the presence of heparin. Both also possess a second, specific site that bound only A2RE and was unaffected by heparin, hnRNP A2 bound A2RE in the latter site with a K-d near 50 nM, whereas the K-d for hnRNP A1 was above 10 muM. UV cross-linking assays led to a similar conclusion. Mutant A2RE sequences, that in earlier qualitative studies appeared not to bind hnRNP A2 or support RNA trafficking in oligodendrocytes, had dissociation constants above 5 muM for this protein. The two concatenated RNA recognition motifs (RRMs), but not the individual RRMs, mimicked the binding behavior of hnRNP A2. These data highlight the specificity of the interaction of A2RE with these hnRNPs and suggest that the sequence-specific A2RE-binding site on hnRNP A2 is formed by both RRMs acting in cis.

Forced retroevolution of an RNA bacteriophage

The operator hairpin ahead of the replicase gene in RNA bacteriophage MS2 contains overlapping signals for binding the coat protein and ribosomes. Coat protein binding inhibits further translation of the gene and forms the first step in capsid formation. The hairpin sequence was partially randomized to assess the importance of this structure element for the bacteriophage and to monitor alternative solutions that would evolve on the passaging of mutant phages. The evolutionary reconstruction of the operator failed in the majority of mutants. Instead, a poor imitation developed containing only some of the recognition signals for the coat protein. Three mutants were of particular interest in that they contained double nonsense codons in the lysis reading frame that runs through the operator hairpin. The simultaneous reversion of two stop codons into sense codons has a very low probability of occurring. Therefore the phage solved the problem by deleting the nonsense signals and, in fact, the complete operator, except for the initiation codon of the replicase gene. Several revertants were isolated with activities ranging from 1% to 20% of wild type. The operator, long thought to be a critical regulator, now appears to be a dispensable element. In addition, the results indicate how RNA viruses can be forced to step back to an attenuated form.

Long-range translational coupling in single-stranded RNA bacteriophages: an evolutionary analysis

In coliphage MS2 RNA a long-distance interaction (LDI) between an internal segment of the upstream coat gene and the start region of the replicase gene prevents initiation of replicase synthesis in the absence of coat gene translation. Elongating ribosomes break up the repressor LDI and thus activate the hidden initiation site. Expression studies on partial MS2 cDNA clones identified base pairing between 1427-1433 and 1738-1744, the so-called Min Jou (MJ) interaction, as the molecular basis for the long-range coupling mechanism. Here, we examine the biological significance of this interaction for the control of replicase gene translation. The LDI was disrupted by mutations in the 3'-side and the evolutionary adaptation was monitored upon phage passaging. Two categories of pseudorevertants emerged. The first type had restored the MJ interaction but not necessarily the native sequence. The pseudorevertants of the second type acquired a compensatory substitution some 80 nt downstream of the MJ interaction that stabilizes an adjacent LDI. In one examined case we confirmed that the second site mutations had restored coat-replicase translational coupling. Our results show the importance of translational control for fitness of the phage. They also reveal that the structure that buries the replicase start extends to structure elements bordering the MJ interaction.

LacI-mediated sequence-specific affinity purification of plasmid DNA for therapeutic applications

Affinity purification of plasmid DNA is an attractive option for the biomanufacture of therapeutic plasmids, which are strictly controlled for levels of host protein, DNA, RNA, and endotoxin. Plasmid vectors are considered to be a safer alternative than viruses for gene therapy, but milligram quantities of DNA are required per dose. Previous affinity approaches have involved triplex DNA formation and a sequence-specific zinc finger protein. We present a more generically applicable protein-based approach, which exploits the lac operator, present in a wide diversity of plasmids, as a target sequence. We used a GFP/His-tagged Lacl protein, which is precomplexed with the plasmid, and the resulting complex was immobilized on a solid support (TALON resin). Ensuing elution gives plasmid DNA, in good yield (>80% based on recovered starting material, 35-50% overall process), free from detectable RNA and protein and with minimal genomic DNA contamination. Such an affinity-based process should enhance plasmid purity and ultimately, after appropriate development, may simplify the biomanufacturing process of therapeutic plasmids.

Evaluation of DNA enzymes targeted against the RNA component of human telomerase

Glioblastoma Multiforme (GBM) is a highly malignant form of brain cancer for which there is currently no effective cure. Consequently, developing new therapies and elucidating effective targets is crucial for this fatal disease. In recent years, DNA enzymes, deoxyribonucleic acid molecules with enzymatic activity, have emerged. In the same manner as ribozymes, DNA enzymes are able to effect cleavage of RNA in a sequence-specific manner, and operate with catalytic efficiency. In this study, two DNA enzymes were designed to target the template region of human telomerase RNA (hTR), utilising the 10-23 and 8-17 catalytic motifs elucidated by Santoro and Joyce (1997). Telomerase is an RNA-dependent DNA polymerase, which stabilises telomere lengths by adding hexameric repeats (TTAGGG in humans) to chromosome termini, thus preventing the telomere shortening that usually occurs during mitotic cell division. Telomerase activity, whilst absent in normal somatic tissues, is present in almost 90% of all tumours. Thus, there is speculation that telomerase may be the much sought universal target for therapeutic intervention in cancer. In vitro cleavage assays showed both DNA enzymes to be catalytically competent. Unmodified phosphodiester (PO) backbone DNA enzymes were rapidly degraded in the presence of serum, with a half-life of 10 minutes. The common approach of introducing phosphorothioate (PS) linkages was used in an effort to overcome this instability. As a result of concurrent activity and stability studies on the DNA enzymes with various numbers of PS linkages, the DNA enzymes with a PO core and PS arms were chosen for use in further cell work. The cleavage activity of both was shown to be specific and affected by temperature, pH, MgCI2 concentration and enzyme concentration. Both DNA enzyme motifs reduced telomerase activity in cell lysates, as assessed by the telomerase repeat amplification protocol (TRAP) with an IC50 of 100nM. DNA enzymes being polyanionic molecules do not readily cross biological barriers. Cellular association of naked DNA enzyme was inefficient at less than 2%. Cellular delivery of the DNA enzymes was effectively improved using commercial cationic lipid formulations. However, the lipid-mediated delivery of DNA enzymes to U87-MG cells over a 4-hour period did not significantly inhibit cell proliferation compared to controls. This is possibly due to an expected lag period between the inhibition of telomere maintenance and cell death. Therefore, biodegradable polymer microspheres were investigated as a potential delivery option for prolonged and sustained delivery. In vitro release profiles showed that after an initial burst, sustained release of DNA enzymes was observed over 35 days. Finally, the efficacy and specificity of the DNA enzymes were demonstrated in a luciferase based reporter assay. Specific inhibition of luciferase expression was displayed at 10nM. Thus DNA enzymes have potential against endogenous cellular targets.

On the structural repertoire of pools of short, random RNA sequences

A detailed knowledge of the mapping between sequence and structure spaces in populations of RNA molecules is essential to better understand their present-day functional properties, to envisage a plausible early evolution of RNA in a prebiotic chemical environment and to improve the design of in vitro evolution experiments, among others. Analysis of natural RNAs, as well as in vitro and computational studies, show that certain RNA structural motifs are much more abundant than others, pointing out a complex relation between sequence and structure. Within this framework, we have investigated computationally the structural properties of a large pool (10 molecules) of single-stranded, 35 nt-long, random RNA sequences. The secondary structures obtained are ranked and classified into structure families. The number of structures in main families is analytically calculated and compared with the numerical results. This permits a quantification of the fraction of structure space covered by a large pool of sequences. We further show that the number of structural motifs and their frequency is highly unbalanced with respect to the nucleotide composition: simple structures such as stem-loops and hairpins arise from sequences depleted in G, while more complex structures require an enrichment of G. In general, we observe a strong correlation between subfamilies-characterized by a fixed number of paired nucleotides-and nucleotide composition. Our results are compared to the structural repertoire obtained in a second pool where isolated base pairs are prohibited. © 2008 Elsevier Ltd. All rights reserved.

Efficient HIV-1 inhibition by a 16 nt-long RNA aptamer designed by combining in vitro selection and in silico optimisation strategies

The human immunodeficiency virus type-1 (HIV-1) genome contains multiple, highly conserved structural RNA domains that play key roles in essential viral processes. Interference with the function of these RNA domains either by disrupting their structures or by blocking their interaction with viral or cellular factors may seriously compromise HIV-1 viability. RNA aptamers are amongst the most promising synthetic molecules able to interact with structural domains of viral genomes. However, aptamer shortening up to their minimal active domain is usually necessary for scaling up production, what requires very time-consuming, trial-and-error approaches. Here we report on the in vitro selection of 64 nt-long specific aptamers against the complete 5' -untranslated region of HIV-1 genome, which inhibit more than 75% of HIV-1 production in a human cell line. The analysis of the selected sequences and structures allowed for the identification of a highly conserved 16 nt-long stem-loop motif containing a common 8 nt-long apical loop. Based on this result, an in silico designed 16 nt-long RNA aptamer, termed RNApt16, was synthesized, with sequence 5'-CCCCGGCAAGGAGGGG-3-'. The HIV-1 inhibition efficiency of such an aptamer was close to 85%, thus constituting the shortest RNA molecule so far described that efficiently interferes with HIV-1 replication.

Sequence analysis of complementary DNAs encoding C3 in the nurse shark ({\it Ginglymostoma cirratum\/})

Mammalian C3 is a pivotal complement protein, encoded for by a single gene. In some vertebrate species multiple C3 isoforms are products of different C3 genes. The goal of this study was to determine whether multiple genes encode for shark C3. A protocol was developed for the isolation of mRNA from shark blood for the isolation of C3 cDNA clones. RT-PCR amplification of mRNA, using sense (GCGEQNM) and antisense (TWLTAYV) primers encoding conserved regions of human C3, yielded 21 clones. The C3-like clones isolated shared 97% similarity with each other and 40% similarity to human C3. RACE-PCR amplification of shark liver RNA, using gene specific primers, yielded products ranging from 1800bp to 3000bp. Deduced amino acid sequence, corresponding to 408bp of the 1800bp fragment, was obtained which showed 51% similarity to human C3. These results suggest that nurse shark C3 might be encoded for by more than one gene. ^

Sequence-related structural DNA anomalies affect restriction enzyme activity and DNA polymerase I binding

DNA serves as a target molecule for several types of enzymes and may assume a wide variety of structural motifs depending upon the local sequence. The BssHII restriction site (GC)3 resides in a 9bp region of alternating pyrimidine and purine residues within the &phis;X174 genome. Such sequences are known to demonstrate non-canonical helical behavior under the appropriate conditions. The kinetics of BssHII cleavage was investigated in supercoiled and linear plasmid DNA, and in a 323bp DNA fragment obtained via amplification of &phis;X174. The rate of enzyme cleavage was enhanced in the supercoiled form and in the presence of 50μM cobalt hexamine. Similarly, cobalt hexamine was also found to enhance TaqI activity directly adjacent to the (GC)3 region. ^ Initial DNA polymerase I binding studies (including a gel mobility shift assay and a protection assay) indicated a notable interaction between DNA polymerase I and the BssHII site. An in-depth study revealed that equilibrium binding of DNA polymerase I to the T7 RNA polymerase promoter was comparable to that of the (GC)3 site, however the strongest interaction was observed with a cruciform containing region. Increasing the ionic strength of the solution environment, including the addition of DNA polymerase I reaction buffer significantly decreased the equilibrium dissociation constant values. ^ It is suggested that the region within or around the BssHII site experiences a conformational change generating a novel structure under the influence of supercoiled tension or 50μM cobalt hexamine. It is proposed that this transition may enhance enzyme activity and binding by providing an initial enzyme-docking site—the rate-limiting step in restriction enzyme kinetics. The high binding potential of DNA polymerase I for each of the motifs described, is hypothesized to be due to recognition of the structural DNA anomalies by the 3′–5′ exonuclease domain. ^