Skip to main content
  • Research article
  • Open access
  • Published:

A molecular signature of epithelial host defense: comparative gene expression analysis of cultured bronchial epithelial cells and keratinocytes



Epithelia are barrier-forming tissues that protect the organism against external noxious stimuli. Despite the similarity in function of epithelia, only few common protective mechanisms that are employed by these tissues have been systematically studied. Comparative analysis of genome-wide expression profiles generated by means of Serial Analysis of Gene Expression (SAGE) is a powerful approach to yield further insight into epithelial host defense mechanisms. We performed an extensive comparative analysis of previously published SAGE data sets of two types of epithelial cells, namely bronchial epithelial cells and keratinocytes, in which the response to pro-inflammatory cytokines was assessed. These data sets were used to elucidate a common denominator in epithelial host defense.


Bronchial epithelial cells and keratinocytes were found to have a high degree of overlap in gene expression. Using an in silico approach, an epithelial-specific molecular signature of gene expression was identified in bronchial epithelial cells and keratinocytes comprising of family members of keratins, small proline-rich proteins and proteinase inhibitors. Whereas some of the identified genes were known to be involved in inflammation, the majority of the signature represented genes that were previously not associated with host defense. Using polymerase chain reaction, presence of expression of selected tissue-specific genes was validated.


Our comparative analysis of gene transcription reveals that bronchial epithelial cells and keratinocytes both express a subset of genes that is likely to be essential in epithelial barrier formation in these cell types. The expression of these genes is specific for bronchial epithelial cells and keratinocytes and is not seen in non-epithelial cells. We show that bronchial epithelial cells, similar to keratinocytes, express components that are able to form a cross-linked protein envelope that may contribute to an effective barrier against noxious stimuli and pathogens.


Epithelial tissues in the mammalian airways and skin are among the largest organs and form the interface between the internal milieu of the host and the outside world. They not only protect the host against invading pathogens but also provide an effective barrier to noxious external (chemical and physical) stimuli and dehydration [1, 2]. The effectiveness of the epithelial barrier is demonstrated by the rare incidence of severe infections to the lung or skin in healthy individuals. It has become clear that epithelia also play an active role in innate and adaptive immunity [3, 4]. Epithelial tissues display three main mechanisms to protect the organism from infection. First, epithelial cells form an impermeable physical barrier which both prevents pathogen entry and minimizes dehydration. Second, epithelial cells are capable of producing defense molecules such as antimicrobial peptides and proteinase inhibitors. Finally, these cells are able to produce signaling molecules such as cytokines and chemokines. These molecules may attract or activate cells of the innate and adaptive immune system [5, 6]. Interaction between cells of the immune system is mediated by adhesion molecules and cytokine receptors [7, 8] that are present on epithelial cells.

Host defense mechanisms in epithelial cells are coordinated by a complex program of gene expression. Very powerful and sophisticated laboratory techniques such as Serial Analysis of Gene Expression (SAGE) [9] and DNA microarrays [10] have been developed to assess the expression of thousands of genes at the mRNA level in a single experiment. To delineate the barrier function of epithelial cells, the transcriptional change induced by pro-inflammatory cytokines was recently assessed by means of SAGE in two well-established culture models of epithelial inflammation using subcultures of primary bronchial epithelial cells [11] and primary keratinocytes [12]. These independent studies showed a marked overlap in gene families expressed in response to pro-inflammatory cytokines in both cell types. Upon cytokine exposure, in particular genes associated with cytoskeletal architecture and epidermal barrier function such as keratins, S100 calcium-binding proteins and various antimicrobial proteinase inhibitors were differentially expressed. These studies indicated that bronchial epithelial cells and keratinocytes might respond similarly to external influences to ultimately provide effective host protection. This is especially of interest because the epithelia of the skin and conducting airways are markedly different in morphology. The potential functional resemblance of these types of epithelia is also demonstrated by comparative analysis of genetic studies in patients with asthma and atopic dermatitis showing that similar patterns of gene expression may contribute to susceptibility to these diseases [13]. This prompted us to a conduct a comparative analysis of our previously generated gene expression in culture models of epithelial inflammation. The aim was to test the hypothesis whether bronchial epithelial cells and keratinocytes employ similar mechanisms for providing effective host defense at these epithelia.

Therefore, in the present study, our previously generated SAGE data sets derived from bronchial epithelial cells [11] and keratinocytes [12] that were exposed to pro-inflammatory cytokines were compared to identify a common denominator in host defense in the different types of epithelial cells. SAGE libraries of resting and IL1β/TNFα-exposed primary bronchial epithelial cells (~28.000 tags in each library) were compared to SAGE libraries of resting and TNFα-exposed human primary keratinocytes (~13.000 tags in each library). The in silico method Tissue Preferential Expression (TPE) [14] was used for the recognition of putative cell-specific gene expression in these SAGE libraries. Previously, this method has been successfully applied to identify novel specific markers for disease [14, 15]. To verify the in silico prediction analysis of tissue specific gene expression, polymerase chain reaction was performed on seven target genes that were identified by the TPE algorithm in a panel of nine different cell types of which seven are normally present in the airways or lungs. The airway- and lung-derived NCI-H292 and A549 cell lines were included since these cell lines are frequently used to study epithelial cell function. We have identified and validated a signature of specific gene expression for bronchial epithelial cells and keratinocytes. The majority of genes in this signature was previously not associated with host defense or inflammation. These results indicate that epithelia of the airways and skin exploit unified host defense strategies to protect the host, despite their morphological differences.


Transcriptional overlap between PBEC and KC and epithelial-specific gene expression upon cytokine exposure was characterized. By comparing the four SAGE libraries of primary bronchial epithelial cells (PBEC) and keratinocytes (KC), an overlap in tags of approximately 80% was observed indicating a high similarity in the repertoire of genes expressed by these types of epithelial cells. Although remarkable commonalities were found in gene families found to be expressed by PBEC and KC, the repertoire of transcribed family members differed among the two cell types (table 2). To extract a pattern of genes that is specifically expressed in epithelial cells that could likely be involved in epithelial host defense we explored which of the genes are preferentially expressed by PBEC and KC using the TPE algorithm. The scatter plot in figure 1 displays the individual tags observed in the cytokine-exposed PBEC and KC libraries. Each dot represents a single tag with the corresponding TPE values for PBEC and KC. In this analysis, four groups of tags were identified: epithelial non-specific tags (i), tags preferentially expressed by either PBEC (ii) or KC (iii) and tags that were preferentially expressed by both PBEC and KC (iv). The expression of the 30 tags observed in the latter group represents putative epithelial-specific genes because a TPE score ≥ 9 was observed in both PBEC and KC (table 3). Almost half of these tags corresponded to genes encoding for keratins, small proline-rich proteins, kallikreins and proteinase inhibitors (table 3). Interestingly, the expression of a large proportion of these genes was found to be affected by cytokine exposure in PBEC or KC (or both) as observed in the initial SAGE studies (as indicated by underlined tag numbers in table 3). A similar picture in preferential tag expression was obtained when using the libraries of resting PBEC and KC since the majority of genes do not show an on/off expression profile upon stimulation with cytokines (data not shown).

Table 1 Primer sequences and conditions for RT-PCR. Gene sequences used for primer design were retrieved from the Ensembl website [36], annealing temperature and MgCl2 concentration used in PCR reactions are listed.
Table 2 Inventory of gene families and their members expressed by PBEC and/or KC.
Table 3 30 Epithelial-specific genes as identified by the TPE analysis.
Figure 1
figure 1

TPE scatter plot of SAGE tags of PBEC and KC libraries after cytokine exposure. Reliable 3'-end tags with TPE>4 and tag frequency of = 2 in at least one library were plotted. Tags with corresponding TPE values = 9 in both libraries were considered to be potential epithelial cell-specific tags as indicated by the threshold lines in the figure. A similar picture was obtained when TPE values of tags from the resting libraries were plotted.

To validate this in silico TPE prediction analysis, expression of seven putative epithelial-specific genes by PBEC and KC was assessed by reverse transcriptase polymerase chain reaction (RT-PCR) in nine different cell types. Each cell type in the panel was exposed to medium alone or to IL1β/TNFα. KC were exposed to medium or TNFα alone to maintain comparability with the original SAGE experiment. In concordance with the SAGE data and TPE analysis, expression for SPRR2A was only observed in PBEC. On the other hand, CALML5 was expected to be expressed by KC alone. However, PBEC were shown to be positive for this transcript as well and weak expression was observed in NCI-H292 cells. As demonstrated by the TPE analysis KRT6A, SPRR1A, SPRR1B, IL1F9, S100A2 all showed TPE values of ≥ 9 in both PBEC and KC libraries. The RT-PCR results in figure 2 demonstrates that preferential expression of SPRR1B was found in PBEC, KC and NCI-H292 cells, whereas moderate to weak expression was also detected in fibroblasts, HUVEC, HASM and monocytes. Expression of KRT6A is restricted to PBEC, KC and the bronchial epithelial cell line NCI-H292, whereas expression of this transcript was negative in all other cell types. Transcription of SPRR1A, IL1F9 and S100A2 was only detected in primary cultures of PBEC and KC and was completely absent in all other cell types.

Figure 2
figure 2

PCR verification of 7 potential preferentially expressed tags identified using the TPE algorithm. The expression of these seven genes was assessed both under resting conditions and after cytokine exposure in all cell types. On the right, the TPE values of the gene are listed for both PBEC and KC after cytokine exposure. Tags for which no TPE value could be calculated because of absence of the tag in the particular library are indicated by "not available" (N/A). The predicted preferential expression could be verified for all genes. The expression of six of these genes seems to be selective for epithelial cells only. Whereas SPRR1B is preferentially expressed by epithelial cells, moderate to low levels of expression were also detected in other cell types as well. SPRR2A is preferentially expressed by PBEC only. No tags for SPRR2A were found in the KC libraries after TNFα exposure whereas expression for CALML5 was observed by RT-PCR in both PBEC and KC, while no tags for this gene were found in PBEC libraries.


Comparative genomics approaches have the potential to gain additional insight into a biological process at the mRNA expression level by integrating and combining data obtained from similar model systems. Particularly, SAGE is excellent for this purpose since digital, scalable expression data is generated that allows comparison without the need for complex mathematical normalization methods. Although the SAGE libraries used in the present analysis were not initially intended for comparative genomic research, remarkable commonalities in epithelial-specific gene expression were found that related to host defense.

The tissue preferential expression (TPE) algorithm was employed to recognize specific tag expression by PBEC and KC under inflammatory conditions (pane iv; figure 1, table 3). Experimental verification of selected epithelial-specific genes by RT-PCR showed a good correlation between the in silico approach and RT-PCR (figure 2). The PCR setup was designed to detect true presence or absence of validation genes and was not intended to be quantitative. The observed discrepancies between SAGE and PCR results can be explained by the difference in detection sensitivity between techniques: RT-PCR is far more sensitive than SAGE in detecting low abundant gene expression.

The majority of tags of the molecular signature corresponded to genes encoding structural components of the cytoskeleton (keratins, small proline-rich proteins, elafin) and for proteins that are involved in the assembly/disassembly (transglutaminase 1, kallikreins and matrixmetalloproteinases) of the cornified cell envelope in keratinocytes (reviewed in [16]). Components of the cross-linked or cornified envelope are linked by transglutaminases (reviewed in [17]). The observation that bronchial epithelial cells express components of and assembly/disassembly enzymes forming a cross-linked envelope is relevant to our understanding of epithelial host defense in the airways. Additional support for this observation is provided by abundant transcription of genes that are known to be involved in cornification in skin, including the S100 calcium-binding proteins [18], annexins [18] and cystatins [19, 20] (table 2). So far, only few studies provide evidence for the existence of a protein envelope in bronchial epithelial cells. Components such as small proline-rich proteins (SPRR) have been suggested to be associated with squamous differentiation [21, 22]. Low SPRR expression has been associated with squamous cell carcinoma [23], whereas high expression of SPRR1B enhanced G0-arrest resulting in growth arrest [24]. Interestingly, the families of small proline-rich proteins and S100 calcium-binding proteins are encoded in the epidermal differentiation complex (EDC) [25, 26]. Proteins encoded in this region share significant sequence similarities, particularly in the glutamine- and lysine-rich regions that are involved in the cross-linking by transglutaminases [17]. This indicates that PBEC and KC not only share structural characteristics, but may also share functional characteristics.

A disadvantage of the present study might be the differences in type of cytokine-exposure and duration of the treatment. The opposite directional changes in expression in gene families (table 2) observed could be explained either by dissimilarities in the initial model systems or by the inherent differences between PBEC and KC. By using the TPE algorithm, highly cell-specific tag expression can be predicted largely independently from transcriptional levels because the more unique a tag is to a particular tissue, the less important is its level of expression. Therefore, we are confident that the signature of epithelial host defense that was extracted is representative for bronchial epithelial cells and keratinocytes.

Computational subtraction methods such as the TPE algorithm allow functional clustering of genes derived from large and complex genome-wide expression profiles without having full knowledge of the repertoire of genes involved in biological processes of interest. Although the identified molecular signature of host defense is characteristic for bronchial epithelial cells and keratinocytes, it would be of great interest to study whether this gene expression pattern is also applicable to other types of epithelial cells, a finding that would greatly enhance our understanding of epithelial defense strategies.


In summary, our comprehensive comparison of overlapping genes across bronchial epithelial cells and keratinocytes provides novel insights in epithelial host defense strategies, in particular of the airway epithelium. Combining in silico and experimental approaches is very valuable in accelerating the interpretation of genomics data and defining follow-up research. We identified an expression signature of genes that were specifically expressed by bronchial epithelial cells and keratinocytes. These genes are likely to fulfill an eminent function in epithelial host defense. Based on the present findings we propose that formation of a cross-linked protein envelope by bronchial epithelial cells is an effective host defense strategy of the mucosal epithelium in the human airways. This function would be analogous to the host defense function of cornifying keratinocytes. Finally, a better understanding of unified host defense strategies in different epithelia may lead to the identification of novel therapeutic targets for epithelial inflammatory disorders such as asthma and atopic dermatitis.


SAGE data

The previously published SAGE libraries that were compared in this study were derived from two models of epithelial inflammation using primary bronchial epithelial cells [11] and primary keratinocytes [12]. The SAGE data from the original studies is accessible through NCBI's Gene Expression Omnibus [27] with GEO accessions GSM37337 (PBEC_unstimulated), GSM37339 (PBEC_IL1beta/TNFalpha), GSM1121 (NormCultKC_Diff) and GSM1122 (TNF_AlphaCultKC).

For tag mapping, after discarding tags occurring only once, the libraries were compared with NCBI's "reliable Unigene cluster to SAGE tag map" [28] and with SAGEgenie of the Cancer Genome Anatomy Project [29]. Both maps were based on Unigene build#171. Additionally, to enhance the reliability of tag identity we included the virtual tag classification as used in SAGEgenie to assess the location of each tag within the corresponding transcript. Reliable tags can be discriminated from tags that are not isolated from the 3'-end such as internally primed transcripts and tags derived from internal NlaIII restriction sites [29, 30]. In the data set comparisons, only tags were included that were derived from the most 3'- restriction site of NlaIII, tags that matched to undefined 3'-end transcripts and tags for which no additional information was available. The last category may contain tags that correspond to novel transcripts.

TPE analysis

Epithelial-specific gene expression in PBEC and KC was identified using the Tissue Preferential Expression (TPE) algorithm [14]. The calculated Tissue Preferential Expression (TPE) value is based both on the presence of a particular tags and its level of expression in the SAGE library of interest in comparison to a panel of reference SAGE libraries derived from a range of different whole tissues. To allow calculation of TPE values, each of the PBEC and KC SAGE libraries as well as the reference libraries were normalized to a equal number of tags per library. TPE values were determined for each tag in the individual PBEC and KC libraries by selecting each of the libraries as library of interest before applying the TPE algorithm. After calculation, the TPE values were ranked according to their value. Large positive TPE values represent tissue-specific genes that are overexpressed in the PBEC and/or KC libraries. Tags with TPE values of <4 were excluded from further analysis since these tags occur very frequently in other cell types as well (See supplement for a detailed description of the TPE algorithm). The threshold value for the TPE analysis that is indicative for tissue-specific expression was chosen very high to prevent possible false positives. Tags with a corresponding TPE value of ≥ 9 is indicative for tissue preferential expression in the epithelial cells used in this analysis. See additional file for a detailed description of the TPE algorithm.

Cell culture

Primary Bronchial Epithelial Cells (PBEC) and Primary Keratinocytes (KC)

Subcultures of human primary bronchial epithelial cells and human primary keratinocytes were derived and cultured as described previously [31, 32].

A549 and NCI-H292 cells

the lung derived epithelial cell lines A549 (CCL-185, American Type Culture Collection, Rockville, MD, USA) and NCI-H292 (CRL-1848, American Type Culture Collection) were cultured according to the supplier's recommendation. Prior to stimulation, cells were cultured overnight in serum free medium.

Human airway smooth muscle cells (HASM)

Human airway smooth muscle cells (HASM) from two donors were purchased from Stratagene (La Jolla, CA) and were cultured as described previously [33].

Human mast cells (HMC-1)

HMC-1 were kindly provided by J.H. Butterfield [34] and were cultured in IMDM medium containing 25 mM Hepes, 2 mM L-glutamine, 20 U/ml penicillin, 20 μg/ml streptomycin (all from Bio Whittaker, Walkersville, MD), 5 μg/ml apo-transferrin, 0,36% β-mercaptoethanol and 10% heat-inactivated Fetal Calf Serum (FCS; GibcoBRL/Life Technologies, Breda, The Netherlands). Prior to stimulation, cells were cultured overnight in serum free medium (same as above, without heat-inactivated FCS)

Human lung fibroblasts (HFL-1)

HFL-1 (CCL-153, American Type Culture Collection) were cultured according to the supplier's recommendations. Prior to stimulation, cells were cultured overnight in serum free medium.


CD14-purified monocytes were kindly provided by the department of Nephrology (Leiden University Medical Center, Leiden, The Netherlands), and were resuspended in RPMI 1640 medium supplemented with 20 U/ml penicillin, 20 μg/ml streptomycin, 2 mM glutamine (all from Bio Whittaker) and 10% heat-inactivated FCS and were seeded in 6-wells plates to allow adherence. After 2 hours, the medium was replaced by serum free medium (same as above, without heat-inactivated FCS) for overnight incubation prior to stimulation.

Human umbilical vein endothelial cells (HUVEC)

HUVEC were kindly provided by the department of Nephrology, Leiden University Medical Center, Leiden, The Netherlands. All cell cultures were performed at 37°C, 5% CO2 and 95% relative humidity.

Stimulation of cell cultures

All cell types, except for the keratinocytes, were stimulated for 6 hours with either medium alone, a mixture of the pro-inflammatory cytokines IL1β (20 ng/ml; PeproTech, Rocky Hill, NJ) and TNFα (20 ng/ml; PeproTech). Keratinocytes were stimulated with TNFα (25 ng/ml) for 48 hours as described previously [12].

Reverse Transcription PCR

RT-PCR was used to verify preferential expression of genes identified in the TPE analysis. Total RNA from the nine different cell cultures was extracted using the RNeasy mini kit (Qiagen, Westburg, Leusden, The Netherlands) and on-column DNA digestion was performed with DNase I (Qiagen, Westburg), all according to the manufacturer's instructions. Single-stranded cDNA was synthesized of total RNA using M-MLV Reverse Transcriptase primed with Oligo-dT (both from Invitrogen/Life Technologies, Breda, The Netherlands) in the presence of a RNase inhibitor (RNaseOUT; Invitrogen/Life Technologies) according to the manufacturer's instructions. Gene-specific primers were designed for keratin 6A (KRT6A), small proline-rich protein and 1B (SPRR1B), the calcium-binding protein S100A2, IL1 family member 9 (IL1F9), calmodulin-like 5 (CALML5) and β-actin (ACTB) as internal control (table 1). Primers for the small proline-rich protein 1A and 2A [35] were kindly provided by Claude Backendorf. Primers were synthesized by Isogen (Maarssen, The Netherlands). All PCR reactions were carried out according to the following PCR conditions: initial denaturation of 3 minutes at 95°C, then 35 cycles of 15 seconds denaturing at 95°C, 15 seconds primer annealing, 30 seconds elongation at 72°C, and a final extension of 3 minutes at 72°C in the last cycle. For β-actin, the cycle number was limited to 25. PCR products were visualized on 1% agarose gels using ethidium bromide. Both the reverse transcription and PCR reactions were performed on a Biometra T-Gradient thermocycler (Biometra GmbH, Goettingen, Germany).


  1. Presland RB, Jurevic RJ: Making sense of the epithelial barrier: what molecular biology and genetics tell us about the functions of oral mucosal and epidermal tissues. J Dent Educ. 2002, 66: 564-574.

    PubMed  Google Scholar 

  2. Thompson AB, Robbins RA, Romberger DJ, Sisson JH, Spurzem JR, Teschler H, Rennard SI: Immunological functions of the pulmonary epithelium. Eur Respir J. 1995, 8: 127-149. 10.1183/09031936.95.08010127.

    Article  PubMed  CAS  Google Scholar 

  3. Bals R, Hiemstra PS: Innate immunity in the lung: how epithelial cells fight against respiratory pathogens. Eur Respir J. 2004, 23: 327-333. 10.1183/09031936.03.00098803.

    Article  PubMed  CAS  Google Scholar 

  4. Diamond G, Legarda D, Ryan LK: The innate immune response of the respiratory epithelium. Immunol Rev. 2000, 173: 27-38. 10.1034/j.1600-065X.2000.917304.x.

    Article  PubMed  CAS  Google Scholar 

  5. Takizawa H: Airway epithelial cells as regulators of airway inflammation (Review). Int J Mol Med. 1998, 1: 367-378.

    PubMed  CAS  Google Scholar 

  6. Laberge S, El Bassam S: Cytokines, structural cells of the lungs and airway inflammation. Paediatr Respir Rev. 2004, 5 Suppl A: S41-S45. 10.1016/S1526-0542(04)90009-7.

    Article  PubMed  Google Scholar 

  7. Pilewski JM, Albelda SM: Adhesion molecules in the lung. An overview. Am Rev Respir Dis. 1993, 148: S31-S37.

    Article  PubMed  CAS  Google Scholar 

  8. Barker JN, Mitra RS, Griffiths CE, Dixit VM, Nickoloff BJ: Keratinocytes as initiators of inflammation. Lancet. 1991, 337: 211-214. 10.1016/0140-6736(91)92168-2.

    Article  PubMed  CAS  Google Scholar 

  9. Velculescu VE, Zhang L, Vogelstein B, Kinzler KW: Serial analysis of gene expression. Science. 1995, 270: 484-487.

    Article  PubMed  CAS  Google Scholar 

  10. Schena M, Shalon D, Davis RW, Brown PO: Quantitative Monitoring of Gene-Expression Patterns with A Complementary-Dna Microarray. Science. 1995, 270: 467-470.

    Article  PubMed  CAS  Google Scholar 

  11. Vos JB, van Sterkenburg MA, Rabe KF, Schalkwijk J, Hiemstra PS, Datson NA: Transcriptional response of bronchial epithelial cells to Pseudomonas aeruginosa: identification of early mediators of host defense. Physiol Genomics. 2005, 21: 324-336. 10.1152/physiolgenomics.00289.2004.

    Article  PubMed  CAS  Google Scholar 

  12. Jansen BJ, van Ruissen F, de Jongh G, Zeeuwen PL, Schalkwijk J: Serial analysis of gene expression in differentiated cultures of human epidermal keratinocytes. J Invest Dermatol. 2001, 116: 12-22. 10.1046/j.1523-1747.2001.00218.x.

    Article  PubMed  CAS  Google Scholar 

  13. Cookson W: The immunogenetics of asthma and eczema: a new focus on the epithelium. Nat Rev Immunol. 2004, 4: 978-988. 10.1038/nri1500.

    Article  PubMed  CAS  Google Scholar 

  14. Moreno JC, Pauws E, Van Kampen AH, Jedlickova M, de Vijlder JJ, Ris-Stalpers C: Cloning of tissue-specific genes using serial analysis of gene expression and a novel computational substraction approach. Genomics. 2001, 75: 70-76. 10.1006/geno.2001.6586.

    Article  PubMed  CAS  Google Scholar 

  15. Pauws E, Veenboer GJ, Smit JW, de Vijlder JJ, Morreau H, Ris-Stalpers C: Genes differentially expressed in thyroid carcinoma identified by comparison of SAGE expression profiles. FASEB J. 2004, 18: 560-561.

    PubMed  CAS  Google Scholar 

  16. Nemes Z, Steinert PM: Bricks and mortar of the epidermal barrier. Exp Mol Med. 1999, 31: 5-19.

    Article  PubMed  CAS  Google Scholar 

  17. Kalinin AE, Kajava AV, Steinert PM: Epithelial barrier function: assembly and structural features of the cornified cell envelope. Bioessays. 2002, 24: 789-800. 10.1002/bies.10144.

    Article  PubMed  CAS  Google Scholar 

  18. Robinson NA, Lapic S, Welter JF, Eckert RL: S100A11, S100A10, annexin I, desmosomal proteins, small proline-rich proteins, plasminogen activator inhibitor-2, and involucrin are components of the cornified envelope of cultured human epidermal keratinocytes. J Biol Chem. 1997, 272: 12035-12046. 10.1074/jbc.272.18.12035.

    Article  PubMed  CAS  Google Scholar 

  19. Abrahamson M, Alvarez-Fernandez M, Nathanson CM: Cystatins. Biochem Soc Symp. 2003, 179-199.

    Google Scholar 

  20. Zettergren JG, Peterson LL, Wuepper KD: Keratolinin: the soluble substrate of epidermal transglutaminase from human and bovine tissue. Proc Natl Acad Sci U S A. 1984, 81: 238-242.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  21. An G, Huang TH, Tesfaigzi J, Garcia-Heras J, Ledbetter DH, Carlson DM, Wu R: An unusual expression of a squamous cell marker, small proline-rich protein gene, in tracheobronchial epithelium: differential regulation and gene mapping. Am J Respir Cell Mol Biol. 1992, 7: 104-111.

    Article  PubMed  CAS  Google Scholar 

  22. Tesfaigzi J, Th'ng J, Hotchkiss JA, Harkema JR, Wright PS: A small proline-rich protein, SPRR1, is upregulated early during tobacco smoke-induced squamous metaplasia in rat nasal epithelia. Am J Respir Cell Mol Biol. 1996, 14: 478-486.

    Article  PubMed  CAS  Google Scholar 

  23. Luo A, Kong J, Hu G, Liew CC, Xiong M, Wang X, Ji J, Wang T, Zhi H, Wu M, Liu Z: Discovery of Ca2+-relevant and differentiation-associated genes downregulated in esophageal squamous cell carcinoma using cDNA microarray. Oncogene. 2004, 23: 1291-1299. 10.1038/sj.onc.1207218.

    Article  PubMed  CAS  Google Scholar 

  24. Tesfaigzi Y, Wright PS, Belinsky SA: SPRR1B overexpression enhances entry of cells into the G0 phase of the cell cycle. Am J Physiol Lung Cell Mol Physiol. 2003, 285: L889-L898.

    Article  PubMed  CAS  Google Scholar 

  25. Gibbs S, Fijneman R, Wiegant J, van Kessel AG, van de PP, Backendorf C: Molecular characterization and evolution of the SPRR family of keratinocyte differentiation markers encoding small proline-rich proteins. Genomics. 1993, 16: 630-637. 10.1006/geno.1993.1240.

    Article  PubMed  CAS  Google Scholar 

  26. Mischke D, Korge BP, Marenholz I, Volz A, Ziegler A: Genes encoding structural proteins of epidermal cornification and S100 calcium-binding proteins form a gene complex ("epidermal differentiation complex") on human chromosome 1q21. J Invest Dermatol. 1996, 106: 989-992. 10.1111/1523-1747.ep12338501.

    Article  PubMed  CAS  Google Scholar 

  27. Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. Lash AE, Tolstoshev CM, Wagner L, Schuler GD, Strausberg RL, Riggins GJ, Altschul SF: SAGEmap: a public gene expression resource. Genome Res. 2000, 10: 1051-1060. 10.1101/gr.10.7.1051.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  29. Boon K, Osorio EC, Greenhut SF, Schaefer CF, Shoemaker J, Polyak K, Morin PJ, Buetow KH, Strausberg RL, De Souza SJ, Riggins GJ: An anatomy of normal and malignant gene expression. Proc Natl Acad Sci U S A. 2002, 99: 11287-11292. 10.1073/pnas.152324199.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  30. Liang P: SAGE Genie: a suite with panoramic view of gene expression. Proc Natl Acad Sci U S A. 2002, 99: 11547-11548. 10.1073/pnas.192436299.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  31. van Wetering S, van der Linden AC, van Sterkenburg MA, de Boer WI, Kuijpers AL, Schalkwijk J, Hiemstra PS: Regulation of SLPI and elafin release from bronchial epithelial cells by neutrophil defensins. Am J Physiol Lung Cell Mol Physiol. 2000, 278: L51-L58.

    PubMed  CAS  Google Scholar 

  32. van Ruissen F, de Jongh GJ, Zeeuwen PL, Van Erp PE, Madsen P, Schalkwijk J: Induction of normal and psoriatic phenotypes in submerged keratinocyte cultures. J Cell Physiol. 1996, 168: 442-452. 10.1002/(SICI)1097-4652(199608)168:2<442::AID-JCP23>3.0.CO;2-3.

    Article  PubMed  CAS  Google Scholar 

  33. Zuyderduyn S, Hiemstra PS, Rabe KF: TGF-beta differentially regulates TH2 cytokine-induced eotaxin and eotaxin-3 release by human airway smooth muscle cells. J Allergy Clin Immunol. 2004, 114: 791-798. 10.1016/j.jaci.2004.06.037.

    Article  PubMed  CAS  Google Scholar 

  34. Butterfield JH, Weiler DA: In vitro sensitivity of immature human mast cells to chemotherapeutic agents. Int Arch Allergy Appl Immunol. 1989, 89: 297-300.

    Article  PubMed  CAS  Google Scholar 

  35. Cabral A, Voskamp P, Cleton-Jansen AM, South A, Nizetic D, Backendorf C: Structural organization and regulation of the small proline-rich family of cornified envelope precursors suggest a role in adaptive barrier function. J Biol Chem. 2001, 276: 19231-19237. 10.1074/jbc.M100336200.

    Article  PubMed  CAS  Google Scholar 

  36. Ensembl Genome Browser []. 2005

  37. NCBI Gene Expression Omnibus []. 2005

Download references


We would like to thank Marianne van Sterkenburg (Department of Pulmonology, LUMC, Leiden, The Netherlands) for technical support, Melinda Oroszlàn for providing the HUVEC, Sandra van der Kooy (Department of Nephrology, LUMC, Leiden, The Netherlands) for providing the CD14+ Monocytes and René Lutter (Department of Pulmonology, AMC, Amsterdam, The Netherlands) for providing the HMC-1 mast cells. We also would like to thank Claude Backendorf (Department of Molecular Genetics, Leiden Institute of Chemistry, Leiden, The Netherlands) for providing the SPRR1A and 2A primers and for critically reading the manuscript. This research was supported by grant 902-11-092 and 903-42-085 from the Dutch Organisation for Scientific Research (NWO)

Author information

Authors and Affiliations


Corresponding author

Correspondence to Joost B Vos.

Additional information

Authors' contributions

JBV was primarily responsible for study design and coordination, building comparative SAGE expression databases, performing computational analyses and drafting the manuscript. NAD and PSH participated in study design, coordination and manuscript editing. AHvK and ACL advised in study design and performed the TPE computational analyses. DO, JS and PLZ provided the keratinocyte SAGE data, supplied samples and reagents and assisted in editing the manuscript. RMV performed cell culture experiments, PCR analyses and provided comments and discussion. KFR provided useful discussion and manuscript editing. All authors read and approved the final manuscript.

Electronic supplementary material


Additional File 1: A detailed description of the Tissue Preferential Algorithm has been provided as supplemental material in a Word-document entitled "Materials supplement.doc". (DOC 26 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

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 ( ), 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

Vos, J.B., Datson, N.A., van Kampen, A.H. et al. A molecular signature of epithelial host defense: comparative gene expression analysis of cultured bronchial epithelial cells and keratinocytes. BMC Genomics 7, 9 (2006).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: