TCF1-LEF1 co-expression identifies a multipotent progenitor cell (TH2-MPP) across human allergic diseases.

Repetitive exposure to antigen in chronic infection and cancer drives T cell exhaustion, limiting adaptive immunity. In contrast, aberrant, sustained T cell responses can persist over decades in human allergic disease. To understand these divergent outcomes, we employed bioinformatic, immunophenotyping and functional approaches with human diseased tissues, identifying an abundant population of type 2 helper T (TH2) cells with co-expression of TCF7 and LEF1, and features of chronic activation. These cells, which we termed TH2-multipotent progenitors (TH2-MPP) could self-renew and differentiate into cytokine-producing effector cells, regulatory T (Treg) cells and follicular helper T (TFH) cells. Single-cell T-cell-receptor lineage tracing confirmed lineage relationships between TH2-MPP, TH2 effectors, Treg cells and TFH cells. TH2-MPP persisted despite in vivo IL-4 receptor blockade, while thymic stromal lymphopoietin (TSLP) drove selective expansion of progenitor cells and rendered them insensitive to glucocorticoid-induced apoptosis in vitro. Together, our data identify TH2-MPP as an aberrant T cell population with the potential to sustain type 2 inflammation and support the paradigm that chronic T cell responses can be coordinated over time by progenitor cells.

Increased glycolysis and cellular crosstalk in eosinophilic chronic rhinosinusitis with nasal polyps.

Introduction: Chronic rhinosinusitis (CRS) is a chronic inflammatory disease of the sinonasal mucosa with distinct endotypes including type 2 (T2) high eosinophilic CRS with nasal polyps (eCRSwNP), T2 low non-eosinophilic CRS with nasal polyps (neCRSwNP), and CRS without nasal polyps (CRSsNP). Methods: Given the heterogeneity of disease, we hypothesized that assessment of single cell RNA sequencing (scRNA-seq) across this spectrum of disease would reveal connections between infiltrating and activated immune cells and the epithelial and stromal populations that reside in sinonasal tissue. Results: Here we find increased expression of genes encoding glycolytic enzymes in epithelial cells (EpCs), stromal cells, and memory T-cell subsets from patients with eCRSwNP, as compared to healthy controls. In basal EpCs, this is associated with a program of cell motility and Rho GTPase effector expression. Across both stromal and immune subsets, glycolytic programming was associated with extracellular matrix interactions, proteoglycan generation, and collagen formation. Furthermore, we report increased cell-cell interactions between EpCs and stromal/immune cells in eCRSwNP compared to healthy control tissue, and we nominate candidate receptor-ligand pairs that may drive tissue remodeling. Discussion: These findings support a role for glycolytic reprograming in T2-elicited tissue remodeling and implicate increased cellular crosstalk in eCRSwNP.

Integrating circulating T follicular memory cells and autoantibody repertoires for characterization of autoimmune disorders.

Introduction: Autoimmune diseases are heterogeneous and often lack specific or sensitive diagnostic tests. Increased percentages of CD4+CXCR5+PD1+ circulating T follicular helper (cTfh) cells and skewed distributions of cTfh subtypes have been associated with autoimmunity. However, cTfh cell percentages can normalize with immunomodulatory treatment despite persistent disease activity, indicating the need for identifying additional cellular and/or serologic features correlating with autoimmunity. Methods: The cohort included 50 controls and 56 patients with autoimmune cytopenias, gastrointestinal, pulmonary, and/or neurologic autoimmune disease. Flow cytometry was used to measure CD4+CXCR5+ T cell subsets expressing the chemokine receptors CXCR3 and/or CCR6: CXCR3+CCR6- Type 1, CXCR3-CCR6- Type 2, CXCR3+CCR6+ Type 1/17, and CXCR3- CCR6+ Type 17 T cells. IgG and IgA autoantibodies were quantified using a microarray featuring 1616 full-length, conformationally intact protein antigens. The 97.5th percentile in the control cohort defined normal limits for T cell subset percentages and total number (burden) of autoantibodies. Results: This study focused on CD4+CXCR5+ T cells because CXCR5 upregulation occurs after cognate T-B cell interactions characteristic of autoimmune diseases. We refer to these cells as circulating T follicular memory (cTfm) cells to acknowledge the dynamic nature of antigen-experienced CXCR5+ T cells, which encompass progenitors of cTfh or Tfh cells as well as early effector memory T cells that have not yet lost CXCR5. Compared to controls, 57.1% of patients had increased CXCR5+CXCR3+CCR6+ cTfm1/17 and 25% had increased CXCR5+CXCR3-CCR6+ cTfm17 cell percentages. Patients had significantly more diverse IgG and IgA autoantibodies than controls and 44.6% had an increased burden of autoantibodies of either isotype. Unsupervised autoantibody clustering identified three clusters of patients with IgG autoantibody profiles distinct from those of controls, enriched for patients with active autoimmunity and monogenic diseases. An increased percentage of cTfm17 cells was most closely associated with an increased burden of high-titer IgG and IgA autoantibodies. A composite measure integrating increased cTfm1/17, cTfm17, and high-titer IgG and/or IgA autoantibodies had 91.1% sensitivity and 90.9% specificity for identifying patients with autoimmunity. Percentages of cTfm1/17 and cTfm17 percentages and numbers of high-titer autoantibodies in patients receiving immunomodulatory treatment did not differ from those in untreated patients, thus suggesting that measurements of cTfm can complement measurements of other cellular markers affected by treatment. Conclusions: This study highlights two new approaches for assessing autoimmunity: measuring CD4+CXCR5+ cTfm subsets as well as total burden of autoantibodies. Our findings suggest that these approaches are particularly relevant to patients with rare autoimmune disorders for whom target antigens and prognosis are often unknown.

Non-coding autoimmune risk variant defines role for ICOS in T peripheral helper cell development.

Fine-mapping and functional studies implicate rs117701653, a non-coding single nucleotide polymorphism in the CD28/CTLA4/ICOS locus, as a risk variant for rheumatoid arthritis and type 1 diabetes. Here, using DNA pulldown, mass spectrometry, genome editing and eQTL analysis, we establish that the disease-associated risk allele is functional, reducing affinity for the inhibitory chromosomal regulator SMCHD1 to enhance expression of inducible T-cell costimulator (ICOS) in memory CD4+ T cells from healthy donors. Higher ICOS expression is paralleled by an increase in circulating T peripheral helper (Tph) cells and, in rheumatoid arthritis patients, of blood and joint fluid Tph cells as well as circulating plasmablasts. Correspondingly, ICOS ligation and carriage of the rs117701653 risk allele accelerate T cell differentiation into CXCR5-PD-1high Tph cells producing IL-21 and CXCL13. Thus, mechanistic dissection of a functional non-coding variant in human autoimmunity discloses a previously undefined pathway through which ICOS regulates Tph development and abundance.

Rhinovirus infection of airway epithelial cells uncovers the non-ciliated subset as a likely driver of genetic susceptibility to childhood-onset asthma.

Asthma is a complex disease caused by genetic and environmental factors. Epidemiological studies have shown that in children, wheezing during rhinovirus infection (a cause of the common cold) is associated with asthma development during childhood. This has led scientists to hypothesize there could be a causal relationship between rhinovirus infection and asthma or that RV-induced wheezing identifies individuals at increased risk for asthma development. However, not all children who wheeze when they have a cold develop asthma. Genome-wide association studies (GWAS) have identified hundreds of genetic variants contributing to asthma susceptibility, with the vast majority of likely causal variants being non-coding. Integrative analyses with transcriptomic and epigenomic datasets have indicated that T cells drive asthma risk, which has been supported by mouse studies. However, the datasets ascertained in these integrative analyses lack airway epithelial cells. Furthermore, large-scale transcriptomic T cell studies have not identified the regulatory effects of most non-coding risk variants in asthma GWAS, indicating there could be additional cell types harboring these "missing regulatory effects". Given that airway epithelial cells are the first line of defense against rhinovirus, we hypothesized they could be mediators of genetic susceptibility to asthma. Here we integrate GWAS data with transcriptomic datasets of airway epithelial cells subject to stimuli that could induce activation states relevant to asthma. We demonstrate that epithelial cultures infected with rhinovirus significantly upregulate childhood-onset asthma-associated genes. We show that this upregulation occurs specifically in non-ciliated epithelial cells. This enrichment for genes in asthma risk loci, or 'asthma heritability enrichment' is also significant for epithelial genes upregulated with influenza infection, but not with SARS-CoV-2 infection or cytokine activation. Additionally, cells from patients with asthma showed a stronger heritability enrichment compared to cells from healthy individuals. Overall, our results suggest that rhinovirus infection is an environmental factor that interacts with genetic risk factors through non-ciliated airway epithelial cells to drive childhood-onset asthma.

Understanding Spondyloarthritis Pathogenesis: The Promise of Single-Cell Profiling.

PURPOSE OF REVIEW: Single-cell profiling, either in suspension or within the tissue context, is a rapidly evolving field. The purpose of this review is to outline recent advancements and emerging trends with a specific focus on studies in spondyloarthritis. RECENT FINDINGS: The introduction of sequencing-based approaches for the quantification of RNA, protein, or epigenetic modifications at single-cell resolution has provided a major boost to discovery-driven research. Fluorescent flow cytometry, mass cytometry, and image-based cytometry continue to evolve. Spatial transcriptomics and imaging mass cytometry have extended high-dimensional analysis to cells in tissues. Applications in spondyloarthritis include the indexing and functional characterization of cells, discovery of disease-associated cell states, and identification of signatures associated with therapeutic responses. Single-cell TCR-seq has provided evidence for clonal expansion of CD8+ T cells in spondyloarthritis. The use of single-cell profiling approaches in spondyloarthritis research is still in its early stages. Challenges include high cost and limited availability of diseased tissue samples. To harness the full potential of the rapidly expanding technical capabilities, large-scale collaborative efforts are imperative.

FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases.

The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

Blood immunophenotyping identifies distinct kidney histopathology and outcomes in patients with lupus nephritis.

Lupus nephritis (LN) is a frequent manifestation of systemic lupus erythematosus, and fewer than half of patients achieve complete renal response with standard immunosuppressants. Identifying non-invasive, blood-based pathologic immune alterations associated with renal injury could aid therapeutic decisions. Here, we used mass cytometry immunophenotyping of peripheral blood mononuclear cells in 145 patients with biopsy-proven LN and 40 healthy controls to evaluate the heterogeneity of immune activation in patients with LN and to identify correlates of renal parameters and treatment response. Unbiased analysis identified 3 immunologically distinct groups of patients with LN that were associated with different patterns of histopathology, renal cell infiltrates, urine proteomic profiles, and treatment response at one year. Patients with enriched circulating granzyme B+ T cells at baseline showed more severe disease and increased numbers of activated CD8 T cells in the kidney, yet they had the highest likelihood of treatment response. A second group characterized primarily by a high type I interferon signature had a lower likelihood of response to therapy, while a third group appeared immunologically inactive by immunophenotyping at enrollment but with chronic renal injuries. Main immune profiles could be distilled down to 5 simple cytometric parameters that recapitulate several of the associations, highlighting the potential for blood immune profiling to translate to clinically useful non-invasive metrics to assess immune-mediated disease in LN.

High-throughput identification of functional regulatory SNPs in systemic lupus erythematosus.

Genome-wide association studies implicate multiple loci in risk for systemic lupus erythematosus (SLE), but few contain exonic variants, rendering systematic identification of non-coding variants essential to decoding SLE genetics. We utilized SNP-seq and bioinformatic enrichment to interrogate 2180 single-nucleotide polymorphisms (SNPs) from 87 SLE risk loci for potential binding of transcription factors and related proteins from B cells. 52 SNPs that passed initial screening were tested by electrophoretic mobility shift and luciferase reporter assays. To validate the approach, we studied rs2297550 in detail, finding that the risk allele enhanced binding to the transcription factor Ikaros (IKZF1), thereby modulating expression of IKBKE. Correspondingly, primary cells from genotyped healthy donors bearing the risk allele expressed higher levels of the interferon / NF-κB regulator IKKϵ. Together, these findings define a set of likely functional non-coding lupus risk variants and identify a new regulatory pathway involving rs2297550, Ikaros, and IKKϵ implicated by human genetics in risk for SLE.

PPARγ agonist treatment reduces fibroadipose tissue in secondary lymphedema by exhausting fibroadipogenic PDGFRα+ mesenchymal cells.

Secondary lymphedema occurs in up to 20% of patients after lymphadenectomy performed for the surgical management of tumors involving the breast, prostate, uterus, and skin. Patients develop progressive edema of the affected extremity due to retention of protein-rich lymphatic fluid. Despite compression therapy, patients progress to chronic lymphedema in which noncompressible fibrosis and adipose tissue are deposited within the extremity. The presence of fibrosis led to our hypothesis that rosiglitazone, a PPARγ agonist that inhibits fibrosis, would reduce fibrosis in a mouse model of secondary lymphedema after hind limb lymphadenectomy. In vivo, rosiglitazone reduced fibrosis in the hind limb after lymphadenectomy. Our findings verified that rosiglitazone reestablished the adipogenic features of TGF-ß1-treated mesenchymal cells in vitro. Despite this, rosiglitazone led to a reduction in adipose tissue deposition. Single-cell RNA-Seq data obtained from human tissues and flow cytometric and histological evaluation of mouse tissues demonstrated increased presence of PDGFRα+ cells in lymphedema; human tissue analysis verified these cells have the capacity for adipogenic and fibrogenic differentiation. Upon treatment with rosiglitazone, we noted a reduction in the overall quantity of PDGFRα+ cells and LipidTOX+ cells. Our findings provide a framework for treating secondary lymphedema as a condition of fibrosis and adipose tissue deposition, both of which, paradoxically, can be prevented with a pro-adipogenic agent.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution.

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues. To mitigate technical confounding, we developed scHLApers, a pipeline to accurately quantify single-cell HLA expression using personalized reference genomes. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B and T cells. For example, a T cell HLA-DQA1 eQTL ( rs3104371 ) is strongest in cytotoxic cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Identification of a regulatory pathway governing TRAF1 via an arthritis-associated non-coding variant.

TRAF1/C5 was among the first loci shown to confer risk for inflammatory arthritis in the absence of an associated coding variant, but its genetic mechanism remains undefined. Using Immunochip data from 3,939 patients with juvenile idiopathic arthritis (JIA) and 14,412 control individuals, we identified 132 plausible common non-coding variants, reduced serially by single-nucleotide polymorphism sequencing (SNP-seq), electrophoretic mobility shift, and luciferase studies to the single variant rs7034653 in the third intron of TRAF1. Genetically manipulated experimental cells and primary monocytes from genotyped donors establish that the risk G allele reduces binding of Fos-related antigen 2 (FRA2), encoded by FOSL2, resulting in reduced TRAF1 expression and enhanced tumor necrosis factor (TNF) production. Conditioning on this JIA variant eliminated attributable risk for rheumatoid arthritis, implicating a mechanism shared across the arthritis spectrum. These findings reveal that rs7034653, FRA2, and TRAF1 mediate a pathway through which a non-coding functional variant drives risk of inflammatory arthritis in children and adults.

Deconstruction of rheumatoid arthritis synovium defines inflammatory subtypes.

Rheumatoid arthritis is a prototypical autoimmune disease that causes joint inflammation and destruction1. There is currently no cure for rheumatoid arthritis, and the effectiveness of treatments varies across patients, suggesting an undefined pathogenic diversity1,2. Here, to deconstruct the cell states and pathways that characterize this pathogenic heterogeneity, we profiled the full spectrum of cells in inflamed synovium from patients with rheumatoid arthritis. We used multi-modal single-cell RNA-sequencing and surface protein data coupled with histology of synovial tissue from 79 donors to build single-cell atlas of rheumatoid arthritis synovial tissue that includes more than 314,000 cells. We stratified tissues into six groups, referred to as cell-type abundance phenotypes (CTAPs), each characterized by selectively enriched cell states. These CTAPs demonstrate the diversity of synovial inflammation in rheumatoid arthritis, ranging from samples enriched for T and B cells to those largely lacking lymphocytes. Disease-relevant cell states, cytokines, risk genes, histology and serology metrics are associated with particular CTAPs. CTAPs are dynamic and can predict treatment response, highlighting the clinical utility of classifying rheumatoid arthritis synovial phenotypes. This comprehensive atlas and molecular, tissue-based stratification of rheumatoid arthritis synovial tissue reveal new insights into rheumatoid arthritis pathology and heterogeneity that could inform novel targeted treatments.

Increased epithelial mTORC1 activity in chronic rhinosinusitis with nasal polyps.

Background: The airway epithelium plays a central role in the pathogenesis of chronic respiratory diseases such as asthma and chronic rhinosinusitis with nasal polyps (CRSwNP), but the mechanisms by which airway epithelial cells (EpCs) maintain inflammation are poorly understood. Objective: We hypothesized that transcriptomic assessment of sorted airway EpCs across the spectrum of differentiation would allow us to define mechanisms by which EpCs perpetuate airway inflammation. Methods: Ethmoid sinus EpCs from adult patients with CRS were sorted into 3 subsets, bulk RNA sequenced, and analyzed for differentially expressed genes and pathways. Single cell RNA-seq (scRNA-seq) datasets from eosinophilic and non-eosinophilic CRSwNP and bulk RNA-seq of EpCs from mild/moderate and severe asthma were assessed. Immunofluorescent staining and ex vivo functional analysis of sinus EpCs were used to validate our findings. Results: Analysis within and across purified EpC subsets revealed an enrichment in glycolytic programming in CRSwNP vs CRSsNP. Correlation analysis identified mammalian target of rapamycin complex 1 (mTORC1) as a potential regulator of the glycolytic program and identified EpC expression of cytokines and wound healing genes as potential sequelae. mTORC1 activity was upregulated in CRSwNP, and ex vivo inhibition demonstrated that mTOR is critical for EpC generation of CXCL8, IL-33, and CXCL2. Across patient samples, the degree of glycolytic activity was associated with T2 inflammation in CRSwNP, and with both T2 and non-T2 inflammation in severe asthma. Conclusions: Together, these findings highlight a metabolic axis required to support epithelial generation of cytokines critical to both chronic T2 and non-T2 inflammation in CRSwNP and asthma.

Functional genomics implicates natural killer cells as potential key drivers in the pathogenesis of ankylosing spondylitis.

Objective: Multiple lines of evidence indicate that ankylosing spondylitis (AS) is a lymphocyte-driven disease. However, which lymphocyte populations are critical in AS pathogenesis is not known. In this study, we aimed to identify the key cell types mediating the genetic risk in AS using an unbiased integrative functional genomics approach. Methods: We integrated GWAS data with epigenomic and transcriptomic datasets of immune cells in healthy humans. To quantify enrichment of cell type-specific open chromatin regions or gene expression in AS risk loci, we used three published methods which have identified cell types for other diseases. Additionally, we performed co-localization analyses between GWAS risk loci and genetic variants associated with gene expression (eQTL) to find putative target genes of AS risk variants. Results: Natural killer (NK) cell-specific open chromatin regions are significantly enriched in heritability for AS, compared to other immune cell types such as T cells, B cells, and monocytes. This finding was consistent between two AS GWAS. Using RNA-seq data, we validated that genes in AS risk loci are enriched in NK cell-specific gene expression. Expression levels of AS-associated genes, such as RUNX3, TBX21, TNFRSF1A, and NPEPPS, were found to be highest in NK cells compared to five T cell subsets. Using the human Space-Time Gut Cell Atlas we found significant upregulation of AS-associated genes predominantly in NK cells. Co-localization analysis revealed four AS risk loci affecting regulation of candidate target genes in NK cells: two known loci, ERAP1 and TNFRSF1A, and two under-studied loci, ENTR1 (aka SDCCAG3) and B3GNT2. Conclusion: Our results point to NK cells as potential key drivers in the development of AS and highlight four putative target genes for functional follow-up in NK cells.

History of tuberculosis disease is associated with genetic regulatory variation in Peruvians.

A quarter of humanity is estimated to be latently infected with Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n=63) or did not progress to TB (controls, n=63). Transcriptomic profiling of monocyte-derived dendritic cells (DCs) and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Five genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution.

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation, and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here, we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues, using personalized reference genomes to mitigate technical confounding. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B, and T cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Type 2 inflammation drives an airway basal stem cell program through insulin receptor substrate signaling.

BACKGROUND: Chronic rhinosinusitis with nasal polyposis (CRSwNP) is a type 2 (T2) inflammatory disease associated with an increased number of airway basal cells (BCs). Recent studies have identified transcriptionally distinct BCs, but the molecular pathways that support or inhibit human BC proliferation and differentiation are largely unknown. OBJECTIVE: We sought to determine the role of T2 cytokines in regulating airway BCs. METHODS: Single-cell and bulk RNA sequencing of sinus and lung airway epithelial cells was analyzed. Human sinus BCs were stimulated with IL-4 and IL-13 in the presence and absence of inhibitors of IL-4R signaling. Confocal analysis of human sinus tissue and murine airway was performed. Murine BC subsets were sorted for RNA sequencing and functional assays. Fate labeling was performed in a murine model of tracheal injury and regeneration. RESULTS: Two subsets of BCs were found in human and murine respiratory mucosa distinguished by the expression of basal cell adhesion molecule (BCAM). BCAM expression identifies airway stem cells among P63+KRT5+NGFR+ BCs. In the sinonasal mucosa, BCAMhi BCs expressing TSLP, IL33, CCL26, and the canonical BC transcription factor TP63 are increased in patients with CRSwNP. In cultured BCs, IL-4/IL-13 increases the expression of BCAM and TP63 through an insulin receptor substrate-dependent signaling pathway that is increased in CRSwNP. CONCLUSIONS: These findings establish BCAM as a marker of airway stem cells among the BC pool and demonstrate that airway epithelial remodeling in T2 inflammation extends beyond goblet cell metaplasia to the support of a BC stem state poised to perpetuate inflammation.

Comparison between qPCR and RNA-seq reveals challenges of quantifying HLA expression.

Human leukocyte antigen (HLA) class I and II loci are essential elements of innate and acquired immunity. Their functions include antigen presentation to T cells leading to cellular and humoral immune responses, and modulation of NK cells. Their exceptional influence on disease outcome has now been made clear by genome-wide association studies. The exons encoding the peptide-binding groove have been the main focus for determining HLA effects on disease susceptibility/pathogenesis. However, HLA expression levels have also been implicated in disease outcome, adding another dimension to the extreme diversity of HLA that impacts variability in immune responses across individuals. To estimate HLA expression, immunogenetic studies traditionally rely on quantitative PCR (qPCR). Adoption of alternative high-throughput technologies such as RNA-seq has been hampered by technical issues due to the extreme polymorphism at HLA genes. Recently, however, multiple bioinformatic methods have been developed to accurately estimate HLA expression from RNA-seq data. This opens an exciting opportunity to quantify HLA expression in large datasets but also brings questions on whether RNA-seq results are comparable to those by qPCR. In this study, we analyze three classes of expression data for HLA class I genes for a matched set of individuals: (a) RNA-seq, (b) qPCR, and (c) cell surface HLA-C expression. We observed a moderate correlation between expression estimates from qPCR and RNA-seq for HLA-A, -B, and -C (0.2 ≤ rho ≤ 0.53). We discuss technical and biological factors which need to be accounted for when comparing quantifications for different molecular phenotypes or using different techniques.