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.

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.

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.

An optimized protocol for quality control of gene therapy vectors using nanopore direct RNA sequencing.

Despite recent advances made toward improving the efficacy of lentiviral gene therapies, a sizeable proportion of produced vector contains an incomplete and thus potentially nonfunctional RNA genome. This can undermine gene delivery by the lentivirus as well as increase manufacturing costs and must be improved to facilitate the widespread clinical implementation of lentiviral gene therapies. Here, we compare three long-read sequencing technologies for their ability to detect issues in vector design and determine nanopore direct RNA sequencing to be the most powerful. We show how this approach identifies and quantifies incomplete RNA caused by cryptic splicing and polyadenylation sites, including a potential cryptic polyadenylation site in the widely used Woodchuck Hepatitis Virus Posttranscriptional Regulatory Element (WPRE). Using artificial polyadenylation of the lentiviral RNA, we also identify multiple hairpin-associated truncations in the analyzed lentiviral vectors (LVs), which account for most of the detected RNA fragments. Finally, we show that these insights can be used for the optimization of LV design. In summary, nanopore direct RNA sequencing is a powerful tool for the quality control and optimization of LVs, which may help to improve lentivirus manufacturing and thus the development of higher quality lentiviral gene therapies.

Measuring X-Chromosome inactivation skew for X-linked diseases with adaptive nanopore sequencing.

X-linked genetic disorders typically affect females less severely than males owing to the presence of a second X Chromosome not carrying the deleterious variant. However, the phenotypic expression in females is highly variable, which may be explained by an allelic skew in X-Chromosome inactivation. Accurate measurement of X inactivation skew is crucial to understand and predict disease phenotype in carrier females, with prediction especially relevant for degenerative conditions. We propose a novel approach using nanopore sequencing to quantify skewed X inactivation accurately. By phasing sequence variants and methylation patterns, this single assay reveals the disease variant, X inactivation skew, and its directionality and is applicable to all patients and X-linked variants. Enrichment of X Chromosome reads through adaptive sampling enhances cost-efficiency. Our study includes a cohort of 16 X-linked variant carrier females affected by two X-linked inherited retinal diseases: choroideremia and RPGR-associated retinitis pigmentosa. As retinal DNA cannot be readily obtained, we instead determine the skew from peripheral samples (blood, saliva, and buccal mucosa) and correlate it to phenotypic outcomes. This reveals a strong correlation between X inactivation skew and disease presentation, confirming the value in performing this assay and its potential as a way to prioritize patients for early intervention, such as gene therapy currently in clinical trials for these conditions. Our method of assessing skewed X inactivation is applicable to all long-read genomic data sets, providing insights into disease risk and severity and aiding in the development of individualized strategies for X-linked variant carrier females.

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.

Venetoclax Dose Escalation Rapidly Activates a BAFF/BCL-2 Survival Axis in Chronic Lymphocytic Leukemia.

Venetoclax, a first-in-class BH3 mimetic drug targeting BCL-2, has improved outcomes for patients with chronic lymphocytic leukemia (CLL). Early measurements of the depth of the venetoclax treatment response, assessed by minimal residual disease, are strong predictors of long-term clinical outcomes. Yet, there are limited data concerning the early changes induced by venetoclax treatment that might inform strategies to improve responses. To address this gap, we conducted longitudinal mass cytometric profiling of blood cells from patients with CLL during the first five weeks of venetoclax monotherapy. At baseline, we resolved CLL heterogeneity at the single-cell level to define multiple subpopulations in all patients distinguished by proliferative, metabolic and cell survival proteins. Venetoclax induced significant reduction in all CLL subpopulations coincident with rapid upregulation of pro-survival BCL-2, BCL-XL and MCL-1 proteins in surviving cells, which had reduced sensitivity to the drug. Mouse models recapitulated the venetoclax-induced elevation of survival proteins in B cells and CLL-like cells that persisted in vivo, with genetic models demonstrating that extensive apoptosis and access to the B cell cytokine, BAFF, were essential. Accordingly, analysis of patients with CLL that were treated with venetoclax or the anti-CD20 antibody obinutuzumab exhibited marked elevation of BAFF and increased pro-survival proteins in leukemic cells that persisted. Overall, these data highlight the rapid adaptation of CLL cells to targeted therapies via homeostatic factors and support co-targeting of cytokine signals to achieve deeper and more durable long-term responses.

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.

VEGF-dependent testicular vascularisation involves MEK1/2 signalling and the essential angiogenesis factors, SOX7 and SOX17.

BACKGROUND: Abnormalities of in utero testis development are strongly associated with reproductive health conditions, including male infertility and testis cancer. In mouse testes, SOX9 and FGF9 support Sertoli cell development, while VEGF signalling is essential for the establishment of vasculature. The mitogen-activated protein kinase (MAPK) pathway is a major signalling cascade, essential for cell proliferation, differentiation and activation of Sry during primary sex-determination, but little is known about its function during fetal testis morphogenesis. We explored potential functions of MAPK signalling immediately after the establishment of testis cords in embryonic day (E)12.5 Oct4-eGFP transgenic mouse testes cultured using a MEK1/2 inhibitor. RESULTS: RNA sequencing in isolated gonadal somatic cells identified 116 and 114 differentially expressed genes after 24 and 72 h of MEK1/2 inhibition, respectively. Ingenuity Pathway Analysis revealed an association of MEK1/2 signalling with biological functions such as angiogenesis, vasculogenesis and cell migration. This included a failure to upregulate the master transcriptional regulators of vascular development, Sox7 and Sox17, VEGF receptor genes, the cell adhesion factor gene Cd31 and a range of other endothelial cell markers such as Cdh5 (encoding VE-cadherin) and gap junction genes Gja4 and Gja5. In contrast, only a small number of Sertoli cell enriched genes were affected. Immunofluorescent analyses of control testes revealed that the MEK1/2 downstream target, ERK1/2 was phosphorylated in endothelial cells and Sertoli cells. Inhibition of MEK1/2 eliminated pERK1/2 in fetal testes, and CD31, VE-cadherin, SOX7 and SOX17 and endothelial cells were lost. Consistent with a role for VEGF in driving endothelial cell development in the testis, inhibition of VEGFR also abrogated pERK1/2 and SOX7 and SOX17 expressing endothelial cells. Moreover, while Sertoli cell proliferation and localisation to the testis cord basement membrane was disrupted by inhibition of MEK1/2, it was unaffected by VEGFR inhibition. Instead, inhibition of FGF signalling compromised Sertoli cell proliferation and localisation to the testis cord basement membrane. CONCLUSIONS: Together, our data highlight an essential role for VEGF-dependent MEK1/2 signalling in promoting vasculature and indicate that FGF signalling through MEK1/2 regulates Sertoli cell organisation in the developing mouse testis.

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.