Chimeric Reads Research Articles

A Genomic Language Model for Chimera Artifact Detection in Nanopore Direct RNA Sequencing.

Chimera artifacts in nanopore direct RNA sequencing (dRNA-seq) can significantly distort transcriptome analyses, yet their detection and removal remain challenging due to limitations in existing basecalling models. We present Deep-Chopper, a genomic language model that precisely identifies and removes adapter sequences from base-called dRNA-seq long reads at single-base resolution, operating independently of raw signal or alignment information to effectively eliminate chimeric read artifacts. By removing these artifacts, DeepChopper substantially improves the accuracy of critical downstream analyses, such as transcript annotation and gene fusion detection, thereby enhancing the reliability and utility of nanopore dRNA-seq for transcriptomics research.

RNA architecture of porcine deltacoronavirus genome inside virions detected by vRIC-seq

Porcine deltacoronavirus (PDCoV) is a newly emerging and special delta coronavirus, which infect mammals such as pigs, cattle and humans, as well as chickens and birds. Exploring RNA structures in the viral genome benefits the understanding of the role of RNA in the lifecycle of viruses. In this study, vRIC-seq is employed to analyze the RNA-RNA interaction in the whole genome structure of PDCoV in virions. About 12.87 and 13.52 million paired reads are obtained in two biological replicates, respectively, with 17.9% and 14.8% of them are identified as valid chimeric reads. These are employed to predict the RNA secondary structure, which is compact and highly structured. A twisted-cyclized conformation is observed in the RNA-RNA interaction map of PDCoV for the first time. 77 multi-way junctions are evenly distributed in the PDCoV genome. Our work provides fundamental structural insights that are essential for understanding the genomic structure and function, genetic evolution, and packaging characteristics of PDCoV.

Global analysis of the RNA-RNA interactome in Acinetobacter baumannii AB5075 uncovers a small regulatory RNA repressing the virulence-related outer membrane protein CarO.

Acinetobacter baumannii is an opportunistic Gram-negative pathogen that infects critically ill patients. The emergence of antimicrobial resistant A. baumannii has exacerbated the need to characterize environmental adaptation, antibiotic resistance and pathogenicity and their genetic regulators to inform intervention strategies. Critical to adaptation to changing environments in bacteria are small regulatory RNAs (sRNAs), however, the role that sRNAs play in the biology of A. baumannii is poorly understood. To assess the regulatory function of sRNAs and to uncover their RNA interaction partners, we employed an RNA proximity ligation and sequencing method (Hi-GRIL-seq) in three different environmental conditions. Forty sRNAs were ligated to sRNA-RNA chimeric sequencing reads, suggesting that sRNA-mediated gene regulation is pervasive in A. baumannii. In-depth characterization uncovered the sRNA Aar to be a post-transcriptional regulator of four mRNA targets including the transcript encoding outer membrane protein CarO. Aar initiates base-pairing with these mRNAs using a conserved seed region of nine nucleotides, sequestering the ribosome binding sites and inhibiting translation. Aar is differentially expressed in multiple stress conditions suggesting a role in fine-tuning translation of the Aar-target molecules. Our study provides mechanistic insights into sRNA-mediated gene regulation in A. baumannii and represents a valuable resource for future RNA-centric research endeavours.

Inference of genomic lesions from single-cell RNA-seq in myeloma improves functional intraclonal and interclonal analysis

AbstractSmoldering multiple myeloma (SMM) is an asymptomatic plasma cell (PC) neoplasm that may evolve with variable frequency into multiple myeloma (MM). SMM is initiated by chromosomal translocations involving the immunoglobulin heavy-chain locus or by hyperdiploidy and evolves through acquisition of additional genetic lesions. In this scenario, we aimed at establishing a reliable analysis pipeline to infer genomic lesions from transcriptomic analysis, by combining single-cell RNA sequencing (scRNA-seq) with B-cell receptor sequencing and copy number abnormality (CNA) analysis to identify clonal PCs at the genetic level along their specific transcriptional landscape. We profiled 20 465 bone marrow PCs derived from 5 patients with SMM/MM and unbiasedly identified clonal and polyclonal PCs. Hyperdiploidy, t(11;14), and t(6;14) were identified at the scRNA level by analysis of chimeric reads. Subclone functional analysis was improved by combining transcriptome with CNA analysis. As examples, we illustrate the different functional properties of a light-chain escape subclone in SMM and of different B-cell and PC subclones in a patient affected by Wäldenstrom macroglobulinemia and SMM. Overall, our data provide a proof of principle for inference of clinically relevant genotypic data from scRNA-seq, which in turn will refine functional annotation of the clonal architecture of PC dyscrasias.

Identification of transcriptionally-active human papillomavirus integrants through nanopore sequencing reveals viable targets for gene therapy against cervical cancer.

Integration of the human papillomavirus (HPV) genome into the cellular genome is a key event that leads to constitutive expression of viral oncoprotein E6/E7 and drives the progression of cervical cancer. However, HPV integration patterns differ on a case-by-case basis among related malignancies. Next-generation sequencing technologies still face challenges for interrogating HPV integration sites. In this study, utilizing Nanopore long-read sequencing, we identified 452 and 108 potential integration sites from the cervical cancer cell lines (CaSki and HeLa) and five tissue samples, respectively. Based on long Nanopore chimeric reads, we were able to analyze the methylation status of the HPV long control region (LCR), which controls oncogene E6/E7 expression, and to identify transcriptionally-active integrants among the numerous integrants. As a proof of concept, we identified an active HPV integrant in between RUNX2 and CLIC5 on chromosome 6 in the CaSki cell line, which was supported by ATAC-seq, H3K27Ac ChIP-seq, and RNA-seq analysis. Knockout of the active HPV integrant, by the CRISPR/Cas9 system, dramatically crippled cell proliferation and induced cell senescence. In conclusion, identifying transcriptionally-active HPV integrants with Nanopore sequencing can provide viable targets for gene therapy against HPV-associated cancers.

Variations in the persistence of 5′-end genomic and subgenomic SARS-CoV-2 RNAs in wastewater from aircraft, airports and wastewater treatment plants

Wastewater sequencing has become a powerful supplement to clinical testing in monitoring SARS-CoV-2 infections in the post-COVID-19 pandemic era. While its applications in measuring the viral burden and main circulating lineages in the community have proved their efficacy, the variations in sequencing quality and coverage across the different regions of the SARS-CoV-2 genome are not well understood. Furthermore, it is unclear how different sample origins, viral extraction and concentration methods and environmental factors impact the reads sequenced from wastewater. Using high-coverage, amplicon-based, paired-end read sequencing of viral RNA extracted from wastewater collected directly from aircraft, pooled from different aircraft and airport buildings or from regular wastewater plants, we assessed the genome coverage across the sample groups with a focus on the 5′-end region covering the leader sequence and investigated whether it was possible to detect subgenomic RNA from viral material recovered from wastewater. We identified distinct patterns in the persistence of the different genomic regions across the different types of wastewaters and the existence of chimeric reads mapping to non-amplified regions. Our findings suggest that preservation of the 5′-end of the genome and the ability to detect subgenomic RNA reads, though highly susceptible to environment and sample processing conditions, may be indicative of the quality and amount of the viral RNA present in wastewater.

Spatial transcriptomics reveals a low extent of transcriptionally active hepatitis B virus integration in patients with HBsAg loss

ObjectiveHepatitis B virus (HBV) can integrate into the chromosomes of infected hepatocytes, contributing to the production of hepatitis B surface antigen (HBsAg) and to hepatocarcinogenesis. In this study, we aimed...

Characterization and mitigation of artifacts derived from NGS library preparation due to structure-specific sequences in the human genome

BackgroundHybridization capture-based targeted next generation sequencing (NGS) is gaining importance in routine cancer clinical practice. DNA library preparation is a fundamental step to produce high-quality sequencing data. Numerous unexpected, low variant allele frequency calls were observed in libraries using sonication fragmentation and enzymatic fragmentation. In this study, we investigated the characteristics of the artifact reads induced by sonication and enzymatic fragmentation. We also developed a bioinformatic algorithm to filter these sequencing errors.ResultsWe used pairwise comparisons of somatic single nucleotide variants (SNVs) and insertions and deletions (indels) of the same tumor DNA samples prepared using both ultrasonic and enzymatic fragmentation protocols. Our analysis revealed that the number of artifact variants was significantly greater in the samples generated using enzymatic fragmentation than using sonication. Most of the artifacts derived from the sonication-treated libraries were chimeric artifact reads containing both cis- and trans-inverted repeat sequences of the genomic DNA. In contrast, chimeric artifact reads of endonuclease-treated libraries contained palindromic sequences with mismatched bases. Based on these distinctive features, we proposed a mechanistic hypothesis model, PDSM (pairing of partial single strands derived from a similar molecule), by which these sequencing errors derive from ultrasonication and enzymatic fragmentation library preparation. We developed a bioinformatic algorithm to generate a custom mutation “blacklist” in the BED region to reduce errors in downstream analyses.ConclusionsWe first proposed a mechanistic hypothesis model (PDSM) of sequencing errors caused by specific structures of inverted repeat sequences and palindromic sequences in the natural genome. This new hypothesis predicts the existence of chimeric reads that could not be explained by previous models, and provides a new direction for further improving NGS analysis accuracy. A bioinformatic algorithm, ArtifactsFinder, was developed and used to reduce the sequencing errors in libraries produced using sonication and enzymatic fragmentation.

Enrichment Reveals Extensive Integration of Hepatitis B Virus DNA in Hepatitis Delta Virus-Infected Patients.

Hepatitis B virus (HBV) DNA may become integrated into the human genome of infected human hepatocytes. Expression of integrations can produce the surface antigen (HBsAg) that is required for synthesis of hepatitis D virus (HDV) particles and the abundant subviral particles in the blood of HBV- and HDV-infected subjects. Knowledge about the extent and variation of HBV integrations and impact on chronic HDV is still limited. We investigated 50 pieces of liver explant tissue from 5 patients with hepatitis D-induced cirrhosis, using a deep-sequencing strategy targeting HBV RNA. We found that integrations were abundant and highly expressed, with large variation in the number of integration-derived (HBV/human chimeric) reads, both between and within patients. The median number of unique integrations for each patient correlated with serum levels of HBsAg. However, most of the HBV reads represented a few predominant integrations. The results suggest that HBV DNA integrates in a large proportion of hepatocytes, and that the HBsAg output from these integrations vary >100-fold depending on clone size and expression rate. A small proportion of the integrations seems to determine the serum levels of HBsAg and HDV RNA in HBV/HDV coinfected patients with liver cirrhosis.

Analysis of chimeric reads characterises the diverse targetome of AGO2-mediated regulation

Argonaute proteins are instrumental in regulating RNA stability and translation. AGO2, the major mammalian Argonaute protein, is known to primarily associate with microRNAs, a family of small RNA ‘guide’ sequences, and identifies its targets primarily via a ‘seed’ mediated partial complementarity process. Despite numerous studies, a definitive experimental dataset of AGO2 ‘guide’–’target’ interactions remains elusive. Our study employs two experimental methods—AGO2 CLASH and AGO2 eCLIP, to generate thousands of AGO2 target sites verified by chimeric reads. These chimeric reads contain both the AGO2 loaded small RNA ‘guide’ and the target sequence, providing a robust resource for modeling AGO2 binding preferences. Our novel analysis pipeline reveals thousands of AGO2 target sites driven by microRNAs and a significant number of AGO2 ‘guides’ derived from fragments of other small RNAs such as tRNAs, YRNAs, snoRNAs, rRNAs, and more. We utilize convolutional neural networks to train machine learning models that accurately predict the binding potential for each ‘guide’ class and experimentally validate several interactions. In conclusion, our comprehensive analysis of the AGO2 targetome broadens our understanding of its ‘guide’ repertoire and potential function in development and disease. Moreover, we offer practical bioinformatic tools for future experiments and the prediction of AGO2 targets. All data and code from this study are freely available at https://github.com/ML-Bioinfo-CEITEC/HybriDetector/.

Heterozygosity alters Msh5 binding to meiotic chromosomes in the baker's yeast.

Meiotic crossovers are initiated from programmed DNA double-strand breaks. The Msh4-Msh5 heterodimer is an evolutionarily conserved mismatch repair-related protein complex that promotes meiotic crossovers by stabilizing strand invasion intermediates and joint molecule structures such as Holliday junctions. In vivo studies using homozygous strains of the baker's yeast Saccharomyces cerevisiae (SK1) show that the Msh4-Msh5 complex associates with double-strand break hotspots, chromosome axes, and centromeres. Many organisms have heterozygous genomes that can affect the stability of strand invasion intermediates through heteroduplex rejection of mismatch-containing sequences. To examine Msh4-Msh5 function in a heterozygous context, we performed chromatin immunoprecipitation and sequencing (ChIP-seq) analysis in a rapidly sporulating hybrid S. cerevisiae strain (S288c-sp/YJM789, containing sporulation-enhancing QTLs from SK1), using SNP information to distinguish reads from homologous chromosomes. Overall, Msh5 localization in this hybrid strain was similar to that determined in the homozygous strain (SK1). However, relative Msh5 levels were reduced in regions of high heterozygosity, suggesting that high mismatch densities reduce levels of recombination intermediates to which Msh4-Msh5 binds. Msh5 peaks were also wider in the hybrid background compared to the homozygous strain (SK1). We determined regions containing heteroduplex DNA by detecting chimeric sequence reads with SNPs from both parents. Msh5-bound double-strand break hotspots overlap with regions that have chimeric DNA, consistent with Msh5 binding to heteroduplex-containing recombination intermediates.

A streamlined pipeline based on HmmUFOtu for microbial community profiling using 16S rRNA amplicon sequencing.

Microbial community profiling using 16S rRNA amplicon sequencing allows for taxonomic characterization of diverse microorganisms. While amplicon sequence variant (ASV) methods are increasingly favored for their fine-grained resolution of sequence variants, they often discard substantial portions of sequencing reads during quality control, particularly in datasets with large number samples. We present a streamlined pipeline that integrates FastP for read trimming, HmmUFOtu for operational taxonomic units (OTU) clustering, Vsearch for chimera checking, and Kraken2 for taxonomic assignment. To assess the pipeline's performance, we reprocessed two published stool datasets of normal Korean populations: one with 890 and the other with 1,462 independent samples. In the first dataset, HmmUFOtu retained 93.2% of over 104 million read pairs after quality trimming, discarding chimeric or unclassifiable reads, while DADA2, a commonly used ASV method, retained only 44.6% of the reads. Nonetheless, both methods yielded qualitatively similar β-diversity plots. For the second dataset, HmmUFOtu retained 89.2% of read pairs, while DADA2 retained a mere 18.4% of the reads. HmmUFOtu, being a closed-reference clustering method, facilitates merging separately processed datasets, with shared OTUs between the two datasets exhibiting a correlation coefficient of 0.92 in total abundance (log scale). While the first two dimensions of the β-diversity plot exhibited a cohesive mixture of the two datasets, the third dimension revealed the presence of a batch effect. Our comparative evaluation of ASV and OTU methods within this streamlined pipeline provides valuable insights into their performance when processing large-scale microbial 16S rRNA amplicon sequencing data. The strengths of HmmUFOtu and its potential for dataset merging are highlighted.

RegCloser: a robust regression approach to closing genome gaps

BackgroundClosing gaps in draft genomes leads to more complete and continuous genome assemblies. The ubiquitous genomic repeats are challenges to the existing gap-closing methods, based on either the k-mer representation by the de Bruijn graph or the overlap-layout-consensus paradigm. Besides, chimeric reads will cause erroneous k-mers in the former and false overlaps of reads in the latter.ResultsWe propose a novel local assembly approach to gap closing, called RegCloser. It represents read coordinates and their overlaps respectively by parameters and observations in a linear regression model. The optimal overlap is searched only in the restricted range consistent with insert sizes. Under this linear regression framework, the local DNA assembly becomes a robust parameter estimation problem. We solved the problem by a customized robust regression procedure that resists the influence of false overlaps by optimizing a convex global Huber loss function. The global optimum is obtained by iteratively solving the sparse system of linear equations. On both simulated and real datasets, RegCloser outperformed other popular methods in accurately resolving the copy number of tandem repeats, and achieved superior completeness and contiguity. Applying RegCloser to a plateau zokor draft genome that had been improved by long reads further increased contig N50 to 3-fold long. We also tested the robust regression approach on layout generation of long reads.ConclusionsRegCloser is a competitive gap-closing tool. The software is available at https://github.com/csh3/RegCloser. The robust regression approach has a prospect to be incorporated into the layout module of long read assemblers.

Long-read Pore-C shows the 3D structure of the cattle genome

Context Recent advances in molecular technology have allowed us to examine the cattle genome with an accuracy never before possible. Genetic variations, both small and large, as well as the transcriptional landscape of the bovine genome, have both been explored in many studies. However, the topological configuration of the genome has not been extensively investigated, largely due to the cost of the assays required. Such assays can both identify topologically associated domains and be used for genome scaffolding. Aims This study aimed to implement a chromatin conformation capture together with long-read nanopore sequencing (Pore-C) pipeline for scaffolding a draft assembly and identifying topologically associating domains (TADs) of a Bos indicus Brahman cow. Methods Genomic DNA from a liver sample was first cross-linked to proteins, preserving the spatial proximity of loci. Restriction digestion and proximity ligation were then used to join cross-linked fragments, followed by nucleic isolation. The Pore-C DNA extracts were then prepped and sequenced on a PromethION device. Two genome assemblies were used to analyse the data, namely, one generated from sequencing of the same Brahman cow, and the other is the ARS-UCD1.2 Bos taurus assembly. The Pore-C snakemake pipeline was used to map, assign bins and scaffold the draft and current annotated bovine assemblies. The contact matrices were then used to identify TADs. Key results The study scaffolded a chromosome-level Bos indicus assembly representing 30 chromosomes. The scaffolded assembly showed a total of 215 contigs (2.6 Gbp) with N50 of 44.8 Mb. The maximum contig length was 156.8 Mb. The GC content of the scaffold assembly is 41 ± 0.02%. Over 50% of mapped chimeric reads identified for both assemblies had three or more contacts. This is the first experimental study to identify TADs in bovine species. In total, 3036 and 3094 TADs across 30 chromosomes were identified for input Brahman and ARS-UCD1.2 assemblies respectively. Conclusions The Pore-C pipeline presented herein will be a valuable approach to scaffold draft assemblies for agricultural species and understand the chromatin structure at different scales. Implications The Pore-C approach will open a new era of 3D genome-organisation studies across agriculture species.

A draft chromosome-scale genome assembly of a commercial sugarcane

Sugarcane accounts for a large portion of the worlds sugar production. Modern commercial cultivars are complex hybrids of S. officinarum, S. spontaneum, and several other Saccharum species, resulting in an auto-allopolyploid with 8–12 copies of each chromosome. The current genome assembly gold standard is to generate a long read assembly followed by chromatin conformation capture sequencing to scaffold. We used the PacBio RSII and chromatin conformation capture sequencing to sequence and assemble the genome of a South East Asian commercial sugarcane cultivar, known as Khon Kaen 3. The Khon Kaen 3 genome assembled into 104,477 contigs totalling 7 Gb, which scaffolded into 56 pseudochromosomes containing 5.2 Gb of sequence. Genome annotation produced 242,406 genes from 30,927 orthogroups. Aligning the Khon Kaen 3 genome sequence to S. officinarum and S. spontaneum revealed a high level of apparent recombination, indicating a chimeric assembly. This assembly error is explained by high nucleotide identity between S. officinarum and S. spontaneum, where 91.8% of S. spontaneum aligns to S. officinarum at 94% identity. Thus, the subgenomes of commercial sugarcane are so similar that using short reads to correct long PacBio reads produced chimeric long reads. Future attempts to sequence sugarcane must take this information into account.

Identification of Global miRNA Targetomes within Human Hematopoietic Progenitors and Acute Myeloid Leukemia Uncovers Unique Molecular Networks Centered on Mitochondria Function in Leukemic Stem Cells

Previous studies comparing transcriptome and proteome profiles of the cells that constitute both the normal and malignant hematopoietic hierarchies have revealed substantial uncoupling of transcript and protein expression levels. Despite the widely established role of miRNAs in regulating many cellular functions, the precise contribution of miRNA-mediated post-transcriptional regulation linking gene expression and proteome signatures of hematopoietic stem cells (HSC) and leukemic stem cells (LSC) remains largely unknown. Previous efforts to characterize miRNA targets have largely relied on computational prediction algorithms, which suffer from high rates of false-positive and negative target predictions and ultimately preclude discovery of non-canonical miRNA-target interactions. In this study, we perform Argonaute 2 (AGO2) crosslinking and immunoprecipitation (CLIP)-seq to experimentally map endogenous miRNA-target interactions in human CD34+ hematopoietic stem and progenitor cells (HSPC) from cord blood and six primary AML patient samples in an unbiased, transcriptome-wide manner (Manakov et al. 2022). Using this approach, we identify hundreds of miRNAs directly ligated to thousands of unique targets from miRNA-target chimeric sequencing reads. We find limited overlap (~30%) between miRNA targets captured by CLIP-seq and computationally predicted targets, underscoring the importance of functional validation of miRNA targets in specific cell types. Surprisingly, we find remarkable similarity between abundant miRNAs from HSPC and AML samples, including miR-142 and miRNAs encoded by the miR-17-92 cluster. However, the targets and the gene regulatory networks controlled by these highly abundant miRNAs are vastly different between HSPC and AML. In HSPC, we find that targets of abundant miRNAs are enriched in cell cycle regulation, histone modification and megakaryocyte/myeloid differentiation. In contrast, targets of the same miRNAs are enriched in protein localization/translation, neutrophil function, regulation of apoptosis and electron transport chain (ETC) pathways in AML. Collectively, we find that the composition of the abundant miRNAs is similar in normal HSPC and AML tumors; however they control divergent gene regulatory networks. To identify miRNA-target networks that play functional roles in LSC, we performed small RNA sequencing in 23 AML patient samples where each sorted fraction was assessed for LSC content by xenotransplantation. We derived an LSC-specific miRNA signature consisting of 5 miRNAs (LSC miR-5 score) predictive of overall survival and capable of stratifying AML patients independent of mRNA expression and other established indicators of outcome. Comparison of CLIP-seq targets of the LSC-specific miRNAs from HSPC and AML patient blasts revealed the existence of AML-specific gene regulatory networks centred on oxidative phosphorylation, ETC and cell-matrix adhesion. We further interrogated the functional significance of one miRNA from the LSC miR-5 score, miR-let-7c. We find that miR-let-7c is highly expressed in LSC-enriched populations of two AML in vitro models that recapitulate the leukemic hierarchy. Enforced expression of miR-let-7c in OCI-AML22 cells resulted in expansion of the LSC-enriched CD34+CD38- population and reduction in the CD34- terminally differentiated blast population. Comparison of miR-let-7c targets from HSPC and AML patient samples revealed the existence of distinct gene regulatory networks. In HSPC, miR-let-7c targets are enriched in cell cycle and regulation of transcription, whereas in AML, mitochondrial ETC, oxidative phosphorylation, and cell-matrix adhesion were the most enriched pathways. Moreover, analysis of miR-let-7c targets identified a high prevalence of transcripts encoded by the mitochondrial (mt) DNA (MT-ATP6, MT-ATP8, MT-ND1, -2, -4, -5, MT-CO1), pointing towards a largely uncharacterized function of miRNAs in regulating transcripts encoded by mtDNA in AML. Collectively, miRNA sequencing in LSC-enriched AML samples and identification of miRNA targetomes in HSPC and AML patient samples have revealed unexpected mechanisms of post-transcriptional regulation in normal and malignant hematopoietic stem cells.

TatDB: a database of Ago1-mediated targets of transfer RNA fragments.

tRNA-derived fragments (tRFs) are a class of emerging post-transcriptional regulators of gene expression likely binding to the transcripts of target genes. However, only a few tRFs targets have been experimentally validated, making it hard to extrapolate the functions or binding mechanisms of tRFs. The paucity of resources supporting the identification of the targets of tRFs creates a bottleneck in the fast-developing field. We have previously analyzed chimeric reads in crosslinked Argonaute1-RNA complexes to help infer the guide-target pairs and binding mechanisms of multiple tRFs based on experimental data in human HEK293 cells. To efficiently disseminate these results to the research community, we designed a web-based database tatDB (targets of tRFs DataBase) populated with close to 250 000 experimentally determined guide-target pairs with ∼23 000 tRF isoforms. tatDB has a user-friendly interface with flexible query options/filters allowing one to obtain comprehensive information on given tRFs (or targets). Modes of interactions are supported by secondary structures of potential guide-target hybrids and binding motifs, essential for understanding the targeting mechanisms of tRFs. Further, we illustrate the value of the database on an example of hypothesis-building for a tRFs potentially involved in the lifecycle of the SARS-CoV-2 virus. tatDB is freely accessible at https://grigoriev-lab.camden.rutgers.edu/tatdb.

SCRAP: a bioinformatic pipeline for the analysis of small chimeric RNA-seq data.

MicroRNAs (miRNAs) are small non-coding RNAs (sncRNAs) that function in post-transcriptional gene regulation through imperfect base pairing with mRNA targets which results in inhibition of translation and typically destabilization of bound transcripts. Sequence-based algorithms historically used to predict miRNA targets face inherent challenges in reliably reflecting in vivo interactions. Recent strategies have directly profiled miRNA-target interactions by crosslinking and ligation of sncRNAs to their targets within the RNA-induced silencing complex (RISC), followed by high throughput sequencing of the chimeric sncRNA:target RNAs. Despite the strength of these direct profiling approaches, standardized pipelines for effectively analyzing the resulting chimeric sncRNA:target RNA sequencing data are not readily available. Here we present SCRAP, a robust Small Chimeric RNA Analysis Pipeline for the bioinformatic processing of chimeric sncRNA:target RNA sequencing data. SCRAP consists of two parts, each of which are specifically optimized for the distinctive characteristics of chimeric small RNA sequencing reads: first, read processing and alignment and second, peak calling and annotation. We apply SCRAP to benchmark chimeric sncRNA:target RNA sequencing datasets generated by distinct molecular approaches, and compare SCRAP to existing chimeric RNA analysis pipelines. SCRAP has minimal hardware requirements, is cross-platform, and contains extensive annotation to broaden accessibility for processing small chimeric RNA sequencing data and enable insights about the targets of small non-coding RNAs in regulating diverse biological systems.

Chimeric DNA byproducts in strand displacement amplification using the T7 replisome.

Recent advances in next generation sequencing technologies enable reading DNA molecules hundreds of kilobases in length and motivate development of DNA amplification methods capable of producing long amplicons. In vivo, DNA replication is performed not by a single polymerase enzyme, but multiprotein complexes called replisomes. Here, we investigate strand-displacement amplification reactions using the T7 replisome, a macromolecular complex of a helicase, a single-stranded DNA binding protein, and a DNA polymerase. The T7 replisome may initiate processive DNA synthesis from DNA nicks, and the reaction of a 48 kilobase linear double stranded DNA substrate with the T7 replisome and nicking endonucleases is shown to produce discrete DNA amplicons. To gain a mechanistic understanding of this reaction, we utilized Oxford Nanopore long-read sequencing technology. Sequence analysis of the amplicons revealed chimeric DNA reads and uncovered a connection between template switching and polymerase exonuclease activity. Nanopore sequencing provides insight to guide the further development of isothermal amplification methods for long DNA, and our results highlight the need for high-specificity, high-turnover nicking endonucleases to initiate DNA amplification without thermal denaturation.

Identification and quantification of chimeric sequencing reads in a highly multiplexed RAD-seq protocol.

Highly multiplexed approaches have become common in genomic studies. They have improved the cost-effectiveness of genotyping hundreds of individuals using combinatorially barcoded adapters. These strategies, however, can potentially misassigned reads to incorrect samples. Here, we used a modified quaddRAD protocol to analyse the occurrence of index hopping and PCR chimeras in a series of experiments with up to 100 multiplexed samples per sequencing lane (639 samples in total). We created two types of sequencing libraries: four libraries of type A, where PCRs were run on individual samples before multiplexing, and three libraries of type B, where PCRs were run on pooled samples. We used fixed pairs of inner barcodes to identify chimeric reads. Type B libraries show a higher percentage of misassigned reads (1.15%) than type A libraries (0.65%). We also quantify the commonly undetectable chimeric sequences that occur whenever multiplexed groups of samples with different outer barcodes are sequenced together on a single flow cell. Our results suggest that these types of chimeric sequences represent up to 1.56% and 1.29% of reads in type A and B libraries, respectively. We also show that increasing the number of mismatches allowed for barcode rescue to above 2 dramatically increases the number of recovered chimeric reads. We provide recommendations for developing highly multiplexed RAD-seq protocols and analysing the resulting data to minimize the generation of chimeric sequences, allowing their quantification and a finer control on the number of PCR cycles necessary to generate enough input DNA for library preparation.