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.

Deconvoluting multi-person biological mixtures and accurate characterization and identification of separated contributors using non-targeted single-cell DNA sequencing.

The genetic characterization and identification of individuals who contributed to biological mixtures are complex and mostly unresolved tasks. These tasks are relevant in various fields, particularly in forensic investigations, which frequently encounters crime scene stains generated by more than one person. Currently, forensic mixture deconvolution is mostly performed subsequent to forensic DNA profiling at the level of the mixed DNA profiles, which comes with several limitations. Some previous studies attempted at separating single cells prior to forensic DNA profiling. However, these approaches are biased at selection of the cells and, due to their targeted DNA analysis on low template DNA, provide incomplete and unreliable forensic DNA profiles. We recently demonstrated the feasibility of performing mixture deconvolution prior to forensic DNA profiling through the utilization of a non-targeted single-cell transcriptome sequencing (scRNA-seq). In addition to individual-specific mixture deconvolution, this approach also allowed accurate characterisation of biological sex, biogeographic ancestry and individual identification of the separated mixture contributors. However, RNA has the forensic disadvantage of being prone to degradation, and sequencing RNA - focussing on coding regions - limits the number of single nucleotide polymorphisms (SNPs) utilized for genetic mixture deconvolution, characterization, and identification. These limitations can be overcome by performing single-cell sequencing on the level of DNA instead of RNA. Here, for the first time, we applied non-targeted single-cell DNA sequencing (scDNA-seq) by applying the scATAC-seq (Assay for Transposase-Accessible Chromatin with sequencing) technique to address the challenges of mixture deconvolution in the forensic context. We demonstrated that scATAC-seq, together with our recently developed De-goulash data analysis pipeline, is capable of deconvoluting complex blood mixtures of five individuals from both sexes with varying biogeographic ancestries. We further showed that our approach achieved correct genetic characterization of the biological sex and the biogeographic ancestry of each of the separated mixture contributors and established their identity. Furthermore, by analysing in-silico generated scATAC-seq data mixtures, we demonstrated successful individual-specific mixture deconvolution of i) highly complex mixtures of 11 individuals, ii) balanced mixtures containing as few as 20 cells (10 per each individual), and iii) imbalanced mixtures with a ratio as low as 1:80. Overall, our proof-of-principle study demonstrates the general feasibility of scDNA-seq in general, and scATAC-seq in particular, for mixture deconvolution, genetic characterization and individual identification of the separated mixture contributors. Furthermore, it shows that compared to scRNA-seq, scDNA-seq detects more SNPs from fewer cells, providing higher sensitivity, that is valuable in forensic genetics.

Single-cell transcriptome sequencing allows genetic separation, characterization and identification of individuals in multi-person biological mixtures.

Identifying individuals from biological mixtures to which they contributed is highly relevant in crime scene investigation and various biomedical research fields, but despite previous attempts, remains nearly impossible. Here we investigated the potential of using single-cell transcriptome sequencing (scRNA-seq), coupled with a dedicated bioinformatics pipeline (De-goulash), to solve this long-standing problem. We developed a novel approach and tested it with scRNA-seq data that we de-novo generated from multi-person blood mixtures, and also in-silico mixtures we assembled from public single individual scRNA-seq datasets, involving different numbers, ratios, and bio-geographic ancestries of contributors. For all 2 up to 9-person balanced and imbalanced blood mixtures with ratios up to 1:60, we achieved a clear single-cell separation according to the contributing individuals. For all separated mixture contributors, sex and bio-geographic ancestry (maternal, paternal, and bi-parental) were correctly determined. All separated contributors were correctly individually identified with court-acceptable statistical certainty using de-novo generated whole exome sequencing reference data. In this proof-of-concept study, we demonstrate the feasibility of single-cell approaches to deconvolute biological mixtures and subsequently genetically characterise, and individually identify the separated mixture contributors. With further optimisation and implementation, this approach may eventually allow moving to challenging biological mixtures, including those found at crime scenes.

Retrospective analysis of enhancer activity and transcriptome history.

Cell state changes in development and disease are controlled by gene regulatory networks, the dynamics of which are difficult to track in real time. In this study, we used an inducible DCM-RNA polymerase subunit b fusion protein which labels active genes and enhancers with a bacterial methylation mark that does not affect gene transcription and is propagated in S-phase. This DCM-RNA polymerase fusion protein enables transcribed genes and active enhancers to be tagged and then examined at later stages of development or differentiation. We apply this DCM-time machine (DCM-TM) technology to study intestinal homeostasis, revealing rapid and coordinated activation of enhancers and nearby genes during enterocyte differentiation. We provide new insights in absorptive-secretory lineage decision-making in intestinal stem cell (ISC) differentiation and show that ISCs retain a unique chromatin landscape required to maintain ISC identity and delineate future expression of differentiation-associated genes. DCM-TM has wide applicability in tracking cell states, providing new insights in the regulatory networks underlying cell state changes.

DNA methylation and hydroxymethylation characterize the identity of D1 and D2 striatal projection neurons.

Neuronal DNA modifications differ from those in other cells, including methylation outside CpG context and abundant 5-hydroxymethylation whose relevance for neuronal identities are unclear. Striatal projection neurons expressing D1 or D2 dopamine receptors allow addressing this question, as they share many characteristics but differ in their gene expression profiles, connections, and functional roles. We compare translating mRNAs and DNA modifications in these two populations. DNA methylation differences occur predominantly in large genomic clusters including differentially expressed genes, potentially important for D1 and D2 neurons. Decreased gene body methylation is associated with higher gene expression. Hydroxymethylation differences are more scattered and affect transcription factor binding sites, which can influence gene expression. We also find a strong genome-wide hydroxymethylation asymmetry between the two DNA strands, particularly pronounced at expressed genes and retrotransposons. These results identify novel properties of neuronal DNA modifications and unveil epigenetic characteristics of striatal projection neurons heterogeneity.

Endothelial Zeb2 preserves the hepatic angioarchitecture and protects against liver fibrosis.

AIMS: Hepatic capillaries are lined with specialized liver sinusoidal endothelial cells (LSECs) which support macromolecule passage to hepatocytes and prevent fibrosis by keeping hepatic stellate cells (HSCs) quiescent. LSEC specialization is co-determined by transcription factors. The zinc-finger E-box-binding homeobox (Zeb)2 transcription factor is enriched in LSECs. Here, we aimed to elucidate the endothelium-specific role of Zeb2 during maintenance of the liver and in liver fibrosis. METHODS AND RESULTS: To study the role of Zeb2 in liver endothelium we generated EC-specific Zeb2 knock-out (ECKO) mice. Sequencing of liver EC RNA revealed that deficiency of Zeb2 results in prominent expression changes in angiogenesis-related genes. Accordingly, the vascular area was expanded and the presence of pillars inside ECKO liver vessels indicated that this was likely due to increased intussusceptive angiogenesis. LSEC marker expression was not profoundly affected and fenestrations were preserved upon Zeb2 deficiency. However, an increase in continuous EC markers suggested that Zeb2-deficient LSECs are more prone to dedifferentiation, a process called 'capillarization'. Changes in the endothelial expression of ligands that may be involved in HSC quiescence together with significant changes in the expression profile of HSCs showed that Zeb2 regulates LSEC-HSC communication and HSC activation. Accordingly, upon exposure to the hepatotoxin carbon tetrachloride (CCl4), livers of ECKO mice showed increased capillarization, HSC activation, and fibrosis compared to livers from wild-type littermates. The vascular maintenance and anti-fibrotic role of endothelial Zeb2 was confirmed in mice with EC-specific overexpression of Zeb2, as the latter resulted in reduced vascularity and attenuated CCl4-induced liver fibrosis. CONCLUSION: Endothelial Zeb2 preserves liver angioarchitecture and protects against liver fibrosis. Zeb2 and Zeb2-dependent genes in liver ECs may be exploited to design novel therapeutic strategies to attenuate hepatic fibrosis.

Comprehensive multi-omics integration identifies differentially active enhancers during human brain development with clinical relevance.

BACKGROUND: Non-coding regulatory elements (NCREs), such as enhancers, play a crucial role in gene regulation, and genetic aberrations in NCREs can lead to human disease, including brain disorders. The human brain is a complex organ that is susceptible to numerous disorders; many of these are caused by genetic changes, but a multitude remain currently unexplained. Understanding NCREs acting during brain development has the potential to shed light on previously unrecognized genetic causes of human brain disease. Despite immense community-wide efforts to understand the role of the non-coding genome and NCREs, annotating functional NCREs remains challenging. METHODS: Here we performed an integrative computational analysis of virtually all currently available epigenome data sets related to human fetal brain. RESULTS: Our in-depth analysis unravels 39,709 differentially active enhancers (DAEs) that show dynamic epigenomic rearrangement during early stages of human brain development, indicating likely biological function. Many of these DAEs are linked to clinically relevant genes, and functional validation of selected DAEs in cell models and zebrafish confirms their role in gene regulation. Compared to enhancers without dynamic epigenomic rearrangement, DAEs are subjected to higher sequence constraints in humans, have distinct sequence characteristics and are bound by a distinct transcription factor landscape. DAEs are enriched for GWAS loci for brain-related traits and for genetic variation found in individuals with neurodevelopmental disorders, including autism. CONCLUSION: This compendium of high-confidence enhancers will assist in deciphering the mechanism behind developmental genetics of human brain and will be relevant to uncover missing heritability in human genetic brain disorders.

The Bone-Forming Properties of Periosteum-Derived Cells Differ Between Harvest Sites.

The development of alternatives for autologous bone grafts is a major focus of bone tissue engineering. To produce living bone-forming implants, skeletal stem and progenitor cells (SSPCs) are envisioned as key ingredients. SSPCs can be obtained from different tissues including bone marrow, adipose tissue, dental pulp, and periosteum. Human periosteum-derived cells (hPDCs) exhibit progenitor cell characteristics and have well-documented in vivo bone formation potency. Here, we have characterized and compared hPDCs derived from tibia with craniofacial hPDCs, from maxilla and mandible, respectively, each representing a potential source for cell-based tissue engineered implants for craniofacial applications. Maxilla and mandible-derived hPDCs display similar growth curves as tibial hPDCs, with equal trilineage differentiation potential toward chondrogenic, osteogenic, and adipogenic cells. These craniofacial hPDCs are positive for SSPC-markers CD73, CD164, and Podoplanin (PDPN), and negative for CD146, hematopoietic and endothelial lineage markers. Bulk RNA-sequencing identified genes that are differentially expressed between the three sources of hPDC. In particular, differential expression was found for genes of the HOX and DLX family, for SOX9 and genes involved in skeletal system development. The in vivo bone formation, 8 weeks after ectopic implantation in nude mice, was observed in constructs seeded with tibial and mandibular hPDCs. Taken together, we provide evidence that hPDCs show different profiles and properties according to their anatomical origin, and that craniofacial hPDCs are potential sources for cell-based bone tissue engineering strategies. The mandible-derived hPDCs display - both in vitro and in vivo - chondrogenic and osteogenic differentiation potential, which supports their future testing for use in craniofacial bone regeneration applications.

An evolutionarily ancient mechanism for regulation of hemoglobin expression in vertebrate red cells.

The oxygen transport function of hemoglobin (HB) is thought to have arisen ∼500 million years ago, roughly coinciding with the divergence between jawless (Agnatha) and jawed (Gnathostomata) vertebrates. Intriguingly, extant HBs of jawless and jawed vertebrates were shown to have evolved twice, and independently, from different ancestral globin proteins. This raises the question of whether erythroid-specific expression of HB also evolved twice independently. In all jawed vertebrates studied to date, one of the HB gene clusters is linked to the widely expressed NPRL3 gene. Here we show that the nprl3-linked hb locus of a jawless vertebrate, the river lamprey (Lampetra fluviatilis), shares a range of structural and functional properties with the equivalent jawed vertebrate HB locus. Functional analysis demonstrates that an erythroid-specific enhancer is located in intron 7 of lamprey nprl3, which corresponds to the NPRL3 intron 7 MCS-R1 enhancer of jawed vertebrates. Collectively, our findings signify the presence of an nprl3-linked multiglobin gene locus, which contains a remote enhancer that drives globin expression in erythroid cells, before the divergence of jawless and jawed vertebrates. Different globin genes from this ancestral cluster evolved in the current NPRL3-linked HB genes in jawless and jawed vertebrates. This provides an explanation of the enigma of how, in different species, globin genes linked to the same adjacent gene could undergo convergent evolution.

MAU2 and NIPBL Variants Impair the Heterodimerization of the Cohesin Loader Subunits and Cause Cornelia de Lange Syndrome.

The NIPBL/MAU2 heterodimer loads cohesin onto chromatin. Mutations in NIPBL account for most cases of the rare developmental disorder Cornelia de Lange syndrome (CdLS). Here we report a MAU2 variant causing CdLS, a deletion of seven amino acids that impairs the interaction between MAU2 and the NIPBL N terminus. Investigating this interaction, we discovered that MAU2 and the NIPBL N terminus are largely dispensable for normal cohesin and NIPBL function in cells with a NIPBL early truncating mutation. Despite a predicted fatal outcome of an out-of-frame single nucleotide duplication in NIPBL, engineered in two different cell lines, alternative translation initiation yields a form of NIPBL missing N-terminal residues. This form cannot interact with MAU2, but binds DNA and mediates cohesin loading. Altogether, our work reveals that cohesin loading can occur independently of functional NIPBL/MAU2 complexes and highlights a novel mechanism protective against out-of-frame mutations that is potentially relevant for other genetic conditions.

Integrative and perturbation-based analysis of the transcriptional dynamics of TGFß/BMP system components in transition from embryonic stem cells to neural progenitors.

Cooperative actions of extrinsic signals and cell-intrinsic transcription factors alter gene regulatory networks enabling cells to respond appropriately to environmental cues. Signaling by transforming growth factor type ß (TGFß) family ligands (eg, bone morphogenetic proteins [BMPs] and Activin/Nodal) exerts cell-type specific and context-dependent transcriptional changes, thereby steering cellular transitions throughout embryogenesis. Little is known about coordinated regulation and transcriptional interplay of the TGFß system. To understand intrafamily transcriptional regulation as part of this system's actions during development, we selected 95 of its components and investigated their mRNA-expression dynamics, gene-gene interactions, and single-cell expression heterogeneity in mouse embryonic stem cells transiting to neural progenitors. Interrogation at 24 hour intervals identified four types of temporal gene transcription profiles that capture all stages, that is, pluripotency, epiblast formation, and neural commitment. Then, between each stage we performed esiRNA-based perturbation of each individual component and documented the effect on steady-state mRNA levels of the remaining 94 components. This exposed an intricate system of multilevel regulation whereby the majority of gene-gene interactions display a marked cell-stage specific behavior. Furthermore, single-cell RNA-profiling at individual stages demonstrated the presence of detailed co-expression modules and subpopulations showing stable co-expression modules such as that of the core pluripotency genes at all stages. Our combinatorial experimental approach demonstrates how intrinsically complex transcriptional regulation within a given pathway is during cell fate/state transitions.

Id4 promotes the elimination of the pro-activation factor Ascl1 to maintain quiescence of adult hippocampal stem cells.

Quiescence is essential for the long-term maintenance of adult stem cells but how stem cells maintain quiescence is poorly understood. Here, we show that neural stem cells (NSCs) in the adult mouse hippocampus actively transcribe the pro-activation factor Ascl1 regardless of their activated or quiescent states. We found that the inhibitor of DNA binding protein Id4 is enriched in quiescent NSCs and that elimination of Id4 results in abnormal accumulation of Ascl1 protein and premature stem cell activation. Accordingly, Id4 and other Id proteins promote elimination of Ascl1 protein in NSC cultures. Id4 sequesters Ascl1 heterodimerization partner E47, promoting Ascl1 protein degradation and stem cell quiescence. Our results highlight the importance of non-transcriptional mechanisms for the maintenance of NSC quiescence and reveal a role for Id4 as a quiescence-inducing factor, in contrast with its role of promoting the proliferation of embryonic neural progenitors.

DamC reveals principles of chromatin folding in vivo without crosslinking and ligation.

Current understanding of chromosome folding is largely reliant on chromosome conformation capture (3C)-based experiments, where chromosomal interactions are detected as ligation products after chromatin crosslinking. To measure chromosome structure in vivo, quantitatively and without crosslinking and ligation, we implemented a modified version of DNA adenine methyltransferase identification (DamID) named DamC, which combines DNA methylation-based detection of chromosomal interactions with next-generation sequencing and biophysical modeling of methylation kinetics. DamC performed in mouse embryonic stem cells provides the first in vivo validation of the existence of topologically associating domains (TADs), CTCF loops and confirms 3C-based measurements of the scaling of contact probabilities. Combining DamC with transposon-mediated genomic engineering shows that new loops can be formed between ectopic and endogenous CTCF sites, which redistributes physical interactions within TADs. DamC provides the first crosslinking- and ligation-free demonstration of the existence of key structural features of chromosomes and provides novel insights into how chromosome structure within TADs can be manipulated.

Genomes of Ellobius species provide insight into the evolutionary dynamics of mammalian sex chromosomes.

The X and Y sex chromosomes of placental mammals show hallmarks of a tumultuous evolutionary past. The X Chromosome has a rich and conserved gene content, while the Y Chromosome has lost most of its genes. In the Transcaucasian mole vole Ellobius lutescens, the Y Chromosome including Sry has been lost, and both females and males have a 17,X diploid karyotype. Similarly, the closely related Ellobius talpinus, has a 54,XX karyotype in both females and males. Here, we report the sequencing and assembly of the E. lutescens and E. talpinus genomes. The results indicate that the loss of the Y Chromosome in E. lutescens and E. talpinus occurred in two independent events. Four functional homologs of mouse Y-Chromosomal genes were detected in both female and male E. lutescens, of which three were also detected in the E. talpinus genome. One of these is Eif2s3y, known as the only Y-derived gene that is crucial for successful male meiosis. Female and male E. lutescens can carry one and the same X Chromosome with a largely conserved gene content, including all genes known to function in X Chromosome inactivation. The availability of the genomes of these mole vole species provides unique models to study the dynamics of sex chromosome evolution.

A whole-genome sequence and transcriptome perspective on HER2-positive breast cancers.

HER2-positive breast cancer has long proven to be a clinically distinct class of breast cancers for which several targeted therapies are now available. However, resistance to the treatment associated with specific gene expressions or mutations has been observed, revealing the underlying diversity of these cancers. Therefore, understanding the full extent of the HER2-positive disease heterogeneity still remains challenging. Here we carry out an in-depth genomic characterization of 64 HER2-positive breast tumour genomes that exhibit four subgroups, based on the expression data, with distinctive genomic features in terms of somatic mutations, copy-number changes or structural variations. The results suggest that, despite being clinically defined by a specific gene amplification, HER2-positive tumours melt into the whole luminal-basal breast cancer spectrum rather than standing apart. The results also lead to a refined ERBB2 amplicon of 106 kb and show that several cases of amplifications are compatible with a breakage-fusion-bridge mechanism.

Incomplete meiotic sex chromosome inactivation in the domestic dog.

BACKGROUND: In mammalian meiotic prophase, homologous chromosome recognition is aided by formation and repair of programmed DNA double-strand breaks (DSBs). Subsequently, stable associations form through homologous chromosome synapsis. In male mouse meiosis, the largely heterologous X and Y chromosomes synapse only in their short pseudoautosomal regions (PARs), and DSBs persist along the unsynapsed non-homologous arms of these sex chromosomes. Asynapsis of these arms and the persistent DSBs then trigger transcriptional silencing through meiotic sex chromosome inactivation (MSCI), resulting in formation of the XY body. This inactive state is partially maintained in post-meiotic haploid spermatids (postmeiotic sex chromatin repression, PSCR). For the human, establishment of MSCI and PSCR have also been reported, but X-linked gene silencing appears to be more variable compared to mouse. To gain more insight into the regulation and significance of MSCI and PSCR among different eutherian species, we have performed a global analysis of XY pairing dynamics, DSB repair, MSCI and PSCR in the domestic dog (Canis lupus familiaris), for which the complete genome sequence has recently become available, allowing a thorough comparative analyses. RESULTS: In addition to PAR synapsis between X and Y, we observed extensive self-synapsis of part of the dog X chromosome, and rapid loss of known markers of DSB repair from that part of the X. Sequencing of RNA from purified spermatocytes and spermatids revealed establishment of MSCI. However, the self-synapsing region of the X displayed higher X-linked gene expression compared to the unsynapsed area in spermatocytes, and was post-meiotically reactivated in spermatids. In contrast, genes in the PAR, which are expected to escape MSCI, were expressed at very low levels in both spermatocytes and spermatids. Our comparative analysis was then used to identify two X-linked genes that may escape MSCI in spermatocytes, and 21 that are specifically re-activated in spermatids of human, mouse and dog. CONCLUSIONS: Our data indicate that MSCI is incomplete in the dog. This may be partially explained by extensive, but transient, self-synapsis of the X chromosome, in association with rapid completion of meiotic DSB repair. In addition, our comparative analysis identifies novel candidate male fertility genes.