
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

72348
10.1038/s41598-024-72348-5
Article
Phenotype wide association study links bronchopulmonary dysplasia with eosinophilia in children
Kelchtermans Jelte kelchtermj@chop.edu

123
March Michael E. 2
Hakonarson Hakon 123
McGrath-Morrow Sharon A. 13
1 grid.25879.31 0000 0004 1936 8972 Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA USA
2 https://ror.org/01z7r7q48 grid.239552.a 0000 0001 0680 8770 The Center of Applied Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA USA
3 https://ror.org/01z7r7q48 grid.239552.a 0000 0001 0680 8770 Division of Pulmonary and Sleep Medicine, The Children’s Hospital of Philadelphia, 3615 Civic Center Blvd, Abramson Building, Philadelphia, PA 19104 USA
13 9 2024
13 9 2024
2024
14 213912 2 2024
5 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Bronchopulmonary dysplasia (BPD) is a common complication of preterm birth. Despite this, genetic drivers of BPD are poorly understood. The objective of this study is to better understand the impact of single nucleotide polymorphisms (SNPs) previously associated with BPD by examining associations with other phenotypes. We drew pediatric subjects from the biorepository at the Center for Applied Genomics to identify associations between these SNPs and 2,146 imputed phenotypes. Methylation data, external cohorts, and in silico validation methods were used to corroborate significant associations. We identified 60 SNPs that were previously associated with BPD. We found a significant association between rs3771150 and rs3771171 and mean eosinophil percentage in a European cohort of 6,999 patients and replicated this in external cohorts. Both SNPs were also associated with asthma, COPD and FEV1/FVC ratio. These SNPs displayed associations with methylation probes and were functionally linked to ST2 (IL1RL1) levels in blood and lung tissue. Our findings support a genetic justification for the epidemiological link between BPD and asthma. Given the well-established link between ST2 and type 2 inflammation in asthma, these findings provide a rationale for future studies exploring the role of type 2 inflammation in the pathogenesis of BPD.

Subject terms

Clinical genetics
Immunogenetics
Parker B Francis Foundationhttp://dx.doi.org/10.13039/100000066 National Institute of Environmental Health Sciences P30 ES013508 issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Bronchopulmonary dysplasia (BPD) is a common complication of preterm birth and remains a global problem1,2. BPD is broadly defined as a need for supplemental oxygen support and/or respiratory support around the time that a preterm infant is close to term2. Given the observational nature of this diagnosis, it comes as no surprise that it has offered limited prognostic and mechanistic homogeneity2. In pursuit of a deeper understanding of BPD pathophysiology, twin studies, candidate gene studies, and whole genome association studies have been conducted to identify and validate genetic risk factors for BPD3.

Initial twin studies suggested a substantial genetic influence on BPD prevalence, but recent research contradicts this claim for extremely premature infants3,4. Similarly, candidate gene studies and genome wide association studies (GWAS) have struggled to identify reproducible relationships between variants and BPD3. However, some notable exceptions exist. First, Hadchouel et al. implicated SPOCK2 in BPD in a discovery series and a replication population; they also demonstrated increased SPOCK2 expression in hyperoxia exposed rats5. Secondly, through exome sequencing, Li et al. implicated 258 genes, with significant enrichment in collagen fibril organization, morphogenesis of embryonic epithelium, and regulation of Wnt signaling pathways6. Finally, Blume et al., validated links between genes that regulate immune cell adhesion and BPD7.

While these inconsistent findings may be due to limited sample size and lack of power, we echo the sentiment articulated by Lal et al., that BPD likely comprises distinct sub-phenotypes driven by varied pathophysiological processes and genetic factors8. Furthermore, how genetic factors affect BPD risk is likely highly dependent on the specific conditions during pregnancy, birth, and the postnatal period, implying significant gene environment interactions.

Overall, the heterogeneity in BPD diagnosis and the varying effects of specific genetic variants influenced by environmental factors may considerably complicate the genetic study of BPD. Therefore, the complex web of environmental variables surrounding premature infants may obscure the potential genetic signal contributing to BPD risk. To address this issue, we look beyond the neonatal intensive care unit context. Considering the anticipated impact of BPD-linked genetic variants on lung growth, inflammation, and injury response, we propose that these variants may also manifest in phenotypes beyond BPD, including that of asthma. To test this hypothesis, we completed a phenotype wide association screen in children enrolled in the biobank at the Center for Applied Genomics (CAG), exploring potential links between SNPs previously associated with BPD and other phenotypes.

Results

We found 60 SNPs that were previously associated with bronchopulmonary dysplasia (Supplemental Table I)9–42. Of these, 3 were significantly associated with at least one phenotype in the current study after Bonferroni correction (Table 1). First, the effect alleles of SNP rs3771150 and rs3771171 were found to be significantly associated with higher mean eosinophil percentage (P = 1.06 × 10−5; coefficient 0.03; 95% CI 0.02–0.04 and P = 6.63 × 10−5; coefficient 0.03; 95% CI 0.02–0.04, respectively) amongst European patients (Fig. 1). Second, rs2077079 was significantly associated with higher odds of chronic liver disease in meta-analysis of all ancestries (P = 1.37 × 10−5; OR 1.39; 95% CI 1.20–1.61).Table 1 Significant associations revealed by PheWAS by ancestry.

Ancestry	Phecode	Description	SNP	p-value	CAG records*	Cases	Controls	
European	2000.33	Mean eosinophil percent (Log10 transformed)	rs3771150	1.06E−05	6999			
2000.33	Mean eosinophil percent (Log10 transformed)	rs3771171	6.63E−06	6999			
Meta analysis	571	Chronic liver disease and cirrhosis	rs2077079	1.37E−05	72,313	383	71,930	
Results for phenotypes with sample sizes of 200 cases or more for binary traits and 1000 individuals or more for quantitative traits. P-values from logistic regression model for binary traits and linear regression model for quantitative traits. * ‘CAG records’ refers to number of patients within the CAG biobank with information available for the relevant phenotype.

Fig. 1 Manhattan plot with PheWAS results for rs3771150 and rs3771171 in a European cohort. Significance threshold of P = 2.3 × 10−5 is represented by a dotted horizontal line.

Of these associations, the associations between rs3771150 and rs3771171 and eosinophil percentage were externally validated (Table 2). Our European cohort in whom this association was found included 3736 (53%) males and 3263 (47%) females. For the 6972 patients that had age at time of blood draw available, the mean age per patient at time of blood draw was 10.0 years with a standard deviation of 4.6 years. Using the Ensenbl VEP tool, these variants were categorized as intron variants for IL18R1 and IL18RAP respectively. Using the OpenTargets interface, the effect alleles of rs3771150 and rs3771171, previously associated with BPD, were found to have been strongly associated with higher “Eosinophil percentage of white cells” (P-value 2.4 × 10−195 and 3.1 × 10−198 respectively) in a study combining UKBiobank with Blood Cell Consortium data43. Furthermore, both are strongly associated with higher odds of asthma in the UKBiobank (P-value 6.0 × 10−26 and 1.7 × 10−26 respectively)44. The AtFKP database implements a Bottom-Line meta-analysis of all datasets in the database to identify significant variant-phenotype associations. This tool again replicated the relationship between rs3771150 and rs3771171 and higher eosinophil percentage (P-value 1.2 × 10−232 and 1.03 × 10−234 respectively) as well as the relationship between both and higher odds of asthma (P-value 1.86 × 10−71 and 2.58 × 10−72 respectively). Beyond these relationships, this tool also identified significant associations with higher odds of chronic obstructive pulmonary disease (COPD) (P-value 6.22 × 10−9 and 2.89 × 10−8 respectively) and lower forced expiratory volume 1 s over forced vital capacity (FEV1/FVC) (P-value 3.48 × 10−8 and 6.22 × 10–9 respectively). Using the AtFKP Signal Sifter tool we can identify the specific region on chromosome two that has variants associated with eosinophil percentage, asthma, COPD, and FEV1/FVC Ratio, we can then visualize the Bottom-Line meta-analysis results for these different phenotypes using Locus Zoom (Fig. 2)45.Table 2 External validation of associations of interest.

Data source	Phenotype	rs3771150	rs3771171	
Coefficient/odds ratio	P-value	Coefficient/odds ratio	P-value	
OpenTargets	Eosinophil percentage	0.070	2.4 × 10−195	0.070	3.1 × 10−198	
Odds of asthma	1.080	6.0 × 10−26	1.080	1.7 × 10−26	
AtFKP	Eosinophil percentage	0.069	1.2 × 10−232	0.070	1.03 × 10−234	
Odds of asthma	1.080	1.86 × 10−71	1.080	2.58 × 10−72	
Odds of COPD	1.040	6.22 × 10−9	1.040	2.89 × 10−8	
FEV1/FVC	− 0.014	3.48 × 10−8	− 0.015	6.22 × 10−9	
This table presents the significant associations between specific SNPs (rs3771150 and rs3771171) and various phenotypes, validated across two large databases: Open Targets Genetics and the Association to Function Knowledge Portal (AtFKP).

Fig. 2 LocusZoom plots of the region surrounding rs3771150 (2:103060851_G/A) and rs3771171 (2:102985950_T/C). These plots visualize the Bottom-Line meta-analysis results for associations between variants in the region and eosinophil percentage, asthma, COPD, and FEV1/FVC Ratio. Lead variants per phenotype are annotated. LD Ref Var = Linkage Disequilibrium Reference variants. Note variation in p-value scale between phenotypes.

After Bonferroni correction, both rs3771150 and rs3771171 remained significantly associated with decreased methylation at probes cg03938978, cg05295703, and cg19795292 and increased methylation at probes cg00652772, cg13315345, cg02648057, cg05168012, and cg08023416 (Table 3). These probes are located on chromosome 2 and relate to both IL18RAP and IL1RL2.Table 3 Significant association between rs3771150 and rs3771171 and methylation status at select methylation probes.

Probe	Gene	rs3771150	rs3771171	
Beta	SE*	p-values	Samples	Beta	SE*	p-values	Samples	
cg03938978	IL18RAP	− 1.04	0.06	1.05E−52	868	− 0.98	0.06	1.19E−46	868	
cg00652772	FAM183DP	0.20	0.02	3.43E−22	858	0.21	0.02	1.02E−24	859	
cg13315345	IL1RL2	0.18	0.02	1.24E−18	863	0.19	0.02	4.04E−20	863	
cg02648057	FAM183DP	0.23	0.03	7.03E−14	861	0.24	0.03	4.82E−15	861	
cg05295703	NA	− 0.12	0.02	1.04E−09	864	− 0.13	0.02	2.08E−11	864	
cg08023416	IL1RL2	0.15	0.02	3.54E−09	861	0.16	0.02	3.26E−10	861	
cg05168012	FAM183DP	0.12	0.02	3.80E−09	845	0.11	0.02	4.83E−08	845	
cg19795292	IL18R1;IL1RL1	− 0.11	0.02	1.02E−08	850	− 0.13	0.02	6.13E−11	850	
*SE: Standard error.

Both variants are significant protein QTLs (Table 4). In Blood plasma, both SNPs are QTLs for IL1RL1, IL18R1, and IL1R246. Furthermore, in Haploreg rs3771150 and rs3771171 were found to be associated with IL18RAP expression in whole blood47,48.Table 4 Expression QTL associations for rs3771150 and rs3771171 across tissues.

Tissue	Gene	rs3771150	rs3771171	
Beta/NES	P-value	Beta/NES	P-value	
Blood Plasma	IL1RL1	− 0.566	1.0 × 10−111	− 0.584	2.0 × 10−118	
IL18R1	− 0.170	2.3 × 10−10	− 0.190	2.2 × 10−12	
IL1R2	− 0.147	4.8 × 10−8	− 0.139	2.8 × 10−7	
Whole blood	IL18RAP	N/A	6.0 × 10−22	N/A	7.3 × 10−84	
AC007278.3	N/A	N/A	− 0.280	9.7 × 10−32	
Lung tissue	IL1RL1	− 0.290	3.2 × 10−15	− 0.300	7.9 × 10−16	
AC007278.3	− 0.270	4.3 × 10−17	− 0.270	4.3 × 10−17	
This table presents the quantitative trait loci (QTL) associations of rs3771150 and rs3771171 across tissues. Beta coefficients/normalized effect sizes (Beta/NES) and corresponding P-values are reported for each gene-tissue pair. N/A indicates that the association was not available or not significant for that tissue-gene pair.

In the GTEx database, rs3771150 was found to be an expression QTL for AC007278.3 and IL1RL1 in lung tissue. Rs3771171 is an expression QTL for AC007278.3 in whole blood and for both AC007278.3 and IL1RL1 in lung tissue. Both Rs3771150 and rs3771171 are also splicing QTLs for IL1RL1 in lung tissue (NES 0.33 and 0.34, p-value 3.8 × 10−16 and 1.6 × 10−16, respectively). Since these variants are in high linkage disequilibrium (LD) with an r2 of 0.81 in the HaploReg database, we sought to refine the signal by performing a formal QTL-colocalization analysis for IL1RL1 in lung tissue49. For this, we utilized the CAVIAR high confidence set derived from the GTEx database, leveraging the UCSC Genome Browser50,51. The variant with the highest causal posterior probability (CPP) to affect IL1RL1 expression in lung tissue was rs953934 with a CPP of 0.77. In the GTEx database, rs953934 has a NES of − 0.5 with associated p-value 4.7 × 10−72 for IL1RL1 in lung tissue.

Discussion

In this study, we present a novel and reproducible genetic association between BPD, eosinophilia, and asthma. These findings are consistent with epidemiological data showing increased asthma incidence in children with BPD, as well as a family history of asthma being a risk factor for asthma diagnosis in BPD children52. Specifically, our findings highlight the significance of two intronic SNPs, rs3771150 and rs3771171, within the IL18R1 and IL18RAP genes. The effect allele of these SNPs, initially identified by Floros et al. in an African American cohort, were successfully linked to higher eosinophil percentage in our pediatric cohort, a relationship validated externally32. Moreover, these SNPs were significantly associated with higher odds of asthma in the UKBiobank and with higher odds of asthma and COPD, and lower FEV1/FVC ratio in the AtKFP database. Both variants are known eQTLs for IL1RL1 in lung tissue and in our own data these SNPs were connected with methylation patterns previously associated with protein level changes in IL18R1 and IL1RL153,54. In this context, the high LD between both variants should be noted. Both variants are associated with methylation status in our control data, are expression and splicing QTLs in lung tissue, and have pleiotrophic effects on eosinophil percentage, asthma, COPD, and FEV1/FVC. While QTL colocalization analysis revealed that rs953934 on chromosome two appears to be the most likely variant driving IL1RL1 expression in lung tissue, future investigations using functional assays are required to fully asses if rs953934 is the primary driver between the observed associations. It is intriguing that these findings implicate a single genomic area in pulmonary diseases across the lifespan. This observation appears to support our hypothesis that variants impacting lung inflammation have the potential to impact various pulmonary phenotypes. The identified region on chromosome 2 contains the Il1RL1 gene. This gene encodes an IL-33 receptor known as ST2. ST2 has two transcriptional variants—ST2L and sST2—the first is a transmembrane receptor that activates the NF-KB pathway while the second is a soluble protein that acts as a decoy receptor for IL-33 and decreases ST2L mediated signaling55,56. Aligning with our genetic findings, decreased sST2 levels have been implicated in asthma with type 2 inflammation56,57. These observations have led to ongoing efforts to develop pharmacological interventions targeting IL-33-ST2L signaling for use in asthma58.

Paradoxically, increased as opposed to decreased sST2 expression has been found in patients with asthma and type 2 inflammation56. Similarly, a high as opposed to low sST2 on day of life 14 has been associated with BPD59. Two possible explanations for these observations exist. First, during episodes of acute injury the body may initially rely on a preformed pool of sST2 to prevent runaway inflammation, placing patients with lower baseline sST2 at a disadvantage56. However, this explanation seems less likely, given the presence of a negative feedback loop that significantly induces sST2 secretion60. An alternative explanation that could be explored is that variants within IL1RL1 may shift the balance between sST2 and ST2L, decreasing the efficacy of negative feedback mechanisms. In this context, it is worth noting that both rs3771150 and rs3771171 are sQTLs in lung tissue located between IL1RL1 exon 5 and 661.

We note that we found a significant association with eosinophil percentage only in patients of European ancestry, whereas the initial study implicating these variants in BPD reported an effect limited to their African ancestry cohort. There are several potential explanations for this discrepancy. First, it is possible that this variation could be attributed to random sampling effects. Second, the higher minor allele frequency (MAF) in patients of European ancestry compared to African ancestry (0.28 vs. 0.10 for rs3771150 and 0.29 vs. 0.10 for rs3771171) may have provided us with greater statistical power in the European cohort for the phenotype of interest (higher eosinophil percentage). Finally, variations in LD across populations could contribute to discrepancies in findings among studies.

As outlined in the introduction, the lack of reproducible relationships between genetic variants and BPD is a major challenge facing the field3. By extension, the use of suggestive disease association SNPs is a weakness of the current study, until fully validated. However, validating relationships between genetic variants and BPD is complicated by BPD compromising distinct sub-phenotypes driven by varied pathophysiological processes and genetic sensitivities8. As such, rather than immediately executing a validation study in an inherently heterogeneous BPD phenotype, we designed the current study to inform future validation efforts. Specifically, by assessing if variants implicated in BPD are associated with non-BPD phenotypes in the absence of premature birth, our aim was to better understand the potential physiologic impact of these variants as to identify potential BPD sub-populations. In this light, the association between rs3771150 and rs3771171 and higher eosinophil percentage creates a rational for future studies assessing if a BPD sub-phenotype with more prominent eosinophilic inflammation can be identified. If this is the case, attempts at validating the effect of rs3771150 and rs3771171 in this sub-phenotype may be a powerful approach. Of note, despite the current lack of direct genomic validation, recent mouse models demonstrating a BPD phenotype in mice with IL-33 knockdown support the potential role of IL1RL1 related variants in BPD62. Another limitation is the use of genetic ancestry as a categorical variable as opposed to using quantitative measures of genetic ancestry. This approach was selected to limit the number of confounders included in the PheWAS analysis and by extension maintain power in this pediatric cohort.

Finally, the previously published link between methylation at cg08023416 and particulate matter exposure may hint at environmental interactions48. Given that our findings support that a subset of patients with BPD is predisposed to develop asthma with significant type 2 inflammation, further characterization of the role of ambient air pollution in long term BPD morbidity and early intervention may be important in optimizing long term outcomes.

Conclusion

In conclusion, our study highlights a genetic links between BPD, eosinophilia, and asthma. While the mechanism behind this relationship remains poorly understood, our findings stress the importance of asthma vigilance and early preventative interventions to optimize lung health. The central role of the IL-33-ST2L pathway in pulmonary health across the lifespan opens avenues for future research and therapeutic exploration, pending validation in subsequent studies.

Materials and methods

Single nucleotide polymorphisms selection

SNPs were selected from enriched literature review starting with SNPs previously summarized by Blume et al. and Solaligue et al.7,63 This list was complemented by reviewing all studies in PubMed identified with MeSH terms ‘Bronchopulmonary Dysplasia’ and ‘Polymorphism, Single Nucleotide’. Any SNPs found to be significantly associated with the development of bronchopulmonary dysplasia, either by meeting genome wide significance in GWAS studies or significance by published p-value in candidate gene studies was included in the current study.

Population

Subjects were drawn from The Children's Hospital of Philadelphia biorepository at the Center for Applied Genomics (CAG). The pediatric samples included in this biorepository are linked to subjects’ EMRs. All subjects have consented to both genomic analysis and EMR mining64. This study does not involve human participants and is limited to secondary analysis of existing de-identified data. As such, it does not constitute Human Subjects Research. This study does not involve human participants and is limited to secondary analysis of existing de-identified data. As such, as defined by 45CFR 46 from the United States Department of Health and Human Services, this study does not constitute Human Subjects Research. All methods were carried out in accordance with relevant guidelines and regulations.

Genotype imputation

Genotype data were generated on 31 different genotyping arrays, with 96.8% of the data coming from four major families of arrays (HumanHapMap550/610Q, OmniExpress, OMNI2.5M, and the GSA array). Array versions from revisions of these arrays were merged on SNPs present on all arrays and filtered for genotype missingness (geno 0.1), individual missingness (mind 0.02), and minor allele frequency (MAF ≥ 0.01), in that order, using PLINK v1.965. Data were imputed using the TOPMed v2 reference panel on the TOPMed Imputation Server66–68. Each imputed file set was filtered for imputation quality on a combination of r-squared (R2) and MAF (for SNPs with MAF ≥ 0.05, R2 ≥ 0.3 were kept; for MAF < 0.05, R2 ≥ 0.5 were kept). File sets were merged, and variants present in 95% of samples were retained. Ancestry was assigned based on the results of principal component analysis (PCA). PCA was performed using flashpca on approximately 2.4 million imputed SNPs with MAF > 0.05 that had been pruned for linkage disequilibrium (LD) using PLINK v1.9, leaving 577,224 variants65,69. The first three principle components were plotted, and patients were assigned to either Asian, African/African American, European, or Hispanic/Latino ancestry by smallest distance to reference genotypes from the HapMap consortium70. No data regarding self-reported race or ethnicity was collected for or used in the current study.

Phenotype wide association study (PheWAS)

International Classification of Diseases 9 (ICD-9) codes were obtained from an anonymized extraction of the CHOP diagnosis database that contained subjects that had been recruited into the patient collection of the CAG. Counts of the occurrence of each ICD-9 code for each subject were generated, and the resulting table was converted into the PheWAS phenotype table. Subjects were included in the case group for each PheWAS phenotype if they possessed two or more occurrences of any of the ICD-9 codes that composed the phenotype in question. Subjects were listed as controls for the PheWAS phenotype if they lacked the case-defining ICD-9 codes, as well as ICD-9 codes corresponding to closely related phenotypes. Phenotypes were analyzed in the PheWAS if they were represented by 20 or more cases in the cohort. Quantitative traits were added to the phenotype table where lab test data were available. For measurements reported multiple times, mean values were calculated and used. Each quantitative trait was examined for normal distribution, and where skewing was in evidence, a log transformation was performed. The subject’s sex and age were included as covariates in the analysis, as were the 10 flashpca generated principle components and a variable representing the group in which genotyping array had been imputed. Genotypes were extracted from the imputed data as allele dose information to preserve some information regarding genotype probability, and the allele doses were used as the genotype inputs to the PheWAS. The PheWAS analysis was performed individually on each PCA-defined ancestry, and then a meta-analysis was performed combining all four ancestries. For the association tests, a logistic regression model was used for binary traits and a linear regression model was used for quantitative traits using p-value significance threshold of P = 2.3 × 10−5 after Bonferroni correction. To ensure appropriate statistical power with variable sample sizes, number of cases, case–control ratio, and minor allele frequency we only report sample sizes of 200 cases or more for binary traits and 1000 individuals or more for quantitative traits as calculated by Verma et al71.

Validation

The Ensembl variant effect predictor (VEP) was then used to assess the likely impact of these variants72. Furthermore, associations of interest were validated in other cohorts by querying the publicly available Open Target Genetics and Association to Function Knowledge Portal (AtFKP) databases73,74. The functional significance of SNPs of interest was validated by investigating their relationship to previously published methylation data obtained from CAG participants75.

The effect of the variants on gene expression was assessed by querying the GTEx database and refined using the USCF Genome Browser and the CAVIAR high confidence set50,51. (See acknowledgement section) Specifically, we conducted a colocalization analysis to identify the most likely causal SNP influencing IL1RL1 expression in lung tissue using GTEx v8 data. SNPs within 1MB of the IL1RL1 transcription start site (chr2:102,311,502–102,352,037) were analyzed using the CAVIAR software, which calculates the causal posterior probability (CPP) for each SNP. Given our focus on respiratory phenotypes, lung tissue was selected for testing.

Software used

R-4.4.1 for Windows available from https://cran.r-project.org/bin/windows/base/

CAVIAR track UCSF genome browser as implemented on https://genome.ucsc.edu/cgi-bin/hgTrackUi?db=hg38&c=chrX&g=gtexEqtlHighConf

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-72348-5.

Acknowledgements

Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The data used for the analyses described in this manuscript were obtained from the GTEx Portal on 08/17/23.

Author contributions

J.K. conceptualized and designed the study, completed data analysis and drafted the initial manuscript. M.M. designed the analysis pipeline and critically reviewed and revised the manuscript. H.H. and S.M. supervised study design and data collection, and critically reviewed and revised the manuscript.

Funding

This work was supported by the Parker B. Francis Fellowship Program. The study was funded by an Institute Development Fund and a K-readiness pilot grant from The Children´s Hospital of Philadelphia. Research reported in this publication was supported by the National Institute of Environmental Health Sciences of the National Institutes of Health under grant number P30 ES013508. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Data availability

The data that support the findings of this study are available from the corresponding author, JK, upon reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors jointly supervised this work and contributed equally: Hakon Hakonarson and Sharon A. McGrath-Morrow.
==== Refs
References

1. Siffel C Kistler KD Lewis JFM Sarda SP Global incidence of bronchopulmonary dysplasia among extremely preterm infants: A systematic literature review J. Matern. Fetal Neonatal Med. 2021 34 1721 1731 10.1080/14767058.2019.1646240 31397199
Siffel, C., Kistler, K. D., Lewis, J. F. M. & Sarda, S. P. Global incidence of bronchopulmonary dysplasia among extremely preterm infants: A systematic literature review. J. Matern. Fetal Neonatal Med. 34, 1721–1731. 10.1080/14767058.2019.1646240 (2021).31397199 10.1080/14767058.2019.1646240
2. Thébaud B Bronchopulmonary dysplasia Nat. Rev. Dis. Primers 2019 5 78 10.1038/s41572-019-0127-7 31727986
Thébaud, B. et al. Bronchopulmonary dysplasia. Nat. Rev. Dis. Primers 5, 78. 10.1038/s41572-019-0127-7 (2019).31727986 10.1038/s41572-019-0127-7
3. Yu K-H Li J Snyder M Shaw GM O’Brodovich HM The genetic predisposition to bronchopulmonary dysplasia Curr. Opin. Pediatr. 2016 28 318 323 10.1097/mop.0000000000000344 26963946
Yu, K.-H., Li, J., Snyder, M., Shaw, G. M. & O’Brodovich, H. M. The genetic predisposition to bronchopulmonary dysplasia. Curr. Opin. Pediatr. 28, 318–323. 10.1097/mop.0000000000000344 (2016).26963946 10.1097/mop.0000000000000344
4. Parad RB Role of genetic susceptibility in the development of bronchopulmonary dysplasia J. Pediatr. 2018 203 234 241 10.1016/j.jpeds.2018.07.099 30287068
Parad, R. B. et al. Role of genetic susceptibility in the development of bronchopulmonary dysplasia. J. Pediatr. 203, 234–241. 10.1016/j.jpeds.2018.07.099 (2018).30287068 10.1016/j.jpeds.2018.07.099
5. Hadchouel A Identification of SPOCK2 as a susceptibility gene for bronchopulmonary dysplasia Am. J. Respir. Crit. Care Med. 2011 184 1164 1170 10.1164/rccm.201103-0548OC 21836138
Hadchouel, A. et al. Identification of SPOCK2 as a susceptibility gene for bronchopulmonary dysplasia. Am. J. Respir. Crit. Care Med. 184, 1164–1170. 10.1164/rccm.201103-0548OC (2011).21836138 10.1164/rccm.201103-0548OC
6. Li J Exome sequencing of neonatal blood spots and the identification of genes implicated in bronchopulmonary dysplasia Am. J. Respir. Crit. Care Med. 2015 192 589 596 10.1164/rccm.201501-0168OC 26030808
Li, J. et al. Exome sequencing of neonatal blood spots and the identification of genes implicated in bronchopulmonary dysplasia. Am. J. Respir. Crit. Care Med. 192, 589–596. 10.1164/rccm.201501-0168OC (2015).26030808 10.1164/rccm.201501-0168OC
7. Blume F Verification of immunology-related genetic associations in BPD supports ABCA3 and five other genes Pediatr. Res. 2022 92 190 198 10.1038/s41390-021-01689-y 34465876
Blume, F. et al. Verification of immunology-related genetic associations in BPD supports ABCA3 and five other genes. Pediatr. Res. 92, 190–198. 10.1038/s41390-021-01689-y (2022).34465876 10.1038/s41390-021-01689-y
8. Lal CV Ambalavanan N Genetic predisposition to bronchopulmonary dysplasia Semin. Perinatol. 2015 39 584 591 10.1053/j.semperi.2015.09.004 26471063
Lal, C. V. & Ambalavanan, N. Genetic predisposition to bronchopulmonary dysplasia. Semin. Perinatol. 39, 584–591. 10.1053/j.semperi.2015.09.004 (2015).26471063 10.1053/j.semperi.2015.09.004
9. Hamvas A Exome sequencing identifies gene variants and networks associated with extreme respiratory outcomes following preterm birth BMC Genet. 2018 19 94 10.1186/s12863-018-0679-7 30342483
Hamvas, A. et al. Exome sequencing identifies gene variants and networks associated with extreme respiratory outcomes following preterm birth. BMC Genet. 19, 94. 10.1186/s12863-018-0679-7 (2018).30342483 10.1186/s12863-018-0679-7
10. Huusko JM Polymorphisms of the gene encoding Kit ligand are associated with bronchopulmonary dysplasia Pediatr. Pulmonol. 2015 50 260 270 10.1002/ppul.23018 24610823
Huusko, J. M. et al. Polymorphisms of the gene encoding Kit ligand are associated with bronchopulmonary dysplasia. Pediatr. Pulmonol. 50, 260–270. 10.1002/ppul.23018 (2015).24610823 10.1002/ppul.23018
11. Capoluongo E Mannose-binding lectin polymorphisms and pulmonary outcome in premature neonates: A pilot study Intensive Care Med. 2007 33 1787 1794 10.1007/s00134-007-0793-x 17653692
Capoluongo, E. et al. Mannose-binding lectin polymorphisms and pulmonary outcome in premature neonates: A pilot study. Intensive Care Med. 33, 1787–1794. 10.1007/s00134-007-0793-x (2007).17653692 10.1007/s00134-007-0793-x
12. Zhang S Surfactant protein B gene polymorphisms is associated with risk of bronchopulmonary dysplasia in Chinese Han population Int. J. Clin. Exp. Pathol. 2015 8 2971 2978 26045806
Zhang, S. et al. Surfactant protein B gene polymorphisms is associated with risk of bronchopulmonary dysplasia in Chinese Han population. Int. J. Clin. Exp. Pathol. 8, 2971–2978 (2015).26045806
13. Winters AH Single nucleotide polymorphism in toll-like receptor 6 is associated with a decreased risk for ureaplasma respiratory tract colonization and bronchopulmonary dysplasia in preterm infants Pediatr. Infect. Dis. J. 2013 32 898 904 10.1097/INF.0b013e31828fc693 23518821
Winters, A. H. et al. Single nucleotide polymorphism in toll-like receptor 6 is associated with a decreased risk for ureaplasma respiratory tract colonization and bronchopulmonary dysplasia in preterm infants. Pediatr. Infect. Dis. J. 32, 898–904. 10.1097/INF.0b013e31828fc693 (2013).23518821 10.1097/INF.0b013e31828fc693
14. Karjalainen MK Haataja R Hallman M Haplotype analysis of ABCA3: Association with respiratory distress in very premature infants Ann. Med. 2008 40 56 65 10.1080/07853890701611094 18246475
Karjalainen, M. K., Haataja, R. & Hallman, M. Haplotype analysis of ABCA3: Association with respiratory distress in very premature infants. Ann. Med. 40, 56–65. 10.1080/07853890701611094 (2008).18246475 10.1080/07853890701611094
15. Ambalavanan N Integrated genomic analyses in bronchopulmonary dysplasia J. Pediatr. 2015 166 531 537 10.1016/j.jpeds.2014.09.052 25449221
Ambalavanan, N. et al. Integrated genomic analyses in bronchopulmonary dysplasia. J. Pediatr. 166, 531–537. 10.1016/j.jpeds.2014.09.052 (2015).25449221 10.1016/j.jpeds.2014.09.052
16. Mailaparambil B Genetic and epidemiological risk factors in the development of bronchopulmonary dysplasia Dis. Markers 2010 29 1 9 10.3233/dma-2010-0720 20826912
Mailaparambil, B. et al. Genetic and epidemiological risk factors in the development of bronchopulmonary dysplasia. Dis. Markers 29, 1–9. 10.3233/dma-2010-0720 (2010).20826912 10.3233/dma-2010-0720
17. Poggi C Genetic contributions to the development of complications in preterm newborns PLoS One 2015 10 e0131741 10.1371/journal.pone.0131741 26172140
Poggi, C. et al. Genetic contributions to the development of complications in preterm newborns. PLoS One 10, e0131741. 10.1371/journal.pone.0131741 (2015).26172140 10.1371/journal.pone.0131741
18. Cakmak BC Association between bronchopulmonary dysplasia and MBL2 and IL1-RN polymorphisms Pediatr. Int. 2012 54 863 868 10.1111/j.1442-200X.2012.03714.x 22882323
Cakmak, B. C. et al. Association between bronchopulmonary dysplasia and MBL2 and IL1-RN polymorphisms. Pediatr. Int. 54, 863–868. 10.1111/j.1442-200X.2012.03714.x (2012).22882323 10.1111/j.1442-200X.2012.03714.x
19. Sampath V Antioxidant response genes sequence variants and BPD susceptibility in VLBW infants Pediatr. Res. 2015 77 477 483 10.1038/pr.2014.200 25518008
Sampath, V. et al. Antioxidant response genes sequence variants and BPD susceptibility in VLBW infants. Pediatr. Res. 77, 477–483. 10.1038/pr.2014.200 (2015).25518008 10.1038/pr.2014.200
20. Rezvani M Association of a FGFR-4 gene polymorphism with bronchopulmonary dysplasia and neonatal respiratory distress Dis. Markers 2013 35 633 640 10.1155/2013/932356 24288432
Rezvani, M. et al. Association of a FGFR-4 gene polymorphism with bronchopulmonary dysplasia and neonatal respiratory distress. Dis. Markers 35, 633–640. 10.1155/2013/932356 (2013).24288432 10.1155/2013/932356
21. Fujioka K Association of a vascular endothelial growth factor polymorphism with the development of bronchopulmonary dysplasia in Japanese premature newborns Sci. Rep. 2014 4 4459 10.1038/srep04459 24662923
Fujioka, K. et al. Association of a vascular endothelial growth factor polymorphism with the development of bronchopulmonary dysplasia in Japanese premature newborns. Sci. Rep. 4, 4459. 10.1038/srep04459 (2014).24662923 10.1038/srep04459
22. Pavlovic J Genetic variants of surfactant proteins A, B, C, and D in bronchopulmonary dysplasia Dis. Markers 2006 22 277 291 10.1155/2006/817805 17264398
Pavlovic, J. et al. Genetic variants of surfactant proteins A, B, C, and D in bronchopulmonary dysplasia. Dis. Markers 22, 277–291. 10.1155/2006/817805 (2006).17264398 10.1155/2006/817805
23. Koroglu OA Association of vitamin D receptor gene polymorphisms and bronchopulmonary dysplasia Pediatr. Res. 2014 76 171 176 10.1038/pr.2014.63 24796371
Koroglu, O. A. et al. Association of vitamin D receptor gene polymorphisms and bronchopulmonary dysplasia. Pediatr. Res. 76, 171–176. 10.1038/pr.2014.63 (2014).24796371 10.1038/pr.2014.63
24. Derzbach L Selectin polymorphisms and perinatal morbidity in low-birthweight infants Acta Paediatr. 2006 95 1213 1217 10.1080/08035250600575404 16982492
Derzbach, L. et al. Selectin polymorphisms and perinatal morbidity in low-birthweight infants. Acta Paediatr. 95, 1213–1217. 10.1080/08035250600575404 (2006).16982492 10.1080/08035250600575404
25. Bokodi G Derzbach L Bányász I Tulassay T Vásárhelyi B Association of interferon gamma T+874A and interleukin 12 p40 promoter CTCTAA/GC polymorphism with the need for respiratory support and perinatal complications in low birthweight neonates Arch. Dis. Child. Fetal Neonatal Ed. 2007 92 F25 29 10.1136/adc.2005.086421 16754651
Bokodi, G., Derzbach, L., Bányász, I., Tulassay, T. & Vásárhelyi, B. Association of interferon gamma T+874A and interleukin 12 p40 promoter CTCTAA/GC polymorphism with the need for respiratory support and perinatal complications in low birthweight neonates. Arch. Dis. Child. Fetal Neonatal Ed. 92, F25-29. 10.1136/adc.2005.086421 (2007).16754651 10.1136/adc.2005.086421
26. Hadchouel A Matrix metalloproteinase gene polymorphisms and bronchopulmonary dysplasia: Identification of MMP16 as a new player in lung development PLoS One 2008 3 e3188 10.1371/journal.pone.0003188 18784838
Hadchouel, A. et al. Matrix metalloproteinase gene polymorphisms and bronchopulmonary dysplasia: Identification of MMP16 as a new player in lung development. PLoS One 3, e3188. 10.1371/journal.pone.0003188 (2008).18784838 10.1371/journal.pone.0003188
27. Trittmann JK Arginase I gene single-nucleotide polymorphism is associated with decreased risk of pulmonary hypertension in bronchopulmonary dysplasia Acta Paediatr. 2014 103 e439 443 10.1111/apa.12717 24919409
Trittmann, J. K. et al. Arginase I gene single-nucleotide polymorphism is associated with decreased risk of pulmonary hypertension in bronchopulmonary dysplasia. Acta Paediatr. 103, e439-443. 10.1111/apa.12717 (2014).24919409 10.1111/apa.12717
28. Ali S Functional genetic variation in NFKBIA and susceptibility to childhood asthma, bronchiolitis, and bronchopulmonary dysplasia J. Immunol. 2013 190 3949 3958 10.4049/jimmunol.1201015 23487427
Ali, S. et al. Functional genetic variation in NFKBIA and susceptibility to childhood asthma, bronchiolitis, and bronchopulmonary dysplasia. J. Immunol. 190, 3949–3958. 10.4049/jimmunol.1201015 (2013).23487427 10.4049/jimmunol.1201015
29. Krueger M Heinzmann A Mailaparambil B Härtel C Göpel W Polymorphisms of interleukin 18 in the genetics of preterm birth and bronchopulmonary dysplasia Arch. Dis. Child. Fetal Neonatal Ed. 2011 96 F299 300 10.1136/adc.2009.174862 20971720
Krueger, M., Heinzmann, A., Mailaparambil, B., Härtel, C. & Göpel, W. Polymorphisms of interleukin 18 in the genetics of preterm birth and bronchopulmonary dysplasia. Arch. Dis. Child. Fetal Neonatal Ed. 96, F299-300. 10.1136/adc.2009.174862 (2011).20971720 10.1136/adc.2009.174862
30. Kazzi SN Kim UO Quasney MW Buhimschi I Polymorphism of tumor necrosis factor-alpha and risk and severity of bronchopulmonary dysplasia among very low birth weight infants Pediatrics 2004 114 e243 248 10.1542/peds.114.2.e243 15286263
Kazzi, S. N., Kim, U. O., Quasney, M. W. & Buhimschi, I. Polymorphism of tumor necrosis factor-alpha and risk and severity of bronchopulmonary dysplasia among very low birth weight infants. Pediatrics 114, e243-248. 10.1542/peds.114.2.e243 (2004).15286263 10.1542/peds.114.2.e243
31. Zachaki S GSTP1 and CYP2B6 genetic polymorphisms and the risk of bronchopulmonary dysplasia in preterm neonates Am. J. Perinatol. 2017 34 729 734 10.1055/s-0036-1597994 28081574
Zachaki, S. et al. GSTP1 and CYP2B6 genetic polymorphisms and the risk of bronchopulmonary dysplasia in preterm neonates. Am. J. Perinatol. 34, 729–734. 10.1055/s-0036-1597994 (2017).28081574 10.1055/s-0036-1597994
32. Floros J IL-18R1 and IL-18RAP SNPs may be associated with bronchopulmonary dysplasia in African-American infants Pediatr. Res. 2012 71 107 114 10.1038/pr.2011.14 22289858
Floros, J. et al. IL-18R1 and IL-18RAP SNPs may be associated with bronchopulmonary dysplasia in African-American infants. Pediatr. Res. 71, 107–114. 10.1038/pr.2011.14 (2012).22289858 10.1038/pr.2011.14
33. Chen LL Dual-specificity phosphatase (DUSP) genetic variants predict pulmonary hypertension in patients with bronchopulmonary dysplasia Pediatr. Res. 2020 87 81 87 10.1038/s41390-019-0502-9 31330530
Chen, L. L. et al. Dual-specificity phosphatase (DUSP) genetic variants predict pulmonary hypertension in patients with bronchopulmonary dysplasia. Pediatr. Res. 87, 81–87. 10.1038/s41390-019-0502-9 (2020).31330530 10.1038/s41390-019-0502-9
34. Luo X Identification of genetic susceptibility in preterm newborns with bronchopulmonary dysplasia by whole-exome sequencing: BIVM gene may play a role Eur. J. Pediatr. 2023 182 1707 1718 10.1007/s00431-022-04779-z 36757497
Luo, X. et al. Identification of genetic susceptibility in preterm newborns with bronchopulmonary dysplasia by whole-exome sequencing: BIVM gene may play a role. Eur. J. Pediatr. 182, 1707–1718. 10.1007/s00431-022-04779-z (2023).36757497 10.1007/s00431-022-04779-z
35. Ryckman KK Dagle JM Kelsey K Momany AM Murray JC Genetic associations of surfactant protein D and angiotensin-converting enzyme with lung disease in preterm neonates J. Perinatol. 2012 32 349 355 10.1038/jp.2011.104 21960125
Ryckman, K. K., Dagle, J. M., Kelsey, K., Momany, A. M. & Murray, J. C. Genetic associations of surfactant protein D and angiotensin-converting enzyme with lung disease in preterm neonates. J. Perinatol. 32, 349–355. 10.1038/jp.2011.104 (2012).21960125 10.1038/jp.2011.104
36. Haas DM The impact of glucocorticoid polymorphisms on markers of neonatal respiratory disease after antenatal betamethasone administration Am. J. Obstet. Gynecol. 2013 208 215 e211 216 10.1016/j.ajog.2012.12.031
Haas, D. M. et al. The impact of glucocorticoid polymorphisms on markers of neonatal respiratory disease after antenatal betamethasone administration. Am. J. Obstet. Gynecol. 208(215), e211-216. 10.1016/j.ajog.2012.12.031 (2013).10.1016/j.ajog.2012.12.031
37. Mahlman M Genes encoding vascular endothelial growth factor A (VEGF-A) and VEGF receptor 2 (VEGFR-2) and risk for bronchopulmonary dysplasia Neonatology 2015 108 53 59 10.1159/000381279 25998098
Mahlman, M. et al. Genes encoding vascular endothelial growth factor A (VEGF-A) and VEGF receptor 2 (VEGFR-2) and risk for bronchopulmonary dysplasia. Neonatology 108, 53–59. 10.1159/000381279 (2015).25998098 10.1159/000381279
38. Akat A Bronchopulmonary dysplasia and wnt pathway-associated single nucleotide polymorphisms Pediatr. Res. 2022 92 888 898 10.1038/s41390-021-01851-6 34853430
Akat, A. et al. Bronchopulmonary dysplasia and wnt pathway-associated single nucleotide polymorphisms. Pediatr. Res. 92, 888–898. 10.1038/s41390-021-01851-6 (2022).34853430 10.1038/s41390-021-01851-6
39. Sampath V A TLR5 (g.1174C > T) variant that encodes a stop codon (R392X) is associated with bronchopulmonary dysplasia Pediatr. Pulmonol. 2012 47 460 468 10.1002/ppul.21568 22058078
Sampath, V. et al. A TLR5 (g.1174C > T) variant that encodes a stop codon (R392X) is associated with bronchopulmonary dysplasia. Pediatr. Pulmonol. 47, 460–468. 10.1002/ppul.21568 (2012).22058078 10.1002/ppul.21568
40. Hilgendorff A Association of polymorphisms in the mannose-binding lectin gene and pulmonary morbidity in preterm infants Genes Immun. 2007 8 671 677 10.1038/sj.gene.6364432 17898783
Hilgendorff, A. et al. Association of polymorphisms in the mannose-binding lectin gene and pulmonary morbidity in preterm infants. Genes Immun. 8, 671–677. 10.1038/sj.gene.6364432 (2007).17898783 10.1038/sj.gene.6364432
41. Prencipe G A polymorphism in the macrophage migration inhibitory factor promoter is associated with bronchopulmonary dysplasia Pediatr. Res. 2011 69 142 147 10.1203/PDR.0b013e3182042496 21045753
Prencipe, G. et al. A polymorphism in the macrophage migration inhibitory factor promoter is associated with bronchopulmonary dysplasia. Pediatr. Res. 69, 142–147. 10.1203/PDR.0b013e3182042496 (2011).21045753 10.1203/PDR.0b013e3182042496
42. Kwinta P Genetic risk factors of bronchopulmonary dysplasia Pediatr. Res. 2008 64 682 688 10.1203/PDR.0b013e318184edeb 18614962
Kwinta, P. et al. Genetic risk factors of bronchopulmonary dysplasia. Pediatr. Res. 64, 682–688. 10.1203/PDR.0b013e318184edeb (2008).18614962 10.1203/PDR.0b013e318184edeb
43. Vuckovic D The polygenic and monogenic basis of blood traits and diseases Cell 2020 182 1214 1231 10.1016/j.cell.2020.08.008 32888494
Vuckovic, D. et al. The polygenic and monogenic basis of blood traits and diseases. Cell 182, 1214–1231. 10.1016/j.cell.2020.08.008 (2020).32888494 10.1016/j.cell.2020.08.008
44. Valette K Prioritization of candidate causal genes for asthma in susceptibility loci derived from UK Biobank Commun. Biol. 2021 4 700 10.1038/s42003-021-02227-6 34103634
Valette, K. et al. Prioritization of candidate causal genes for asthma in susceptibility loci derived from UK Biobank. Commun. Biol. 4, 700. 10.1038/s42003-021-02227-6 (2021).34103634 10.1038/s42003-021-02227-6
45. Boughton AP LocusZoom.js: Interactive and embeddable visualization of genetic association study results Bioinformatics (Oxf, Engl) 2021 37 3017 3018 10.1093/bioinformatics/btab186
Boughton, A. P. et al. LocusZoom.js: Interactive and embeddable visualization of genetic association study results. Bioinformatics (Oxf, Engl) 37, 3017–3018. 10.1093/bioinformatics/btab186 (2021).10.1093/bioinformatics/btab186
46. Sun BB Genomic atlas of the human plasma proteome Nature 2018 558 73 79 10.1038/s41586-018-0175-2 29875488
Sun, B. B. et al. Genomic atlas of the human plasma proteome. Nature 558, 73–79. 10.1038/s41586-018-0175-2 (2018).29875488 10.1038/s41586-018-0175-2
47. Fehrmann RS Trans-eQTLs reveal that independent genetic variants associated with a complex phenotype converge on intermediate genes, with a major role for the HLA PLoS Genet. 2011 7 e1002197 10.1371/journal.pgen.1002197 21829388
Fehrmann, R. S. et al. Trans-eQTLs reveal that independent genetic variants associated with a complex phenotype converge on intermediate genes, with a major role for the HLA. PLoS Genet. 7, e1002197. 10.1371/journal.pgen.1002197 (2011).21829388 10.1371/journal.pgen.1002197
48. Westra HJ Systematic identification of trans eQTLs as putative drivers of known disease associations Nat. Genet. 2013 45 1238 1243 10.1038/ng.2756 24013639
Westra, H. J. et al. Systematic identification of trans eQTLs as putative drivers of known disease associations. Nat. Genet. 45, 1238–1243. 10.1038/ng.2756 (2013).24013639 10.1038/ng.2756
49. Ward LD Kellis M HaploReg v4: Systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease Nucleic Acids Res. 2016 44 D877 881 10.1093/nar/gkv1340 26657631
Ward, L. D. & Kellis, M. HaploReg v4: Systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease. Nucleic Acids Res. 44, D877-881. 10.1093/nar/gkv1340 (2016).26657631 10.1093/nar/gkv1340
50. Nassar LR The UCSC genome browser database: 2023 update Nucleic Acids Res. 2023 51 D1188 d1195 10.1093/nar/gkac1072 36420891
Nassar, L. R. et al. The UCSC genome browser database: 2023 update. Nucleic Acids Res. 51, D1188-d1195. 10.1093/nar/gkac1072 (2023).36420891 10.1093/nar/gkac1072
51. Hormozdiari F Kostem E Kang EY Pasaniuc B Eskin E Identifying causal variants at loci with multiple signals of association Genetics 2014 198 497 508 10.1534/genetics.114.167908 25104515
Hormozdiari, F., Kostem, E., Kang, E. Y., Pasaniuc, B. & Eskin, E. Identifying causal variants at loci with multiple signals of association. Genetics 198, 497–508. 10.1534/genetics.114.167908 (2014).25104515 10.1534/genetics.114.167908
52. Sun T Yu HY Yang M Song YF Fu JH Risk of asthma in preterm infants with bronchopulmonary dysplasia: A systematic review and meta-analysis World J. Pediatr. 2023 19 549 556 10.1007/s12519-023-00701-1 36857022
Sun, T., Yu, H. Y., Yang, M., Song, Y. F. & Fu, J. H. Risk of asthma in preterm infants with bronchopulmonary dysplasia: A systematic review and meta-analysis. World J. Pediatr. 19, 549–556. 10.1007/s12519-023-00701-1 (2023).36857022 10.1007/s12519-023-00701-1
53. Hillary RF Multi-method genome- and epigenome-wide studies of inflammatory protein levels in healthy older adults Genome Med. 2020 10.1186/s13073-020-00754-1 32641083
Hillary, R. F. et al. Multi-method genome- and epigenome-wide studies of inflammatory protein levels in healthy older adults. Genome Med.10.1186/s13073-020-00754-1 (2020).32641083 10.1186/s13073-020-00754-1
54. Dijk FN Genetic regulation of <i>IL1RL1</i> methylation and IL1RL1-a protein levels in asthma Eur. Respir. J. 2018 51 1701377 10.1183/13993003.01377-2017 29519908
Dijk, F. N. et al. Genetic regulation of <i>IL1RL1</i> methylation and IL1RL1-a protein levels in asthma. Eur. Respir. J. 51, 1701377. 10.1183/13993003.01377-2017 (2018).29519908 10.1183/13993003.01377-2017
55. Tunc T Predictive value of soluble urokinase plasminogen activator receptor, soluble ST2, and IL-33 in bronchopulmonary dysplasia Pediatr. Res. 2014 75 788 792 10.1038/pr.2014.28 24603291
Tunc, T. et al. Predictive value of soluble urokinase plasminogen activator receptor, soluble ST2, and IL-33 in bronchopulmonary dysplasia. Pediatr. Res. 75, 788–792. 10.1038/pr.2014.28 (2014).24603291 10.1038/pr.2014.28
56. Gordon ED IL1RL1 asthma risk variants regulate airway type 2 inflammation JCI Insight 2016 1 e87871 10.1172/jci.insight.87871 27699235
Gordon, E. D. et al. IL1RL1 asthma risk variants regulate airway type 2 inflammation. JCI Insight 1, e87871. 10.1172/jci.insight.87871 (2016).27699235 10.1172/jci.insight.87871
57. Rojo-Tolosa S Influence of genetics on the response to Omalizumab in patients with severe uncontrolled asthma with an allergic phenotype Int. J. Mol. Sci. 2023 24 7029 10.3390/ijms24087029 37108192
Rojo-Tolosa, S. et al. Influence of genetics on the response to Omalizumab in patients with severe uncontrolled asthma with an allergic phenotype. Int. J. Mol. Sci. 24, 7029. 10.3390/ijms24087029 (2023).37108192 10.3390/ijms24087029
58. Gaurav R Poole JA Interleukin (IL)-33 immunobiology in asthma and airway inflammatory diseases J. Asthma 2022 59 2530 2538 10.1080/02770903.2021.2020815 34928757
Gaurav, R. & Poole, J. A. Interleukin (IL)-33 immunobiology in asthma and airway inflammatory diseases. J. Asthma 59, 2530–2538. 10.1080/02770903.2021.2020815 (2022).34928757 10.1080/02770903.2021.2020815
59. Tunc T Predictive value of soluble urokinase plasminogen activator receptor, soluble ST2, and IL-33 in bronchopulmonary dysplasia Pediatr. Res. 2014 75 788 792 10.1038/pr.2014.28 24603291
Tunc, T. et al. Predictive value of soluble urokinase plasminogen activator receptor, soluble ST2, and IL-33 in bronchopulmonary dysplasia. Pediatr. Res. 75, 788–792. 10.1038/pr.2014.28 (2014).24603291 10.1038/pr.2014.28
60. Cayrol C IL-33, an alarmin of the IL-1 family involved in allergic and non allergic inflammation: Focus on the mechanisms of regulation of its activity Cells 2021 11 107 10.3390/cells11010107 35011670
Cayrol, C. IL-33, an alarmin of the IL-1 family involved in allergic and non allergic inflammation: Focus on the mechanisms of regulation of its activity. Cells 11, 107. 10.3390/cells11010107 (2021).35011670 10.3390/cells11010107
61. Cunningham F Ensembl 2022 Nucleic Acids Res. 2021 50 D988 D995 10.1093/nar/gkab1049
Cunningham, F. et al. Ensembl 2022. Nucleic Acids Res. 50, D988–D995. 10.1093/nar/gkab1049 (2021).10.1093/nar/gkab1049
62. Zhu Y Yao HC Lu HY Hao XB Xu SQ <scp>IL-33-ST2</scp> pathway regulates <scp>AECII</scp> transdifferentiation by targeting alveolar macrophage in a bronchopulmonary dysplasia mouse model J. Cell. Mol. Med. 2023 27 304 308 10.1111/jcmm.17654 36573439
Zhu, Y., Yao, H. C., Lu, H. Y., Hao, X. B. & Xu, S. Q. <scp>IL-33-ST2</scp> pathway regulates <scp>AECII</scp> transdifferentiation by targeting alveolar macrophage in a bronchopulmonary dysplasia mouse model. J. Cell. Mol. Med. 27, 304–308. 10.1111/jcmm.17654 (2023).36573439 10.1111/jcmm.17654
63. Surate Solaligue DE Rodríguez-Castillo JA Ahlbrecht K Morty RE Recent advances in our understanding of the mechanisms of late lung development and bronchopulmonary dysplasia Am. J. Physiol. Lung Cell Mol. Physiol. 2017 313 L1101 l1153 10.1152/ajplung.00343.2017 28971976
Surate Solaligue, D. E., Rodríguez-Castillo, J. A., Ahlbrecht, K. & Morty, R. E. Recent advances in our understanding of the mechanisms of late lung development and bronchopulmonary dysplasia. Am. J. Physiol. Lung Cell Mol. Physiol. 313, L1101-l1153. 10.1152/ajplung.00343.2017 (2017).28971976 10.1152/ajplung.00343.2017
64. Gottesman O The electronic medical records and genomics (eMERGE) network: Past, present, and future Genet. Med. 2013 15 761 771 10.1038/gim.2013.72 23743551
Gottesman, O. et al. The electronic medical records and genomics (eMERGE) network: Past, present, and future. Genet. Med. 15, 761–771. 10.1038/gim.2013.72 (2013).23743551 10.1038/gim.2013.72
65. Chang CC Second-generation PLINK: Rising to the challenge of larger and richer datasets GigaScience 2015 4 7 10.1186/s13742-015-0047-8 25722852
Chang, C. C. et al. Second-generation PLINK: Rising to the challenge of larger and richer datasets. GigaScience 4, 7. 10.1186/s13742-015-0047-8 (2015).25722852 10.1186/s13742-015-0047-8
66. Das S Next-generation genotype imputation service and methods Nat. Genet. 2016 48 1284 1287 10.1038/ng.3656 27571263
Das, S. et al. Next-generation genotype imputation service and methods. Nat. Genet. 48, 1284–1287. 10.1038/ng.3656 (2016).27571263 10.1038/ng.3656
67. Fuchsberger C Abecasis GR Hinds DA minimac2: Faster genotype imputation Bioinformatics (Oxf., Engl.) 2015 31 782 784 10.1093/bioinformatics/btu704
Fuchsberger, C., Abecasis, G. R. & Hinds, D. A. minimac2: Faster genotype imputation. Bioinformatics (Oxf., Engl.) 31, 782–784. 10.1093/bioinformatics/btu704 (2015).10.1093/bioinformatics/btu704
68. Taliun D Sequencing of 53,831 diverse genomes from the NHLBI TOPMed program Nature 2021 590 290 299 10.1038/s41586-021-03205-y 33568819
Taliun, D. et al. Sequencing of 53,831 diverse genomes from the NHLBI TOPMed program. Nature 590, 290–299. 10.1038/s41586-021-03205-y (2021).33568819 10.1038/s41586-021-03205-y
69. Abraham G Inouye M Fast principal component analysis of large-scale genome-wide data PLoS One 2014 9 e93766 10.1371/journal.pone.0093766 24718290
Abraham, G. & Inouye, M. Fast principal component analysis of large-scale genome-wide data. PLoS One 9, e93766. 10.1371/journal.pone.0093766 (2014).24718290 10.1371/journal.pone.0093766
70. Altshuler DM Integrating common and rare genetic variation in diverse human populations Nature 2010 467 52 58 10.1038/nature09298 20811451
Altshuler, D. M. et al. Integrating common and rare genetic variation in diverse human populations. Nature 467, 52–58. 10.1038/nature09298 (2010).20811451 10.1038/nature09298
71. Verma A A simulation study investigating power estimates in phenome-wide association studies BMC Bioinform. 2018 10.1186/s12859-018-2135-0
Verma, A. et al. A simulation study investigating power estimates in phenome-wide association studies. BMC Bioinform.10.1186/s12859-018-2135-0 (2018).10.1186/s12859-018-2135-0
72. McLaren W The ensembl variant effect predictor Genome Biol. 2016 17 122 10.1186/s13059-016-0974-4 27268795
McLaren, W. et al. The ensembl variant effect predictor. Genome Biol. 17, 122. 10.1186/s13059-016-0974-4 (2016).27268795 10.1186/s13059-016-0974-4
73. Ghoussaini M Open targets genetics: Systematic identification of trait-associated genes using large-scale genetics and functional genomics Nucleic Acids Res. 2021 49 D1311 D1320 10.1093/nar/gkaa840 33045747
Ghoussaini, M. et al. Open targets genetics: Systematic identification of trait-associated genes using large-scale genetics and functional genomics. Nucleic Acids Res. 49, D1311–D1320. 10.1093/nar/gkaa840 (2021).33045747 10.1093/nar/gkaa840
74. Association to Function Knowledge Portal, https://a2f.hugeamp.org/
75. Spindola LM Detecting multiple differentially methylated CpG sites and regions related to dimensional psychopathology in youths Clin. Epigenetics 2019 11 146 10.1186/s13148-019-0740-z 31639064
Spindola, L. M. et al. Detecting multiple differentially methylated CpG sites and regions related to dimensional psychopathology in youths. Clin. Epigenetics 11, 146. 10.1186/s13148-019-0740-z (2019).31639064 10.1186/s13148-019-0740-z
