Enhancer retargeting of CDX2 and UBTF::ATXN7L3 define a subtype of high-risk B-progenitor acute lymphoblastic leukemia.

Transcriptome sequencing has identified multiple subtypes of B-progenitor acute lymphoblastic leukemia (B-ALL) of prognostic significance, but a minority of cases lack a known genetic driver. Here, we used integrated whole-genome (WGS) and -transcriptome sequencing (RNA-seq), enhancer mapping, and chromatin topology analysis to identify previously unrecognized genomic drivers in B-ALL. Newly diagnosed (n = 3221) and relapsed (n = 177) B-ALL cases with tumor RNA-seq were studied. WGS was performed to detect mutations, structural variants, and copy number alterations. Integrated analysis of histone 3 lysine 27 acetylation and chromatin looping was performed using HiChIP. We identified a subset of 17 newly diagnosed and 5 relapsed B-ALL cases with a distinct gene expression profile and 2 universal and unique genomic alterations resulting from aberrant recombination-activating gene activation: a focal deletion downstream of PAN3 at 13q12.2 resulting in CDX2 deregulation by the PAN3 enhancer and a focal deletion of exons 18-21 of UBTF at 17q21.31 resulting in a chimeric fusion, UBTF::ATXN7L3. A subset of cases also had rearrangement and increased expression of the PAX5 gene, which is otherwise uncommon in B-ALL. Patients were more commonly female and young adult with median age 35 (range,12-70 years). The immunophenotype was characterized by CD10 negativity and immunoglobulin M positivity. Among 16 patients with known clinical response, 9 (56.3%) had high-risk features including relapse (n = 4) or minimal residual disease >1% at the end of remission induction (n = 5). CDX2-deregulated, UBTF::ATXN7L3 rearranged (CDX2/UBTF) B-ALL is a high-risk subtype of leukemia in young adults for which novel therapeutic approaches are required.

Human yolk sac-like haematopoiesis generates RUNX1-, GFI1- and/or GFI1B-dependent blood and SOX17-positive endothelium.

The genetic regulatory network controlling early fate choices during human blood cell development are not well understood. We used human pluripotent stem cell reporter lines to track the development of endothelial and haematopoietic populations in an in vitro model of human yolk-sac development. We identified SOX17-CD34+CD43- endothelial cells at day 2 of blast colony development, as a haemangioblast-like branch point from which SOX17-CD34+CD43+ blood cells and SOX17+CD34+CD43- endothelium subsequently arose. Most human blood cell development was dependent on RUNX1. Deletion of RUNX1 only permitted a single wave of yolk sac-like primitive erythropoiesis, but no yolk sac myelopoiesis or aorta-gonad-mesonephros (AGM)-like haematopoiesis. Blocking GFI1 and/or GFI1B activity with a small molecule inhibitor abrogated all blood cell development, even in cell lines with an intact RUNX1 gene. Together, our data define the hierarchical requirements for RUNX1, GFI1 and/or GFI1B during early human haematopoiesis arising from a yolk sac-like SOX17-negative haemogenic endothelial intermediate.

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.

STRipy: A graphical application for enhanced genotyping of pathogenic short tandem repeats in sequencing data.

Expansions of short tandem repeats (STRs) have been implicated as the causal variant in over 50 diseases known to date. There are several tools which can genotype STRs from high-throughput sequencing (HTS) data. However, running these tools out of the box only allows around half of the known disease-causing loci to be genotyped. Furthermore, the genotypes estimated at these loci are often underestimated with maximum lengths limited to either the read or fragment length, which is less than the pathogenic cutoff for some diseases. Although analysis tools can be customized to genotype extra loci, this requires proficiency in bioinformatics to set up, limiting their widespread usage by other researchers and clinicians. To address these issues, we have developed a new software called STRipy, which is able to target all known disease-causing STRs from HTS data. We created an intuitive graphical interface for STRipy and significantly simplified the detection of STRs expansions. Moreover, we genotyped all disease loci for over two and half thousand samples to provide population-wide distributions to assist with interpretation of results. We believe the simplicity and breadth of STRipy will increase the genotyping of STRs in sequencing data resulting in further diagnoses of rare STR diseases.

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.

Kidney micro-organoids in suspension culture as a scalable source of human pluripotent stem cell-derived kidney cells.

Kidney organoids have potential uses in disease modelling, drug screening and regenerative medicine. However, novel cost-effective techniques are needed to enable scaled-up production of kidney cell types in vitro We describe here a modified suspension culture method for the generation of kidney micro-organoids from human pluripotent stem cells. Optimisation of differentiation conditions allowed the formation of micro-organoids, each containing six to ten nephrons that were surrounded by endothelial and stromal populations. Single cell transcriptional profiling confirmed the presence and transcriptional equivalence of all anticipated renal cell types consistent with a previous organoid culture method. This suspension culture micro-organoid methodology resulted in a three- to fourfold increase in final cell yield compared with static culture, thereby representing an economical approach to the production of kidney cells for various biological applications.

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.

Detecting copy number alterations in RNA-Seq using SuperFreq.

MOTIVATION: Calling copy number alterations (CNAs) from RNA sequencing (RNA-Seq) is challenging, because of the marked variability in coverage across genes and paucity of single nucleotide polymorphisms (SNPs). We have adapted SuperFreq to call absolute and allele sensitive CNAs from RNA-Seq. SuperFreq uses an error-propagation framework to combine and maximize information from read counts and B-allele frequencies. RESULTS: We used datasets from The Cancer Genome Atlas (TCGA) to assess the validity of CNA calls from RNA-Seq. When ploidy estimates were consistent, we found agreement with DNA SNP-arrays for over 98% of the genome for acute myeloid leukaemia (TCGA-AML, n = 116) and 87% for colorectal cancer (TCGA-CRC, n = 377). The sensitivity of CNA calling from RNA-Seq was dependent on gene density. Using RNA-Seq, SuperFreq detected 78% of CNA calls covering 100 or more genes with a precision of 94%. Recall dropped for focal events, but this also depended on signal intensity. For example, in the CRC cohort SuperFreq identified all cases (7/7) with high-level amplification of ERBB2, where the copy number was typically >20, but identified only 6% of cases (1/17) with moderate amplification of IGF2, which occurs over a smaller interval. SuperFreq offers an integrated platform for identification of CNAs and point mutations. As evidence of how SuperFreq can be applied, we used it to reproduce the established relationship between somatic mutation load and CNA profile in CRC using RNA-Seq alone. AVAILABILITY AND IMPLEMENTATION: SuperFreq is implemented in R and the code is available through GitHub: https://github.com/ChristofferFlensburg/SuperFreq/. Data and code to reproduce the figures are available at: https://gitlab.wehi.edu.au/flensburg.c/SuperFreq_RNA_paper. Data from TCGA (phs000178) was accessed from GDC following completion of a data access request through the database of Genotypes and Phenotypes (dbGaP). Data from the Leucegene consortium was downloaded from GEO (AML samples: GSE67040; normal CD34+ cells: GSE48846). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

RNA-Seq of amniotic fluid cell-free RNA: a discovery phase study of the pathophysiology of congenital cytomegalovirus infection.

BACKGROUND: Congenital cytomegalovirus infection is the most common perinatal infection and a significant cause of sensorineural hearing loss, cerebral palsy, and neurodevelopmental disability. There is a paucity of human gene expression studies examining the pathophysiology of cytomegalovirus infection. OBJECTIVE: This study aimed to perform a whole transcriptomic assessment of amniotic fluid from pregnancies with live fetuses to identify differentially expressed genes and enriched Gene Ontology categories associated with congenital cytomegalovirus infection. STUDY DESIGN: Amniotic fluid supernatant was prospectively collected from pregnant women undergoing amniocentesis for suspected congenital cytomegalovirus infection because of first-trimester maternal primary infection or ultrasound features suggestive of fetal infection. Women who had received therapy to prevent fetal infection were excluded. Congenital cytomegalovirus infection was diagnosed via viral polymerase chain reaction of amniotic fluid; cytomegalovirus-infected fetuses were paired with noninfected controls, matched for gestational age and fetal sex. Paired-end RNA sequencing was performed on amniotic fluid cell-free RNA with the Novaseq 6000 at a depth of 30 million reads per sample. Following quality control and filtering, reads were mapped to the human genome and counts summarized across genes. Differentially expressed genes were identified using 2 approaches: voomWithQualityWeights in conjunction with limma and RUVSeq with edgeR. Genes with a false discovery rate <0.05 were considered statistically significant. Differential exon use was analyzed using DEXSeq. Functional analysis was performed using gene set enrichment analysis and Ingenuity Pathway Analysis. Manual curation of differentially regulated genes was also performed. RESULTS: Amniotic fluid samples were collected from 50 women; 16 (32%) had congenital cytomegalovirus infection confirmed by polymerase chain reaction. After excluding 3 samples without matched controls, 13 cytomegalovirus-infected samples collected at 18 to 23 weeks and 13 cytomegalovirus-negative gestation-matched controls were submitted for RNA sequencing and analysis (N=26). Ten of the 13 pregnancies with cytomegalovirus-infected fetuses had amniocentesis because of serologic evidence of maternal primary infection with normal fetal ultrasound, and 3 had amniocentesis because of ultrasound abnormality suggestive of cytomegalovirus infection. Four cytomegalovirus-infected pregnancies ended in termination (n=3) or fetal death (n=1), and 9 resulted in live births. Pregnancy outcomes were available for 11 of the 13 cytomegalovirus-negative controls; all resulted in live births of clinically-well infants. Differential gene expression analysis revealed 309 up-regulated and 32 down-regulated genes in the cytomegalovirus-infected group compared with the cytomegalovirus-negative group. Gene set enrichment analysis showed significant enrichment of multiple Gene Ontology categories involving the innate immune response to viral infection and interferon signaling. Of the 32 significantly down-regulated genes, 8 were known to be involved in neurodevelopment and preferentially expressed by the brain. Six specific cellular restriction factors involved in host defense to cytomegalovirus infection were up-regulated in the cytomegalovirus-infected group. Ingenuity Pathway Analysis predicted the activation of pathways involved in progressive neurologic disease and inflammatory neurologic disease. CONCLUSION: In this next-generation sequencing study, we revealed new insights into the pathophysiology of congenital cytomegalovirus infection. These data on the up-regulation of the intraamniotic innate immune response to cytomegalovirus infection and the dysregulation of neurodevelopmental genes may inform future approaches to developing prognostic markers and assessing fetal responses to in utero therapy.

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.

SuperFreq: Integrated mutation detection and clonal tracking in cancer.

Analysing multiple cancer samples from an individual patient can provide insight into the way the disease evolves. Monitoring the expansion and contraction of distinct clones helps to reveal the mutations that initiate the disease and those that drive progression. Existing approaches for clonal tracking from sequencing data typically require the user to combine multiple tools that are not purpose-built for this task. Furthermore, most methods require a matched normal (non-tumour) sample, which limits the scope of application. We developed SuperFreq, a cancer exome sequencing analysis pipeline that integrates identification of somatic single nucleotide variants (SNVs) and copy number alterations (CNAs) and clonal tracking for both. SuperFreq does not require a matched normal and instead relies on unrelated controls. When analysing multiple samples from a single patient, SuperFreq cross checks variant calls to improve clonal tracking, which helps to separate somatic from germline variants, and to resolve overlapping CNA calls. To demonstrate our software we analysed 304 cancer-normal exome samples across 33 cancer types in The Cancer Genome Atlas (TCGA) and evaluated the quality of the SNV and CNA calls. We simulated clonal evolution through in silico mixing of cancer and normal samples in known proportion. We found that SuperFreq identified 93% of clones with a cellular fraction of at least 50% and mutations were assigned to the correct clone with high recall and precision. In addition, SuperFreq maintained a similar level of performance for most aspects of the analysis when run without a matched normal. SuperFreq is highly versatile and can be applied in many different experimental settings for the analysis of exomes and other capture libraries. We demonstrate an application of SuperFreq to leukaemia patients with diagnosis and relapse samples.

Direct reprogramming to human nephron progenitor-like cells using inducible piggyBac transposon expression of SNAI2-EYA1-SIX1.

All nephrons in the mammalian kidney arise from a transient nephron progenitor population that is lost close to the time of birth. The generation of new nephron progenitors and their maintenance in culture are central to the success of kidney regenerative strategies. Using a lentiviral screening approach, we previously generated a human induced nephron progenitor-like state in vitro using a pool of six transcription factors. Here, we sought to develop a more efficient approach for direct reprogramming of human cells that could be applied in vivo. PiggyBac transposons are a non-viral integrating gene delivery system that is suitable for in vivo use and allows for simultaneous delivery of multiple genes. Using an inducible piggyBac transposon system, we optimized a protocol for the direct reprogramming of HK2 cells to induced nephron progenitor-like cells with expression of only 3 transcription factors (SNAI2, EYA1, and SIX1). Culture in conditions supportive of the nephron progenitor state further increased the expression of nephron progenitor genes. The refined protocol was then applied to primary human renal epithelial cells, which integrated into developing nephron structures in vitro and in vivo. Such inducible reprogramming to nephron progenitor-like cells could facilitate direct cellular reprogramming for kidney regeneration.

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.

Targeted therapy and disease monitoring in CNTRL-FGFR1-driven leukaemia.

We report two patients with leukaemia driven by the rare CNTRL-FGFR1 fusion oncogene. This fusion arises from a t(8;9)(p12;q33) translocation, and is a rare driver of biphenotypic leukaemia in children. We used RNA sequencing to report novel features of expressed CNTRL-FGFR1, including CNTRL-FGFR1 fusion alternative splicing. From this knowledge, we designed and tested a Droplet Digital PCR assay that detects CNTRL-FGFR1 expression to approximately one cell in 100 000 using fusion breakpoint-specific primers and probes. We also utilised cell-line models to show that effective tyrosine kinase inhibitors, which may be included in treatment regimens for this disease, are only those that block FGFR1 phosphorylation.

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.

MCM3AP in recessive Charcot-Marie-Tooth neuropathy and mild intellectual disability.

Defects in mRNA export from the nucleus have been linked to various neurodegenerative disorders. We report mutations in the gene MCM3AP, encoding the germinal center associated nuclear protein (GANP), in nine affected individuals from five unrelated families. The variants were associated with severe childhood onset primarily axonal (four families) or demyelinating (one family) Charcot-Marie-Tooth neuropathy. Mild to moderate intellectual disability was present in seven of nine affected individuals. The affected individuals were either compound heterozygous or homozygous for different MCM3AP variants, which were predicted to cause depletion of GANP or affect conserved amino acids with likely importance for its function. Accordingly, fibroblasts of affected individuals from one family demonstrated severe depletion of GANP. GANP has been described to function as an mRNA export factor, and to suppress TDP-43-mediated motor neuron degeneration in flies. Thus our results suggest defective mRNA export from nucleus as a potential pathogenic mechanism of axonal degeneration in these patients. The identification of MCM3AP variants in affected individuals from multiple centres establishes it as a disease gene for childhood-onset recessively inherited Charcot-Marie-Tooth neuropathy with intellectual disability.

Transcriptional profiles for distinct aggregation states of mutant Huntingtin exon 1 protein unmask new Huntington's disease pathways.

Huntington's disease is caused by polyglutamine (polyQ)-expansion mutations in the CAG tandem repeat of the Huntingtin gene. The central feature of Huntington's disease pathology is the aggregation of mutant Huntingtin (Htt) protein into micrometer-sized inclusion bodies. Soluble mutant Htt states are most proteotoxic and trigger an enhanced risk of death whereas inclusions confer different changes to cellular health, and may even provide adaptive responses to stress. Yet the molecular mechanisms underpinning these changes remain unclear. Using the flow cytometry method of pulse-shape analysis (PulSA) to sort neuroblastoma (Neuro2a) cells enriched with mutant or wild-type Htt into different aggregation states, we clarified which transcriptional signatures were specifically attributable to cells before versus after inclusion assembly. Dampened CREB signalling was the most striking change overall and invoked specifically by soluble mutant Httex1 states. Toxicity could be rescued by stimulation of CREB signalling. Other biological processes mapped to different changes before and after aggregation included NF-kB signalling, autophagy, SUMOylation, transcription regulation by histone deacetylases and BRD4, NAD+ biosynthesis, ribosome biogenesis and altered HIF-1 signalling. These findings open the path for therapeutic strategies targeting key molecular changes invoked prior to, and subsequently to, Httex1 aggregation.