Search PubMedSearch

SEARCH · Search PubMed

Results for “Direct 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 37 records · Page 2Linked to original sources

Analysis of Leishbuviridae from Trypanosomatids.

Over the last decade, considerable progress has been made in unraveling RNA virus diversity. This has contributed to our understanding of the evolution of these viruses, which include emerging zoonotic human pathogens. Current success has been greatly facilitated by the development of next-generation sequencing platforms instrumental for meta-transcriptomic studies. However, due to the rapid evolution of RNA viruses, there are numerous "blind spots" waiting to be explored; one of those is the RNA virome of unicellular eukaryotes. Here, we present the pipeline, which has been successfully used to characterize various types of RNA viruses, including Leishbuviridae (Bunyaviricetes, Hareavirales) in the parasitic flagellates of the family Trypanosomatidae. The pipeline relies on axenic in vitro cell culture and double-stranded RNA enrichment, followed by direct RNA-sequencing. A detailed procedure description starting from the initial total RNA preparation to the final assembly of the viral segments is provided.

High-Throughput Nucleotide Sequencing

A novel reusable transcriptome-wide association study workflow used to map key genes linked to important cattle traits.

Transcriptome-wide association studies (TWAS) are a powerful approach for studying the genes underlying complex traits by directly integrating GWAS and gene expression datasets. In cattle, they have been previously applied to identify genes driving fertility, milk production, and health. However, these studies have also highlighted several challenges, from difficulties in reproducing these complex analyses to limitations from poor genotype calls, especially when called directly from RNA sequencing data. To address these and other challenges, for the H2020 BovReg Project, we have developed a streamlined, species-agnostic, and reusable Nextflow TWAS workflow to integrate transcriptomic and GWAS summary statistic datasets. Our workflow first generates accurate genotype calls and gene expression prediction models from transcriptomic datasets and then applies these tools to impute gene expression levels into GWAS cohorts, enabling the association of genes with traits of interest. We explore optimal strategies for calling genetic variants directly from transcriptomic data and illustrate that using imputation approaches specifically designed for low-pass sequencing data can improve variant calling over previously adopted methods. We demonstrate the utility of our TWAS workflow by applying it to both novel and publicly available GWAS cohorts for cattle, detecting novel gene-trait associations for complex traits. Using a new transcriptome annotation of the cattle genome generated for the BovReg project we also illustrate how previously un-assayable associations can be detected. The results and the workflow we present, provide a new resource for the community and contribute to a better understanding of the molecular drivers of complex traits in cattle with the goal of eventually leveraging this information in future breeding decisions.

Animals

Enzymes in high-throughput RNA sequencing: Applications and challenges.

High-throughput RNA sequencing provides genome-wide information on the dynamics of RNA in each cell and how the dynamics responds to environmental changes. Next-generation sequencing by the Illumina platform currently provides the highest information output as compared to other platforms. A key component of next generation sequencing of each RNA is the successful end-to-end reverse-transcription into a cDNA strand. This can be highly challenging given the propensity of each RNA to adopt ordered structures and to contain post-transcriptional modifications. While many reverse transcriptase (RT) enzymes have been developed over the years to maximize read-through of an RNA, their processivity and efficiency varies, raising the question of how to select the RT for the experiment at hand. Here, we use tRNA as a model for genome-wide sequencing, as tRNA has a stable secondary and tertiary structure and has a high density and wide variety of post-transcriptional modifications, presenting one of the most challenging problems of sequencing RNA. We compare the efficiency of end-to-end cDNA synthesis of tRNA among several recent RT enzymes and provide a general sequencing workflow that is applicable to most of these enzymes.

High-Throughput Nucleotide Sequencing

The use of nuclease P1 in sequence analysis of end group labeled RNA.

A method is described for the direct sequence analysis of 20-25 nucleotides from the termini of 5'- or 3'-end-group [32P] labeled RNA. The method involves partial endonucleolytic digestion of the labeled RNA with nuclease P1 (from Penicillium citrinum) followed by separation of the partial digestion products by two-dimensional homochromatography, the nucleotide sequence being determined by mobility shift analysis. This procedure has been applied to the sequence analysis of the terminal regions of tRNAs and of high molecular weight RNA, such as messenger RNA or viral RNA. A further application involves its use in conjunction with snake venom phosphodiesterase to determine the sequence of 5'-end group labeled oligonucleotides, containing modified bases, derived from T1 or pancreatic RNase digestion of tRNA.

Base Sequence

The structure of a transcriptional unit on colicin E1 plasmid.

In an RNA-synthesizing system in vitro, a low-molecular-weight RNA consisting of about 110 residues (RNA-I) was efficiently synthesized on DNA of colicin E 1 plasmid (ColE1) and its deletion derivatives. The promoter site for RNA-I was analysed by testing the RNA polymerase-binding ability and template activity of restriction fragments; it was mapped in the region between the replication initiation site and the colicin immunity gene of ColE1. The direction of transcription was determined by hybridization tests to the separated strands of the template. The DNA region directing RNA-I was sequenced, and RNA-I was assigned on the sequence based on the nearest-neighbour data of RNA. The sequences of its promoter and terminator regions were also deduced. Although the function of this small RNA species is unknown, a unique secondary structure could be constructed from its sequence and sensitivity to RNase.

Bacteriocin Plasmids

Genome sequencing reveals the impact of pseudoexons in rare genetic disease.

PURPOSE: Advancements in sequencing technologies have significantly improved clinical genetic testing; yet, the diagnostic yield remains around 30% to 40%. Emerging technologies are now being deployed to address the remaining diagnostic gap. METHODS: We tested whether short-read genome sequencing could increase the diagnostic yield in individuals enrolled into the UCI-GREGoR research study, who had suspected Mendelian conditions and prior inconclusive testing. Two other collaborative research cohorts, focused on aortopathy and dilated cardiomyopathy, consisted of individuals who were undiagnosed but had not undergone harmonized prior testing. RESULTS: We sequenced 353 families (754 participants) and found a molecular diagnosis in 54 (15.3%) of them. Of these diagnoses, 55.5% were previously missed because the causative variants were in regions not originally interrogated. In 5 cases, they were deep intronic variants, all of which led to abnormal splicing and pseudoexons, as directly shown by RNA sequencing. All 5 of these variants had inconclusive spliceAI scores. In 26% of newly diagnosed cases, the causal variant could have been detected by exome sequencing reanalysis. CONCLUSION: Genome sequencing can overcome limitations of clinical genetic testing, such as the inability to call intronic variants. Our findings highlight pseudoexons as a common mechanism via which deep intronic variants cause Mendelian disease.

Humans

Model-directed generation of artificial CRISPR-Cas13a guide RNA sequences improves nucleic acid detection.

CRISPR guide RNA sequences deriving exactly from natural sequences may not perform optimally in every application. Here we implement and evaluate algorithms for designing maximally fit, artificial CRISPR-Cas13a guides with multiple mismatches to natural sequences that are tailored for diagnostic applications. These guides offer more sensitive detection of diverse pathogens and discrimination of pathogen variants compared with guides derived directly from natural sequences and illuminate design principles that broaden Cas13a targeting.

CRISPR-Cas Systems

Site specific enzymatic cleavage of RNA.

The hybridization of a DNA oligonucleotide a specific tetramer or longer) will direct a cleavage by RNase H (EC 3.1.4.34) to a specific site in RNA. The resulting fragments can then be labeled at their 5' or 3' ends, purified, and sequenced directly. This procedure is demonstrated with two RNA molecules of known sequence: 5.8S rRNA from yeast (158 nucleotides) and satellite tobacco necrosis virus (STNV) RNA (1240 nucleotides).

Animals

Mapping the spliced and unspliced late lytic SV40 RNAs.

The sizes and map positions of the major late lytic SV40 cytoplasmic mRNAs and the abundant nuclear RNA species have been determined by the technique of Berk and Sharp (1977, 1978). From these experiments, the coding sequences (bodies) of the 16S and 19S late cytoplasmic SV40 RNAs have been located at 0.935-0.17 and 0.765-0.17 map units, respectively. The cytoplasmic 16S RNA molecules contain a leader sequence of approximately 210 nucleotides, corresponding to SV40 map positions 0.72-0.76 units, spliced to the coding sequences. In a population of the late 19S RNA molecules, there are several different leader segments, each spliced to the same coding sequences. The size of these 19S leader RNA segments was estimated to be 50-70, 100-120 and 200-210 nucleotides in length. The 5' ends of the 19S leader RNA segments were located at 0.72, 0.71, 0.695 and 0.69 map units. An analysis of the nuclear viral RNAs has provided insight into the biogenesis of the cytoplasmic messages. Poly (A)-containing nuclear RNA has a number of species in addition to those found in the cytoplasm. The 3' ends of the poly (A)-containing RNAs map at 0.17 SV40 units. The 5' ends of the more abundant nuclear molecules map approximately at 0.72, 0.70, 0.67, 0.64 and 0.59 units. Since these nuclear SV40 RNA molecules are both colinear with the viral DNA and larger than the cytoplasmic nRNAs, they may represent intermediates in a stepwise processing system. Alternatively, the variation in 5' ends of nuclear SV40 transcripts may represent a number of separate initiation sites for transcription. The presence of the intervening RNA sequences (between the leader and the coding sequences of the mature mRNAs) in these nuclear RNA molecules suggests that the synthesis of "spliced" SV40 RNA involves the direct transcription of the DNA sequences and the subsequent splicing out of the intervening seqment of RNA. Evaluation of the more abundant nonpolyadenylated nuclear RNA molecules showed that they have the same 5' ends as the poly (A)-containing nuclear RNAs. The 3' ends of the nonpolyadenylated RNA molecules map heterogeneously in a broad region extending beyond 0.28 map units. The presence of these long nuclear viral transcripts suggests that transcription of late SV40 RNA does not terminate at 0.17 map units. The location of poly (A) in mature cytoplasmic viral RNA at 0.17 map units suggest that poly (A) addition to RNA molecules may occur by a specific cleavage of the longer transcripts. Based on these analyses, we propose that the longer nonpolyadenylated viral RNA molecules in the nuclei of SV40-infected cells may represent the primary transcripts. While their 5' termini are being processed, the specific addition of poly (A) at 0.17 map units takes place. The polyadenylation of RNA is followed by splicing events to generate the cytoplasmic forms of SV40 mRNA.

Base Sequence

Contacts between Escherichia coli RNA polymerase and thymines in the lac UV5 promoter.

I have identified those 5 positions of thymines in the lac UV5 promoter that lie close to bound Escherichia coli RNA polymerase (nucleosidetriphosphate:RNA nucleotidyltransferase, EC 2.7.7.6). Although ultraviolet irradiation of DNA with 5-bromouracil substituted in place of thymine normally cleaves the DNA at the bromouracils, a protein bound to the DNA can perturb these cleavages at those locations at which the protein lies close to the bromine. In the lac promoter most of these contacts lie in three regions. Four contacts lie in the region where transcription initiates; four lie in the "Pribnow box," which is located about 10 base pairs upstream from the initiation site; and three more lie in the "-35 region," located about 35 base pairs upstream from the initiation site. The "Pribnow box" and the "-35 region" are regions whose sequences are partially conserved between promoters and in which most promoter mutations are located; thus, contacts in these two regions probably represent sites of sequence-specific recognition by RNA polymerase.

Base Sequence

Nucleotide sequence at the 5' terminus of the avian sarcoma virus genome.

Transcription of DNA from the RNA genome of avian sarcoma virus by RNA-directed DNA polymerase in vitro initiates on a primer (tRNATrp) located near the 5'-terminus of the viral genome. One of the major products of transcription is a single-stranded DNA chain complementary to a sequence of 101 nucleotides immediately distal to the site of initiation of DNA synthesis. We have determined the complete nucleotide sequence of this transcribed chain for the Prague strain of avian sarcoma virus, a partial sequence of the transcribed chain for the Bratislava 77 strain of avian sarcoma virus, and the sequence of a DNA transcript that is shorter than the transcribed single-stranded chain. Our data define the location of tRNATrp on the genome of avian sarcoma virus and provide the sequence of 119 nucleotides at the 5'-terminus of the genome. Portions of this sequence may be involved in the binding of RNA-directed DNA polymerase, the initiation of translation from viral messenger RNA, the extension of RNA-directed DNA synthesis from the 5'- to the 3'-terminus of viral RNA, and the integration of viral DNA into the host genome.

Avian Sarcoma Viruses

Targeted loss of heterozygosity in Candida albicans using CRISPR-Cas9 reveals the functional impact of allelic variation.

The diploid genome of the fungal pathogen Candida albicans is highly heterozygous, with most allele pairs diverging at either the coding or regulatory level. When faced with selection pressure like antifungal exposure, this hidden genetic diversity can provide a reservoir of adaptive mutations through loss of heterozygosity (LOH) events. Validating the potential phenotypic impact of LOH events observed in clinical or experimentally evolved strains can be difficult due to the challenge of precisely targeting one allele over the other. Here, we show that a CRISPR-Cas9 system can be used to overcome this challenge. By designing allele-specific guide RNA sequences, we can induce targeted, directed LOH events, which we validate by whole-genome long-read sequencing. Using this approach, we efficiently recapitulate a recently described LOH event that increases resistance to the antifungal fluconazole. Additionally, we find that the recombination tracts of these induced LOH events have similar lengths to those observed naturally. To facilitate future use of this method, we provide a database of allele-specific sgRNA sequences for Cas9 that provide near genome-wide coverage of heterozygous sites through either direct or indirect targeting. This approach will be useful in probing the adaptive role of LOH events in this important human pathogen.

Candida albicans

Sequence of a RNA templated by the 3'-OH RNA terminus of defective interfering particles of vesicular stomatitis virus.

We have sequenced the endogenous RNA polymerase product produced by disrupted purified virions of vesicular stomatitis virus defective interfering particles by using the newer one-dimensional rapid gel sequencing techniques and confirming this with a modified two-dimensional gel vectoring technique. The sequence of this 46-nucleotide RNA is: 5'(pp)pACGAAGACCACAAAACCA-GAUAAAAAAUAAAAACCACAAGAGGG(U)COH3'. We infer that this sequence is identical to the sequence at the 5' end of infectious vesicular stomatitis virus RNA and is complementary to the sequence of the 3'-OH terminus of this defective interfering particle genome RNA.

Base Sequence

Distinctive nucleotide sequences of promoters recognized by RNA polymerase containing a phage-coded "sigma-like" protein.

We report the nucleotide sequences of two promoters for bacteriophage SP01 "middle" genes. These promoters are recognized by a modified form of Bacillus subtilis RNA polymerase that contains a phage-coded "sigma-like" regulatory protein (gp28) in place of the bacterial sigma factor. Both promoters shared the identical hexanucleotide 5'A-G-G-A-G-A at about 35 base pairs preceding the start point of transcription and the identical heptanucleotide 5'-T-T-T-A-T-T-T (T is the thymine analog 5-hydroxymethyluracil in SP01 DNA) located about 10 base pairs preceding the transcriptional start point. The significance of these sequences in comparison with nucleotide sequences of promoters recognized by sigma-containing RNA polymerases is discussed.

Bacillus subtilis

Cap and internal nucleotides of reovirus mRNA primers are incorporated into influenza viral complementary RNA during transcription in vitro.

Reovirus mRNA's containing a 5'-terminal methylated cap structure (m(7)GpppG(m)) were shown to be effective primers for influenza viral RNA transcription in vitro catalyzed by the influenza virion transcriptase. Priming activity required the presence of methyl groups in the cap since reovirus mRNA's with 5'-terminal GpppG were inactive as primers. Both the cap and internal nucleotides were physically transferred from radiolabeled reovirus mRNA to influenza viral complementary RNA (cRNA) during transcription in vitro. By using reovirus mRNA's with methyl-(3)H-labeled caps as primers, we showed that the influenza viral cRNA synthesized in the presence of unlabeled nucleoside triphosphates contained [methyl-(3)H]m(7)GpppG(m), identical to that found in the reovirus mRNA primer. To demonstrate transfer of internal residues, reovirus mRNA's synthesized in the presence of all four alpha-(32)P-labeled ribonucleoside triphosphates were used as primers. The resulting influenza viral cRNA was (32)P-labeled. Diethyl-aminoethyl-Sephadex chromatography of the RNase T2 digest of this cRNA demonstrated (32)P radiolabel in both internal residues (charge -2) and the cap (charge -4.6). Approximately 25 internal nucleotides along with the cap of reovirus mRNA were transferred to each chain of influenza viral cRNA. Gel electrophoretic analysis indicated that the segments of influenza viral cRNA primed by reovirus mRNA were approximately the same size as those primed by a different mRNA, globin mRNA, strongly suggesting that the influenza virion transcriptase complex transfers approximately the same number of nucleotides plus the cap from different mRNA primers to the 5' end of influenza viral RNA transcripts.

Base Sequence

Deletion mutant of the Bratislava-77 strain of Rous sarcoma virus containing a fusion of the group-specific antigen and envelope genes.

The genetic compositions of two independently derived preparations of the Bratislava-77 strain (B77) of Rous sarcoma virus were analyzed after each was passaged seven or more times in duck embryo fibroblasts. RNase, T1-resistant oligonucleotide fingerprint analysis of virion RNA from both preparations of duck-passaged B77 revealed the presence of two large noncontiguous deletions. Approximately 75% of the RNAs contained a deletion which spans oligonucleotides 304 to 4 on the viral genome (about 3,500 nucleotides) and encompasses all of the B77 polymerase gene. More than 90% of the RNAs also contained a deletion which spans src-specific oligonucleotides 6 and 5(about 2,200 nucleotides) and is identical to the deletion observed in transformation-defective B77. Virion RNA from duck-passaged B77 also contained two oligonucleotides (D1 and D2) not observed in the RNA of B77 virus grown on chicken embryo fibroblasts. Analysis of the virion RNA of duck-passaged B77 by denaturing agarose gel electrophoresis revealed four major subunits with molecular weights of 3.40 x 10(6), 2.65 x 10(6), 2.25 x 10(6), and 1.55 x 10(6). Whereas the 3.40- and 2.65-megadalton (Mdal) RNA species comigrated with the nondefective and transformation-defective RNAs of B77 propagated on chicken embryo fibroblasts, no counterparts to the 2.25- and 1.55-Mdal RNAs were observed in the RNA of B77 grown on chicken embryo fibroblasts. Oligonucleotide fingerprint analysis of these RNA species revealed that the 2.65-Mdal RNA contains the src-specific deletion and that 2.25-Mdal RNA contains the polymerase region deletion; both of these deletions were observed in the 1.55-Mdal RNA, which was the major RNA subunit species detected in duck-passaged B77. The new oligonucleotides (D1 and D2) observed in the duck-passaged virus were present in the 2.25- and 1.55-Mdal RNA species in vitro and in vivo and directs the synthesis of a 130,000-dalton protein (p130). p130 contains antigenic determinants specific for p27 (gag gene) and gp85 (env gene) but does not contain sequences which cross-react with antisera directed against the alpha beta form of RNA-dependent DNA polymerase (pol gene). This RNA, therefore, is generated by a fusion of the gag and env genes of Rous sarcoma virus B77.

Antigens, Viral