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.
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.
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.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
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.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
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. In this study, we utilized CRISPR/Cas9 genome editing to repair AXIN1-truncated mutations in five HCC cell lines. For each cell line we successfully obtained 2-4 correctly repaired clones, which all show reduced beta-catenin signaling accompanied with reduced cell viability and colony formation. Although exposure of repaired clones to Wnt3A-conditioned medium restored beta-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 beta-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. Our study provides insights into the effects of repairing AXIN1 mutations on beta-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.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
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
itself. The homozygous deletion of this site demonstrates that autoregulation of
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
deficiency, and also other monogenic syndromes.
Mesenchymal stromal cells (MSC) are increasingly used as an investigative therapeutic product for immune disorders and degenerative disease. Typically, MSC are isolated from human tissue, expanded in ...culture, and cryopreserved until usage. The safety and efficacy of MSC therapy will depend on the phenotypical and functional characteristics of MSC. The freeze-thawing procedure may change these characteristics. Furthermore, the cells encounter a microenvironment after administration that may impact their properties. It has been demonstrated that the majority of MSC localize to the lungs after intravenous infusion, making this the site to study the effects of the in vivo milieu on administered MSC. In this study, we investigated the effect of freeze-thawing and the mouse lung microenvironment on human adipose tissue-derived MSC. There were effects of freeze-thawing on the whole genome expression profile of MSC, although the effects did not exceed interdonor differences. There were no major changes in the expression of hemostatic regulators on transcriptional level, but significantly increased expression of procoagulant tissue factor on the surface of thawed adipose MSC, correlating with increased procoagulant activity of thawed cells. Exposure for 2 h to the lung microenvironment had a major effect on MSC gene expression and affected several immunological pathways. This indicates that MSC undergo functional changes shortly after infusion and this may influence the efficacy of MSC to modulate inflammatory responses. The results of this study demonstrate that MSC rapidly alter in response to the local milieu and disease-specific conditions may shape MSC after administration.
Noncoding RNAs have been widely recognized as essential mediators of gene regulation. However, in contrast to protein‐coding genes, much less is known about the influence of noncoding RNAs on human ...diseases. Here we examined the association of genetic variants located in primary microRNA sequences and long noncoding RNAs (lncRNAs) with Alzheimer disease (AD) by leveraging data from the largest genome‐wide association meta‐analysis of late‐onset AD. Variants annotated to 5 miRNAs and 10 lncRNAs (in seven distinct loci) exceeded the Bonferroni‐corrected significance threshold (p < 1.02 × 10−6). Among these, a leading variant (rs2526377:A>G) at the 17q22 locus annotated to two noncoding RNAs (MIR142 and BZRAP1‐AS) was significantly associated with a reduced risk of AD and fulfilled predefined criteria for being a functional variant. Our functional genomic analyses revealed that rs2526377 affects the promoter activity and decreases the expression of miR‐142. Moreover, differential expression analysis by RNA‐Seq in human iPSC‐derived neural progenitor cells and the hippocampus of miR‐142 knockout mice demonstrated multiple target genes of miR‐142 in the brain that are likely to be involved in the inflammatory and neurodegenerative manifestations of AD. These include TGFBR1 and PICALM, of which their derepression in the brain due to reduced expression levels of miR‐142‐3p may reduce the risk of AD.
We examined the association of genetic variants located in primary microRNA sequences and long noncoding RNAs with Alzheimer disease (AD) by leveraging data from the largest genome‐wide association meta‐analysis of late‐onset AD. Five miRNAs and ten lncRNAs in seven distinct loci exceeded the Bonferroni‐corrected significance threshold, including a new susceptibility locus at 17q22. The top SNP annotated to MIR142 and BZRAP1‐AS was significantly associated with a reduced risk of AD and fulfilled predefined criteria for being functional variant. Our functional genomic analyses revealed that rs2526377 affects promoter activity and decreases the expression of miR‐142. Moreover, differential expression analysis by RNA‐Seq in human iPSC‐derived neural progenitor cells and the hippocampus of miR‐142 knockout mice demonstrated multiple miR‐142‐3p target genes (e.g., TGFBR1 and PICALM) in the brain that are likely to be involved in the inflammatory and neurodegenerative manifestations of AD.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
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.
Cohesin is crucial for genome stability, cell division, transcription and chromatin organization. Its functions critically depend on NIPBL, the cohesin-loader protein that is found to be mutated in ...>60% of the cases of Cornelia de Lange syndrome (CdLS). Other mutations are described in the cohesin subunits SMC1A, RAD21, SMC3 and the HDAC8 protein. In 25-30% of CdLS cases no mutation in the known CdLS genes is detected. Until now, functional elements in the noncoding genome were not characterized in the molecular etiology of CdLS and therefore are excluded from mutation screening, although the impact of such mutations has now been recognized for a wide range of diseases. We have identified different elements of the noncoding genome involved in regulation of the NIPBL gene. NIPBL-AS1 is a long non-coding RNA transcribed upstream and antisense to NIPBL. By knockdown and transcription blocking experiments, we could show that not the NIPBL-AS1 gene product, but its actual transcription is important to regulate NIPBL expression levels. This reveals a possibility to boost the transcriptional activity of the NIPBL gene by interfering with the NIPBL-AS1 lncRNA. Further, we have identified a novel distal enhancer regulating both NIPBL and NIPBL-AS1. Deletion of the enhancer using CRISPR genome editing in HEK293T cells reduces expression of NIPBL, NIPBL-AS1 as well as genes found to be dysregulated in CdLS.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The domestic ferret,
Mustela putorius furo
, is an important mammalian animal model to study human respiratory infection. However, insufficient genomic annotation hampers detailed studies of ferret T ...cell responses. In this study, we analyzed the published T cell receptor beta (TRB) locus and performed high-throughput sequencing (HTS) of peripheral blood of four healthy adult ferrets to identify expressed V, D, J, and C genes. The HTS data is used as a guide to manually curate the expressed V, D, J, and C genes. The ferret locus appears to be most similar to that of the dog. Like other mammalian TRB loci, the ferret TRB locus contains a library of variable genes located upstream of two D-J-C gene clusters, followed by a (in the ferret non-functional) V gene with an inverted transcriptional orientation. All TRB genes (expressed or not) reported here have been approved by the IMGT/WHO-IUIS nomenclature committee.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OBVAL, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ