propeller: testing for differences in cell type proportions in single cell data.

MOTIVATION: Single cell RNA-Sequencing (scRNA-seq) has rapidly gained popularity over the last few years for profiling the transcriptomes of thousands to millions of single cells. This technology is now being used to analyse experiments with complex designs including biological replication. One question that can be asked from single cell experiments, which has been difficult to directly address with bulk RNA-seq data, is whether the cell type proportions are different between two or more experimental conditions. As well as gene expression changes, the relative depletion or enrichment of a particular cell type can be the functional consequence of disease or treatment. However, cell type proportion estimates from scRNA-seq data are variable and statistical methods that can correctly account for different sources of variability are needed to confidently identify statistically significant shifts in cell type composition between experimental conditions. RESULTS: We have developed propeller, a robust and flexible method that leverages biological replication to find statistically significant differences in cell type proportions between groups. Using simulated cell type proportions data, we show that propeller performs well under a variety of scenarios. We applied propeller to test for significant changes in cell type proportions related to human heart development, ageing and COVID-19 disease severity. AVAILABILITY AND IMPLEMENTATION: The propeller method is publicly available in the open source speckle R package (https://github.com/phipsonlab/speckle). All the analysis code for the article is available at the associated analysis website: https://phipsonlab.github.io/propeller-paper-analysis/. The speckle package, analysis scripts and datasets have been deposited at https://doi.org/10.5281/zenodo.7009042. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Sex-Specific Control of Human Heart Maturation by the Progesterone Receptor.

BACKGROUND: Despite in-depth knowledge of the molecular mechanisms controlling embryonic heart development, little is known about the signals governing postnatal maturation of the human heart. METHODS: Single-nucleus RNA sequencing of 54 140 nuclei from 9 human donors was used to profile transcriptional changes in diverse cardiac cell types during maturation from fetal stages to adulthood. Bulk RNA sequencing and the Assay for Transposase-Accessible Chromatin using sequencing were used to further validate transcriptional changes and to profile alterations in the chromatin accessibility landscape in purified cardiomyocyte nuclei from 21 human donors. Functional validation studies of sex steroids implicated in cardiac maturation were performed in human pluripotent stem cell-derived cardiac organoids and mice. RESULTS: Our data identify the progesterone receptor as a key mediator of sex-dependent transcriptional programs during cardiomyocyte maturation. Functional validation studies in human cardiac organoids and mice demonstrate that the progesterone receptor drives sex-specific metabolic programs and maturation of cardiac contractile properties. CONCLUSIONS: These data provide a blueprint for understanding human heart maturation in both sexes and reveal an important role for the progesterone receptor in human heart development.

Single cell analysis of the developing mouse kidney provides deeper insight into marker gene expression and ligand-receptor crosstalk.

Recent advances in the generation of kidney organoids and the culture of primary nephron progenitors from mouse and human have been based on knowledge of the molecular basis of kidney development in mice. Although gene expression during kidney development has been intensely investigated, single cell profiling provides new opportunities to further subsect component cell types and the signalling networks at play. Here, we describe the generation and analysis of 6732 single cell transcriptomes from the fetal mouse kidney [embryonic day (E)18.5] and 7853 sorted nephron progenitor cells (E14.5). These datasets provide improved resolution of cell types and specific markers, including subdivision of the renal stroma and heterogeneity within the nephron progenitor population. Ligand-receptor interaction and pathway analysis reveals novel crosstalk between cellular compartments and associates new pathways with differentiation of nephron and ureteric epithelium cell types. We identify transcriptional congruence between the distal nephron and ureteric epithelium, showing that most markers previously used to identify ureteric epithelium are not specific. Together, this work improves our understanding of metanephric kidney development and provides a template to guide the regeneration of renal tissue.

Evaluation of variability in human kidney organoids.

The utility of human pluripotent stem cell-derived kidney organoids relies implicitly on the robustness and transferability of the protocol. Here we analyze the sources of transcriptional variation in a specific kidney organoid protocol. Although individual organoids within a differentiation batch showed strong transcriptional correlation, we noted significant variation between experimental batches, particularly in genes associated with temporal maturation. Single-cell profiling revealed shifts in nephron patterning and proportions of component cells. Distinct induced pluripotent stem cell clones showed congruent transcriptional programs, with interexperimental and interclonal variation also strongly associated with nephron patterning. Epithelial cells isolated from organoids aligned with total organoids at the same day of differentiation, again implicating relative maturation as a confounder. This understanding of experimental variation facilitated an optimized analysis of organoid-based disease modeling, thereby increasing the utility of kidney organoids for personalized medicine and functional genomics.

Patient-iPSC-Derived Kidney Organoids Show Functional Validation of a Ciliopathic Renal Phenotype and Reveal Underlying Pathogenetic Mechanisms.

Despite the increasing diagnostic rate of genomic sequencing, the genetic basis of more than 50% of heritable kidney disease remains unresolved. Kidney organoids differentiated from induced pluripotent stem cells (iPSCs) of individuals affected by inherited renal disease represent a potential, but unvalidated, platform for the functional validation of novel gene variants and investigation of underlying pathogenetic mechanisms. In this study, trio whole-exome sequencing of a prospectively identified nephronophthisis (NPHP) proband and her parents identified compound-heterozygous variants in IFT140, a gene previously associated with NPHP-related ciliopathies. IFT140 plays a key role in retrograde intraflagellar transport, but the precise downstream cellular mechanisms responsible for disease presentation remain unknown. A one-step reprogramming and gene-editing protocol was used to derive both uncorrected proband iPSCs and isogenic gene-corrected iPSCs, which were differentiated to kidney organoids. Proband organoid tubules demonstrated shortened, club-shaped primary cilia, whereas gene correction rescued this phenotype. Differential expression analysis of epithelial cells isolated from organoids suggested downregulation of genes associated with apicobasal polarity, cell-cell junctions, and dynein motor assembly in proband epithelial cells. Matrigel cyst cultures confirmed a polarization defect in proband versus gene-corrected renal epithelium. As such, this study represents a "proof of concept" for using proband-derived iPSCs to model renal disease and illustrates dysfunctional cellular pathways beyond the primary cilium in the setting of IFT140 mutations, which are established for other NPHP genotypes.

Exploring the single-cell RNA-seq analysis landscape with the scRNA-tools database.

As single-cell RNA-sequencing (scRNA-seq) datasets have become more widespread the number of tools designed to analyse these data has dramatically increased. Navigating the vast sea of tools now available is becoming increasingly challenging for researchers. In order to better facilitate selection of appropriate analysis tools we have created the scRNA-tools database (www.scRNA-tools.org) to catalogue and curate analysis tools as they become available. Our database collects a range of information on each scRNA-seq analysis tool and categorises them according to the analysis tasks they perform. Exploration of this database gives insights into the areas of rapid development of analysis methods for scRNA-seq data. We see that many tools perform tasks specific to scRNA-seq analysis, particularly clustering and ordering of cells. We also find that the scRNA-seq community embraces an open-source and open-science approach, with most tools available under open-source licenses and preprints being extensively used as a means to describe methods. The scRNA-tools database provides a valuable resource for researchers embarking on scRNA-seq analysis and records the growth of the field over time.

Haploinsufficiency for the Six2 gene increases nephron progenitor proliferation promoting branching and nephron number.

The regulation of final nephron number in the kidney is poorly understood. Cessation of nephron formation occurs when the self-renewing nephron progenitor population commits to differentiation. Transcription factors within this progenitor population, such as SIX2, are assumed to control expression of genes promoting self-renewal such that homozygous Six2 deletion results in premature commitment and an early halt to kidney development. In contrast, Six2 heterozygotes were assumed to be unaffected. Using quantitative morphometry, we found a paradoxical 18% increase in ureteric branching and final nephron number in Six2 heterozygotes, despite evidence for reduced levels of SIX2 protein and transcript. This was accompanied by a clear shift in nephron progenitor identity with a distinct subset of downregulated progenitor genes such as Cited1 and Meox1 while other genes were unaffected. The net result was an increase in nephron progenitor proliferation, as assessed by elevated EdU (5-ethynyl-2'-deoxyuridine) labeling, an increase in MYC protein, and transcriptional upregulation of MYC target genes. Heterozygosity for Six2 on an Fgf20-/- background resulted in premature differentiation of the progenitor population, confirming that progenitor regulation is compromised in Six2 heterozygotes. Overall, our studies reveal a unique dose response of nephron progenitors to the level of SIX2 protein in which the role of SIX2 in progenitor proliferation versus self-renewal is separable.

MOZ and BMI1 play opposing roles during Hox gene activation in ES cells and in body segment identity specification in vivo.

Hox genes underlie the specification of body segment identity in the anterior-posterior axis. They are activated during gastrulation and undergo a dynamic shift from a transcriptionally repressed to an active chromatin state in a sequence that reflects their chromosomal location. Nevertheless, the precise role of chromatin modifying complexes during the initial activation phase remains unclear. In the current study, we examined the role of chromatin regulators during Hox gene activation. Using embryonic stem cell lines lacking the transcriptional activator MOZ and the polycomb-family repressor BMI1, we showed that MOZ and BMI1, respectively, promoted and repressed Hox genes during the shift from the transcriptionally repressed to the active state. Strikingly however, MOZ but not BMI1 was required to regulate Hox mRNA levels after the initial activation phase. To determine the interaction of MOZ and BMI1 in vivo, we interrogated their role in regulating Hox genes and body segment identity using Moz;Bmi1 double deficient mice. We found that the homeotic transformations and shifts in Hox gene expression boundaries observed in single Moz and Bmi1 mutant mice were rescued to a wild type identity in Moz;Bmi1 double knockout animals. Together, our findings establish that MOZ and BMI1 play opposing roles during the onset of Hox gene expression in the ES cell model and during body segment identity specification in vivo. We propose that chromatin-modifying complexes have a previously unappreciated role during the initiation phase of Hox gene expression, which is critical for the correct specification of body segment identity.

Early lineage priming by trisomy of Erg leads to myeloproliferation in a Down syndrome model.

Down syndrome (DS), with trisomy of chromosome 21 (HSA21), is the commonest human aneuploidy. Pre-leukemic myeloproliferative changes in DS foetal livers precede the acquisition of GATA1 mutations, transient myeloproliferative disorder (DS-TMD) and acute megakaryocytic leukemia (DS-AMKL). Trisomy of the Erg gene is required for myeloproliferation in the Ts(1716)65Dn DS mouse model. We demonstrate here that genetic changes specifically attributable to trisomy of Erg lead to lineage priming of primitive and early multipotential progenitor cells in Ts(1716)65Dn mice, excess megakaryocyte-erythroid progenitors, and malignant myeloproliferation. Gene expression changes dependent on trisomy of Erg in Ts(1716)65Dn multilineage progenitor cells were correlated with those associated with trisomy of HSA21 in human DS hematopoietic stem and primitive progenitor cells. These data suggest a role for ERG as a regulator of hematopoietic lineage potential, and that trisomy of ERG in the context of DS foetal liver hemopoiesis drives the pre-leukemic changes that predispose to subsequent DS-TMD and DS-AMKL.

missMethyl: an R package for analyzing data from Illumina's HumanMethylation450 platform.

UNLABELLED: DNA methylation is one of the most commonly studied epigenetic modifications due to its role in both disease and development. The Illumina HumanMethylation450 BeadChip is a cost-effective way to profile >450 000 CpGs across the human genome, making it a popular platform for profiling DNA methylation. Here we introduce missMethyl, an R package with a suite of tools for performing normalization, removal of unwanted variation in differential methylation analysis, differential variability testing and gene set analysis for the 450K array. AVAILABILITY AND IMPLEMENTATION: missMethyl is an R package available from the Bioconductor project at www.bioconductor.org. CONTACT: alicia.oshlack@mcri.edu.au SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

limma powers differential expression analyses for RNA-sequencing and microarray studies.

limma is an R/Bioconductor software package that provides an integrated solution for analysing data from gene expression experiments. It contains rich features for handling complex experimental designs and for information borrowing to overcome the problem of small sample sizes. Over the past decade, limma has been a popular choice for gene discovery through differential expression analyses of microarray and high-throughput PCR data. The package contains particularly strong facilities for reading, normalizing and exploring such data. Recently, the capabilities of limma have been significantly expanded in two important directions. First, the package can now perform both differential expression and differential splicing analyses of RNA sequencing (RNA-seq) data. All the downstream analysis tools previously restricted to microarray data are now available for RNA-seq as well. These capabilities allow users to analyse both RNA-seq and microarray data with very similar pipelines. Second, the package is now able to go past the traditional gene-wise expression analyses in a variety of ways, analysing expression profiles in terms of co-regulated sets of genes or in terms of higher-order expression signatures. This provides enhanced possibilities for biological interpretation of gene expression differences. This article reviews the philosophy and design of the limma package, summarizing both new and historical features, with an emphasis on recent enhancements and features that have not been previously described.

Polycomb repressive complex 2 (PRC2) suppresses Eµ-myc lymphoma.

Deregulation of polycomb group complexes polycomb repressive complex 1 (PRC1) and 2 (PRC2) is associated with human cancers. Although inactivating mutations in PRC2-encoding genes EZH2, EED, and SUZ12 are present in T-cell acute lymphoblastic leukemia and in myeloid malignancies, gain-of-function mutations in EZH2 are frequently observed in B-cell lymphoma, implying disease-dependent effects of individual mutations. We show that, in contrast to PRC1, PRC2 is a tumor suppressor in Eµ-myc lymphomagenesis, because disease onset was accelerated by heterozygosity for Suz12 or by short hairpin RNA-mediated knockdown of Suz12 or Ezh2. Accelerated lymphomagenesis was associated with increased accumulation of B-lymphoid cells in the absence of effects on apoptosis or cell cycling. However, Suz12-deficient B-lymphoid progenitors exhibit enhanced serial clonogenicity. Thus, PRC2 normally restricts the self-renewal of B-lymphoid progenitors, the disruption of which contributes to lymphomagenesis. This finding provides new insight regarding the functional contribution of mutations in PRC2 in a range of leukemias.

Bcl-2, Bcl-x(L), and Bcl-w are not equivalent targets of ABT-737 and navitoclax (ABT-263) in lymphoid and leukemic cells.

The BH3-mimetic ABT-737 and an orally bioavailable compound of the same class, navitoclax (ABT-263), have shown promising antitumor efficacy in preclinical and early clinical studies. Although both drugs avidly bind Bcl-2, Bcl-x(L), and Bcl-w in vitro, we find that Bcl-2 is the critical target in vivo, suggesting that patients with tumors overexpressing Bcl-2 will probably benefit. In human non-Hodgkin lymphomas, high expression of Bcl-2 but not Bcl-x(L) predicted sensitivity to ABT-263. Moreover, we show that increasing Bcl-2 sensitized normal and transformed lymphoid cells to ABT-737 by elevating proapoptotic Bim. In striking contrast, increasing Bcl-x(L) or Bcl-w conferred robust resistance to ABT-737, despite also increasing Bim. Cell-based protein redistribution assays unexpectedly revealed that ABT-737 disrupts Bcl-2/Bim complexes more readily than Bcl-x(L)/Bim or Bcl-w/Bim complexes. These results have profound implications for how BH3-mimetics induce apoptosis and how the use of these compounds can be optimized for treating lymphoid malignancies.

SuperCellCyto: enabling efficient analysis of large scale cytometry datasets.

Advancements in cytometry technologies have enabled quantification of up to 50 proteins across millions of cells at single cell resolution. Analysis of cytometry data routinely involves tasks such as data integration, clustering, and dimensionality reduction. While numerous tools exist, many require extensive run times when processing large cytometry data containing millions of cells. Existing solutions, such as random subsampling, are inadequate as they risk excluding rare cell subsets. To address this, we propose SuperCellCyto, an R package that builds on the SuperCell tool which groups highly similar cells into supercells. SuperCellCyto is available on GitHub ( https://github.com/phipsonlab/SuperCellCyto ) and Zenodo ( https://doi.org/10.5281/zenodo.10521294 ).

Library size confounds biology in spatial transcriptomics data.

Spatial molecular data has transformed the study of disease microenvironments, though, larger datasets pose an analytics challenge prompting the direct adoption of single-cell RNA-sequencing tools including normalization methods. Here, we demonstrate that library size is associated with tissue structure and that normalizing these effects out using commonly applied scRNA-seq normalization methods will negatively affect spatial domain identification. Spatial data should not be specifically corrected for library size prior to analysis, and algorithms designed for scRNA-seq data should be adopted with caution.

In situ single-cell profiling sheds light on IFI27 localisation during SARS-CoV-2 infection.

The utilization of single-cell resolved spatial transcriptomics to delineate immune responses during SARS-CoV-2 infection was able to identify M1 macrophages to have elevated expression of IFI27 in areas of infection.

Proteomic and metabolomic analyses of mitochondrial complex I-deficient mouse model generated by spontaneous B2 short interspersed nuclear element (SINE) insertion into NADH dehydrogenase (ubiquinone) Fe-S protein 4 (Ndufs4) gene.

Eukaryotic cells generate energy in the form of ATP, through a network of mitochondrial complexes and electron carriers known as the oxidative phosphorylation system. In mammals, mitochondrial complex I (CI) is the largest component of this system, comprising 45 different subunits encoded by mitochondrial and nuclear DNA. Humans diagnosed with mutations in the gene NDUFS4, encoding a nuclear DNA-encoded subunit of CI (NADH dehydrogenase ubiquinone Fe-S protein 4), typically suffer from Leigh syndrome, a neurodegenerative disease with onset in infancy or early childhood. Mitochondria from NDUFS4 patients usually lack detectable NDUFS4 protein and show a CI stability/assembly defect. Here, we describe a recessive mouse phenotype caused by the insertion of a transposable element into Ndufs4, identified by a novel combined linkage and expression analysis. Designated Ndufs4(fky), the mutation leads to aberrant transcript splicing and absence of NDUFS4 protein in all tissues tested of homozygous mice. Physical and behavioral symptoms displayed by Ndufs4(fky/fky) mice include temporary fur loss, growth retardation, unsteady gait, and abnormal body posture when suspended by the tail. Analysis of CI in Ndufs4(fky/fky) mice using blue native PAGE revealed the presence of a faster migrating crippled complex. This crippled CI was shown to lack subunits of the "N assembly module", which contains the NADH binding site, but contained two assembly factors not present in intact CI. Metabolomic analysis of the blood by tandem mass spectrometry showed increased hydroxyacylcarnitine species, implying that the CI defect leads to an imbalanced NADH/NAD(+) ratio that inhibits mitochondrial fatty acid ß-oxidation.

A recombination hotspot leads to sequence variability within a novel gene (AK005651) and contributes to type 1 diabetes susceptibility.

More than 25 loci have been linked to type 1 diabetes (T1D) in the nonobese diabetic (NOD) mouse, but identification of the underlying genes remains challenging. We describe here the positional cloning of a T1D susceptibility locus, Idd11, located on mouse chromosome 4. Sequence analysis of a series of congenic NOD mouse strains over a critical 6.9-kb interval in these mice and in 25 inbred strains identified several haplotypes, including a unique NOD haplotype, associated with varying levels of T1D susceptibility. Haplotype diversity within this interval between congenic NOD mouse strains was due to a recombination hotspot that generated four crossover breakpoints, including one with a complex conversion tract. The Idd11 haplotype and recombination hotspot are located within a predicted gene of unknown function, which exhibits decreased expression in relevant tissues of NOD mice. Notably, it was the recombination hotspot that aided our mapping of Idd11 and confirms that recombination hotspots can create genetic variation affecting a common polygenic disease. This finding has implications for human genetic association studies, which may be affected by the approximately 33,000 estimated hotspots in the genome.