Reversal of pulmonary veno-occlusive disease phenotypes by inhibition of the integrated stress response.

Pulmonary veno-occlusive disease (PVOD) is a rare form of pulmonary hypertension arising from EIF2AK4 gene mutations or mitomycin C (MMC) administration. The lack of effective PVOD therapies is compounded by a limited understanding of the mechanisms driving vascular remodeling in PVOD. Here we show that administration of MMC in rats mediates activation of protein kinase R (PKR) and the integrated stress response (ISR), which leads to the release of the endothelial adhesion molecule vascular endothelial (VE) cadherin (VE-Cad) in complex with RAD51 to the circulation, disruption of endothelial barrier and vascular remodeling. Pharmacological inhibition of PKR or ISR attenuates VE-Cad depletion, elevation of vascular permeability and vascular remodeling instigated by MMC, suggesting potential clinical intervention for PVOD. Finally, the severity of PVOD phenotypes was increased by a heterozygous BMPR2 mutation that truncates the carboxyl tail of the receptor BMPR2, underscoring the role of deregulated bone morphogenetic protein signaling in the development of PVOD.

Genetics and precision genomics approaches to pulmonary hypertension.

Considerable progress has been made in the genomics of pulmonary arterial hypertension (PAH) since the 6th World Symposium on Pulmonary Hypertension, with the identification of rare variants in several novel genes, as well as common variants that confer a modest increase in PAH risk. Gene and variant curation by an expert panel now provides a robust framework for knowing which genes to test and how to interpret variants in clinical practice. We recommend that genetic testing be offered to specific subgroups of symptomatic patients with PAH, and to children with certain types of group 3 pulmonary hypertension (PH). Testing of asymptomatic family members and the use of genetics in reproductive decision-making require the involvement of genetics experts. Large cohorts of PAH patients with biospecimens now exist and extension to non-group 1 PH has begun. However, these cohorts are largely of European origin; greater diversity will be essential to characterise the full extent of genomic variation contributing to PH risk and treatment responses. Other types of omics data are also being incorporated. Furthermore, to advance gene- and pathway-specific care and targeted therapies, gene-specific registries will be essential to support patients and their families and to lay the foundation for genetically informed clinical trials. This will require international outreach and collaboration between patients/families, clinicians and researchers. Ultimately, harmonisation of patient-derived biospecimens, clinical and omic information, and analytic approaches will advance the field.

Integrative Multiomics in the Lung Reveals a Protective Role of Asporin in Pulmonary Arterial Hypertension.

BACKGROUND: Integrative multiomics can elucidate pulmonary arterial hypertension (PAH) pathobiology, but procuring human PAH lung samples is rare. METHODS: We leveraged transcriptomic profiling and deep phenotyping of the largest multicenter PAH lung biobank to date (96 disease and 52 control) by integration with clinicopathologic data, genome-wide association studies, Bayesian regulatory networks, single-cell transcriptomics, and pharmacotranscriptomics. RESULTS: We identified 2 potentially protective gene network modules associated with vascular cells, and we validated ASPN, coding for asporin, as a key hub gene that is upregulated as a compensatory response to counteract PAH. We found that asporin is upregulated in lungs and plasma of multiple independent PAH cohorts and correlates with reduced PAH severity. We show that asporin inhibits proliferation and transforming growth factor-ß/phosphorylated SMAD2/3 signaling in pulmonary artery smooth muscle cells from PAH lungs. We demonstrate in Sugen-hypoxia rats that ASPN knockdown exacerbated PAH and recombinant asporin attenuated PAH. CONCLUSIONS: Our integrative systems biology approach to dissect the PAH lung transcriptome uncovered asporin as a novel protective target with therapeutic potential in PAH.

Rare Sequence Variation Underlying Suspected Familial Cerebral Small-Vessel Disease.

BACKGROUND: Cerebral small-vessel disease (cSVD) is the leading monogenic cause of stroke. Despite genetic screening in routine diagnosis, many cases remain without a known causative variant. Using a cohort with suspected familial cSVD and whole-genome sequencing, we screened for variants in genes associated with monogenic cSVD and searched for novel variants associated with the disease. METHODS AND RESULTS: Rare variants were identified in whole-genome sequencing data from the NBR (National Institute for Health Research BioResource Rare Disease) study. Pathogenic variants in known monogenic cSVD genes were identified. Gene-based burden tests and family analysis were performed to identify novel variants associated with familial cSVD. A total of 257 suspected cSVD cases (mean ± SD age, 56.2 ± 16.1 years), and 13 086 controls with other nonstroke diseases (5874 [44.9%] men) were studied. A total of 8.9% of the cases carried a variant in known cSVD genes. Excluding these known causes, 23.6% of unrelated subjects with cSVD carried predicted deleterious variants in the Genomics England gene panel, but no association was found with cSVD in burden tests. We identified potential associations with cSVD in noncoding genes, including RP4-568F9.3 (adjusted P = 7.1 × 10-25), RP3-466I7.1 (adjusted P = 8.9 × 10-16), and ZNF209P (adjusted P = 1.0 × 10-15), and matrisomal genes (adjusted P = 5.1 × 10-6), including FAM20C, INHA, LAMC1, and VWA5B2. CONCLUSIONS: Predicted deleterious variants in known cSVD genes were present in 23.6% of unrelated cases with cSVD, but none of the genes were associated with the disease. Rare variants in noncoding and matrisomal genes could potentially contribute to cSVD development. These genes could play a role in tissue development and brain endothelial cell function. However, further studies are needed to confirm their pathophysiological roles.

SEMA3G regulates BMP9 inhibition of VEGF-mediated migration and network formation in pulmonary endothelial cells.

AIMS: Bone morphogenetic protein-9 (BMP9) is critical for bone morphogenetic protein receptor type-2 (BMPR2) signalling in pulmonary vascular endothelial cells. Furthermore, human genetics studies support the central role of disrupted BMPR2 mediated BMP9 signalling in vascular endothelial cells in the initiation of pulmonary arterial hypertension (PAH). In addition, loss-of-function mutations in BMP9 have been identified in PAH patients. BMP9 is considered to play an important role in vascular homeostasis and quiescence. METHODS AND RESULTS: We identified a novel BMP9 target as the class-3 semaphorin, SEMA3G. Although originally identified as playing a role in neuronal development, class-3 semaphorins may have important roles in endothelial function. Here we show that BMP9 transcriptional regulation of SEMA3G occurs via ALK1 and the canonical Smad pathway, requiring both Smad1 and Smad5. Knockdown studies demonstrated redundancy between type-2 receptors in that BMPR2 and ACTR2A were compensatory. Increased SEMA3G expression by BMP9 was found to be regulated by the transcription factor, SOX17. Moreover, we observed that SEMA3G regulates VEGF signalling by inhibiting VEGFR2 phosphorylation and that VEGF, in contrast to BMP9, negatively regulated SEMA3G transcription. Functional endothelial cell assays of VEGF-mediated migration and network formation revealed that BMP9 inhibition of VEGF was abrogated by SEMA3G knockdown. Conversely, treatment with recombinant SEMA3G partially mimicked the inhibitory action of BMP9 in these assays. CONCLUSIONS: This study provides further evidence for the anti-angiogenic role of BMP9 in microvascular endothelial cells and these functions are mediated at least in part via SOX17 and SEMA3G induction.

Shared and Distinct Genomics of Chronic Thromboembolic Pulmonary Hypertension and Pulmonary Embolism.

Rationale: Chronic thromboembolic pulmonary hypertension involves the formation and nonresolution of thrombus, dysregulated inflammation, angiogenesis, and the development of a small-vessel vasculopathy. Objectives: We aimed to establish the genetic basis of chronic thromboembolic pulmonary hypertension to gain insight into its pathophysiological contributors. Methods: We conducted a genome-wide association study on 1,907 European cases and 10,363 European control subjects. We coanalyzed our results with existing results from genome-wide association studies on deep vein thrombosis, pulmonary embolism, and idiopathic pulmonary arterial hypertension. Measurements and Main Results: Our primary association study revealed genetic associations at the ABO, FGG, F11, MYH7B, and HLA-DRA loci. Through our coanalysis, we demonstrate further associations with chronic thromboembolic pulmonary hypertension at the F2, TSPAN15, SLC44A2, and F5 loci but find no statistically significant associations shared with idiopathic pulmonary arterial hypertension. Conclusions: Chronic thromboembolic pulmonary hypertension is a partially heritable polygenic disease, with related though distinct genetic associations with pulmonary embolism and deep vein thrombosis.

Mendelian Randomization Study With Clinical Follow-Up Links Metabolites to Risk and Severity of Pulmonary Arterial Hypertension.

BACKGROUND: Pulmonary arterial hypertension (PAH) exhibits phenotypic heterogeneity and variable response to therapy. The metabolome has been implicated in the pathogenesis of PAH, but previous works have lacked power to implicate specific metabolites. Mendelian randomization (MR) is a method for causal inference between exposures and outcomes. METHODS AND RESULTS: Using genome-wide association study summary statistics, we implemented MR analysis to test for potential causal relationships between serum concentration of 575 metabolites and PAH. Five metabolites were causally associated with the risk of PAH after multiple testing correction. Next, we measured serum concentration of candidate metabolites in an independent clinical cohort of 449 patients with PAH to check whether metabolite concentrations are correlated with markers of disease severity. Of the 5 candidates nominated by our MR work, serine was negatively associated and homostachydrine was positively associated with clinical severity of PAH via direct measurement in this independent clinical cohort. Finally we used conditional and orthogonal approaches to explore the biology underlying our lead metabolites. Rare variant burden testing was carried out using whole exome sequencing data from 578 PAH cases and 361 675 controls. Multivariable MR is an extension of MR that uses a single set of instrumental single-nucleotide polymorphisms to measure multiple exposures; multivariable MR is used to determine interdependence between the effects of different exposures on a single outcome. Rare variant analysis demonstrated that loss-of-function mutations within activating transcription factor 4, a transcription factor responsible for upregulation of serine synthesis under conditions of serine starvation, are associated with higher risk for PAH. Homostachydrine is a xenobiotic metabolite that is structurally related to l-proline betaine, which has previously been linked to modulation of inflammation and tissue remodeling in PAH. Our multivariable MR analysis suggests that the effect of l-proline betaine is actually mediated indirectly via homostachydrine. CONCLUSIONS: Our data present a method for study of the metabolome in the context of PAH, and suggests several candidates for further evaluation and translational research.

Blood DNA methylation profiling identifies cathepsin Z dysregulation in pulmonary arterial hypertension.

Pulmonary arterial hypertension (PAH) is characterised by pulmonary vascular remodelling causing premature death from right heart failure. Established DNA variants influence PAH risk, but susceptibility from epigenetic changes is unknown. We addressed this through epigenome-wide association study (EWAS), testing 865,848 CpG sites for association with PAH in 429 individuals with PAH and 1226 controls. Three loci, at Cathepsin Z (CTSZ, cg04917472), Conserved oligomeric Golgi complex 6 (COG6, cg27396197), and Zinc Finger Protein 678 (ZNF678, cg03144189), reached epigenome-wide significance (p < 10-7) and are hypermethylated in PAH, including in individuals with PAH at 1-year follow-up. Of 16 established PAH genes, only cg10976975 in BMP10 shows hypermethylation in PAH. Hypermethylation at CTSZ is associated with decreased blood cathepsin Z mRNA levels. Knockdown of CTSZ expression in human pulmonary artery endothelial cells increases caspase-3/7 activity (p < 10-4). DNA methylation profiles are altered in PAH, exemplified by the pulmonary endothelial function modifier CTSZ, encoding protease cathepsin Z.

Allele-specific control of rodent and human lncRNA KMT2E-AS1 promotes hypoxic endothelial pathology in pulmonary hypertension.

Hypoxic reprogramming of vasculature relies on genetic, epigenetic, and metabolic circuitry, but the control points are unknown. In pulmonary arterial hypertension (PAH), a disease driven by hypoxia inducible factor (HIF)-dependent vascular dysfunction, HIF-2α promoted expression of neighboring genes, long noncoding RNA (lncRNA) histone lysine N-methyltransferase 2E-antisense 1 (KMT2E-AS1) and histone lysine N-methyltransferase 2E (KMT2E). KMT2E-AS1 stabilized KMT2E protein to increase epigenetic histone 3 lysine 4 trimethylation (H3K4me3), driving HIF-2α-dependent metabolic and pathogenic endothelial activity. This lncRNA axis also increased HIF-2α expression across epigenetic, transcriptional, and posttranscriptional contexts, thus promoting a positive feedback loop to further augment HIF-2α activity. We identified a genetic association between rs73184087, a single-nucleotide variant (SNV) within a KMT2E intron, and disease risk in PAH discovery and replication patient cohorts and in a global meta-analysis. This SNV displayed allele (G)-specific association with HIF-2α, engaged in long-range chromatin interactions, and induced the lncRNA-KMT2E tandem in hypoxic (G/G) cells. In vivo, KMT2E-AS1 deficiency protected against PAH in mice, as did pharmacologic inhibition of histone methylation in rats. Conversely, forced lncRNA expression promoted more severe PH. Thus, the KMT2E-AS1/KMT2E pair orchestrates across convergent multi-ome landscapes to mediate HIF-2α pathobiology and represents a key clinical target in pulmonary hypertension.

Reversal of pulmonary veno-occlusive disease phenotypes by inhibition of the integrated stress response.

Pulmonary veno-occlusive disease (PVOD) is a rare form of pulmonary hypertension arising from EIF2AK4 gene mutations or mitomycin C (MMC) administration. The lack of effective PVOD therapies is compounded by a limited understanding of the mechanisms driving the vascular remodeling in PVOD. We show that the administration of MMC in rats mediates the activation of protein kinase R (PKR) and the integrated stress response (ISR), which lead to the release of the endothelial adhesion molecule VE-Cadherin in the complex with Rad51 to the circulation, disruption of endothelial barrier, and vascular remodeling. Pharmacological inhibition of PKR or ISR attenuates the depletion of VE-Cadherin, elevation of vascular permeability, and vascular remodeling instigated by MMC, suggesting potential clinical intervention for PVOD. Finally, the severity of PVOD phenotypes was increased by a heterozygous BMPR2 mutation that truncates the carboxyl tail of BMPR2, underscoring the role of deregulated BMP signal in the development of PVOD.

Defining the clinical validity of genes reported to cause pulmonary arterial hypertension.

PURPOSE: Pulmonary arterial hypertension (PAH) is a rare, progressive vasculopathy with significant cardiopulmonary morbidity and mortality. Genetic testing is currently recommended for adults diagnosed with heritable, idiopathic, anorexigen-, hereditary hemorrhagic telangiectasia-, and congenital heart disease-associated PAH, PAH with overt features of venous/capillary involvement, and all children diagnosed with PAH. Variants in at least 27 genes have putative evidence for PAH causality. Rigorous assessment of the evidence is needed to inform genetic testing. METHODS: An international panel of experts in PAH applied a semi-quantitative scoring system developed by the NIH Clinical Genome Resource to classify the relative strength of evidence supporting PAH gene-disease relationships based on genetic and experimental evidence. RESULTS: Twelve genes (BMPR2, ACVRL1, ATP13A3, CAV1, EIF2AK4, ENG, GDF2, KCNK3, KDR, SMAD9, SOX17, and TBX4) were classified as having definitive evidence and 3 genes (ABCC8, GGCX, and TET2) with moderate evidence. Six genes (AQP1, BMP10, FBLN2, KLF2, KLK1, and PDGFD) were classified as having limited evidence for causal effects of variants. TOPBP1 was classified as having no known PAH relationship. Five genes (BMPR1A, BMPR1B, NOTCH3, SMAD1, and SMAD4) were disputed because of a paucity of genetic evidence over time. CONCLUSION: We recommend that genetic testing includes all genes with definitive evidence and that caution be taken in the interpretation of variants identified in genes with moderate or limited evidence. Genes with no known evidence for PAH or disputed genes should not be included in genetic testing.

Risk Factors for Longer-Term Mortality in Discharged Patients with Dementia and SARS-CoV-2 Infection: A Matched Case-Control Study.

BACKGROUND: Persisting symptoms and increased mortality after SARS-CoV-2 infection has been described in COVID-19 survivors. OBJECTIVE: We examined longer-term mortality in patients with dementia and SARS-CoV-2 infection. METHODS: A retrospective matched case-control study of 165 patients with dementia who survived an acute hospital admission with COVID-19 infection, and 1325 patients with dementia who survived a hospital admission but without SARS-CoV-2 infection. Potential risk factors investigated included socio-demographic factors, clinical features, and results of investigations. Data were fitted using a Cox proportional hazard model. RESULTS: Compared to patients with dementia but without SARS-CoV-2 infection, people with dementia and SARS-CoV-2 infection had a 4.4-fold risk of death (adjusted hazard ratio [aHR] = 4.44, 95% confidence interval [CI] 3.13-6.30) even beyond the acute phase of infection. This excess mortality could be seen up to 125 days after initial recovery but was not elevated beyond this time. Risk factors for COVID-19-associated mortality included prescription of antipsychotics (aHR = 3.06, 95% CI 1.40-6.69) and benzodiazepines (aHR = 3.00, 95% CI 1.28-7.03). Abnormalities on investigation associated with increased mortality included high white cell count (aHR = 1.21, 95% CI 1.04-1.39), higher absolute neutrophil count (aHR = 1.28, 95% CI 1.12-1.46), higher C-reactive protein (aHR = 1.01, 95% CI 1.00-1.02), higher serum sodium (aHR = 1.09, 95% CI 1.01-1.19), and higher ionized calcium (aHR = 1.03, 95% CI 1.00-1.06). The post-acute COVID mortality could be modeled for the first 120 days after recovery with a balanced accuracy of 87.2%. CONCLUSION: We found an increased mortality in patients with dementia beyond the acute phase of illness. We identified several investigation results associated with increased mortality, and increased mortality in patients prescribed antipsychotics or benzodiazepines.

Reduced circulating BMP9 and pBMP10 in hospitalized COVID-19 patients.

Similar to other causes of acute respiratory distress syndrome, coronavirus disease 2019 (COVID-19) is characterized by the aberrant expression of vascular injury biomarkers. We present the first report that circulating plasma bone morphogenetic proteins (BMPs), BMP9 and pBMP10, involved in vascular protection, are reduced in hospitalized patients with COVID-19.

Integrative Multiomics to Dissect the Lung Transcriptional Landscape of Pulmonary Arterial Hypertension.

Pulmonary arterial hypertension (PAH) remains an incurable and often fatal disease despite currently available therapies. Multiomics systems biology analysis can shed new light on PAH pathobiology and inform translational research efforts. Using RNA sequencing on the largest PAH lung biobank to date (96 disease and 52 control), we aim to identify gene co-expression network modules associated with PAH and potential therapeutic targets. Co-expression network analysis was performed to identify modules of co-expressed genes which were then assessed for and prioritized by importance in PAH, regulatory role, and therapeutic potential via integration with clinicopathologic data, human genome-wide association studies (GWAS) of PAH, lung Bayesian regulatory networks, single-cell RNA-sequencing data, and pharmacotranscriptomic profiles. We identified a co-expression module of 266 genes, called the pink module, which may be a response to the underlying disease process to counteract disease progression in PAH. This module was associated not only with PAH severity such as increased PVR and intimal thickness, but also with compensated PAH such as lower number of hospitalizations, WHO functional class and NT-proBNP. GWAS integration demonstrated the pink module is enriched for PAH-associated genetic variation in multiple cohorts. Regulatory network analysis revealed that BMPR2 regulates the main target of FDA-approved riociguat, GUCY1A2, in the pink module. Analysis of pathway enrichment and pink hub genes (i.e. ANTXR1 and SFRP4) suggests the pink module inhibits Wnt signaling and epithelial-mesenchymal transition. Cell type deconvolution showed the pink module correlates with higher vascular cell fractions (i.e. myofibroblasts). A pharmacotranscriptomic screen discovered ubiquitin-specific peptidases (USPs) as potential therapeutic targets to mimic the pink module signature. Our multiomics integrative study uncovered a novel gene subnetwork associated with clinicopathologic severity, genetic risk, specific vascular cell types, and new therapeutic targets in PAH. Future studies are warranted to investigate the role and therapeutic potential of the pink module and targeting USPs in PAH.

Genetic counselling and testing in pulmonary arterial hypertension: a consensus statement on behalf of the International Consortium for Genetic Studies in PAH.

Pulmonary arterial hypertension (PAH) is a rare disease that can be caused by (likely) pathogenic germline genomic variants. In addition to the most prevalent disease gene, BMPR2 (bone morphogenetic protein receptor 2), several genes, some belonging to distinct functional classes, are also now known to predispose to the development of PAH. As a consequence, specialist and non-specialist clinicians and healthcare professionals are increasingly faced with a range of questions regarding the need for, approaches to and benefits/risks of genetic testing for PAH patients and/or related family members. We provide a consensus-based approach to recommendations for genetic counselling and assessment of current best practice for disease gene testing. We provide a framework and the type of information to be provided to patients and relatives through the process of genetic counselling, and describe the presently known disease causal genes to be analysed. Benefits of including molecular genetic testing within the management protocol of patients with PAH include the identification of individuals misclassified by other diagnostic approaches, the optimisation of phenotypic characterisation for aggregation of outcome data, including in clinical trials, and importantly through cascade screening, the detection of healthy causal variant carriers, to whom regular assessment should be offered.

Molecular Function and Contribution of TBX4 in Development and Disease.

Over the past decade, recognition of the profound impact of the TBX4 (T-box 4) gene, which encodes a member of the evolutionarily conserved family of T-box-containing transcription factors, on respiratory diseases has emerged. The developmental importance of TBX4 is emphasized by the association of TBX4 variants with congenital disorders involving respiratory and skeletal structures; however, the exact role of TBX4 in human development remains incompletely understood. Here, we discuss the developmental, tissue-specific, and pathological TBX4 functions identified through human and animal studies and review the published TBX4 variants resulting in variable disease phenotypes. We also outline future research directions to fill the gaps in our understanding of TBX4 function and of how TBX4 disruption affects development.

Biallelic variants in the calpain regulatory subunit CAPNS1 cause pulmonary arterial hypertension.

Purpose: The aim of this study was to identify the monogenic cause of pulmonary arterial hypertension (PAH), a multifactorial and often fatal disease, in 2 unrelated consanguine families. Methods: We performed exome sequencing and validated variant pathogenicity by whole-blood RNA and protein expression analysis in both families. Further RNA sequencing of preserved lung tissue was performed to investigate the consequences on selected genes that are involved in angiogenesis, proliferation, and apoptosis. Results: We identified 2 rare biallelic variants in CAPNS1, encoding the regulatory subunit of calpain. The variants cosegregated with PAH in the families. Both variants lead to loss of function (LoF), which is demonstrated by aberrant splicing resulting in the complete absence of the CAPNS1 protein in affected patients. No other LoF CAPNS1 variant was identified in the genome data of more than 1000 patients with unresolved PAH. Conclusion: The calpain holoenzyme was previously linked to pulmonary vascular development and progression of PAH in patients. We demonstrated that biallelic LoF variants in CAPNS1 can cause idiopathic PAH by the complete absence of CAPNS1 protein. Screening of this gene in patients who are affected by PAH, especially with suspected autosomal recessive inheritance, should be considered.

A pilot study to examine association of BMI with functional class and 6 min walk distance in idiopathic and heritable PAH: Possible association with estrogen metabolism.

The hypothesis that a relationship exists between body mass index (BMI), functional class, and 6 min walk distance (6MWD) in Group 1-pulmonary arterial hypertension (PAH) was examined. Analysis of data from the UK National Cohort Study for heritable pulmonary arterial/idiopathic PAH suggests increased BMI is a predictor of worse functional class and shorter 6MWD; increased body-weight in mice and man may be associated with increased estrogen metabolism.

First Genotype-Phenotype Study in TBX4 Syndrome: Gain-of-Function Mutations Causative for Lung Disease.

Rationale: Despite the increased recognition of TBX4 (T-BOX transcription factor 4)-associated pulmonary arterial hypertension (PAH), genotype-phenotype associations are lacking and may provide important insights. Objectives: To compile and functionally characterize all TBX4 variants reported to date and undertake a comprehensive genotype-phenotype analysis. Methods: We assembled a multicenter cohort of 137 patients harboring monoallelic TBX4 variants and assessed the pathogenicity of missense variation (n = 42) using a novel luciferase reporter assay containing T-BOX binding motifs. We sought genotype-phenotype correlations and undertook a comparative analysis with patients with PAH with BMPR2 (Bone Morphogenetic Protein Receptor type 2) causal variants (n = 162) or no identified variants in PAH-associated genes (n = 741) genotyped via the National Institute for Health Research BioResource-Rare Diseases. Measurements and Main Results: Functional assessment of TBX4 missense variants led to the novel finding of gain-of-function effects associated with older age at diagnosis of lung disease compared with loss-of-function effects (P = 0.038). Variants located in the T-BOX and nuclear localization domains were associated with earlier presentation (P = 0.005) and increased incidence of interstitial lung disease (P = 0.003). Event-free survival (death or transplantation) was shorter in the T-BOX group (P = 0.022), although age had a significant effect in the hazard model (P = 0.0461). Carriers of TBX4 variants were diagnosed at a younger age (P < 0.001) and had worse baseline lung function (FEV1, FVC) (P = 0.009) than the BMPR2 and no identified causal variant groups. Conclusions: We demonstrated that TBX4 syndrome is not strictly the result of haploinsufficiency but can also be caused by gain of function. The pleiotropic effects of TBX4 in lung disease may be in part explained by the differential effect of pathogenic mutations located in critical protein domains.