Search

SUMMARYThis study establishes the homeodomain only protein, HOPX, as a determinant controlling the molecular switch between cardiomyocyte progenitor and maturation gene programs. Time-course single-cell gene expression with genome-wide footprinting reveal that HOPX interacts with and controls core cardiac networks by regulating the activity of mutually exclusive developmental gene programs. Upstream hypertrophy and proliferation pathways compete to regulate HOPX transcription. Mitogenic signals override hypertrophic growth signals to suppress HOPX and maintain cardiomyocyte progenitor gene programs. Physiological studies show HOPX directly governs genetic control of cardiomyocyte cell stress responses, electro-mechanical coupling, proliferation, and contractility. We use human genome-wide association studies (GWAS) to show that genetic variation in the HOPX-regulome is significantly associated with complex traits affecting cardiac structure and function. Collectively, this study provides a mechanistic link situating HOPX between competing upstream pathways where HOPX acts as a molecular switch controlling gene regulatory programs underpinning metabolic, signaling, and functional maturation of cardiomyocytes.

Summary Here, we outline detailed protocols to isolate and profile murine splenic dendritic cells (DCs) through advanced flow cytometry of the myeloid compartment and single-cell transcriptomic profiling with integrated cell surface protein expression through CITE-seq. This protocol provides a general transferrable road map for different tissues and species. For complete details on the use and execution of this protocol, please refer to Lukowski et al. (2021)., Graphical abstract, Highlights • Protocol to obtain integrated single-cell gene and protein expression data • Optimized flow cytometry panel for confident delineation of six main myeloid lineages • Gating strategy identifies large cell state heterogeneity within each lineage, Here, we outline detailed protocols to isolate and profile murine splenic dendritic cells (DCs) through advanced flow cytometry of the myeloid compartment and single-cell transcriptomic profiling with integrated cell surface protein expression through CITE-seq. This protocol provides a general transferrable road map for different tissues and species.

We present Nanopore-DamID, a method to simultaneously detect cytosine methylation and DNA-protein interactions from single molecules, via selective sequencing of adenine-labelled DNA. Assaying LaminB1 and CTCF binding with Nanopore-DamID, we identify escape from LAD-associated repression of hypomethylated promoters amidst generalised hypermethylation of LaminB1-associated regulatory elements. We detect novel CTCF binding sites in highly repetitive regions, and allele-specific CTCF binding to imprinted genes and the active X chromosome. Nanopore-DamID highlights the importance of DNA methylation to transcription factor activity.

With the exception of a few master transcription factors, regulators of neutrophil maturation are poorly annotated in the intermediate phenotypes between the granulocyte-macrophage progenitor (GMP) and the mature neutrophil phenotype. Additional challenges in identifying gene expression regulators in differentiation pathways relate to challenges wherein starting cell populations are heterogeneous in lineage potential and development, are spread across various states of quiescence, as well as sample quality and input limitations. These factors contribute to data variability make it difficult to draw simple regulatory inferences. In response we have applied a multi-omics approach using primary blood progenitor cells primed for homogeneous proliferation and granulocyte differentiation states which combines whole transcriptome resequencing (Ampliseq RNA) supported by droplet digital PCR (ddPCR) validation and mass spectrometry-based proteomics in a hypothesis-generation study of neutrophil differentiation pathways. Primary CD34+ cells isolated from human cord blood were first precultured in non-lineage driving medium to achieve an active, proliferating phenotype from which a neutrophil primed progenitor was isolated and cultured in neutrophil lineage supportive medium. Samples were then taken at 24-hour intervals over 9 days and analysed by Ampliseq RNA and mass spectrometry. The Ampliseq dataset depth, breadth and quality allowed for several unexplored transcriptional regulators and ncRNAs to be identified using a combinatorial approach of hierarchical clustering, enriched transcription factor binding motifs, and network mapping. Network mapping in particular increased comprehension of neutrophil differentiation regulatory relationships by implicating ARNT, NHLH1, PLAG1, and 6 non-coding RNAs associated with PU.1 regulation as cell-engineering targets with the potential to increase total neutrophil culture output. Overall, this study develops and demonstrates an effective new hypothesis generation methodology for transcriptome profiling during differentiation, thereby enabling identification of novel gene targets for editing interventions.

Summary: Conventional dendritic cells (cDCs) are traditionally subdivided into cDC1 and cDC2 lineages. Batf3 is a cDC1-required transcription factor, and we observed that Batf3−/− mice harbor a population of cDC1-like cells co-expressing cDC2-associated surface molecules. Using single-cell RNA sequencing with integrated cell surface protein expression (CITE-seq), we found that Batf3−/− mitotic immature cDC1-like cells showed reduced expression of cDC1 features and increased levels of cDC2 features. In wild type, we also observed a proportion of mature cDC1 cells expressing surface features characteristic to cDC2 and found that overall cDC cell state heterogeneity was mainly driven by developmental stage, proliferation, and maturity. We detected population diversity within Sirpa+ cDC2 cells, including a Cd33+ cell state expressing high levels of Sox4 and lineage-mixed features characteristic to cDC1, cDC2, pDCs, and monocytes. In conclusion, these data suggest that multiple cDC cell states can co-express lineage-overlapping features, revealing a level of previously unappreciated cDC plasticity.

Understanding the molecular mechanisms underlying mammalian kidney function requires transcriptome profiling of the interplay between cells comprising nephron segments. Traditional transcriptomics requires cell dissociation, resulting in loss of the spatial context of gene expression within native tissue. To address this problem, we performed spatial transcriptomics (ST) to retain the spatial context of the transcriptome in human and mouse kidneys. The generated ST data allowed spatially resolved differential gene expression analysis, spatial identification of functional nephron segments, cell-to-cell interaction analysis, and chronic kidney disease-associated genetic variant calling. Novel ST thus provides an opportunity to enhance kidney diagnostics and knowledge, by retaining the spatial context of gene expression within intact tissue.

Acute gastrointestinal (GI) graft-versus-host disease (GVHD) is a primary determinant of mortality after allogeneic hematopoietic stem cell transplantation (alloSCT). The condition is mediated by alloreactive donor CD4(+) T cells that differentiate into pathogenic subsets expressing IFN-γ, IL-17A, or GM-CSF and is regulated by subsets expressing IL-10 and/or Foxp3. Developmental relationships between Th cell states during priming in mesenteric lymph nodes (mLNs) and effector function in the GI tract remain undefined at genome scale. We applied scRNA-Seq and computational modeling to a mouse model of donor DC-mediated GVHD exacerbation, creating an atlas of putative CD4(+) T cell differentiation pathways in vivo. Computational trajectory inference suggested emergence of pathogenic and regulatory states along a single developmental trajectory in mLNs. Importantly, we inferred an unexpected second trajectory, categorized by little proliferation or cytokine expression, reduced glycolysis, and high tcf7 expression. TCF1(hi) cells upregulated α4β7 before gut migration and failed to express cytokines. These cells exhibited recall potential and plasticity following secondary transplantation, including cytokine or Foxp3 expression, but reduced T cell factor 1 (TCF1). Thus, scRNA-Seq suggested divergence of alloreactive CD4(+) T cells into quiescent and effector states during gut GVHD exacerbation by donor DC, reflecting putative heterogeneous priming in vivo. These findings, which are potentially the first at a single-cell level during GVHD over time, may assist in examination of T cell differentiation in patients undergoing alloSCT.

Acute gastrointestinal Graft-versus-Host-Disease (GVHD) is a primary determinant of mortality after allogeneic hematopoietic stem-cell transplantation (alloSCT). It is mediated by alloreactive donor CD4+ T cells that differentiate into pathogenic subsets expressing IFNγ, IL-17A or GM-CSF, and is regulated by subsets expressing IL-10 and/or Foxp3. Developmental relationships between T-helper states during priming in mesenteric lymph nodes (mLN) and effector function in the GI tract remain undefined at genome-scale. We used scRNA-seq and computational modelling to create an atlas of putative differentiation pathways during GVHD. Computational trajectory inference suggested emergence of pathogenic and regulatory states along a single developmental trajectory in mLN. Importantly, we identified an unexpected second trajectory, categorised by little proliferation or cytokine expression, reduced glycolysis, and high TCF1 expression. TCF1hi cells upregulated α4β7 prior to gut migration and failed to express cytokines therein. Nevertheless, they demonstrated recall potential and plasticity following secondary transplantation, including cytokine or Foxp3 expression, but reduced TCF1. Thus, scRNA-seq revealed divergence of allo-reactive CD4+ T cells into quiescent and effector states during gut GVHD, reflecting putative heterogenous priming in vivo. These findings, the first at a single-cell level during GVHD over time, can now be used to interrogate T cell differentiation in patients undergoing alloSCT.

A variety of methods have been developed to demultiplex pooled samples in a single cell RNA sequencing (scRNA-seq) experiment which either require hashtag barcodes or sample genotypes prior to pooling. We introduce scSplit which utilizes genetic differences inferred from scRNA-seq data alone to demultiplex pooled samples. scSplit also enables mapping clusters to original samples. Using simulated, merged, and pooled multi-individual datasets, we show that scSplit prediction is highly concordant with demuxlet predictions and is highly consistent with the known truth in cell-hashing dataset. scSplit is ideally suited to samples without external genotype information and is available at: https://github.com/jon-xu/scSplit

Background Recent developments in single-cell RNA sequencing (scRNA-seq) platforms have vastly increased the number of cells typically assayed in an experiment. Analysis of scRNA-seq data is multidisciplinary in nature, requiring careful consideration of the application of statistical methods with respect to the underlying biology. Few analysis packages exist that are at once robust, are computationally fast, and allow flexible integration with other bioinformatics tools and methods. Findings ascend is an R package comprising tools designed to simplify and streamline the preliminary analysis of scRNA-seq data, while addressing the statistical challenges of scRNA-seq analysis and enabling flexible integration with genomics packages and native R functions, including fast parallel computation and efficient memory management. The package incorporates both novel and established methods to provide a framework to perform cell and gene filtering, quality control, normalization, dimension reduction, clustering, differential expression, and a wide range of visualization functions. Conclusions ascend is designed to work with scRNA-seq data generated by any high-throughput platform and includes functions to convert data objects between software packages. The ascend workflow is simple and interactive, as well as suitable for implementation by a broad range of users, including those with little programming experience.

The libraries generated by high-throughput single cell RNA-sequencing (scRNA-seq) platforms such as the Chromium from 10× Genomics require considerable amounts of sequencing, typically due to the large number of cells. The ability to use these data to address biological questions is directly impacted by the quality of the sequence data. Here we have compared the performance of the Illumina NextSeq 500 and NovaSeq 6000 against the BGI MGISEQ-2000 platform using identical Single Cell 3′ libraries consisting of over 70 000 cells generated on the 10× Genomics Chromium platform. Our results demonstrate a highly comparable performance between the NovaSeq 6000 and MGISEQ-2000 in sequencing quality, and the detection of genes, cell barcodes, Unique Molecular Identifiers. The performance of the NextSeq 500 was also similarly comparable to the MGISEQ-2000 based on the same metrics. Data generated by both sequencing platforms yielded similar analytical outcomes for general single-cell analysis. The performance of the NextSeq 500 and MGISEQ-2000 were also comparable for the deconvolution of multiplexed cell pools via variant calling, and detection of guide RNA (gRNA) from a pooled CRISPR single-cell screen. Our study provides a benchmark for high-capacity sequencing platforms applied to high-throughput scRNA-seq libraries.

The identification of genetic variation with next-generation sequencing is confounded by the complexity of the human genome sequence and by biases that arise during library preparation, sequencing and analysis. We have developed a set of synthetic DNA standards, termed 'sequins', that emulate human genetic features and constitute qualitative and quantitative spike-in controls for genome sequencing. Sequencing reads derived from sequins align exclusively to an artificial in silico reference chromosome, rather than the human reference genome, which allows them them to be partitioned for parallel analysis. Here we use this approach to represent common and clinically relevant genetic variation, ranging from single nucleotide variants to large structural rearrangements and copy-number variation. We validate the design and performance of sequin standards by comparison to examples in the NA12878 reference genome, and we demonstrate their utility during the detection and quantification of variants. We provide sequins as a standardized, quantitative resource against which human genetic variation can be measured and diagnostic performance assessed.

Summary: We assessed the pluripotency of human induced pluripotent stem cells (iPSCs) maintained on an automated platform using StemFlex and TeSR-E8 media. Analysis of transcriptome of single cells revealed similar expression of core pluripotency genes, as well as genes associated with naive and primed states of pluripotency. Analysis of individual cells from four samples consisting of two different iPSC lines each grown in the two culture media revealed a shared subpopulation structure with three main subpopulations different in pluripotency states. By implementing a machine learning approach, we estimated that most cells within each subpopulation are very similar between all four samples. The single-cell RNA sequencing analysis of iPSC lines grown in both media reports the molecular signature in StemFlex medium and how it compares to that observed in the TeSR-E8 medium. : Stem Cells Research; Transcriptomics; Automation Subject Areas: Stem Cells Research, Transcriptomics, Automation

Persistent human papillomavirus (HPV) infection is responsible for at least 5% of human malignancies. Most HPV-associated cancers are initiated by the HPV16 genotype, as confirmed by detection of integrated HPV DNA in cells of oral and anogenital epithelial cancers. However, single-cell RNA sequencing may enable prediction of HPV involvement in carcinogenesis at other sites. We conducted single-cell RNA sequencing on keratinocytes from a mouse transgenic for the E7 gene of HPV16 and showed sensitive and specific detection of HPV16-E7 mRNA, predominantly in basal keratinocytes. We showed that increased E7 mRNA copy number per cell was associated with increased expression of E7 induced genes. This technique enhances detection of active viral transcription in solid tissue and may clarify possible linkage of HPV infection to development of squamous cell carcinoma.

We assessed the pluripotency of human induced pluripotent stem cells (iPSCs) maintained on an automated platform using StemFlexTM and TeSRTM-E8TM media. Analysis of transcriptome of single cells revealed similar expression of core pluripotency genes, as well as genes associated with naive and primed states of pluripotency. Analysis of individual cells from four samples consisting of two different iPSC lines each grown in the two culture media revealed a shared subpopulation structure with three main subpopulations different in pluripotency states. By implementing a machine-learning approach, we estimated that most cells within each subpopulation are very similar between all four samples. The single cell RNA sequencing analysis of iPSC lines grown in both media is the first report of the molecular signature of StemFlex medium and how it compares to that observed in TeSR-E8 medium.

Flow cytometry based electrostatic cell sorting is an important tool in the separation of cell populations. Existing instruments can sort single cells into multi-well collection plates, and keep track of cell of origin and sorted well location. However currently single sorted cell results reflect the population distribution and fail to capture the population diversity. Software was designed that implements a novel sorting approach, "Slice and Dice Sorting," that links a graphical representation of a multi-well plate to logic that ensures that single cells are sampled and sorted from all areas defined by the sort region/s. Therefore the diversity of the total population is captured, and the more frequently occurring or rarer cell types are all sampled. The sorting approach was tested computationally, and using functional cell based assays. Computationally we demonstrate that conventional single cell sorting can sample as little as 50% of the population diversity dependant on the population distribution, and that Slice and Dice sorting samples much more of the variety present within a cell population. We then show by sorting single cells into wells using the Slice and Dice sorting method that there are cells sorted using this method that would be either rarely sorted, or not sorted at all using conventional single cell sorting approaches. The present study demonstrates a novel single cell sorting method that samples much more of the population diversity than current methods. It has implications in clonal selection, stem cell sorting, single cell sequencing and any areas where population heterogeneity is of importance.

Summaryascend is an R package comprised of fast, streamlined analysis functions optimized to address the statistical challenges of single cell RNA-seq. The package incorporates novel and established methods to provide a flexible framework to perform filtering, quality control, normalization, dimension reduction, clustering, differential expression and a wide-range of plotting. ascend is designed to work with scRNA-seq data generated by any high-throughput platform, and includes functions to convert data objects between software packages.AvailabilityThe R package and associated vignettes are freely available at https://github.com/IMB-Computational-Genomics-Lab/ascend.Contactjoseph.powell@uq.edu.auSupplementary informationAn example dataset is available at ArrayExpress, accession number E-MTAB-6108

Absolute neutrophil count (ANC) is used clinically to monitor physiological dysfunctions such as myelosuppression or infection. In the research laboratory, ANC is a valuable measure to monitor the evolution of a wide range of disease states in disease models. Flow cytometry (FCM) is a fast, widely used approach to confidently identify thousands of cells within minutes. FCM can be optimised for absolute counting using spiked-in beads or by measuring the sample volume analysed. Here we combine the 1A8 antibody, specific for the mouse granulocyte protein Ly6G, with flow cytometric counting in straightforward FCM assays for mouse ANC, easily implementable in the research laboratory. Volumetric and Trucount™ bead assays were optimized for mouse neutrophils, and ANC values obtained with these protocols were compared to ANC measured by a dual-platform assay using the Orphee Mythic 18 veterinary haematology analyser. The single platform assays were more precise with decreased intra-assay variability compared with ANC obtained using the dual protocol. Defining ANC based on Ly6G expression produces a 15% higher estimate than the dual protocol. Allowing for this difference in ANC definition, the flow cytometry counting assays using Ly6G can be used reliably in the research laboratory to quantify mouse ANC from a small volume of blood. We demonstrate the utility of the volumetric protocol in a time-course study of chemotherapy induced neutropenia using four drug regimens.

RNA sequencing (RNA-seq) can be used to assemble spliced isoforms, quantify expressed genes and provide a global profile of the transcriptome. However, the size and diversity of the transcriptome, the wide dynamic range in gene expression and inherent technical biases confound RNA-seq analysis. We have developed a set of spike-in RNA standards, termed 'sequins' (sequencing spike-ins), that represent full-length spliced mRNA isoforms. Sequins have an entirely artificial sequence with no homology to natural reference genomes, but they align to gene loci encoded on an artificial in silico chromosome. The combination of multiple sequins across a range of concentrations emulates alternative splicing and differential gene expression, and it provides scaling factors for normalization between samples. We demonstrate the use of sequins in RNA-seq experiments to measure sample-specific biases and determine the limits of reliable transcript assembly and quantification in accompanying human RNA samples. In addition, we have designed a complementary set of sequins that represent fusion genes arising from rearrangements of the in silico chromosome to aid in cancer diagnosis. RNA sequins provide a qualitative and quantitative reference with which to navigate the complexity of the human transcriptome.

Flow cytometry based electrostatic cell sorting is an important tool in the separation of cell populations. Existing instruments can sort single cells into multi-well collection plates, and keep track of cell of origin and sorted well location. However currently single sorted cell results reflect the population distribution and fail to capture the population diversity. Software was designed that implements a novel sorting approach, "Slice and Dice Sorting," that links a graphical representation of a multi-well plate to logic that ensures that single cells are sampled and sorted from all areas defined by the sort region/s. Therefore the diversity of the total population is captured, and the more frequently occurring or rarer cell types are all sampled. The sorting approach was tested computationally, and using functional cell based assays. Computationally we demonstrate that conventional single cell sorting can sample as little as 50% of the population diversity dependant on the population distribution, and that Slice and Dice sorting samples much more of the variety present within a cell population. We then show by sorting single cells into wells using the Slice and Dice sorting method that there are cells sorted using this method that would be either rarely sorted, or not sorted at all using conventional single cell sorting approaches. The present study demonstrates a novel single cell sorting method that samples much more of the population diversity than current methods. It has implications in clonal selection, stem cell sorting, single cell sequencing and any areas where population heterogeneity is of importance.

During the splicing reaction, the 5′ intron end is joined to the branchpoint nucleotide, selecting the next exon to incorporate into the mature RNA and forming an intron lariat, which is excised. Despite a critical role in gene splicing, the locations and features of human splicing branchpoints are largely unknown. We use exoribonuclease digestion and targeted RNA-sequencing to enrich for sequences that traverse the lariat junction and, by split and inverted alignment, reveal the branchpoint. We identify 59,359 high-confidence human branchpoints in >10,000 genes, providing a first map of splicing branchpoints in the human genome. Branchpoints are predominantly adenosine, highly conserved, and closely distributed to the 3′ splice site. Analysis of human branchpoints reveals numerous novel features, including distinct features of branchpoints for alternatively spliced exons and a family of conserved sequence motifs overlapping branchpoints we term B-boxes, which exhibit maximal nucleotide diversity while maintaining interactions with the keto-rich U2 snRNA. Different B-box motifs exhibit divergent usage in vertebrate lineages and associate with other splicing elements and distinct intron–exon architectures, suggesting integration within a broader regulatory splicing code. Lastly, although branchpoints are refractory to common mutational processes and genetic variation, mutations occurring at branchpoint nucleotides are enriched for disease associations.