BART-Seq: cost-effective massively parallelized targeted sequencing for genomics, transcriptomics, and single-cell analysis

CRISPR screening in hematology research: from bulk to single-cell level

The CRISPR genome editing technology has revolutionized the way gene function is studied. Genome editing can be achieved in single genes or for thousands of genes simultaneously in sensitive genetic screens. While conventional genetic screens are limited to bulk measurements of cell behavior, recent developments in single-cell technologies make it possible to combine CRISPR screening with single-cell profiling. In this way, cell behavior and gene expression can be monitored simultaneously, with the additional possibility of including data on chromatin accessibility and protein levels. Moreover, the availability of various Cas proteins leading to inactivation, activation, or other effects on gene function further broadens the scope of such screens. The integration of single-cell multi-omics approaches with CRISPR screening open the path to high-content information on the impact of genetic perturbations at single-cell resolution. Current limitations in cell throughput and data density need to be taken into consideration, but new technologies are rapidly evolving and are likely to easily overcome these limitations. In this review, we discuss the use of bulk CRISPR screening in hematology research, as well as the emergence of single-cell CRISPR screening and its added value to the field.

New approach methodologies to facilitate and improve the hazard assessment of non-genotoxic carcinogens-a PARC project

Carcinogenic chemicals, or their metabolites, can be classified as genotoxic or non-genotoxic carcinogens (NGTxCs). Genotoxic compounds induce DNA damage, which can be detected by an established in vitro and in vivo battery of genotoxicity assays. For NGTxCs, DNA is not the primary target, and the possible modes of action (MoA) of NGTxCs are much more diverse than those of genotoxic compounds, and there is no specific in vitro assay for detecting NGTxCs. Therefore, the evaluation of the carcinogenic potential is still dependent on long-term studies in rodents. This 2-year bioassay, mainly applied for testing agrochemicals and pharmaceuticals, is time-consuming, costly and requires very high numbers of animals. More importantly, its relevance for human risk assessment is questionable due to the limited predictivity for human cancer risk, especially with regard to NGTxCs. Thus, there is an urgent need for a transition to new approach methodologies (NAMs), integrating human-relevant in vitro assays and in silico tools that better exploit the current knowledge of the multiple processes involved in carcinogenesis into a modern safety assessment toolbox. Here, we describe an integrative project that aims to use a variety of novel approaches to detect the carcinogenic potential of NGTxCs based on different mechanisms and pathways involved in carcinogenesis. The aim of this project is to contribute suitable assays for the safety assessment toolbox for an efficient and improved, internationally recognized hazard assessment of NGTxCs, and ultimately to contribute to reliable mechanism-based next-generation risk assessment for chemical carcinogens.

Adding gene transcripts into genomic prediction improves accuracy and reveals sampling time dependence

Recent developments allowed generating multiple high-quality 'omics' data that could increase the predictive performance of genomic prediction for phenotypes and genetic merit in animals and plants. Here, we have assessed the performance of parametric and nonparametric models that leverage transcriptomics in genomic prediction for 13 complex traits recorded in 478 animals from an outbred mouse population. Parametric models were implemented using the best linear unbiased prediction, while nonparametric models were implemented using the gradient boosting machine algorithm. We also propose a new model named GTCBLUP that aims to remove between-omics-layer covariance from predictors, whereas its counterpart GTBLUP does not do that. While gradient boosting machine models captured more phenotypic variation, their predictive performance did not exceed the best linear unbiased prediction models for most traits. Models leveraging gene transcripts captured higher proportions of the phenotypic variance for almost all traits when these were measured closer to the moment of measuring gene transcripts in the liver. In most cases, the combination of layers was not able to outperform the best single-omics models to predict phenotypes. Using only gene transcripts, the gradient boosting machine model was able to outperform best linear unbiased prediction for most traits except body weight, but the same pattern was not observed when using both single nucleotide polymorphism genotypes and gene transcripts. Although the GTCBLUP model was not able to produce the most accurate phenotypic predictions, it showed the highest accuracies for breeding values for 9 out of 13 traits. We recommend using the GTBLUP model for prediction of phenotypes and using the GTCBLUP for prediction of breeding values.

Sample-multiplexing approaches for single-cell sequencing

Single-cell sequencing is widely used in biological and medical studies. However, its application with multiple samples is hindered by inefficient sample processing, high experimental costs, ambiguous identification of true single cells, and technical batch effects. Here, we introduce sample-multiplexing approaches for single-cell sequencing in transcriptomics, epigenomics, genomics, and multiomics. In single-cell transcriptomics, sample multiplexing uses variants of native or artificial features as sample markers, enabling sample pooling and decoding. Such features include: (1) natural genetic variation, (2) nucleotide-barcode anchoring on cellular or nuclear membranes, (3) nucleotide-barcode internalization to the cytoplasm or nucleus, (4) vector-based barcode expression in cells, and (5) nucleotide-barcode incorporation during library construction. Other single-cell omics methods are based on similar concepts, particularly single-cell combinatorial indexing. These methods overcome current challenges, while enabling super-loading of single cells. Finally, selection guidelines are presented that can accelerate technological application.

Prime-seq, efficient and powerful bulk RNA sequencing

Cost-efficient library generation by early barcoding has been central in propelling single-cell RNA sequencing. Here, we optimize and validate prime-seq, an early barcoding bulk RNA-seq method. We show that it performs equivalently to TruSeq, a standard bulk RNA-seq method, but is fourfold more cost-efficient due to almost 50-fold cheaper library costs. We also validate a direct RNA isolation step, show that intronic reads are derived from RNA, and compare cost-efficiencies of available protocols. We conclude that prime-seq is currently one of the best options to set up an early barcoding bulk RNA-seq protocol from which many labs would profit.

Dissecting intercellular and intracellular signaling networks with barcoded genetic tools

The power of next-generation sequencing has stimulated the development of many analysis techniques for transcriptomics and genomics. More recently, the concept of 'molecular barcoding' has broadened the spectrum of sequencing-based applications to dissect different aspects of intracellular and intercellular signaling. In these assay formats, barcode reporters replace standard reporter genes. The virtually infinitive number of expressed barcode sequences allows high levels of multiplexing, hence accelerating experimental progress. Furthermore, reporter barcodes are used to quantitatively monitor a variety of biological events in living cells which has already provided much insight into complex cellular signaling and will further increase our knowledge in the future.

Multiplexing Methods for Simultaneous Large-Scale Transcriptomic Profiling of Samples at Single-Cell Resolution

Barcoding technology has greatly improved the throughput of cells and genes detected in single-cell RNA sequencing (scRNA-seq) studies. Recently, increasing studies have paid more attention to the use of this technology to increase the throughput of samples, as it has greatly reduced the processing time, technical batch effects, and library preparation costs, and lowered the per-sample cost. In this review, the various DNA-based barcoding methods for sample multiplexing are focused on, specifically, on the four major barcoding strategies. A detailed comparison of the barcoding methods is also presented, focusing on aspects such as sample/cell throughput and gene detection, and guidelines for choosing the most appropriate barcoding technique according to the personalized requirements are developed. Finally, the critical applications of sample multiplexing and technical challenges in combinatorial labeling, barcoding in vivo, and multimodal tagging at the spatially resolved resolution, as well as, the future prospects of multiplexed scRNA-seq, for example, prioritizing and predicting the severity of coronavirus disease 2019 (COVID-19) in patients of different gender and age are highlighted.

scPNMF: sparse gene encoding of single cells to facilitate gene selection for targeted gene profiling

Motivation

Single-cell RNA sequencing (scRNA-seq) captures whole transcriptome information of individual cells. While scRNA-seq measures thousands of genes, researchers are often interested in only dozens to hundreds of genes for a closer study. Then, a question is how to select those informative genes from scRNA-seq data. Moreover, single-cell targeted gene profiling technologies are gaining popularity for their low costs, high sensitivity and extra (e.g. spatial) information; however, they typically can only measure up to a few hundred genes. Then another challenging question is how to select genes for targeted gene profiling based on existing scRNA-seq data.

Results

Here, we develop the single-cell Projective Non-negative Matrix Factorization (scPNMF) method to select informative genes from scRNA-seq data in an unsupervised way. Compared with existing gene selection methods, scPNMF has two advantages. First, its selected informative genes can better distinguish cell types. Second, it enables the alignment of new targeted gene profiling data with reference data in a low-dimensional space to facilitate the prediction of cell types in the new data. Technically, scPNMF modifies the PNMF algorithm for gene selection by changing the initialization and adding a basis selection step, which selects informative bases to distinguish cell types. We demonstrate that scPNMF outperforms the state-of-the-art gene selection methods on diverse scRNA-seq datasets. Moreover, we show that scPNMF can guide the design of targeted gene profiling experiments and the cell-type annotation on targeted gene profiling data.

Availability and implementation

The R package is open-access and available at https://github.com/JSB-UCLA/scPNMF. The data used in this work are available at Zenodo: https://doi.org/10.5281/zenodo.4797997.

Supplementary information

Supplementary data are available at Bioinformatics online.

Merging Genomics and Transcriptomics for Predicting Fusarium Head Blight Resistance in Wheat

Genomic selection with genome-wide distributed molecular markers has evolved into a well-implemented tool in many breeding programs during the last decade. The resistance against Fusarium head blight (FHB) in wheat is probably one of the most thoroughly studied systems within this framework. Aside from the genome, other biological strata like the transcriptome have likewise shown some potential in predictive breeding strategies but have not yet been investigated for the FHB-wheat pathosystem. The aims of this study were thus to compare the potential of genomic with transcriptomic prediction, and to assess the merit of blending incomplete transcriptomic with complete genomic data by the single-step method. A substantial advantage of gene expression data over molecular markers has been observed for the prediction of FHB resistance in the studied diversity panel of breeding lines and released cultivars. An increase in prediction ability was likewise found for the single-step predictions, although this can mostly be attributed to an increased accuracy among the RNA-sequenced genotypes. The usage of transcriptomics can thus be seen as a complement to already established predictive breeding pipelines with pedigree and genomic data, particularly when more cost-efficient multiplexing techniques for RNA-sequencing will become more accessible in the future.

HyPR-seq: Single-cell quantification of chosen RNAs via hybridization and sequencing of DNA probes

Single-cell quantification of RNAs is important for understanding cellular heterogeneity and gene regulation, yet current approaches suffer from low sensitivity for individual transcripts, limiting their utility for many applications. Here we present Hybridization of Probes to RNA for sequencing (HyPR-seq), a method to sensitively quantify the expression of hundreds of chosen genes in single cells. HyPR-seq involves hybridizing DNA probes to RNA, distributing cells into nanoliter droplets, amplifying the probes with PCR, and sequencing the amplicons to quantify the expression of chosen genes. HyPR-seq achieves high sensitivity for individual transcripts, detects nonpolyadenylated and low-abundance transcripts, and can profile more than 100,000 single cells. We demonstrate how HyPR-seq can profile the effects of CRISPR perturbations in pooled screens, detect time-resolved changes in gene expression via measurements of gene introns, and detect rare transcripts and quantify cell-type frequencies in tissue using low-abundance marker genes. By directing sequencing power to genes of interest and sensitively quantifying individual transcripts, HyPR-seq reduces costs by up to 100-fold compared to whole-transcriptome single-cell RNA-sequencing, making HyPR-seq a powerful method for targeted RNA profiling in single cells.

Oligonucleotide Probe for Transcriptome in Vivo Analysis (TIVA) of Single Neurons with Minimal Background

Messenger RNA (mRNA) isolated from single cells can generate powerful biological insights, including the discovery of new cell types with unique functions as well as markers potentially predicting a cell's response to various therapeutic agents. We previously introduced an oligonucleotide-based technique for site-selective, photoinduced biotinylation and capture of mRNA within a living cell called transcriptome in vivo analysis (TIVA). Successful application of the TIVA technique hinges upon its oligonucleotide probe remaining completely inert (or "caged") to mRNA unless photoactivated. To improve the reliability of TIVA probe caging in diverse and challenging biological conditions, we applied a rational design process involving iterative modifications to the oligonucleotide construct. In this work, we discuss these design motivations and present an optimized probe with minimal background binding to mRNA prior to photolysis. We assess its caging performance through multiple in vitro assays including FRET analysis, native gel comigration, and pull down with model mRNA transcripts. Finally, we demonstrate that this improved probe can also isolate mRNA from single living neurons in brain tissue slices with excellent caging control.

Cytologically targeted next-generation sequencing: a synergy for diagnosing urothelial carcinoma

INTRODUCTION

Cytology and cystoscopy are used to detect urothelial carcinoma (UC), but together they still fail to detect some UC cases and are not suitable for screening asymptomatic individuals. Mutations are present in more than 98% of UC, mutations have therapeutic significance, and they can be detected by next generation sequencing (NGS) in urine samples. We review the role of NGS in UC detection.

MATERIALS AND METHODS

Comprehensive literature review on UC genetics, economics of NGS, and previous reports of UC detection by NGS.

RESULTS

The raw costs of NGS have decreased to about 14,000 base pairs per penny, making it appear economically feasible to use NGS widely. Reported NGS assays fall short of predicted sensitivity. Decreased sensitivity is attributed to a low frequency of mutant alleles in many urine samples. Attempts to increase the percentage of mutant alleles, by using cell-free urinary DNA, or by using cell sorting and microfluidics, have been unsuccessful or remain unproven. However, cytologic examination can immediately enable NGS: Urine cytologies with sufficient proportions of abnormal cells could be directly triaged to NGS with high sensitivity for UC detection. For cases with a low proportion of abnormal cells, cytologically targeted microdissection of cells for NGS should maintain sensitivity and decrease sequencing costs. Cytologically targeted urothelial cells for NGS could allow a screening test for low grade UC.

CONCLUSIONS

Cytology is immediately poised to allow NGS to improve the diagnosis of UC, allowing NGS to be an ancillary test for atypical cytologies, and potentially allowing a screening test for low-grade UC.

Targeted Perturb-seq enables genome-scale genetic screens in single cells

The transcriptome contains rich information on molecular, cellular and organismal phenotypes. However, experimental and statistical limitations constrain sensitivity and throughput of genetic screening with single-cell transcriptomics readout. To overcome these limitations, we introduce targeted Perturb-seq (TAP-seq), a sensitive, inexpensive and platform-independent method focusing single-cell RNA-seq coverage on genes of interest, thereby increasing the sensitivity and scale of genetic screens by orders of magnitude. TAP-seq permits routine analysis of thousands of CRISPR-mediated perturbations within a single experiment, detects weak effects and lowly expressed genes, and decreases sequencing requirements by up to 50-fold. We apply TAP-seq to generate perturbation-based enhancer-target gene maps for 1,778 enhancers within 2.5% of the human genome. We thereby show that enhancer-target association is jointly determined by three-dimensional contact frequency and epigenetic states, allowing accurate prediction of enhancer targets throughout the genome. In addition, we demonstrate that TAP-seq can identify cell subtypes with only 100 sequencing reads per cell.

Targeted Perturb-seq enables genome-scale genetic screens in single cells

The transcriptome contains rich information on molecular, cellular and organismal phenotypes. However, experimental and statistical limitations constrain sensitivity and throughput of genetic screening with single-cell transcriptomics readout. To overcome these limitations, we introduce targeted Perturb-seq (TAP-seq), a sensitive, inexpensive and platform-independent method focusing single-cell RNA-seq coverage on genes of interest, thereby increasing the sensitivity and scale of genetic screens by orders of magnitude. TAP-seq permits routine analysis of thousands of CRISPR-mediated perturbations within a single experiment, detects weak effects and lowly expressed genes, and decreases sequencing requirements by up to 50-fold. We apply TAP-seq to generate perturbation-based enhancer-target gene maps for 1,778 enhancers within 2.5% of the human genome. We thereby show that enhancer-target association is jointly determined by three-dimensional contact frequency and epigenetic states, allowing accurate prediction of enhancer targets throughout the genome. In addition, we demonstrate that TAP-seq can identify cell subtypes with only 100 sequencing reads per cell.

Multiparametric Assays for Accelerating Early Drug Discovery

Drug discovery campaigns are hampered by substantial attrition rates largely due to a lack of efficacy and safety reasons associated with candidate drugs. This is true in particular for genetically complex diseases, where insufficient knowledge of the modulatory actions of candidate drugs on targets and entire target pathways further adds to the problem of attrition. To better profile compound actions on targets, potential off-targets, and disease-linked pathways, new innovative technologies need to be developed that can elucidate the complex cellular signaling networks in health and disease. Here, we discuss progress in genetically encoded multiparametric assays and mass spectrometry (MS)-based proteomics, which both represent promising toolkits to profile multifactorial actions of drug candidates in disease-relevant cellular systems to promote drug discovery and personalized medicine.

BART-Seq: cost-effective massively parallelized targeted sequencing for genomics, transcriptomics, and single-cell analysis

We describe a highly sensitive, quantitative, and inexpensive technique for targeted sequencing of transcript cohorts or genomic regions from thousands of bulk samples or single cells in parallel. Multiplexing is based on a simple method that produces extensive matrices of diverse DNA barcodes attached to invariant primer sets, which are all pre-selected and optimized in silico. By applying the matrices in a novel workflow named Barcode Assembly foR Targeted Sequencing (BART-Seq), we analyze developmental states of thousands of single human pluripotent stem cells, either in different maintenance media or upon Wnt/β-catenin pathway activation, which identifies the mechanisms of differentiation induction. Moreover, we apply BART-Seq to the genetic screening of breast cancer patients and identify BRCA mutations with very high precision. The processing of thousands of samples and dynamic range measurements that outperform global transcriptomics techniques makes BART-Seq first targeted sequencing technique suitable for numerous research applications.