Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “RNA sequencing”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44Linked to original sources

Direct RNA nanopore sequencing of full-length coronavirus genomes provides novel insights into structural variants and enables modification analysis.

Sequence analyses of RNA virus genomes remain challenging owing to the exceptional genetic plasticity of these viruses. Because of high mutation and recombination rates, genome replication by viral RNA-dependent RNA polymerases leads to populations of closely related viruses, so-called "quasispecies." Standard (short-read) sequencing technologies are ill-suited to reconstruct large numbers of full-length haplotypes of (1) RNA virus genomes and (2) subgenome-length (sg) RNAs composed of noncontiguous genome regions. Here, we used a full-length, direct RNA sequencing (DRS) approach based on nanopores to characterize viral RNAs produced in cells infected with a human coronavirus. By using DRS, we were able to map the longest (∼26-kb) contiguous read to the viral reference genome. By combining Illumina and Oxford Nanopore sequencing, we reconstructed a highly accurate consensus sequence of the human coronavirus (HCoV)-229E genome (27.3 kb). Furthermore, by using long reads that did not require an assembly step, we were able to identify, in infected cells, diverse and novel HCoV-229E sg RNAs that remain to be characterized. Also, the DRS approach, which circumvents reverse transcription and amplification of RNA, allowed us to detect methylation sites in viral RNAs. Our work paves the way for haplotype-based analyses of viral quasispecies by showing the feasibility of intra-sample haplotype separation. Even though several technical challenges remain to be addressed to exploit the potential of the nanopore technology fully, our work illustrates that DRS may significantly advance genomic studies of complex virus populations, including predictions on long-range interactions in individual full-length viral RNA haplotypes.

Cell Line↗

Pattern generation in molecular evolution: exploitation of the variation in RNA landscapes.

Evolution of RNA secondary structure is studied using simulation techniques and statistical analysis of fitness landscapes. The transition from RNA sequence to RNA secondary structure leads to fitness landscapes that have local variations in their "ruggedness." Evolution exploits these variations. In stable environments it moves the quasispecies toward relatively "flat" peaks, where not only the master sequence but also its mutants have a high fitness. In a rapidly changing environment, the situation is reversed; evolution moves the quasispecies to a region where the correlation between secondary structures of "neighboring" RNA sequences is relatively low. In selection for simple secondary structures the movement toward flat peaks leads to pattern generation in the RNA sequences. Patterns are generated at the level of polynucleotide frequencies and the distribution of purines and pyrimidines. The patterns increase the modularity of the sequence. They thereby prevent the formation of alternative secondary structures after mutations. The movement of the quasispecies toward relatively rugged parts of the landscape results in pattern generation at the level of the RNA secondary structure. The base-pairing frequency of the sequences increases. The patterns that are generated in the RNA sequences and the RNA secondary structures are not directly selected for and can be regarded as a side effect of the evolutionary dynamics of the system.

Animals↗

Complete nucleotide sequence of RNA 3 from alfalfa mosaic virus, strain S.

We report the sequence of RNA 3 from strain S of Alfalfa mosaic virus (2,055 nucleotides). This RNA codes for a 32.4 kd protein (P3) and for the 24 kd coat protein (P4). The largest part of the sequence was established using RNA sequencing methods. The completion of the sequence in the region coding for P3 was achieved with cloned cDNA synthesized after priming at internal sites of RNA 3. Comparison of the RNA sequences coding P3 and P4 proteins in strain S with those reported in the literature for strain 425 revealed a higher amino acid substitution rate (3%) for P3 than for P4 (congruent to 1%) despite a similar average base substitution of 3-4% in these regions. In P3, two out of nine amino acid changes occur in hydrophilic regions. The amino acid changes in P4 do not modify the local hydrophilicity distribution. The intercistronic region displays a low degree of base substitution (2%) when compared with the untranslated 3'-end region (3.6%) or the 5'-end leader region (8%), the average substitution rate being 3.2%.

Amino Acid Sequence↗

U1-snRNP-A protein selects a ten nucleotide consensus sequence from a degenerate RNA pool presented in various structural contexts.

The U1snRNP-A (U1-A) protein was used to select specific RNA sequences from a degenerate pool of transcripts using direct RNA binding and polymerase chain reaction amplification (PCR). Sequences were randomized in loops of 10 or 13 nucleotides or as a linear stretch of 25 nucleotides. From all three structural contexts, an unpaired ten nucleotide consensus sequence was obtained. A selected stem-loop structure that resembled the natural U1-A protein binding site on loop II of U1 RNA demonstrated the highest affinity of binding in comparison with the other structural contexts. A data profile of selected sequences identified U1 RNA upon searching the GenBank database. Thus, this method was useful in determining the sequence specificity of an RNA binding protein and may complement the use of phylogenetic comparisons to predict conserved recognition elements. These findings also suggest that the evolutionary conservation of loop II of U1 RNA results from constraints imposed by protein binding.

Base Sequence↗

Sex-lethal, a Drosophila sex determination switch gene, exhibits sex-specific RNA splicing and sequence similarity to RNA binding proteins.

The switch gene, Sex-lethal (Sxl), controls sexual development and dosage compensation. It must be active in females and inactive in males throughout development. Analysis of Sxl cDNAs shows that this on/off regulation may be explained by differential RNA splicing; only female transcripts appear to encode functional products, whereas all male transcripts contain an exon that truncates the open reading frame. The functional female product shows sequence similarities with ribonucleoproteins, suggesting that it is an RNA binding protein. Thus, we propose that Sxl encodes a factor that interacts with both its own pre-mRNA (accounting for positive autoregulation) and that of downstream genes to confer female-specific splicing. In this way, a single, simple mechanism could account for both the maintenance and expression of the sexually determined state.

Amino Acid Sequence↗

Protocol to decode the role of transcriptionally active microbes in SARS-CoV-2-positive patients using an RNA-seq-based approach.

The elucidation of the role of microorganisms in human infections has been hindered by difficulties using conventional culture-based techniques. Here, we present a protocol for the investigation of transcriptionally active microbes (TAMs) using an RNA sequencing (RNA-seq)-based approach. We describe the steps for RNA isolation, viral genome sequencing, RNA-seq library preparation, and metatranscriptomic and transcriptomic analysis. This protocol permits a comprehensive evaluation of TAMs' contributions to the differential severity of infectious diseases, with a particular focus on diseases such as COVID-19. For complete details on the use and execution of this protocol, please refer to Devi et al.1.

Humans↗

Conserved and variable domains within divergent rnase P RNA gene sequences of Prochlorococcus strains.

RNase P RNA gene (rnpB) sequences were PCR-amplified from different members of the Prochlorococcus group. Aligned nucleotide sequences revealed a variance of up to 27% for rnpB. Comparative secondary structure analysis showed that domains P12, P18 and P19 of these novel ribozyme sequences in particular are highly divergent. Thus, these regions in RNase P RNA might serve as potential targets for deoxyoligonucleotide primers for the identification of specific genotypes of Prochlorococcus and for probing field populations. Phylogenetic trees constructed from RNase P RNA sequences were similar to, but not fully congruent with, those derived previously using sequences of the 16S rRNA gene. However, the application of rnpB sequences allowed a better resolution within clades of very closely related genotypes. As is known from 16S rRNA-based phylogenetic trees, sequences from individual strains clustered according to their physiology and the conditions at the original site of isolation, rather than their geographical origin. All sequences obtained from high-light-adapted strains formed a single coherent clade, as did the four sequences from low-light-adapted strains that were previously isolated from the North Atlantic and the subtropical North Pacific. This suggests a remarkable genetic stability of Prochlorococcus genotypes that thrive under identical ecological conditions.

Base Sequence↗

RNA processing of beta-globin transcripts containing 5' flanking and structural gene sequences.

RNA species have been identified in murine erythroid cells which contain both 5' flanking and structural gene sequences from the beta maj globin gene. Two nonpolyadenylated RNA transcripts, average 3700 and 1800-1900 nucleotides long, were identified by denaturing agarose gel electrophoresis and were found to hybridize to both 5' and 3' beta maj globin flanking sequences. This finding suggests that transcription initiated in the 5' flanking region does proceed past the polyadenylation site. The apparent higher concentration of the 5' flanking sequences in precursor RNA molecules, as compared to mature polyadenylated globin mRNA, suggests increased precursor stability of beta globin transcripts initiated in the 5' flanking region.

Animals↗

Molecular phylogenetics of the four Schistosoma species groups determined with partial 28S ribosomal RNA gene sequences.

Partial 28S ribosomal RNA (rRNA) gene sequences, including the variable domains D1, D2 and D3, were determined for representative species from the 4 Schistosoma species groups. On an alignment of 1345 bp from S. mansoni, S. haematobium, S. spindale and S. japonicum (with Heterobilharzia americana chosen as an outgroup), both maximum likelihood and maximum parsimony analyses provide a robust molecular phylogeny for the genus; ((((S. haematobium, S. spindale), S. mansoni), S. japonicum), H. americana). When analysed separately, both domain D1 and domain D2 yielded similarly informative data whereas D3 failed to resolve the phylogeny. These results confirm a phylogeny previously suggested by 18S rRNA gene sequences, corroborating the status of S. spindale as a sister taxon to S. haematobium, and demonstrate the utility of 28S rRNA gene sequence data for resolving phylogenies within the Schistosomatidae.

Animals↗

Cloning and sequencing of the S RNA from a Bulgarian isolate of tomato spotted wilt virus.

Libraries of cloned cDNA were prepared from complete genomic RNA and isolated S RNA of the Bulgarian L3 isolate of tomato spotted wilt virus (TSWV-L3). Northern blotting of TSWV genomic RNA detected clones specific for the L, M and S RNAs in the library from complete RNA. S RNA-specific clones selected from both libraries covered approximately 2.8 kb (about 95%) of the S RNA. Sequencing of these clones showed TSWV-L3 S RNA to be ambisense. It contains two open reading frames (ORFs); one of 1401 nucleotides located on the viral RNA encodes an Mr 52,400 (52K) protein, and the other of 774 nucleotides on the complementary strand encodes an Mr 28,900 (29K) protein. Expression of the 29K ORF in bacteria and immunological analysis of the fusion protein synthesized confirmed that the 29K protein is the N protein of TSWV-L3. Comparison with the published sequence for the S RNA of a Brazilian TSWV isolate, CNPH1, revealed almost complete identity in the amino acid sequences for the 29K protein, but several clustered amino acid exchanges in the putative 52K protein. In addition, the separating non-translated intergenic region of the S RNA of the Bulgarian isolate is 81 nucleotides longer than that of CNPH1.

Amino Acid Sequence↗

Classification of acidophilic, neutrotolerant and neutrophilic streptomycetes by nucleotide sequencing of 5S ribosomal RNA.

Complete 5S ribosomal RNA sequences were obtained for four acidophilic actinomycetes, seven neutrophilic streptomycetes and a strain of Streptoverticillium baldaccii. All of the organisms contained RNAs belonging to the 120 nucleotide type. An evolutionary tree was generated after combining the test data with results from similar studies on representative Gram-positive bacteria. The acidophilic, neutrotolerant and neutrophilic actinomycetes were recovered in a distinct cluster that was equated with the genus Streptomyces. The sequence data support the view that the genera Chainia, Elytrosporangium, Kitasatoa and Microellobosporia should be considered as synonyms of the genus Streptomyces. The recovery of the Streptoverticillium baldaccii strain on the fringe of the Streptomyces cluster is also consistent with current trends in the taxonomy of these organisms. Further work is needed to determine the taxonomic status of the two streptomycete subgroups that comprised the streptomycete cluster.

Actinomycetales↗

Nucleotide sequences of Escherichia coli 16-S RNA associated with ribosomal proteins S7, S9, S10, S14 and S19.

1. A ribonucleoprotein fragment containing proteins S7, S9, S10, S14, and S19 was isolated in high yield from Escherichia coli 30-S ribosomal sub-particles. The same fragment was obtained whether ribosomes from E. coli strain A19 or MRE 600 were used, despite the fact that protein S7 differs widely between the two strains. RNA was extracted from this fragment and fractionated on gels containing 7 M urea, to reveal "hidden breaks". A well-defined and reporducible pattern of RNA fragments was obtained, with the main components being approximately 300, 240, 130, 115 and 75 nucleotides in length, respectively. The pattern of RNA fragments obtained was also independent of the strain of E. coli used. 2. Two-dimensional fingerprints were made from ribonuclease T1 hydrolysates of these RNA fragments, labelled with 32P, and the oligonucleotides were further analysed by digestion with either ribonuclease A or U2. The data obtained were fitted in detail to the new 16-S RNA sequence map of Ehresmann et al. (1975). Again no significant differences were observed between the RNA from E. coli A19 or MRE 600. The RNA sequences found lay in the region O'-D-E'-K-P-P'-E-A of the 16-S RNA, with a clear excision of several nucleotides in section E'-K. The total sequence length covered was approximately 430 nucleotides.

Base Sequence↗

A systematic search for RNA editing sites in pea chloroplasts: an editing event causes diversification from the evolutionarily conserved amino acid sequence.

RNA editing in higher plant chloroplasts involves C-to-U conversion at specific sites in the transcripts. To examine whether pea shares editing sites with other angiosperms, a systematic search for editing sites in pea chloroplast transcripts was performed. Based on amino acid sequence alignment, 451 RNA editing sites were predicted from 60 transcripts. Sequence analysis of amplified cDNAs for these potential editing sites revealed 19 true editing sites from 13 transcripts. Together with those reported previously, the total number of editing sites is 27 from 16 transcripts in pea chloroplasts. Twenty-two sites are conserved among other plant species, whereas five sites are unique to pea. Among the 27 editing sites, seven are partially edited. The most interesting is the ndhG site 1, which has led to the diversification of the evolutionarily conserved amino acid sequence. This observation suggests that some of the editing events cause the diversity of amino acid sequences, and hence, that prediction of editing sites based on amino acid sequence alignment has its own limitations.

Amino Acid Sequence↗

RNA editing of wheat mitochondrial ATP synthase subunit 9: direct protein and cDNA sequencing.

RNA editing of subunit 9 of the wheat mitochondrial ATP synthase has been studied by cDNA and protein sequence analysis. Most of the cDNA clones sequenced (95%) showed that editing by C-to-U transitions occurred at eight positions in the coding region. Consequently, 5 amino acids were changed in the protein when compared with the sequence predicted from the gene. Two edited codons gave no changes (silent editing). One of the C-to-U transitions generated a stop codon by modifying the arginine codon CGA to UGA. Thus, the protein produced is 6 amino acids shorter than that deduced from the genomic sequence. Minor forms of cDNA with partial or overedited sequences were also found. Protein sequence and amino acid composition analyses confirmed the results obtained by cDNA sequencing and showed that the major form of edited atp9 mRNA is translated.

Adenosine Triphosphatases↗

Suprageneric classification of peptidoglycan group B actinomycetes by nucleotide sequencing of 5S ribosomal RNA.

5S ribosomal RNA sequences were determined for thirteen actinomycetes mainly representatives with the rare group B type peptidoglycan. The primary and secondary structure of the resultant sequences were of the type characteristic of Gram-positive bacteria with DNA rich in guanine plus cytosine. The sequencing and associated chemotaxonomic data provide compelling grounds for classifying actinomycetes with a group B type peptidoglycan in a single family. The family Microbacteriaceae fam. nov. is proposed to accommodate actinomycetes classified in the genera Agromyces, Aureobacterium, Clavibacter, Curtobacterium and Microbacterium.

Actinomyces↗

Sequence conservation and divergence of hepatitis delta virus RNA.

The complete RNA sequence of the hepatitis delta virus (HDV) obtained from the Nauru Island in the Pacific was determined by cDNA cloning and amplification by polymerase chain reaction (PCR). The sequence showed 14-17% divergence from the two known HDV RNA sequences. There are three highly conserved domains: the region around the autocatalytic cleavage site of the genomic RNA (nucleotides 659 to 772), the region around the autocatalytic cleavage site of the antigenomic-sense RNA (nucleotides 847 to 966), and the region around the middle one-third domain of the open reading frame (ORF) encoding the hepatitis delta antigen on the antigenomic RNA (nucleotides 1267 to 1347). The two autocatalytic activities are required for the cleavage and ligation of HDV RNA during RNA replication. The third conserved domain codes for the RNA-binding domain of HDAg, which specifically interacts with HDV RNA. Three nucleotide changes within the genomic catalytic sequence are present but did not alter the catalytic cleavage activity of the HDV RNA. Microheterogeneity of the RNA sequences was also detected. One of these occurred within the coding region of the delta antigen, creating an amber termination codon in some of the RNA species. Thus, this HDV strain contains two different RNA species, one of which encodes a delta antigen of 214 amino acids and the other 195 amino acids. These two protein species were detected by immunoblotting of the patient's plasma. In contrast to other HDV strains, only three ORFs capable of encoding more than 100 amino acids each are present in this HDV RNA. We recommend that oligonucleotides complementary to the highly conserved sequences should be used as primers for PCR in clinical detection assays of hepatitis delta virus infection.

Adult↗