Author Correction: The neuropeptide VIP confers anticipatory mucosal immunity by regulating ILC3 activity.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The neuropeptide VIP confers anticipatory mucosal immunity by regulating ILC3 activity.

Group 3 innate lymphoid cell (ILC3)-mediated production of the cytokine interleukin-22 (IL-22) is critical for the maintenance of immune homeostasis in the gastrointestinal tract. Here, we find that the function of ILC3s is not constant across the day, but instead oscillates between active phases and resting phases. Coordinate responsiveness of ILC3s in the intestine depended on the food-induced expression of the neuropeptide vasoactive intestinal peptide (VIP). Intestinal ILC3s had high expression of the G protein-coupled receptor vasoactive intestinal peptide receptor 2 (VIPR2), and activation by VIP markedly enhanced the production of IL-22 and the barrier function of the epithelium. Conversely, deficiency in signaling through VIPR2 led to impaired production of IL-22 by ILC3s and increased susceptibility to inflammation-induced gut injury. Thus, intrinsic cellular rhythms acted in synergy with the cyclic patterns of food intake to drive the production of IL-22 and synchronize protection of the intestinal epithelium through a VIP-VIPR2 pathway in ILC3s.

A new lymphoid-primed progenitor marked by Dach1 downregulation identified with single cell multi-omics.

A classical view of blood cell development is that multipotent hematopoietic stem and progenitor cells (HSPCs) become lineage-restricted at defined stages. Lin-c-Kit+Sca-1+Flt3+ cells, termed lymphoid-primed multipotent progenitors (LMPPs), have lost megakaryocyte and erythroid potential but are heterogeneous in their fate. Here, through single-cell RNA sequencing, we identify the expression of Dach1 and associated genes in this fraction as being coexpressed with myeloid/stem genes but inversely correlated with lymphoid genes. Through generation of Dach1-GFP reporter mice, we identify a transcriptionally and functionally unique Dach1-GFP- subpopulation within LMPPs with lymphoid potential with low to negligible classic myeloid potential. We term these 'lymphoid-primed progenitors' (LPPs). These findings define an early definitive branch point of lymphoid development in hematopoiesis and a means for prospective isolation of LPPs.

Clonal multi-omics reveals Bcor as a negative regulator of emergency dendritic cell development.

Despite advances in single-cell multi-omics, a single stem or progenitor cell can only be tested once. We developed clonal multi-omics, in which daughters of a clone act as surrogates of the founder, thereby allowing multiple independent assays per clone. With SIS-seq, clonal siblings in parallel "sister" assays are examined either for gene expression by RNA sequencing (RNA-seq) or for fate in culture. We identified, and then validated using CRISPR, genes that controlled fate bias for different dendritic cell (DC) subtypes. This included Bcor as a suppressor of plasmacytoid DC (pDC) and conventional DC type 2 (cDC2) numbers during Flt3 ligand-mediated emergency DC development. We then developed SIS-skew to examine development of wild-type and Bcor-deficient siblings of the same clone in parallel. We found Bcor restricted clonal expansion, especially for cDC2s, and suppressed clonal fate potential, especially for pDCs. Therefore, SIS-seq and SIS-skew can reveal the molecular and cellular mechanisms governing clonal fate.

Apoptotic caspases suppress mtDNA-induced STING-mediated type I IFN production.

Activated caspases are a hallmark of apoptosis induced by the intrinsic pathway, but they are dispensable for cell death and the apoptotic clearance of cells in vivo. This has led to the suggestion that caspases are activated not just to kill but to prevent dying cells from triggering a host immune response. Here, we show that the caspase cascade suppresses type I interferon (IFN) production by cells undergoing Bak/Bax-mediated apoptosis. Bak and Bax trigger the release of mitochondrial DNA. This is recognized by the cGAS/STING-dependent DNA sensing pathway, which initiates IFN production. Activated caspases attenuate this response. Pharmacological caspase inhibition or genetic deletion of caspase-9, Apaf-1, or caspase-3/7 causes dying cells to secrete IFN-ß. In vivo, this precipitates an elevation in IFN-ß levels and consequent hematopoietic stem cell dysfunction, which is corrected by loss of Bak and Bax. Thus, the apoptotic caspase cascade functions to render mitochondrial apoptosis immunologically silent.

Systematic comparison of sequencing-based spatial transcriptomic methods.

Recent developments of sequencing-based spatial transcriptomics (sST) have catalyzed important advancements by facilitating transcriptome-scale spatial gene expression measurement. Despite this progress, efforts to comprehensively benchmark different platforms are currently lacking. The extant variability across technologies and datasets poses challenges in formulating standardized evaluation metrics. In this study, we established a collection of reference tissues and regions characterized by well-defined histological architectures, and used them to generate data to compare 11 sST methods. We highlighted molecular diffusion as a variable parameter across different methods and tissues, significantly affecting the effective resolutions. Furthermore, we observed that spatial transcriptomic data demonstrate unique attributes beyond merely adding a spatial axis to single-cell data, including an enhanced ability to capture patterned rare cell states along with specific markers, albeit being influenced by multiple factors including sequencing depth and resolution. Our study assists biologists in sST platform selection, and helps foster a consensus on evaluation standards and establish a framework for future benchmarking efforts that can be used as a gold standard for the development and benchmarking of computational tools for spatial transcriptomic analysis.

Systematic assessment of long-read RNA-seq methods for transcript identification and quantification.

The Long-read RNA-Seq Genome Annotation Assessment Project Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. Using different protocols and sequencing platforms, the consortium generated over 427 million long-read sequences from complementary DNA and direct RNA datasets, encompassing human, mouse and manatee species. Developers utilized these data to address challenges in transcript isoform detection, quantification and de novo transcript detection. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. Incorporating additional orthogonal data and replicate samples is advised when aiming to detect rare and novel transcripts or using reference-free approaches. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis.

Benchmarking long-read RNA-sequencing analysis tools using in silico mixtures.

The lack of benchmark data sets with inbuilt ground-truth makes it challenging to compare the performance of existing long-read isoform detection and differential expression analysis workflows. Here, we present a benchmark experiment using two human lung adenocarcinoma cell lines that were each profiled in triplicate together with synthetic, spliced, spike-in RNAs (sequins). Samples were deeply sequenced on both Illumina short-read and Oxford Nanopore Technologies long-read platforms. Alongside the ground-truth available via the sequins, we created in silico mixture samples to allow performance assessment in the absence of true positives or true negatives. Our results show that StringTie2 and bambu outperformed other tools from the six isoform detection tools tested, DESeq2, edgeR and limma-voom were best among the five differential transcript expression tools tested and there was no clear front-runner for performing differential transcript usage analysis between the five tools compared, which suggests further methods development is needed for this application.

Dividing out quantification uncertainty allows efficient assessment of differential transcript expression with edgeR.

Differential expression analysis of RNA-seq is one of the most commonly performed bioinformatics analyses. Transcript-level quantifications are inherently more uncertain than gene-level read counts because of ambiguous assignment of sequence reads to transcripts. While sequence reads can usually be assigned unambiguously to a gene, reads are very often compatible with multiple transcripts for that gene, particularly for genes with many isoforms. Software tools designed for gene-level differential expression do not perform optimally on transcript counts because the read-to-transcript ambiguity (RTA) disrupts the mean-variance relationship normally observed for gene level RNA-seq data and interferes with the efficiency of the empirical Bayes dispersion estimation procedures. The pseudoaligners kallisto and Salmon provide bootstrap samples from which quantification uncertainty can be assessed. We show that the overdispersion arising from RTA can be elegantly estimated by fitting a quasi-Poisson model to the bootstrap counts for each transcript. The technical overdispersion arising from RTA can then be divided out of the transcript counts, leading to scaled counts that can be input for analysis by established gene-level software tools with full statistical efficiency. Comprehensive simulations and test data show that an edgeR analysis of the scaled counts is more powerful and efficient than previous differential transcript expression pipelines while providing correct control of the false discovery rate. Simulations explore a wide range of scenarios including the effects of paired vs single-end reads, different read lengths and different numbers of replicates.

A method for stabilising the XX karyotype in female mESC cultures.

Female mouse embryonic stem cells (mESCs) present differently from male mESCs in several fundamental ways; however, complications with their in vitro culture have resulted in an under-representation of female mESCs in the literature. Recent studies show that the second X chromosome in female, and more specifically the transcriptional activity from both of these chromosomes due to absent X chromosome inactivation, sets female and male mESCs apart. To avoid this undesirable state, female mESCs in culture preferentially adopt an XO karyotype, with this adaption leading to loss of their unique properties in favour of a state that is near indistinguishable from male mESCs. If female pluripotency is to be studied effectively in this system, it is crucial that high-quality cultures of XX mESCs are available. Here, we report a method for better maintaining XX female mESCs in culture that also stabilises the male karyotype and makes study of female-specific pluripotency more feasible.

Flexiplex: a versatile demultiplexer and search tool for omics data.

MOTIVATION: The process of analyzing high throughput sequencing data often requires the identification and extraction of specific target sequences. This could include tasks, such as identifying cellular barcodes and UMIs in single-cell data, and specific genetic variants for genotyping. However, existing tools, which perform these functions are often task-specific, such as only demultiplexing barcodes for a dedicated type of experiment, or are not tolerant to noise in the sequencing data. RESULTS: To overcome these limitations, we developed Flexiplex, a versatile and fast sequence searching and demultiplexing tool for omics data, which is based on the Levenshtein distance and thus allows imperfect matches. We demonstrate Flexiplex's application on three use cases, identifying cell-line-specific sequences in Illumina short-read single-cell data, and discovering and demultiplexing cellular barcodes from noisy long-read single-cell RNA-seq data. We show that Flexiplex achieves an excellent balance of accuracy and computational efficiency compared to leading task-specific tools. AVAILABILITY AND IMPLEMENTATION: Flexiplex is available at https://davidsongroup.github.io/flexiplex/.

Hookworms dynamically respond to loss of Type 2 immune pressure.

The impact of the host immune environment on parasite transcription and fitness is currently unknown. It is widely held that hookworm infections have an immunomodulatory impact on the host, but whether the converse is true remains unclear. Immunity against adult-stage hookworms is largely mediated by Type 2 immune responses driven by the transcription factor Signal Transducer and Activator of Transcription 6 (STAT6). This study investigated whether serial passage of the rodent hookworm Nippostrongylus brasiliensis in STAT6-deficient mice (STAT6 KO) caused changes in parasites over time. After adaptation to STAT6 KO hosts, N. brasiliensis increased their reproductive output, feeding capacity, energy content, and body size. Using an improved N. brasiliensis genome, we found that these physiological changes corresponded with a dramatic shift in the transcriptional landscape, including increased expression of gene pathways associated with egg production, but a decrease in genes encoding neuropeptides, proteases, SCP/TAPS proteins, and transthyretin-like proteins; the latter three categories have been repeatedly observed in hookworm excreted/secreted proteins (ESPs) implicated in immunosuppression. Although transcriptional changes started to appear in the first generation of passage in STAT6 KO hosts for both immature and mature adult stages, downregulation of the genes putatively involved in immunosuppression was only observed after multiple generations in this immunodeficient environment. When STAT6 KO-adapted N. brasiliensis were reintroduced to a naive WT host after up to 26 generations, this progressive change in host-adaptation corresponded to increased production of inflammatory cytokines by the WT host. Surprisingly, however, this single exposure of STAT6 KO-adapted N. brasiliensis to WT hosts resulted in worms that were morphologically and transcriptionally indistinguishable from WT-adapted parasites. This work uncovers remarkable plasticity in the ability of hookworms to adapt to their hosts, which may present a general feature of parasitic nematodes.

Actinobacillus pleuropneumoniae encodes multiple phase-variable DNA methyltransferases that control distinct phasevarions.

Actinobacillus pleuropneumoniae is the cause of porcine pleuropneumonia, a severe respiratory tract infection that is responsible for major economic losses to the swine industry. Many host-adapted bacterial pathogens encode systems known as phasevarions (phase-variable regulons). Phasevarions result from variable expression of cytoplasmic DNA methyltransferases. Variable expression results in genome-wide methylation differences within a bacterial population, leading to altered expression of multiple genes via epigenetic mechanisms. Our examination of a diverse population of A. pleuropneumoniae strains determined that Type I and Type III DNA methyltransferases with the hallmarks of phase variation were present in this species. We demonstrate that phase variation is occurring in these methyltransferases, and show associations between particular Type III methyltransferase alleles and serovar. Using Pacific BioSciences Single-Molecule, Real-Time (SMRT) sequencing and Oxford Nanopore sequencing, we demonstrate the presence of the first ever characterised phase-variable, cytosine-specific Type III DNA methyltransferase. Phase variation of distinct Type III DNA methyltransferase in A. pleuropneumoniae results in the regulation of distinct phasevarions, and in multiple phenotypic differences relevant to pathobiology. Our characterisation of these newly described phasevarions in A. pleuropneumoniae will aid in the selection of stably expressed antigens, and direct and inform development of a rationally designed subunit vaccine against this major veterinary pathogen.

Arginine-rich C9ORF72 ALS proteins stall ribosomes in a manner distinct from a canonical ribosome-associated quality control substrate.

Hexanucleotide expansion mutations in C9ORF72 are a frequent cause of amyotrophic lateral sclerosis. We previously reported that long arginine-rich dipeptide repeats (DPRs), mimicking abnormal proteins expressed from the hexanucleotide expansion, caused translation stalling when expressed in cell culture models. Whether this stalling provides a mechanism of pathogenicity remains to be determined. Here, we explored the molecular features of DPR-induced stalling and examined whether known mechanisms such as ribosome quality control (RQC) regulate translation elongation on sequences that encode arginine-rich DPRs. We demonstrate that arginine-rich DPRs lead to stalling in a length-dependent manner, with lengths longer than 40 repeats invoking severe translation arrest. Mutational screening of 40×Gly-Xxx DPRs shows that stalling is most pronounced when Xxx is a charged amino acid (Arg, Lys, Glu, or Asp). Through a genome-wide knockout screen, we find that genes regulating stalling on polyadenosine mRNA coding for poly-Lys, a canonical RQC substrate, act differently in the case of arginine-rich DPRs. Indeed, these findings point to a limited scope for natural regulatory responses to resolve the arginine-rich DPR stalls, even though the stalls may be sensed, as evidenced by an upregulation of RQC gene expression. These findings therefore implicate arginine-rich DPR-mediated stalled ribosomes as a source of stress and toxicity and may be a crucial component in pathomechanisms.

A functional genetic screen identifies aurora kinase b as an essential regulator of Sox9-positive mouse embryonic lung progenitor cells.

Development of a branching tree in the embryonic lung is crucial for the formation of a fully mature functional lung at birth. Sox9+ cells present at the tip of the primary embryonic lung endoderm are multipotent cells responsible for branch formation and elongation. We performed a genetic screen in murine primary cells and identified aurora kinase b (Aurkb) as an essential regulator of Sox9+ cells ex vivo. In vivo conditional knockout studies confirmed that Aurkb was required for lung development but was not necessary for postnatal growth and the repair of the adult lung after injury. Deletion of Aurkb in embryonic Sox9+ cells led to the formation of a stunted lung that retained the expression of Sox2 in the proximal airways, as well as Sox9 in the distal tips. Although we found no change in cell polarity, we showed that loss of Aurkb or chemical inhibition of Aurkb caused Sox9+ cells to arrest at G2/M, likely responsible for the lack of branch bifurcation. This work demonstrates the power of genetic screens in identifying novel regulators of Sox9+ progenitor cells and lung branching morphogenesis.

Single-cell multiomics reveal the scale of multilayered adaptations enabling CLL relapse during venetoclax therapy.

Venetoclax (VEN) inhibits the prosurvival protein BCL2 to induce apoptosis and is a standard therapy for chronic lymphocytic leukemia (CLL), delivering high complete remission rates and prolonged progression-free survival in relapsed CLL but with eventual loss of efficacy. A spectrum of subclonal genetic changes associated with VEN resistance has now been described. To fully understand clinical resistance to VEN, we combined single-cell short- and long-read RNA-sequencing to reveal the previously unappreciated scale of genetic and epigenetic changes underpinning acquired VEN resistance. These appear to be multilayered. One layer comprises changes in the BCL2 family of apoptosis regulators, especially the prosurvival family members. This includes previously described mutations in BCL2 and amplification of the MCL1 gene but is heterogeneous across and within individual patient leukemias. Changes in the proapoptotic genes are notably uncommon, except for single cases with subclonal losses of BAX or NOXA. Much more prominent was universal MCL1 gene upregulation. This was driven by an overlying layer of emergent NF-κB (nuclear factor kappa B) activation, which persisted in circulating cells during VEN therapy. We discovered that MCL1 could be a direct transcriptional target of NF-κB. Both the switch to alternative prosurvival factors and NF-κB activation largely dissipate following VEN discontinuation. Our studies reveal the extent of plasticity of CLL cells in their ability to evade VEN-induced apoptosis. Importantly, these findings pinpoint new approaches to circumvent VEN resistance and provide a specific biological justification for the strategy of VEN discontinuation once a maximal response is achieved rather than maintaining long-term selective pressure with the drug.

FGF-independent MEK1/2 signalling in the developing foetal testis is essential for male germline differentiation in mice.

BACKGROUND: Disrupted germline differentiation or compromised testis development can lead to subfertility or infertility and are strongly associated with testis cancer in humans. In mice, SRY and SOX9 induce expression of Fgf9, which promotes Sertoli cell differentiation and testis development. FGF9 is also thought to promote male germline differentiation but the mechanism is unknown. FGFs typically signal through mitogen-activated protein kinases (MAPKs) to phosphorylate ERK1/2 (pERK1/2). We explored whether FGF9 regulates male germline development through MAPK by inhibiting either FGF or MEK1/2 signalling in the foetal testis immediately after gonadal sex determination and testis cord formation, but prior to male germline commitment. RESULTS: pERK1/2 was detected in Sertoli cells and inhibition of MEK1/2 reduced Sertoli cell proliferation and organisation and resulted in some germ cells localised outside of the testis cords. While pERK1/2 was not detected in germ cells, inhibition of MEK1/2 after somatic sex determination profoundly disrupted germ cell mitotic arrest, dysregulated a broad range of male germline development genes and prevented the upregulation of key male germline markers, DPPA4 and DNMT3L. In contrast, while FGF inhibition reduced Sertoli cell proliferation, expression of male germline markers was unaffected and germ cells entered mitotic arrest normally. While male germline differentiation was not disrupted by FGF inhibition, a range of stem cell and cancer-associated genes were commonly altered after 24 h of FGF or MEK1/2 inhibition, including genes involved in the maintenance of germline stem cells, Nodal signalling, proliferation, and germline cancer. CONCLUSIONS: Together, these data demonstrate a novel role for MEK1/2 signalling during testis development that is essential for male germline differentiation, but indicate a more limited role for FGF signalling. Our data indicate that additional ligands are likely to act through MEK1/2 to promote male germline differentiation and highlight a need for further mechanistic understanding of male germline development.

Benchmarking single cell RNA-sequencing analysis pipelines using mixture control experiments.

Single cell RNA-sequencing (scRNA-seq) technology has undergone rapid development in recent years, leading to an explosion in the number of tailored data analysis methods. However, the current lack of gold-standard benchmark datasets makes it difficult for researchers to systematically compare the performance of the many methods available. Here, we generated a realistic benchmark experiment that included single cells and admixtures of cells or RNA to create 'pseudo cells' from up to five distinct cancer cell lines. In total, 14 datasets were generated using both droplet and plate-based scRNA-seq protocols. We compared 3,913 combinations of data analysis methods for tasks ranging from normalization and imputation to clustering, trajectory analysis and data integration. Evaluation revealed pipelines suited to different types of data for different tasks. Our data and analysis provide a comprehensive framework for benchmarking most common scRNA-seq analysis steps.

The EMT modulator SNAI1 contributes to AML pathogenesis via its interaction with LSD1.

Modulators of epithelial-to-mesenchymal transition (EMT) have recently emerged as novel players in the field of leukemia biology. The mechanisms by which EMT modulators contribute to leukemia pathogenesis, however, remain to be elucidated. Here we show that overexpression of SNAI1, a key modulator of EMT, is a pathologically relevant event in human acute myeloid leukemia (AML) that contributes to impaired differentiation, enhanced self-renewal, and proliferation of immature myeloid cells. We demonstrate that ectopic expression of Snai1 in hematopoietic cells predisposes mice to AML development. This effect is mediated by interaction with the histone demethylase KDM1A/LSD1. Our data shed new light on the role of SNAI1 in leukemia development and identify a novel mechanism of LSD1 corruption in cancer. This is particularly pertinent given the current interest surrounding the use of LSD1 inhibitors in the treatment of multiple different malignancies, including AML.

The impact of influenza pulmonary infection and inflammation on vagal bronchopulmonary sensory neurons.

Influenza A virus (IAV) is rapidly detected in the airways by the immune system, with resident parenchymal cells and leukocytes orchestrating viral sensing and the induction of antiviral inflammatory responses. The airways are innervated by heterogeneous populations of vagal sensory neurons which also play an important role in pulmonary defense. How these neurons respond to IAV respiratory infection remains unclear. Here, we use a murine model to provide the first evidence that vagal sensory neurons undergo significant transcriptional changes following a respiratory IAV infection. RNA sequencing on vagal sensory ganglia showed that IAV infection induced the expression of many genes associated with an antiviral and pro-inflammatory response and this was accompanied by a significant increase in inflammatory cell recruitment into the vagal ganglia. Assessment of gene expression in single-vagal sensory neurons confirmed that IAV infection induced a neuronal inflammatory phenotype, which was most prominent in bronchopulmonary neurons, and also evident in some neurons innervating other organs. The altered transcriptome could be mimicked by intranasal treatment with cytokines and the lung homogenates of infected mice, in the absence of infectious virus. These data argue that IAV pulmonary infection and subsequent inflammation induces vagal sensory ganglia neuroinflammation and this may have important implications for IAV-induced morbidity.