Developmental Origin Governs CD8+ T Cell Fate Decisions during Infection.

Heterogeneity is a hallmark feature of the adaptive immune system in vertebrates. Following infection, naive T cells differentiate into various subsets of effector and memory T cells, which help to eliminate pathogens and maintain long-term immunity. The current model suggests there is a single lineage of naive T cells that give rise to different populations of effector and memory T cells depending on the type and amounts of stimulation they encounter during infection. Here, we have discovered that multiple sub-populations of cells exist in the naive CD8+ T cell pool that are distinguished by their developmental origin, unique transcriptional profiles, distinct chromatin landscapes, and different kinetics and phenotypes after microbial challenge. These data demonstrate that the naive CD8+ T cell pool is not as homogeneous as previously thought and offers a new framework for explaining the remarkable heterogeneity in the effector and memory T cell subsets that arise after infection.

CDK2 kinase activity is a regulator of male germ cell fate.

The ability of men to remain fertile throughout their lives depends upon establishment of a spermatogonial stem cell (SSC) pool from gonocyte progenitors, and thereafter balancing SSC renewal versus terminal differentiation. Here, we report that precise regulation of the cell cycle is crucial for this balance. Whereas cyclin-dependent kinase 2 (Cdk2) is not necessary for mouse viability or gametogenesis stages prior to meiotic prophase I, mice bearing a deregulated allele (Cdk2Y15S ) are severely deficient in spermatogonial differentiation. This allele disrupts an inhibitory phosphorylation site (Tyr15) for the kinase WEE1. Remarkably, Cdk2Y15S/Y15S mice possess abnormal clusters of mitotically active SSC-like cells, but these are eventually removed by apoptosis after failing to differentiate properly. Analyses of lineage markers, germ cell proliferation over time, and single cell RNA-seq data revealed delayed and defective differentiation of gonocytes into SSCs. Biochemical and genetic data demonstrated that Cdk2Y15S is a gain-of-function allele causing elevated kinase activity, which underlies these differentiation defects. Our results demonstrate that precise regulation of CDK2 kinase activity in male germ cell development is crucial for the gonocyte-to-spermatogonia transition and long-term spermatogenic homeostasis.

Robust partitioning of microRNA targets from downstream regulatory changes.

The biological impact of microRNAs (miRNAs) is determined by their targets, and robustly identifying direct miRNA targets remains challenging. Existing methods suffer from high false-positive rates and are unable to effectively differentiate direct miRNA targets from downstream regulatory changes. Here, we present an experimental and computational framework to deconvolute post-transcriptional and transcriptional changes using a combination of RNA-seq and PRO-seq. This novel approach allows us to systematically profile the regulatory impact of a miRNA. We refer to this approach as CARP: Combined Analysis of RNA-seq and PRO-seq. We apply CARP to multiple miRNAs and show that it robustly distinguishes direct targets from downstream changes, while greatly reducing false positives. We validate our approach using Argonaute eCLIP-seq and ribosome profiling, demonstrating that CARP defines a comprehensive repertoire of targets. Using this approach, we identify miRNA-specific activity of target sites within the open reading frame. Additionally, we show that CARP facilitates the dissection of complex changes in gene regulatory networks triggered by miRNAs and identification of transcription factors that mediate downstream regulatory changes. Given the robustness of the approach, CARP would be particularly suitable for dissecting miRNA regulatory networks in vivo.

A draft genome sequence of the pulse crop chickpea (Cicer arietinum L.).

Cicer arietinum L. (chickpea) is the third most important food legume crop. We have generated the draft sequence of a desi-type chickpea genome using next-generation sequencing platforms, bacterial artificial chromosome end sequences and a genetic map. The 520-Mb assembly covers 70% of the predicted 740-Mb genome length, and more than 80% of the gene space. Genome analysis predicts the presence of 27,571 genes and 210 Mb as repeat elements. The gene expression analysis performed using 274 million RNA-Seq reads identified several tissue-specific and stress-responsive genes. Although segmental duplicated blocks are observed, the chickpea genome does not exhibit any indication of recent whole-genome duplication. Nucleotide diversity analysis provides an assessment of a narrow genetic base within the chickpea cultivars. We have developed a resource for genetic markers by comparing the genome sequences of one wild and three cultivated chickpea genotypes. The draft genome sequence is expected to facilitate genetic enhancement and breeding to develop improved chickpea varieties.

Distinct pulmonary and systemic effects of dexamethasone in severe COVID-19.

Dexamethasone is the standard of care for critically ill patients with COVID-19, but the mechanisms by which it decreases mortality and its immunological effects in this setting are not understood. Here we perform bulk and single-cell RNA sequencing of samples from the lower respiratory tract and blood, and assess plasma cytokine profiling to study the effects of dexamethasone on both systemic and pulmonary immune cell compartments. In blood samples, dexamethasone is associated with decreased expression of genes associated with T cell activation, including TNFSFR4 and IL21R. We also identify decreased expression of several immune pathways, including major histocompatibility complex-II signaling, selectin P ligand signaling, and T cell recruitment by intercellular adhesion molecule and integrin activation, suggesting these are potential mechanisms of the therapeutic benefit of steroids in COVID-19. We identify additional compartment- and cell- specific differences in the effect of dexamethasone that are reproducible in publicly available datasets, including steroid-resistant interferon pathway expression in the respiratory tract, which may be additional therapeutic targets. In summary, we demonstrate compartment-specific effects of dexamethasone in critically ill COVID-19 patients, providing mechanistic insights with potential therapeutic relevance. Our results highlight the importance of studying compartmentalized inflammation in critically ill patients.

The gene regulatory basis of bystander activation in CD8+ T cells.

CD8+ T cells are classically recognized as adaptive lymphocytes based on their ability to recognize specific foreign antigens and mount memory responses. However, recent studies indicate that some antigen-inexperienced CD8+ T cells can respond to innate cytokines alone in the absence of cognate T cell receptor stimulation, a phenomenon referred to as bystander activation. Here, we demonstrate that neonatal CD8+ T cells undergo a robust and diverse program of bystander activation, which corresponds to enhanced innate-like protection against unrelated pathogens. Using a multi-omics approach, we found that the ability of neonatal CD8+ T cells to respond to innate cytokines derives from their capacity to undergo rapid chromatin remodeling, resulting in the usage of a distinct set of enhancers and transcription factors typically found in innate-like T cells. We observed that the switch between innate and adaptive functions in the CD8+ T cell compartment is mediated by changes in the abundance of distinct subsets of cells. The innate CD8+ T cell subset that predominates in early life was also present in adult mice and humans. Our findings provide support for the layered immune hypothesis and indicate that the CD8+ T cell compartment is more functionally diverse than previously thought.

Single-cell and spatial multi-omics highlight effects of anti-integrin therapy across cellular compartments in ulcerative colitis.

Ulcerative colitis (UC) is driven by immune and stromal subsets, culminating in epithelial injury. Vedolizumab (VDZ) is an anti-integrin antibody that is effective for treating UC. VDZ is known to inhibit lymphocyte trafficking to the intestine, but its broader effects on other cell subsets are less defined. To identify the inflammatory cells that contribute to colitis and are affected by VDZ, we performed single-cell transcriptomic and proteomic analyses of peripheral blood and colonic biopsies in healthy controls and patients with UC on VDZ or other therapies. Here we show that VDZ treatment is associated with alterations in circulating and tissue mononuclear phagocyte (MNP) subsets, along with modest shifts in lymphocytes. Spatial multi-omics of formalin-fixed biopsies demonstrates trends towards increased abundance and proximity of MNP and fibroblast subsets in active colitis. Spatial transcriptomics of archived specimens pre-treatment identifies epithelial-, MNP-, and fibroblast-enriched genes related to VDZ responsiveness, highlighting important roles for these subsets in UC.

Single-cell and spatial multi-omics highlight effects of anti-integrin therapy across cellular compartments in ulcerative colitis.

Ulcerative colitis (UC) is driven by immune and stromal subsets, culminating in epithelial injury. Vedolizumab (VDZ) is an anti-integrin antibody that is effective for treating UC. VDZ is known to inhibit lymphocyte trafficking to the intestine, but its broader effects on other cell subsets are less defined. To identify the inflammatory cells that contribute to colitis and are affected by VDZ, we perform single-cell transcriptomic and proteomic analyses of peripheral blood and colonic biopsies in healthy controls and patients with UC on VDZ or other therapies. Here we show that VDZ treatment is associated with alterations in circulating and tissue mononuclear phagocyte (MNP) subsets, along with modest shifts in lymphocytes. Spatial multi-omics of formalin-fixed biopsies demonstrates trends towards increased abundance and proximity of MNP and fibroblast subsets in active colitis. Spatial transcriptomics of archived specimens pre-treatment identifies epithelial-, MNP-, and fibroblast-enriched genes related to VDZ responsiveness, highlighting important roles for these subsets in UC.

Integrated longitudinal multiomics study identifies immune programs associated with acute COVID-19 severity and mortality.

BACKGROUNDPatients hospitalized for COVID-19 exhibit diverse clinical outcomes, with outcomes for some individuals diverging over time even though their initial disease severity appears similar to that of other patients. A systematic evaluation of molecular and cellular profiles over the full disease course can link immune programs and their coordination with progression heterogeneity.METHODSWe performed deep immunophenotyping and conducted longitudinal multiomics modeling, integrating 10 assays for 1,152 Immunophenotyping Assessment in a COVID-19 Cohort (IMPACC) study participants and identifying several immune cascades that were significant drivers of differential clinical outcomes.RESULTSIncreasing disease severity was driven by a temporal pattern that began with the early upregulation of immunosuppressive metabolites and then elevated levels of inflammatory cytokines, signatures of coagulation, formation of neutrophil extracellular traps, and T cell functional dysregulation. A second immune cascade, predictive of 28-day mortality among critically ill patients, was characterized by reduced total plasma Igs and B cells and dysregulated IFN responsiveness. We demonstrated that the balance disruption between IFN-stimulated genes and IFN inhibitors is a crucial biomarker of COVID-19 mortality, potentially contributing to failure of viral clearance in patients with fatal illness.CONCLUSIONOur longitudinal multiomics profiling study revealed temporal coordination across diverse omics that potentially explain the disease progression, providing insights that can inform the targeted development of therapies for patients hospitalized with COVID-19, especially those who are critically ill.TRIAL REGISTRATIONClinicalTrials.gov NCT04378777.FUNDINGNIH (5R01AI135803-03, 5U19AI118608-04, 5U19AI128910-04, 4U19AI090023-11, 4U19AI118610-06, R01AI145835-01A1S1, 5U19AI062629-17, 5U19AI057229-17, 5U19AI125357-05, 5U19AI128913-03, 3U19AI077439-13, 5U54AI142766-03, 5R01AI104870-07, 3U19AI089992-09, 3U19AI128913-03, and 5T32DA018926-18); NIAID, NIH (3U19AI1289130, U19AI128913-04S1, and R01AI122220); and National Science Foundation (DMS2310836).

Astroglial toxicity promotes synaptic degeneration in the thalamocortical circuit in frontotemporal dementia with GRN mutations.

Mutations in the human progranulin (GRN) gene are a leading cause of frontotemporal lobar degeneration (FTLD). While previous studies implicate aberrant microglial activation as a disease-driving factor in neurodegeneration in the thalamocortical circuit in Grn-/- mice, the exact mechanism for neurodegeneration in FTLD-GRN remains unclear. By performing comparative single-cell transcriptomics in the thalamus and frontal cortex of Grn-/- mice and patients with FTLD-GRN, we have uncovered a highly conserved astroglial pathology characterized by upregulation of gap junction protein GJA1, water channel AQP4, and lipid-binding protein APOE, and downregulation of glutamate transporter SLC1A2 that promoted profound synaptic degeneration across the two species. This astroglial toxicity could be recapitulated in mouse astrocyte-neuron cocultures and by transplanting induced pluripotent stem cell-derived astrocytes to cortical organoids, where progranulin-deficient astrocytes promoted synaptic degeneration, neuronal stress, and TDP-43 proteinopathy. Together, these results reveal a previously unappreciated astroglial pathology as a potential key mechanism in neurodegeneration in FTLD-GRN.

Distinct pulmonary and systemic effects of dexamethasone in severe COVID-19.

Dexamethasone is the standard of care for critically ill patients with COVID-19, but the mechanisms by which it decreases mortality and its immunological effects in this setting are not understood. We performed bulk and single-cell RNA sequencing of the lower respiratory tract and blood, and plasma cytokine profiling to study the effect of dexamethasone on systemic and pulmonary immune cells. We find decreased signatures of antigen presentation, T cell recruitment, and viral injury in patients treated with dexamethasone. We identify compartment- and cell- specific differences in the effect of dexamethasone in patients with severe COVID-19 that are reproducible in publicly available datasets. Our results highlight the importance of studying compartmentalized inflammation in critically ill patients.

Cyclone: an accessible pipeline to analyze, evaluate and optimize multiparametric cytometry data.

In the past decade, high-dimensional single cell technologies have revolutionized basic and translational immunology research and are now a key element of the toolbox used by scientists to study the immune system. However, analysis of the data generated by these approaches often requires clustering algorithms and dimensionality reduction representation which are computationally intense and difficult to evaluate and optimize. Here we present Cyclone, an analysis pipeline integrating dimensionality reduction, clustering, evaluation and optimization of clustering resolution, and downstream visualization tools facilitating the analysis of a wide range of cytometry data. We benchmarked and validated Cyclone on mass cytometry (CyTOF), full spectrum fluorescence-based cytometry, and multiplexed immunofluorescence (IF) in a variety of biological contexts, including infectious diseases and cancer. In each instance, Cyclone not only recapitulates gold standard immune cell identification, but also enables the unsupervised identification of lymphocytes and mononuclear phagocytes subsets that are associated with distinct biological features. Altogether, the Cyclone pipeline is a versatile and accessible pipeline for performing, optimizing, and evaluating clustering on variety of cytometry datasets which will further power immunology research and provide a scaffold for biological discovery.

Cyclone: an accessible pipeline to analyze, evaluate, and optimize multiparametric cytometry data.

In the past decade, high-dimensional single-cell technologies have revolutionized basic and translational immunology research and are now a key element of the toolbox used by scientists to study the immune system. However, analysis of the data generated by these approaches often requires clustering algorithms and dimensionality reduction representation, which are computationally intense and difficult to evaluate and optimize. Here, we present Cytometry Clustering Optimization and Evaluation (Cyclone), an analysis pipeline integrating dimensionality reduction, clustering, evaluation, and optimization of clustering resolution, and downstream visualization tools facilitating the analysis of a wide range of cytometry data. We benchmarked and validated Cyclone on mass cytometry (CyTOF), full-spectrum fluorescence-based cytometry, and multiplexed immunofluorescence (IF) in a variety of biological contexts, including infectious diseases and cancer. In each instance, Cyclone not only recapitulates gold standard immune cell identification but also enables the unsupervised identification of lymphocytes and mononuclear phagocyte subsets that are associated with distinct biological features. Altogether, the Cyclone pipeline is a versatile and accessible pipeline for performing, optimizing, and evaluating clustering on a variety of cytometry datasets, which will further power immunology research and provide a scaffold for biological discovery.

Integrated longitudinal multi-omics study identifies immune programs associated with COVID-19 severity and mortality in 1152 hospitalized participants.

Hospitalized COVID-19 patients exhibit diverse clinical outcomes, with some individuals diverging over time even though their initial disease severity appears similar. A systematic evaluation of molecular and cellular profiles over the full disease course can link immune programs and their coordination with progression heterogeneity. In this study, we carried out deep immunophenotyping and conducted longitudinal multi-omics modeling integrating ten distinct assays on a total of 1,152 IMPACC participants and identified several immune cascades that were significant drivers of differential clinical outcomes. Increasing disease severity was driven by a temporal pattern that began with the early upregulation of immunosuppressive metabolites and then elevated levels of inflammatory cytokines, signatures of coagulation, NETosis, and T-cell functional dysregulation. A second immune cascade, predictive of 28-day mortality among critically ill patients, was characterized by reduced total plasma immunoglobulins and B cells, as well as dysregulated IFN responsiveness. We demonstrated that the balance disruption between IFN-stimulated genes and IFN inhibitors is a crucial biomarker of COVID-19 mortality, potentially contributing to the failure of viral clearance in patients with fatal illness. Our longitudinal multi-omics profiling study revealed novel temporal coordination across diverse omics that potentially explain disease progression, providing insights that inform the targeted development of therapies for hospitalized COVID-19 patients, especially those critically ill.

Gene discovery and tissue-specific transcriptome analysis in chickpea with massively parallel pyrosequencing and web resource development.

Chickpea (Cicer arietinum) is an important food legume crop but lags in the availability of genomic resources. In this study, we have generated about 2 million high-quality sequences of average length of 372 bp using pyrosequencing technology. The optimization of de novo assembly clearly indicated that hybrid assembly of long-read and short-read primary assemblies gave better results. The hybrid assembly generated a set of 34,760 transcripts with an average length of 1,020 bp representing about 4.8% (35.5 Mb) of the total chickpea genome. We identified more than 4,000 simple sequence repeats, which can be developed as functional molecular markers in chickpea. Putative function and Gene Ontology terms were assigned to at least 73.2% and 71.0% of chickpea transcripts, respectively. We have also identified several chickpea transcripts that showed tissue-specific expression and validated the results using real-time polymerase chain reaction analysis. Based on sequence comparison with other species within the plant kingdom, we identified two sets of lineage-specific genes, including those conserved in the Fabaceae family (legume specific) and those lacking significant similarity with any non chickpea species (chickpea specific). Finally, we have developed a Web resource, Chickpea Transcriptome Database, which provides public access to the data and results reported in this study. The strategy for optimization of de novo assembly presented here may further facilitate the transcriptome sequencing and characterization in other organisms. Most importantly, the data and results reported in this study will help to accelerate research in various areas of genomics and implementing breeding programs in chickpea.

Inhibition of MET Signaling with Ficlatuzumab in Combination with Chemotherapy in Refractory AML: Clinical Outcomes and High-Dimensional Analysis.

Acute myeloid leukemia patients refractory to induction therapy or relapsed within one year have poor outcomes. Autocrine production of hepatocyte growth factor by myeloid blasts drives leukemogenesis in pre-clinical models. A phase Ib trial evaluated ficlatuzumab, a first-in-class anti-HGF antibody, in combination with cytarabine in this high-risk population. Dose-limiting toxicities were not observed, and 20 mg/kg was established as the recommended phase II dose. The most frequent treatment-related adverse event was febrile neutropenia. Among 17 evaluable patients, the overall response rate was 53%, all complete remissions. Phospho-proteomic mass cytometry showed potent on-target suppression of p-MET after ficlatuzumab treatment and that attenuation of p-S6 was associated with clinical response. Multiplexed single cell RNA sequencing using prospectively acquired patient specimens identified interferon response genes as adverse predictive factors. The ficlatuzumab and cytarabine combination is well-tolerated with favorable efficacy. High-dimensional analyses at single-cell resolution represent promising approaches for identifying biomarkers of response and mechanisms of resistance in prospective clinical studies.

In Vitro Evolution Reveals a Single Mutation as Sole Source of Src-Family Kinase C-Helix-out Inhibitor Resistance.

Understanding cancer cell drug resistance to protein-tyrosine kinase inhibitors, which often arises from acquired mutations in the target kinase, is central to the development of more durable therapies. Experimental systems that reveal potential paths to resistance for a given inhibitor and kinase target have an important role in preclinical development of kinase inhibitor drugs. Here, we employed a codon mutagenesis strategy to define the mutational landscape of acquired resistance in HCK, a member of the SRC tyrosine kinase family and therapeutic target in acute myeloid leukemia (AML). Using PCR-based saturation mutagenesis, we created a cDNA library designed to replace each codon in the HCK open reading frame with all possible codons. This HCK mutant library was used to transform Rat-2 fibroblasts, followed by selection for resistant colonies with A-419259, a pyrrolopyrimidine HCK inhibitor and drug lead for AML. X-ray crystallography has shown that A-419259 binding induces outward rotation of the kinase domain αC-helix, a conformation incompatible with phosphotransfer. Remarkably, only a single resistance mutation evolved during A-419259 selection: histidine substitution for threonine at the gatekeeper position in the kinase domain. Deep sequencing confirmed representation of nearly all other missense mutations across the entire HCK open reading frame. This observation suggests that A-419259 and other C-helix-out Src-family kinase inhibitors may have a narrow path to acquired resistance in the context of AML cases where Hck is an oncogenic driver.

Expression of myeloid Src-family kinases is associated with poor prognosis in AML and influences Flt3-ITD kinase inhibitor acquired resistance.

Unregulated protein-tyrosine kinase signaling is a common feature of AML, often involving mutations in Flt3 and overexpression of myeloid Src-family kinases (Hck, Fgr, Lyn). Here we show that high-level expression of these Src kinases predicts poor survival in a large cohort of AML patients. To test the therapeutic benefit of Flt3 and Src-family kinase inhibition, we used the pyrrolopyrimidine kinase inhibitor A-419259. This compound potently inhibits Hck, Fgr, and Lyn as well as Flt3 bearing an activating internal tandem duplication (ITD). Flt3-ITD expression sensitized human TF-1 myeloid cells to growth arrest by A-419259, supporting direct action on the Flt3-ITD kinase domain. Cells transformed with the Flt3-ITD mutants D835Y and F691L were resistant to A-419259, while co-expression of Hck or Fgr restored inhibitor sensitivity to Flt3-ITD D835Y. Conversely, Hck and Fgr mutants with engineered A-419259 resistance mutations decreased sensitivity of TF-1/Flt3-ITD cells. To investigate de novo resistance mechanisms, A-419259-resistant Flt3-ITD+ AML cell populations were derived via long-term dose escalation. Whole exome sequencing identified a distinct Flt3-ITD kinase domain mutation (N676S/T) among all A-419259 target kinases in each of six independent resistant cell populations. These studies show that Hck and Fgr expression influences inhibitor sensitivity and the pathway to acquired resistance in Flt3-ITD+ AML.

MiR-146a wild-type 3' sequence identity is dispensable for proper innate immune function in vivo.

The prevailing model of microRNA function is that the "seed region" (nt 2-8) is sufficient to mediate target recognition and repression. However, numerous recent studies have challenged this model, either by demonstrating extensive 3' pairing between physically defined miRNA-mRNA pairs or by showing in Caenorhabditis elegans that disrupted 3' pairing can result in impaired function in vivo. To test the importance of miRNA 3' pairing in a mammalian system in vivo, we engineered a mutant murine mir-146a allele in which the 5' half of the mature microRNA retains its wild-type sequence, but the 3' half's sequence has been altered to robustly disrupt predicted pairing to this latter region. Mice homozygous or hemizygous for this mutant allele are phenotypically indistinguishable from wild-type controls and do not recapitulate any of the immunopathology previously described for mir-146a-null mice. Our results indicate that 3' pairing is dispensable for the established myeloid function of this key mammalian microRNA.

The Notch signaling pathway promotes basophil responses during helminth-induced type 2 inflammation.

Type 2 inflammation drives the clearance of gastrointestinal helminth parasites, which infect over two billion people worldwide. Basophils are innate immune cells that support host-protective type 2 inflammation during murine infection with the helminth Trichuris muris However, the mechanisms required for basophil function and gene expression regulation in this context remain unclear. We show that during T. muris infection, basophils localized to the intestine and up-regulated Notch receptor expression, rendering them sensitive to Notch signals that rapidly regulate gene expression programs. In vitro, Notch inhibition limited basophil cytokine production in response to cytokine stimulation. Basophil-intrinsic Notch signaling was required for T. muris-elicited changes in genome-wide basophil transcriptional programs. Mice lacking basophil-intrinsic functional Notch signaling had impaired worm clearance, decreased intestinal type 2 inflammation, altered basophil localization in the intestine, and decreased CD4+ T helper 2 cell responses following infection. These findings demonstrate that Notch is required for basophil gene expression and effector function associated with helminth expulsion during type 2 inflammation.