|Home | About | Journals | Submit | Contact Us | Français|
Modern functional genomic approaches may help to better understand the molecular events involved in tissue morphogenesis and to identify molecular signatures and pathways. We have recently applied transcriptomic profiling to evidence molecular signatures in the development of the normal chicken chorioallantoic membrane (CAM) and in tumor engrafted on the CAM. We have now extended our studies by performing a transcriptome analysis in the "wound model" of the chicken CAM, which is another relevant model of tissue morphogenesis.
To induce granulation tissue (GT) formation, we performed wounding of the chicken CAM and compared gene expression to normal CAM at the same stage of development. Matched control samples from the same individual were used. We observed a total of 282 genes up-regulated and 44 genes down-regulated assuming a false-discovery rate at 5% and a fold change > 2. Furthermore, bioinformatics analysis lead to the identification of several categories that are associated to organismal injury, tissue morphology, cellular movement, inflammatory disease, development and immune system. Endothelial cell data filtering leads to the identification of several new genes with an endothelial cell signature.
The chick chorioallantoic wound model allows the identification of gene signatures and pathways involved in GT formation and neoangiogenesis. This may constitute a fertile ground for further studies.
Different physiological as well as pathological conditions trigger tissue remodeling including surgery, infection, chemical or physical burns, ischemia or immunological reaction . The restoration of tissue integrity involves alteration in tissue elasticity, interstitial fluid pressure and oxygen tension, which is normalized by vascularization of the affected region . Revascularization is accomplished by the ingrown of the granulation tissue (GT) that is composed of a dense network of enlarged vessels forming specific and leaky temporary vasculature . When not disturbed, GT vasculature is normalized during course of scarification. The healing process proceeds according to that general pattern e.g. in the skin but also during regenerative healing after brain or myocardium stroke . Wound healing can be perturbed by pathological changes that include ulceration, hypertrophic scaring or keloids formation and fibrosis . Modern therapy requires the targeting of drugs directly to the site of interest and to accomplish that goal in systemic treatment, the molecular signatures distinguishing the expanding vasculature of the GT from the normal vessels need to be known.
The chicken embryo model has been widely used in developmental biology to understand vascular development and to test the effect of molecules predicted to interfere with the angiogenic process or lymphangiogenesis . For example, the effect of flow on vessel ontology such as venous or arterial patterning has been elucidated using the chicken chorioallantoic membrane (CAM) . Furthermore, the effect of different angiogenesis stimulators such as VEGF-A, VEGF-C or inhibitors has been tested in the chick embryo. Adult wound healing involves movement from the epidermis and connective tissue and the recruitment of inflammatory and immunocompetent cells. Embryonic wound healing also involves wound contraction, followed by re-epithelialization but without recruitment of immunocompetent cells. The inflammatory response in wound healing is crucial for fighting infection so that tissue damage does not lead to death through septicaemia. But, aside from this role, recruitment of leukocytes may more negatively impact wound healing. Indeed, knockout and knockdown studies suggest that immun cells do not promote wound healing and their depletion can even enhance it [1,8]. Thus, models of embryonic wound healing will evidence gene regulations that are crucial for the healing process and independent from the perturbation induced by immunocompetent cells. However, neutrophile-like inflammatory cells and monocyte-like cells are accumulated in growth factor-stimulated CAM which may participate in the wound healing process.
Kilarski et al.  have developed a method to investigate GT formation in the CAM. This model has allowed a better understanding of the formation of the vasculature during GT formation . The major advantages of CAM wound healing model is that the CAM is composed of blood vessels and enclosed within 2 layers of epithelium and a fibroblast matrix. This is in contrast to a skin model, for example, where there are multiple cell and tissue types (epidermis, dermis and subdermis). Wound healing in the CAM model reflects primary changes in vasculature and in stromal fibroblasts that is not affected by "noise" from other cell types. Furthermore, a matched control can be obtained from the same CAM.
Modern functional genomics approaches may facilitate a better understanding of the molecular events involved in tissue morphogenesis and allow the identification of molecular signatures and pathways. We have recently applied transcriptomic profiling to elucidate the molecular signatures involved in the development of the normal chicken chorioallantoic membrane and in tumors engrafted onto the CAM [11,12]. Furthermore, engraftment of human tumour tissue onto the CAM, followed by transcriptomic analyses with both human and chicken microarrays, enables the gene signatures of both the host stroma and the human tumour to be distinguished. We have now extended our earlier studies by performing a transcriptome analysis in the "wound model" of the chicken CAM. This has allowed us to identify gene signatures involved in GT formation and neoangiogenesis. These results further indicate that the chicken embryo model is an excellent tool for discovering networks that are associated with granulation-tissue formation and tissue repair.
Chick embryos were cultured for 10 days and CAMs were inflicted by parallel scalpel superficial cuts of 1 cm area, followed with a subsequent scarping off of the epithelium of the injured chorioallantoic membrane. The wound area was then covered with 1.5 cm square nylon grid (figure (figure1A).1A). As seen in figure figure1B,1B, significant GT was formed and grew through the nylon grid in response to injury. This GT contained a significant amount of blood vessels with abnormal morphology (tortuous blood vessels). Newly formed blood vessels within the CAM became functional since injection of Indian ink clearly evidenced newly perfused blood vessels in the GT (figure (figure1C1C).
The CAM wound model has been established to analyze GT formation and the role of invading fibroblasts and blood vessels in this process . It has been found that tissue tension generated by activated fibroblasts or myofibroblasts during wound contraction, mediated and directed translocation of the vasculature. This vasculature can be expanded secondarily by elongation and vessel enlargement, and finally, through splitting and sprouting. To verify the presence of myofibroblasts in our experimental set up, we performed immunohistology using anti-α2-somooth muscle actin antibodies. The invasion of α2-smooth muscle (α2-SM) actin positive myofibroblasts and of blood vessels in the wounded area was clearly visible and is depicted in figure figure22.
Other chick wound models that have been proposed such as epithelial regeneration models at the surface of the embryo such as at the wing bud or the midbrain region [8,13,14]. These models have helped to characterize some of the morphological and molecular events occurring during embryonic tissue repair involving actin cable assembly and the Rho kinases[8,14]. Tissue wound contraction is present in this model to some extent, however without the presence of α2-SM positive myofibroblasts. Another model is characterized by the removal of only the peridermal layer. Wound closure in this model is essentially driven by the conversion of the basal layer, from monolayer to multilayer .
The advantage of our model is that it clearly distinguishes between preexisting and newly formed tissue and vasculature and that an ingrowth of α2-SM actin positive myofibroblasts is observed. Furthermore, recruitment and translocation of the vasculature in the wound area can be clearly envisioned. This allows us to perform transcriptomic analysis after wounding to establish which genes are important players in this process.
We therefore extracted RNA from eggs 6 days post-wounding. We also extracted RNA from areas where no injury was inflicted to the CAM from the same eggs. Thus, each wound had a matched control from the same egg and three eggs were used for hybridization on each chip with their own controls (6 samples in total). Data were normalized using the robust multi-averaging method, part of the Affymetrix library in the statistical programming language R. Hierarchical clustering, performed on the three wounded and non-wounded samples, indicated separate clustering of wounded and non-wound CAM tissue (figure (figure3A).3A). This indicates good quality of the samples collected and of hybridizations. We observed a total of 282 genes up-regulated and 44 genes down-regulated, applying a FDR < 5% and a fold change > 2 (figure (figure3B3B).
Up-regulation as high as 79 fold was observed. Among the most up-regulated genes, we observed: fatty acid binding protein 4 FABP4 (79 fold), retinol binding protein 7 (RBP7) (74 fold), transthyretin (TTR) (37fold), osteopontin (SPP1) (30 fold), neutrophile cytosolic factor (NCF2) (25 fold), chemokine ah221 (LOC417536) (10 fold) and cysteine-rich secretory protein (CRISP3 (8 fold)) (Table (Table1).1). Among the most down-regulated genes, we observed: inter-alpha (globulin) inhibitor H5 (ITIH5) (0.11 fold), Collagen type VIII a1 (0.15 fold), testican (SPOCK1) (0.29 fold), Atonal homolog 8 (ATOH8) (0.29 fold), laminin a1 (LAMA1) (0.32 fold), CXCL12 (0.32 fold), C1q and TNF-related protein 1 (C1QTNF1) (0.32 fold) and plexin A2 (PLXNA2) (0.33 fold) (Table (Table22).
In order to analyze potential functional trends in the gene lists identified, we performed a functional analysis of the lists of differentially expressed genes. With the purpose to describe the overall representation of functions differentially modulated in response to wounding, we first computed the number of genes represented in the most common gene ontology (GO) categories. As seen in figure figure4,4, the most abundant categories found in this analysis were "binding", "catalytic activity", "cell part", "cellular processes" and "metabolic processes". There was a slight preponderance of up-regulated genes within the first three categories. However, within "metabolic processes" the majority of genes were up-regulated rather than down-regulated.
Furthermore, differences in four other categories related to extracellular matrix and adhesion ("biological adhesion", "extracellular matrix", "extracellular region" and "extracellular region part") were quite striking. In these categories, down-regulated genes were much more abundant. In addition to this analysis, and in order to identify functional categories enriched in our gene lists, we also performed a functional analysis using the web based tool, DAVID [16,17]. This procedure clusters genes with similar functions (defined by similar Gene ontology or KEGG pathway terms) into annotation clusters (Additional file 1). It then assesses whether these "functional clusters" are over-represented in the list of genes differentially regulated. We submitted the up and down regulated genes separately and summarized the results in table table3.3. Due to the relatively low number of significant genes at 5% FDR level, we have also verified that results of the analysis were similar at less stringent FDR thresholds (10% and 15%FDR). Generally no new functions emerged when the FDR was increased but as expected, each function increased in significance and enrichment.
The expression of a number of cytokines or chemokines was up-regulated after wounding. Some of them are represented in the GO term "Cytokines-Chemokines". Among the others clearly associated with this category are IL10R, IL4R, chemokine ligand 20 (CCL20), CX3CR1 and IL1β. Components of the major histocompatibility complex (MHC) were also found up-regulated after wounding. The latter may represent antigen presentation as a consequence to injury. This transcriptional response is not likely to be part of the immune system but may be associated to other cellular types of the GT such as endothelial cells or fibroblasts. Interestingly, the expression of some components of the extracellular matrix and genes involved in developmental processes were decreased. These include Wingless-type MMTV integration site family member 2B (WNT2B), fibrillin-1(FBN1), laminin-α1 (LAMA1), collagen VIII (COL8A1), FRAS1-related extracellular matrix-1 (FREM1), Cysteine-rich transmembrane BMP regulator 1 (CRIM1), Semaphorin 3G (SEMA3G) or Eph receptor A7 (EphA7). This indicates that synthesis of these components is no more required when significant GT formation has occurred. For the list of individual genes represented in our gene lists and belonging to these enriched categories, see additional files 2 and 3.
In order to identify the structure of regulatory networks underlying response to wounding we performed Ingenuity Pathway Analysis (IPA). IPA identifies gene interaction networks representing potential regulatory pathways by integrating lists of differentially expressed genes with a vast public domain literature database, representing several types of gene-gene interactions. Contrary to gene ontology analysis IPA networks represented gene interactions linked to specific mechanisms (e.g. transcriptional activation, protein-protein interactions, etc).
Additional file 4 summarizes the results of this analysis by listing the most significant networks identified. The Ingenuity category with the best scores (> 20) were "Cardiovascular disease, organismal injury, tissue morphology" (score: 50), "Free Radical Scavenging, Cellular Movement, Hematological System Development and Function» (score: 45), « Inflammatory Disease, Respiratory Disease, Carbohydrate Metabolism » (score: 37), « Cellular Development, Hematological System Development and Function, Immune and Lymphatic System Development and Function » (score: 31), « Lipid Metabolism, Small Molecule Biochemistry, Vitamin and Mineral Metabolism » (score: 24), « Cancer, Immune and Lymphatic System Development and Function, Gene Expression » (score: 22).
The most significant networks identified with the best IPA scores are Network 1, 2, 3 and 4 (figures 5A-D). Network 1 and 2 rank top in the IPA score and reflect gene regulations involved in organismal injury, tissue morphology and cellular movement such as IL1, IL1b, ITGB2, PDGFBB, Rac, Ras, MMP9, PPAR, SERPINB2&B5, IL1RAP or CXC3CR1. Network 3 depicts several genes that are involved in cell cycle regulation (cyclines A, E), cell movement (CD44, caspase). Network 4 involves regulations implicated in cellular development and cell-cell interactions (IL8, IL2R, IFN-γ etc). The represented networks are associated with the production and response of cytokine components. More specifically, network 1 and 2 show the importance of IL-1 in the response to wounding. Our results suggest that role for this cytokine in the development of tissue repair and possibly in the onset of angiogenesis. It is known that IL-1β interacts with endothelial cells and induces VEGF and iNOS expression . IL-1β also signals through IRAKs in vascular endothelial cells and induces genes in the VEGF pathway [19,20]. Furthermore, IL-1β induces fibrosis after radiation in experimental models, which indicate a possible interaction of IL-1β with fibroblasts. In addition, IL1β concentrations are elevated in chronic wounds in patients . Thus, IL-1β may regulate fibroblast proliferation and migration during GT formation. Among the other central cytokines are IL-8 and interferon-γ (network 4 and 6, see also additional file 5). The interaction of these cytokines with the vasculature and the stroma is well established [23,24]. Furthermore, IL-8 expression is elevated in chronic wounds but not acute wounds, which supports a role of this cytokine in GT formation . IL-8 is a powerful angiogenic chemokine implicated in vessel formation in cancer, tissue repair and inflammation. It stimulates VEGF expression and autocrine activation of VEGFR2 in vascular endothelial cells . Interferon-γ may interact with stroma fibroblasts and paradoxically negatively impact on angiogenesis . It has been described that the induction of the CXC chemokine interferon-γ-inducible protein-10 (IP10) regulates the reparative response following myocardial infarction . Interferon-γ may, thus, modify the cellular composition of the healing tissue and promotes wound contraction, attenuating adverse remodelling.
The expression of 9 of the differentially regulated genes in the Affymetrix analysis was then investigated by quantitative real-time polymerase chain reaction (qRT-PCR) in wound tissue in comparison to its own control CAM (figure (figure6).6). We chose to analyze these genes because of their expression levels and because not much is known about their involvement in angiogenesis or tissue repair. Transcript levels, as determined by qPCR, were in accordance with the results obtained by the Affymetrix analysis. The respective expression levels were: FABP4 (100.5 fold), ah221 (35.4 fold), HOPX (3.3 fold), TTR (109 fold), CCL20 (37.5 fold), MSLN (5.5 fold), TCF21 (- 3 fold), ITH5 (-8.4 fold), and SPOCK1 (-4.5 fold). Control and wound samples from 9 eggs were analyzed.
We then analyzed expression of the different genes using GenePaint http://www.genepaint.org, ProteinAtlas http://www.proteinatlas.org and Geisha http://www.geisha.arizona.edu. In GenePaint, FABP4 was expressed in the vasculature in the stage 14.5 mouse embryo. Significant staining was seen in blood vessels (figure (figure7A7A to to7C),7C), especially in the kidney (figure (figure7C).7C). For comparison, the pattern of staining as seen in the kidney is identical to that of VE-Cadherin (CDH5), which is a specific marker of blood vessels (figure (figure7D).7D). Furthermore, in the ProteinAtlas, staining of FABP4 was also seen in the vasculature of several organs such as the urinary bladder (figure (figure7E)7E) but it is of note that not all blood vessels are stained for FABP4. There is currently only one publication describing the expression of FABP4 in the vasculature . Our results are in agreement with these findings and indicate a strong expression of FABP4 after wounding. Chemokine ah221 is also significantly expressed after wounding. The human ortholog of ah221 is not yet identified but a similarity of 43% to human CCL3 is observed. CCL3 has leukocyte chemotactic activity and is involved in the recruitment of leucocytes and fibroblasts into neoangiogenic sites such as tumors . CCL3 has also been shown to be produced by endothelial cells and to have an autocrine function . It is likely that ah221 or its human ortholog has similar functions and plays a significant role in wound repair. HOPX is up-regulated 8.41 fold in our transcriptomic analysis. This gene is expressed during cardiac development, in various tumors, such as choriocarcinoma, and it is likely it may have a function in tissue repair as well [31,32].
Transthyretin (TTR) has been found highly up-regulated in our analysis. Transthyretin (TTR) is a plasma protein mostly known for being the transporter of thyroxine and retinol [33,34]. When mutated, TTR is also well-described as the cause of familial amyloid polyneuropathy, a neurodegenerative lethal disorder characterized by systemic deposition of TTR amyloid fibrils. A potential role in tissue repair has not yet been described for this gene but it may have an indirect effect on tissue morphogenesis through retinol. This is supported by our IPA (Additional file 5, Network 6) where 15 genes up-regulated in the wound are known to be directly induced by retinoic acid.
CCL20 (MIP-3 α) is a CC chemokine family member that is highly expressed in our wound transcriptomic analysis. CCL20 activates CCR6 and leads to calcium mobilization and elevated active RhoA, phosphorylated myosin light chain, and F-actin accumulation and stimulation of epithelial cell migration . Furthermore, in rat models of oral wound healing, CCL20 is up-regulated during the peak phase of wound healing . These findings, together with ours, support an important role of this chemokine in wound healing.
Mesothelin (MSLN) is also significantly up-regulated after wounding. This is surprising since MSLN up-regulation is mainly found in mesothelioma, pancreatic, breast and ovarian carcinoma, and tumors of the GI tract http://www.proteinatlas.org. In normal tissue, significant expression is only observed in the fallopian tube http://www.proteinatlas.org. It has been shown that MSLN can interfere with cell cycle regulators by activating ERK kinase and decreasing BIM . Furthermore, an increase in Stat3 activation and cyclin E in MSLN transfected pancreatic tumor cells is observed . As in tumors, MSLN may promote GT formation by promoting proliferation of stromal fibroblasts and vascular cells.
TCF21, which is decreased by 3 fold in our qPCR analysis, has been implicated in kidney and lung organogenesis . TCF21 -/- die in the early perinatal period because of multiple renal defects. TCF21 has also been described as a tumor suppressor gene that undergoes epigenetic modifications . TCF21 has been implicated in myofibroblast differentiation and control of proliferation in mesenchymal progenitor cells . During development, TCF21 is expressed at various sites in the chick embryo such as the pericardium or the allantois (figure (figure7H).7H). During wound healing, TCF21 may be implicated in the regulation of fibroblast proliferation and differentiation in the GT. Inter-α globulin inhibitor 5 (ITIH5) decreased by 8.4 fold in our qPCR analysis, encodes one of the heavy chains of ITI, It is a protease inhibitor associated with the extracellular matrix and contributes to matrix stability by covalent linkage to hyaluronan. Loss of expression has been observed in various human solid tumors . Furthermore, its loss by promoter hyper-methylation is associated with poor prognosis in mammary carcinoma . During wound healing, this molecule is possibly involved in matrix remodeling. SPOCK1/SPARC/Osteonectin/testican-1 is decreased by 4.5 fold in our qPCR analysis. This molecule is a proteoglycan able to inhibit proteases such as MMP2 or Catepsin L . SPOCK1 is also expressed in cancer associated fibroblasts where it reflects EMT . However, SPOCK1 by itself inhibits angiogenesis, enhances tumor stroma formation and prevents fibroblast activation . This may further explain its anti-tumor effect. Down-regulation of SPOCK1, as observed in our study, may contribute to an increase in vessel ingrowth and fibroblast activation during GT formation.
Herbert and coll. have developed a method to assign endothelial cell signatures to genes identified after gene profiling studies . Endothelial cell data filtering evidenced several genes with significant endothelial cell signature (Table (Table4).4). For FABP4, despite the fact that endothelial transcript counts were not high, differences were significant because no counts were observed in the non-endothelial cell pool. As indicated previously, vascular endothelial cell expression of this molecule is very pronounced. For example, FABP4 is found in the vasculature at E14.5 (figure (figure7B7B).
The molecule with the highest endothelial transcript counts was epithelial membrane protein-1 (EMP1). EMP1 was induced in our transcriptomic analysis by about 3.3 fold in the GT after wounding. In GenePaint, expression of EMP1 was seen in vascular cells and possibly in the mesenchyme (figures (figures7F7F and and7G).7G). Expression is much more diffuse than that of FABP4. EMP-1 is however, highly expressed in the endothelial cell lining of the endocard (figure (figure7G).7G). There is only one single publication that reports EMP1 to be present at tight junctions in vascular endothelial cells . This observation fits well with a role of EMP1 in tissue repair, since endothelial cell junctions are remodeled during endothelial cell migration in the GT.
Another gene highly expressed in endothelial cells was MYCT1. MYCT1 is a direct target of c-myc and phenocopies many of the effects of c-myc [49,50]. It has been described to be over-expressed in gastric carcinoma . No publications with regard to endothelial cell expression of this gene have been reported, which is reflected by a surprisingly low Angioscore.
CCL20 exhibited a significant Endofactor and a high Angioscore. It has been reported that endothelial cells in culture express CCL20 upon thrombin stimulation . CCL20 seems to be implicated in endothelial cell-lymphocyte interaction through CCR6 . Lymphatic endothelial cells have been reported to express CCL20 upon induction by lipoteichoic acid (LTA) . However, there have been reports where expression of CCL20 is outside the vasculature, such as in tumor cells . The reason for these differences is not known but genome instability of tumor cells, leading to aberrant CCL20 expression, could be the reason.
From our results, some hypothesis can be formulated of how identified genes may fit into a general scheme of wound repair and GT formation. Tissue repair- is driven by positively and negatively acting factors. Early ingrowth of vessels and fibroblasts are driven through growth factors and cytokines such as IL8, IL1 and PDGFRB. Some additional chemokines including CCL20 or ah221 may also contribute to endothelial cell activation, vessel remodeling during wound healing and fibroblast recruitment. FABP4 may participate in wound repair by promoting endothelial cell proliferation. EMP1 may be involved in the modulation of intercellular adhesion in vessels after endothelial cell activation and participate in the mobility and sprouting of vessels in the GT.
The expression of several matrix or matrix-associated proteins (CL8A1, FBN1, Laminin-α1, FRAS-1, ITIH5 etc.) is likely to be modulated during wound repair and is decreased when significant GT formation had occurred. Down-regulation of SPOCK1, as observed in our study, may contribute to an increase in vessel in growth and fibroblast activation during GT formation. During wound healing, TCF21 may be implicated in the regulation of fibroblast proliferation and differentiation in the GT, and ITIH5 is possibly involved in matrix remodeling. Interferon-γ may interact with stroma fibroblasts and modify the cellular composition of the healing tissue, thus, promoting wound contraction, attenuating adverse effects on remodeling.
It is possible that vitamin A has a role in wound healing  as it is interesting that three genes out of the top four up-regulated in the wounded CAM (from 36 fold to 79 fold up-regulation) from this study are potentially retinol related (FABP4, RBP7 and TTR). It has been shown that TTR (Transthyretin) forms a complex with Retinol Binding Protein (potentially RBP7 here) for transport of retinol around the circulation . In addition, retinol binding proteins have already been shown to be differentially expressed in GT. Retinol is a precursor to retinoic acid, which acts as a steroid hormone, targeting nuclear receptors of genes involved in tissue morphogenesis . It is possible this steroid hormone could be delivered to cells bound to a RBP7-TTR complex and be transported through the cell membrane by FABP4. FABP4 could also deliver retinoic acid to signalling molecules such as Retinoic acid receptors (RARs), Peroxisome Proliferator-Activated Receptors (PPARs) and nuclear response elements. Fatty acid binding proteins have been previously shown to do this . This hypothesis is supported by previous work that found the topical addition of retinoic acid, derived from retinol, to genetically diabetic mice improves wound healing  and that corneal endothelial healing rates increase in the presence of retinoic acid. The participation of other molecules such as Mesothelin is more difficult to envision because of lack of sufficient functional data.
There have been several studies that report transcriptomic profiling in wound tissue in different experimental settings. These include, for example, the transcriptome-wide analysis in excisional murine cutaneous wound inflammation or in chronic ischemic wounds in the pig model . These studies are different to ours because they are performed in an immunocompetent setting and, thus, do not address exclusively the role of the stromal fibroblasts and blood vessels. There has been one study that performed a transcriptome-wide analysis of blood vessels laser captured from human skin and chronic wound-edge tissue . However, in this case contamination by circulating mononuclear cells can also not not be excluded. There have also been transcriptomic profiling studies using models of tissue regeneration such as regeneration of Xenopus laevis hindlimbsor fin regeneration in the medaka fish . These studies are different from ours because myofibroblast invasion does not occur in these models. Our study is complementary to these existing transcriptome profling studies and provides additional informations on the gene networks implicated in wound repair and GT formation.
The CAM wound model has been established to analyze GT formation and the role of invading fibroblasts and blood vessels in this process . It has been found that tissue tension generated by activated fibroblasts or myofibroblasts during wound contraction, mediated and directed translocation of the vasculature. This vasculature can be expanded, secondarily by elongation and vessel enlargement, and finally through splitting and sprouting. We report herein a complete transcriptome analysis of the "wound model" in the chicken CAM, which allowed the identification of gene signatures involved in GT formation and neoangiogenesis. Cytokines and chemokines clearly play a central role as evidenced in our analysis. The limitation of our work is that, contrary to the adult organism, our model is devoid of immunocompetent cells. However, it has been described that MMP-9 positive neutrophile-like inflammatory cells and MMP-13 positive monocyte-like cells are accumulated in growth factor-stimulated CAM. Thus, these cells may also participate, besides blood vessels and stromal fibroblasts, in GT formation after wounding in the CAM.
Another possible limitation is the relevance of our findings for the mammalian setting. Indeed, it is known that some of the regulators identified in the mammalian system that are involved in vascular development are not present in the chick such as the VE-statins . However, this is a general problem for every model organism including murin models. As an example, CXCL4L1 is only expressed in man, mouse and chimpanzee. CXCL4L1 is a potent angio-inhibitory chemokine that has potent inhibitory activity across species. Furthermore, our laboratory has performed molecular profiling studies using human xenograft tissue in the chick CAM and identified gene regulatory mechanisms relevant for the mammalian setting [12,70]. Thus, we believe, that our results are of importance to the general understanding of GT formation and tissue repair.
Brown Leghorn eggs were cultured at 38°C for 3 days. The shells were then cracked and the contents transferred to 10 cm cell-culture Petri dishes. Embryo culture was continued for a another 7 days, when two injuries to the CAM were inflicted by parallel scalpel superficial cuts of 1 cm area, with a subsequent scarping off the epithelium of the injured chorioallantoic membrane. The wound area was then covered with 1.5 cm square nylon grid and after 6 days, the CAM tissue (control and wound) were excised and processed for subsequent analysis.
Total RNA from cells or snap-frozen tissues were extracted by using RNeasy mini kit (Qiagen, Courtaboeuf, France). RNA quality and quantity were assessed by agarose gel electrophoresis and optical density measurement. First strand cDNA was prepared from 1 μg of total RNA with Quantitect Reverse Transcription kit (Qiagen). For all samples, a negative control was realized with mRNA without reverse transcriptase in the reaction mixture.
RNA was isolated from control and "wound" CAM 6 days after injury. Three eggs were used for independent transcriptomic profiling. It is important to note that from each egg, unwounded control CAM and wound were analyzed. RNA was isolated according to standard procedures and hybridized to Affymetrix chicken GeneChips using the Affymetrix standard protocol (Affymetrix UK Ltd, High Wycombe, UK).
RNAs hybridized to Affymetrix chicken GeneChips using the Affymetrix standard protocol (Affymetrix UK Ltd, High Wycombe, UK). The chicken GeneChip covers 32773 transcripts, corresponding to > 28000 chicken genes, and has a probe set oligonucleotide length of 25 and a detection sensitivity of 1:100000 http://www.affymetrix.com. Data were analyzed with the GCOS 1.2 software (Affymetrix), using the default analysis settings; global scaling as first normalization method, with a trimmed mean target intensity value (TGT) of each array arbitrarily set to 100.
Gene expression profiles were identified using two-class Significance Analysis of Microarrays (SAM) method  (http://www-stat.stanford.edu/~tibs/SAM/, which utilizes a Wilcoxon-test statistic and sample-label permutation to evaluate statistical significance between sample groups. SAM provides mean fold change values (FC) (mean fold-change > 2) and a false discovery rate (FDR) confidence percentage based on data permutation (n = 200). The False Discovery Rate (FDR), an estimate of the fraction of selective genes, was kept below 5% in all statistical analyses.
Data analysis was done using the Gene ontology database included in the statistical environment R library GOstat http://www.geneontology.org and Ingenuity Pathway Analysis (IPA) (Ingenuity Systems, Redwood City, CA 94063) software. The functional clustering was performed using the method implemented on the DAVID website http://david.abcc.ncifcrf.gov/.
Annotation of genes was performed using NetAffx http://www.affymetrix.com. The microarray data files have be submitted to the US National Center for Biotechnology Information, Gene Expression Omnibus (GEO), and released in May 6th, 2010 (GEO accession number: GSE21679).
Real-time PCR was carried out in an Mx3000P thermocycler (Stratagene, La Jolla, CA) by using SYBR Green dye (ABgene, Courtaboeuf, France). Chicken-specific primers were designed and respectively evaluated for amplification efficiency using total RNA isolated from a total chicken embryo on embryonic day 5. Only primers pairs with amplification efficacy between 90 and 100% were used. The PCR specificity was verified by dissociation curve analysis and agarose gel electrophoresis of the amplification product.
The primer sequences are: gah221 forward CTGGCCCTCTGCTCCTCA and reverse GGACGGGACGTTGAACATAG, gCCL20 forward CGGAAGGTCATTAAGGGC and reverse AAACCATATCACATTGACATCCTC, gFABP4 forward AGACTGCTACCTGGCCTGAC and reverse GCCATCTTCCTGGTAGCAAA, gHOPX forward GCAGTCACGCTGGCTATAAA and reverse CCATTTCTCCTGGATGGTG, gITIH5 forward TCTTGTTGCCCTTGGAAATC and reverse TTCTTTCCTCCCACCTCCTT, gMsln forward AAAATGAACAGGCTGCTGCT and reverse TCAGGCTGTTGGGGTCTATC, gSPOCK1 forward AAAGCAGGGGACCGTTAGTT and reverse TTCCAAATCATCCAGCAACA, gTCF21 forward CCATCCAGTCAACCTGACCT and reverse AGCGGTTTGTGTTCACCACT, gTTR forward TTGATTCCAAATGCCCTCTC and reverse TAGCAAAGTCCTGCCAGGTT and the house keeping gene gHNRPH1 forward GCTGTGTCTGCCACGAGTTA and reverse GCTTTCGGCTGAGAGACAAT.
To concentrate on genes of importance to human pathology and physiology, human orthologs of the chicken genes present on the Affymetrix chicken GeneChips were predicted using a Reciprocal Best Hit (RBH) approach. In this work, both human and chicken Refseq nucleotide and protein sequences were downloaded from the NCBI on 30th January 2008 ftp://ftp.ncbi.nih.gov/refseq. Likewise, cDNA accession numbers of sequences used to design the microarray probes were extracted from Affymetrix chicken chip file "Chicken.na22.annot.csv" http://www.affymetrix.com. Each cDNA sequence, depending on the source of the probe design, was downloaded from ENSEMBL or the NCBI, and, as most of the cDNA sources were Expressed Sequence Tags, full length chicken mRNAs were sought by BLAST searching each of them against the Refseq chicken nucleotide database. The resulting matches were ranked as good, reasonable or bad, depending on the alignment quality (Good: sequence alignment > 100 base pairs with a percent identity > = 96%. Reasonable: a sequence alignment > 100 bases and a percent identity > 90% and < = 95%. Bad: all other hits). The full length chicken mRNA sequences ranked good and reasonable were then used in a RBH analysis.
To identify which chicken genes could have an endothelial cell expression signature, the human orthologs were compared with the results from a novel in-silico bioinformatics screen, where an accurate EST-to-gene assignment and a new likelihood ratio statistic were used to find genes preferentially expressed in endothelial cells using cDNA library analyses (see Herbert et al. 2008 for a full description). The intersection of the comparison were endothelial genes and only those genes  with a q-value < = 0.01 were considered. An "Endofactor" describes how significant a gene was endothelial as found with the q-value.
For all the genes found differentially expressed on the chicken chip, a literature search of the human orthologs of chicken genes were carried out to find those having literature relating to relevant pathologies and physiologies. To accomplish this, Perl scripts were written that searched article abstracts for the following keywords typical of angiogenic research. They were "angiogenic", "angiogenesis", "neovascularis(z)ation", "vasculogenesis", "vascular", "VEGF", "hypoxia" and "endoth" (for endothelial or endothelium).
In situ hybridization and immunostaining data were retrieved from the Genepaint data base (Max-Planck-Institute of Biophysical Chemistry, Dept. Genes and Behavior, 37077 Goettingen
Germany; http://www.genepaint.org; figure figure7A7A to 7D, F and and7G),7G), Proteinatlas (AlbaNova University Center at the Royal Institute of Technology, Stockholm, Sweden, the Rudbeck Laboratory, Uppsala University, Uppsala, Sweden and Lab Surgpath, Mumbai, India; http://www.proteinatlas.org; figure figure7E)7E) and the GEISHA database (University of Arizona, Tucson, AZ 8572; http://geisha.arizona.edu; figure figure7F7F).
Wound areas were fixed for at least 24 hours in a Zn-fixative . Tissues were embedded in paraplast, and section of 10 μm with a MICROM HM325 were performed and placed onto Super Frost slides. Dewaxed slides were either stained with Weigert's hematoxylin and eosin or incubated 40 min at 95°C in a citrate buffer for the antigen recovery before immunohistochemistry. These slides were then fixed with 4% paraformaldehyde, permeabilized with Triton-X100 (0.1%), saturated with 5% BSA in PBS (pH 7.4) and incubated with the primary antibody over-night at 4°C (anti-human Smooth Muscle actin, DAKO, IR611) This antibody also recognizes very well the chicken protein. Secondary fluorescent anti-mouse antibody was from Molecular Probes (used at1:1,000, Invitrogen). Chick blood vessels were visualized by using fluorescein-coupled Sambucus nigra lectin (SNA-1 lectin, 1:100, Vector Laboratories). Cell nuclei were visualized by DAPI (Invitrogen). Microphotographs were taken with a Nikon eclipse E600 microscope equipped with a digital camera Nikon DS-Ri1.
ATOH8: Atonal hololog 8; CCL: chemokine (CC motif) ligand; CXCL: chemokine (C-X-C motif) ligand; CAM: chorioallantoic membrane; C1q: Complement factor 1q; CRISP3: cysteine-riche secretory protein; EMP1: epithelial membrane protein-1; FABP4: fatty acid binding protein-4; GT: granulation tissue; HOPX: Homeobox PX; ITIH5: inter-a (globulin) inhibitor H5; LAMA1: laminin a1; LTA: lymphotoxin alpha; MSLN: mesothelin; MYCT1: myc target-1; NCF2: neutrophile cytosolic factor; PLXNA2: plexin A2; RBH: reciprocal Best Hit; RTBP7: retinol binding protein 7; SPP1: osteopontin; SPOCK1: sparc/osteonectin, cwcv and kazal-like domains proteoglycan (testican); C1QTNF1: TNF related factor-1; TCF21: transcription factor 21; TTR: transthyretin; TGT: trimmed mean target intensity value.
WK did CAM experiments, FS did expression analysis, PA and FF did bioinformatics analysis, JH and RB did the endothelial and Angioscore screens, AB supervised the work and wrote the manuscript. All authors read and approved the final manuscript.
WK and FS have been post-doctoral fellows in the AB laboratory, PA is PhD student in the FF laboratory, JH is bioinformatics officer in the RB laboratory, RB is Professor in the division of Immunity and Infection at the Institute for Biomedical Research at the University of Birmingham Medical School (UK), FF is senior lecturer at the School of Biosciences at the University of Birmingham (UK), AB is Professor in cell and molecular biology at the university Bordeaux and director of the molecular angiogenesis laboratory of the National Institute for Health and Medical Research (INSERM, France).
(additional Table S1): Clustering of genes with similar function according to Gene Ontology terms.
(additional Table S2): List of all up-regulated genes (cut off > 2).
(additional Table S3): List of all down-regulated genes (cut off > 0.1).
(additional Table S4): Ingenuity pathway analysis. Regulated genes were analyzed by Ingenuity pathway analysis (IPA) in order to determine functional categories. Functional categories are listed in the tables that depict categories for up-regulated (A) or down-regulated genes (B). As it can been seen in the tables, categories related to tissue morphological processes have the highest Ingenuity scores.
(additional Figure S1A-J): Representation of networks from the IPA for up-(A-H) or down-regulated genes (I-J).
This work was supported by grant from the Agence Nationale de la Recherche (ANR, EGC) to AB. WK was supported by a fellowship from the Lefoulon-Delalande Foundation.