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.

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.

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.

A single-cell atlas of mouse lung development.

Lung organogenesis requires precise timing and coordination to effect spatial organization and function of the parenchymal cells. To provide a systematic broad-based view of the mechanisms governing the dynamic alterations in parenchymal cells over crucial periods of development, we performed a single-cell RNA-sequencing time-series yielding 102,571 epithelial, endothelial and mesenchymal cells across nine time points from embryonic day 12 to postnatal day 14 in mice. Combining computational fate-likelihood prediction with RNA in situ hybridization and immunofluorescence, we explore lineage relationships during the saccular to alveolar stage transition. The utility of this publicly searchable atlas resource (www.sucrelab.org/lungcells) is exemplified by discoveries of the complexity of type 1 pneumocyte function and characterization of mesenchymal Wnt expression patterns during the saccular and alveolar stages - wherein major expansion of the gas-exchange surface occurs. We provide an integrated view of cellular dynamics in epithelial, endothelial and mesenchymal cell populations during lung organogenesis.

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.

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.

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.

A comparison of gene expression and DNA methylation patterns across tissues and species.

Previously published comparative functional genomic data sets from primates using frozen tissue samples, including many data sets from our own group, were often collected and analyzed using nonoptimal study designs and analysis approaches. In addition, when samples from multiple tissues were studied in a comparative framework, individuals and tissues were confounded. We designed a multitissue comparative study of gene expression and DNA methylation in primates that minimizes confounding effects by using a balanced design with respect to species, tissues, and individuals. We also developed a comparative analysis pipeline that minimizes biases attributable to sequence divergence. Thus, we present the most comprehensive catalog of similarities and differences in gene expression and DNA methylation levels between livers, kidneys, hearts, and lungs, in humans, chimpanzees, and rhesus macaques. We estimate that overall, interspecies and inter-tissue differences in gene expression levels can only modestly be accounted for by corresponding differences in promoter DNA methylation. However, the expression pattern of genes with conserved inter-tissue expression differences can be explained by corresponding interspecies methylation changes more often. Finally, we show that genes whose tissue-specific regulatory patterns are consistent with the action of natural selection are highly connected in both gene regulatory and protein-protein interaction networks.

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.

Genome-wide DNA methylation and gene expression patterns reflect genetic ancestry and environmental differences across the Indonesian archipelago.

Indonesia is the world's fourth most populous country, host to striking levels of human diversity, regional patterns of admixture, and varying degrees of introgression from both Neanderthals and Denisovans. However, it has been largely excluded from the human genomics sequencing boom of the last decade. To serve as a benchmark dataset of molecular phenotypes across the region, we generated genome-wide CpG methylation and gene expression measurements in over 100 individuals from three locations that capture the major genomic and geographical axes of diversity across the Indonesian archipelago. Investigating between- and within-island differences, we find up to 10.55% of tested genes are differentially expressed between the islands of Sumba and New Guinea. Variation in gene expression is closely associated with DNA methylation, with expression levels of 9.80% of genes correlating with nearby promoter CpG methylation, and many of these genes being differentially expressed between islands. Genes identified in our differential expression and methylation analyses are enriched in pathways involved in immunity, highlighting Indonesia's tropical role as a source of infectious disease diversity and the strong selective pressures these diseases have exerted on humans. Finally, we identify robust within-island variation in DNA methylation and gene expression, likely driven by fine-scale environmental differences across sampling sites. Together, these results strongly suggest complex relationships between DNA methylation, transcription, archaic hominin introgression and immunity, all jointly shaped by the environment. This has implications for the application of genomic medicine, both in critically understudied Indonesia and globally, and will allow a better understanding of the interacting roles of genomic and environmental factors shaping molecular and complex phenotypes.

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.

Integrated Single-Cell Atlas of Endothelial Cells of the Human Lung.

BACKGROUND: Cellular diversity of the lung endothelium has not been systematically characterized in humans. We provide a reference atlas of human lung endothelial cells (ECs) to facilitate a better understanding of the phenotypic diversity and composition of cells comprising the lung endothelium. METHODS: We reprocessed human control single-cell RNA sequencing (scRNAseq) data from 6 datasets. EC populations were characterized through iterative clustering with subsequent differential expression analysis. Marker genes were validated by fluorescent microscopy and in situ hybridization. scRNAseq of primary lung ECs cultured in vitro was performed. The signaling network between different lung cell types was studied. For cross-species analysis or disease relevance, we applied the same methods to scRNAseq data obtained from mouse lungs or from human lungs with pulmonary hypertension. RESULTS: Six lung scRNAseq datasets were reanalyzed and annotated to identify >15 000 vascular EC cells from 73 individuals. Differential expression analysis of EC revealed signatures corresponding to endothelial lineage, including panendothelial, panvascular, and subpopulation-specific marker gene sets. Beyond the broad cellular categories of lymphatic, capillary, arterial, and venous ECs, we found previously indistinguishable subpopulations; among venous EC, we identified 2 previously indistinguishable populations: pulmonary-venous ECs (COL15A1neg) localized to the lung parenchyma and systemic-venous ECs (COL15A1pos) localized to the airways and the visceral pleura; among capillary ECs, we confirmed their subclassification into recently discovered aerocytes characterized by EDNRB, SOSTDC1, and TBX2 and general capillary EC. We confirmed that all 6 endothelial cell types, including the systemic-venous ECs and aerocytes, are present in mice and identified endothelial marker genes conserved in humans and mice. Ligand-receptor connectome analysis revealed important homeostatic crosstalk of EC with other lung resident cell types. scRNAseq of commercially available primary lung ECs demonstrated a loss of their native lung phenotype in culture. scRNAseq revealed that endothelial diversity is maintained in pulmonary hypertension. Our article is accompanied by an online data mining tool (www.LungEndothelialCellAtlas.com). CONCLUSIONS: Our integrated analysis provides a comprehensive and well-crafted reference atlas of ECs in the normal lung and confirms and describes in detail previously unrecognized endothelial populations across a large number of humans and mice.

Impact of regulatory variation across human iPSCs and differentiated cells.

Induced pluripotent stem cells (iPSCs) are an essential tool for studying cellular differentiation and cell types that are otherwise difficult to access. We investigated the use of iPSCs and iPSC-derived cells to study the impact of genetic variation on gene regulation across different cell types and as models for studies of complex disease. To do so, we established a panel of iPSCs from 58 well-studied Yoruba lymphoblastoid cell lines (LCLs); 14 of these lines were further differentiated into cardiomyocytes. We characterized regulatory variation across individuals and cell types by measuring gene expression levels, chromatin accessibility, and DNA methylation. Our analysis focused on a comparison of inter-individual regulatory variation across cell types. While most cell-type-specific regulatory quantitative trait loci (QTLs) lie in chromatin that is open only in the affected cell types, we found that 20% of cell-type-specific regulatory QTLs are in shared open chromatin. This observation motivated us to develop a deep neural network to predict open chromatin regions from DNA sequence alone. Using this approach, we were able to use the sequences of segregating haplotypes to predict the effects of common SNPs on cell-type-specific chromatin accessibility.

Genetic Variation, Not Cell Type of Origin, Underlies the Majority of Identifiable Regulatory Differences in iPSCs.

The advent of induced pluripotent stem cells (iPSCs) revolutionized human genetics by allowing us to generate pluripotent cells from easily accessible somatic tissues. This technology can have immense implications for regenerative medicine, but iPSCs also represent a paradigm shift in the study of complex human phenotypes, including gene regulation and disease. Yet, an unresolved caveat of the iPSC model system is the extent to which reprogrammed iPSCs retain residual phenotypes from their precursor somatic cells. To directly address this issue, we used an effective study design to compare regulatory phenotypes between iPSCs derived from two types of commonly used somatic precursor cells. We find a remarkably small number of differences in DNA methylation and gene expression levels between iPSCs derived from different somatic precursors. Instead, we demonstrate genetic variation is associated with the majority of identifiable variation in DNA methylation and gene expression levels. We show that the cell type of origin only minimally affects gene expression levels and DNA methylation in iPSCs, and that genetic variation is the main driver of regulatory differences between iPSCs of different donors. Our findings suggest that studies using iPSCs should focus on additional individuals rather than clones from the same individual.

Methylation QTLs are associated with coordinated changes in transcription factor binding, histone modifications, and gene expression levels.

DNA methylation is an important epigenetic regulator of gene expression. Recent studies have revealed widespread associations between genetic variation and methylation levels. However, the mechanistic links between genetic variation and methylation remain unclear. To begin addressing this gap, we collected methylation data at ∼300,000 loci in lymphoblastoid cell lines (LCLs) from 64 HapMap Yoruba individuals, and genome-wide bisulfite sequence data in ten of these individuals. We identified (at an FDR of 10%) 13,915 cis methylation QTLs (meQTLs)-i.e., CpG sites in which changes in DNA methylation are associated with genetic variation at proximal loci. We found that meQTLs are frequently associated with changes in methylation at multiple CpGs across regions of up to 3 kb. Interestingly, meQTLs are also frequently associated with variation in other properties of gene regulation, including histone modifications, DNase I accessibility, chromatin accessibility, and expression levels of nearby genes. These observations suggest that genetic variants may lead to coordinated molecular changes in all of these regulatory phenotypes. One plausible driver of coordinated changes in different regulatory mechanisms is variation in transcription factor (TF) binding. Indeed, we found that SNPs that change predicted TF binding affinities are significantly enriched for associations with DNA methylation at nearby CpGs.

Association between mosaic loss of chromosome Y and pulmonary fibrosis susceptibility and severity.

Background Pulmonary fibrosis (PF) is a rare lung disease with diverse pathogenesis and multiple interconnected underlying biological mechanisms. Mosaic loss of chromosome Y (mLOY) is one of the most common forms of acquired chromosome abnormality in men, which has been reported to be associated with increased risk of various chronic progressive diseases including fibrotic diseases. However, the exact role of mLOY in the development of PF remains elusive and to be elucidated. Methods: We adopted three complementary approaches to explore the role of mLOY in the pathogenesis of PF. We used copy number on chromosome Y to estimate mLOY comparing patients in PROFILE and gnomAD cohorts and between cases and control patients from the GE100KGP cohort. Correlation of mLOY with demographic and clinical variables was tested using patients from PROFILE cohort. Lung single-cell transcriptomic data were analysed to assess the cell types implicated in mLOY. We performed Mendelian randomisation to examine the causal relationship between mLOY, IPF, and telomere length. Results: The genetic analysis suggests that mLOY is found in PF from both case cohorts but when compared with an age matched population the effect is minimal (P = 0.0032). mLOY is related to age (P = 0.00021) and shorter telomere length (P = 0.0081) rather than PF severity or progression. Single-cell analysis indicates that mLOY appears to be found primarily in immune cells and appears to be related to presence and severity of fibrosis. Mendelian randomisation demonstrates that mLOY is not on the causal pathway for IPF, but partial evidence supports that telomere shortening is on the causal pathway for mLOY. Conclusion: Our study confirms the existence of mLOY in PF patients and suggests that mLOY is not a major driver of IPF. The combined evidence suggests a triangulation model where telomere shortening leads to both IPF and mLOY.

Bronchopulmonary Dysplasia with Pulmonary Hypertension Associates with Loss of Semaphorin Signaling and Functional Decrease in FOXF1 Expression.

Lung injury in preterm infants leads to structural and functional respiratory deficits, with a risk for bronchopulmonary dysplasia (BPD) that in its most severe form is accompanied by pulmonary hypertension (PH). To examine cellular and molecular dynamics driving evolving BPD in humans, we performed single-cell RNA sequencing of preterm infant lungs in early stages of BPD and BPD+PH compared to term infants. Analysis of the endothelium revealed a unique aberrant capillary cell-state primarily in BPD+PH marked by ANKRD1 expression. Predictive signaling analysis identified deficits in the semaphorin guidance-cue signaling pathway and decreased expression of pro-angiogenic transcription factor FOXF1 within the alveolar parenchyma in neonatal lung samples with BPD/BPD+PH. Loss of semaphorin signaling was replicated in a murine BPD model and in humans with alveolar capillary dysplasia (ACDMPV), suggesting a mechanistic link between the developmental programs underlying BPD and ACDMPV and a critical role for semaphorin signaling in normal lung development.

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.

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.