Tissue signals imprint ILC2 identity with anticipatory function.

Group 2 innate lymphoid cells (ILC2s) are distributed systemically and produce type 2 cytokines in response to a variety of stimuli, including the epithelial cytokines interleukin (IL)-25, IL-33, and thymic stromal lymphopoietin (TSLP). Transcriptional profiling of ILC2s from different tissues, however, grouped ILC2s according to their tissue of origin, even in the setting of combined IL-25-, IL-33-receptor-, and TSLP-receptor-deficiency. Single-cell profiling confirmed a tissue-organizing transcriptome and identified ILC2 subsets expressing distinct activating receptors, including the major subset of skin ILC2s, which were activated preferentially by IL-18. Tissue ILC2 subsets were unaltered in number and expression in germ-free mice, suggesting that endogenous, tissue-derived signals drive the maturation of ILC2 subsets by controlling expression of distinct patterns of activating receptors, thus anticipating tissue-specific perturbations occurring later in life.

Genomic characterization of the adolescent idiopathic scoliosis-associated transcriptome and regulome.

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is the most common pediatric musculoskeletal disorder, affecting ~3% of the population worldwide. However, its genetic bases and tissues of origin remain largely unknown. Several genome-wide association studies (GWAS) have implicated nucleotide variants in non-coding sequences that control genes with important roles in cartilage, muscle, bone, connective tissue and intervertebral disks (IVDs) as drivers of AIS susceptibility. Here, we set out to define the expression of AIS-associated genes and active regulatory elements by performing RNA-seq and chromatin immunoprecipitation-sequencing against H3 lysine 27 acetylation in these tissues in mouse and human. Our study highlights genetic pathways involving AIS-associated loci that regulate chondrogenesis, IVD development and connective tissue maintenance and homeostasis. In addition, we identify thousands of putative AIS-associated regulatory elements which may orchestrate tissue-specific expression in musculoskeletal tissues of the spine. Quantification of enhancer activity of several candidate regulatory elements from our study identifies three functional enhancers carrying AIS-associated GWAS SNPs at the ADGRG6 and BNC2 loci. Our findings provide a novel genome-wide catalog of AIS-relevant genes and regulatory elements and aid in the identification of novel targets for AIS causality and treatment.

On the cross-population generalizability of gene expression prediction models.

The genetic control of gene expression is a core component of human physiology. For the past several years, transcriptome-wide association studies have leveraged large datasets of linked genotype and RNA sequencing information to create a powerful gene-based test of association that has been used in dozens of studies. While numerous discoveries have been made, the populations in the training data are overwhelmingly of European descent, and little is known about the generalizability of these models to other populations. Here, we test for cross-population generalizability of gene expression prediction models using a dataset of African American individuals with RNA-Seq data in whole blood. We find that the default models trained in large datasets such as GTEx and DGN fare poorly in African Americans, with a notable reduction in prediction accuracy when compared to European Americans. We replicate these limitations in cross-population generalizability using the five populations in the GEUVADIS dataset. Via realistic simulations of both populations and gene expression, we show that accurate cross-population generalizability of transcriptome prediction only arises when eQTL architecture is substantially shared across populations. In contrast, models with non-identical eQTLs showed patterns similar to real-world data. Therefore, generating RNA-Seq data in diverse populations is a critical step towards multi-ethnic utility of gene expression prediction.

The Type 2 Asthma Mediator IL-13 Inhibits Severe Acute Respiratory Syndrome Coronavirus 2 Infection of Bronchial Epithelium.

Asthma is associated with chronic changes in the airway epithelium, a key target of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Many epithelial changes, including goblet cell metaplasia, are driven by the type 2 cytokine IL-13, but the effects of IL-13 on SARS-CoV-2 infection are unknown. We found that IL-13 stimulation of differentiated human bronchial epithelial cells (HBECs) cultured at air-liquid interface reduced viral RNA recovered from SARS-CoV-2-infected cells and decreased double-stranded RNA, a marker of viral replication, to below the limit of detection in our assay. An intact mucus gel reduced SARS-CoV-2 infection of unstimulated cells, but neither a mucus gel nor SPDEF, which is required for goblet cell metaplasia, were required for the antiviral effects of IL-13. Bulk RNA sequencing revealed that IL-13 regulated 41 of 332 (12%) mRNAs encoding SARS-CoV-2-associated proteins that were detected in HBECs (>1.5-fold change; false discovery rate < 0.05). Although both IL-13 and IFN-α each inhibit SARS-CoV-2 infection, their transcriptional effects differed markedly. Single-cell RNA sequencing revealed cell type-specific differences in SARS-CoV-2-associated gene expression and IL-13 responses. Many IL-13-induced gene expression changes were seen in airway epithelium from individuals with type 2 asthma and chronic obstructive pulmonary disease. IL-13 effects on airway epithelial cells may protect individuals with type 2 asthma from COVID-19 and could lead to identification of novel strategies for reducing SARS-CoV-2 infection.

Flow-Cytometric Analysis and Purification of Airway Epithelial-Cell Subsets.

The human airway epithelium is essential in homeostasis, and epithelial dysfunction contributes to chronic airway disease. Development of flow-cytometric methods to characterize subsets of airway epithelial cells will enable further dissection of airway epithelial biology. Leveraging single-cell RNA-sequencing data in combination with known cell type-specific markers, we developed panels of antibodies to characterize and isolate the major airway epithelial subsets (basal, ciliated, and secretory cells) from human bronchial epithelial-cell cultures. We also identified molecularly distinct subpopulations of secretory cells and demonstrated cell subset-specific expression of low-abundance transcripts and microRNAs that are challenging to analyze with current single-cell RNA-sequencing methods. These new tools will be valuable for analyzing and separating airway epithelial subsets and interrogating airway epithelial biology.

Whole-Genome Sequencing Identifies Novel Functional Loci Associated with Lung Function in Puerto Rican Youth.

Rationale: Puerto Ricans have the highest childhood asthma prevalence in the United States (23.6%); however, the etiology is uncertain.Objectives: In this study, we sought to uncover the genetic architecture of lung function in Puerto Rican youth with and without asthma who were recruited from the island (n = 836).Methods: We used admixture-mapping and whole-genome sequencing data to discover genomic regions associated with lung function. Functional roles of the prioritized candidate SNPs were examined with chromatin immunoprecipitation sequencing, RNA sequencing, and expression quantitative trait loci data.Measurements and Main Results: We discovered a genomic region at 1q32 that was significantly associated with a 0.12-L decrease in the lung volume of exhaled air (95% confidence interval, -0.17 to -0.07; P = 6.62 × 10-8) with each allele of African ancestry. Within this region, two SNPs were expression quantitative trait loci of TMEM9 in nasal airway epithelial cells and MROH3P in esophagus mucosa. The minor alleles of these SNPs were associated with significantly decreased lung function and decreased TMEM9 gene expression. Another admixture-mapping peak was observed on chromosome 5q35.1, indicating that each Native American ancestry allele was associated with a 0.15-L increase in lung function (95% confidence interval, 0.08-0.21; P = 5.03 × 10-6). The region-based association tests identified four suggestive windows that harbored candidate rare variants associated with lung function.Conclusions: We identified common and rare genetic variants that may play a critical role in lung function among Puerto Rican youth. We independently validated an inflammatory pathway that could potentially be used to develop more targeted treatments and interventions for patients with asthma.

Comparative Genomic Characterization of the Multimammate Mouse Mastomys coucha.

Mastomys are the most widespread African rodent and carriers of various diseases such as the plague or Lassa virus. In addition, mastomys have rapidly gained a large number of mammary glands. Here, we generated a genome, variome, and transcriptomes for Mastomys coucha. As mastomys diverged at similar times from mouse and rat, we demonstrate their utility as a comparative genomic tool for these commonly used animal models. Furthermore, we identified over 500 mastomys accelerated regions, often residing near important mammary developmental genes or within their exons leading to protein sequence changes. Functional characterization of a noncoding mastomys accelerated region, located in the HoxD locus, showed enhancer activity in mouse developing mammary glands. Combined, our results provide genomic resources for mastomys and highlight their potential both as a comparative genomic tool and for the identification of mammary gland number determining factors.

Integrative approach identifies corticosteroid response variant in diverse populations with asthma.

BACKGROUND: Although inhaled corticosteroid (ICS) medication is considered the cornerstone treatment for patients with persistent asthma, few ICS pharmacogenomic studies have involved nonwhite populations. OBJECTIVE: We sought to identify genetic predictors of ICS response in multiple population groups with asthma. METHODS: The discovery group comprised African American participants from the Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-Ethnicity (SAPPHIRE) who underwent 6 weeks of monitored ICS therapy (n = 244). A genome-wide scan was performed to identify single nucleotide polymorphism (SNP) variants jointly associated (ie, the combined effect of the SNP and SNP × ICS treatment interaction) with changes in asthma control. Top associations were validated by assessing the joint association with asthma exacerbations in 3 additional groups: African Americans (n = 803 and n = 563) and Latinos (n = 1461). RNA sequencing data from 408 asthmatic patients and 405 control subjects were used to examine whether genotype was associated with gene expression. RESULTS: One variant, rs3827907, was significantly associated with ICS-mediated changes in asthma control in the discovery set (P = 7.79 × 10-8) and was jointly associated with asthma exacerbations in 3 validation cohorts (P = .023, P = .029, and P = .041). RNA sequencing analysis found the rs3827907 C-allele to be associated with lower RNASE2 expression (P = 6.10 × 10-4). RNASE2 encodes eosinophil-derived neurotoxin, and the rs3827907 C-allele appeared to particularly influence ICS treatment response in the presence of eosinophilic inflammation (ie, high pretreatment eosinophil-derived neurotoxin levels or blood eosinophil counts). CONCLUSION: We identified a variant, rs3827907, that appears to influence response to ICS treatment in multiple population groups and likely mediates its effect through eosinophils.

Cell fate analysis in fetal mouse lung reveals distinct pathways for TI and TII cell development.

Alveolar type I (TI) cells are large squamous cells that cover >95% of the internal surface area of the lung; type II (TII) cells are small cuboidal cells with distinctive intracellular surfactant storage organelles. Based on autoradiographic studies in the 1970s, the long-held paradigm of alveolar epithelial development has been a linear progression from undifferentiated progenitor cells through TII cells to TI cells. Subsequent data support the existence of more complex pathways. Recently, a bipotent TI/TII progenitor cell at embryonic day E18 has been inferred both from marker expression in developing airways and from statistical analyses of gene expression data obtained from single-lung embryonic cells. To study cell lineage directly by fate mapping, we developed new transgenic mouse models in which rtTA is driven either by the rat podoplanin or the mouse Sftpc gene to mark cells irreversibly in development. Using these models, we found two distinct lineage pathways. One pathway, evident as early as E12-15, is devoted almost exclusively to TI cell development; a second pathway gives rise predominantly to TII cells but also a subpopulation of TI cells. We have defined the molecular phenotypes of these distinct progenitor populations and have identified potential regulatory factors in TI and TII cell differentiation. By analyzing gene pathways in mature TI and TII cells, we identified potential novel functions of each cell type. These results provide novel insights into lung development and suggest a basis for testing strategies to promote alveolar differentiation and repair, including potential transplantation of lineage-specific progenitor cells.

Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma.

RATIONALE: Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response. OBJECTIVES: To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children. METHODS: We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR. MEASUREMENTS AND MAIN RESULTS: We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P < 3.53 × 10-7) and suggestive (P < 7.06 × 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and ß-adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR-associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole-genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age-matched and population-matched cohorts of sufficient sample size for replication of our common variant findings. CONCLUSIONS: The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations.

Bat Accelerated Regions Identify a Bat Forelimb Specific Enhancer in the HoxD Locus.

The molecular events leading to the development of the bat wing remain largely unknown, and are thought to be caused, in part, by changes in gene expression during limb development. These expression changes could be instigated by variations in gene regulatory enhancers. Here, we used a comparative genomics approach to identify regions that evolved rapidly in the bat ancestor, but are highly conserved in other vertebrates. We discovered 166 bat accelerated regions (BARs) that overlap H3K27ac and p300 ChIP-seq peaks in developing mouse limbs. Using a mouse enhancer assay, we show that five Myotis lucifugus BARs drive gene expression in the developing mouse limb, with the majority showing differential enhancer activity compared to the mouse orthologous BAR sequences. These include BAR116, which is located telomeric to the HoxD cluster and had robust forelimb expression for the M. lucifugus sequence and no activity for the mouse sequence at embryonic day 12.5. Developing limb expression analysis of Hoxd10-Hoxd13 in Miniopterus natalensis bats showed a high-forelimb weak-hindlimb expression for Hoxd10-Hoxd11, similar to the expression trend observed for M. lucifugus BAR116 in mice, suggesting that it could be involved in the regulation of the bat HoxD complex. Combined, our results highlight novel regulatory regions that could be instrumental for the morphological differences leading to the development of the bat wing.

Genomic Characterization of Metformin Hepatic Response.

Metformin is used as a first-line therapy for type 2 diabetes (T2D) and prescribed for numerous other diseases. However, its mechanism of action in the liver has yet to be characterized in a systematic manner. To comprehensively identify genes and regulatory elements associated with metformin treatment, we carried out RNA-seq and ChIP-seq (H3K27ac, H3K27me3) on primary human hepatocytes from the same donor treated with vehicle control, metformin or metformin and compound C, an AMP-activated protein kinase (AMPK) inhibitor (allowing to identify AMPK-independent pathways). We identified thousands of metformin responsive AMPK-dependent and AMPK-independent differentially expressed genes and regulatory elements. We functionally validated several elements for metformin-induced promoter and enhancer activity. These include an enhancer in an ataxia telangiectasia mutated (ATM) intron that has SNPs in linkage disequilibrium with a metformin treatment response GWAS lead SNP (rs11212617) that showed increased enhancer activity for the associated haplotype. Expression quantitative trait locus (eQTL) liver analysis and CRISPR activation suggest that this enhancer could be regulating ATM, which has a known role in AMPK activation, and potentially also EXPH5 and DDX10, its neighboring genes. Using ChIP-seq and siRNA knockdown, we further show that activating transcription factor 3 (ATF3), our top metformin upregulated AMPK-dependent gene, could have an important role in gluconeogenesis repression. Our findings provide a genome-wide representation of metformin hepatic response, highlight important sequences that could be associated with interindividual variability in glycemic response to metformin and identify novel T2D treatment candidates.

Genome-wide discovery of drug-dependent human liver regulatory elements.

Inter-individual variation in gene regulatory elements is hypothesized to play a causative role in adverse drug reactions and reduced drug activity. However, relatively little is known about the location and function of drug-dependent elements. To uncover drug-associated elements in a genome-wide manner, we performed RNA-seq and ChIP-seq using antibodies against the pregnane X receptor (PXR) and three active regulatory marks (p300, H3K4me1, H3K27ac) on primary human hepatocytes treated with rifampin or vehicle control. Rifampin and PXR were chosen since they are part of the CYP3A4 pathway, which is known to account for the metabolism of more than 50% of all prescribed drugs. We selected 227 proximal promoters for genes with rifampin-dependent expression or nearby PXR/p300 occupancy sites and assayed their ability to induce luciferase in rifampin-treated HepG2 cells, finding only 10 (4.4%) that exhibited drug-dependent activity. As this result suggested a role for distal enhancer modules, we searched more broadly to identify 1,297 genomic regions bearing a conditional PXR occupancy as well as all three active regulatory marks. These regions are enriched near genes that function in the metabolism of xenobiotics, specifically members of the cytochrome P450 family. We performed enhancer assays in rifampin-treated HepG2 cells for 42 of these sequences as well as 7 sequences that overlap linkage-disequilibrium blocks defined by lead SNPs from pharmacogenomic GWAS studies, revealing 15/42 and 4/7 to be functional enhancers, respectively. A common African haplotype in one of these enhancers in the GSTA locus was found to exhibit potential rifampin hypersensitivity. Combined, our results further suggest that enhancers are the predominant targets of rifampin-induced PXR activation, provide a genome-wide catalog of PXR targets and serve as a model for the identification of drug-responsive regulatory elements.

Differential expression of conserved and novel microRNAs during tail regeneration in the lizard Anolis carolinensis.

BACKGROUND: Lizards are evolutionarily the most closely related vertebrates to humans that can lose and regrow an entire appendage. Regeneration in lizards involves differential expression of hundreds of genes that regulate wound healing, musculoskeletal development, hormonal response, and embryonic morphogenesis. While microRNAs are able to regulate large groups of genes, their role in lizard regeneration has not been investigated. RESULTS: MicroRNA sequencing of green anole lizard (Anolis carolinensis) regenerating tail and associated tissues revealed 350 putative novel and 196 known microRNA precursors. Eleven microRNAs were differentially expressed between the regenerating tail tip and base during maximum outgrowth (25 days post autotomy), including miR-133a, miR-133b, and miR-206, which have been reported to regulate regeneration and stem cell proliferation in other model systems. Three putative novel differentially expressed microRNAs were identified in the regenerating tail tip. CONCLUSIONS: Differentially expressed microRNAs were identified in the regenerating lizard tail, including known regulators of stem cell proliferation. The identification of 3 putative novel microRNAs suggests that regulatory networks, either conserved in vertebrates and previously uncharacterized or specific to lizards, are involved in regeneration. These findings suggest that differential regulation of microRNAs may play a role in coordinating the timing and expression of hundreds of genes involved in regeneration.

Comparative proteomic analysis of tail regeneration in the green anole lizard, Anolis carolinensis.

As amniote vertebrates, lizards are the most closely related organisms to humans capable of appendage regeneration. Lizards can autotomize, or release their tails as a means of predator evasion, and subsequently regenerate a functional replacement. Green anoles (Anolis carolinensis) can regenerate their tails through a process that involves differential expression of hundreds of genes, which has previously been analyzed by transcriptomic and microRNA analysis. To investigate protein expression in regenerating tissue, we performed whole proteomic analysis of regenerating tail tip and base. This is the first proteomic data set available for any anole lizard. We identified a total of 2,646 proteins - 976 proteins only in the regenerating tail base, 796 only in the tail tip, and 874 in both tip and base. For over 90% of these proteins in these tissues, we were able to assign a clear orthology to gene models in either the Ensembl or NCBI databases. For 13 proteins in the tail base, 9 proteins in the tail tip, and 10 proteins in both regions, the gene model in Ensembl and NCBI matched an uncharacterized protein, confirming that these predictions are present in the proteome. Ontology and pathways analysis of proteins expressed in the regenerating tail base identified categories including actin filament-based process, ncRNA metabolism, regulation of phosphatase activity, small GTPase mediated signal transduction, and cellular component organization or biogenesis. Analysis of proteins expressed in the tail tip identified categories including regulation of organelle organization, regulation of protein localization, ubiquitin-dependent protein catabolism, small GTPase mediated signal transduction, morphogenesis of epithelium, and regulation of biological quality. These proteomic findings confirm pathways and gene families activated in tail regeneration in the green anole as well as identify uncharacterized proteins whose role in regrowth remains to be revealed. This study demonstrates the insights that are possible from the integration of proteomic and transcriptomic data in tail regrowth in the green anole, with potentially broader application to studies in other regenerative models.

Defective interferon signaling in the circulating monocytes of type 2 diabetic mice.

Type 2 diabetes mellitus (T2DM) is associated with poor outcome after stroke. Peripheral monocytes play a critical role in the secondary injury and recovery of damaged brain tissue after stroke, but the underlying mechanisms are largely unclear. To investigate transcriptome changes and molecular networks across monocyte subsets in response to T2DM and stroke, we performed single-cell RNA-sequencing (scRNAseq) from peripheral blood mononuclear cells and bulk RNA-sequencing from blood monocytes from four groups of adult mice, consisting of T2DM model db/db and normoglycemic control db/+ mice with or without ischemic stroke. Via scRNAseq we found that T2DM expands the monocyte population at the expense of lymphocytes, which was validated by flow cytometry. Among the monocytes, T2DM also disproportionally increased the inflammatory subsets with Ly6C+ and negative MHC class II expression (MO.6C+II-). Conversely, monocytes from control mice without stroke are enriched with steady-state classical monocyte subset of MO.6C+II+ but with the least percentage of MO.6C+II- subtype. Apart from enhancing inflammation and coagulation, enrichment analysis from both scRNAseq and bulk RNAseq revealed that T2DM specifically suppressed type-1 and type-2 interferon signaling pathways crucial for antigen presentation and the induction of ischemia tolerance. Preconditioning by lipopolysaccharide conferred neuroprotection against ischemic brain injury in db/+ but not in db/db mice and coincided with a lesser induction of brain Interferon-regulatory-factor-3 in the brains of the latter mice. Our results suggest that the increased diversity and altered transcriptome in the monocytes of T2DM mice underlie the worse stroke outcome by exacerbating secondary injury and potentiating stroke-induced immunosuppression. Significance Statement: The mechanisms involved in the detrimental diabetic effect on stroke are largely unclear. We show here, for the first time, that peripheral monocytes have disproportionally altered the subsets and changed transcriptome under diabetes and/or stroke conditions. Moreover, genes in the IFN-related signaling pathways are suppressed in the diabetic monocytes, which underscores the immunosuppression and impaired ischemic tolerance under the T2DM condition. Our data raise a possibility that malfunctioned monocytes may systemically and focally affect the host, leading to the poor outcome of diabetes in the setting of stroke. The results yield important clues to molecular mechanisms involved in the detrimental diabetic effect on stroke outcome.

Regulatory T cell-derived enkephalin imparts pregnancy-induced analgesia.

T cells have emerged as sex-dependent orchestrators of pain chronification but the sexually dimorphic mechanisms by which T cells control pain sensitivity is not resolved. Here, we demonstrate an influence of regulatory T cells (Tregs) on pain processing that is distinct from their canonical functions of immune regulation and tissue repair. Specifically, meningeal Tregs (mTregs) express the endogenous opioid, enkephalin, and mTreg-derived enkephalin exerts an antinociceptive action through a presynaptic opioid receptor signaling mechanism that is dispensable for immunosuppression. mTregs are both necessary and sufficient for suppressing mechanical pain sensitivity in female but not male mice. Notably, the mTreg modulation of pain thresholds depends on sex-hormones and expansion of enkephalinergic mTregs during gestation imparts a remarkable pregnancy-induced analgesia in a pre-existing, chronic, unremitting neuropathic pain model. These results uncover a fundamental sex-specific, pregnancy-pronounced, and immunologically-derived endogenous opioid circuit for nociceptive regulation with critical implications for pain biology and maternal health.

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.

Host-microbe multiomic profiling reveals age-dependent immune dysregulation associated with COVID-19 immunopathology.

Age is a major risk factor for severe coronavirus disease 2019 (COVID-19), yet the mechanisms behind this relationship have remained incompletely understood. To address this, we evaluated the impact of aging on host immune response in the blood and the upper airway, as well as the nasal microbiome in a prospective, multicenter cohort of 1031 vaccine-naïve patients hospitalized for COVID-19 between 18 and 96 years old. We performed mass cytometry, serum protein profiling, anti-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) antibody assays, and blood and nasal transcriptomics. We found that older age correlated with increased SARS-CoV-2 viral abundance upon hospital admission, delayed viral clearance, and increased type I interferon gene expression in both the blood and upper airway. We also observed age-dependent up-regulation of innate immune signaling pathways and down-regulation of adaptive immune signaling pathways. Older adults had lower naïve T and B cell populations and higher monocyte populations. Over time, older adults demonstrated a sustained induction of pro-inflammatory genes and serum chemokines compared with younger individuals, suggesting an age-dependent impairment in inflammation resolution. Transcriptional and protein biomarkers of disease severity differed with age, with the oldest adults exhibiting greater expression of pro-inflammatory genes and proteins in severe disease. Together, our study finds that aging is associated with impaired viral clearance, dysregulated immune signaling, and persistent and potentially pathologic activation of pro-inflammatory genes and proteins.

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.