Split Pool Ligation-based Single-cell Transcriptome sequencing (SPLiT-seq) data processing pipeline comparison

Background: Single-cell sequencing techniques are revolutionizing every field of biology by providing the ability to measure the abundance of biological molecules at a single-cell resolution. Although single-cell sequencing approaches have been developed for several molecular modalities, single-cell transcriptome sequencing is the most prevalent and widely applied technique. SPLiT-seq (split-pool ligation-based transcriptome sequencing) is one of these single-cell transcriptome techniques that applies a unique combinatorial-barcoding approach by splitting and pooling cells into multi-well plates containing barcodes. This unique approach required the development of dedicated computational tools to preprocess the data and extract the count matrices. Here we compare eight bioinformatic pipelines (alevin-fry splitp, LR-splitpipe, SCSit, splitpipe, splitpipeline, SPLiTseq-demultiplex, STARsolo and zUMI) that have been developed to process SPLiT-seq data. We provide an overview of the tools, their computational performance, functionality and impact on downstream processing of the single-cell data, which vary greatly depending on the tool used. Results: We show that STARsolo, splitpipe and alevin-fry splitp can all handle large amount of data within reasonable time. In contrast, the other five pipelines are slow when handling large datasets. When using smaller dataset, cell barcode results are similar with the exception of SPLiTseq-demultiplex and splitpipeline. LR-splitpipe that is originally designed for processing long-read sequencing data is the slowest of all pipelines. Alevin-fry produced different down-stream results that are difficult to interpret. STARsolo functions nearly identical to splitpipe and produce results that are highly similar to each other. However, STARsolo lacks the function to collapse random hexamer reads for which some additional coding is required. Conclusion: Our comprehensive comparative analysis aids users in selecting the most suitable analysis tool for efficient SPLiT-seq data processing, while also detailing the specific prerequisites for each of these pipelines. From the available pipelines, we recommend splitpipe or STARSolo for SPLiT-seq data analysis.</p

Species-specific responses during Seoul orthohantavirus infection in human and rat lung microvascular endothelial cells

Seoul orthohantavirus (SEOV) is a rat-borne zoonotic virus that is transmitted via inhalation of aerosolized infectious excreta, and can cause hemorrhagic fever with renal syndrome (HFRS) in humans worldwide. In rats, SEOV predominantly exists as a persistent infection in the absence of overt clinical signs. Lack of disease in rats is attributed to downregulation of pro-inflammatory and upregulation of regulatory host responses. As lung microvascular endothelial cells (LMECs) represent a primary target of infection in both human and rats, infections in these cells provide a unique opportunity to study the central role of LMECs in the dichotomy between pathogenicity in both species. In this study, host responses to SEOV infection in primary human and rat LMECs were directly compared on a transcriptional level. As infection of rat LMECs was more efficient than human LMECs, the majority of anti-viral defense responses were observed earlier in rat LMECs. Most prominently, SEOV-induced processes in both species included responses to cytokine stimulus, negative regulation of innate immune responses, responses to type I and II interferons, regulation of pattern recognition receptor signaling and MHC-I signaling. However, over time, in the rat LMECs, responses shifted from an anti-viral state towards a more immunotolerant state displayed by a PD-L1, B2M-, JAK2-focused interaction network aiding in negative regulation of cytotoxic CD8-positive T cell activation. This suggests a novel mechanism by which species-specific orthohantavirus-induced endothelium and T cell crosstalk may play a crucial role in the development of acute disease in humans and persistence in rodents.</p

Species-specific responses during Seoul orthohantavirus infection in human and rat lung microvascular endothelial cells

Seoul orthohantavirus (SEOV) is a rat-borne zoonotic virus that is transmitted via inhalation of aerosolized infectious excreta, and can cause hemorrhagic fever with renal syndrome (HFRS) in humans worldwide. In rats, SEOV predominantly exists as a persistent infection in the absence of overt clinical signs. Lack of disease in rats is attributed to downregulation of pro-inflammatory and upregulation of regulatory host responses. As lung microvascular endothelial cells (LMECs) represent a primary target of infection in both human and rats, infections in these cells provide a unique opportunity to study the central role of LMECs in the dichotomy between pathogenicity in both species. In this study, host responses to SEOV infection in primary human and rat LMECs were directly compared on a transcriptional level. As infection of rat LMECs was more efficient than human LMECs, the majority of anti-viral defense responses were observed earlier in rat LMECs. Most prominently, SEOV-induced processes in both species included responses to cytokine stimulus, negative regulation of innate immune responses, responses to type I and II interferons, regulation of pattern recognition receptor signaling and MHC-I signaling. However, over time, in the rat LMECs, responses shifted from an anti-viral state towards a more immunotolerant state displayed by a PD-L1, B2M-, JAK2-focused interaction network aiding in negative regulation of cytotoxic CD8-positive T cell activation. This suggests a novel mechanism by which species-specific orthohantavirus-induced endothelium and T cell crosstalk may play a crucial role in the development of acute disease in humans and persistence in rodents.</p

Comparison of Single Cell Transcriptome Sequencing Methods:Of Mice and Men

Single cell RNAseq has been a big leap in many areas of biology. Rather than investigating gene expression on a whole organism level, this technology enables scientists to get a detailed look at rare single cells or within their cell population of interest. The field is growing, and many new methods appear each year. We compared methods utilized in our core facility: Smart-seq3, PlexWell, FLASH-seq, VASA-seq, SORT-seq, 10X, Evercode, and HIVE. We characterized the equipment requirements for each method. We evaluated the performances of these methods based on detected features, transcriptome diversity, mitochondrial RNA abundance and multiplets, among others and benchmarked them against bulk RNA sequencing. Here, we show that bulk transcriptome detects more unique transcripts than any single cell method. While most methods are comparable in many regards, FLASH-seq and VASA-seq yielded the best metrics, e.g., in number of features. If no equipment for automation is available or many cells are desired, then HIVE or 10X yield good results. In general, more recently developed methods perform better. This also leads to the conclusion that older methods should be phased out, and that the development of single cell RNAseq methods is still progressing considerably.</p

Zeb2 DNA-Binding Sites in Neuroprogenitor Cells Reveal Autoregulation and Affirm Neurodevelopmental Defects, Including in Mowat-Wilson Syndrome

Functional perturbation and action mechanism studies have shown that the transcription factor Zeb2 controls cell fate decisions, differentiation, and/or maturation in multiple cell lineages in embryos and after birth. In cultured embryonic stem cells (ESCs), Zeb2’s mRNA/protein upregulation is necessary for the exit from primed pluripotency and for entering general and neural differentiation. We edited mouse ESCs to produce Flag-V5 epitope-tagged Zeb2 protein from one endogenous allele. Using chromatin immunoprecipitation coupled with sequencing (ChIP-seq), we mapped 2432 DNA-binding sites for this tagged Zeb2 in ESC-derived neuroprogenitor cells (NPCs). A new, major binding site maps promoter-proximal to Zeb2 itself. The homozygous deletion of this site demonstrates that autoregulation of Zeb2 is necessary to elicit the appropriate Zeb2-dependent effects in ESC-to-NPC differentiation. We have also cross-referenced all the mapped Zeb2 binding sites with previously obtained transcriptome data from Zeb2 perturbations in ESC-derived NPCs, GABAergic interneurons from the ventral forebrain of mouse embryos, and stem/progenitor cells from the post-natal ventricular-subventricular zone (V-SVZ) in mouse forebrain, respectively. Despite the different characteristics of each of these neurogenic systems, we found interesting target gene overlaps. In addition, our study also contributes to explaining developmental disorders, including Mowat-Wilson syndrome caused by ZEB2 deficiency, and also other monogenic syndromes.</p

Nimbus: A design-driven analyses suite for amplicon-based NGS data

Motivation: PCR-based DNA enrichment followed by massively parallel sequencing is a straightforward and cost effective method to sequence genes up to high depth. The full potential of ampliconbased sequencing assays is currently not achieved as analysis methods do not take into account the source amplicons of the detected variants. Tracking the source amplicons has the potential to identify systematic biases, enhance variant calling and improve the designs of future assays. Results: We present Nimbus, a software suite for the analysis of amplicon-based sequencing data. Nimbus includes tools for data pre-processing, alignment, single nucleotide polymorphism (SNP), insertion and deletion calling, quality control and visualization. Nimbus can detect SNPs in its alignment seeds and reduces alignment issues by the usage of decoy amplicons. Tracking the amplicons throughout analysis allows easy and fast design optimization by amplicon performance comparison. It enables detection of probable false positive variants present in a single amplicon from real variants present in multiple amplicons and provides multiple sample visualization. Nimbus was tested using HaloPlex Exome datasets and outperforms other callers for low-frequency variants. The variants called by Nimbus were highly concordant between twin samples and SNP-arrays. The Nimbus suite provides an end-to-end solution for variant calling, design optimization and visualization of amplicon-derived next-generation sequencing datasets

Targeted chromatin conformation analysis identifies novel distal neural enhancers of ZEB2 in pluripotent stem cell differentiation

The transcription factor zinc finger E-box binding protein 2 (ZEB2) controls embryonic and adult cell fate decisions and cellular maturation in many stem/progenitor cell types. Defects in these processes in specific cell types underlie several aspects of Mowat-Wilson syndrome (MOWS), which is caused by ZEB2 haplo-insufficiency. Human ZEB2, like mouse Zeb2, is located on chromosome 2 downstream of a ±3.5 Mb-long gene-desert, lacking any protein-coding gene. Using temporal targeted chromatin capture (T2C), we show major chromatin structural changes based on mapping in-cis proximities between the ZEB2 promoter and this gene desert during neural differentiation of human-induced pluripotent stem cells, including at early neuroprogenitor cell (NPC)/rosette state, where ZEB2 mRNA levels increase significantly. Combining T2C with histone-3 acetylation mapping, we identified three novel candidate enhancers about 500 kb upstream of the ZEB2 transcription start site. Functional luciferase-based assays in heterologous cells and NPCs reveal co-operation between these three enhancers. This study is the first to document in-cis Regulatory Elements located in ZEB2's gene desert. The results further show the usability of T2C for future studies of ZEB2 REs in differentiation and maturation of multiple cell types and the molecular characterization of newly identified MOWS patients that lack mutations in ZEB2 protein-coding exons

Exome sequencing and functional analyses suggest that SIX6 is a gene involved in an altered proliferation-differentiation balance early in life and optic nerve degeneration at old age

Primary open-angle glaucoma (POAG) is a hereditary neurodegenerative disease, characterized by optic nerve changes including increased excavation, notching and optic disc hemorrhages. The excavation can be described by the vertical cup-disc ratio (VCDR). Previously, genome-wide significant evidence for the association of rs10483727 in SIX1-SIX6 locus with VCDR and subsequent POAG was found. Using 1000 genomes-based imputation of four independent population-based cohorts in the Netherlands, we identified a missense variant rs33912345 (His141Asn) in SIX6 associated with VCDR (Pmeta = 7.74 × 10-7, n = 11 473) and POAG (Pmeta = 6.09 × 10-3, n = 292). Exome sequencing analysis revealed another missense variant rs146737847 (Glu129Lys) also in SIX6 associated with VCDR (P = 5.09 × 10-3, n = 1208). These two findings point to SIX6 as the responsible gene for the previously reported association signal. Functional characterization of SIX6 in zebrafish revealed that knockdown of six6b led to a small eye phenotype. Histological analysis showed retinal lamination, implying an apparent normal development of the eye, but an underdeveloped lens, and reduced optic nerve diameter. Expression analysis of morphants at 3 dpf showed a 5.5-fold up-regulation of cdkn2b, a cyclin-dependent kinase inhibitor, involved in cell cycle regulation and previously associated with VCDR and POAG in genome-wide association studies (GWASs). Since both six6b and cdkn2b play a key role in cell proliferation, we assessed the proliferative activity in the eye of morphants and found an alteration in the proliferative pattern of retinal cells. Our findings in humans and zebrafish suggest a functional involvement of six6b in early eye development, and open new insights into the genetic architecture of POAG

Exome-wide meta-analysis identifies rare 3'-UTR variant in ERCC1/CD3EAP associated with symptoms of sleep apnea

Obstructive sleep apnea (OSA) is a common sleep breathing disorder associated with an increased risk of cardiovascular and cerebrovascular diseases and mortality. Although OSA is fairly heritable (~40%), there have been only few studies looking into the genetics of OSA. In the present study, we aimed to identify genetic variants associated with symptoms of sleep apnea by performing a whole-exome sequence meta-analysis of symptoms of sleep apnea in 1,475 individuals of European descent. We identified 17 rare genetic variants with at least suggestive evidence of significance. Replication in an independent dataset confirmed the association of a rare genetic variant (rs2229918; minor allele frequency = 0.3%) with symptoms of sleep apnea (p-valuemeta = 6.98 × 10-9, ßmeta = 0.99). Rs2229918 overlaps with the 3' untranslated regions of ERCC1 and CD3EAP genes on chromosome 19q13. Both genes are expressed in tissues in the neck area, such as the tongue, muscles, cartilage and the trachea. Further, CD3EAP is localized in the nucleus and mitochondria and involved in the tumor necrosis factor-alpha/nuclear factor kappa B signaling pathway. Our results and biological functions of CD3EAP/ERCC1 genes suggest that the 19q13 locus is interesting for further OSA research

Deciphering the RNA landscape by RNAome sequencing

Current RNA expression profiling methods rely on enrichment steps for specific RNA classes, thereby not detecting all RNA species in an unperturbed manner. We report strand-specific RNAome sequencing that determines expression of small and large RNAs from rRNA-depleted total RNA in a single sequence run. Since current analysis pipelines cannot reliably analyze small and large RNAs simultaneously, we developed TRAP, Total Rna Analysis Pipeline, a robust interface that is also compatible with existing RNA sequencing protocols. RNAome sequencing quantitatively preserved all RNA classes, allowing cross-class comparisons that facilitates the identification of relationships between different RNA classes. We demonstrate the strength of RNAome sequencing in mouse embryonic stem cells treated with cisplatin. MicroRNA and mRNA expression in RNAome sequencing significantly correlated between replicates and was in concordance with both existing RNA sequencing methods and gene expression arrays generated from the same samples. Moreover, RNAome sequencing also detected additional RNA classes such as enhancer RNAs, anti-sense RNAs, novel RNA species and numerous differentially expressed RNAs undetectable by other methods. At the level of complete RNA classes, RNAome sequencing also identified a specific global repression of the microRNA and microRNA isoform classes after cisplatin treatment whereas all other classes such as mRNAs were unchanged. These characteristics of RNAome sequencing will significantly improve expression analysis as well as studies on RNA biology not covered by existing methods