Cell-type-specific and disease-associated expression quantitative trait loci in the human lung.

Common genetic variants confer substantial risk for chronic lung diseases, including pulmonary fibrosis. Defining the genetic control of gene expression in a cell-type-specific and context-dependent manner is critical for understanding the mechanisms through which genetic variation influences complex traits and disease pathobiology. To this end, we performed single-cell RNA sequencing of lung tissue from 66 individuals with pulmonary fibrosis and 48 unaffected donors. Using a pseudobulk approach, we mapped expression quantitative trait loci (eQTLs) across 38 cell types, observing both shared and cell-type-specific regulatory effects. Furthermore, we identified disease interaction eQTLs and demonstrated that this class of associations is more likely to be cell-type-specific and linked to cellular dysregulation in pulmonary fibrosis. Finally, we connected lung disease risk variants to their regulatory targets in disease-relevant cell types. These results indicate that cellular context determines the impact of genetic variation on gene expression and implicates context-specific eQTLs as key regulators of lung homeostasis and disease.

Locoregional delivery of IL-13Rα2-targeting CAR-T cells in recurrent high-grade glioma: a phase 1 trial.

Chimeric antigen receptor T cell (CAR-T) therapy is an emerging strategy to improve treatment outcomes for recurrent high-grade glioma, a cancer that responds poorly to current therapies. Here we report a completed phase I trial evaluating IL-13Rα2-targeted CAR-T cells in 65 patients with recurrent high-grade glioma, the majority being recurrent glioblastoma (rGBM). Primary objectives were safety and feasibility, maximum tolerated dose/maximum feasible dose and a recommended phase 2 dose plan. Secondary objectives included overall survival, disease response, cytokine dynamics and tumor immune contexture biomarkers. This trial evolved to evaluate three routes of locoregional T cell administration (intratumoral (ICT), intraventricular (ICV) and dual ICT/ICV) and two manufacturing platforms, culminating in arm 5, which utilized dual ICT/ICV delivery and an optimized manufacturing process. Locoregional CAR-T cell administration was feasible and well tolerated, and as there were no dose-limiting toxicities across all arms, a maximum tolerated dose was not determined. Probable treatment-related grade 3+ toxicities were one grade 3 encephalopathy and one grade 3 ataxia. A clinical maximum feasible dose of 200 × 106 CAR-T cells per infusion cycle was achieved for arm 5; however, other arms either did not test or achieve this dose due to manufacturing feasibility. A recommended phase 2 dose will be refined in future studies based on data from this trial. Stable disease or better was achieved in 50% (29/58) of patients, with two partial responses, one complete response and a second complete response after additional CAR-T cycles off protocol. For rGBM, median overall survival for all patients was 7.7 months and for arm 5 was 10.2 months. Central nervous system increases in inflammatory cytokines, including IFNγ, CXCL9 and CXCL10, were associated with CAR-T cell administration and bioactivity. Pretreatment intratumoral CD3 T cell levels were positively associated with survival. These findings demonstrate that locoregional IL-13Rα2-targeted CAR-T therapy is safe with promising clinical activity in a subset of patients. ClinicalTrials.gov Identifier: NCT02208362 .

Targeting Wnt signaling for improved glioma immunotherapy.

Introduction: Despite aggressive standard-of-care therapy, including surgery, radiation, and chemotherapy, glioblastoma recurrence is almost inevitable and uniformly lethal. Activation of glioma-intrinsic Wnt/ß-catenin signaling is associated with a poor prognosis and the proliferation of glioma stem-like cells, leading to malignant transformation and tumor progression. Impressive results in a subset of cancers have been obtained using immunotherapies including anti-CTLA4, anti-PD-1, and anti-PD-L1 or chimeric antigen receptor (CAR) T cell therapies. However, the heterogeneity of tumors, low mutational burden, single antigen targeting, and associated antigen escape contribute to non-responsiveness and potential tumor recurrence despite these therapeutic efforts. In the current study, we determined the effects of the small molecule, highly specific Wnt/CBP (CREB Binding Protein)/ß-catenin antagonist ICG-001, on glioma tumor cells and the tumor microenvironment (TME)-including its effect on immune cell infiltration, blood vessel decompression, and metabolic changes. Methods: Using multiple glioma patient-derived xenografts cell lines and murine tumors (GL261, K-Luc), we demonstrated in vitro cytostatic effects and a switch from proliferation to differentiation after treatment with ICG-001. Results: In these glioma cell lines, we further demonstrated that ICG-001 downregulated the CBP/ß-catenin target gene Survivin/BIRC5-a hallmark of Wnt/CBP/ß-catenin inhibition. We found that in a syngeneic mouse model of glioma (K-luc), ICG-001 treatment enhanced tumor infiltration by CD3+ and CD8+ cells with increased expression of the vascular endothelial marker CD31 (PECAM-1). We also observed differential gene expression and induced immune cell infiltration in tumors pretreated with ICG-001 and then treated with CAR T cells as compared with single treatment groups or when ICG-001 treatment was administered after CAR T cell therapy. Discussion: We conclude that specific Wnt/CBP/ß-catenin antagonism results in pleotropic changes in the glioma TME, including glioma stem cell differentiation, modulation of the stroma, and immune cell activation and recruitment, thereby suggesting a possible role for enhancing immunotherapy in glioma patients.

A prognostic molecular signature of hepatic steatosis is spatially heterogeneous and dynamic in human liver.

Metabolic dysfunction-associated steatotic liver disease (MASLD) prevalence is increasing in parallel with an obesity pandemic, calling for novel strategies for prevention and treatment. We defined a circulating proteome of human MASLD across ≈7000 proteins in ≈5000 individuals from diverse, at-risk populations across the metabolic health spectrum, demonstrating reproducible diagnostic performance and specifying both known and novel metabolic pathways relevant to MASLD (central carbon and amino acid metabolism, hepatocyte regeneration, inflammation, fibrosis, insulin sensitivity). A parsimonious proteomic signature of MASLD was associated with a protection from MASLD and its related multi-system metabolic consequences in >26000 free-living individuals, with an additive effect to polygenic risk. The MASLD proteome was encoded by genes that demonstrated transcriptional enrichment in liver, with spatial transcriptional activity in areas of steatosis in human liver biopsy and dynamicity for select targets in human liver across stages of steatosis. We replicated several top relations from proteomics and spatial tissue transcriptomics in a humanized "liver-on-a-chip" model of MASLD, highlighting the power of a full translational approach to discovery in MASLD. Collectively, these results underscore utility of blood-based proteomics as a dynamic "liquid biopsy" of human liver relevant to clinical biomarker and mechanistic applications.

Proteomic Biomarkers of Quantitative Interstitial Abnormalities in COPDGene and CARDIA Lung Study.

Rationale: Quantitative interstitial abnormalities (QIAs) are early measures of lung injury automatically detected on chest computed tomography scans. QIAs are associated with impaired respiratory health and share features with advanced lung diseases, but their biological underpinnings are not well understood. Objectives: To identify novel protein biomarkers of QIAs using high-throughput plasma proteomic panels within two multicenter cohorts. Methods: We measured the plasma proteomics of 4,383 participants in an older, ever-smoker cohort (COPDGene [Genetic Epidemiology of Chronic Obstructive Pulmonary Disease]) and 2,925 participants in a younger population cohort (CARDIA [Coronary Artery Disease Risk in Young Adults]) using the SomaLogic SomaScan assays. We measured QIAs using a local density histogram method. We assessed the associations between proteomic biomarker concentrations and QIAs using multivariable linear regression models adjusted for age, sex, body mass index, smoking status, and study center (Benjamini-Hochberg false discovery rate-corrected P ⩽ 0.05). Measurements and Main Results: In total, 852 proteins were significantly associated with QIAs in COPDGene and 185 in CARDIA. Of the 144 proteins that overlapped between COPDGene and CARDIA, all but one shared directionalities and magnitudes. These proteins were enriched for 49 Gene Ontology pathways, including biological processes in inflammatory response, cell adhesion, immune response, ERK1/2 regulation, and signaling; cellular components in extracellular regions; and molecular functions including calcium ion and heparin binding. Conclusions: We identified the proteomic biomarkers of QIAs in an older, smoking population with a higher prevalence of pulmonary disease and in a younger, healthier community cohort. These proteomics features may be markers of early precursors of advanced lung diseases.

Spatial Transcriptomics Resolve an Emphysema-Specific Lymphoid Follicle B Cell Signature in Chronic Obstructive Pulmonary Disease.

Rationale: Within chronic obstructive pulmonary disease (COPD), emphysema is characterized by a significant yet partially understood B cell immune component. Objectives: To characterize the transcriptomic signatures from lymphoid follicles (LFs) in ever-smokers without COPD and patients with COPD with varying degrees of emphysema. Methods: Lung sections from 40 patients with COPD and ever-smokers were used for LF proteomic and transcriptomic spatial profiling. Formalin- and O.C.T.-fixed lung samples obtained from biopsies or lung explants were assessed for LF presence. Emphysema measurements were obtained from clinical chest computed tomographic scans. High-confidence transcriptional target intersection analyses were conducted to resolve emphysema-induced transcriptional networks. Measurements and Main Results: Overall, 115 LFs from ever-smokers and Global Initiative for Chronic Obstructive Lung Disease (GOLD) 1-2 and GOLD 3-4 patients were analyzed. No LFs were found in never-smokers. Differential gene expression analysis revealed significantly increased expression of LF assembly and B cell marker genes in subjects with severe emphysema. High-confidence transcriptional analysis revealed activation of an abnormal B cell activity signature in LFs (q-value = 2.56E-111). LFs from patients with GOLD 1-2 COPD with emphysema showed significantly increased expression of genes associated with antigen presentation, inflammation, and B cell activation and proliferation. LFs from patients with GOLD 1-2 COPD without emphysema showed an antiinflammatory profile. The extent of centrilobular emphysema was significantly associated with genes involved in B cell maturation and antibody production. Protein-RNA network analysis showed that LFs in emphysema have a unique signature skewed toward chronic B cell activation. Conclusions: An off-targeted B cell activation within LFs is associated with autoimmune-mediated emphysema pathogenesis.

Spatial Transcriptomics Resolve an Emphysema-specific Lymphoid Follicle B Cell Signature in COPD.

RATIONALE: Within chronic obstructive pulmonary disease (COPD), emphysema is characterized by a significant yet partially understood B cell immune component. OBJECTIVE: To characterize the transcriptomic signatures from lymphoid follicles (LFs) in ever-smokers without COPD and COPD patients with varying degrees of emphysema. METHODS: Lung sections from 40 COPD patients and ever-smokers were used for LF proteomic and transcriptomic spatial profiling. Formalin and OCT-fixed lung samples obtained from biopsies or lung explants, were assessed for LF presence. Emphysema measurements were obtained from clinical chest CT scans. High confidence transcriptional (HCT) target intersection analyses were conducted to resolve emphysema-induced transcriptional networks. MEASUREMENTS AND MAIN RESULTS: Overall, 115 LFs from ever-smokers and GOLD 1-2 and GOLD 3-4 patients were analyzed. No LFs were found in never-smokers. Differential gene expression analysis revealed significantly increased expression of LF assembly and B cell markers genes in subjects with severe emphysema. HCT analysis revealed activation of abnormal B cell activity signature in LFs (q-value: 2.56E-111). LFs from GOLD 1-2 COPD patients with emphysema showed significantly increased expression of genes associated with antigen presentation, inflammation, and B cell activation and proliferation. LFs from GOLD 1-2 COPD patients without emphysema showed an anti-inflammatory profile. The extent of centrilobular emphysema was significantly associated with genes involved in B cell maturation and antibody production. Protein-RNA network analysis showed that LFs in emphysema have a unique signature skewed towards chronic B cell activation. CONCLUSIONS: An off-targeted B cell activation within LFs is associated with autoimmune-mediated emphysema pathogenesis.

Cell type-specific and disease-associated eQTL in the human lung.

Common genetic variants confer substantial risk for chronic lung diseases, including pulmonary fibrosis (PF). Defining the genetic control of gene expression in a cell-type-specific and context-dependent manner is critical for understanding the mechanisms through which genetic variation influences complex traits and disease pathobiology. To this end, we performed single-cell RNA-sequencing of lung tissue from 67 PF and 49 unaffected donors. Employing a pseudo-bulk approach, we mapped expression quantitative trait loci (eQTL) across 38 cell types, observing both shared and cell type-specific regulatory effects. Further, we identified disease-interaction eQTL and demonstrated that this class of associations is more likely to be cell-type specific and linked to cellular dysregulation in PF. Finally, we connected PF risk variants to their regulatory targets in disease-relevant cell types. These results indicate that cellular context determines the impact of genetic variation on gene expression, and implicates context-specific eQTL as key regulators of lung homeostasis and disease.

PCSK6 and Survival in Idiopathic Pulmonary Fibrosis.

Rationale: Idiopathic pulmonary fibrosis (IPF) is a devastating disease characterized by limited treatment options and high mortality. A better understanding of the molecular drivers of IPF progression is needed. Objectives: To identify and validate molecular determinants of IPF survival. Methods: A staged genome-wide association study was performed using paired genomic and survival data. Stage I cases were drawn from centers across the United States and Europe and stage II cases from Vanderbilt University. Cox proportional hazards regression was used to identify gene variants associated with differential transplantation-free survival (TFS). Stage I variants with nominal significance (P < 5 × 10-5) were advanced for stage II testing and meta-analyzed to identify those reaching genome-wide significance (P < 5 × 10-8). Downstream analyses were performed for genes and proteins associated with variants reaching genome-wide significance. Measurements and Main Results: After quality controls, 1,481 stage I cases and 397 stage II cases were included in the analysis. After filtering, 9,075,629 variants were tested in stage I, with 158 meeting advancement criteria. Four variants associated with TFS with consistent effect direction were identified in stage II, including one in an intron of PCSK6 (proprotein convertase subtilisin/kexin type 6) reaching genome-wide significance (hazard ratio, 4.11 [95% confidence interval, 2.54-6.67]; P = 9.45 × 10-9). PCSK6 protein was highly expressed in IPF lung parenchyma. PCSK6 lung staining intensity, peripheral blood gene expression, and plasma concentration were associated with reduced TFS. Conclusions: We identified four novel variants associated with IPF survival, including one in PCSK6 that reached genome-wide significance. Downstream analyses suggested that PCSK6 protein plays a potentially important role in IPF progression.

The Genetic Landscape of Familial Pulmonary Fibrosis.

Rationale and Objectives: Up to 20% of idiopathic interstitial lung disease is familial, referred to as familial pulmonary fibrosis (FPF). An integrated analysis of FPF genetic risk was performed by comprehensively evaluating for genetic rare variants (RVs) in a large cohort of FPF kindreds. Methods: Whole-exome sequencing and/or candidate gene sequencing from affected individuals in 569 FPF kindreds was performed, followed by cosegregation analysis in large kindreds, gene burden analysis, gene-based risk scoring, cell-type enrichment analysis, and coexpression network construction. Measurements and Main Results: It was found that 14.9-23.4% of genetic risk in kindreds could be explained by RVs in genes previously linked to FPF, predominantly telomere-related genes. New candidate genes were identified in a small number of families-including SYDE1, SERPINB8, GPR87, and NETO1-and tools were developed for evaluation and prioritization of RV-containing genes across kindreds. Several pathways were enriched for RV-containing genes in FPF, including focal adhesion and mitochondrial complex I assembly. By combining single-cell transcriptomics with prioritized candidate genes, expression of RV-containing genes was discovered to be enriched in smooth muscle cells, type II alveolar epithelial cells, and endothelial cells. Conclusions: In the most comprehensive FPF genetic study to date, the prevalence of RVs in known FPF-related genes was defined, and new candidate genes and pathways relevant to FPF were identified. However, new RV-containing genes shared across multiple kindreds were not identified, thereby suggesting that heterogeneous genetic variants involving a variety of genes and pathways mediate genetic risk in most FPF kindreds.

Image-based spatial transcriptomics identifies molecular niche dysregulation associated with distal lung remodeling in pulmonary fibrosis.

The human lung is structurally complex, with a diversity of specialized epithelial, stromal and immune cells playing specific functional roles in anatomically distinct locations, and large-scale changes in the structure and cellular makeup of this distal lung is a hallmark of pulmonary fibrosis (PF) and other progressive chronic lung diseases. Single-cell transcriptomic studies have revealed numerous disease-emergent/enriched cell types/states in PF lungs, but the spatial contexts wherein these cells contribute to disease pathogenesis has remained uncertain. Using sub-cellular resolution image-based spatial transcriptomics, we analyzed the gene expression of more than 1 million cells from 19 unique lungs. Through complementary cell-based and innovative cell-agnostic analyses, we characterized the localization of PF-emergent cell-types, established the cellular and molecular basis of classical PF histopathologic disease features, and identified a diversity of distinct molecularly-defined spatial niches in control and PF lungs. Using machine-learning and trajectory analysis methods to segment and rank airspaces on a gradient from normal to most severely remodeled, we identified a sequence of compositional and molecular changes that associate with progressive distal lung pathology, beginning with alveolar epithelial dysregulation and culminating with changes in macrophage polarization. Together, these results provide a unique, spatially-resolved characterization of the cellular and molecular programs of PF and control lungs, provide new insights into the heterogeneous pathobiology of PF, and establish analytical approaches which should be broadly applicable to other imaging-based spatial transcriptomic studies.

Thromboxane-Prostanoid Receptor Signaling Drives Persistent Fibroblast Activation in Pulmonary Fibrosis.

Rationale: Although persistent fibroblast activation is a hallmark of idiopathic pulmonary fibrosis (IPF), mechanisms regulating persistent fibroblast activation in the lungs have not been fully elucidated. Objectives: On the basis of our observation that lung fibroblasts express TBXA2R (thromboxane-prostanoid receptor) during fibrosis, we investigated the role of TBXA2R signaling in fibrotic remodeling. Methods: We identified TBXA2R expression in lungs of patients with IPF and mice and studied primary mouse and human lung fibroblasts to determine the impact of TBXA2R signaling on fibroblast activation. We used TBXA2R-deficient mice and small-molecule inhibitors to investigate TBXA2R signaling in preclinical lung fibrosis models. Measurements and Main Results: TBXA2R expression was upregulated in fibroblasts in the lungs of patients with IPF and in mouse lungs during experimental lung fibrosis. Genetic deletion of TBXA2R, but not inhibition of thromboxane synthase, protected mice from bleomycin-induced lung fibrosis, thereby suggesting that an alternative ligand activates profibrotic TBXA2R signaling. In contrast to thromboxane, F2-isoprostanes, which are nonenzymatic products of arachidonic acid induced by reactive oxygen species, were persistently elevated during fibrosis. F2-isoprostanes induced TBXA2R signaling in fibroblasts and mediated a myofibroblast activation profile due, at least in part, to potentiation of TGF-ß (transforming growth factor-ß) signaling. In vivo treatment with the TBXA2R antagonist ifetroban reduced profibrotic signaling in the lungs, protected mice from lung fibrosis in three preclinical models (bleomycin, Hermansky-Pudlak mice, and radiation-induced fibrosis), and markedly enhanced fibrotic resolution after bleomycin treatment. Conclusions: TBXA2R links oxidative stress to fibroblast activation during lung fibrosis. TBXA2R antagonists could have utility in treating pulmonary fibrosis.

Secretory Cells Are the Primary Source of pIgR in Small Airways.

Loss of secretory IgA (SIgA) is common in chronic obstructive pulmonary disease (COPD) small airways and likely contributes to disease progression. We hypothesized that loss of SIgA results from reduced expression of pIgR (polymeric immunoglobulin receptor), a chaperone protein needed for SIgA transcytosis, in the COPD small airway epithelium. pIgR-expressing cells were defined and quantified at single-cell resolution in human airways using RNA in situ hybridization, immunostaining, and single-cell RNA sequencing. Complementary studies in mice used immunostaining, primary murine tracheal epithelial cell culture, and transgenic mice with secretory or ciliated cell-specific knockout of pIgR. SIgA degradation by human neutrophil elastase or secreted bacterial proteases from nontypeable Haemophilus influenzae was evaluated in vitro. We found that secretory cells are the predominant cell type responsible for pIgR expression in human and murine airways. Loss of SIgA in small airways was not associated with a reduction in secretory cells but rather a reduction in pIgR protein expression despite intact PIGR mRNA expression. Neutrophil elastase and nontypeable H. influenzae-secreted proteases are both capable of degrading SIgA in vitro and may also contribute to a deficient SIgA immunobarrier in COPD. Loss of the SIgA immunobarrier in small airways of patients with severe COPD is complex and likely results from both pIgR-dependent defects in IgA transcytosis and SIgA degradation.

Multiplatform Single-Cell Analysis Identifies Immune Cell Types Enhanced in Pulmonary Fibrosis.

Immune cells have been implicated in idiopathic pulmonary fibrosis (IPF), but the phenotypes and effector mechanisms of these cells remain incompletely characterized. We performed mass cytometry to quantify immune cell subsets in lungs of 12 patients with IPF and 15 organ donors without chronic lung disease and used existing single-cell RNA-sequencing data to investigate transcriptional profiles of immune cells overrepresented in IPF. Among myeloid cells, we found increased numbers of alveolar macrophages (AMØs) and dendritic cells (DCs) in IPF, as well as a subset of monocyte-derived DCs. In contrast, monocyte-like cells and interstitial macrophages were reduced in IPF. Transcriptomic profiling identified an enrichment for IFN-γ response pathways in AMØs and DCs from IPF, as well as antigen processing in DCs and phagocytosis in AMØs. Among T cells, we identified three subsets of memory T cells that were increased in IPF, including CD4+ and CD8+ resident memory T cells (TRM) and CD8+ effector memory cells. The response to the IFN-γ pathway was enriched in CD4 TRM and CD8 TRM cells in IPF, together with T cell activation and immune response-regulating signaling pathways. Increased AMØs, DCs, and memory T cells were present in IPF lungs compared with control subjects. In IPF, these cells possess an activation profile indicating increased IFN-γ signaling and upregulation of adaptive immunity in the lungs. Together, these studies highlight critical features of the immunopathogenesis of IPF.

Genetic architecture of gene regulation in Indonesian populations identifies QTLs associated with global and local ancestries.

Lack of diversity in human genomics limits our understanding of the genetic underpinnings of complex traits, hinders precision medicine, and contributes to health disparities. To map genetic effects on gene regulation in the underrepresented Indonesian population, we have integrated genotype, gene expression, and CpG methylation data from 115 participants across three island populations that capture the major sources of genomic diversity in the region. In a comparison with European datasets, we identify eQTLs shared between Indonesia and Europe as well as population-specific eQTLs that exhibit differences in allele frequencies and/or overall expression levels between populations. By combining local ancestry and archaic introgression inference with eQTLs and methylQTLs, we identify regulatory loci driven by modern Papuan ancestry as well as introgressed Denisovan and Neanderthal variation. GWAS colocalization connects QTLs detected here to hematological traits, and further comparison with European datasets reflects the poor overall transferability of GWAS statistics across diverse populations. Our findings illustrate how population-specific genetic architecture, local ancestry, and archaic introgression drive variation in gene regulation across genetically distinct and in admixed populations and highlight the need for performing association studies on non-European populations.