Skip to main content

Identification of novel functional sequence variants in the gene for peptidase inhibitor 3

Abstract

Background

Peptidase inhibitor 3 (PI3) inhibits neutrophil elastase and proteinase-3, and has a potential role in skin and lung diseases as well as in cancer. Genome-wide expression profiling of chorioamniotic membranes revealed decreased expression of PI3 in women with preterm premature rupture of membranes. To elucidate the molecular mechanisms contributing to the decreased expression in amniotic membranes, the PI3 gene was searched for sequence variations and the functional significance of the identified promoter variants was studied.

Methods

Single nucleotide polymorphisms (SNPs) were identified by direct sequencing of PCR products spanning a region from 1,173 bp upstream to 1,266 bp downstream of the translation start site. Fourteen SNPs were genotyped from 112 and nine SNPs from 24 unrelated individuals. Putative transcription factor binding sites as detected by in silico search were verified by electrophoretic mobility shift assay (EMSA) using nuclear extract from Hela and amnion cell nuclear extract. Deviation from Hardy-Weinberg equilibrium (HWE) was tested by χ2 goodness-of-fit test. Haplotypes were estimated using expectation maximization (EM) algorithm.

Results

Twenty-three sequence variations were identified by direct sequencing of polymerase chain reaction (PCR) products covering 2,439 nt of the PI3 gene (-1,173 nt of promoter sequences and all three exons). Analysis of 112 unrelated individuals showed that 20 variants had minor allele frequencies (MAF) ranging from 0.02 to 0.46 representing "true polymorphisms", while three had MAF ≤ 0.01. Eleven variants were in the promoter region; several putative transcription factor binding sites were found at these sites by database searches. Differential binding of transcription factors was demonstrated at two polymorphic sites by electrophoretic mobility shift assays, both in amniotic and HeLa cell nuclear extracts. Differential binding of the transcription factor GATA1 at -689C>G site was confirmed by a supershift.

Conclusion

The promoter sequences of PI3 have a high degree of variability. Functional promoter variants provide a possible mechanism for explaining the differences in PI3 mRNA expression levels in the chorioamniotic membranes, and are also likely to be useful in elucidating the role of PI3 in other diseases.

Peer Review reports

Background

PI3 [Gene ID: 5266] is a member of the 'trappin' gene family [1]. The trappin gene family members are defined by an amino-terminal transglutaminase substrate domain consisting of hexapeptide repeats with the consensus sequence of GQDPVK and a carboxy-terminal four-disulphide bond core. PI3, also known as trappin-2, elafin, elastase specific inhibitor and skin-derived antileukoproteinase (SKALP), is a low-molecular weight, 6 kDa serine protease inhibitor [2], that is capable of inhibiting neutrophil elastase (also known as elastase 2; ELA2; [GeneID: 1991]) and proteinase 3 (PRTN3; [GeneID: 5657]; also known as the Wegener autoantigen, P29). PI3 has been mapped to chromosome 20q12-13.1 [3], and this locus contains 14 genes expressing protease inhibitor domains with homology to whey acidic protein (WAP). Human PI3 gene spans about 11,620 bp and consists of three exons [2, 4]. The gene has multiple transcription start sites and the mRNA has been reported to have an unusually short 5'-UTR (5'-untranslated region) [5].

Initially, PI3 was identified in human epidermis of psoriatic patients [6], and later in bronchial secretions from patients with bronchial carcinoma [7] and chronic obstructive pulmonary disease [2], as well as in epidermal [8] and breast tumors [9]. In addition to its antipeptidase role, PI3 has antimicrobial activity and is a component of the innate immune system to protect epithelial surfaces from infection [1013]. Expression of PI3 can be induced by inflammatory mediators such as tumor necrosis factor (TNF) and interleukin 1 beta (IL1B) [14, 15].

In our previous report we identified PI3 as a down-regulated gene in the chorioamnionitic membranes of patients with preterm premature rupture of membranes (PPROM) [16]. In this study, we investigated the possible molecular mechanisms that control the expression of PI3 by carrying out a detailed analysis of the PI3 gene sequences.

Methods

Genomic DNA isolation

Blood samples were obtained from 112 healthy unrelated African-American individuals after written informed consent. The collection of samples, and their utilization for research purposes, was approved by the Institutional Review Boards of Wayne State University and the National Institute of Child Health and Human Development, NIH. Genomic DNA was extracted from blood samples using QIAGEN® DNA Blood BioRobot® 9604 kit (QIAGEN Inc., Valencia, CA.).

Direct sequencing of PCR products

Genomic DNA was used as a template to generate three overlapping PCR products of 724 bp, 717 bp and 1,328 bp in size extending from 1,173 bp upstream to 1,266 bp downstream of the translation start site of the PI3 gene [GenBank: NT_011362]. Primers are listed in Table 1. All PCRs were carried out in 100-μl volumes containing 1.5 mM of MgCl2, 0.2 mM dNTPs, 0.4 μM of each primer, 3 U of Taq DNA polymerase (Roche Molecular Systems, Inc., Branchburg, NJ) and 100 ng of genomic DNA. A 10 minute initial denaturation at 94°C was followed by 40 cycles consisting of 30 s denaturation at 94°C, 30 s annealing at 50°C to 55°C, and 1 minute extension at 72°C. PCR products were analyzed on 2% agarose gels. PCR products were purified by ultrafiltration (Centricon Centrifugal Filter Devices, Millipore, Bedford, MA), and sequenced by cycle sequencing and dye terminator labeling (ABI® BigDye™ Terminator v1.1 Cycle Sequencing kit, Applied Biosystems, Foster City, CA). Sequencing reactions were purified using gel filtration columns (CENTRI-SEP, Princeton Separation, Adelphia, NJ) and run on 310 or 3700 Genetic Analyzer (Applied Biosystems). Sequences were edited using BioEdit [17]. Fourteen SNPs were genotyped from 112 unrelated individuals and nine SNPs from 24 unrelated individuals (Table 2).

Table 1 Oligonucleotide primers used in the study.
Table 2 Minor allele frequencies of the 23 SNPs detected in the PI3 gene.

In silicosearch for transcription factor binding sites

The sequences in and around the SNP sites in the promoter region were searched for putative transcription factor binding sites using three different computer programs: TESS [18, 19], Alibaba 2.1 [20, 21], and MatInspector [22, 23]. Default parameters were used as search criteria.

Electrophoretic mobility shift assays (EMSA)

Oligonucleotides (Table 1) and their complementary strands were designed and purchased as gel purified (IDT, Coraville, IA). Complementary oligonucleotides were annealed to each other to generate double-stranded probes. EMSAs were performed using commercially available HeLa cell nuclear extracts (Promega, Madison, WI) and nuclear extracts prepared from primary amnion cell cultures as previously described [24] since we had previously demonstrated that PI3 protein was produced by a variety of chorioamniotic membrane cell types with the highest amount produced by the amniotic epithelial cells [16]. Primary amnion cell cultures were established using amniotic membranes obtained from women not in labor at term who underwent elective cesarean deliveries for obstetrical indications. All other reagents were purchased from a commercial source and used according to the manufacturer's protocol (Promega, Madison, WI). The concentration of poly(dI-dC) (Amersham Biosciences Corp., Piscataway, NJ) in the reaction was optimized to 0.05 μg/μl to minimize non-specific binding. The concentrations of components in 10 μl reaction mixtures were as follows: 1X binding buffer [without poly(dI-dC)], 3.75 μg of HeLa or amnion cell nuclear extract, 0.05 μg/μl poly(dI-dC) and 50 fmole of 32P-labeled double-stranded probe (>50,000 cpm). All the components, except 32P-probe, were added to the reaction and incubated for 15 min on ice and 10 min at 20°C, followed by the addition of the 32P-labelled probe, and incubation for 20 min at 20°C. For competition experiments, a 100-fold molar excess of unlabeled double-stranded oligonucleotides was added to the reaction mixture prior to the addition of the labeled probe. For supershift experiments, polyclonal antibodies against AP1 (Cat. No. sc-253X and sc-44X; Santa Cruz Biotechnology, Santa Cruz, CA), and GATA1 (Active Motif, Carlsbad, CA) were used. For AP1, after 20 min of incubation at 20°C with 32P-labelled probe, 400 ng of corresponding antibodies were added to the reaction and incubated for another 15 min at 20°C. For the GATA1 assay, antibodies were added and incubated for 20 min at 20°C before adding the labeled probe [25]. Samples were run on non-denaturing 6% polyacrylamide gels in 0.5X TBE buffer, at 100 V for 80 min. X-ray film (Kodak, Rochester, NY) was exposed to dried gels 2 to 5 h at -80°C depending on signal intensity.

Nomenclature for sequence variants and genes

The variants and nucleotides are described following the guidelines of the Human Genome Variation Society (HGVS) [26]. SNPs are described using the genomic sequence AL049767.12 as a reference and numbered relative to the translation start site. Official gene symbols provided by Human Genome Organization (HUGO) Nomenclature Committee (HGNC) were used [27].

Statistical analyses

Tests for deviations from HWE were performed by using the χ2 goodness-of-fit test. Haplotypes were estimated following expectation maximization (EM) algorithm as implemented in the software Arlequin [28].

Results

SNP genotyping and haplotype construction

When this study was initiated, only two polymorphisms were known to exist in the PI3 gene, neither of them in the promoter region. We identified 23 SNPs (Table 2 and Fig 1) in the PI3 gene sequences of 24 unrelated individuals by direct sequencing of PCR products that spanned the region from 1,173 bp upstream to 1,266 bp downstream of the translation start site. Of the 23 SNPs, nine are in the dbSNP database [29] (Table 2). Eleven SNPs were located in the promoter region, one in exon 1, seven in intron 1, two in exon 2 and two in intron 2. To obtain more reliable allele frequency estimates, a larger sample of 112 unrelated individuals was genotyped for 14 SNPs. Of the 23 SNPs, three (-911G>A, -453T>G and -258A>G) had minor allele frequencies (MAF) of <0.01 and 20 SNPs had a MAF between 0.02 and 0.46 (Table 2). The genotype frequencies of all SNPs were in HWE. Three SNPs (+50C>T, +959A>C and +1123C>T) altered the codons (Table 2). However, only two SNPs (+50C>T and +959A>C) altered amino acid (Table 2). The T allele at +50 altered the 17th amino acid from threonine to methionine and C allele at +959 altered the 34th amino acid from threonine to proline. The 17th amino acid is part of the signal peptide sequence, whereas the 34th amino acid is part of the proprotein.

Figure 1
figure 1

Locations of the 23 SNPs detected in the region from 1,173 bp upstream to 1,266 bp downstream of the translation start site of the PI3 gene. Dark shaded boxes represent the three exons and light shaded boxes represent introns. For more information on the SNPs, see Table 2.

Thirteen of the 23 SNPs were in complete linkage disequilibrium (Table 3). Altogether 16 haplotypes were identified (Table 3). PI3_F was the most common haplotype followed by PI3_H and PI3_K.

Table 3 Haplotypes constructed with 23 SNPs in the PI3 gene.

Effect of SNPs on protein binding

We performed in silico searches for putative transcription factor binding sites at the 11 SNP sites in the promoter region of the PI3 gene. Except for one SNP (-675C>T), all other sites showed potential differential binding for at least one transcription factor (Table 4). In other words, the transcription factor was predicted to bind to one of the alleles, but not the other for these 10 SNPs.

Table 4 Results from in silico searches for putative transcription factor binding sites.

To verify experimentally the differential binding of transcription factors, we conducted electrophoretic mobility shift assays (EMSA) for 10 SNP sites located in the promoter region. Due to the presence of a long stretch of AC-repeats, EMSA was not carried out for the SNP -868C>G. Of the 10 putative sites, six (-1077A>G, -1067A>G, -1063G>A, -960T>Del, -689C>G, -338G>A) showed differential binding by transcription factors in nuclear extracts derived from HeLa cells (not shown), while only two (-1063G>A and -689A>G) showed differential binding using amnion cell nuclear extract (Fig 2). Those SNP sites that did not show differential binding with HeLa cell nuclear extract, also did not show differential binding with amnion cell nuclear extracts. A transcription factor in HeLa cell nuclear extract bound to -960T>Del. There was however, no binding by a transcription factor in amnion cell nuclear extract to this same site. No transcription factor in either HeLa or amnion cell nuclear extract bound to the -258A>G site. Our interest was the transcriptional regulation of the PI3 gene in amnion cells [16]. We, therefore, focused on the two SNPs (-1063G>A and -689C>G, Fig 2) that showed differential binding by transcription factors derived from the amnion cell nuclear extract. For -1063A>G and -689C>G, the banding patterns representing the protein-DNA complexes were similar when using HeLa and amnion cell nuclear extracts, although the band intensities were lower with the latter, probably due to a lower concentration of functional proteins (Fig 2). To determine the specificity of the binding, we used a competition assay (Fig 3). The differential binding that persisted after cross-competition (100-fold) was considered to be due to the SNP. For -1063A>G, one protein-DNA complex persisted after a labeled double-stranded A-probe was competed with double-stranded sequence differing only at the SNP (G instead of A). For -689C>G, two protein-DNA complexes persisted after labeled double-stranded G-probe was competed with a double-stranded sequence differing only at the SNP (C instead of G) (Fig 3).

Figure 2
figure 2

EMSA showing the banding patterns with HeLa and amniotic cell nuclear extracts for -1063A>G and -689C>G sites. The arrows indicate protein-DNA complexes formed when transcription factors bind to their target sites.

Figure 3
figure 3

EMSA showing self- and cross-competition for differential binding for -1063A>G and -689C>G sites with amniotic or HeLa cell nuclear extracts. The arrows indicate the differential binding that consistently persisted after cross-competition.

Our in silico search predicted that AP1 was the transcription factor that would bind differentially at the -1063 SNP. To investigate this we used 100-fold excess of a competitor with the consensus sequence for AP1 binding or the anti-AP1 antibody in the reaction. No change in the banding pattern was observed in the competition assay (Fig 4). Similarly, no supershift with anti-c-jun or anti-c-fos antibody was observed for -1063G>A polymorphism using amnion or HeLa cell nuclear extracts (Fig 4). Since a positive control, using consensus AP1 binding sequence, demonstrated a supershift against anti-c-jun and anti-c-fos antibodies (Fig 4), a failure in the supershift was unlikely to be due to technical problems. We, therefore, concluded that the protein that binds to the A probe at nt -1063 does not contain the AP1 epitope.

Figure 4
figure 4

Results of competition and supershift experiments for the -1063A>G site using HeLa cell nuclear extract. S, shift; SS, supershift. Arrow on the left indicates a protein-DNA complex specific to the transcription factor binding to the A-allele.

For the SNP at nt -689, the transcription factor, whose binding was predicted to change due to the SNP, was GATA1 (Table 5). As shown in Fig 5, a consensus sequence containing the GATA1 binding site was able to compete with the -689G probe (Fig 5A) and a supershift was observed with anti-GATA1 antibody when using amniotic cell nuclear extract (Fig 5B) indicating that GATA1 binds to the G-allele of the -689C>G polymorphism in the promoter region of PI3 gene.

Table 5 Previously reported transcription factor binding sites in the PI3 promoter.
Figure 5
figure 5

EMSA with amniotic cell nuclear extract for the nt -689C>G SNP site. (A) Competition with -689C, -689G and GATA consensus sequences. The arrows indicate protein-DNA complexes which consistently persisted after cross-competition. (B) Supershift experiment. S, shift; SS, supershift.

Discussion

We observed a high degree of polymorphism within the PI3 gene with 23 SNPs detected, 11 of which were located in the promoter region. We found an amino acid substitution, T34P, in the 4th amino acid of the amino terminal-transglutaminase substrate domain, GQDPVK, of PI3. To determine if this SNP has a significant effect on the function of this domain, we searched for the consensus sequence of the transglutaminase substrate in other mammals. A similar sequence domain was identified in seminal vesicle protein I (Semg1) repeats in guinea pig (PROSITE documentation PDOC000282). Semg1 is a clotting protein that serves as the substrate in the formation of the copulatory plug [30]. Covalent clotting of this protein is catalyzed by a transglutaminase and involves the formation of γ-glutamyl-ε-lysine crosslinks. The consensus signature of this consensus sequence was [IVM]-X-G-Q-D-X-V-K-X5- [KN]-G-X3- [STLV]. The 6th amino acid from the left of this pattern, X, corresponds to the 34th amino acid of the PI3 protein, and appears not to be conserved.

Putative binding of several transcription factors in the promoter region of PI3 gene have been previously reported (Table 5) [5, 9, 14, 15, 31]. Except for two, none of these sites were polymorphic in our study. The two consensus sequences of NFKB1 at nt -964 to -956 and nt -340 to -331 identified in previous studies [5, 9] contained the -960T>Del and -338G>A SNP sites, respectively, in our study. Based on our in silico analysis, the binding of transcription factors NFATC2 and AP1 could be altered by the sequence changes at -960T>Del and -338G>A (Table 4), but with EMSA we did not observe any differences between the two alleles using amnion cell nuclear extract. For -960T>Del we did not observe any binding to a transcription factor using amnion cell extract.

We identified 10 SNPs with alleles that were predicted to have different binding sites for one or more transcription factors by in silico searches (Table 4). We tested the predicted differential binding at these sites by EMSA with both HeLa and amnion cell nuclear extracts. HeLa cells are used widely for studying the functionality of promoter polymorphisms. Since PI3 mRNA is down-regulated in chorioamniotic membranes of patients with PPROM [16], we also used a nuclear extract derived from amnion cells. Six of the 10 sites exhibited differential binding to transcription factors with the HeLa extract, in contrast to only two with amnion cell nuclear extract. Since we did not observe supershift at -1063G>A using antibody against AP1 and nuclear extract from either HeLa or amniotic cell line, it is likely that the differential binding was to a transcription factor other than AP1 or the antibody did not have the specific epitope. The presence of AP1 in the nuclear extract of both cell lines was confirmed by the supershift seen when using AP1 consensus probe. We demonstrated the binding of GATA1 to the G allele at the -689C>G site by supershift with an antibody against GATA1.

These findings suggest the involvement of GATA1 in the transcriptional regulation of PI3 gene in amnion cells and provide a possible genetic explanation for the downregulation of PI3 in chorioamniotic membranes from PPROM cases. We have previously demonstrated that the levels of neutrophil elastase (ELA2, [LocusID: 1991]) are increased in the amniotic fluid of patients with PPROM and acute chorioamnionitis [32]. It is therefore plausible to speculate that the production of PI3 in the fetal membranes is to protect the tissue from the damage that could be caused by increased amounts of neutrophil elastase. Our recent study [16] showing decreased expression of PI3 in the chorioamniotic membranes from patients with PPROM supports our hypothesis that patients who are not capable of producing adequate amounts of PI3 may be predisposed to PPROM.

It has been suggested that PI3 is involved in the pathophysiology of many clinical conditions. For example, PI3 was found in the epidermis of patients with psoriasis, but not in normal human epidermis [33]. Higher levels of PI3 were also observed in bronchial secretions from patients with chronic obstructive pulmonary disease [2] and bronchial carcinoma [7], and the expression of PI3 was decreased in breast [9] and in epidermal tumors [8, 9]. The SNPs identified here will likely be useful for studying the molecular mechanisms of these diseases.

Conclusion

A high degree of polymorphism was detected in the PI3 gene with 23 SNPs, 11 of which are in the promoter region. Two SNP sites (-1063G>A and -689C>G) showed differential binding of transcription factors in nuclear extracts derived from both amnion and HeLa cells suggesting possible involvement of these two SNPs in the expression of PI3 gene. As the SNP site at -1063G>A did not bind to the transcription factor AP1 as suggested by in silico search, the bound transcription factor may not be in current database and needs to be characterized. Binding of GATA1 to the G allele at the -689C>G site suggests the involvement of GATA1 in the transcriptional regulation of PI3 gene in amnion cells. We have performed a genetic association study with PI3 variants, including the -689C>G variant, and found that it is associated with PPROM [manuscript in preparation]. We also previously demonstrated by immunohistochemistry that many cell types of the chorioamniotic membranes produce PI3 and that PI3 protein is decreased in chorioamniotic membranes from PPROM cases [16]. Together, these lines of evidence provide a plausible genetic explanation for the down regulation of PI3 in chorioamniotic membranes from PPROM cases. Previously the involvement of PI3 in the pathophysiology of many clinical conditions was suggested. The SNPs identified here provide the tools for studying the molecular mechanism of these diseases.

References

  1. Schalkwijk J, Wiedow O, Hirose S: The trappin gene family: proteins defined by an N-terminal transglutaminase substrate domain and a C-terminal four-disulphide core. Biochem J. 1999, 340: 569-577. 10.1042/0264-6021:3400569.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Sallenave JM, Silva A: Characterization and gene sequence of the precursor of elafin, an elastase-specific inhibitor in bronchial secretions. Am J Respir Cell Mol Biol. 1993, 8 (4): 439-445.

    Article  CAS  PubMed  Google Scholar 

  3. Clauss A, Lilja H, Lundwall A: A locus on human chromosome 20 contains several genes expressing protease inhibitor domains with homology to whey acidic protein. Biochem J. 2002, 368 (Pt 1): 233-242. 10.1042/BJ20020869.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Saheki T, Ito F, Hagiwara H, Saito Y, Kuroki J, Tachibana S, Hirose S: Primary structure of the human elafin precursor preproelafin deduced from the nucleotide sequence of its gene and the presence of unique repetitive sequences in the prosegment. Biochem Biophys Res Commun. 1992, 185 (1): 240-245. 10.1016/S0006-291X(05)80981-7.

    Article  CAS  PubMed  Google Scholar 

  5. Pol A, Pfundt R, Zeeuwen P, Molhuizen H, Schalkwijk J: Transcriptional regulation of the elafin gene in human keratinocytes. J Invest Dermatol. 2003, 120 (2): 301-307. 10.1046/j.1523-1747.2003.12043.x.

    Article  CAS  PubMed  Google Scholar 

  6. Wiedow O, Schroder JM, Gregory H, Young JA, Christophers E: Elafin: an elastase-specific inhibitor of human skin. Purification, characterization, and complete amino acid sequence. J Biol Chem. 1990, 265 (25): 14791-14795.

    CAS  PubMed  Google Scholar 

  7. Sallenave JM, Marsden MD, Ryle AP: Isolation of elafin and elastase-specific inhibitor (ESI) from bronchial secretions. Evidence of sequence homology and immunological cross-reactivity. Biol Chem Hoppe Seyler. 1992, 373 (1): 27-33.

    Article  CAS  PubMed  Google Scholar 

  8. Alkemade HA, Molhuizen HO, van Vlijmen-Willems IM, van Haelst UJ, Schalkwijk J: Differential expression of SKALP/Elafin in human epidermal tumors. Am J Pathol. 1993, 143 (6): 1679-1687.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Zhang M, Zou Z, Maass N, Sager R: Differential expression of elafin in human normal mammary epithelial cells and carcinomas is regulated at the transcriptional level. Cancer Res. 1995, 55 (12): 2537-2541.

    CAS  PubMed  Google Scholar 

  10. Hiemstra PS, Maassen RJ, Stolk J, Heinzel-Wieland R, Steffens GJ, Dijkman JH: Antibacterial activity of antileukoprotease. Infect Immun. 1996, 64 (11): 4520-4524.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Sallenave JM: Antimicrobial activity of antiproteinases. Biochem Soc Trans. 2002, 30 (2): 111-115. 10.1042/BST0300111.

    Article  CAS  PubMed  Google Scholar 

  12. Simpson AJ, Maxwell AI, Govan JR, Haslett C, Sallenave JM: Elafin (elastase-specific inhibitor) has anti-microbial activity against gram-positive and gram-negative respiratory pathogens. FEBS Lett. 1999, 452 (3): 309-313. 10.1016/S0014-5793(99)00670-5.

    Article  CAS  PubMed  Google Scholar 

  13. Wiedow O, Harder J, Bartels J, Streit V, Christophers E: Antileukoprotease in human skin: an antibiotic peptide constitutively produced by keratinocytes. Biochem Biophys Res Commun. 1998, 248 (3): 904-909. 10.1006/bbrc.1998.9069.

    Article  CAS  PubMed  Google Scholar 

  14. Sallenave JM, Shulmann J, Crossley J, Jordana M, Gauldie J: Regulation of secretory leukocyte proteinase inhibitor (SLPI) and elastase-specific inhibitor (ESI/elafin) in human airway epithelial cells by cytokines and neutrophilic enzymes. Am J Respir Cell Mol Biol. 1994, 11 (6): 733-741.

    Article  CAS  PubMed  Google Scholar 

  15. Bingle L, Tetley TD, Bingle CD: Cytokine-mediated induction of the human elafin gene in pulmonary epithelial cells is regulated by nuclear factor-kappaB. Am J Respir Cell Mol Biol. 2001, 25 (1): 84-91.

    Article  CAS  PubMed  Google Scholar 

  16. Tromp G, Kuivaniemi H, Romero R, Chaiworapongsa T, Kim YM, Kim MR, Maymon E, Edwin S: Genome-wide expression profiling of fetal membranes reveals a deficient expression of protease inhibitor 3 in premature rupture of membranes. American Journal of Obstetrics and Gynecology. 2004, 191: 1331-1338. 10.1016/j.ajog.2004.07.010.

    Article  CAS  PubMed  Google Scholar 

  17. Hall T: BioEdit Sequence Alignment Editor for Windows 95/98/NT/XP. [http://www.mbio.ncsu.edu/BioEdit/bioedit.html]

  18. Schug J, Overton GC: TESS : Transcription Element Search System. [http://www.cbil.upenn.edu/cgi-bin/tess/tess]

  19. Schug J, Overton GC: Modeling transcription factor binding sites with Gibbs Sampling and Minimum Description Length encoding. Proc Int Conf Intell Syst Mol Biol. 1997, 5: 268-271.

    CAS  PubMed  Google Scholar 

  20. Grabe N: AliBaba 2.1. [http://darwin.nmsu.edu/~molb470/fall2003/Projects/solorz/aliBaba_2_1.htm]

  21. Grabe N: AliBaba2: context specific identification of transcription factor binding sites. In Silico Biol. 2002, 2 (1): S1-1.

    PubMed  Google Scholar 

  22. Cartharius K, Frech K, Grote K, Klocke B, Haltmeier M, Klingenhoff A, Frisch M, Bayerlein M, Werner T: Genomatix Software GmbH - GEMS Launcher. [http://www.genomatix.de/products/MatInspector/index.html]

  23. Quandt K, Frech K, Karas H, Wingender E, Werner T: MatInd and MatInspector: new fast and versatile tools for detection of consensus matches in nucleotide sequence data. Nucleic Acids Res. 1995, 23: 4878-4884.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Dignam JD, Lebovitz RM, Roeder RG: Accurate transcription initiation by RNA polymerase II in a soluble extract from isolated mammalian nuclei. Nucleic Acids Res. 1983, 11 (5): 1475-1489.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Komatsu BN, Kirito K, Kashii Y, Furukawa Y, Kikuchi J, Suwabe N, Yamamoto M, Miura Y: Cell-cycle-dependent regulation of erythropoietin receptor gene. Blood. 1997, 89 (4): 1182-1188.

    CAS  PubMed  Google Scholar 

  26. Mutation Nomenclature. [http://www.hgvs.org/mutnomen/]

  27. Povey S: HGNC Home Page. [http://www.gene.ucl.ac.uk/nomenclature]

  28. Excoffier L: Arlequin's Home on the Web. [http://anthro.unige.ch/arlequin/]

  29. dbSNP Home Page. [http://www.ncbi.nlm.nih.gov/SNP/]

  30. Moore JT, Hagstrom J, McCormick DJ, Harvey S, Madden B, Holicky E, Stanford DR, Wieben ED: The major clotting protein from guinea pig seminal vesicle contains eight repeats of a 24-amino acid domain. Proc Natl Acad Sci U S A. 1987, 84 (19): 6712-6714. 10.1073/pnas.84.19.6712.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Zhang M, Magit D, Pardee AB, Sager R: Re-expression of elafin in 21MT2 breast carcinomas by phorbol 12-myristate 13-acetate is mediated by the Ap1 site in the elafin promoter. Cancer Res. 1997, 57 (20): 4631-4636.

    CAS  PubMed  Google Scholar 

  32. Helmig BR, Romero R, Espinoza J, Chaiworapongsa T, Bujold E, Gomez R, Ohlsson K, Uldbjerg N: Neutrophil elastase and secretory leukocyte protease inhibitor in prelabor rupture of membranes, parturition and intra-amniotic infection. J Matern Fetal Neonatal Med. 2002, 12 (4): 237-246.

    Article  CAS  PubMed  Google Scholar 

  33. Schalkwijk J, Chang A, Janssen P, de Jongh GJ, Mier PD: Skin-derived antileucoproteinase (SKALP): characterization of two new elastase inhibitors from psoriatic epidermis. Br J Dermatol. 1990, 122: 631-641.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

We thank Ms. Magdalena Skunca, Ms. Irene Broner, and Ms. Tamara Briazova for expert technical assistance in DNA isolations and the personnel under the directorship of Dr. Susan Land at the Applied Genomics Technology Center for DNA sequencing. This research was supported by the Intramural Research Program of the National Institute of Child Health and Human Development, NIH, DHHS.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Roberto Romero or Gerard Tromp.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Authors' contributions

Mahboob A. Chowdhury designed and interpreted experiments; carried out sequencing, EMSAs and in silico searches; drafted the manuscript and approved the final version of it.

Helena Kuivaniemi designed and interpreted experiments; carried out PCRs; drafted the manuscript and approved the final version of it.

Roberto Romero provided overall direction to the project; was responsible for the clinical data, revised manuscript and approved the final version of it.

Samuel Edwin established amniotic cell cultures; revised manuscript and approved the final version of it.

Tinnakorn Chaiworapongsa provided clinical data and approved the final version of the manuscript.

Gerard Tromp designed and interpreted experiments; carried out all statistical analyses; provided overall direction to the project; drafted and revised the manuscript as well as approved the final version of it.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Chowdhury, M.A., Kuivaniemi, H., Romero, R. et al. Identification of novel functional sequence variants in the gene for peptidase inhibitor 3. BMC Med Genet 7, 49 (2006). https://doi.org/10.1186/1471-2350-7-49

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2350-7-49

Keywords