Association of a vascular endothelial growth factor polymorphism with the development of bronchopulmonary dysplasia in Japanese premature newborns.

Our objective was to correlate vascular endothelial growth factor (VEGF) genetic polymorphisms with the risk of bronchopulmonary dysplasia (BPD) development in premature newborns. Fifty-five newborns with BPD (BPD: median gestational age [GA]: 27 weeks, birthweight [BW]: 786 g) and 42 newborns without BPD (non-BPD: median GA: 29 weeks, BW: 1,165 g), who were born at <32 weeks gestational age and were admitted to Kobe University Hospital, were included. BPD was defined as oxygen dependency at 36 weeks postmenstrual age. Genomic DNA was extracted from the umbilical cord, cord blood, or buccal mucosa. Six VEGF genotypes (-1498T > C, -1154G > A, -634C > G, -7C > T, 936C > T, and 1612G > A) were determined by DNA sequencing. Clinical characteristics, and allele and genotype frequencies of VEGF in the BPD and non-BPD groups were analyzed. G allele frequencies in -634C > G of the BPD group were significantly higher than in the non-BPD group (66.4% vs. 50%, P = 0.02). -634C > G genotype distributions differed significantly between the BPD and non-BPD groups (BPD: CC 7%/CG 53%/GG 40%; non-BPD: CC 24%/CG 52%/GG 24%; P = 0.04). Multivariate logistic regression showed that duration of ventilation, VEGF-634G > C G alleles, and male gender were independent risk factors for BPD. In conclusion, polymorphism VEGF -634C > G may influence the risk of BPD.

The incidence of bronchopulmonary dysplasia (BPD), a form of chronic lung disease, is approximately 30% for infants with birth weights < 1,000 g1, and is associated with poor neurodevelopmental and medical outcomes2. Moreover, the absolute numbers of extremely premature infants surviving with BPD is rising because of increasing numbers of preterm births2.

BPD is characterized by impaired alveolarization and vascularization in developing lungs3, and its pathogenesis is multi-factorial, including fetal infection or inflammation, absence of antenatal steroids, oxidant stress, ventilator-induced lung injury, postnatal inflammation/infection, poor nutrition, abnormal growth-factor signaling, and genetic factors45.

Among the neonatal morbidities, significant genetic susceptibility for BPD in preterm infants has previously been identified in two separate multicenter retrospective twin studies67. In the study comparing the intrapair occurrence of BPD within monozygotic and dizygotic twins, sharing 100% and approximately 50% of their genome, respectively, Bhandari et al. found that the observed concordance for developing BPD was significantly higher than the expected concordance in monozygotic twins6. Similarly, Lavoie et al. reported that intrapair correlation coefficients in monozygotic pairs were greater than in dizygotic pairs for the co-occurrence of 36-week oxygen need-based BPD7. According to these studies, the heritability of moderate-severe BPD is estimated at 53–79%8.

Vascular endothelial growth factor (VEGF), a major mediator of vascular permeability, endothelial cell proliferation, and migration, is important in vasculogenesis and angiogenesis910. It also plays a crucial role in fetal development. For example, mice deficient in VEGF showed embryonic lethality because of defects in the development of the vascular system1112. In humans, VEGF mRNA can be detected in fetal tissues from 16 weeks of gestation13. During the fetal period, VEGF is expressed in villous and extravillous trophoblasts14, and the expression level alters with adverse pregnancy outcomes15. Previously, a developmental increase in VEGF levels was reported in umbilical cord plasma, with premature fetuses having significantly lower mean VEGF concentrations than term fetuses16.

VEGF expression levels are particularly high in the lung, where it is essential in lung development and structural maintenance17. A previous immunohistochemical study of aborted human fetuses showed that VEGF was strongly expressed in the epithelial cells of most distal airways, suggesting that it is involved in the maturation as well as proliferation of capillary endothelial cells18. Interestingly, abnormal VEGF signaling was associated with the pathophysiology of pulmonary hypertension1920, which is strongly linked to BPD development2122. Infants with BPD demonstrate decreased VEGF levels in their lungs, suggesting an association between disrupted alveolar vascular development and impaired VEGF expression in BPD patients23.

The VEGF gene is located on chromosome 6q21.3, and consists of eight exons and seven introns. Several VEGF polymorphisms have been identified that influence VEGF protein production through the regulation of VEGF expression24252627. For example, some are associated with neonatal disease morbidities, especially in retinopathy of prematurity2829. However, no studies have investigated the correlation between VEGF polymorphisms and BPD risk in the Japanese population, so the present study set out to achieve this aim.

Results

Clinical characteristics of patients in BPD and non-BPD groups

Univariate and multivariate analyses were performed of the 97 patients. The median GA of preterm infants in the BPD group was significantly younger than in the non-BPD group (27 [range, 23–31] vs. 29 [range, 25–31] weeks). The median BWs of the preterm infants also differed significantly between the groups (786 [range, 430–1,626] vs. 1,165 [range, 800–1,810] g). Furthermore, the BPD group had a significantly lower Apgar score at 5 min, and a significantly higher duration of ventilation and oxygen than the non-BPD group. Univariate analysis of categorized factors revealed that gender, maternal smoking status, SGA, respiratory distress syndrome, and use of surfactant were significant risk factors for the development of BPD. The risk of BPD was not significantly influenced by other tested factors (Table 1).

Table 1

Patient clinical characteristics of BPD and non-BPD groups

	BPD, n = 55	non-BPD, n = 42	P-value
Gestational age (weeks)	27 (23–31)	29 (25–31)	<0.01
Birth weight (g)	786 (430–1626)	1165 (800–1810)	<0.01
Apgar score at 5 min	8 (1–10)	9 (3–10)	<0.01
Male (%)	35/55 (63.6)	17/42 (40.5)	<0.05
Multiple pregnancy (%)	15/55 (27.3)	10/42 (23.8)	0.70
Cesarean section (%)	52/55 (94.6)	37/42 (88.1)	0.25
Antenatal steroids (%)	22/55 (40)	21/41 (51.2)	0.27
Maternal age (years)	32 (19–42)	33 (21–43)	0.13
Primiparity (%)	25/55 (45.5)	12/42 (28.6)	0.09
Histological chorioamnionitis (%)	18/54 (33.3)	12/38 (31.6)	0.86
Maternal smoking (%)	6/54 (11.1)	11/41 (26.8)	<0.05
Assisted reproductive technology (%)	12/54 (22.2)	11/42 (26.2)	0.65
Pregnancy-induced hypertension (%)	13/55 (23.6)	9/42 (21.4)	0.80
Premature rupture of the membranes (%)	13/55 (23.6)	7/42 (16.7)	0.40
Small-for-gestational age (%)	31/55 (56.4)	12/42 (28.6)	<0.01
Symptomatic patent ductus arteriosus (%)	17/48 (35.4)	9/42 (21.4)	0.14
Respiratory distress syndrome (%)	51/55 (92.7)	32/42 (76.2)	<0.05
Use of surfactant (%)	52/55 (94.5)	30/42 (71.4)	<0.01
Duration of ventilation (days)	61 (2–129)	25 (0–60)	<0.01
Duration of oxygen (days)	75 (33–143)	33 (7–67)	<0.01
Culture-proven sepsis (%)	5/55 (9.1)	2/42 (4.8)	0.41
Serum IgM at birth ≥ 20 mg/dL (%)	1/54 (1.9)	2/41 (4.9)	0.40

Data are expressed as medians (range) or numbers (%).

VEGF allele and genotype frequencies of the BPD and non-BPD groups

The observed genotype frequencies did not deviate from Hardy-Weinberg equilibrium. The VEGF polymorphisms -1498 T > C, -1154 G > A, -634 C > G, -7 C > T, 936 C > T, and 1612 G > A in the BPD and non-BPD groups are shown in Table 2. The distribution of allele (BPD: C 33.6%/G 66.4%; non-BPD: C 50%/G 50%; P = 0.02) and genotype (BPD: CC 7%/CG 53%/GG 40%; non-BPD: CC 24%/CG 52%/GG 24%; P = 0.04) frequencies of the VEGF -634C > G polymorphism differed significantly between the BPD and non-BPD groups (Table 2).

Table 2

VEGF allele and genotype frequencies of BPD and non-BPD groups

Polymorphism	Genotype	BPD, n = 55	non-BPD, n = 42	P-value
-1498T > C	TT/TC/CC	24/23/8	22/14/6	0.66
-1154G > A	GG/GA/AA	38/15/2	25/16/1	0.51
-634C > G	CC/CG/GG	4/29/22	10/22/10	0.04
-7C > T	CC/CT/TT	34/20/1	33/9/0	0.17
936C > T	CC/CT/TT	27/27/1	26/15/1	0.42
1612G > A	GG/GA/AA	40/15/0	33/9/0	0.51

VEGF haplotypes

Linkage disequilibrium (LD) plots from the analysis of VEGF SNPs are shown in Figure 1. Strong LD was observed between the -634C > G (rs2010963) polymorphism and the three SNPs in the promoter and 5′ UTR region. Haploview analysis based on the four-gamete rule44 identified one block (<1 kb) containing -634C > G (rs2010963) and -7C > T (rs25648) in the 5′ UTR. Table 3 shows the two-locus [rs2010963 and rs25648] VEGF haplotype. As shown, the -634G/-7C haplotype was significantly decreased in patients with BPD compared with those without (P = 0.02).

Figure 1

Haploview linkage disequilibrium (LD) plots of the six VEGF SNPs analyzed.

Red blocks, D′ (normalized linkage disequilibrium measure or D) ≤ 1.0, with logarithm of odds (LOD) score ≥ 2.0; white blocks, D′ < 1.0 with LOD < 2.0; blue blocks, D′ = 1.0 with LOD < 2.0. Numbers in blocks denote D′ values. The genomic organization is described above the LD plot. LOD was defined as log10(L1/L0), where L1 = likelihood of the data under linkage disequilibrium, and L0 = likelihood of the data under linkage equilibrium. D′ was calculated as follows: D′ = (D) divided by the theoretical maximum for the observed allele frequencies.

Table 3

VEGF haplotype frequencies of BPD and non-BPD groups

	Frequency	BPD, n = 110	non-BPD, n = 84	χ2	P-value
-634G/-7C	0.433	51 (46.4%)	33 (39.3%)	0.972	0.32
-634C/-7C	0.407	37 (33.6%)	42 (50.0%)	5.283	0.02
-634G/-7T	0.16	22 (20.0%)	9 (10.7%)	3.059	0.08

Multivariate stepwise logistic regression analysis of factors associated with BPD development

Stepwise logistic regression analysis was performed to determine whether G alleles in VEGF-634C > G were independently associated with the development of BPD based on factors that showed a significant difference between groups, including gestational age, birth weight, Apgar score at 5 min, gender, maternal smoking status, SGA, respiratory distress syndrome, use of surfactant, and duration of ventilation. Based on the finding by Awata et al. that healthy Japanese subjects without VEGF-634C > G G alleles showed significantly higher serum VEGF levels than those with G alleles35, we assumed a dominant model in this analysis. Because BPD was defined in terms of oxygen requirement, the duration of oxygen supplementation was not included as an independent variable in this analysis45. This revealed that duration of ventilation (P < 0.01), the carrying of G alleles in VEGF-634C > G (P = 0.03), and male gender (P = 0.02) were significant independent risk factors for the development of BPD (Table 4).

Table 4

Multivariate stepwise logistic regression analysis of BPD risk factors

	Likelihood ratio chi2	R2	P-value
Duration of ventilation	43.57	0.34	<0.01
Any G -634C > G	4.53	0.37	0.03
Male gender	5.20	0.41	0.02

Discussion

Because BPD has a multifactorial etiology5, it is important to collect major causative factors to explore the effects of genetic factors on its development6. We collected details of antenatal and postnatal causative factors from our population that was treated uniformly at a single institution, and used multivariate analysis to show that VEGF-634C > G, the duration of mechanical ventilation, and male gender were independent risk factors. Mechanical ventilation4647 and male gender4849 have previously been described as risk factors of BPD; however, the causative effect of VEGF-634C > G on the development of BPD was elucidated for the first time in this study.

Recently, an imbalance of pro-inflammatory and anti-inflammatory cytokines leading to the activation of cellular death pathways in the lung, followed by healing and repair was proposed to be a pathologic hallmark of BPD. Interestingly, each of these steps is closely associated with genetic susceptibility, and BPD is considered to develop from gene-environment interactions, not as secondary to environmental factors650. Additionally, many growth factors under genetic regulation were found to be associated with the disturbance of late lung development51. Thus, we believe that a genetic predisposition is intriguing in consideration of the development of BPD.

Lung angiogenesis is largely regulated by VEGF52, and the importance of VEGF in the development of BPD has been reported in several previous studies. Bhatt et al. carried out a lung comparison study between infants dying with BPD or nonpulmonary causes, and found that the former group had decreased VEGF mRNA levels and decreased VEGF immunostaining compared with the latter group23. Another study investigated VEGF in tracheal aspirates, and revealed that preterm infants who developed BPD had lower VEGF levels during the early postnatal days than those without the disease, suggesting a prolonged and more severe respiratory distress53. Meanwhile, Levesque et al. reported that low VEGF protein levels in infant urine and the need for mechanical ventilation were both associated with BPD development; however, the VEGF protein association was no longer significant once the model was adjusted for respiratory support54. Therefore, the correlation between BPD development and systemic VEGF levels remains unclear. Nevertheless, we speculate that pathogenic lung angiogenesis in BPD is mostly driven by local VEGF synthesis.

Various VEGF genetic polymorphisms have been identified that influence the levels of VEGF expression2427. Our results showed that there was a significant difference in distribution of the -634C > G allele and genotype in premature Japanese infants with BPD. Additionally, one haploblock that included the -634C > G allele was identified, while the frequency of haplotype -634C/-7C that contains only C alleles showed a significant decrease in the BPD compared with the non-BPD population. The -634 C > G polymorphism is located in the 5′ UTR of the VEGF gene, and is predicted to lie within a potential myeloid zinc finger protein (MZF1) binding site, which controls the binding specificity of this motif24. Adult carriers of the -634 C > G polymorphism CC genotype were found to have significantly higher levels of plasma VEGF35, while Lambrechts et al. showed that the G allele reduces internal ribosome entry site-mediated VEGF expression and translation of the large L-VEGF isoform55. However, no data are available on the genetic effects on VEGF levels in immature and developing lungs. We predict that Japanese premature newborns with the VEGF -634 G allele will show decreased VEGF levels in their lungs and will be susceptible to BPD.

Multiple testing procedures usually increase the probability of obtaining false-positives unless statistical adjustments such as Sidak or Bonferroni corrections are employed. However, such corrections are based on the hypothesis that each testing is completely independent, and would remarkably overcorrect for the false-positive rate if each test correlated, resulting in a reduction in power56. Indeed, the types of correction to be adapted for a genetic polymorphism study remain controversial. In the current study, we observed remarkable LD in our six VEGF SNPs studied, so we attempted to assess the statistical significance of haplotype frequency differences by performing permutation tests using Haploview instead of conservative correction, following previous reports575859. We observed a single χ2 value greater than 5.283 392 times (empirical P = 0.0392) from 10,000 permutation tests.

In contrast with our results and others, Kwinta et al. found no association between the allele state of VEGF -634 C > G and risk of BPD in a Polish population60. The first genome-wide association study (GWAS) for BPD that included only Caucasian and African ethnic infants identified SPOCK2 SNPs as genetic risks, but not VEGF -634 C > G61. Moreover, the recent GWAS for BPD could not identify SNPs associated with BPD at a genome-wide significance level8. However, we cannot compare our results with these directly because of the study differences with respect to ethnic background, since ethnicity greatly influences the distribution of gene polymorphisms. We were interested in studying the VEGF -634 C > G polymorphism in the Japanese population because it has a homogeneous ethnic origin compared with the more heterogeneous characteristics of ethnic groups examined in previous studies860.

One limitation of the present study is its relatively small sample size as it was a retrospective study conducted at a single center. Although this sample size is much smaller than those of previous studies861, the first GWAS did not include an Asian population61 while the more recent GWAS included a population of 167 Asian/Pacific islanders that was only 9.6% of the total cohort8. We therefore believe that our results based on 97 infants of single Japanese ethnicity that were uniformly treated in a single institution are meaningful considering the genetic risks of BPD in Asian populations. Further studies of larger sample sizes are therefore necessary to clarify the observed VEGF genotype-specific effects on BPD development.

In conclusion, we identified VEGF -634 C > G as an independent risk factor for the development of BPD using a single center study of Japanese infants less than 32 weeks' gestational age; those infants with the -634 G allele were shown to have a higher susceptibility for BPD.

Methods

Study design and patient groups

This retrospective study was conducted under the approval of the ethical committee of Kobe University Graduate School of Medicine with written informed consent acquired from the parents of the patients. This study was also performed in accordance with the ethics guidelines for human genome/gene analysis research from the Ministries of Education, Culture, Sports, Science and Technology, Health, Labour and Welfare, and Economy, Trade and Industry in Japan.

Ninety-seven infants who were born at our center between 2002 and 2011 with a gestational age (GA) less than 32 weeks were included in the study. All infants were treated at the neonatal intensive care unit of Kobe University Hospital. Patients with congenital or chromosomal anomalies were excluded. Study patients were classified into newborns with BPD (BPD group, n = 55) and newborns without BPD (non-BPD group, n = 42), according to their level of oxygen dependency at a postnatal gestational age of 36 weeks. Clinical, laboratory, and VEGF genetic polymorphism data between BPD and non-BPD groups were analyzed.

Patient clinical data

Clinical data were collected from the patient's record and included GA, birth weight (BW), Apgar score (a numerical expression of an infant's condition after birth, usually determined from heart rate, respiratory efforts, muscle tone, reflex irritability, and color) at 5 min30, gender, antenatal steroid uses, maternal age, primiparity, chorioamnionitis histology, maternal smoking status, the presence or absence of assisted reproductive technology, pregnancy-induced hypertension (maternal systolic blood pressure > 140 mmHg and/or diastolic pressure > 90 mmHg during pregnancy), premature rupture of the membranes (rupture of the membranes more than 24 h before delivery), small for gestational age (SGA; BW less than the 10th percentile of mean BW at the same gestational age in Japanese newborns)31, symptomatic patent ductus arteriosus (therapeutic use of indomethacin or surgical ligation), respiratory distress syndrome diagnosed by chest X-ray findings and microbubble test, use of surfactant, duration of ventilation (regardless of the mode of ventilation), duration of oxygen, culture-proven sepsis and increased serum IgM concentration at birth (≥20 mg/dL), and details of multiple pregnancies and birth by Cesarean section28. BPD was diagnosed when the infants required oxygen therapy at a corrected postnatal gestational age of 36 weeks in association with characteristic radiologic findings32.

VEGF polymorphism genotyping

Genomic DNA was extracted from the umbilical cord, cord blood, or buccal mucosa using a DNeasy Tissue Kit (Qiagen, Hilden, Germany) or QIAamp DNA Mini Kit (Qiagen), according to the manufacturer's instructions. Six VEGF candidate polymorphisms were selected: -1498T > C (rs833061) and -1154G > A (rs1570360) in the promoter region, -634C > G (rs2010963) and -7C > T (rs25648) in the 5′ UTR, and 936C > T (rs3025039) and 1612G > A (rs10434) in the 3′ UTR. Among the polymorphisms that showed a minor allele frequency of at least 0.10 in the Japanese population, we chose these six single nucleotide polymorphisms (SNPs) because they have been shown to influence gene expression2733, protein production2425, and plasma levels26, and were also previously investigated in terms of their association with disease or clinical situations in Japanese populations3435363738, according to previous reports2839. They were identified by PCR, and restriction fragment length polymorphism analysis or direct sequencing28.

Statistical analysis

Data are expressed as the median (range) or number (%). Univariate analyses were performed using the Mann-Whitney nonparametric rank test, chi-square test, or Fisher's exact test as appropriate to compare data between the two groups. To prevent the possible effects of confounding factors, significant factors identified in univariate analysis underwent further analyses, including stepwise logistic regression. Differences were deemed statistically significant when P < 0.05. Analyses were performed using GraphPad Prism 5.0 software (Graphpad Software, Inc., San Diego, CA, USA) and JMP version 8.0.2 (SAS Institute Japan., Tokyo, Japan).

To test whether the genotypes were in Hardy-Weinberg equilibrium, we compared the observed genotype frequencies with expected frequencies at equilibrium based on the chi-square test, according to previous reports4041. LD analysis and haplotype reconstruction were performed using Haploview 4.2 (http://www.broadinstitute.org/scientific-community/science/programs/medical-and-population-genetics/haploview/haploview)42. Sample size calculations were performed using G*Power 3.1 (http://www.gpower.hhu.de/)43, based on our previous report that showed the frequency of common alleles of the six VEGF SNPs to range between 60 and 90% and the prevalence of BPD to be approximately 50% in the studied population28. If a common allele frequency was calculated to differ by 20% between BPD and non-BPD groups, this would result in an 80% chance of detecting this difference with a type I risk error of 5% with 50 infants (100 alleles) in each group.

Author Contributions

K.F. and I.M. designed the study, carried out the statistical analyses and wrote the manuscript. All authors contributed to the intellectual content of this manuscript. K.F., A.S., T.Y., T.K., M.N., K.I. and I.M., were pediatricians who treated the patients and were collected the clinical data and samples. H.Y. were obstetricians and collected the cord blood samples of the patients. K.F., M.Y. and Y.T. analyzed the VEGF polymorphism genotyping.