Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “RNA sequencing analysis”

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 919 records · Page 51Linked to original sources

Large-scale gene expression profiling of discrete brain regions: potential, limitations, and application in genetics of aggressive behavior.

Many behavioral geneticists are interested in unraveling the molecular mechanisms underlying aggressive behavior. So far, most scientists have based their search for aggression-related genes on a preliminary functional hypothesis. Large-scale gene expression profiling techniques, such as serial analysis of gene expression (SAGE) and DNA microarrays, now enable the screening of expression levels of thousands of genes simultaneously, allowing the identification of new candidate aggression-related genes expressed in brain and thus the generation of new hypotheses. However, expression profiling in the brain is challenging, as brain is a complex heterogeneous tissue where large numbers of genes are expressed and relatively small changes in gene expression occur. In this special issue, we focus on the principles of SAGE and DNA microarrays, as well as their advantages and disadvantages and application to analysis in brain tissue in order to identify aggression-related genes.

Aggression↗

Phylogenetic analysis of Histoplasma capsulatum based on partial sequence of the D1/D2 region of the 28S rRNA gene.

In order to confirm the phylogenetic relationships of Histoplasma capsulatum, the partial sequences of large subunit (28S) ribosomal gene (D1/D2 region) of 49 isolates were studied. The similarity values of the 49 isolates were more than 99.0% across 617 base pairs, however, the 49 isolates were divided into 9 groups. These 9 groups were independent of 3 varieties, var. capsulatum, var. farciminosum and var. duboisii. These results showed that analysis of the nucleotide sequence of the 28S rRNA gene was very effective for identification of H. capsulatum and that three varieties of H. capsulatum should be reclassified according to the phylogenetic relationship established from analysis of the D1/D2 region sequences.

Animals↗

Ribosomal RNA analysis of Babesia odocoilei isolates from farmed reindeer (Rangifer tarandus tarandus) and elk (Cervus elaphus canadensis) in Wisconsin.

Piroplasms isolated from a farmed reindeer and elk in Wisconsin were determined to be Babesia odocoilei, based on morphology and ribosomal RNA (rRNA) analysis. Different clinical manifestations were observed in the two host species. The reindeer was parasitemic and exhibited acute babesiosis resulting in death, while the elk showed no parasites in blood smears and no overt clinical signs of babesiosis. B. odocoilei was, however, readily cultured from elk erythrocytes. Small subunit rRNA gene sequences from the two isolates were identical to that previously reported for B. odocoilei. Internal transcribed spacers 1 and 2 and 5.8S rRNA sequence analysis showed an overall identity range of 94.3-98.1% to corresponding sequences from three previously reported B. odocoilei isolates, but the Wisconsin reindeer B. odocoilei shared only 87.3% identity with a previously reported Babesia sp. isolated from a reindeer in California (RD61).

Animals↗

Comparison of the beet necrotic yellow vein virus P75 nucleotide sequences of Belgian type A and type B sources.

RNA-2 P75 open reading frame sequencing was conducted on BNYVV-infected sugar beet plants originating from 10 Belgian sources and was compared with sequences originating from France, Kazakhstan, Japan and China. This allowed the characterization of types A and B in Belgium. Sequence analysis confirmed that type A is intermediate between types P and B. The analysis of nucleotide sequences of the capsid protein epitopes and of different motifs involved in the transmission of the virus to the plant by Polymyxa betae reveals that they are well conserved among the different sequences analysed. In silico analysis showed that a potential dimerization/dsRNA-binding domain is present on P75. Finally, in order to compare the different BNYVV types, a real-time quantitative (RT)-PCR targeting the readthrough domain of P75 was developed and assessed over a period of 40 days on BNYVV susceptible sugar beet cv. Cadyx.

Amino Acid Sequence↗

A deep learning framework for denoising and ordering scRNA-seq data using adversarial autoencoder with dynamic batching.

Single-cell RNA sequencing (scRNA-seq) provides high resolution of cell-to-cell variation in gene expression and offers insights into cell heterogeneity, differentiating dynamics, and disease mechanisms. However, technical challenges such as low capture rates and dropout events can introduce noise in data analysis. Here, we present a deep learning framework, called the dynamic batching adversarial autoencoder (DB-AAE), for denoising scRNA-seq datasets. First, we describe steps to set up the computing environment, training, and tuning. Then, we depict the visualization of the denoising results. For complete details on the use and execution of this protocol, please refer to Ko et al.1.

Deep Learning↗

Gill-associated virus of Penaeus monodon prawns: an invertebrate virus with ORF1a and ORF1b genes related to arteri- and coronaviruses.

A 20089 nucleotide (nt) sequence was determined for the 5' end of the (+)-ssRNA genome of gill-associated virus (GAV), a yellow head-like virus infecting Penaeus monodon prawns. Clones were generated from a approximately 22 kb dsRNA purified from lymphoid organ total RNA of GAV-infected prawns. The region contains a single gene comprising two long overlapping open reading frames, ORF1a and ORF1b, of 4060 and 2646 amino acids, respectively. The ORFs are structurally related to the ORF1a and ORF1ab polyproteins of coronaviruses and arteriviruses. The 99 nt overlap between ORF1a and ORF1b contains a putative AAAUUUU 'slippery' sequence associated with -1 ribosomal frameshifting. A 131 nt stem-loop with the potential to form a complex pseudoknot resides 3 nt downstream of this sequence. Although different to the G/UUUAAAC frameshift sites and 'H-type' pseudoknots of nidoviruses, in vitro transcription/translation analysis demonstrated that the GAV element also facilitates read-through of the ORF1a/1b junction. As in coronaviruses, GAV ORF1a encodes a 3C-like cysteine protease domain located between two hydrophobic regions. However, its sequence suggests some structural relationship to the chymotrypsin-like serine proteases of arteriviruses. ORF1b encodes homologues of the 'SDD' polymerase, which among (+)-RNA viruses is unique to nidoviruses, as well as metal-ion-binding and helicase domains. The presence of a dsRNA replicative intermediate and ORF1a and ORF1ab polyproteins translated by a-1 frameshift suggests that GAV represents the first invertebrate member of the Order NIDOVIRALES:

Amino Acid Sequence↗

Isolation and characterization of cDNA clones for carrot extensin and a proline-rich 33-kDa protein.

Extensins are hydroxyproline-rich glycoproteins associated with most dicotyledonous plant cell walls. To isolate cDNA clones encoding extensin, we started by isolating poly(A)(+) RNA from carrot root tissue, and then translating the RNA in vitro, in the presence of tritiated leucine or proline. A 33-kDa peptide was identified in the translation products as a putative extensin precursor because: (i) it is rich in proline and poor in leucine, and (ii) the message appears to be more abundant when carrot tissue is wounded. From a cDNA library constructed with poly(A)(+) RNA from wounded carrots, one cDNA clone (pDC5) was identified that specifically hybridized to poly(A)(+) RNA encoding this 33-kDa peptide. We isolated three cDNA clones (pDC11, pDC12, and pDC16) from another cDNA library using pDC5 as a probe. DNA sequence data, RNA hybridization analysis, and hybrid released in vitro translation indicate that the cDNA clone pDC11 encodes extensin and that cDNA clones pDC12 and pDC16 encode the 33-kDa peptide, which as yet has an unknown identity and function. The assumption that the 33-kDa peptide was an extensin precursor was invalid. RNA hybridization and DNA sequence analysis indicate that pDC5 is a hybrid clone corresponding to two RNA species. RNA hybridization analysis showed that RNA encoded by both clone types is accumulated upon wounding.

Journal Article↗

Structure clustering features on the Sfold Web server.

UNLABELLED: The energy landscape of RNA secondary structures is often complex, and the Boltzmann-weighted ensemble usually contains distinct clusters. Furthermore, the minimum free energy structure often lies outside of the cluster containing the structure determined by comparative sequence analysis. We have developed procedures to characterize and visualize the Boltzmann-weighted ensemble, and have made them available on the Sfold Web server. The new features on the Web server include clustering statistics, ensemble and cluster centroids, multi-dimensional scaling display and energy landscape representation of the Boltzmann-weighted ensemble. AVAILABILITY: http://sfold.wadsworth.org; http://www.bioinfo.rpi.edu/applications/sfold CONTACT: chanc@wadsworth.org.

Algorithms↗

Characterization of RNA binding specificity of the Drosophila sex-lethal protein by in vitro ligand selection.

The Drosophila sex-lethal (Sxl) protein, a regulator of somatic sexual differentiation, is an RNA binding protein with two potential RNA recognition motifs (RRMs). It is thought to exert its function on splicing by binding to specific RNA sequences within Sxl and transformer (tra) pre-mRNAs. To examine the Sxl RNA binding specificity in detail, we performed in vitro selection and amplification of ligand RNAs from a random sequence pool on the basis of affinity with Sxl protein. After three cycles of selection and amplification, we cloned and sequenced 17 cDNAs corresponding to the RNAs selected in vitro. Sequencing showed that most of the RNAs selected contain polyuridine stretches surrounded by purine residues. In vitro binding analysis revealed that the sequences of the in vitro selected RNAs with relatively high affinity for Sxl show similarity to that of the Sxl- and tra-regulated acceptor regions, including the invariant AG sequence for splicing. These results suggest that Sxl recognizes and preferentially binds to a polyuridine stretch with a downstream AG sequence.

Animals↗

MGAlignIt: A web service for the alignment of mRNA/EST and genomic sequences.

Splicing is a biological phenomenon that removes the non-coding sequence from the transcripts to produce a mature transcript suitable for translation. To study this phenomenon, information on the intron-exon arrangement of a gene is essential, usually obtained by aligning mRNA/EST sequences to their cognate genomic sequences. MGAlign is a novel, rapid, memory efficient and practical method for aligning mRNA/EST and genome sequences. We present here a freely available web service, MGAlignIt (http://origin.bic.nus.edu.sg/mgalign/mgalignit), based on MGAlign. Besides the alignment itself, this web service allows users to effectively visualize the alignment in a graphical manner and to perform limited analysis on the alignment output. The server also permits the alignment to be saved in several forms, both graphical and text, suitable for further processing and analysis by other programs.

Algorithms↗

Identification of Freezing-Responsive microRNAs and Their Targets in Chinese Jujube by Small RNA and Degradome Sequencing.

The jujube tree fruit remains a primary fruit in northern China, yet its geographical distribution and yield are significantly constrained by freezing stress during winter. Numerous studies have highlighted the pivotal regulatory function of microRNAs (miRNAs) in plant responses to low-temperature stress. Nevertheless, the specific miRNAs involved in the response to low temperatures and their associated gene networks in Ziziphus jujuba Mill are not well understood. In this investigation, we utilized high-throughput sequencing to analyze small RNA libraries from branches subjected to temperatures of 4 °C and -30 °C. Our analysis identified a total of 342 miRNAs, comprising 123 known miRNAs and 219 novel miRNAs. The differential expression analysis revealed that under low-temperature conditions, 177 miRNAs underwent significant changes. Among them, specific upregulation of miR319 in the less cold-resistant variety and miR6483 in sensitive variety was observed. By employing degradome sequencing, we identified a total of 1551 target genes corresponding to 3059 unique miRNA target interaction pairs involving 299 miRNAs. Functional analysis using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways indicated that these target genes are primarily associated with transcriptional regulation, metabolic pathways, and genetic information processing. Through a comprehensive analysis, we pinpointed 11 genes corresponding to 9 miRNAs that are implicated in jujube tree cold stress, and 7 target genes of 7 miRNAs were confirmed by 5'-RACE analysis. These miRNAs are likely to exert crucial regulatory functions in the context of jujube tree cold stress. This study is the first to systematically identify miRNAs and their target genes in the response of Ziziphus jujuba Mill to low-temperature stress, which provides important resources for in-depth analysis of the molecular mechanism of jujube tree cold resistance and for cold-resistant breeding.

Ziziphus↗

Phylogenetic position of Blastocystis hominis and of stramenopiles inferred from multiple molecular sequence data.

Blastocystis hominis, a parasite of the human intestine, has recently been positioned within stramenopiles by the small subunit rRNA phylogeny. To further confirm its phylogenetic position using multiple molecular sequence data, we determined the nucleotide sequences putatively encoding small subunit ribosomal RNA, cytosolic-type 70-kDa heat shock protein, translation elongation factor 2, and the non-catalytic 'B' subunit of vacuolar ATPase of B. hominis (HE87-1 strain). Moreover, we determined the translation elongation factor 2 sequence of an apicomplexan parasite, Plasmodium falciparum, that belongs to alveolates. The maximum likelihood analyses of small subunit rRNA and cytosolic-type 70-kDa heat shock protein clearly demonstrated that B. hominis (HE87-1 strain) is positioned within stramenopiles, being congruent with the previous small subunit rRNA analysis, including the sequences of B. hominis (Nand strain) and a Blastocystis isolate from guinea pig. Although no clear resolution among major eukaryotic groups was obtained by the individual phylogenies based on the four molecules analyzed here, a combined analysis of various molecules, including these, clearly indicated that Blastocystis/stramenopiles are the closest relatives of alveolates.

Animals↗

RNA-RNA interactions in the binding site of protein L24 on 23S ribosomal RNA of E. coli. II. Sequence analysis of the interacting fragments.

A ribonucleoprotein complex containing several RNA subfragments from the 5' part of 23S RNA was recovered after digestion of the reconstituted complex between 23S RNA and protein L24. It was suggested in the preceding paper that the RNA subfragments 4B, 10A and 9, which are widely separated in the sequence, strongly interact. These subfragments were previously partially sequenced by the classical fingerprinting methods. Their sequences have now been completed with rapid new RNA sequencing methods. We propose here a base-pairing model showing how these subfragments may interact with one another.

Base Sequence↗

Obstacles in quantifying A-to-I RNA editing by Sanger sequencing.

Adenosine-to-Inosine (A-to-I) RNA editing is the most prevalent type of RNA editing, in which adenosine within a completely or largely double-stranded RNA (dsRNA) is converted to inosine by deamination. RNA editing was shown to be involved in many neurological diseases and cancer; therefore, detection of A-to-I RNA editing and quantitation of editing levels are necessary for both basic and clinical biomedical research. While high-throughput sequencing (HTS) is widely used for global detection of editing events, Sanger sequencing is the method of choice for precise characterization of editing site clusters (hyper-editing) and for comparing levels of editing at a particular site under different environmental conditions, developmental stages, genetic backgrounds, or disease states. To detect A-to-I editing events and quantify them using Sanger sequencing, RNA samples are reverse transcribed, cDNA is amplified using gene-specific primers, and then sequenced. The chromatogram outputs are then compared to the genomic DNA sequence. As editing occurs in the context of dsRNA, the reverse transcription step is performed at a temperature as high as 65 °C, using thermostable reverse transcriptase to open double-stranded structures. However, this measure alone is insufficient for transcripts possessing long stems comprised of hundreds of nucleotide pairs. Consequently, the editing levels detected by Sanger sequencing are significantly lower than those obtained by HTS, and the amplification yield is low. We suggest that the reverse transcription is biased towards unedited transcripts, and the severity of the bias is dependent on the transcript's secondary structure. Here, we show how this bias can be significantly reduced to allow reliable detection of editing levels and sufficient product yield.

RNA Editing↗

The gene for the U6 small nuclear RNA in fission yeast has an intron.

The small nuclear RNAs (snRNAs) are a class of metabolically stable small RNAs present in the nuclei of eukaryotic cells. In mammalian cells, there are six major molecular species (U1 to U6 snRNA), which are complexed with proteins, forming small nuclear ribonucleoprotein particles, snRNPs. Of these, the U1, U2, U4, U5, U6 snRNPs are thought to participate in pre-mRNA splicing as part of the spliceosome. Here, we describe the characterization of the gene coding for the Schizosaccharomyces pombe U6 snRNA. Unexpectedly, the Schiz. pombe U6 RNA gene was found to contain an intron-like sequence of 50 base pairs. Northern blot analysis and RNA sequencing revealed that this intron-like sequence is precisely removed from the transcript. The mature U6 RNA of Schiz. pombe has 77% sequence homology with the mammalian U6 RNA. In Schiz. pombe, it is possible that U6 RNA is not only involved in pre-mRNA splicing, but is also a splicing substrate. This is the first report of an intron in a snRNA gene.

Base Sequence↗

OligoMatcher: analysis and selection of specific oligonucleotide sequences for gene silencing by antisense or siRNA.

UNLABELLED: OligoMatcher is a web-based tool for analysis and selection of unique oligonucleotide sequences for gene silencing by antisense oligonucleotides (ASOs) or small interfering RNA (siRNA). A specific BLAST server was built for analysing sequences of ASOs that target pre-mRNA in the cell nucleus. Tissue- and cell-specific expression data of potential cross-reactive genes are integrated in the OligoMatcher program, which allows biologists to select unique oligonucleotide sequences for their target genes in specific experimental systems. AVAILABILITY: The OligoMatcher web server is available at http://shelob.cs.iupui.edu:18081/oligomatch.php. The source code is freely available for non-profit use on request to the authors. CONTACT: Mathew Palakal (mpalakal@cs.iupui.edu) or Shuyu Li (li_shuyu_dan@lilly.com).

Algorithms↗

Bacterial phylogeny based on comparative sequence analysis.

Comparative sequence analysis of small subunit rRNA is currently one of the most important methods for the elucidation of bacterial phylogeny as well as bacterial identification. Phylogenetic investigations targeting alternative phylogenetic markers such as large subunit rRNA, elongation factors, and ATPases have shown that 16S rRNA-based trees reflect the history of the corresponding organisms globally. However, in comparison with three to four billion years of evolution the phylogenetic information content of these markers is limited. Consequently, the limited resolution power of the marker molecules allows only a spot check of the evolutionary history of microorganisms. This is often indicated by locally different topologies of trees based on different markers, data sets or the application of different treeing approaches. Sequence peculiarities as well as methods and parameters for data analysis were studied with respect to their effects on the results of phylogenetic investigations. It is shown that only careful data analysis starting with a proper alignment, followed by the analysis of positional variability, rates and character of change, testing various data selections, applying alternative treeing methods and, finally, performing confidence tests, allows reasonable utilization of the limited phylogenetic information.

Bacteria↗

Fragmentation heterogeneity of 23S ribosomal RNA in Haemophilus species.

The fragmentation of 23S rRNA of 23 Haemophilus influenzae strains and eight strains belonging to other Haemophilus species was investigated. Instead of intact molecules, the 23S rRNA molecules were found to be cleaved into two to five smaller conserved fragments in most strains examined, especially in H. influenzae type b (5/6) and nontypeable strains (5/5). One or two conserved potential cleavage sites were identified by PCR analysis of the strains showing a fragmented 23S rRNA pattern. The relevant nucleotide sequences were determined and compared to H. influenzae Rd, which contains intact 23S rRNA molecules. An identical 112bp long intervening sequence (IVS) at position 542 and a conserved 121-123bp IVS sequence at position 1171 were found in two H. influenzae type b strains and one nontypeable strain. Among the strains with fragmented 23S rRNA, nearly half showed a heterogeneous cleavage pattern due to the dispersion of IVSs among different 23S rRNA operons. The localization of the conserved H. influenzae IVSs coincided well with the extensively studied IVSs among other bacteria, but differed in nucleotide sequence from any other reported IVSs. Therefore, the IVSs of Haemophilus 23S rRNA may originate from a common source that is independent of other bacteria.

Base Sequence↗