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 1,747 records · Page 97Linked to original sources

A comparative study of the ribosomal RNA operons of Streptomyces coelicolor A3(2) and sequence analysis of rrnA.

S. coelicolor A3(2) contains six ribosomal RNA operons. Here we describe the cloning of rrnA, rrnC and rrnE, thereby completing the cloning of all operons. Southern hybridisation of genomic DNA with a heterologous probe from the E.coli rrnB 16S rRNA gene showed differences in hybridisation among the six rRNA operon-containing bands. The nucleotide sequence of the 16S rRNA gene and the upstream region of rrnA was determined and compared with the corresponding sequence of rrnD, showing that the 16S rRNA genes are 99% identical. Substantial differences were found, however, in the upstream regions corresponding to the P1 and P2 promoters of rrnD. Southern analysis showed that some of the other rRNA operons of S.coelicolor A3(2) also differed in this part of the upstream region.

Amino Acid Sequence↗

Reverse transcription using random pentadecamer primers increases yield and quality of resulting cDNA.

Reverse transcription of RNA is an invaluable method for gene expression analysis by real-time PCR or microarray methods. Random primers of varying lengths were compared with respect to their efficiency of priming reverse transcription reactions. The results showed that 15-nucleotide-long random oligonucleotides (pentadecamers) consistently yielded at least 2-fold as much cDNA as did random hexamers using either poly(A) RNA or an amplified version of messenger RNA (aRNA) as a template. The cDNA generated using pentadecamers did not differ in size distribution or the amount of incorporated label compared with cDNA generated with random hexamers. The increased efficiency of priming using random pentadecamers resulted in reverse transcription of >80% of the template aRNA, while random hexamers induced reverse transcription of only 40% of the template aRNA. This suggests a better coverage of the transcriptome when using random pentadecamers over random hexamers. Using the same amount of aRNA as starting material, random pentadecamer-primed reactions resulted in 11-fold more genes being detected in whole transcriptome DNA microarray experiments than random hexamer-primed reactions. The results indicate that random pentadecamers can replace random hexamers in reverse transcription reactions on both poly(A) RNA and amplified RNA, resulting in higher cDNA yields and quality.

DNA Primers↗

Amaranth: enhanced single-cell transcript assembly via discriminative modelling of UMI reads and internal reads.

MOTIVATION: Single-cell RNA sequencing (scRNA-seq) has transformed transcriptome profiling at cellular resolution, yet accurate reconstruction of full-length transcripts for individual cells remains a central challenge. Emerging scRNA-seq protocols can produce reads that span entire transcripts, enabling isoform-level expression analysis. For example, Smart-seq protocols combine unique molecular identifier (UMI)-linked reads that index and stitch together multiple reads from the same molecule, with internal reads filling coverage gaps. We demonstrate that these read types exhibit markedly different biological and statistical properties in strandness, 5'/3' coverage bias, and genomic locality. Existing assemblers fail to leverage these distinctions, yielding suboptimal assembly. RESULTS: We developed Amaranth, a novel single-cell assembler that discriminatively models UMI and internal reads. Amaranth implements heuristics specifically designed to address the distinct biases of UMI-linked and internal reads, enabling accurate strandness assignment for internal reads, reliable splicing graph refinement, and precise transcript start site determination. We also developed Amaranth-meta, which integrates information across cells to enhance individual cell assemblies. Benchmarked on Smart-seq3 datasets from human HEK293T and mouse fibroblast cells, Amaranth outperformed other state-of-the-art assemblers in assembling individual cells and in meta-assembly. Amaranth advances isoform-level analysis in single-cell transcriptomics, facilitating detailed studies at cellular resolution. AVAILABILITY AND IMPLEMENTATION: Amaranth is implemented in C++ and is freely available at https://github.com/Shao-Group/amaranth under the BSD-3-Clause license. Scripts, documentation, and data for reproducing experiments in this manuscript are available at https://github.com/Shao-Group/amaranth-test.

Single-Cell Gene Expression Analysis↗

A pluripotent stem cell atlas of multilineage differentiation.

Human pluripotent stem cells offer a scalable platform to study genetic and signalling mechanisms governing cell lineage decisions during differentiation. Genome-wide and single-cell transcriptomics technologies likewise offer high-throughput analysis of heterogeneous cell differentiation states. While in vivo development has been extensively characterised using these technologies, there remains a need for comprehensive single-cell transcriptomic profiling of stem cell differentiation from pluripotency. Understanding gene expression changes governing differentiation in vitro is key to developing high fidelity differentiation protocols and understanding fundamental mechanisms of development. We generated a single-cell RNA sequencing time course to study the role of developmental signalling pathways on multilineage diversification from pluripotency in vitro. The combined dataset of over 60,000 cells spans cell types from a time course of differentiation across all germ layers, ranging from gastrulation cell states to progenitor and committed cell types. These data provide a diverse benchmarking reference point to compare against in vivo development and advance understanding of signalling regulation of differentiation, providing insights into protocol development, drug screening, and regenerative medicine applications.

Pluripotent Stem Cells↗

Diving Deeper Into Mechanisms of Acrylamide-Induced Toxicity: RNA Sequencing Reveals Transcriptomic Alteration and Retrotransposon Expression in Drosophila melanogaster.

Given the inevitability of human and animal exposure to acrylamide, there is increasing concern regarding its potential health risks. While a number of molecular mechanisms have been proposed, the complexity of acrylamide toxicological pathways and interactions remains incompletely characterized. In this study, we employed a transcriptomic approach to investigate the transcriptional responses of Drosophila melanogaster following exposure to acrylamide (100 mg/kg). Our analysis identified 634 differentially expressed genes (DEGs), with 362 upregulated and 272 downregulated. Functional analysis revealed these DEGs are enriched in pathways related to reproduction, detoxification, cellular and metabolic processes, signaling, synaptic formation and organization. Notably, acrylamide exposure upregulated the expression of tau and beta-amyloid protein precursor-like genes, both implicated in Alzheimer's disease pathology. An aversive memory test further demonstrated that acrylamide impaired the short-term memory of treated flies. Additionally, acrylamide-induced toxicity altered the expression of nine long terminal repeat retrotransposons, belonging to the gypsy and pao superfamilies. By exploring the potential role of transposable element activity in acrylamide-mediated toxicity, this study provides novel insights into the molecular mechanisms underlying its effects. Collectively, these findings offer a more comprehensive understanding of the mechanisms and pathways associated with the toxic action and detoxification of acrylamide in D. melanogaster.

Animals↗

Reliable amplification method for bacterial RNA.

DNA microarray technology has been increasingly applied for studies of clinical samples. Frequently, RNA probes from clinical samples are available in limited amounts. We describe a reliable amplification method for bacterial RNA. We verified this method on mycobacterial RNA applying mycobacterial genome-directed primers (mtGDPs). Glass slide-based oligoarrays were employed to assess the quality of the amplification method. We observed a relatively small bias in amplified RNA pool when compared to the unamplified one. Up to 1000-fold linear RNA amplification in a single amplification round was obtained. To our knowledge, this study describes the first amplification method for mycobacterial RNA.

Moloney murine leukemia virus↗

Benchmarking the CATMA microarray. A novel tool for Arabidopsis transcriptome analysis.

Transcript profiling is crucial to study biological systems, and various platforms have been implemented to survey mRNAs at the genome scale. We have assessed the performance of the CATMA microarray designed for Arabidopsis (Arabidopsis thaliana) transcriptome analysis and compared it with the Agilent and Affymetrix commercial platforms. The CATMA array consists of gene-specific sequence tags of 150 to 500 bp, the Agilent (Arabidopsis 2) array of 60mer oligonucleotides, and the Affymetrix gene chip (ATH1) of 25mer oligonucleotide sets. We have matched each probe repertoire with the Arabidopsis genome annotation (The Institute for Genomic Research release 5.0) and determined the correspondence between them. Array performance was analyzed by hybridization with labeled targets derived from eight RNA samples made of shoot total RNA spiked with a calibrated series of 14 control transcripts. CATMA arrays showed the largest dynamic range extending over three to four logs. Agilent and Affymetrix arrays displayed a narrower range, presumably because signal saturation occurred for transcripts at concentrations beyond 1,000 copies per cell. Sensitivity was comparable for all three platforms. For Affymetrix GeneChip data, the RMA software package outperformed Microarray Suite 5.0 for all investigated criteria, confirming that the information provided by the mismatch oligonucleotides has no added value. In addition, taking advantage of replicates in our dataset, we conducted a robust statistical analysis of the platform propensity to yield false positive and false negative differentially expressed genes, and all gave satisfactory results. The results establish the CATMA array as a mature alternative to the Affymetrix and Agilent platforms.

Arabidopsis↗

Substrate-dependent regulation of anaerobic degradation pathways for toluene and ethylbenzene in a denitrifying bacterium, strain EbN1.

Anaerobic biodegradation of toluene and ethylbenzene is of environmental concern and biochemical interest due to toxicity and novel reactions, respectively. The denitrifying strain EbN1 is unique in anaerobically degrading both alkylbenzenes via different pathways which converge at benzoyl coenzyme A. The organization of genes involved in both pathways was only recently determined for strain EbN1. In the present study, global expression analysis (DNA microarray and proteomics) indicated involvement of several thus-far-unknown proteins in the degradation of both alkylbenzenes. For example, orf68 and orf57, framing the ebd operon, are implicated in ethylbenzene degradation, and the ebA1932 and ebA1936 genes, located 7.2 kb upstream of the bbs operon, are implicated in toluene degradation. In addition, expression studies were now possible on the level of the complete pathways. Growth experiments demonstrated that degradative capacities for toluene and ethylbenzene could be simultaneously induced, regardless of the substrate used for adaptation. Regulation was studied at the RNA (real-time reverse transcription-PCR and DNA microarray) and protein (two-dimensional-difference gel electrophoresis) level by using cells adapted to anaerobic growth with benzoate, toluene, ethylbenzene, or a mixture of toluene and ethylbenzene. Expression of the two toluene-related operons (bss and bbs) was specifically induced in toluene-adapted cells. In contrast, genes involved in anaerobic ethylbenzene degradation were induced in ethylbenzene- and toluene-adapted cells, suggesting that toluene may act as a gratuitous inducer. In agreement with the predicted sequential regulation of the ethylbenzene pathway, Ebd proteins (encoding subunits of ethylbenzene dehydrogenase) were formed in ethylbenzene- but not in acetophenone-adapted cells, while Apc proteins (subunits of predicted acetophenone carboxylase) were formed under both conditions.

Anaerobiosis↗

Automated detection and quantitation of bacterial RNA by using electrical microarrays.

Low-density electrical 16S rRNA specific oligonucleotide microarrays and an automated analysis system have been developed for the identification and quantitation of pathogens. The pathogens are Escherichia coli, Pseudomonas aeruginosa, Enterococcus faecalis, Staphylococcus aureus, and Staphylococcus epidermidis, which are typically involved in urinary tract infections. Interdigitated gold array electrodes (IDA-electrodes), which have structures in the nanometer range, have been used for very sensitive analysis. Thiol-modified oligonucleotides are immobilized on the gold IDA as capture probes. They mediate the specific recognition of the target 16S rRNA by hybridization. Additionally three unlabeled oligonucleotides are hybridized in close proximity to the capturing site. They are supporting molecules, because they improve the RNA hybridization at the capturing site. A biotin labeled detector oligonucleotide is also allowed to hybridize to the captured RNA sequence. The biotin labels enable the binding of avidin alkaline phophatase conjugates. The phosphatase liberates the electrochemical mediator p-aminophenol from its electrically inactive phosphate derivative. The electrical signals were generated by amperometric redox cycling and detected by a unique multipotentiostat. The read out signals of the microarray are position specific current and change over time in proportion to the analyte concentration. If two additional biotins are introduced into the affinity binding complex via the supporting oligonucleotides, the sensitivity of the assays increase more than 60%. The limit of detection of Escherichia coli total RNA has been determined to be 0.5 ng/microL. The control of fluidics for variable assay formats as well as the multichannel electrical read out and data handling have all been fully automated. The fast and easy procedure does not require any amplification of the targeted nucleic acids by PCR.

Automation↗

Genome-directed primers for selective labeling of bacterial transcripts for DNA microarray analysis.

DNA microarrays have the ability to analyze the expression of thousands of the same set of genes under at least two different experimental conditions. However, DNA microarrays require substantial amounts of RNA to generate the probes, especially when bacterial RNA is used for hybridization (50 microg of bacterial total RNA contains approximately 2 microg of mRNA). We have developed a computer-based algorithm for prediction of the minimal number of primers to specifically anneal to all genes in a given genome. The algorithm predicts, for example, that 37 oligonucleotides should prime all genes in the Mycobacterium tuberculosis genome. We tested the usefulness of the genome-directed primers (GDPs) in comparison to random primers for gene expression profiling using DNA microarrays. Both types of primers were used to generate fluorescent-labeled probes and to hybridize to an array of 960 mycobacterial genes. Compared to random-primer probes, the GDP probes were more sensitive and more specific, especially when mammalian RNA samples were spiked with mycobacterial RNA. The GDPs were used for gene expression profiling of mycobacterial cultures grown to early log or stationary growth phases. This approach could be useful for accurate genome-wide expression analysis, especially for in vivo gene expression profiling, as well as directed amplification of sequenced genomes.

Algorithms↗

Gene expression levels in small specimens from patients detected using oligonucleotide arrays.

Large-scale gene expression profiling using microarray technology is often limited by the amount of tissue or cells available. A number of RNA amplification protocols have been published to overcome this problem. However, additional amplification steps can result in both a 3' bias and poor reproducibility for low abundance transcripts. We performed microarray experiments using HG-U133A GeneChip arrays to ascertain whether less than the recommended amount of RNA can be used, thus avoiding additional amplification steps. In a titration experiment, 2-10 microg of total RNA from a single cryopreserved patient specimen was used to prepare biotinylated cRNA, and the recommended standard amount of 15 microg of each preparation was used for hybridization. Statistical analysis using box plots, correlation coefficients, MvA plots, and concordance percentages revealed almost identical levels of gene expression, independent of the amount of RNA used for target preparation. Most importantly, there was no statistically significant difference when the concordance percentages for low abundance genes were compared, demonstrating that as little as 2 microg of total RNA is sufficient to perform GeneChip analysis.

Child↗

Identification and analysis of gene families from the duplicated genome of soybean using EST sequences.

BACKGROUND: Large scale gene analysis of most organisms is hampered by incomplete genomic sequences. In many organisms, such as soybean, the best source of sequence information is the existence of expressed sequence tag (EST) libraries. Soybean has a large (1115 Mbp) genome that has yet to be fully sequenced. However it does have the 6th largest EST collection comprised of ESTs from a variety of soybean genotypes. Many EST libraries were constructed from RNA extracted from various genetic backgrounds, thus gene identification from these sources is complicated by the existence of both gene and allele sequence differences. We used the ESTminer suite of programs to identify potential soybean gene transcripts from a single genetic background allowing us to observe functional classifications between gene families as well as structural differences between genes and gene paralogs within families. The identification of potential gene sequences (pHaps) from soybean allows us to begin to get a picture of the genomic history of the organism as well as begin to observe the evolutionary fates of gene copies in this highly duplicated genome. RESULTS: We identified approximately 45,000 potential gene sequences (pHaps) from EST sequences of Williams/Williams82, an inbred genotype of soybean (Glycine max L. Merr.) using a redundancy criterion to identify reproducible sequence differences between related genes within gene families. Analysis of these sequences revealed single base substitutions and single base indels are the most frequently observed form of sequence variation between genes within families in the dataset. Genomic sequencing of selected loci indicate that intron-like intervening sequences are numerous and are approximately 220 bp in length. Functional annotation of gene sequences indicate functional classifications are not randomly distributed among gene families containing few or many genes. CONCLUSION: The predominance of single nucleotide insertion/deletions and substitution events between genes within families (individual genes and gene paralogs) is consistent with a model of gene amplification followed by single base random mutational events expected under the classical model of duplicated gene evolution. Molecular functions of small and large gene families appear to be non-randomly distributed possibly indicating a difference in retention of duplicates or local expansion.

Evolution, Molecular↗

The role of the omega subunit of RNA polymerase in expression of the relA gene in Escherichia coli.

The rpoZ gene for the omega subunit of Escherichia coli RNA polymerase constitutes single operon with the spoT gene, which is responsible for the maintenance of stringent response under nutrient starvation conditions. To identify the physiological role of the omega subunit, we compared the gene expression profile of wild-type Escherichia coli with that of an rpoZ deleted strain by microarray analysis using an E. coli DNA chip. Here we report on a set of genes which show changes in expression profile following the removal of rpoZ. We have seen that relA, which is responsible for the synthesis of the stringent factor ppGpp and many ribosomal proteins, exhibited noticeable changes in mRNA levels and were therefore further analyzed for their expression using a GFP/RFP two-fluorescent protein promoter assay vector. In the absence of rpoZ, the promoter for the relA gene was severely impaired, but the promoters from the ribosomal protein genes were not affected as much. Taking these results together we propose that the omega subunit is involved in regulation of the relA gene, but induction of the stringently controlled genes in the absence of rpoZ is, at least in part, attributable to a decrease in ppGpp level.

Artificial Gene Fusion↗

Analysis of the Leishmania donovani transcriptome reveals an ordered progression of transient and permanent changes in gene expression during differentiation.

Leishmania donovani is an intracellular protozoan parasite that causes kala-azar in humans. During infection the extracellular insect forms (promastigotes) undergo rapid differentiation to intracellular amastigotes that proliferates in phagolysosomes of mammalian macrophages. We used microarray-based expression profiling to investigate the time-course of changes in RNA abundance during promastigote-to-amastigote differentiation in a host-free system that mimics this process. These studies revealed that several hundred genes underwent an ordered progression of transient or permanent up- and down-regulation during differentiation. Genes that were permanently up-regulated in amastigotes were enriched for transporters and surface proteins, but under-represented in genes involved in protein and other metabolism. Most of these changes occurred late in the differentiation process, when morphological differentiation was essentially complete. Down-regulated genes were over-represented in those involved in cell motility, growth and/or maintenance, and these changes generally occurred earlier in the process. Genes that were transiently up- or down-regulated during differentiation included those encoding heat shock proteins, ubiquitin hydrolases, RNA binding proteins, protein kinases, a protein phosphatase, and a histone deacetylase. These results suggest that changes in mRNA abundance may be important in signal transduction, as well as protein and mRNA turnover, during differentiation. In addition to these mRNA changes, other transcripts including one or more rRNAs and snoRNAs, and non-coding RNAs from several telomeres, also showed substantial changes in abundance during the differentiation process. This paper provides the first genome-scale quantitative analysis of gene expression during the transition from promastigotes to amastigotes and demonstrates the utility of the host-free differentiation system.

Animals↗

Evaluation of differential gene expression in Leishmania major Friedlin procyclics and metacyclics using DNA microarray analysis.

The various Leishmania species are flagellated protozoans, responsible for a wide spectrum of human diseases. The sequence of the L. major genome is nearing completion and a large proportion of the identified genes have yet to be ascribed functions. DNA microarrays containing PCR-amplified DNA from a random amplified genomic library of L. major Friedlin (LmjF) [Mol. Biochem. Parasitol. 113 (2001) 337] were hybridized with fluorescent probes made from L. major Friedlin RNA from five time-points during differentiation from procyclics to metacyclics. The data were normalized for background and probe intensity and the relative abundance of RNA for each spot was calculated. Almost 15% (1387/9282) of the DNAs showed statistically significant (P<0.01) changes in expression (1.1-5-fold) during the transition, with 1.16% (108) showing up-regulation at two or more time-points and 0.14% (13) showing down-regulation. Northern blot analyses of selected genes confirmed these results. These studies confirmed the stage-specific expression of several known genes, as well as identifying a number of novel genes that are up-regulated in either procyclics or metacyclics.

Animals↗

Label-free detection of 16S ribosomal RNA hybridization on reusable DNA arrays using surface plasmon resonance imaging.

In this paper, we describe the detection of bacterial cell-extracted 16S ribosomal RNA (rRNA) using an emerging technology, surface plasmon resonance (SPR) imaging of DNA arrays. Surface plasmon resonance enables detection of molecular interactions on surfaces in response to changes in the index of refraction, therefore eliminating the need for a fluorescent or radioactive label. A variation of the more common SPR techniques, SPR imaging enables detection from multiple probes in a reusable array format. The arrays developed here contain DNA probes (15-21 bases) designed to be complementary to 16S rRNA gene sequences of Escherichia coli and Bacillus subtilis as well as to a highly conserved sequence found in rRNAs from most members of the domain Bacteria. We report species-specific hybridization of cell-extracted total RNA and in vitro transcribed 16S rRNA to oligonucleotide probes on SPR arrays. We tested multiple probe sequences for each species, and found that success or failure of hybridization was dependent upon probe position in the 16S rRNA molecule. It was also determined that one of the probes intended to bind 16S rRNA also bound an unknown protein. The amount of binding to these probes was quantified with SPR imaging. A detection limit of 2 micro g ml-1 was determined for fragmented E. coli total cellular RNA under the experimental conditions used. These results indicate the feasibility of using SPR imaging for 16S rRNA identification and encourage further development of this method for direct detection of other RNA molecules.

DNA Probes↗

NextLongIso: a comprehensive Nextflow pipeline for multi-dimensional long-read RNA-seq analysis.

SUMMARY: Long-read RNA sequencing technologies, including Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT), enable direct characterization of full-length transcripts and transcriptome complexity. However, analysis of long-read RNA-seq data remains fragmented across multiple tools, limiting the ability to obtain a unified view of transcript structure, expression, and regulatory variation in long-read transcriptomes. We present NextLongIso, a scalable and reproducible Nextflow pipeline that enables coordinated analysis of multiple layers of transcript regulation. Rather than focusing solely on transcript reconstruction, NextLongIso integrates transcript discovery with downstream regulatory analyses to jointly characterize alternative splicing, isoform switching, transcript boundary dynamics (including alternative promoters and polyadenylation), and transposable element-associated transcription from both PacBio and ONT datasets. By eliminating complex cross-tool data harmonization, this unified framework facilitates the transition from transcript identification to functional interpretation of transcriptomic variation. AVAILABILITY AND IMPLEMENTATION: NextLongIso is implemented in Nextflow and is freely available at github: https://github.com/YidanSunResearchLab/nf-LongIso.git and Zenodo: https://doi.org/10.5281/zenodo.21049837.

Software↗

Identification of differentially expressed genes in keratoconus epithelium analyzed on microarrays.

PURPOSE: The present study was conducted to investigate differential gene expression in individual samples derived from fresh-frozen human keratoconus and normal corneal epithelium, using gene microarrays. METHODS: Total RNA was extracted (11 keratoconus and 8 normal samples), labeled, and hybridized to microarrays (GeneChip; Affymetrix, Inc., Santa Clara, CA). GeneChip data were validated by verifying the expression profiles of 10 genes by real-time PCR and by recalculation using dChip software (Wong Laboratory, Department of Biostatistics, Harvard School of Public Health, Boston, MA). Furthermore, 3 of the 10 encoded proteins were stained by immunohistochemistry. RESULTS: In comparison with normal cornea, the expression of 471 of the 5600 genes on the microarrays was changed in the keratoconus samples. This number was reduced to 47 with increased expression and 9 with decreased expression when more stringent selection parameters were applied. These genes are believed to be involved in keratoconus. Two of the candidate genes, lysyl oxidase and tissue inhibitor of metalloproteinase 3, are known to be involved in other eye diseases. Expression profiles were reproduced with the software dChip (Wong Laboratory) and real-time PCR. Increases in keratin 6 and 13 were also detected at the protein level. CONCLUSIONS: Keratoconus epithelium appears to be characterized by massive changes of the cytoskeleton, reduced extracellular matrix remodeling, altered transmembrane signaling, and modified cell-to-cell and cell-to-matrix interactions. Validation of gene expression with dChip analysis and real-time PCR indicates GeneChip to be a valid technique for investigation of epithelium from single dissected corneal samples. Association between alterations at RNA and protein levels was observed for some of the tested candidates.

Adult↗