Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program.

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV1] and its ratio to forced vital capacity [FEV1/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV1 and FEV1/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Metabolomic Signatures Associated With Pulmonary Arterial Hypertension Outcomes.

BACKGROUND: Pulmonary arterial hypertension (PAH) is a complex disease characterized by progressive right ventricular (RV) failure leading to significant morbidity and mortality. Investigating metabolic features and pathways associated with RV dilation, mortality, and measures of disease severity can provide insight into molecular mechanisms, identify subphenotypes, and suggest potential therapeutic targets. METHODS: We collected data from a prospective cohort of PAH participants and performed untargeted metabolomic profiling on 1045 metabolites from circulating blood. Analyses were intended to identify metabolomic differences across a range of common metrics in PAH (eg, dilated versus nondilated RV). Partial least squares discriminant analysis was first applied to assess the distinguishability of relevant outcomes. Significantly altered metabolites were then identified using linear regression, and Cox regression models (as appropriate for the specific outcome) with adjustments for age, sex, body mass index, and PAH cause. Models exploring RV maladaptation were further adjusted for pulmonary vascular resistance. Pathway enrichment analysis was performed to identify significantly dysregulated processes. RESULTS: A total of 117 participants with PAH were included. Partial least squares discriminant analysis showed cluster differentiation between participants with dilated versus nondilated RVs, survivors versus nonsurvivors, and across a range of NT-proBNP (N-terminal pro-B-type natriuretic peptide) levels, REVEAL 2.0 composite scores, and 6-minute-walk distances. Polyamine and histidine pathways were associated with differences in RV dilation, mortality, NT-proBNP, REVEAL score, and 6-minute walk distance. Acylcarnitine pathways were associated with NT-proBNP, REVEAL score, and 6-minute walk distance. Sphingomyelin pathways were associated with RV dilation and NT-proBNP after adjustment for pulmonary vascular resistance. CONCLUSIONS: Distinct plasma metabolomic profiles are associated with RV dilation, mortality, and measures of disease severity in PAH. Polyamine, histidine, and sphingomyelin metabolic pathways represent promising candidates for identifying patients at high risk for poor outcomes and investigation into their roles as markers or mediators of disease progression and RV adaptation.

Plasma Metabolomics Reveals Distinct Biological and Diagnostic Signatures for Melioidosis.

Rationale: The global burden of sepsis is greatest in low-resource settings. Melioidosis, infection with the gram-negative bacterium Burkholderia pseudomallei, is a frequent cause of fatal sepsis in endemic tropical regions such as Southeast Asia. Objectives: To investigate whether plasma metabolomics would identify biological pathways specific to melioidosis and yield clinically meaningful biomarkers. Methods: Using a comprehensive approach, differential enrichment of plasma metabolites and pathways was systematically evaluated in individuals selected from a prospective cohort of patients hospitalized in rural Thailand with infection. Statistical and bioinformatics methods were used to distinguish metabolomic features and processes specific to patients with melioidosis and between fatal and nonfatal cases. Measurements and Main Results: Metabolomic profiling and pathway enrichment analysis of plasma samples from patients with melioidosis (n = 175) and nonmelioidosis infections (n = 75) revealed a distinct immuno-metabolic state among patients with melioidosis, as suggested by excessive tryptophan catabolism in the kynurenine pathway and significantly increased levels of sphingomyelins and ceramide species. We derived a 12-metabolite classifier to distinguish melioidosis from other infections, yielding an area under the receiver operating characteristic curve of 0.87 in a second validation set of patients. Melioidosis nonsurvivors (n = 94) had a significantly disturbed metabolome compared with survivors (n = 81), with increased leucine, isoleucine, and valine metabolism, and elevated circulating free fatty acids and acylcarnitines. A limited eight-metabolite panel showed promise as an early prognosticator of mortality in melioidosis. Conclusions: Melioidosis induces a distinct metabolomic state that can be examined to distinguish underlying pathophysiological mechanisms associated with death. A 12-metabolite signature accurately differentiates melioidosis from other infections and may have diagnostic applications.

Leveraging pleiotropy to discover and interpret GWAS results for sleep-associated traits.

Genetic association studies of many heritable traits resulting from physiological testing often have modest sample sizes due to the cost and burden of the required phenotyping. This reduces statistical power and limits discovery of multiple genetic associations. We present a strategy to leverage pleiotropy between traits to both discover new loci and to provide mechanistic hypotheses of the underlying pathophysiology. Specifically, we combine a colocalization test with a locus-level test of pleiotropy. In simulations, we show that this approach is highly selective for identifying true pleiotropy driven by the same causative variant, thereby improves the chance to replicate the associations in underpowered validation cohorts and leads to higher interpretability. Here, as an exemplar, we use Obstructive Sleep Apnea (OSA), a common disorder diagnosed using overnight multi-channel physiological testing. We leverage pleiotropy with relevant cellular and cardio-metabolic phenotypes and gene expression traits to map new risk loci in an underpowered OSA GWAS. We identify several pleiotropic loci harboring suggestive associations to OSA and genome-wide significant associations to other traits, and show that their OSA association replicates in independent cohorts of diverse ancestries. By investigating pleiotropic loci, our strategy allows proposing new hypotheses about OSA pathobiology across many physiological layers. For example, we identify and replicate the pleiotropy across the plateletcrit, OSA and an eQTL of DNA primase subunit 1 (PRIM1) in immune cells. We find suggestive links between OSA, a measure of lung function (FEV1/FVC), and an eQTL of matrix metallopeptidase 15 (MMP15) in lung tissue. We also link a previously known genome-wide significant peak for OSA in the hexokinase 1 (HK1) locus to hematocrit and other red blood cell related traits. Thus, the analysis of pleiotropic associations has the potential to assemble diverse phenotypes into a chain of mechanistic hypotheses that provide insight into the pathogenesis of complex human diseases.

γδ T Cells Mediate Protection Against Neutrophil-associated Lung Inflammation in Pulmonary Melioidosis.

Pulmonary melioidosis is a severe tropical infection caused by Burkholderia pseudomallei and is associated with high mortality despite early antibiotic treatment. γδ T cells have been increasingly implicated as drivers of the host neutrophil response during bacterial pneumonia, but their role in pulmonary melioidosis is unknown. Here, we report that in patients with melioidosis, a lower peripheral blood γδ T cell concentration is associated with higher mortality even when adjusting for severity of illness. γδ T cells were also enriched in the lung and protected against mortality in a mouse model of pulmonary melioidosis. γδ T cell deficiency in infected mice induced an early recruitment of neutrophils to the lung, independent of bacterial burden. Subsequently, γδ T cell deficiency resulted in increased neutrophil-associated inflammation in the lung as well as impaired bacterial clearance. Additionally, γδ T cells influenced neutrophil function and subset diversity in the lung after infection. Our results indicate that γδ T cells serve a novel protective role in the lung during a severe bacterial pneumonia by regulating excessive neutrophil-associated inflammation.

Fatty acid metabolism promotes TRPV4 activity in lung microvascular endothelial cells in pulmonary arterial hypertension.

Pulmonary arterial hypertension (PAH) is a morbid disease characterized by significant lung endothelial cell (EC) dysfunction. Prior work has shown that microvascular endothelial cells (MVECs) isolated from animals with experimental PAH and patients with PAH exhibit significant abnormalities in metabolism and calcium signaling. With regards to metabolism, we and others have shown evidence of increased aerobic glycolysis and evidence of increased utilization of alternate fuel sources (such as fatty acids) in PAH EC. In the realm of calcium signaling, our prior work linked increased activity of the transient receptor potential vanilloid-4 (TRPV4) channel to increased proliferation of MVECs isolated from the Sugen/Hypoxia rat model of PAH (SuHx-MVECs). However, the relationship between metabolic shifts and calcium abnormalities was not clear. Specifically, whether shifts in metabolism were responsible for increasing TRPV4 channel activity in SuHx-MVECs was not known. In this study, using human data, serum samples from SuHx rats, and SuHx-MVECs, we describe the consequences of increased MVEC fatty acid oxidation in PAH. In human samples, we observed an increase in long-chain fatty acid levels that was associated with PAH severity. Next, using SuHx rats and SuHx-MVECs, we observed increased intracellular levels of lipids. We also show that increasing intracellular lipid content increases TRPV4 activity, whereas inhibiting fatty acid oxidation normalizes basal calcium levels in SuHx-MVECs. By exploring the fate of fatty acid-derived carbons, we observed that the metabolite linking increased intracellular lipids to TRPV4 activity was ß-hydroxybutyrate (BOHB), a product of fatty acid oxidation. Finally, we show that BOHB supplementation alone is sufficient to sensitize the TRPV4 channel in rat and mouse MVECs. Returning to humans, we observe a transpulmonary BOHB gradient in human patients with PAH. Thus, we establish a link between fatty acid oxidation, BOHB production, and TRPV4 activity in MVECs in PAH. These data provide new insight into metabolic regulation of calcium signaling in lung MVECs in PAH.NEW & NOTEWORTHY In this paper, we explore the link between metabolism and intracellular calcium levels in microvascular endothelial cells (MVECs) in pulmonary arterial hypertension (PAH). We show that fatty acid oxidation promotes sensitivity of the transient receptor potential vanilloid-4 (TRPV4) calcium channel in MVECs isolated from a rodent model of PAH.

Assessing the contribution of rare genetic variants to phenotypes of chronic obstructive pulmonary disease using whole-genome sequence data.

RATIONALE: Genetic variation has a substantial contribution to chronic obstructive pulmonary disease (COPD) and lung function measurements. Heritability estimates using genome-wide genotyping data can be biased if analyses do not appropriately account for the nonuniform distribution of genetic effects across the allele frequency and linkage disequilibrium (LD) spectrum. In addition, the contribution of rare variants has been unclear. OBJECTIVES: We sought to assess the heritability of COPD and lung function using whole-genome sequence data from the Trans-Omics for Precision Medicine program. METHODS: Using the genome-based restricted maximum likelihood method, we partitioned the genome into bins based on minor allele frequency and LD scores and estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio in 11 051 European ancestry and 5853 African-American participants. MEASUREMENTS AND MAIN RESULTS: In European ancestry participants, the estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio were 35.5%, 55.6% and 32.5%, of which 18.8%, 19.7%, 17.8% were from common variants, and 16.6%, 35.8%, and 14.6% were from rare variants. These estimates had wide confidence intervals, with common variants and some sets of rare variants showing a statistically significant contribution (P-value < 0.05). In African-Americans, common variant heritability was similar to European ancestry participants, but lower sample size precluded calculation of rare variant heritability. CONCLUSIONS: Our study provides updated and unbiased estimates of heritability for COPD and lung function, and suggests an important contribution of rare variants. Larger studies of more diverse ancestry will improve accuracy of these estimates.

Genome-wide association study of preserved ratio impaired spirometry (PRISm).

BACKGROUND: Preserved ratio impaired spirometry (PRISm) is defined as a forced expiratory volume in 1 s (FEV1) <80% predicted and FEV1/forced vital capacity ≥0.70. PRISm is associated with respiratory symptoms and comorbidities. Our objective was to discover novel genetic signals for PRISm and see if they provide insight into the pathogenesis of PRISm and associated comorbidities. METHODS: We undertook a genome-wide association study (GWAS) of PRISm in UK Biobank participants (Stage 1), and selected single nucleotide polymorphisms (SNPs) reaching genome-wide significance for replication in 13 cohorts (Stage 2). A combined meta-analysis of Stage 1 and Stage 2 was done to determine top SNPs. We used cross-trait linkage disequilibrium score regression to estimate genome-wide genetic correlation between PRISm and pulmonary and extrapulmonary traits. Phenome-wide association studies of top SNPs were performed. RESULTS: 22 signals reached significance in the joint meta-analysis, including four signals novel for lung function. A strong genome-wide genetic correlation (rg) between PRISm and spirometric COPD (rg=0.62, p<0.001) was observed, and genetic correlation with type 2 diabetes (rg=0.12, p=0.007). Phenome-wide association studies showed that 18 of 22 signals were associated with diabetic traits and seven with blood pressure traits. CONCLUSION: This is the first GWAS to successfully identify SNPs associated with PRISm. Four of the signals, rs7652391 (nearest gene MECOM), rs9431040 (HLX), rs62018863 (TMEM114) and rs185937162 (HLA-B), have not been described in association with lung function before, demonstrating the utility of using different lung function phenotypes in GWAS. Genetic factors associated with PRISm are strongly correlated with risk of both other lung diseases and extrapulmonary comorbidity.

Association of clonal hematopoiesis with chronic obstructive pulmonary disease.

Chronic obstructive pulmonary disease (COPD) is associated with age and smoking, but other determinants of the disease are incompletely understood. Clonal hematopoiesis of indeterminate potential (CHIP) is a common, age-related state in which somatic mutations in clonal blood populations induce aberrant inflammatory responses. Patients with CHIP have an elevated risk for cardiovascular disease, but the association of CHIP with COPD remains unclear. We analyzed whole-genome sequencing and whole-exome sequencing data to detect CHIP in 48 835 patients, of whom 8444 had moderate to very severe COPD, from four separate cohorts with COPD phenotyping and smoking history. We measured emphysema in murine models in which Tet2 was deleted in hematopoietic cells. In the COPDGene cohort, individuals with CHIP had risks of moderate-to-severe, severe, or very severe COPD that were 1.6 (adjusted 95% confidence interval [CI], 1.1-2.2) and 2.2 (adjusted 95% CI, 1.5-3.2) times greater than those for noncarriers. These findings were consistently observed in three additional cohorts and meta-analyses of all patients. CHIP was also associated with decreased FEV1% predicted in the COPDGene cohort (mean between-group differences, -5.7%; adjusted 95% CI, -8.8% to -2.6%), a finding replicated in additional cohorts. Smoke exposure was associated with a small but significant increased risk of having CHIP (odds ratio, 1.03 per 10 pack-years; 95% CI, 1.01-1.05 per 10 pack-years) in the meta-analysis of all patients. Inactivation of Tet2 in mouse hematopoietic cells exacerbated the development of emphysema and inflammation in models of cigarette smoke exposure. Somatic mutations in blood cells are associated with the development and severity of COPD, independent of age and cumulative smoke exposure.

Distinct Epithelial-Innate Immune Cell Transcriptional Circuits Underlie Airway Hyperresponsiveness in Asthma.

Rationale: Indirect airway hyperresponsiveness (AHR) is a highly specific feature of asthma, but the underlying mechanisms responsible for driving indirect AHR remain incompletely understood. Objectives: To identify differences in gene expression in epithelial brushings obtained from individuals with asthma who were characterized for indirect AHR in the form of exercise-induced bronchoconstriction (EIB). Methods: RNA-sequencing analysis was performed on epithelial brushings obtained from individuals with asthma with EIB (n = 11) and without EIB (n = 9). Differentially expressed genes (DEGs) between the groups were correlated with measures of airway physiology, sputum inflammatory markers, and airway wall immunopathology. On the basis of these relationships, we examined the effects of primary airway epithelial cells (AECs) and specific epithelial cell-derived cytokines on both mast cells (MCs) and eosinophils (EOS). Measurements and Main Results: We identified 120 DEGs in individuals with and without EIB. Network analyses suggested critical roles for IL-33-, IL-18-, and IFN-γ-related signaling among these DEGs. IL1RL1 expression was positively correlated with the density of MCs in the epithelial compartment, and IL1RL1, IL18R1, and IFNG were positively correlated with the density of intraepithelial EOS. Subsequent ex vivo modeling demonstrated that AECs promote sustained type 2 (T2) inflammation in MCs and enhance IL-33-induced T2 gene expression. Furthermore, EOS increase the expression of IFNG and IL13 in response to both IL-18 and IL-33 as well as exposure to AECs. Conclusions: Circuits involving epithelial interactions with MCs and EOS are closely associated with indirect AHR. Ex vivo modeling indicates that epithelial-dependent regulation of these innate cells may be critical in indirect AHR and modulating T2 and non-T2 inflammation in asthma.

Investigating Associations of Omega-3 Fatty Acids, Lung Function Decline, and Airway Obstruction.

Rationale: Inflammation contributes to lung function decline and the development of chronic obstructive pulmonary disease. Omega-3 fatty acids have antiinflammatory properties and may benefit lung health. Objectives: To investigate associations of omega-3 fatty acids with lung function decline and incident airway obstruction in a diverse sample of adults from general-population cohorts. Methods: Complementary study designs: 1) longitudinal study of plasma phospholipid omega-3 fatty acids and repeated FEV1 and FVC measures in the NHLBI Pooled Cohorts Study and 2) two-sample Mendelian randomization (MR) study of genetically predicted omega-3 fatty acids and lung function parameters. Measurements and Main Results: The longitudinal study found that higher omega-3 fatty acid levels were associated with attenuated lung function decline in 15,063 participants, with the largest effect sizes for the most metabolically downstream omega-3 fatty acid, docosahexaenoic acid (DHA). An increase in DHA of 1% of total fatty acids was associated with attenuations of 1.4 ml/yr for FEV1 (95% confidence interval [CI], 1.1-1.8) and 2.0 ml/yr for FVC (95% CI, 1.6-2.4) and a 7% lower incidence of spirometry-defined airway obstruction (95% CI, 0.89-0.97). DHA associations persisted across sexes and smoking histories and in Black, White, and Hispanic participants, with associations of the largest magnitude in former smokers and Hispanic participants. The MR study showed similar trends toward positive associations of genetically predicted downstream omega-3 fatty acids with FEV1 and FVC. Conclusions: The longitudinal and MR studies provide evidence supporting beneficial effects of higher levels of downstream omega-3 fatty acids, especially DHA, on lung health.

Rhinovirus infection of the airway epithelium enhances mast cell immune responses via epithelial-derived interferons.

BACKGROUND: Mast cells (MCs) within the airway epithelium in asthma are closely related to airway dysfunction, but cross talk between airway epithelial cells (AECs) and MCs in asthma remains incompletely understood. Human rhinovirus (RV) infections are key triggers for asthma progression, and AECs from individuals with asthma may have dysregulated antiviral responses. OBJECTIVE: We utilized primary AECs in an ex vivo coculture model system to examine cross talk between AECs and MCs after epithelial rhinovirus infection. METHODS: Primary AECs were obtained from 11 children with asthma and 10 healthy children, differentiated at air-liquid interface, and cultured in the presence of laboratory of allergic diseases 2 (LAD2) MCs. AECs were infected with rhinovirus serogroup A 16 (RV16) for 48 hours. RNA isolated from both AECs and MCs underwent RNA sequencing. Direct effects of epithelial-derived interferons on LAD2 MCs were examined by real-time quantitative PCR. RESULTS: MCs increased expression of proinflammatory and antiviral genes in AECs. AECs demonstrated a robust antiviral response after RV16 infection that resulted in significant changes in MC gene expression, including upregulation of genes involved in antiviral responses, leukocyte activation, and type 2 inflammation. Subsequent ex vivo modeling demonstrated that IFN-ß induces MC type 2 gene expression. The effects of AEC donor phenotype were small relative to the effects of viral infection and the presence of MCs. CONCLUSIONS: There is significant cross talk between AECs and MCs, which are present in the epithelium in asthma. Epithelial-derived interferons not only play a role in viral suppression but also further alter MC immune responses including specific type 2 genes.

Nephronectin is required to maintain right lung lobar separation during embryonic development.

Nephronectin (NPNT) is a basement membrane (BM) protein and high-affinity ligand of integrin α8ß1 that is required for kidney morphogenesis in mice. In the lung, NPNT also localizes to BMs, but its potential role in pulmonary development has not been investigated. Mice with a floxed Npnt allele were used to generate global knockouts (KOs). Staged embryos were obtained by timed matings of heterozygotes and lungs were isolated for analysis. Although primary and secondary lung bud formation was normal in KO embryos, fusion of right lung lobes, primarily the medial and caudal, was first detected at E13.5 and persisted into adulthood. The lung parenchyma of KO mice was indistinguishable from wild-type (WT) and lobe fusion did not alter respiratory mechanics in adult KO mice. Interrogation of an existing single-cell RNA-seq atlas of embryonic and adult mouse lungs identified Npnt transcripts in mesothelial cells at E12.5 and into the early postnatal period, but not in adult lungs. KO embryonic lungs exhibited increased expression of laminin α5 and deposition of collagen IV in the mesothelial BM, accompanied by abnormalities in collagen fibrils in the adjacent stroma. Cranial and accessory lobes extracted from KO embryonic lungs fused ex vivo when cultured in juxtaposition, with the area of fusion showing loss of the mesothelial marker Wilms tumor 1. Because a similar pattern of lobe fusion was previously observed in integrin α8 KO embryos, our results suggest that NPNT signaling through integrin α8, likely in the visceral pleura, maintains right lung lobe separation during embryogenesis.

Benchmarking association analyses of continuous exposures with RNA-seq in observational studies.

Large datasets of hundreds to thousands of individuals measuring RNA-seq in observational studies are becoming available. Many popular software packages for analysis of RNA-seq data were constructed to study differences in expression signatures in an experimental design with well-defined conditions (exposures). In contrast, observational studies may have varying levels of confounding transcript-exposure associations; further, exposure measures may vary from discrete (exposed, yes/no) to continuous (levels of exposure), with non-normal distributions of exposure. We compare popular software for gene expression-DESeq2, edgeR and limma-as well as linear regression-based analyses for studying the association of continuous exposures with RNA-seq. We developed a computation pipeline that includes transformation, filtering and generation of empirical null distribution of association P-values, and we apply the pipeline to compute empirical P-values with multiple testing correction. We employ a resampling approach that allows for assessment of false positive detection across methods, power comparison and the computation of quantile empirical P-values. The results suggest that linear regression methods are substantially faster with better control of false detections than other methods, even with the resampling method to compute empirical P-values. We provide the proposed pipeline with fast algorithms in an R package Olivia, and implemented it to study the associations of measures of sleep disordered breathing with RNA-seq in peripheral blood mononuclear cells in participants from the Multi-Ethnic Study of Atherosclerosis.

Targeted Genome Sequencing Identifies Multiple Rare Variants in Caveolin-1 Associated with Obstructive Sleep Apnea.

Rationale: Obstructive sleep apnea (OSA) is a common disorder associated with increased risk for cardiovascular disease, diabetes, and premature mortality. There is strong clinical and epidemiologic evidence supporting the importance of genetic factors influencing OSA but limited data implicating specific genes. Objectives: To search for rare variants contributing to OSA severity. Methods: Leveraging high-depth genomic sequencing data from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and imputed genotype data from multiple population-based studies, we performed linkage analysis in the CFS (Cleveland Family Study), followed by multistage gene-based association analyses in independent cohorts for apnea-hypopnea index (AHI) in a total of 7,708 individuals of European ancestry. Measurements and Main Results: Linkage analysis in the CFS identified a suggestive linkage peak on chromosome 7q31 (LOD = 2.31). Gene-based analysis identified 21 noncoding rare variants in CAV1 (Caveolin-1) associated with lower AHI after accounting for multiple comparisons (P = 7.4 × 10-8). These noncoding variants together significantly contributed to the linkage evidence (P < 10-3). Follow-up analysis revealed significant associations between these variants and increased CAV1 expression, and increased CAV1 expression in peripheral monocytes was associated with lower AHI (P = 0.024) and higher minimum overnight oxygen saturation (P = 0.007). Conclusions: Rare variants in CAV1, a membrane-scaffolding protein essential in multiple cellular and metabolic functions, are associated with higher CAV1 gene expression and lower OSA severity, suggesting a novel target for modulating OSA severity.

MAP2K2 Delays Recovery in Murine Models of Acute Lung Injury and Associates with Acute Respiratory Distress Syndrome Outcome.

Acute respiratory distress syndrome (ARDS) remains a significant problem in need of new pharmaceutical approaches to improve its resolution. Studies comparing gene expression signatures in rodents and humans with lung injury reveal conserved pathways, including MAPK (mitogen-activated protein kinase)/ERK (extracellular signal-related protein kinase) activation. In preclinical acute lung injury (ALI) models, inhibition of MAP2K1 (MAPK kinase 1)/MAP2K2 (MAPK kinase 2) improves measures of ALI. Myeloid cell deletion of MAP2K1 results in sustained MAP2K2 activation and nonresolving ALI, suggesting that MAP2K2 deactivation may be a key driver of ALI resolution. We used human genomic data from the iSPAAR (Identification of SNPs Predisposing to Altered Acute Lung Injury Risk) Consortium to assess genetic variants in MAP2K1 and MAP2K2 for association with mortality from ARDS. To determine the role of MAP2K2 in ALI recovery, we studied mice deficient in Map2k2 (Mek2-/-) and wild-type control mice in ALI models. We identified a MAP2K2 variant that was associated with death in ARDS and MAP2K2 expression. In Pseudomonas aeruginosa ALI, Mek2-/- mice had similar early alveolar neutrophilic recruitment but faster resolution of alveolar neutrophilia and vascular leak. Gene expression analysis revealed a role for MAP2K2 in promoting and sustaining select proinflammatory pathway activation in ALI. Bone marrow chimera studies indicate that leukocyte MAP2K2 is the key regulator of ALI duration. These studies implicate a role for MAP2K2 in ALI duration via transcriptional regulation of inflammatory programming with potential relevance to ARDS. Targeting leukocyte MAP2K2 may be an effective strategy to promote ALI resolution.

Effect of lung pericyte-like cell ablation on the bleomycin model of injury and repair.

We previously showed that pericyte-like cells derived from the FoxD1-lineage contribute to myofibroblasts following bleomycin-induced lung injury. However, their functional significance in lung fibrosis remains unknown. In this study, we used a model of lung pericyte-like cell ablation to test the hypothesis that pericyte-like cell ablation attenuates lung fibrosis in bleomycin-induced lung injury. Lung fibrosis was induced by intratracheal instillation of bleomycin. To ablate pericyte-like cells in the lung, diphtheria toxin (DT) was administered to Foxd1-Cre;Rosa26-iDTR mice at two different phases of bleomycin-induced lung injury. For early ablation, we coadministered bleomycin with DT and harvested mice at days 7 and 21. To test the effect of ablation after acute injury, we delivered DT 7 days after bleomycin administration. We assessed fibrosis by lung hydroxyproline content and semiquantitative analysis of picrosirius red staining. We performed bronchoalveolar lavage to determine cell count and differential. We also interrogated mRNA expression of fibrosis-related genes in whole lung RNA. Compared with DT-insensitive littermates where pericyte-like cells were not ablated, DT-sensitive animals exhibited no difference in fibrosis at day 21 both in the early and late pericyte ablation models. However, early ablation of pericytes reduced acute lung inflammation, as indicated by decreased inflammatory cells. Our data confirm a role for pericytes in regulating pulmonary inflammation in early lung injury.

Sequencing Analysis at 8p23 Identifies Multiple Rare Variants in DLC1 Associated with Sleep-Related Oxyhemoglobin Saturation Level.

Average arterial oxyhemoglobin saturation during sleep (AvSpO2S) is a clinically relevant measure of physiological stress associated with sleep-disordered breathing, and this measure predicts incident cardiovascular disease and mortality. Using high-depth whole-genome sequencing data from the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) project and focusing on genes with linkage evidence on chromosome 8p23,1,2 we observed that six coding and 51 noncoding variants in a gene that encodes the GTPase-activating protein (DLC1) are significantly associated with AvSpO2S and replicated in independent subjects. The combined DLC1 association evidence of discovery and replication cohorts reaches genome-wide significance in European Americans (p = 7.9 × 10-7). A risk score for these variants, built on an independent dataset, explains 0.97% of the AvSpO2S variation and contributes to the linkage evidence. The 51 noncoding variants are enriched in regulatory features in a human lung fibroblast cell line and contribute to DLC1 expression variation. Mendelian randomization analysis using these variants indicates a significant causal effect of DLC1 expression in fibroblasts on AvSpO2S. Multiple sources of information, including genetic variants, gene expression, and methylation, consistently suggest that DLC1 is a gene associated with AvSpO2S.

Pericyte-like cells undergo transcriptional reprogramming and distinct functional adaptations in acute lung injury.

We previously reported on the role of pericyte-like cells as functional sentinel immune cells in lung injury. However, much about the biological role of pericytes in lung injury remains unknown. Lung pericyte-like cells are well-positioned to sense disruption to the epithelial barrier and coordinate local inflammatory responses due to their anatomic niche within the alveoli. In this report, we characterized transcriptional responses and functional changes in pericyte-like cells following activation by alveolar components from injured and uninjured lungs in a mouse model of acute lung injury (ALI). Purified pericyte-like cells from lung digests using PDGFRß as a selection marker were expanded in culture as previously described (1). We induced sterile acute lung injury in mice with recombinant human Fas ligand (rhFasL) instillation followed by mechanical ventilation (1). We then collected bronchoalveolar lavage fluid (BALF) from injured and uninjured mice. Purified pericyte-like cells in culture were exposed to growth media only (control), BALF from uninjured mice, and BALF from injured mice for 6 and 24 hours. RNA collected from these treatment conditions were processed for RNAseq. Targets of interest identified by pathway analysis were validated using in vitro and in vivo assays. We observed robust global transcriptional changes in pericyte-like cells following treatment with uninjured and injured BALF at 6 hours, but this response persisted for 24 hours only after exposure to injured BALF. Functional enrichment analysis of pericytes treated with injured BALF revealed the activation of pro-inflammatory, cell migration, and angiogenesis-related pathways, whereas processes associated with tissue development and cell differentiation were down-regulated. We validated select upregulated targets in the inflammatory, angiogenic, and cell migratory pathways using functional biological assays in vitro and in vivo. We conclude that lung pericyte-like cells are highly responsive to alveolar compartment content from both uninjured and injured lungs, but injured BALF elicits a more sustained response. The inflammatory, angiogenic, and migratory changes exhibited by activated pericyte-like cells underscore the phenotypic plasticity of these specialized stromal cells in the setting of acute lung injury.

Parallel Murine and Human Plaque Proteomics Reveals Pathways of Plaque Rupture.

RATIONALE: Plaque rupture is the proximate cause of most myocardial infarctions and many strokes. However, the molecular mechanisms that precipitate plaque rupture are unknown. OBJECTIVE: By applying proteomic and bioinformatic approaches in mouse models of protease-induced plaque rupture and in ruptured human plaques, we aimed to illuminate biochemical pathways through which proteolysis causes plaque rupture and identify substrates that are cleaved in ruptured plaques. METHODS AND RESULTS: We performed shotgun proteomics analyses of aortas of transgenic mice with macrophage-specific overexpression of urokinase (SR-uPA+/0 mice) and of SR-uPA+/0 bone marrow transplant recipients, and we used bioinformatic tools to evaluate protein abundance and functional category enrichment in these aortas. In parallel, we performed shotgun proteomics and bioinformatics studies on extracts of ruptured and stable areas of freshly harvested human carotid plaques. We also applied a separate protein-analysis method (protein topography and migration analysis platform) to attempt to identify substrates and proteolytic fragments in mouse and human plaque extracts. Approximately 10% of extracted aortic proteins were reproducibly altered in SR-uPA+/0 aortas. Proteases, inflammatory signaling molecules, as well as proteins involved with cell adhesion, the cytoskeleton, and apoptosis, were increased. ECM (Extracellular matrix) proteins, including basement-membrane proteins, were decreased. Approximately 40% of proteins were altered in ruptured versus stable areas of human carotid plaques, including many of the same functional categories that were altered in SR-uPA+/0 aortas. Collagens were minimally altered in SR-uPA+/0 aortas and ruptured human plaques; however, several basement-membrane proteins were reduced in both SR-uPA+/0 aortas and ruptured human plaques. Protein topography and migration analysis platform did not detect robust increases in proteolytic fragments of ECM proteins in either setting. CONCLUSIONS: Parallel studies of SR-uPA+/0 mouse aortas and human plaques identify mechanisms that connect proteolysis with plaque rupture, including inflammation, basement-membrane protein loss, and apoptosis. Basement-membrane protein loss is a prominent feature of ruptured human plaques, suggesting a major role for basement-membrane proteins in maintaining plaque stability.