Biological heterogeneity in idiopathic pulmonary arterial hypertension identified through unsupervised transcriptomic profiling of whole blood

Understanding pulmonary hypertension: the need for an integrative metabolomics and transcriptomics approach

Pulmonary hypertension (PH), characterised by mean pulmonary arterial pressure (mPAP) >20 mm Hg at rest, is a complex pathophysiological disorder associated with multiple clinical conditions. The high prevalence of the disease along with increased mortality and morbidity makes it a global health burden. Despite major advances in understanding the disease pathophysiology, much of the underlying complex molecular mechanism remains to be elucidated. Lack of a robust diagnostic test and specific therapeutic targets also poses major challenges. This review provides a comprehensive update on the dysregulated pathways and promising candidate markers identified in PH patients using the transcriptomics and metabolomics approach. The review also highlights the need of using an integrative multi-omics approach for obtaining insight into the disease at a molecular level. The integrative multi-omics/pan-omics approach envisaged to help in bridging the gap from genotype to phenotype is outlined. Finally, the challenges commonly encountered while conducting omics-driven studies are also discussed.

Immature reticulocyte fraction: A novel biomarker of hemodynamic severity in pulmonary arterial hypertension

Various erythropoietic abnormalities are highly prevalent among patients with pulmonary arterial hypertension (PAH) and associated with worse disease severity. Given the poorly understood yet important roles of dysregulated erythropoiesis and iron metabolism in PAH, we sought to further characterize the hematologic and iron profiles in PAH and their relationship to PAH severity. We recruited 67 patients with PAH and 13 healthy controls. Hemodynamics attained within 1 year of blood sample collection were available for 36 patients. Multiple hematologic, iron, and inflammatory parameters were evaluated for their association with hemodynamics. The subset with hemodynamic data consisted of 29 females (81%). The most common etiologies were idiopathic PAH (47%) and connective tissue disease-related PAH (33%). 19 (53%) had functional class 3 or 4 symptomatology, and 12 (33%) were on triple pulmonary vasodilator therapy. Immature reticulocyte fraction (IRF) had significant positive correlations with mean pulmonary artery (PA) pressure (mPAP) (0.59, p < 0.001), pulmonary vascular resistance (0.52, p = 0.001), and right atrial pressure (0.46, p = 0.005), and significant negative correlations with cardiac index (-0.43, p = 0.009), PA compliance (PAC) (-0.60, p < 0.001), stroke volume index (SVI) (-0.57, p < 0.001), and mixed venous oxygen saturation (-0.51, p = 0.003). IRF correlated with markers of iron deficiency (ID) and erythropoiesis. On multivariable linear regression, IRF was associated with elevated mPAP and reduced SVI and PAC independent of EPO levels, transferrin saturation, and soluble transferrin receptor levels. We identified IRF as a novel and potent biomarker of PAH hemodynamic severity, possibly related to its associations with erythropoiesis, ID, and tissue hypoxia.

B-cells in pulmonary arterial hypertension: friend, foe or bystander?

There is an unmet need for new therapeutic strategies that target alternative pathways to improve the prognosis of patients with pulmonary arterial hypertension (PAH). As immunity has been involved in the development and progression of vascular lesions in PAH, we review the potential contribution of B-cells in its pathogenesis and evaluate the relevance of B-cell-targeted therapies. Circulating B-cell homeostasis is altered in PAH patients, with total B-cell lymphopenia, abnormal subset distribution (expansion of naïve and antibody-secreting cells, reduction of memory B-cells) and chronic activation. B-cells are recruited to the lungs through local chemokine secretion, and activated by several mechanisms: 1) interaction with lung vascular autoantigens through cognate B-cell receptors; 2) costimulatory signals provided by T follicular helper cells (interleukin (IL)-21), type 2 T helper cells and mast cells (IL-4, IL-6 and IL-13); and 3) increased survival signals provided by B-cell activating factor pathways. This activity results in the formation of germinal centres within perivascular tertiary lymphoid organs and in the local production of pathogenic autoantibodies that target the pulmonary vasculature and vascular stabilisation factors (including angiotensin-II/endothelin-1 receptors and bone morphogenetic protein receptors). B-cells also mediate their effects through enhanced production of pro-inflammatory cytokines, reduced anti-inflammatory properties by regulatory B-cells, immunoglobulin (Ig)G-induced complement activation, and IgE-induced mast cell activation. Precision-medicine approaches targeting B-cell immunity are a promising direction for select PAH conditions, as suggested by the efficacy of anti-CD20 therapy in experimental models and a trial of rituximab in systemic sclerosis-associated PAH.

Elderly Patients with Idiopathic Pulmonary Hypertension: Clinical Characteristics, Survival, and Risk Stratification in a Single-Center Prospective Registry

INTRODUCTION

The predictive value of the risk stratification scales in elderly patients with IPAH might differ from that in younger patients. It is unknown whether young and older IPAH patients have the same survival dependence on PAH-specific therapy numbers. The aim of this study was to evaluate the prognostic relevance of risk stratification scales and PAH medication numbers in elderly IPAH patients in comparison with young IPAH patients.

MATERIALS AND METHODS

A total of 119 patients from a prospective single-center PAH registry were divided into group I < 60 years old (n = 89) and group II ≥ 60 years old (n = 30). ESC/ERS, REVEAL, and REVEAL 2.0 risk stratification scores were assessed at baseline, as well as H2FpEF score and survival at follow-up.

RESULTS

During a mean follow-up period of 2.9 years (1.63; 6.0), 42 (35.3%) patients died; at 1, 2, 3, 5, 7, and 10 years, survival was 95%, 88.6%, 78.5%, 61.7%, 48.5%, and 33.7%, respectively. No survival differences were observed between the two groups, despite the use of monotherapy in the elderly patients. The best predictive REVEAL value in elderly patients (IPAH patients ≥ 60 years) was AUC 0.73 (0.56-0.91), p = 0.03; and in patients with LHD comorbidities in the entire cohort, it was AUC 0.73 (0.59-0.87), p < 0.009. Factors independently associated with death in the entire cohort were CKD (p = 0.01, HR 0.2), the right-to-left ventricle dimension ratio (p = 0.0047, HR 5.97), and NT-proBNP > 1400 pg/mL (p = 0.008, HR 3.18).

CONCLUSION

Risk stratification in the elderly IPAH patients requires a fundamentally different approach than that of younger patients, taking into account the initial limitations in physical performance and comorbidities that interfere with current assessment scores. The REVEAL score reliably stratifies patients at any age and LHD comorbidities. The initial monotherapy seems to be reasonable in patients over 60 years. Selection tools for initial combination PAH therapy in older IPAH patients with comorbidities need to be validated in prospective observational studies.

Omada: robust clustering of transcriptomes through multiple testing

BACKGROUND

Cohort studies increasingly collect biosamples for molecular profiling and are observing molecular heterogeneity. High-throughput RNA sequencing is providing large datasets capable of reflecting disease mechanisms. Clustering approaches have produced a number of tools to help dissect complex heterogeneous datasets, but selecting the appropriate method and parameters to perform exploratory clustering analysis of transcriptomic data requires deep understanding of machine learning and extensive computational experimentation. Tools that assist with such decisions without prior field knowledge are nonexistent. To address this, we have developed Omada, a suite of tools aiming to automate these processes and make robust unsupervised clustering of transcriptomic data more accessible through automated machine learning-based functions.

FINDINGS

The efficiency of each tool was tested with 7 datasets characterized by different expression signal strengths to capture a wide spectrum of RNA expression datasets. Our toolkit's decisions reflected the real number of stable partitions in datasets where the subgroups are discernible. Within datasets with less clear biological distinctions, our tools either formed stable subgroups with different expression profiles and robust clinical associations or revealed signs of problematic data such as biased measurements.

CONCLUSIONS

In conclusion, Omada successfully automates the robust unsupervised clustering of transcriptomic data, making advanced analysis accessible and reliable even for those without extensive machine learning expertise. Implementation of Omada is available at http://bioconductor.org/packages/omada/.

Association Study Identified HLA-DQA1 as a Novel Genetic Risk of Systemic Lupus Erythematosus-Associated Pulmonary Arterial Hypertension

OBJECTIVE

Pulmonary arterial hypertension (PAH) is a severe complication of systemic lupus erythematosus (SLE). However, the genetic signatures of SLE-associated PAH have not been well studied. We aimed to identify genetic variants implicated in SLE-associated PAH susceptibility within the major histocompatibility complex (MHC) region and assess the contribution to clinical outcomes.

METHODS

A total of 172 patients with SLE-associated PAH confirmed by right heart catheterization, 1,303 patients with SLE without PAH, and 9,906 healthy controls were included. Deep sequencing of the MHC region was performed to identify alleles, single-nucleotide polymorphisms, and amino acids. We compared patients with SLE-associated PAH with patients with SLE without PAH and healthy controls. Clinical association study was conducted to explore the contribution to phenotypes.

RESULTS

A total of 19,881 genetic variants were identified within the MHC region. HLA-DQA1*03:02 was identified as a novel genetic variant associated with SLE-associated PAH in the discovery cohort (P = 5.68 × 10-12 ) and authenticated in an independent replication cohort (P = 1.30 × 10-9 ). The strongest associated amino acid position was mapped to HLA-DQα1 in the region affecting MHC/peptide-CD4+ T cell receptor affinity and antigen binding. Clinical association study demonstrated that patients with SLE-associated PAH with HLA-DQA1*03:02 had significantly lower rates of target role achievement (P = 0.005) and survival (P = 0.04).

CONCLUSION

This study, based on the largest cohort of SLE-associated PAH, is the first to investigate how MHC region genetic variants contribute to SLE-associated PAH susceptibility. HLA-DQA1*03:02 is a novel genetic risk factor and a prognostic factor in SLE-associated PAH. Patients with SLE with this allele require regular monitoring and careful follow-up for early diagnosis and interventions for potential PAH.

Risk Stratification in Pulmonary Arterial Hypertension, Update and Perspectives

Risk stratification in pulmonary arterial hypertension (PAH) is crucial in assessing patient prognosis. It serves a prominent role in everyday patient care and can be determined using several validated risk assessment scores worldwide. The recently published 2022 European Society of Cardiology (ESC)/European Respiratory Society (ERS) guidelines underline the importance of risk stratification not only at baseline but also during follow-up. Achieving a low-risk status has now become the therapeutic goal, emphasising the importance of personalised therapy. The application of these guidelines is also important in determining the timing for lung transplantation referral. In this review, we summarise the most relevant prognostic factors of PAH as well as the parameters used in PAH risk scores and their evolution in the guidelines over the last decade. Finally, we describe the central role that risk stratification plays in the current guidelines not only in European countries but also in Asian countries.

Serum and Pulmonary Expression Profiles of the Activin Signaling System in Pulmonary Arterial Hypertension

BACKGROUND

Activins are novel therapeutic targets in pulmonary arterial hypertension (PAH). We therefore studied whether key members of the activin pathway could be used as PAH biomarkers.

METHODS

Serum levels of activin A, activin B, α-subunit of inhibin A and B proteins, and the antagonists follistatin and follistatin-like 3 (FSTL3) were measured in controls and in patients with newly diagnosed idiopathic, heritable, or anorexigen-associated PAH (n=80) at baseline and 3 to 4 months after treatment initiation. The primary outcome was death or lung transplantation. Expression patterns of the inhibin subunits, follistatin, FSTL3, Bambi, Cripto, and the activin receptors type I (ALK), type II (ACTRII), and betaglycan were analyzed in PAH and control lung tissues.

RESULTS

Death or lung transplantation occurred in 26 of 80 patients (32.5%) over a median follow-up of 69 (interquartile range, 50-81) months. Both baseline (hazard ratio, 1.001 [95% CI, 1.000-1.001]; P=0.037 and 1.263 [95% CI, 1.049-1.520]; P=0.014, respectively) and follow-up (hazard ratio, 1.003 [95% CI, 1.001-1.005]; P=0.001 and 1.365 [95% CI, 1.185-1.573]; P<0.001, respectively) serum levels of activin A and FSTL3 were associated with transplant-free survival in a model adjusted for age and sex. Thresholds determined by receiver operating characteristic analyses were 393 pg/mL for activin A and 16.6 ng/mL for FSTL3. When adjusted with New York Heart Association functional class, 6-minute walk distance, and N-terminal pro-B-type natriuretic peptide, the hazard ratios for transplant-free survival for baseline activin A <393 pg/mL and FSTL3 <16.6 ng/mL were, respectively, 0.14 (95% CI, 0.03-0.61; P=0.009) and 0.17 (95% CI, 0.06-0.45; P<0.001), and for follow-up measures, 0.23 (95% CI, 0.07-0.78; P=0.019) and 0.27 (95% CI, 0.09-0.78, P=0.015), respectively. Prognostic values of activin A and FSTL3 were confirmed in an independent external validation cohort. Histological analyses showed a nuclear accumulation of the phosphorylated form of Smad2/3, higher immunoreactivities for ACTRIIB, ALK2, ALK4, ALK5, ALK7, Cripto, and FSTL3 in vascular endothelial and smooth muscle layers, and lower immunostaining for inhibin-α and follistatin.

CONCLUSIONS

These findings offer new insights into the activin signaling system in PAH and show that activin A and FSTL3 are prognostic biomarkers for PAH.

Artificial intelligence-assisted selection and efficacy prediction of antineoplastic strategies for precision cancer therapy

The rapid development of artificial intelligence (AI) technologies in the context of the vast amount of collectable data obtained from high-throughput sequencing has led to an unprecedented understanding of cancer and accelerated the advent of a new era of clinical oncology with a tone of precision treatment and personalized medicine. However, the gains achieved by a variety of AI models in clinical oncology practice are far from what one would expect, and in particular, there are still many uncertainties in the selection of clinical treatment options that pose significant challenges to the application of AI in clinical oncology. In this review, we summarize emerging approaches, relevant datasets and open-source software of AI and show how to integrate them to address problems from clinical oncology and cancer research. We focus on the principles and procedures for identifying different antitumor strategies with the assistance of AI, including targeted cancer therapy, conventional cancer therapy, and cancer immunotherapy. In addition, we also highlight the current challenges and directions of AI in clinical oncology translation. Overall, we hope this article will provide researchers and clinicians with a deeper understanding of the role and implications of AI in precision cancer therapy, and help AI move more quickly into accepted cancer guidelines.

Pursuing functional biomarkers in complex disease: Focus on pulmonary arterial hypertension

A major gap in diagnosis, classification, risk stratification, and prediction of therapeutic response exists in pulmonary arterial hypertension (PAH), driven in part by a lack of functional biomarkers that are also disease-specific. In this regard, leveraging big data-omics analyses using innovative approaches that integrate network medicine and machine learning correlated with clinically useful indices or risk stratification scores is an approach well-positioned to advance PAH precision medicine. For example, machine learning applied to a panel of 48 cytokines, chemokines, and growth factors could prognosticate PAH patients with immune-dominant subphenotypes at elevated or low-risk for mortality. Here, we discuss strengths and weaknesses of the most current studies evaluating omics-derived biomarkers in PAH. Progress in this field is offset by studies with small sample size, pervasive limitations in bioinformatics, and lack of standardized methods for data processing and interpretation. Future success in this field, in turn, is likely to hinge on mechanistic validation of data outputs in order to couple functional biomarker data with target-specific therapeutics in clinical practice.

Tm4sf1-marked Endothelial Subpopulation Is Dysregulated in Pulmonary Arterial Hypertension

The identification and role of endothelial progenitor cells in pulmonary arterial hypertension (PAH) remain controversial. Single-cell omics analysis can shed light on endothelial progenitor cells and their potential contribution to PAH pathobiology. We aim to identify endothelial cells that may have stem/progenitor potential in rat lungs and assess their relevance to PAH. Differential expression, gene set enrichment, cell-cell communication, and trajectory reconstruction analyses were performed on lung endothelial cells from single-cell RNA sequencing of Sugen-hypoxia, monocrotaline, and control rats. Relevance to human PAH was assessed in multiple independent blood and lung transcriptomic data sets. Rat lung endothelial cells were visualized by immunofluorescence in situ, analyzed by flow cytometry, and assessed for tubulogenesis in vitro. A subpopulation of endothelial cells (endothelial arterial type 2 [EA2]) marked by Tm4sf1 (transmembrane 4 L six family member 1), a gene strongly implicated in cancer, harbored a distinct transcriptomic signature enriched for angiogenesis and CXCL12 signaling. Trajectory analysis predicted that EA2 has a less differentiated state compared with other endothelial subpopulations. Analysis of independent data sets revealed that TM4SF1 is downregulated in lungs and endothelial cells from patients and PAH models, is a marker for hematopoietic stem cells, and is upregulated in PAH circulation. TM4SF1+CD31+ rat lung endothelial cells were visualized in distal pulmonary arteries, expressed hematopoietic marker CD45, and formed tubules in coculture with lung fibroblasts. Our study uncovered a novel Tm4sf1-marked subpopulation of rat lung endothelial cells that may have stem/progenitor potential and demonstrated its relevance to PAH. Future studies are warranted to further elucidate the role of EA2 and Tm4sf1 in PAH.

Tensions in Taxonomies: Current Understanding and Future Directions in the Pathobiologic Basis and Treatment of Group 1 and Group 3 Pulmonary Hypertension

In the over 100 years since the recognition of pulmonary hypertension (PH), immense progress and significant achievements have been made with regard to understanding the pathophysiology of the disease and its treatment. These advances have been mostly in idiopathic pulmonary arterial hypertension (IPAH), which was classified as Group 1 Pulmonary Hypertension (PH) at the Second World Symposia on PH in 1998. However, the pathobiology of PH due to chronic lung disease, classified as Group 3 PH, remains poorly understood and its treatments thus remain limited. We review the history of the classification of the five groups of PH and aim to provide a state-of-the-art review of the understanding of the pathogenesis of Group 1 PH and Group 3 PH including insights gained from novel high-throughput omics technologies that have revealed heterogeneities within these categories as well as similarities between them. Leveraging the substantial gains made in understanding the genomics, epigenomics, proteomics, and metabolomics of PAH to understand the full spectrum of the complex, heterogeneous disease of PH is needed. Multimodal omics data as well as supervised and unbiased machine learning approaches after careful consideration of the powerful advantages as well as of the limitations and pitfalls of these technologies could lead to earlier diagnosis, more precise risk stratification, better predictions of disease response, new sub-phenotype groupings within types of PH, and identification of shared pathways between PAH and other types of PH that could lead to new treatment targets. © 2023 American Physiological Society. Compr Physiol 13:4295-4319, 2023.

GSDMD deficiency ameliorates hyperoxia-induced BPD and ROP in neonatal mice

Bronchopulmonary dysplasia (BPD) and retinopathy of prematurity (ROP) are among the most common morbidities affecting extremely premature infants who receive oxygen therapy. Many clinical studies indicate that BPD is associated with advanced ROP. However, the mechanistic link between hyperoxia, BPD, and ROP remains to be explored. Gasdermin D (GSDMD) is a key executor of inflammasome-induced pyroptosis and inflammation. Inhibition of GSDMD has been shown to attenuate hyperoxia-induced BPD and brain injury in neonatal mice. The objective of this study was to further define the mechanistic roles of GSDMD in the pathogenesis of hyperoxia-induced BPD and ROP in mouse models. Here we show that global GSDMD knockout (GSDMD-KO) protects against hyperoxia-induced BPD by reducing macrophage infiltration, improving alveolarization and vascular development, and decreasing cell death. In addition, GSDMD deficiency prevented hyperoxia-induced ROP by reducing vasoobliteration and neovascularization, improving thinning of multiple retinal tissue layers, and decreasing microglial activation. RNA sequencing analyses of lungs and retinas showed that similar genes, including those from inflammatory, cell death, tissue remodeling, and tissue and vascular developmental signaling pathways, were induced by hyperoxia and impacted by GSDMD-KO in both models. These data highlight the importance of GSDMD in the pathogenesis of BPD and ROP and suggest that targeting GSDMD may be beneficial in preventing and treating BPD and ROP in premature infants.

An explorative metabolomic analysis of the endothelium in pulmonary hypertension

Pulmonary hypertension (PH) is classified into five clinical diagnostic groups, including group 1 [idiopathic pulmonary arterial hypertension (IPAH) and connective tissue disease-associated PAH (CTD-aPAH)] and group 4 (chronic thromboembolic pulmonary hypertension (CTEPH)). PH is a progressive, life-threatening, incurable disease. The pathological mechanisms underlying PH remain elusive; recent evidence has revealed that abnormal metabolic activities in the endothelium may play a crucial role. This research introduces a novel approach for studying PH endothelial function, building on the genome-scale metabolic reconstruction of the endothelial cell (EC) to investigate intracellular metabolism. We demonstrate that the intracellular metabolic activities of ECs in PH patients cluster into four phenotypes independent of the PH diagnosis. Notably, the disease severity differs significantly between the metabolic phenotypes, suggesting their clinical relevance. The significant metabolic differences between the PH phenotypes indicate that they may require different therapeutic interventions. In addition, diagnostic capabilities enabling their identification is warranted to investigate whether this opens a novel avenue of precision medicine.

Autoimmunity Is a Significant Feature of Idiopathic Pulmonary Arterial Hypertension

Rationale: Autoimmunity is believed to play a role in idiopathic pulmonary arterial hypertension (IPAH). It is not clear whether this is causative or a bystander of disease and if it carries any prognostic or treatment significance. Objectives: To study autoimmunity in IPAH using a large cross-sectional cohort. Methods: Assessment of the circulating immune cell phenotype was undertaken using flow cytometry, and the profile of serum immunoglobulins was generated using a standardized multiplex array of 19 clinically validated autoantibodies in 473 cases and 946 control subjects. Additional glutathione S-transferase fusion array and ELISA data were used to identify a serum autoantibody to BMPR2 (bone morphogenetic protein receptor type 2). Clustering analyses and clinical correlations were used to determine associations between immunogenicity and clinical outcomes. Measurements and Main Results: Flow cytometric immune profiling demonstrates that IPAH is associated with an altered humoral immune response in addition to raised IgG3. Multiplexed autoantibodies were significantly raised in IPAH, and clustering demonstrated three distinct clusters: "high autoantibody," "low autoantibody," and a small "intermediate" cluster exhibiting high concentrations of ribonucleic protein complex. The high-autoantibody cluster had worse hemodynamics but improved survival. A small subset of patients demonstrated immunoglobulin reactivity to BMPR2. Conclusions: This study establishes aberrant immune regulation and presence of autoantibodies as key features in the profile of a significant proportion of patients with IPAH and is associated with clinical outcomes.

The Landscape of Noncoding RNA in Pulmonary Hypertension

The transcriptome of pulmonary hypertension (PH) is complex and highly genetically heterogeneous, with noncoding RNA transcripts playing crucial roles. The majority of RNAs in the noncoding transcriptome are long noncoding RNAs (lncRNAs) with less circular RNAs (circRNAs), which are two characteristics gaining increasing attention in the forefront of RNA research field. These noncoding transcripts (especially lncRNAs and circRNAs) exert important regulatory functions in PH and emerge as potential disease biomarkers and therapeutic targets. Recent technological advancements have established great momentum for discovery and functional characterization of ncRNAs, which include broad transcriptome sequencing such as bulk RNA-sequence, single-cell and spatial transcriptomics, and RNA-protein/RNA interactions. In this review, we summarize the current research on the classification, biogenesis, and the biological functions and molecular mechanisms of these noncoding RNAs (ncRNAs) involved in the pulmonary vascular remodeling in PH. Furthermore, we highlight the utility and challenges of using these ncRNAs as biomarkers and therapeutics in PH.

Inflammasome Activation in Pulmonary Arterial Hypertension

Inflammasomes are multi-protein complexes that sense both infectious and sterile inflammatory stimuli, launching a cascade of responses to propagate danger signaling throughout an affected tissue. Recent studies have implicated inflammasome activation in a variety of pulmonary diseases, including pulmonary arterial hypertension (PAH). Indeed, the end-products of inflammasome activation, including interleukin (IL)-1β, IL-18, and lytic cell death (“pyroptosis”) are all key biomarkers of PAH, and are potentially therapeutic targets for human disease. This review summarizes current knowledge of inflammasome activation in immune and vascular cells of the lung, with a focus on the role of these pathways in the pathogenesis of PAH. Special emphasis is placed on areas of potential drug development focused on inhibition of inflammasomes and their downstream effectors.