High-throughput and high-accuracy single-cell RNA isoform analysis using PacBio circular consensus sequencing.
Ontology highlight
ABSTRACT:
SUBMITTER: Shi ZX
PROVIDER: S-EPMC10164132 | biostudies-literature | 2023 May
REPOSITORIES: biostudies-literature
Publications
High-throughput and high-accuracy single-cell RNA isoform analysis using PacBio circular consensus sequencing.
Shi Zhuo-Xing ZX Chen Zhi-Chao ZC Zhong Jia-Yong JY Hu Kun-Hua KH Zheng Ying-Feng YF Zheng Ying-Feng YF Chen Ying Y Xie Shang-Qian SQ Bo Xiao-Chen XC Luo Feng F Tang Chong C Xiao Chuan-Le CL Liu Yi-Zhi YZ
Nature communications 20230506 1
Although long-read single-cell RNA isoform sequencing (scISO-Seq) can reveal alternative RNA splicing in individual cells, it suffers from a low read throughput. Here, we introduce HIT-scISOseq, a method that removes most artifact cDNAs and concatenates multiple cDNAs for PacBio circular consensus sequencing (CCS) to achieve high-throughput and high-accuracy single-cell RNA isoform sequencing. HIT-scISOseq can yield >10 million high-accuracy long-reads in a single PacBio Sequel II SMRT Cell 8M. ...[more]
Similar Datasets
Project description:Long single-molecular sequencing technologies, such as PacBio circular consensus sequencing (CCS) and nanopore sequencing, are advantageous in detecting DNA 5-methylcytosine in CpGs (5mCpGs), especially in repetitive genomic regions. However, existing methods for detecting 5mCpGs using PacBio CCS are less accurate and robust. Here, we present ccsmeth, a deep-learning method to detect DNA 5mCpGs using CCS reads. We sequence polymerase-chain-reaction treated and M.SssI-methyltransferase treated DNA of one human sample using PacBio CCS for training ccsmeth. Using long (≥10 Kb) CCS reads, ccsmeth achieves 0.90 accuracy and 0.97 Area Under the Curve on 5mCpG detection at single-molecule resolution. At the genome-wide site level, ccsmeth achieves >0.90 correlations with bisulfite sequencing and nanopore sequencing using only 10× reads. Furthermore, we develop a Nextflow pipeline, ccsmethphase, to detect haplotype-aware methylation using CCS reads, and then sequence a Chinese family trio to validate it. ccsmeth and ccsmethphase can be robust and accurate tools for detecting DNA 5-methylcytosines.
| S-EPMC10329642 | biostudies-literature
Project description:BackgroundLiver cirrhosis is a well-known risk factor for hepatocellular carcinoma (HCC). However, some HCC cases can also originate from non-cirrhotic livers. The aim of this study was to identify key circular RNAs (circRNAs) associated with the tumorigenesis of non-cirrhotic liver disease.MethodsThe differently expressed circRNAs between non-cirrhotic and cirrhotic HCCs were assessed with use of high-throughput circRNAs sequencing and validated with quantitative reverse transcription polymerase chain reaction (qRT-PCR). Potential biological functions of these dysregulated circRNAs were predicted with use of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. A circRNA-miRNA-mRNA regulation network was constructed as achieved with use of miRanda software and visualized using Cytoscape software. Biological functions of the four most prominent dysregulated circRNAs identified were confirmed by in vitro experiments. Moreover, possible translations of these dysregulated circRNAs were also predicted.ResultsA total of 393 dysregulated circRNAs were identified between non-cirrhotic and cirrhotic HCC, including 213 that were significantly up-regulated and 180 significantly down-regulated circRNAs. Expression levels of the six most prominent dysregulated circRNAs were further validated using qRT-PCR. Many tumor related miRNAs were involved in the circRNA-miRNA-mRNA networks, including miR-182-5p, miR-561-3p, miR-125a-5p, miR-145, miR-23b-3p and miR-30e-3p, and downstream mRNAs of dysregulated circRNAs were significantly related with biological processes involved in the progression of tumors, including proliferation, migration, differentiation, and focal adhesion. Results from the in vitro experiments demonstrated that the most prominent dysregulated circRNAs exerted notable effects upon the proliferation and migration of HCC cells. Finally, we also identified 19 dysregulated circRNAs having potential for the coding of functional peptides.ConclusionThe results of this present study indicate that circRNAs may play important roles in tumorigenesis of non-cirrhotic HCC. Such findings provide some novel insights and pave the way for the development of future studies directed at investigating the initiation and treatment of HCC.
| S-EPMC9356431 | biostudies-literature
Project description:It is only recently that the abundant presence of circular RNAs (circRNAs) in all kingdoms of Life, including the hyperthermophilic archaeon Pyrococcus abyssi, has emerged. This led us to investigate the physiologic significance of a previously observed weak intramolecular ligation activity of Pab1020 RNA ligase. Here we demonstrate that this enzyme, despite sharing significant sequence similarity with DNA ligases, is indeed an RNA-specific polynucleotide ligase efficiently acting on physiologically significant substrates. Using a combination of RNA immunoprecipitation assays and RNA-seq, our genome-wide studies revealed 133 individual circRNA loci in P. abyssi. The large majority of these loci interacted with Pab1020 in cells and circularization of selected C/D Box and 5S rRNA transcripts was confirmed biochemically. Altogether these studies revealed that Pab1020 is required for RNA circularization. Our results further suggest the functional speciation of an ancestral NTase domain and/or DNA ligase toward RNA ligase activity and prompt for further characterization of the widespread functions of circular RNAs in prokaryotes. Detailed insight into the cellular substrates of Pab1020 may facilitate the development of new biotechnological applications e.g. in ligation of preadenylated adaptors to RNA molecules.
| S-EPMC5680704 | biostudies-literature
Project description:MotivationToday, the base code of DNA is mostly determined through sequencing by synthesis as provided by the Illumina sequencers. Although highly accurate, resulting reads are short, making their analyses challenging. Recently, a new technology, single molecule real-time (SMRT) sequencing, was developed that could address these challenges, as it generates reads of several thousand bases. But, their broad application has been hampered by a high error rate. Therefore, hybrid approaches that use high-quality short reads to correct erroneous SMRT long reads have been developed. Still, current implementations have great demands on hardware, work only in well-defined computing infrastructures and reject a substantial amount of reads. This limits their usability considerably, especially in the case of large sequencing projects.ResultsHere we present proovread, a hybrid correction pipeline for SMRT reads, which can be flexibly adapted on existing hardware and infrastructure from a laptop to a high-performance computing cluster. On genomic and transcriptomic test cases covering Escherichia coli, Arabidopsis thaliana and human, proovread achieved accuracies up to 99.9% and outperformed the existing hybrid correction programs. Furthermore, proovread-corrected sequences were longer and the throughput was higher. Thus, proovread combines the most accurate correction results with an excellent adaptability to the available hardware. It will therefore increase the applicability and value of SMRT sequencing.Availability and implementationproovread is available at the following URL: http://proovread.bioapps.biozentrum.uni-wuerzburg.de.
| S-EPMC4609002 | biostudies-literature
Project description:Single-cell transcriptomic studies that require intracellular protein staining, rare cell sorting, or inactivation of infectious pathogens are severely limited. This is because current high-throughput single-cell RNA sequencing methods are either incompatible with or necessitate laborious sample preprocessing for paraformaldehyde treatment, a common tissue and cell fixation and preservation technique. Here we present FD-seq (Fixed Droplet RNA sequencing), a high-throughput method for droplet-based RNA sequencing of paraformaldehyde-fixed, permeabilized and sorted single cells. We show that FD-seq preserves the RNA integrity and relative gene expression levels after fixation and permeabilization. Furthermore, FD-seq can detect a higher number of genes and transcripts than methanol fixation. We first apply FD-seq to analyze a rare subpopulation of cells supporting lytic reactivation of the human tumor virus KSHV, and identify TMEM119 as a potential host factor that mediates viral reactivation. Second, we find that infection with the human betacoronavirus OC43 leads to upregulation of pro-inflammatory pathways in cells that are exposed to the virus but fail to express high levels of viral genes. FD-seq thus enables integrating phenotypic with transcriptomic information in rare cell subpopulations, and preserving and inactivating pathogenic samples.
| S-EPMC8463713 | biostudies-literature
Project description:BACKGROUND: Molecular characterization of highly diverse gene families can be time consuming, expensive, and difficult, especially when considering the potential for relatively large numbers of paralogs and/or pseudogenes. Here we investigate the utility of Pacific Biosciences single molecule real-time (SMRT) circular consensus sequencing (CCS) as an alternative to traditional cloning and Sanger sequencing PCR amplicons for gene family characterization. We target vomeronasal gene receptors, one of the most diverse gene families in mammals, with the goal of better understanding intra-specific V1R diversity of the gray mouse lemur (Microcebus murinus). Our study compares intragenomic variation for two V1R subfamilies found in the mouse lemur. Specifically, we compare gene copy variation within and between two individuals of M. murinus as characterized by different methods for nucleotide sequencing. By including the same individual animal from which the M. murinus draft genome was derived, we are able to cross-validate gene copy estimates from Sanger sequencing versus CCS methods. RESULTS: We generated 34,088 high quality circular consensus sequences of two diverse V1R subfamilies (here referred to as V1RI and V1RIX) from two individuals of Microcebus murinus. Using a minimum threshold of 7× coverage, we recovered approximately 90% of V1RI sequences previously identified in the draft M. murinus genome (59% being identical at all nucleotide positions). When low coverage sequences were considered (i.e. < 7× coverage) 100% of V1RI sequences identified in the draft genome were recovered. At least 13 putatively novel V1R loci were also identified using CCS technology. CONCLUSIONS: Recent upgrades to the Pacific Biosciences RS instrument have improved the CCS technology and offer an alternative to traditional sequencing approaches. Our results suggest that the Microcebus murinus V1R repertoire has been underestimated in the draft genome. In addition to providing an improved understanding of V1R diversity in the mouse lemur, this study demonstrates the utility of CCS technology for characterizing complex regions of the genome. We anticipate that long-read sequencing technologies such as PacBio SMRT will allow for the assembly of multigene family clusters and serve to more accurately characterize patterns of gene copy variation in large gene families, thus revealing novel micro-evolutionary patterns within non-model organisms.
| S-EPMC4152597 | biostudies-literature
Project description:BackgroundLarge insert paired-end sequencing technologies are important tools for assembling genomes, delineating associated breakpoints and detecting structural rearrangements. To facilitate the comprehensive detection of inter- and intra-chromosomal structural rearrangements or variants (SVs) and complex genome assembly with long repeats and segmental duplications, we developed a new method based on single-molecule real-time synthesis sequencing technology for generating long paired-end sequences of large insert DNA libraries.ResultsA Fosmid vector, pHZAUFOS3, was developed with the following new features: (1) two 18-bp non-palindromic I-SceI sites flank the cloning site, and another two sites are present in the skeleton of the vector, allowing long DNA inserts (and the long paired-ends in this paper) to be recovered as single fragments and the vector (~ 8 kb) to be fragmented into 2-3 kb fragments by I-SceI digestion and therefore was effectively removed from the long paired-ends (5-10 kb); (2) the chloramphenicol (Cm) resistance gene and replicon (oriV), necessary for colony growth, are located near the two sides of the cloning site, helping to increase the proportion of the paired-end fragments to single-end fragments in the paired-end libraries. Paired-end libraries were constructed by ligating the size-selected, mechanically sheared pooled Fosmid DNA fragments to the Ampicillin (Amp) resistance gene fragment and screening the colonies with Cm and Amp. We tested this method on yeast and Setaria italica Yugu1. Fosmid-size paired-ends with an average length longer than 2 kb for each end were generated. The N50 scaffold lengths of the de novo assemblies of the yeast and S. italica Yugu1 genomes were significantly improved. Five large and five small structural rearrangements or assembly errors spanning tens of bp to tens of kb were identified in S. italica Yugu1 including deletions, inversions, duplications and translocations.ConclusionsWe developed a new method for long paired-end sequencing of large insert libraries, which can efficiently improve the quality of de novo genome assembly and identify large and small structural rearrangements or assembly errors.
| S-EPMC6878638 | biostudies-literature
Project description:ObjectivesAberrantly expressed circular RNAs (circRNAs) have been detected in many types of tumors. Hence, they are currently investigated as candidate biomarkers for diagnostic and potential targets for therapy in cancers. The objective of this study was to assess the expression profile of circRNA in lung adenocarcinoma (LUAD).MethodsThis study included 14 pairs of postoperative lung adenocarcinoma specimens, including cancer tissues and matched adjacent tissues. Second-generation sequencing was applied to the specimens to determine the circRNA expression in them among the 5242 distinct circRNAs detected.ResultsWe identified a total of 18 significantly dysregulated circRNAs in the LUAD tissues: upregulation in four and downregulation in 14. ROC (The receiver operating characteristic curve) further suggested that hsa_circ_0120106, has_circ_0007342, has_circ_0005937, and circRNA_0000826 could potentially be used as biomarkers in the diagnosis of LUAD. Furthermore, study of the circRNA-microRNA (miRNA)-messenger RNA (mRNA) revealed interactions between the 18 dysregulated circRNA and several cancer-related miRNAs. Finally, a further Kyoto Encyclopedia of Genes and Genomes analysis showed that the cell cycle phase transition, p53 signaling pathway, AMP-activated protein kinase (AMPK) relative signaling pathway, and so on were key putative pathways in the process of LUAD.ConclusionsThese findings demonstrated the correlation between abnormality in circRNA expression and LUAD, which lays the foundation of making CircRNAs candidate biomarkers in the diagnosis of LUAD.
| S-EPMC10234779 | biostudies-literature
Project description:Bacteria colonize almost all parts of the human body and can differ significantly. However, the population level transcriptomics measurements can only describe the average bacteria population behaviors, ignoring the heterogeneity among bacteria. Here, we report a droplet-based high-throughput single-microbe RNA-seq assay (smRandom-seq), using random primers for in situ cDNA generation, droplets for single-microbe barcoding, and CRISPR-based rRNA depletion for mRNA enrichment. smRandom-seq showed a high species specificity (99%), a minor doublet rate (1.6%), a reduced rRNA percentage (32%), and a sensitive gene detection (a median of ~1000 genes per single E. coli). Furthermore, smRandom-seq successfully captured transcriptome changes of thousands of individual E. coli and discovered a few antibiotic resistant subpopulations displaying distinct gene expression patterns of SOS response and metabolic pathways in E. coli population upon antibiotic stress. smRandom-seq provides a high-throughput single-microbe transcriptome profiling tool that will facilitate future discoveries in microbial resistance, persistence, microbe-host interaction, and microbiome research.
| S-EPMC10447461 | biostudies-literature
Project description:Most methods for single-cell transcriptome sequencing amplify the termini of polyadenylated transcripts, capturing only a small fraction of the total cellular transcriptome. This precludes the detection of many long non-coding, short non-coding and non-polyadenylated protein-coding transcripts and hinders alternative splicing analysis. We, therefore, developed VASA-seq to detect the total transcriptome in single cells, which is enabled by fragmenting and tailing all RNA molecules subsequent to cell lysis. The method is compatible with both plate-based formats and droplet microfluidics. We applied VASA-seq to more than 30,000 single cells in the developing mouse embryo during gastrulation and early organogenesis. Analyzing the dynamics of the total single-cell transcriptome, we discovered cell type markers, many based on non-coding RNA, and performed in vivo cell cycle analysis via detection of non-polyadenylated histone genes. RNA velocity characterization was improved, accurately retracing blood maturation trajectories. Moreover, our VASA-seq data provide a comprehensive analysis of alternative splicing during mammalian development, which highlighted substantial rearrangements during blood development and heart morphogenesis.
| S-EPMC9750877 | biostudies-literature