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

The ubiquitin-conjugating enzyme HR6B is required for maintenance of X chromosome silencing in mouse spermatocytes and spermatids

<p>Abstract</p> <p>Background</p> <p>The ubiquitin-conjugating enzyme HR6B is required for spermatogenesis in mouse. Loss of HR6B results in aberrant histone modification patterns on the trancriptionally silenced X and Y chromosomes (XY body) and on centromeric chromatin in meiotic prophase. We studied the relationship between these chromatin modifications and their effects on global gene expression patterns, in spermatocytes and spermatids.</p> <p>Results</p> <p>HR6B is enriched on the XY body and on centromeric regions in pachytene spermatocytes. Global gene expression analyses revealed that spermatid-specific single- and multicopy X-linked genes are prematurely expressed in <it>Hr6b </it>knockout spermatocytes. Very few other differences in gene expression were observed in these cells, except for upregulation of major satellite repeat transcription. In contrast, in <it>Hr6b </it>knockout spermatids, 7298 genes were differentially expressed; 65% of these genes was downregulated, but we observed a global upregulation of gene transcription from the X chromosome. In wild type spermatids, approximately 20% of the single-copy X-linked genes reach an average expression level that is similar to the average expression from autosomes.</p> <p>Conclusions</p> <p>Spermatids maintain an enrichment of repressive chromatin marks on the X chromosome, originating from meiotic prophase, but this does not interfere with transcription of the single-copy X-linked genes that are reactivated or specifically activated in spermatids. HR6B represses major satellite repeat transcription in spermatocytes, and functions in the maintenance of X chromosome silencing in spermatocytes and spermatids. It is discussed that these functions involve modification of chromatin structure, possibly including H2B ubiquitylation.</p

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

Abstract 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

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

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

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