# Association of a FGFR-4 Gene Polymorphism with Bronchopulmonary Dysplasia and Neonatal Respiratory Distress

Authors: Milad Rezvani; Juliane Wilde; Patricia Vitt; Beena Mailaparambil; Ruth Grychtol; Marcus Krueger; Andrea Heinzmann
Journal: Disease Markers
 (2013)

---

## Abstract

Background. Bronchopulmonary dysplasia (BPD) is the most common chronic lung disease of premature birth, characterized by impaired alveolar development and inflammation. Pathomechanisms contributing to BPD are poorly understood. However, it is assumed that genetic factors predispose to BPD and other pulmonary diseases of preterm neonates, such as neonatal respiratory distress syndrome (RDS). For association studies, genes upregulated during alveolarization are major candidates for genetic analysis, for example, matrix metalloproteinases (MMPs) and fibroblast growth factors (FGFs) and their receptors (FGFR). Objective. Determining genetic risk variants in a Caucasian population of premature neonates with BPD and RDS. Methods. We genotyped 27 polymorphisms within 14 candidate genes via restriction fragment length polymorphism (RFLP): MMP-1, -2, -9, and -12, -16, FGF receptors 2 and 4, FGF-2, -3, -4, -7, and -18, Signal-Regulatory Proteinα (SIRPA) and Thyroid Transcription Factor-1 (TTF-1). Results. Five single nucleotide polymorphisms (SNPs) in MMP-9, MMP-12, FGFR-4, FGF-3, and FGF-7 are associated (P<0.05) with RDS, defined as surfactant application within the first 24 hours after birth. One of them, in FGFR-4 (rs1966265), is associated with both RDS (P=0.003) and BPD (P=0.023). Conclusion. rs1966265 in FGF receptor 4 is a possible genetic key variant in alveolar diseases of preterm newborns.


---

## Body

## 1. Introduction

Great progress has been made in neonatal care over the last few decades, reflected by improving survival rates and clinical outcomes of preterm infants. Despite these advances, 45 years after its first description [1], bronchopulmonary dysplasia (BPD) remains a major complication of premature birth, causing ongoing morbidity and mortality: it is the most common neonatal chronic lung disease, affecting around 25% to 35% of VLBW neonates (very low birth weight, <1500 g) [2], and is associated with increased risk for rehospitalization [3, 4], cognitive delay, and neurosensory deficits [5].Initially described by Northway et al. in 1967, the “old BPD” mainly affected modestly premature newborns suffering from respiratory distress and therefore mechanically ventilated with high levels of supplemental oxygen [1]. With the introduction of surfactant treatment, prenatal maternal use of glucocorticoids, improved nutrition, and ventilator strategies the clinical course and pathology of BPD have changed considerably. Unlike the original description, today’s “new BPD” is mainly regarded as a disruption of distal lung growth [6, 7]. The underlying etiology is multifactorial. Thus, influenced by both genetic susceptibility [8, 9] and environmental factors [7] on the immature lung, the pathophysiology is characterized by inflammation, abnormal microvascularization, and impaired alveolarization [10].Alveolar formation of the primitive saccules is a complex process of epithelial morphogenesis, capillary growth, and coordinated extracellular matrix (ECM) remodelling. At this, fibroblast growth factor (FGF) signalling and matrix metalloproteinase (MMP) activity play eminent roles.MMPs are zinc-dependent proteolytic enzymes degrading all forms of ECM-components [11]. Some MMPs are upregulated in inflammatory environment and yet are involved in pulmonary host defense [11]. There is evidence for some MMP isoforms being important determinants for alveolarization, especially MMP-2 [12], -9 [13], and -16 [14]: MMP-2 deficient mice showed fewer and larger alveoli with thinner interstitial tissue [12]. Hadchouel et al. demonstrated an increase of MMP-16 activity during the alveolar stage and moreover found two SNPs within the MMP-16 gene being associated with lower tracheal MMP-2 and -16 activity and to protect from BPD [14]. Prospecting further potential biomarkers for BPD, also MMP-9 shows some promise; for example, Harijith et al. highlighted a MMP-9-dependent lung injury pathway in an IFNγ-mediated animal model of BPD. Mice with a partial MMP-9 deficiency showed a reversal of IFNγ-induced lung injury during hyperoxia [13].MMPs, particularly MMP-2 and -9, activate fibroblast growth factors (FGFs) by cleavage in the ECM, especially during angiogenesis [15]. In turn activated FGFs upregulate MMP expression [15]. FGFs are secreted glycoproteins involved in interactions between epithelium and mesenchyme regulating cell migration and proliferation in embryonic development, especially in fetal pulmogenesis [16, 17]. Their signalling depends on membrane-located receptors (FGFRs) with a tyrosine kinase domain, encoded by four different genes (FGFR 1–4) [18–20]. They are all translated in developing lungs and are suggested to play major roles in modifying distal lung patterns during alveolarization [21]; for example, FGFR-3-FGFR-4 double-knockout mice show no alveolarization [22].It has been assumed that heritable determinants contribute significantly to both BPD [8, 23] and RDS [24]. On this account, we were interested in identifying genetic risk factors in a Caucasian population of premature newborn with BPD and RDS. We genotyped 27 polymorphisms within fourteen candidate genes for BPD: MMP-1, -2, -9, -12, and -16, FGF receptors 2 and 4, FGF-2, -3, -4, -7, and -18, signal-regulatory proteinα (SIRPA), and thyroid transcription factor-1 (TTF-1). We also included SIRPA because of the known effect on surfactant proteins [25] and inhibition of macrophages [26], as well as TTF-1 due to its effect on lung differentiation [27].


## 2. Material and Methods

### 2.1. Subjects

We recruited preterm neonates (≤28 weeks of gestation) born between January 1996 and September 2010 at the Centre for Pediatrics and Adolescent Medicine, University Hospital Freiburg, Germany. To provide an ethnically homogenous population, all infants were of Caucasian origin. Twins and siblings were excluded from the study as were children with chromosomal aberrations, congenital heart defects, or other major congenital malformations. DNA was collected by buccal swabs or by routine blood sampling, between 2 weeks up to 2 years of age.Medical charts of all recruited infants were reviewed and clinical data were recorded. This included gestational week, number of days with supplemental oxygen, need of mechanical ventilation and positive airway pressure, and need of surfactant therapy. As described previously [28], the subdivision of our BPD study population was based on the analysis by Lavoie et al. about the heritability of BPD [9] according to the consensus defined by the National Institute of Health [6]: the BPD population included all infants with moderate and severe BPD, that is, supplemental oxygen for at least 28 days plus need of oxygen and/or positive pressure at 36 weeks of gestation, whereas the control population consisted of all preterm neonates with no or mild BPD.Recruiting neonates for the RDS population was targeted on severe cases of respiratory distress by including only newborns depending on surfactant within the first 24 hours after birth (see Supplementary Material available online athttp://dx.doi.org/10.1155/2013/932356). At our Neonatal Intensive Care Unit (NICU) the following approach has been applied regarding the treatment with surfactant: avoiding of intubation independent of the gestational week. Therefore, even very premature infants are only intubated if they show failure of ventilation and/or need of supplemental oxygen above 40%. Once they required intubation during the immediate postnatal period, they receive surfactant within 2 hours. This practice was consistent during the whole study period.


### 2.2. Genotyping

DNA was extracted by standard procedures as previously described [28]. Genotyping was carried out by RFLP following polymerase chain reaction amplification (PCR). Studied SNPs are summarized in Table 1. We included a minority of polymorphisms that already had been tested for other pathologies (rs1799750 in MMP1 [29] and rs2276109 and rs652438 in MMP12[30]).Table 1
Studied polymorphisms. Genotyped polymorphisms, their corresponding reference SNP ID number, nucleotide substitution, and gene region. SNPs in exon regions are indicated as their amino acid substitution.
GenePolymorphismAllelesGene regionMMP-1rs1799750G/GGPromotorrs7125062C/TIntronMMP-2rs7201A/C3′UTRrs243865C/TPromotorrs17301608C/TIntronMMP-9rs20544C/T3′UTRrs3918242C/TPromoterrs17576A/G279 Gln/ArgMMP-12rs2276109A/GPromoterrs652438A/G357 Asn/SerMMP-16rs2664352C/TIntronFGFR-2rs2981579A/GIntronrs1219648A/GIntronFGFR-4rs376618A/G136 Leu/Prors1966265A/G10 Ile/ValFGF-2rs11938826C/GIntronFGF-3rs11263591G/CIntronrs10796856T/CIntronFGF-4rs3740639C/T3′UTRFGF-7rs16962440C/TIntronFGF-7, C15 or f33rs4316697G/AIntronrs10519230C/GIntronFGF-18rs4559013G/AIntronSIRPArs3828016A/GIntronrs4814734A/TIntronTTF-1rs11628131A/CPromoterrs999460C/TDownstreamFor PCR reactions, genomic DNA was initially denatured at 94°C for 5 minutes and underwent 35–40 cycles of denaturation (94°C for 30 seconds), annealing (1 minute, corresponding temperatures displayed in Table2), extension reaction (72°C for 1 minute), and a final extension step at 72°C for 8 minutes. In Table 2 primers, annealing temperatures, and restriction enzymes (RE) are shown. Some primers contain intended single nucleotide mismatches (mutagenic primers) to create sites for restriction enzymes.Table 2
PCR conditions: primers, annealing temperatures, number of cycles and restriction enzymes (RE) for RFLP analysis.
GenePolymorphismForward-primerReverse-primerPCRREMMP-1rs71250625′-ATC TCT GTC GGC AAA TTC GT-3′5′-AGC AAC AAG AAG GAG CTG GA-3′56.5°C, 35 cycl.BseG1MMP-2rs72015′-TGC CCA AGA ATA GAT GCT GA-3′5′-CTG TGA GCC ACA GAA GGT TG-3′59.1°C, 35 cycl.MboIrs2438655′-CTG ACC CCC AGT CCT ATC TG-3′5′-GAG CTG AGA CCT GAA GAG CCA AA-3′59.1°C, 40 cycl.BstXIrs173016085′-GCA GTG TAG ACA CTA GAG GAA GGA T-3′5′-GTC ACC AAG GCC TGT ACC AC-3′59.1°C, 40 cycl.DpnIIMMP-9rs205445′-AGT GCC CTG AGG ACT AGA G-3′5′-CAT ACA TGC ATA CAT ACG TGC ATA C-3′54.2°C, 40 cycl.SacIrs39182425′-GCC TGG CAC ATA GTA GGC CC-3′5′-CTT CCT AGC CAG CCG GCA TC-3′67°C, 40 cycl.XceIrs175765′-AAT TCA CCC TCC CGC ACT CT-3′5′-GTT TTG GGG GCC AAT ACA TGA-3′56.5°C, 35 cycl.SmaIMMP-16rs26643525′-GGT AGG CCT GCT TTG GTT CT-3′5′-GCC TTC CAA ACT CCA CTC ACA T-3′56.5°C, 35 cycl.NlaIIIFGFR-2rs29815795′-CCA CTC AAG TCC CTG CTT GT-3′5′-GTG GGG ACT GAT CAG AGG AA-3′56.5°C, 40 cycl.AccIrs12196485′-GCA CGC CTA TTT TAC TTG ACA CGC-3′5′-CAT TGT GGT GAT CCT TCA CG-3′59.1°C, 40 cycl.HhaIFGFR-4rs3766185′-GGG GGT AAC TGT GCC TAT TCTA-3′5′-TAG CAG GGA GTG AAG GGA TG-3′56.5°C, 40 cycl.XbaIrs19662655′-CAA AGG TGC ACG TGT AGC AG-3′5′-TCC CAC CTC AGA AGC CAT AC-3′56.5°C, 40 cycl.BamHIFGF-2rs119388265′-GAT GAT TTC TCA GGG CCA AA-3′5′-GCT TCA CAT GCC TAT TTG CTT-3′59.1°C, 40 cycl.AluIFGF-3rs112635915′-GGT GTG TGG GCT GAA GAG AT-3′5′-CAC CCT GAC ACT CCT GGT CT-3′61.8°C, 40 cycl.TaqIrs107968565′-CTT CCC CTC CTA GGC TTC AG-3′5′-GAA TGC CAG GTT AGC AAA GG-3′59.1°C, 40 cycl.HhaIFGF-4rs37406395′-AAA AAC ACA CCC GCA GAA CT-3′5′-TGA AGG TCA CCC ACT TCC TC-3′59.1°C, 40 cycl.HinfIFGF-7rs169624405′-TGT CGA ACA CAG TGG TAC CTG-3′5′-GCA AAA GGA AGA TGA GAA AAG AT-3′56.5°C, 40 cycl.EcoRVFGF-7,  C15 or f33rs43166975′-TGG CTC TTT CCA GCA ATT AGA-3′5′-CCC TTT CTC TCC CAG GTT CT-3′59.1°C, 40 cycl.TaiIrs105192305′-AAT GTG ACC GCA TAT TGG TG-3′5′-CAT TTT GGG AAC TGG GGT AA-3′56.5°C, 35 cycl.TaqIFGF-18rs45590135′-GAG GCA CAG AAT GAG GAA GG-3′5′-TAC CTG GGC TCT CAA CAT CC-3′56.5°C, 35 cycl.Alw26ISIRPArs38280165′-TCA AAC CCA GGC TTT CTG AC-3′5′-AAC CAA GGC TGG CAC ATA AG-3′63°C, 40 cycl.NlaIIIrs48147345′-CTG TGT ACA TCC TCT TAA ATG GAC-3′5′-TGC TGT GTG CCA GAG TAT GC-3′59.1°C, 40 cycl.HinfITTF-1rs116281315′-TGC AGC CTC ATT CAG CTT TA-3′5′-GGT TGG AGA GGA GGA AAG GA-3′56.5°C, 40 cycl.TaqIrs9994605′-ACC TGG GCG AGA AGT GAG TA-3′5′-AGA GAG GGG AGA AGG AAT GG-3′56.5°C, 40 cycl.SspI


### 2.3. Sequencing

Accuracy of the RFLP was confirmed by sequencing via dideoxy chain termination method [31], respectively, three controls (homozygous wildtype, heterozygous, and homozygous mutation) for each polymorphism using the Big Dye Terminator cycle sequencing kit on an ABI 310 sequencer (Applied Biosystems). The RFLP analyses included these control references.


### 2.4. Statistical Analysis

Genotyping data of our case-control populations were analysed by using Armitage’s trend test (ATT) for possible association with BPD and RDS as specified previously [28]. Moreover ATT was used to calculate Hardy Weinberg equilibrium (HWE) for each polymorphism.


### 2.5. Approval

The collection of blood/buccal swabs and the experimental procedures were approved by the Ethical Committee of the University of Freiburg. Parents were given written and verbal information about the study and a statement of informed consent was signed by the parents of all enrolled children.


## 2.1. Subjects

We recruited preterm neonates (≤28 weeks of gestation) born between January 1996 and September 2010 at the Centre for Pediatrics and Adolescent Medicine, University Hospital Freiburg, Germany. To provide an ethnically homogenous population, all infants were of Caucasian origin. Twins and siblings were excluded from the study as were children with chromosomal aberrations, congenital heart defects, or other major congenital malformations. DNA was collected by buccal swabs or by routine blood sampling, between 2 weeks up to 2 years of age.Medical charts of all recruited infants were reviewed and clinical data were recorded. This included gestational week, number of days with supplemental oxygen, need of mechanical ventilation and positive airway pressure, and need of surfactant therapy. As described previously [28], the subdivision of our BPD study population was based on the analysis by Lavoie et al. about the heritability of BPD [9] according to the consensus defined by the National Institute of Health [6]: the BPD population included all infants with moderate and severe BPD, that is, supplemental oxygen for at least 28 days plus need of oxygen and/or positive pressure at 36 weeks of gestation, whereas the control population consisted of all preterm neonates with no or mild BPD.Recruiting neonates for the RDS population was targeted on severe cases of respiratory distress by including only newborns depending on surfactant within the first 24 hours after birth (see Supplementary Material available online athttp://dx.doi.org/10.1155/2013/932356). At our Neonatal Intensive Care Unit (NICU) the following approach has been applied regarding the treatment with surfactant: avoiding of intubation independent of the gestational week. Therefore, even very premature infants are only intubated if they show failure of ventilation and/or need of supplemental oxygen above 40%. Once they required intubation during the immediate postnatal period, they receive surfactant within 2 hours. This practice was consistent during the whole study period.


## 2.2. Genotyping

DNA was extracted by standard procedures as previously described [28]. Genotyping was carried out by RFLP following polymerase chain reaction amplification (PCR). Studied SNPs are summarized in Table 1. We included a minority of polymorphisms that already had been tested for other pathologies (rs1799750 in MMP1 [29] and rs2276109 and rs652438 in MMP12[30]).Table 1
Studied polymorphisms. Genotyped polymorphisms, their corresponding reference SNP ID number, nucleotide substitution, and gene region. SNPs in exon regions are indicated as their amino acid substitution.
GenePolymorphismAllelesGene regionMMP-1rs1799750G/GGPromotorrs7125062C/TIntronMMP-2rs7201A/C3′UTRrs243865C/TPromotorrs17301608C/TIntronMMP-9rs20544C/T3′UTRrs3918242C/TPromoterrs17576A/G279 Gln/ArgMMP-12rs2276109A/GPromoterrs652438A/G357 Asn/SerMMP-16rs2664352C/TIntronFGFR-2rs2981579A/GIntronrs1219648A/GIntronFGFR-4rs376618A/G136 Leu/Prors1966265A/G10 Ile/ValFGF-2rs11938826C/GIntronFGF-3rs11263591G/CIntronrs10796856T/CIntronFGF-4rs3740639C/T3′UTRFGF-7rs16962440C/TIntronFGF-7, C15 or f33rs4316697G/AIntronrs10519230C/GIntronFGF-18rs4559013G/AIntronSIRPArs3828016A/GIntronrs4814734A/TIntronTTF-1rs11628131A/CPromoterrs999460C/TDownstreamFor PCR reactions, genomic DNA was initially denatured at 94°C for 5 minutes and underwent 35–40 cycles of denaturation (94°C for 30 seconds), annealing (1 minute, corresponding temperatures displayed in Table2), extension reaction (72°C for 1 minute), and a final extension step at 72°C for 8 minutes. In Table 2 primers, annealing temperatures, and restriction enzymes (RE) are shown. Some primers contain intended single nucleotide mismatches (mutagenic primers) to create sites for restriction enzymes.Table 2
PCR conditions: primers, annealing temperatures, number of cycles and restriction enzymes (RE) for RFLP analysis.
GenePolymorphismForward-primerReverse-primerPCRREMMP-1rs71250625′-ATC TCT GTC GGC AAA TTC GT-3′5′-AGC AAC AAG AAG GAG CTG GA-3′56.5°C, 35 cycl.BseG1MMP-2rs72015′-TGC CCA AGA ATA GAT GCT GA-3′5′-CTG TGA GCC ACA GAA GGT TG-3′59.1°C, 35 cycl.MboIrs2438655′-CTG ACC CCC AGT CCT ATC TG-3′5′-GAG CTG AGA CCT GAA GAG CCA AA-3′59.1°C, 40 cycl.BstXIrs173016085′-GCA GTG TAG ACA CTA GAG GAA GGA T-3′5′-GTC ACC AAG GCC TGT ACC AC-3′59.1°C, 40 cycl.DpnIIMMP-9rs205445′-AGT GCC CTG AGG ACT AGA G-3′5′-CAT ACA TGC ATA CAT ACG TGC ATA C-3′54.2°C, 40 cycl.SacIrs39182425′-GCC TGG CAC ATA GTA GGC CC-3′5′-CTT CCT AGC CAG CCG GCA TC-3′67°C, 40 cycl.XceIrs175765′-AAT TCA CCC TCC CGC ACT CT-3′5′-GTT TTG GGG GCC AAT ACA TGA-3′56.5°C, 35 cycl.SmaIMMP-16rs26643525′-GGT AGG CCT GCT TTG GTT CT-3′5′-GCC TTC CAA ACT CCA CTC ACA T-3′56.5°C, 35 cycl.NlaIIIFGFR-2rs29815795′-CCA CTC AAG TCC CTG CTT GT-3′5′-GTG GGG ACT GAT CAG AGG AA-3′56.5°C, 40 cycl.AccIrs12196485′-GCA CGC CTA TTT TAC TTG ACA CGC-3′5′-CAT TGT GGT GAT CCT TCA CG-3′59.1°C, 40 cycl.HhaIFGFR-4rs3766185′-GGG GGT AAC TGT GCC TAT TCTA-3′5′-TAG CAG GGA GTG AAG GGA TG-3′56.5°C, 40 cycl.XbaIrs19662655′-CAA AGG TGC ACG TGT AGC AG-3′5′-TCC CAC CTC AGA AGC CAT AC-3′56.5°C, 40 cycl.BamHIFGF-2rs119388265′-GAT GAT TTC TCA GGG CCA AA-3′5′-GCT TCA CAT GCC TAT TTG CTT-3′59.1°C, 40 cycl.AluIFGF-3rs112635915′-GGT GTG TGG GCT GAA GAG AT-3′5′-CAC CCT GAC ACT CCT GGT CT-3′61.8°C, 40 cycl.TaqIrs107968565′-CTT CCC CTC CTA GGC TTC AG-3′5′-GAA TGC CAG GTT AGC AAA GG-3′59.1°C, 40 cycl.HhaIFGF-4rs37406395′-AAA AAC ACA CCC GCA GAA CT-3′5′-TGA AGG TCA CCC ACT TCC TC-3′59.1°C, 40 cycl.HinfIFGF-7rs169624405′-TGT CGA ACA CAG TGG TAC CTG-3′5′-GCA AAA GGA AGA TGA GAA AAG AT-3′56.5°C, 40 cycl.EcoRVFGF-7,  C15 or f33rs43166975′-TGG CTC TTT CCA GCA ATT AGA-3′5′-CCC TTT CTC TCC CAG GTT CT-3′59.1°C, 40 cycl.TaiIrs105192305′-AAT GTG ACC GCA TAT TGG TG-3′5′-CAT TTT GGG AAC TGG GGT AA-3′56.5°C, 35 cycl.TaqIFGF-18rs45590135′-GAG GCA CAG AAT GAG GAA GG-3′5′-TAC CTG GGC TCT CAA CAT CC-3′56.5°C, 35 cycl.Alw26ISIRPArs38280165′-TCA AAC CCA GGC TTT CTG AC-3′5′-AAC CAA GGC TGG CAC ATA AG-3′63°C, 40 cycl.NlaIIIrs48147345′-CTG TGT ACA TCC TCT TAA ATG GAC-3′5′-TGC TGT GTG CCA GAG TAT GC-3′59.1°C, 40 cycl.HinfITTF-1rs116281315′-TGC AGC CTC ATT CAG CTT TA-3′5′-GGT TGG AGA GGA GGA AAG GA-3′56.5°C, 40 cycl.TaqIrs9994605′-ACC TGG GCG AGA AGT GAG TA-3′5′-AGA GAG GGG AGA AGG AAT GG-3′56.5°C, 40 cycl.SspI


## 2.3. Sequencing

Accuracy of the RFLP was confirmed by sequencing via dideoxy chain termination method [31], respectively, three controls (homozygous wildtype, heterozygous, and homozygous mutation) for each polymorphism using the Big Dye Terminator cycle sequencing kit on an ABI 310 sequencer (Applied Biosystems). The RFLP analyses included these control references.


## 2.4. Statistical Analysis

Genotyping data of our case-control populations were analysed by using Armitage’s trend test (ATT) for possible association with BPD and RDS as specified previously [28]. Moreover ATT was used to calculate Hardy Weinberg equilibrium (HWE) for each polymorphism.


## 2.5. Approval

The collection of blood/buccal swabs and the experimental procedures were approved by the Ethical Committee of the University of Freiburg. Parents were given written and verbal information about the study and a statement of informed consent was signed by the parents of all enrolled children.


## 3. Results

The results of the 27 studied polymorphisms (Table1) for association with bronchopulmonary dysplasia and neonatal respiratory distress are specified in Table 3 (BPD) and Table 4 (RDS).Table 3
Results for BPD. Genotype distribution as numbers of individuals for each genotype. P values are given in accordance with Hardy-Weinberg equilibrium (HWE) and for association with BPD by Armitage’s trend test.
GenePolymorphismGenotypesBPDHWEBPDGenotypescontrolsHWEcontrolsP values for associationMMP-1rs179975016, 25, 140,59034, 44, 390.0090.532rs712506241, 15, 30,36474, 43, 110.2410.126MMP-2rs720121, 28, 91,00043, 64, 210.8600.747rs24386535, 18, 60,17080, 35, 120.0130.675rs1730160827, 25, 110,28051, 63, 210.8570.769MMP-9rs2054417, 19, 170,05327, 60, 370.8560.508rs391824235, 12, 40,07776, 36, 41.0000.894rs1757624, 18, 120,04257, 46, 220.0320.639MMP-12rs227610944, 13, 20,32886, 38, 20.5220.574rs65243852, 8, 01,00010, 19, 10.5870.615MMP-16rs266435212, 31, 191,00040, 57, 350.1200.185FGFR-2rs298157921, 32, 100,79951, 58, 250.2810.858rs121964819, 33, 120,80452, 61, 200.8540.215FGFR-4rs37661835, 23, 50,74484, 40, 100.1320.432rs196626547, 15, 10,47379, 45, 91.0000.023FGF-2rs1193882648, 14, 11,00096, 38, 10.3050.557FGF-3rs1126359125, 22, 150,03745, 61, 290.3820.710rs1079685617, 26, 220,13537, 59, 400.1250.632FGF-4rs374063933, 23, 70,36574, 46, 150.0860.815FGF-7rs1696244017, 35, 120,61442, 70, 250.7300.656FGF-7, C15 or f33rs431669728, 30, 40,37765, 55, 150.5500.868FGF-7, C15 or f33rs1051923037, 23, 31,00073, 50, 110.5170.421FGF-18rs455901322, 23, 170,04545, 62, 221.0000.387SIRPArs382801620, 31, 121,00048, 59, 270.2880.788rs481473416, 37, 100,20646, 62, 250.5980.553TTF-1rs1162813120, 27, 160,31154, 50, 290.0120.287rs99946034, 19, 110,01960, 51, 240.0400.414Table 4
Results for RDS. Genotype distribution as numbers of individuals for each genotype. P values are given for Hardy-Weinberg equilibrium (HWE) and for association with RDS.
GenePolymorphismGenotypes RDSHWERDSGenotypes controlsHWE controlsP value for associationMMP-1rs179975025, 44, 201,00024, 29, 330.0040.163rs712506251, 36, 100.00366, 21, 90.3300.076MMP-2rs720130, 49, 160.67834, 43, 160.0670.648rs24386558, 26, 110.01258, 27, 70.1390.540rs1730160842, 47, 180,42237, 43, 140.8270.838MMP-9rs2054429, 34, 240,82814, 44, 300.0520.033rs391824255, 23, 60,17857, 25, 21.0000.421rs1757638, 30, 210,00545, 34, 130.1530.164MMP-12rs227610967, 27, 21,00064, 23, 11.0000.595rs65243889, 10, 01,00073, 17, 11.0000.047MMP-16rs266435226, 46, 270,54626, 43, 270.3120.998FGFR-2rs298157933, 40, 200,84040, 52, 150.2890.351rs121964842, 49, 160,83730, 47, 160.8330.347FGFR-4rs37661871, 27, 90,01851, 36, 61.0000.286rs196626576, 28, 21,00050, 35, 80.6050.003FGF-2rs1193882680, 26, 10,69067, 26, 10.6840.586FGF-3rs1126359131, 48, 250,43840, 37, 190.0850.127rs1079685633, 50, 260,44221, 38, 360.0920.036FGF-4rs374063957, 41, 71,00050, 30, 160.0090.209FGF-7rs1696244030, 57, 210,56630, 49, 170.8360.588FGF-7 C15 or f33rs431669756, 43, 71,00038, 44, 121.0000.044FGF-7, C15 or f33rs1051923061, 41, 50,79051, 33, 90.3010.418FGF-18rs455901332, 45, 260,23635, 43, 131.0000.079SIRPArs382801634, 53, 211,00037, 38, 190.1340.487rs481473431, 58, 180,33432, 43, 170.6740.671TTF-1rs1162813136, 47, 230,32641, 30, 220.0020.453rs99946056, 37, 160,02738, 36, 190.0820.121


### 3.1. Matrix Metalloproteinases

Among the 11 genotyped polymorphisms in differentMMP genes (see Table 1) there was no BPD-associated polymorphism (Table 3) but two polymorphisms associated (P<0.05) with RDS (rs20544 in MMP-9:P=0.033; rs652438 in MMP-12: P=0.047, see Table 4). Both SNPs show no significant deviation from Hardy-Weinberg equilibrium, neither in the control nor in the case population. Analysis of rs20544 (C/T) identifies the T allele as protective against respiratory distress. For the genotyping results of the amino acid substitution rs652438 (A/G, Asn357Ser) the complete absence of the G/G homozygous genotype in the respiratory distress case population must be taken in account. The other MMP-SNPs showed no association, inclusively rs2664352 in MMP16, that had been associated with protection from BPD [14].


### 3.2. Fibroblast Growth Factor Receptors

TheFGFR-4 SNP rs1966265, located in the exon region and causing an amino acid substitution of Isoleucine (Ile) for Valine (Val) is associated with both BPD (P=0.023) and RDS (P=0.003). Here the A/A genotype (Ile) could be identified as protective allele variant against our studied lung diseases. The association results from significant differences in allele frequencies: in both BPD and RDS analysis the G allele is more frequent in the disease populations (see Tables 3 and 4). The other SNPs in theFGFRgenes showed no association with neither BPD nor respiratory distress.


### 3.3. Fibroblast Growth Factors

Whereas no association could be detected between the eightFGF-SNPs and BPD, rs10796856 in FGF-3 and rs4316697 in FGF-7showed associations with RDS. Correspondent P values are P=0.036 (rs10796856) and P=0.044 (rs4316697), and no deviations from Hardy-Weinberg equilibrium were detected (see Table 4).


### 3.4.SIRPA and TTF-1

The four SNPs inSIRPAand TTF-1 showed no association with neither BPD nor RDS. Analysis of TTF-1rs999460 unfolds deviation from Hardy-Weinberg equilibrium in both case and one control populations in our Caucasian population (see Tables 3 and 4).


## 3.1. Matrix Metalloproteinases

Among the 11 genotyped polymorphisms in differentMMP genes (see Table 1) there was no BPD-associated polymorphism (Table 3) but two polymorphisms associated (P<0.05) with RDS (rs20544 in MMP-9:P=0.033; rs652438 in MMP-12: P=0.047, see Table 4). Both SNPs show no significant deviation from Hardy-Weinberg equilibrium, neither in the control nor in the case population. Analysis of rs20544 (C/T) identifies the T allele as protective against respiratory distress. For the genotyping results of the amino acid substitution rs652438 (A/G, Asn357Ser) the complete absence of the G/G homozygous genotype in the respiratory distress case population must be taken in account. The other MMP-SNPs showed no association, inclusively rs2664352 in MMP16, that had been associated with protection from BPD [14].


## 3.2. Fibroblast Growth Factor Receptors

TheFGFR-4 SNP rs1966265, located in the exon region and causing an amino acid substitution of Isoleucine (Ile) for Valine (Val) is associated with both BPD (P=0.023) and RDS (P=0.003). Here the A/A genotype (Ile) could be identified as protective allele variant against our studied lung diseases. The association results from significant differences in allele frequencies: in both BPD and RDS analysis the G allele is more frequent in the disease populations (see Tables 3 and 4). The other SNPs in theFGFRgenes showed no association with neither BPD nor respiratory distress.


## 3.3. Fibroblast Growth Factors

Whereas no association could be detected between the eightFGF-SNPs and BPD, rs10796856 in FGF-3 and rs4316697 in FGF-7showed associations with RDS. Correspondent P values are P=0.036 (rs10796856) and P=0.044 (rs4316697), and no deviations from Hardy-Weinberg equilibrium were detected (see Table 4).


## 3.4.SIRPA and TTF-1

The four SNPs inSIRPAand TTF-1 showed no association with neither BPD nor RDS. Analysis of TTF-1rs999460 unfolds deviation from Hardy-Weinberg equilibrium in both case and one control populations in our Caucasian population (see Tables 3 and 4).


## 4. Discussion

Bronchopulmonary dysplasia and respiratory distress syndrome of preterm infants have complex pathogenic mechanisms. The aim of this study has been to identify genetic risk factors in an ethnically homogenous Caucasian population.Genetic contribution to BPD is suggested on the basis of twin studies demonstrating that at least half of the susceptibility is hereditary [8, 9, 23]. Additionally, Lavoie et al. [9] could differentiate in their study that mild BPD (according to the National Institute of Child Health and Human Development consensus definition [6]) had been mainly attributable to shared environmental factors whereas moderate or severe BPD had been attributable to genetic influence. Following these findings, we defined our control population as neonates with no BPD or mild BPD, whereas our BPD population included neonates with moderate or severe BPD. Furthermore, we recruited only preterm neonates ≤28 weeks of gestational age for the BPD population to avoid false associations based on the fact that BPD hardly develops in newborn older than 30 weeks of gestational age.In contrast to BPD, the results of twin studies on RDS susceptibility showed mostly contradictory results [24, 32–35]. A twin study by Levit et al. with 332 twin pairs of a heterogeneous population has been the first one to include and assess the influence of several independent covariates, revealing that 50% of the variance to RDS susceptibility is hereditary [24].Given these lines of evidence for genetic contribution, we have chosen the candidate-gene approach for our association study based on the hypothesis that genes fundamental in lung organogenesis and alveolar remodelling, that is,MMP and FGF, determine susceptibility to BPD and RDS.Known genetic risk factors for RDS are mostly allelic polymorphisms of the genes encoding surfactant proteins SP-A1, SP-A2, and SP-B [36]. Anyhow, other determinants than components of the surfactant system might also affect the liability to RDS. Genes encoding for growth factors or enzymes that account for alveolarization through proper secondary septation and extracellular remodeling might affect the gas-exchange and therefore aggravate respiratory distress at birth.Supposed genetic risk factors for BPD are mostly genes encoding components of innate immunity and antigen-presentation, cytokines, antioxidant defences, and angiogenic growth factors such as: mannose-binding lectin (MBL2) [37], tumor necrosis factor-alpha (TNF-α) [28, 38], human leucocyte antigen (HLA)-A, -B, and -C alleles [39], glutathione-S-transferase-P1 [40], and vascular endothelial growth factor (VEGF) [28]. Some years ago, two MMP-16 gene polymorphisms were demonstrated to protect from BPD and moreover to be associated with lower tracheal MMP-2 and -16 levels [14].Matrix metalloproteinases are a family of zinc-dependent endopeptidases [11], and they degrade extracellular components and play a crucial role in lung development, especially during alveolarization. Particularly MMP-2 and -9 (so-called gelatinases A and B) seem to be relevant in extracellular remodeling and even pulmonary host defense. They degrade type IV collagen, fibronectin, elastin, and denatured collagen (gelatin). MMP-2 deficient mice show an abnormal saccular development with larger and simplified alveoles [12]. In line with this finding, newborns developing BPD showed low MMP-2 tracheal levels at birth [41, 42]. Recently, MMP-9 could be identified as a pathogenic key mediator in a murine model of BPD [13]. On the other hand, increased tracheal levels of MMP-9 early after birth have been associated with resolving RDS, suggesting that increase in MMP-9-activity is a physiologic repair response [43]. Dik et al. demonstrated that increased MMP-9 activity in neonatal lungs early after birth correlated with resolving respiratory distress syndrome, demonstrating a likely role of MMP-9 in pulmonary host defense [43].In our study we identified an SNP (rs20544) in theMMP-9 gene to be associated (P=0.033) with RDS, but not BPD. Respiratory distress syndrome has been defined as need of surfactant (see Supplementary Material). On one hand, ethnically homogenous populations like our Caucasian population are favourable to detect possible pathogenetic determinants, but one must bear in mind that the size of our RDS population is limited and the total numbers of neonates studied for each polymorphism vary slightly according to the recruiting time point. Furthermore, association studies on RDS are prone to confounding factors. Other pulmonary conditions such as a transient tachypnea provoked by wet lung syndrome or pulmonary infection might mimic respiratory distress syndrome caused by surfactant deficiency and thereby hamper the results of our study.In our study, we includedMMP-16 polymorphisms that had been associated with BPD in a French population (rs2664352) [14]. In our population rs2664352 did not show any association.TheFGF-3 polymorphism rs10796856 showed association with RDS (P=0.036). Up to now, the role of FGF3 has been mainly studied in cancer diseases, that is, lung cancer [44], but its exact role in pulmogenesis remains elusive. FGF-3 is encoded by the fgf3/int2 gene. There is evidence for FGF-3 upregulation to be associated with alveolar type 2 cell hyperplasia [45] and downregulation to be associated with an excessive recruitment of free alveolar macrophages [45] which might lead to symptoms of respiratory distress. Furthermore it has been shown that FGF-3 stimulates the secretion of MMP-2 and -9 propeptides in vitro [46].FGFR-4 polymorphism rs1966265 showed association with both respiratory distress (P=0.003) and bronchopulmonary dysplasia (P=0.023). The A/A genotype (encoding for Isoleucin instead of Valine) has been protective in our association study. The exact-test showed no deviation from Hardy Weinberg equilibrium for this SNP in both case and control populations, suggesting that the association does not result from population admixture or genotyping errors.FGFR-1 to FGFR-4 are expressed in the lung and FGFR-3 and -4 signalling, in particular, appears to be fundamental in alveolar formation. Weinstein et al. demonstrated that mice deficient in both FGFR-3 and -4 show a completely blocked alveolarization and fail to show any formation of secondary septae, whereas solely FGFR-4(−/−) animals exhibit no significant abnormalities, revealing a cooperative effect of FGFR-3 and -4 in lung development [22]. Hyperoxia-exposed (FiO2 0.85) mice show a BPD-like lung pattern of enlarged airspaces and furthermore a reduced expression of FGFR-3 and -4, suggesting a pathogenic role in arrested lung development [47]. Srisuma et al. [48] replicated these results in FGFR-3 and -4 deficient mice and demonstrated in addition that FGFR-3/-4 signaling contributes to excessive elastin production and its alveolar accumulation, which is another typical feature of BPD. But these abnormalities have not been due to fibroblast defects but due to increased expression of paracrine factors of alveolar type 2 cell (AT2) [48].If a reduction inFGFR-3 and -4 expression affects distal lung development, a functionally significant polymorphism within the correspondent gene possibly alters the susceptibility to alveolar disease such as BPD and RDS. Powell et al. showed that there is a peak of FGFR-4 expression at the day of birth, when respiratory distress syndrome occurs [21]. This supports the conclusion that defective FGFR-4 signalling possibly results in neonatal lung diseases.Our associations do not justify general interpretation. False-positive results can only be excluded by replications in other study populations.In conclusion, we describe five SNPs inMMP-9, MMP-12, FGFR-4, FGF-3,andFGF-7that are associated (P<0.05) in our Caucasian population with respiratory distress syndrome of the newborn, defined as surfactant application within the first 24 hours after birth. Among these polymorphisms one polymorphism in FGFR-4 (rs1966265) is additionally associated with bronchopulmonary dysplasia, demonstrating its possible role in the pathogenesis of newborn lung diseases on grounds of pulmonal immaturity.



---

*Source: 932356-2013-10-31.xml*