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.

Epigenomic partitioning of a polygenic risk score for asthma reveals distinct genetically driven disease pathways.

BACKGROUND: Individual differences in susceptibility to develop asthma, a heterogeneous chronic inflammatory lung disease, are poorly understood. It remains debated whether genetics can predict asthma risk and how genetic variants modulate the complex pathophysiology of asthma. AIM: To build polygenic risk scores (PRSs) for asthma risk prediction and epigenomically link predictive genetic variants to pathophysiological mechanisms. METHODS: Restricted PRSs were constructed using single nucleotide variants derived from genome-wide association studies and validated using data generated in the Rotterdam Study, a Dutch prospective cohort of 14 926 individuals. Outcomes used were asthma, childhood-onset asthma (COA), adulthood-onset asthma (AOA), eosinophilic asthma, and asthma exacerbations. Genome-wide chromatin analysis data from 19 disease-relevant cell types were used for epigenomic PRS partitioning. RESULTS: PRSs obtained predicted asthma and related outcomes, with the strongest associations observed for COA (2.55 odds ratios per PRS standard deviation, area under the curve of 0.760). PRSs allowed for the classification of individuals into high and low-risk groups. PRS partitioning using epigenomic profiles identified 5 clusters of variants within putative gene regulatory regions linked to specific asthma-relevant cells, genes, and biological pathways. CONCLUSIONS: PRSs were associated with asthma(-related traits) in a Dutch prospective cohort, with substantially higher predictive power observed for COA than for AOA. Importantly, PRS variants could be epigenomically partitioned into clusters of regulatory variants with different pathophysiological association patterns and effect estimates, which likely represent distinct genetically driven disease pathways. Our findings have potential implications for personalized risk mitigation and treatment strategies.

LSD1/KDM1A and GFI1B repress endothelial fate and induce hematopoietic fate in induced pluripotent stem cell-derived hemogenic endotheliums.

Differentiation of induced pluripotent stem cells (iPSCs) into hematopoietic lineages offers great therapeutic potential. During embryogenesis, hemogenic endothelium (HE) gives rise to hematopoietic stem and progenitor cells through the endothelial-to-hematopoietic transition (EHT). Understanding this process using iPSCs is key to generating functional hematopoietic stem cells (HSCs), a currently unmet challenge. In this study, we examined the role of the transcriptional factor GFI1B and its co-factor LSD1/KDM1A in EHT. To this end, we employed patient-derived iPSC lines with a dominant negative dysfunctional GFI1BQ287* and irreversible pharmacological LSD1/KDM1A inhibition in healthy iPSC lines. The formation of HE remained unaffected; however, hematopoietic output was severely reduced in both conditions. Single-cell RNA sequencing (scRNAseq) performed on the CD144+/CD31+ population derived from healthy iPSCs revealed similar expression dynamics of genes associated with in vivo EHT. Interestingly, LSD1/KDM1A inhibition in healthy lines before EHT resulted in a complete absence of hematopoietic output. However, uncommitted HE cells did not display GFI1B expression, suggesting a timed transcriptional program. To test this hypothesis, we ectopically expressed GFI1B in uncommitted HE cells, leading to downregulation of endothelial genes and upregulation of hematopoietic genes, including GATA2, KIT, RUNX1, and SPI1. Thus, we demonstrate that LSD1/KDM1A and GFI1B can function at distinct temporal points in different cellular subsets during EHT. Although GFI1B is not detected in uncommitted HE cells, its ectopic expression allows for partial hematopoietic specification. These data indicate that precisely timed expression of specific transcriptional regulators during EHT is crucial to the eventual outcome of EHT.

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.

Unraveling the impact of AXIN1 mutations on HCC development: Insights from CRISPR/Cas9 repaired AXIN1-mutant liver cancer cell lines.

BACKGROUND: Hepatocellular carcinoma (HCC) is a highly aggressive liver cancer with significant morbidity and mortality rates. AXIN1 is one of the top-mutated genes in HCC, but the mechanism by which AXIN1 mutations contribute to HCC development remains unclear. METHODS: In this study, we utilized CRISPR/Cas9 genome editing to repair AXIN1-truncated mutations in five HCC cell lines. RESULTS: For each cell line we successfully obtained 2-4 correctly repaired clones, which all show reduced ß-catenin signaling accompanied with reduced cell viability and colony formation. Although exposure of repaired clones to Wnt3A-conditioned medium restored ß-catenin signaling, it did not or only partially recover their growth characteristics, indicating the involvement of additional mechanisms. Through RNA-sequencing analysis, we explored the gene expression patterns associated with repaired AXIN1 clones. Except for some highly-responsive ß-catenin target genes, no consistent alteration in gene/pathway expression was observed. This observation also applies to the Notch and YAP/TAZ-Hippo signaling pathways, which have been associated with AXIN1-mutant HCCs previously. The AXIN1-repaired clones also cannot confirm a recent observation that AXIN1 is directly linked to YAP/TAZ protein stability and signaling. CONCLUSIONS: Our study provides insights into the effects of repairing AXIN1 mutations on ß-catenin signaling, cell viability, and colony formation in HCC cell lines. However, further investigations are necessary to understand the complex mechanisms underlying HCC development associated with AXIN1 mutations.

Genome-wide methylation analysis in patients with proximal hypospadias - a pilot study and review of the literature.

In patients with proximal hypospadias, often no genetic cause is identified despite extensive genetic testing. Many genes involved in sex development encode transcription factors with strict timing and dosing of the gene products. We hypothesised that there might be recurrent differences in DNA methylation in boys with hypospadias and that these might differ between patients born small versus appropriate for gestational age. Genome-wide Methylated DNA sequencing (MeD-seq) was performed on 32bp LpnPI restriction enzyme fragments after RE-digestion in leucocytes from 16 XY boys with unexplained proximal hypospadias, one with an unexplained XX testicular disorder/difference of sex development (DSD) and twelve, healthy, sex- and age-matched controls. Five of seven differentially methylated regions (DMRs) between patients and XY controls were in the Long Intergenic Non-Protein Coding RNA 665 (LINC00665; CpG24525). Three patients showed hypermethylation of MAP3K1. Finally, no DMRs in XX testicular DSD associated genes were identified in the XX boy versus XX controls. In conclusion, we observed no recognizable epigenetic signature in 16 boys with XY proximal hypospadias and no difference between children born small versus appropriate for gestational age. Comparison to previous methylation studies in individuals with hypospadias did not show consistent findings, possibly due to the use of different inclusion criteria, tissues and methods.

Distinctive cell-free DNA methylation characterizes presymptomatic genetic frontotemporal dementia.

OBJECTIVE: Methylation of plasma cell-free DNA (cfDNA) has potential as a marker of brain damage in neurodegenerative diseases such as frontotemporal dementia (FTD). Here, we study methylation of cfDNA in presymptomatic and symptomatic carriers of genetic FTD pathogenic variants, next to healthy controls. METHODS: cfDNA was isolated from cross-sectional plasma of 10 presymptomatic carriers (4 C9orf72, 4 GRN, and 2 MAPT), 10 symptomatic carriers (4 C9orf72, 4 GRN, and 2 MAPT), and 9 healthy controls. Genome-wide methylation of cfDNA was determined using a high-resolution sequencing technique (MeD-seq). Cumulative scores based on the identified differentially methylated regions (DMRs) were estimated for presymptomatic carriers (vs. controls and symptomatic carriers), and reevaluated in a validation cohort (8 presymptomatic: 3 C9orf72, 3 GRN, and 2 MAPT; 26 symptomatic: 7 C9orf72, 6 GRN, 12 MAPT, and 1 TARDBP; 13 noncarriers from genetic FTD families). RESULTS: Presymptomatic carriers showed a distinctive methylation profile compared to healthy controls and symptomatic carriers. Cumulative DMR scores in presymptomatic carriers enabled to significantly differentiate presymptomatic carriers from healthy controls (p < 0.001) and symptomatic carriers (p < 0.001). In the validation cohort, these scores differentiated presymptomatic carriers from symptomatic carriers (p ≤ 0.007) only. Transcription-start-site methylation in presymptomatic carriers, generally associated with gene downregulation, was enriched for genes involved in ubiquitin-dependent processes, while gene body methylation, generally associated with gene upregulation, was enriched for genes involved in neuronal cell processes. INTERPRETATION: A distinctive methylation profile of cfDNA characterizes the presymptomatic stage of genetic FTD, and could reflect neuronal death in this stage.

Runx1+ vascular smooth muscle cells are essential for hematopoietic stem and progenitor cell development in vivo.

Hematopoietic stem cells (HSCs) produce all essential cellular components of the blood. Stromal cell lines supporting HSCs follow a vascular smooth muscle cell (vSMC) differentiation pathway, suggesting that some hematopoiesis-supporting cells originate from vSMC precursors. These pericyte-like precursors were recently identified in the aorta-gonad-mesonephros (AGM) region; however, their role in the hematopoietic development in vivo remains unknown. Here, we identify a subpopulation of NG2+Runx1+ perivascular cells that display a sclerotome-derived vSMC transcriptomic profile. We show that deleting Runx1 in NG2+ cells impairs the hematopoietic development in vivo and causes transcriptional changes in pericytes/vSMCs, endothelial cells and hematopoietic cells in the murine AGM. Importantly, this deletion leads also to a significant reduction of HSC reconstitution potential in the bone marrow in vivo. This defect is developmental, as NG2+Runx1+ cells were not detected in the adult bone marrow, demonstrating the existence of a specialised pericyte population in the HSC-generating niche, unique to the embryo.

Epigenetic and Genomic Hallmarks of PARP-Inhibitor Resistance in Ovarian Cancer Patients.

BACKGROUND: Patients with advanced-stage epithelial ovarian cancer (EOC) receive treatment with a poly-ADP ribose-polymerase (PARP) inhibitor (PARPi) as maintenance therapy after surgery and chemotherapy. Unfortunately, many patients experience disease progression because of acquired therapy resistance. This study aims to characterize epigenetic and genomic changes in cell-free DNA (cfDNA) associated with PARPi resistance. MATERIALS AND METHODS: Blood was taken from 31 EOC patients receiving PARPi therapy before treatment and at disease progression during/after treatment. Resistance was defined as disease progression within 6 months after starting PARPi and was seen in fifteen patients, while sixteen patients responded for 6 to 42 months. Blood cfDNA was evaluated via Modified Fast Aneuploidy Screening Test-Sequencing System (mFast-SeqS to detect aneuploidy, via Methylated DNA Sequencing (MeD-seq) to find differentially methylated regions (DMRs), and via shallow whole-genome and -exome sequencing (shWGS, exome-seq) to define tumor fractions and mutational signatures. RESULTS: Aneuploid cfDNA was undetectable pre-treatment but observed in six patients post-treatment, in five resistant and one responding patient. Post-treatment ichorCNA analyses demonstrated in shWGS and exome-seq higher median tumor fractions in resistant (7% and 9%) than in sensitive patients (7% and 5%). SigMiner analyses detected predominantly mutational signatures linked to mismatch repair and chemotherapy. DeSeq2 analyses of MeD-seq data revealed three methylation signatures and more tumor-specific DMRs in resistant than in responding patients in both pre- and post-treatment samples (274 vs. 30 DMRs, 190 vs. 57 DMRs, Χ2-test p < 0.001). CONCLUSION: Our genome-wide Next-Generation Sequencing (NGS) analyses in PARPi-resistant patients identified epigenetic differences in blood before treatment, whereas genomic alterations were more frequently observed after progression. The epigenetic differences at baseline are especially interesting for further exploration as putative predictive biomarkers for PARPi resistance.

A cellular reporter system to evaluate endogenous fetal hemoglobin induction and screen for therapeutic compounds.

Reactivation of fetal hemoglobin expression alleviates the symptoms associated with ß-globinopathies, severe hereditary diseases with significant global health implications due to their high morbidity and mortality rates. The symptoms emerge following the postnatal transition from fetal-to-adult hemoglobin expression. Extensive research has focused on inducing the expression of the fetal γ-globin subunit to reverse this switch and ameliorate these symptoms. Despite decades of research, only one compound, hydroxyurea, found its way to the clinic as an inducer of fetal hemoglobin. Unfortunately, its efficacy varies among patients, highlighting the need for more effective treatments. Erythroid cell lines have been instrumental in the pursuit of both pharmacological and genetic ways to reverse the postnatal hemoglobin switch. Here, we describe the first endogenously tagged fetal hemoglobin reporter cell line based on the adult erythroid progenitor cell line HUDEP2. Utilizing CRISPR-Cas9-mediated knock-in, a bioluminescent tag was integrated at the HBG1 gene. Subsequent extensive characterization confirmed that the resulting reporter cell line closely mirrors the HUDEP2 characteristics and that the cells report fetal hemoglobin induction with high sensitivity and specificity. This novel reporter cell line is therefore highly suitable for evaluating genetic and pharmacologic strategies to induce fetal hemoglobin. Furthermore, it provides an assay compatible with high-throughput drug screening, exemplified by the identification of a cluster of known fetal hemoglobin inducers in a pilot study. This new tool is made available to the research community, with the aspiration that it will accelerate the search for safer and more effective strategies to reverse the hemoglobin switch.

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.