Gene identification and analysis of transcripts differentially regulated in fracture healing by EST sequencing in the domestic sheep
BMC Genomics volume 7, Article number: 172 (2006)
The sheep is an important model animal for testing novel fracture treatments and other medical applications. Despite these medical uses and the well known economic and cultural importance of the sheep, relatively little research has been performed into sheep genetics, and DNA sequences are available for only a small number of sheep genes.
In this work we have sequenced over 47 thousand expressed sequence tags (ESTs) from libraries developed from healing bone in a sheep model of fracture healing. These ESTs were clustered with the previously available 10 thousand sheep ESTs to a total of 19087 contigs with an average length of 603 nucleotides. We used the newly identified sequences to develop RT-PCR assays for 78 sheep genes and measured differential expression during the course of fracture healing between days 7 and 42 postfracture. All genes showed significant shifts at one or more time points. 23 of the genes were differentially expressed between postfracture days 7 and 10, which could reflect an important role for these genes for the initiation of osteogenesis.
The sequences we have identified in this work are a valuable resource for future studies on musculoskeletal healing and regeneration using sheep and represent an important head-start for genomic sequencing projects for Ovis aries, with partial or complete sequences being made available for over 5,800 previously unsequenced sheep genes.
The accumulation of genome-wide sequence data for model organisms is of obvious benefit for many kinds of scientific analysis. Although the sheep is of great economic importance and has been used as a model organism to study such diverse topics as fracture healing , ovarian function , hypernatremic hypertension , and pulmonary function , relatively little sequence information is available for this organism, and DNA sequences are available for only a small number of sheep genes.
2030 loci, 543 of which are designated as genes, are currently listed in the ArkDB sheep database . In UniGene  Build 7 for Ovis aries, a total of 5,730 EST and mRNA sequences form a total of 1714 clusters, and the EBML database  currently contains 10497 ESTs and 2391 mRNA entries for sheep.
Fractures, and especially osteoporotic fractures in the elderly, represent a major public health problem. As populations age, the costs of treating osteoporotic fractures, currently billions of dollars a year in the United States alone, are predicted to dramatically increase . Consequently, a better understanding of the genetic processes underlying fracture healing is of great medical relevance. The process of fracture healing is a regenerative process that in many ways is similar to bone development. The fracture repair process goes through many stages, beginning with an inflammatory response and the recruitment and proliferation of undifferentiated mesenchymal cells. Mesenchymal cells from the periosteum near the fracture site differentiate into osteoblasts and induce intramembranous ossification resulting in the formation of new bone. Simultaneously, cartilage is formed at the fracture site and undergoes enchondral ossification. The newly produced bone then undergoes extensive remodeling until the original shape and structure of the bone is restored. Many of the molecular signals and processes involved in fracture healing closely resemble those involved in bone development [9–11].
In the present work, we have extracted 47 thousand ESTs from libraries developed from healing bone in a sheep model of fracture healing. We clustered our sequences together with previously identified ESTs, and identified 19087 contigs, the majority of which represented previously unsequenced sheep genes. These newly identified sequences were used to perform RT-PCR on 78 genes from the postfracture day 7, 10, 14, and 42 libraries. We confirmed differential expression for genes previously known to be involved in fracture healing and demonstrated differential expression for numerous genes not previously known to be involved in fracture healing.
ESTs and assembly into contigs and creation of an Ovis aries sequence database
In addition to the 47209 sheep ESTs generated from three libraries in the current project (Table 1), we extracted 10497 ESTs and 2391 mRNA entries from the EMBL database (release 82, Sept 2005). 344 sequences were rejected because less than 40 bp remained unmasked following quality control procedures. The lengths of the ESTs ranged from 40 to 928 with an average length of 555 nucleotides (Fig. 1a). All bases beyond 750 bp were masked. Using our sequences as well as the previously available sequences, 19087 contigs were assembled with an average length of 603 nucleotides (Fig. 1b). 397 contigs were shorter than 100 nucleotides and were excluded from further analysis. On average, each contig was assembled from 3.1 ESTs from all libraries, with a range of 0 to 2392 ESTs per contig (some contigs consisted only of previously available mRNA data and contained no ESTs). If the analysis is restricted to ESTs from the libraries described in this work, then each contig contained an average of 2.5 ESTs (Range: 0 – 2387). The contigs comprised of the largest number of ESTs were identified as collagen Iα 1 (2387 ESTs), collagen Iα 2 (1094 ESTs), and osteonectin (1086 ESTs). A total of 8965 contigs contained only a single EST clone.
We used the GeneNest package [12, 13] to analyze all available mRNA and EST sequences and to present the sequences in a user-friendly internet site. The main features of the GeneNest Sheep database are summarized in Fig. 2.
BLAST analysis of Ovis aries contigs to assign putative gene identities based on homology
BLAST searches were used to attempt to assign putative identities to the contig sequences identified in this study. Analysis against the NCBI RefSeq protein sequence database was performed, and 13149 contigs demonstrated a best e-value of 0.01 or lower (Fig. 3a). We used an arbitrary e-value cutoff of 10-25 to identify significant homology to a known gene and allow a probable assignment of a sheep contig. With this cutoff, a total of 10387 contigs demonstrated significant homology to at least one other protein in the RefSeq protein database and could be assigned a putative identity, consisting of a total of 8789 distinct proteins. In cases where hits to a proteins from several species were identified, the most significant match was used to assign a putative identity. 292 contigs blasted to an Ovis aries RefSeq. A total of 2863 contigs (15%) had e-values of less than 10-100 and are thus likely to represent true orthologs. Of the 19087 contigs, 11591 were assembled from ESTs only from the current project. 6577 were assigned to 5839 distinct RefSeq sequences by virtue of blastx e-values of 10-25 or less. These sequences represent previously unidentified sheep genes. The sequences for these genes are available in FASTA format with the online supplemental material as [Additional file 1]
Gene-Ontology analysis of the sheep fracture libraries
GO analysis was performed by assigning annotations to the sheep contigs from the MGI annotations of putative murine homologs (e-value threshold 10-25) in order to get an overview of the subset of genes expressed in healing fracture tissue of the sheep. Fig. 3b provides an overview of the results of this BLAST analysis. Putative homologs were found for 5637 distinct sheep genes, some of which were defined from multiple contig sequences. MGI annotations were available for 4600 of the 5637 genes. The results of the GO analysis using these annotations are summarized in Fig. 4. Many of the most frequently annotated terms in the subontology biological process, such as development, signal transduction, transcription, biosynthesis, and phosphate metabolism, reflect processes known to be involved in fracture healing and bone development. The distribution of annotations to terms in the subontology cellular component shows that the gene products corresponding to the ESTs identified in this work are active in many different cellular and extracellular locations.
A comparison of the sheep fracture libraries with normal human bone libraries
During development and in related processes such as fracture healing, genes exhibit complex patterns of spatial and temporal regulation. The identification of genes that are differentially expressed during fracture healing therefore represents an important initial investigation to discover genes that are essential for fracture healing.
In this work, we generated two EST libraries from early fracture healing stages (PF7 and PF10) with the goal of identifying genes important in the earliest stages of fracture healing. Following fracture, initial hemorrhage into the fractured bone is followed by hemostasis and a phase of inflammation, which lasts up to two weeks in the sheep (Fig 5). Intramembranous ossification begins around day 10. Therefore, we reasoned that genes that are overexpressed in the PF7 and PF10 libraries could represent genes with potential importance for the initiation of new bone formation. To estimate gene expression levels in the two libraries, we compared EST sequence counts using a two-tailed Fisher exact test. This analysis revealed statistically significant differential expression for only three genes; tartrate resistant acid phosphatase 5 (ACP5) increased by a fold change of 10.60 (5 to 53 ESTs), cytochrome c oxidase subunit I (COX1) increased by a fold change of 2.1 (87 to 178 ESTs), and perlecan decreased by a fold change of 3.9 (31 to 8 ESTs). We did not perform such a comparison using the third EST library, because the latter was derived from a mix of time points using PCR amplification of cDNA inserts in order to get a broader overview of the genes involved in the many stages of fracture healing in sheep.
RNA samples from intact sheep bone were not available. Therefore, we chose to compare the EST distribution in the PF7 and PF10 libraries to that of two resting human bone EST libraries. We reasoned that genes that were significantly more highly expressed in the sheep libraries could represent candidate fracture-healing related genes that might warrant closer investigation. A total of 39 genes with Benjamini-Hochberg-corrected p-values less than an arbitrary threshold of 0.05 were identified as putatively differentially expressed. Table 2 gives details on the genes and number of ESTs found (see Methods for abbreviations). Links to corresponding contigs in the GeneNest database are given in [Additional file 2].
It is striking that 24 of the 39 genes (62%) had previously been known to be involved with developmental processes in some way (references are given in the table). Eleven of the 24 genes were previously known to be involved in bone development or fracture healing. Given that fracture healing is thought to recapitulate bone development on a molecular level , these findings are suggestive that our EST analysis captured some essential aspects of fracture healing.
Expression analysis of candidate fracture-healing related genes
We then reasoned that differential expression during the course of fracture healing could be another kind of indicator for genes with potential importance in fracture healing. Therefore, we performed RT-PCR expression analysis using samples from postfracture day 7, 10, 14, and 42. The RNA source for days 7 and 10 was identical to that used for the PF7 and PF10 EST libraries, and the source for day 14 was one of the RNA samples used to construct the pooled EST library. The day 42 RNA sample was not used for any of the EST libraries. Samples from normal sheep bone were not available.
PCR primers were designed for 30 of the genes in Table 2 and for an additional 48 candidate genes (Table 3), that were chosen for biological interest or because of tendencies toward overexpression that did not reach statistical significance.
Almost all genes showed a large drop in expression at the 42-day time point. Early stage callus tissue is a metabolically highly active tissue. In comparison, healthy bone is less metabolically active, which is not surprising given that cells make up only around 2% of the total mass of bone. It has previously been noted that the expression of many genes related to fracture healing drops to very low or undetectable levels at 6 weeks after fracture in a rat model . In our experiments, 77 of the 78 genes (99%) tested by ANOVA over all four timepoints were found to be significantly differentially expressed, a proportion that dropped to 96% if only the first three time points were considered, and 31% for a comparison between only the first two time points. [Additional file 2] presents the normalized expression values determined in these experiments.
We then performed hierarchical clustering to produce a graphic representation of the the expression patterns for the genes investigated in RT-PCR. We used the gap-stat method  to estimate the number of clusters in the data (Fig. 6). The estimated number of clusters is two. However, the gap function rises again with peaks at 4 and 6 clusters, suggesting that there are less well defined groups of four and of six clusters in the data. For the purposes of discussion and visualization, we chose to partition the genes into six clusters (Fig. 7). Clusters A, B, E and F demonstrated high frequencies of genes of specific functional classes (see below).
In the present work, we have performed large-scale EST sequencing and analysis in order to generate sufficient sequence information to study fracture healing at a genomic level in the sheep for the first time. The results of our study are of interest not only in the context of fracture healing and bone development, but will also provide useful information for all fields of research that use the sheep as a model organism. The ESTs generated in the present work represent an increase in the number of available ESTs by a factor of over 4.5, and provide at least partial sequences for over 5800 sheep genes never before sequenced. The sequences can be downloaded from the online supplemental site and are also available in the GeneNest Sheep database.
Genomic approaches to the sheep have great potential to accelerate progress in research areas using the sheep as a model organism. For instance, the first study using an ovine cDNA microarray prepared from unpublished EST libraries was very recently published. A number of genes differentially expressed in sheep resistant and susceptible to gastrointestinal nematode infection were identified .
Although we now have initial concepts of the molecular processes involved in fracture healing, many questions remain to be answered before a comprehensive picture of the genetic programs involved in fracture healing can emerge. An in-depth understanding of the genes involved in fracture healing and the regulatory mechanisms that control them will be prerequisite for developing optimal treatments for disturbances of fracture healing and will have relevance to other disorders such as osteoporosis. At present, we still do not have a comprehensive list of key genes and gene products involved in bone development and fracture healing. In this study, we have performed the first large-scale EST sequencing project using a model for fracture healing and have identified numerous differentially regulated genes not previously known to be involved in fracture healing. EST sequencing gives a nearly unbiased view of highly expressed sequences in a tissue (assuming that the process of producing the ESTs does not introduce a significant amount of bias), and therefore offers some advantages over candidate gene approaches and even over microarray-based approaches that may not cover the entire genome.
In this work, we have generated a profile of postfracture day 7 and day 10 genetic expression using EST libraries produced by automated plasmid purification. No normalization procedures were applied to select ESTs for sequencing. Therefore, we believe that the distribution of ESTs in these libraries gives a nearly unbiased view of genetic expression at these stages of fracture healing. In addition, we generated a third EST library from a pool of later stages of fracture healing and used RT-PCR to examine differential expression of 30 of the 39 genes found to be enriched in the PF7 and PF10 libraries and 48 other genes (Tables 2 and 3).
Previous genomic approaches to fracture healing have involved analysis of a rat femur fracture model based mainly on microarray expression analysis, but different experimental and statistical approaches were used [11, 17, 18]. These studies have collectively identified thousands of genes that are differentially expressed in fracture healing. The genes were grouped by the authors into a number of functional classes, which overlap to a large extent with the classes identified in the present investigation, including extracellular matrix proteins, ribosomal proteins, resorption/remodeling, transcriptional regulation and signal transduction. However, the specific genes belonging to the functional classes are often different among the three studies cited above as well as in comparison to our results. For instance, Hadjiargyrou and coworkers (2002) identified 20 differentially expressed ribosomal genes; in the present study, 16 differentially expressed ribosomal genes were identified. With the exception of ribosomal protein SA, there is no overlap between the two lists. This suggests that we are just beginning to identify to the full complement of genes involved in fracture healing.
Genes differentially expressed in fracture healing
Many of the differentially expressed genes we identified were previously known to be involved in various biological processes involved in bone development or fracture healing (see Tables 2 and 3 for references), suggesting that our approach has captured some essential aspects of fracture healing. Additionally, we have demonstrated differential expression for many genes not previously known to be involved in fracture healing.
The differentially expressed genes can be divided into several different classes, as will be discussed below. Three of these classes, extracellular matrix proteins, resorption/remodeling/inflammation, and transcriptional regulation/signal transduction, have been well studied in the context of fracture healing, and our results have identified both previously known and novel differentially regulated genes with these functional roles.
In the following, we describe three other classes with numerous differentially regulated genes, angiogenesis, free-radical control, and ribosomal genes. These functional classes have received less attention to date, although they are known to be involved in the process of fracture healing. Fracture-related differential expression had not previously been shown for the great majority of genes in these classes (cf. Tables 2 and 3). These observations thus represent starting points for future studies on the role of processes such as angiogenesis and free-radical control in fracture healing.
The single largest group of putatively upregulated genes in our analysis was made up of ribosomal proteins. There are over 75 ribosomal proteins , which are not necessarily expressed in concert . We found 9 genes coding for ribosomal proteins to be significantly overexpressed in our analysis of EST distribution compared to normal human bone. All of these genes, as well as seven additional genes coding for ribosomal proteins (Table 3) showed differential expression in the RT-PCR analysis, with expression levels peaking at postfracture days 7 and 10 (Fig. 7). We additionally identified four translation factors as being overexpressed in our EST analysis (Table 2). RT-PCR for two of these genes confirmed differential expression. To the best of our knowledge, a specific role in fracture healing has not been identified for all but two of these genes to date, although some of them are known to be differentially expressed in certain developmental processes (Tables 2 and 3), and a number of other ribosomal genes are known to be differentially expressed in fracture healing .
Cluster B (Fig. 7) contained 36 genes. Interestingly, 13 of the 15 ribosomal genes analyzed by RT-PCR, and both of the translation factors analyzed by RT-PCR were in this cluster, whose expression levels peaked during the early phases of fracture healing at days 7 and 10 postfracture, suggesting that there is an especially high level of protein synthesis at these stages.
Extracellular matrix proteins
Several matrix proteins were found to be differentially regulated. The gene with the highest overall expression, Collagen Iα 1, is the main protein component of bone, and the fact that it is one of the most highly upregulated genes is not surprising. Two other collagens and osteonectin were also found to be upregulated by EST analysis (Table 2). We investigated these genes and four other collagen genes and two proteoglycans by RT-PCR and showed differential expression for all genes. In general, the expression of these matrix genes was highest at postfracture day 7 or 10 but was still relatively high at day 14. IBSP is an osteoblast marker that is a major noncollagenous structural protein of the bone matrix. IBSP showed an expression pattern different from that of the other matrix genes, with significant expression beginning at day 10, peaking at day 14, and continuing until day 42.
The lysyl-tRNA synthetase KARS has no known involvement in fracture healing. Given the important role of lysine crosslinks in the biosynthesis of collagen, it is interesting that lysyl-tRNA synthetase was the only tRNA synthetase gene that was overexpressed in our data.
SERPINH1 was also found to be significantly overexpressed in the PF7 and PF10 ESTs. It is a molecular chaperone involved in the maturation of collagen molecules . This is presumably important for bone development because SERPINH1 is localized to regions of type I collagen production in developing murine femurs and tibiae [22, 23].
It has long been known that vascular invasion is necessary for bone differentiation , and angiogenesis is a key process for fracture healing . A number of genes related to angiogenesis have previously been shown to be upregulated during bone repair, such as vascular endothelial growth factor .
Six of the eight genes in our dataset with prominent roles in angiogenesis (CSRP1, COL4A1, COL4A2, Perlecan, HDGF, Endoglin) were in cluster A. The smooth muscle marker CSRP1  and the TGFβ-family auxiliary receptor endoglin are required for angiogenesis and heart development [30, 31]. HDGF is a highly expressed vascular endothelial cell protein .
The genes of this cluster show especially strong expression on day 7 postfracture with lesser amounts of expression on days 10 and 14, which is not surprising given that new blood vessel formation is a prerequisite for the formation of new bone from callus tissue. The other two genes, COL3A1 and TAGLN2, were in cluster B, which like cluster A is characterized by high early expression (PF7 and PF10), but shows a stronger drop in expression at the 14 day timepoint.
Resorption, remodeling and inflammation
Several genes that probably play a role in the inflammatory response characteristic of early stages of fracture healing (FCGRT, Basigin, HLA-A, and CD74) were found to be overexpressed in the PF7 and PF10 ESTs (Table 2).
Resorption of cartilage and bone mediated by osteoclasts is an important part of the remodeling processes involved in new bone formation in fracture healing. The main enzymes involved are thought to be the cathepsins and the matrix metalloproteinases (MMPs), which act in a concerted fashion. Although not all details are clear, our current understanding of this process suggests that osteoclasts create an acidic environment in their resorption lacunae, which results in dissolution of the mineral, and then secrete cysteine proteinases that are active at a low pH, to initiate proteolysis of the proteinaceous bone matrix. Finally, MMPs exert their activity, once the pH has increased sufficiently . Our data has demonstrated differential expression for several enzymes from both classes.
Interestingly, four of the five genes in RT-PCR cluster E encode enzymes with roles in osteogenesis (ACP5, CTSK, MMP9, and MMP13). The first three genes are expressed by osteoclasts. The expression of genes in this cluster was minimal at day 7 but peaked at day 10 postfracture, suggesting that these enzymes might be particularly important for the processes involved in the bone formation that becomes visible by day 10.
We showed differential expression for MMP2, MMP9, MMP13, MMP14, and MMP19. All of these MMPs had previously known roles in bone development except MMP19, providing further evidence for the similarity of bone development and fracture healing. MMP2, MMP9, and MMP14 are differentially regulated in scarless fetal wound healing  and during osteogenesis , which is interesting because bone is one of the few adult tissues that can heal without scarring. MMP19 previously had no known role in skeletal development or fracture healing, but it is highly expressed in dermal wounds, suggesting a role in wound repair . The highest expression of MMP19 was observed in the day 7 and 10 stages.
Cluster F contained only two genes, cystatin C and TPM2. These genes were the only ones with an expression profile that peaked at the last time point (42 days postfracture). Cystatin C (CST3) can inhibit bone resorption and osteoclast formation . Osteoclast activity begins early and increases over the course of fracture healing , and one may speculate that CST3 may be important for regulation of osteoclasts at this time point.
Control of free radicals
Vascular invasion of ossified cartilage during enchondral ossification is associated with breakdown of tissue by the release of lytic enzymes by invading cells such as macrophages and endothelial cells. This process is likely to involve the production of reactive oxygen species (ROS) . At present, little is known about the mechanisms by which unwanted tissue damage by ROS is controlled during fracture healing. The glutathione peroxidases (GPX) are enzymes that catalyze the reduction of hydrogen peroxide, organic hydroperoxide, and lipid peroxides by reduced glutathione. They are thus involved in the protection of cells against oxidative damage. Although none of the four GPX genes we tested in RT-PCR were individually statistically significantly overexpressed in the EST analysis, there were no ESTs for any of these GPX genes in normal bone and a total of 85 ESTs for the four GPX genes in the PF7 and PF10 EST libraries. We showed significant differential expression patterns for each of these genes by RT-PCR, although they did not show a uniform expression pattern (Fig. 7).
HMOX1 was also found to be significantly overexpressed in the PF7 and PF10 EST libraries, as has been previously described for postfracture day 3 in a rat femur fracture model . HMOX1 is an inducible heme oxygenase that could conceivably be involved in heme catabolism during resorption of the initial fracture hematoma, although other roles in development have been postulated for this gene .
Transcription regulation and signal transduction
Two membrane-bound proteins with roles in signal transduction were found to be significantly overexpressed in the PF7 and PF10 libraries. ITM2C is highly expressed during chondro-osteogenic differentiation , but its specific function and its role in fracture healing remain unknown. The G protein GNAI2 was also found to be overexpressed, but no clues are available as to its particular role in fracture healing.
Signal transduction and modification of transcriptional programs by transcription factors represent important phenomena in development, and presumably are essential for fracture healing. Nevertheless, the EST counts for the transcription factors examined in this study were too low to reach statistical significance, even though in many cases ESTs for transcription factors were found only in the PF7 and PF10 libraries but not in the human normal bone libraries. For instance, we identified only one EST each for the transcription factor runx-2 and osterix, which are both key players in skeletogenesis [43, 44].
We performed RT-PCR analysis on a number of transcription factors with relatively large numbers of ESTs in our libraries (Table 3). We were able to demonstrate differential expression with the highest expression at days 7 and 10, for JunB, which among other roles is a positive regulator controlling primarily osteoblast as well as osteoclast activity . MORF4L1, which may play a role in mechanotransduction , and TCF4, which may mediate Wnt signaling during limb development  also showed highest expression levels in the first two time points. Additionally, we showed differential expression for NSEP1, PPIB, and RARA, for which to the best of our knowledge no previous information about differential expression during bone development or fracture healing was available.
In this work we present approximately 47,000 ESTs from the domestic sheep. This represents an important resource to researchers interested in the sheep as a model organism and greatly expands the number of gene sequences available for Ovis aries .The sequences are freely available in GeneNest and the public sequence databases. We have used the sequences we generated as a tool to investigate the early stages of fracture healing and have identified numerous genes that are differentially regulated in early callus tissues, many of which were not previously known to have any specific relation to fracture healing or bone development. These genes will be valuable leads for further research into the molecular mechanisms of fracture healing.
All animal experiments were carried out according to the policies and principles established by the Animal Welfare Act, the NIH Guide for Care and Use of Laboratory Animals and the national animal welfare guidelines. The study was approved by the local legal representative (Landesamt für Arbeitsschutz, Gesundheitsschutz und technische Sicherheit, Berlin: G 0224/01 and G 0172/04).
The experimental system has been described in detail elsewhere [48, 49, 38]. Briefly, twelve two-year old female Merino mix sheep with a mean weight of 66 ± 12 kg underwent a standardized midshaft osteotomy of the right tibia under general anesthesia. The osteotomy was distracted to a gap of 3 mm and stabilized with a monolateral external fixator. The fixator consisted of 6 Schanz screws (5 mm, 3 inserted proximally, 3 inserted distally of the osteotomy) and 2 steel tubes (10 mm) and was mounted medially to stabilize the defect.
The osteotomy wound was sutured and the leg, together with the external fixator, was covered with a tube bandage. The animals received an analgetic (Finadyne, Essex, Germany) for 3 days postoperatively. Daily pin care involved cleaning of the insertions of the Schanz screws with ethacridin lactate (Rivanol, 5%, Chinosol, Germany).
After removal of the fixator under general anesthesia the callus tissue was removed and sheep were sacrificed. The tissue was frozen in liquid nitrogen and stored at -80°C until further use.
Isolation of RNA and library construction
For generation of ESTs, three separate libraries were constructed. The first library (PF7) was isolated from four animals 7 days post-fracture, the second library (PF10) was isolated from three animals 10 days post-fracture, and the third library was isolated from a pool of one animal each at 2 weeks, 3 weeks, and 4.5 weeks post-fracture.
PF7 consisted of 3.6 × 105 clones with an average insert length of 1.6 kb. PF10 consisted of 3.4 × 105 clones with an average insert length of 1.4 kb. The pooled library consisted of 4.5 × 106 clones with an average insert size of 1.4 kb. Average insert lengths were determined by screening at least 300 inserts per library.
For isolation of RNA the tissue was ground in a liquid-nitrogen-cooled stainless steel mortar. RNA was isolated with peqGOLD TriFast (peqLab) and subsequently purified with RNeasy Mini columns (Qiagen) followed by poly(A)-RNA isolation using the polyATtract mRNA Isolation System III (Promega).
For construction of cDNA libraries the CloneMiner cDNA library construction Kit (Invitrogen) was used according to the manufacturer's protocol with the following modifications. All precipitations were carried out using PelletPaint Co-Precipitant (Novagen) instead of glycogen. The first strand reaction was carried out in five separate reactions with 1 μ g poly(A)-RNA each. Four reactions were precipitated and pooled with the fifth for a single second strand reaction. Glass fiber filters for radioactive determination of the cDNA yield were processed separately. Clones for further analysis were arrayed into 384-well-microtiter plates by robotic colony picking.
Sequencing and analysis of ESTs
Template DNA for sequencing was prepared by PCR amplification of cDNA inserts with M13 primers or by automated plasmid purification. PCR product quality control was done by 384 well agarose gels, positive products where rearrayed and diluted to appropriate concentrations. For plasmid preparation, cultures were grown in 384 deepwell plates containing 200 μ l of 2YT Media. Incubation took 17–18 h at 1100 rpm.
Plasmids were isolated and purified automatically by a modified alkaline lysis procedure. Quality control of plasmids was done by randomly testing eight samples of each plate on agarose gels. Sequences were determined using M13 -40 primer, ABI Big Dye Terminator 3 chemistry and ABI 3730XL capillary sequencer systems (ABI, Weiterstadt, Germany). All raw sequences were processed by PHRED [50, 51] and screened for vector or Escherichia coli contamination.
Before clustering, sheep EST sequences were subjected to an extensive quality control procedure. In a first step sequences were clipped at both ends as long as PHRED quality values stayed below 20 in a window of size 15 bp. Secondly, for all sequences (ESTs + mRNAs) putative vector and repeat sequences were masked. Clustering and assembly of the clipped sequences was performed using the GeneNest software tools [12, 13].
An Ovis aries EST database
GeneNest is designed to generate and visualize gene indices based on ESTs and mRNA sequences [12, 13]. We have used the GeneNest software to create a database with 47209 sheep ESTs generated in the current project, as well as all previously available Ovis aries sequences. The database is available at http://genenest.molgen.mpg.de/cgi-bin/search_db?db=Oa5.
All ESTs generated in this project were submitted to NCBI's GenBank and are available with the accession numbers [GenBank:DY475559–DY505707] and [GenBank:DY505714–DY522957].
Identification and annotation of genes
In order to identify and annotate novel and previously known sheep genes, similarities to known proteins were identified with a blastx  search against the RefSeq protein database  of 28 September, 2005 using default settings for blastx. Contigs with a best e-value of 10-25 or lower were assigned a putative identity.
For the purpose of Gene Ontology  annotation, the contigs were mapped to mouse mRNAs using blastn with default settings. Putative identities were assigned using an e-value threshold of 10-25. Mouse-Genome Informatics (MGI) annotations  were assigned to the sheep homologs of annotated mouse genes, and GO analysis was performed using the Ontologizer .
Differential expression of early fracture genes
In order to gain insight into the genes that are upregulated in the early stages of fracture healing in the sheep, we compared the distribution of ESTs in the PF7 and PF10 sheep libraries against the distribution of ESTs in the UniGene  human bone libraries ID.655  and ID. 1124, which have a total of 6573 ESTs. Mapping of sheep contigs to human UniGene clusters was performed using blastn with an e-value cutoff of 10-25 as described above. We then compared the EST counts for the PF7/PF10 libraries against those in the ID.655/ID.1124 libraries using a two-tailed Fisher-Exact test. The result of the Fisher-exact test is a p-value that reflects the probability of observing the given differences in EST counts in the respective libraries assuming no differential expression . The expression ratio was calculated based on the total counts of 32032 ESTs in libraries PF7 and PF10 and 6573 ESTs in UniGene libraries ID.655/ID.1124. A similar test was applied to compare EST counts in PF7 against those in PF10. If the p-value of this association test is below a certain threshold, the gene is considered as putatively differentially expressed. Correction for multiple testing was performed using a false discovery rate (FDR)-control method . The FDR represents the expected proportion of errors among the rejected null hypotheses and is thus a different kind of procedure from multiple testing correction procedures such as a Bonferroni correction that control the family-wise error rate, or the probability that one or more rejected null hypotheses are falsely rejected. The Benjamini-Hochberg-adjusted p-values are shown in the Tables for the EST comparisons.
RNA was isolated for quantitative reverse-transcriptase polymerase chain reaction (RT-PCR) from 7-day postfracture (PF7, n = 4 animals), 10-day postfracture (PF10, n = 3), 14 day (PF14, n = 1) and 42-day (PF42, n = 1) tissue. PF7 and PF10 were identical to the RNA sources used for EST library construction, and PF14 was used for the pooled library. For semi-quantitative RT-PCR measurements, single-stranded cDNA was prepared from the different healing stages using random hexamer primers with the Taqman Reverse Transcription Reagents Kit (Applied Biosystems, Forster City, USA). Gene-specific primers were designed using the Primer Express software (Applied Biosystems) and quantitative RT-PCR was performed with SYBR green PCR Master Mix (Applied Biosystems) on an ABI Prism 7900 HT (Applied Biosystems). The mRNA levels were quantified in triplicate by the standard curve method normalized to a exogenous green fluorescent protein (GFP) transcript of which 1 ng was spiked into the sample per 1 μ g of total RNA. See [Additional file 2] for primer sequences and RT-PCR profiles.
Three measurements were made for each gene at each of the four time points. Statistical analysis of the RT-PCR results was performed by a two-sided Student's t test to compare expression between PF7 and PF10, and one-way ANOVA was performed to identify differential expression among the groups PF7, PF10, PF14 or among all four groups. A Bonferroni correction for multiple testing was performed.
Clustering and visualization of RT-PCR results
Complete-linkage hierarchical clustering was used to visualize the results of RT-PCR on the 78 putative fracture-related genes. For each gene, a vector of four values representing the four time points was normalized to have zero mean and unit variance. In order to estimate the optimal number of clusters, the gap statistic method  was used.
- ACP5 = Acid phosphatase 5:
Actin, α 1
Bone morphogenetic protein 1
Collagen (COL1A1, COL1A2, COL3A1, COL4A1, COL4A2, COL6A1, COL6A2)
Cytochrome c oxidase subunit I
Cysteine and Glycine-Rich Protein 1
Cathepsin (CTSD, CTSK, CTSS)
Eukaryotic translation elongation factor (EEF1A1, EEF2) EIF = Eukaryotic translation initiation factor (EIF4A1, EIF5A)
Fc fragment of IgG, receptor, transporter, α
Ferritin, heavy polypeptide 1
Ferritin, light polypeptide
G protein, α inhibiting activity polypeptide 2
Glutathione peroxidase (GPX1, GPX3, GXP4, GPX7)
Hepatoma-derived growth factor
Major histocompatibility complex, class I, A
Heme oxygenase (decycling) 1)
DNA-binding protein inhibitor (ID1, ID3)
Integrin, alpha 5
Integrin, beta 1
Integral membrane protein 2C
Jun B proto-oncogene
α Tubulin, ubiquitous
KDEL (Lys-Asp-Glu-Leu) endoplasmic reticulum protein retention receptor (KDELR2, KDELR3)
Ribosomal protein SA
LDL receptor-related protein (LRP1, LRP10)
Mortality factor 4 like 1
Matrix metalloproteinase (MMP2, MMP9, MMP13, MMP14, MMP19)
Nuclease-sensitive element-binding protein 1
phosphoglycerate kinase 1
Peptidylprolyl isomerase B (cyclophilin B)
Retinoic acid receptor, alpha
Ribosomal protein L (RPL3, RPL6, RPL7, RPL7A, RPL9, RPL10, RPLlOa, RPL17, RPL29, RPL41, RPLPO)
Ribosomal protein S (RPS2, RPS3, RPS4X, RPS5)
S100 calcium binding protein A4
Serpin peptidase inhibitor, clade H, member 1
Solute carrier family 25 member 6
Transcription factor 4
β Tropomyosin 2
Schell H, Epari DR, Kassi JP, Bragulla H, Bail HJ, Duda GN: The course of bone healing is influenced by the initial shear fixation stability. J Orthop Res. 2005, 23 (5): 1022-1028. 10.1016/j.orthres.2005.03.005.
Arav A, Revel A, Nathan Y, Bor A, Gacitua H, Yavin S, Gavish Z, Uri M, Elami A: Oocyte recovery, embryo development and ovarian function after cryopreservation and transplantation of whole sheep ovary. Hum Reprod. 2005, 20 (12): 3554-3559. 10.1093/humrep/dei278.
Ross MG, Desai M, Guerra C, Wang S: Programmed syndrome of hypernatremic hypertension in ovine twin lambs. Am J Obstet Gynecol. 2005, 192 (4): 1196-1204. 10.1016/j.ajog.2005.01.006.
Dukarm RC, Steinhorn RH, Russell JA, Lakshminrusimha S, Swartz D, Cummings JJ: Selective type 5 phosphodiesterase inhibition alters pulmonary hemodynamics and lung liquid production in near-term fetal lambs. J Appl Physiol. 2005, 99 (6): 2331-2336. 10.1152/japplphysiol.00120.2005.
Hu J, Mungall C, Law A, Papworth R, Nelson JP, Brown A, Simpson I, Leckie S, Burt DW, Hillyard AL, Archibald AL: The ARKdb: genome databases for farmed and other animals. Nucleic Acids Res. 2001, 29: 106-110. 10.1093/nar/29.1.106.
Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Church DM, DiCuccio M, Edgar R, Federhen S, Helmberg W, Kenton DL, Khovayko O, Lipman DJ, Madden TL, Maglott DR, Ostell J, Pontius JU, Pruitt KD, Schuler GD, Schriml LM, Sequeira E, Sherry ST, Sirotkin K, Starchenko G, Suzek TO, Tatusov R, Tatusova TA, Wagner L, Yaschenko E: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2005, D39-D45. 33 Database
Cochrane G, Aldebert P, Althorpe N, Andersson M, Baker W, Baldwin A, Bates K, Bhattacharyya S, Browne P, van den Broek A, Castro M, Duggan K, Eberhardt R, Faruque N, Gamble J, Kanz C, Kulikova T, Lee C, Leinonen R, Lin Q, Lombard V, Lopez R, McHale M, McWilliam H, Mukherjee G, Nardone F, Pastor MPG, Sobhany S, Stoehr P, Tzouvara K, Vaughan R, Wu D, Zhu W, Apweiler R: EMBL Nucleotide Sequence Database: developments in 2005. Nucleic Acids Res. 2006, D10-D15. 10.1093/nar/gkj130. 34 Database
Lindsay R: The burden of osteoporosis: cost. Am J Med. 1995, 98 (2A): 9S-11S. 10.1016/S0002-9343(05)80038-5.
Gerstenfeld LC, Einhorn TA: Developmental aspects of fracture healing and the use of pharmacological agents to alter healing. J Musculoskelet Neuronal Interact. 2003, 3 (4): 297-303. discussion 320-1
Vortkamp A, Pathi S, Peretti GM, Caruso EM, Zaleske DJ, Tabin CJ: Recapitulation of signals regulating embryonic bone formation during postnatal growth and in fracture repair. Mech Dev. 1998, 71 (1–2): 65-76. 10.1016/S0925-4773(97)00203-7.
Hadjiargyrou M, Lombardo F, Zhao S, Ahrens W, Joo J, Ahn H, Jurman M, White DW, Rubin CT: Transcriptional profiling of bone regeneration. Insight into the molecular complexity of wound repair. J Biol Chem. 2002, 277 (33): 30177-30182. 10.1074/jbc.M203171200.
Haas SA, Beissbarth T, Rivals E, Krause A, Vingron M: GeneNest: automated generation and visualization of gene indices. Trends Genet. 2000, 16 (11): 521-523. 10.1016/S0168-9525(00)02116-8.
Krause A, Haas SA, Coward E, Vingron M: SYSTERS, GeneNest, SpliceNest: exploring sequence space from genome to protein. Nucleic Acids Res. 2002, 30: 299-300. 10.1093/nar/30.1.299.
Meyer RA, Meyer MH, Tenholder M, Wondracek S, Wasserman R, Garges P: Gene expression in older rats with delayed union of femoral fractures. J Bone Joint Surg Am. 2003, 85-A (7): 1243-1254.
Tibshirani R, Walther G, Hastie T: Estimating the number of clusters in a data set via the gap statistic. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 2001, 63 (2): 411-423. 10.1111/1467-9868.00293.
Keane OM, Zadissa A, Wilson T, Hyndman DL, Greer GJ, Baird DB, McCulloch AF, Crawford AM, McEwan JC: Gene expression profiling of Naive Sheep genetically resistant and susceptible to Gastrointestinal Nematodes. BMC Genomics. 2006, 7: 42-10.1186/1471-2164-7-42.
Li X, Quigg RJ, Zhou J, Ryaby JT, Wang H: Early signals for fracture healing. J Cell Biochem. 2005, 95: 189-205. 10.1002/jcb.20373.
Rundle CH, Wang H, Yu H, Chadwick RB, Davis El, Wergedal JE, Lau KHW, Mohan S, Ryaby JT, Baylink DJ: Microarray analysis of gene expression during the inflammation and endochondral bone formation stages of rat femur fracture repair. Bone. 2006, 38: 521-9. 10.1016/j.bone.2005.09.015.
Kenmochi N, Kawaguchi T, Rozen S, Davis E, Goodman N, Hudson TJ, Tanaka T, Page DC: A map of 75 human ribosomal protein genes. Genome Res. 1998, 8 (5): 509-523.
Bortoluzzi S, d'Alessi F, Romualdi C, Danieli GA: Differential expression of genes coding for ribosomal proteins in different human tissues. Bioinformatics. 2001, 17 (12): 1152-1157. 10.1093/bioinformatics/17.12.1152.
Koide T, Aso A, Yorihuzi T, Nagata K: Conformational requirements of collagenous pep tides for recognition by the chaperone protein HSP47. J Biol Chem. 2000, 275 (36): 27957-27963.
Shroff B, Smith T, Norris K, Pileggi R, Sauk JJ: Hsp 47 is localized to regions of type I collagen production in developing murine femurs and molars. Connect Tissue Res. 1993, 29 (4): 273-286.
Tiffee JC, Griffin JP, Cooper LF: Immunolocalization of stress proteins and extracellular matrix proteins in the rat tibia. Tissue Cell. 2000, 32 (2): 141-147. 10.1054/tice.2000.0097.
Foidart JM, Reddi AH: Immunofluorescent localization of type IV collagen and laminin during endochondral bone differentiation and regulation by pituitary growth hormone. Dev Biol. 1980, 75: 130-136. 10.1016/0012-1606(80)90149-9.
Carano RAD, Filvaroff EH: Angiogenesis and bone repair. Drug Discov Today. 2003, 8 (21): 980-989. 10.1016/S1359-6446(03)02866-6.
Pacicca DM, Patel N, Lee C, Salisbury K, Lehmann W, Carvalho R, Gerstenfeld LC, Einhorn TA: Expression of angiogenic factors during distraction osteogenesis. Bone. 2003, 33 (6): 889-898. 10.1016/j.bone.2003.06.002.
Paralkar VM, Nandedkar AK, Pointer RH, Kleinman HK, Reddi AH: Interaction of osteogenin, a heparin binding bone morphogenetic protein, with type IV collagen. J Biol Chem. 1990, 265 (28): 17281-17284.
Anstrom JA, Thore CR, Moody DM, Challa VR, Block SM, Brown WR: Morphometric assessment of collagen accumulation in germinal matrix vessels of premature human neonates. Neuropathol Appl Neurobiol. 2005, 31 (2): 181-190. 10.1111/j.1365-2990.2004.00626.x.
Chang DF, Belaguli NS, Iyer D, Roberts WB, Wu SP, Dong XR, Marx JG, Moore MS, Beckerle MC, Majesky MW, Schwartz RJ: Cysteine-rich LIM-only proteins CRP1 and CRP2 are potent smooth muscle differentiation cofactors. Dev Cell. 2003, 4: 107-118. 10.1016/S1534-5807(02)00396-9.
Jonker L, Arthur HM: Endoglin expression in early development is associated with vasculogenesis and angiogenesis. Mech Dev. 2002, 110 (1–2): 193-196. 10.1016/S0925-4773(01)00562-7.
Lebrin F, Goumans MJ, Jonker L, Carvalho RLC, Valdimarsdottir G, Thorikay M, Mummery C, Arthur HM, ten Dijke P: Endoglin promotes endothelial cell proliferation and TGF-beta/ALKl signal transduction. EMBO J. 2004, 23 (20): 4018-4028. 10.1038/sj.emboj.7600386.
Everett AD, Natron JV, Stoops T, Nakamura H, Tucker A: Hepatoma-derived growth factor is a pulmonary endothelial cell-expressed angiogenic factor. Am J Physiol Lung Cell Mol Physiol. 2004, 286 (6): L1194-L1201. 10.1152/ajplung.00427.2003.
Uusitalo H, Hiltunen A, Soderström M, Aro HT, Vuorio E: Expression of cathepsins B, H, K, L, and S and matrix metalloproteinases 9 and 13 during chondrocyte hypertrophy and endochondral ossification in mouse fracture callus. Calcif Tissue Int. 2000, 67 (5): 382-390. 10.1007/s002230001152.
Dang CM, Beanes SR, Lee H, Zhang X, Soo C, Ting K: Scarless fetal wounds are associated with an increased matrix metalloproteinase-to-tissue-derived inhibitor of metalloproteinase ratio. Plast Reconstr Surg. 2003, 111 (7): 2273-2285. 10.1097/01.PRS.0000060102.57809.DA.
Filanti C, Dickson GR, Martino DD, Ulivi V, Sanguineti C, Romano P, Palermo C, Manduca P: The expression of metalloproteinase-2, -9, and -14 and of tissue inhibitors-1 and -2 is developmentally modulated during osteogenesis in vitro, the mature osteoblastic phenotype expressing metalloproteinase-14. J Bone Miner Res. 2000, 15 (11): 2154-2168. 10.1359/jbmr.2000.15.11.2154.
Hieta N, Impola U, Lopez-Otfn C, Saarialho-Kere U, Kahari VM: Matrix metalloproteinase-19 expression in dermal wounds and by fibroblasts in culture. J Invest Dermatol. 2003, 121 (5): 997-1004. 10.1046/j.1523-1747.2003.12533.x.
Brage M, Abrahamson M, Lindström V, Grubb A, Lerner UH: Different cysteine proteinases involved in bone resorption and osteoclast formation. Calcif Tissue Int. 2005, 76 (6): 439-447. 10.1007/s00223-004-0043-y.
Schell H, Lienau J, Epari DR, Seebeck P, Exner C, Muchow S, Bragulla H, Haas NP, Duda GN: Osteoclastic activity begins early and increases over the course of bone healing. Bone. 2006, 38: 547-54. 10.1016/j.bone.2005.09.018.
Fragonas E, Pollesello P, Mlinárik V, Toffanin R, Grando C, Godeas C, Vittur F: Sensitivity of chondrocytes of growing cartilage to reactive oxygen species. Biochim Biophys Acta. 1998, 1425: 103-111.
Komatsu DE, Hadjiargyrou M: Activation of the transcription factor HIF-1 and its target genes, VEGF, HO-1, iNOS, during fracture repair. Bone. 2004, 34 (4): 680-688. 10.1016/j.bone.2003.12.024.
Watanabe S, Akagi R, Mori M, Tsuchiya T, Sassa S: Marked developmental changes in heme oxygenase-1 (HO-1) expression in the mouse placenta: correlation between HO-1 expression and placental development. Placenta. 2004, 25 (5): 387-395. 10.1016/j.placenta.2003.10.012.
Deleersnijder W, Hong G, Cortvrindt R, Poirier C, Tylzanowski P, Pittois K, Marck EV, Merregaert J: Isolation of markers for chondro-osteogenic differentiation using cDNA library subtraction. Molecular cloning and characterization of a gene belonging to a novel multigene family of integral membrane proteins. J Biol Chem. 1996, 271 (32): 19475-19482. 10.1074/jbc.271.32.19475.
Mundlos S, Otto F, Mundlos C, Mulliken JB, Aylsworth AS, Albright S, Lindhout D, Cole WG, Henn W, Knoll JH, Owen MJ, Mertelsmann R, Zabel BU, Olsen BR: Mutations involving the transcription factor CBFA1 cause cleidocranial dysplasia. Cell. 1997, 89 (5): 773-779. 10.1016/S0092-8674(00)80260-3.
Nakashima K, Zhou X, Kunkel G, Zhang Z, Deng JM, Behringer RR, de Crombrugghe B: The novel zinc finger-containing transcription factor osterix is required for osteoblast differentiation and bone formation. Cell. 2002, 108: 17-29. 10.1016/S0092-8674(01)00622-5.
Kenner L, Hoebertz A, Beil T, Keon N, Karreth F, Eferl R, Scheuch H, Szremska A, Amling M, Schorpp-Kistner M, Angel P, Wagner EF: Mice lacking JunB are osteopenic due to cell-autonomous osteoblast and osteoclast defects. J Cell Biol. 2004, 164 (4): 613-623. 10.1083/jcb.200308155.
Myokai F, Oyama M, Nishimura F, Ohira T, Yamamoto T, Arai H, Takashiba S, Murayama Y: Unique genes induced by mechanical stress in periodontal ligament cells. J Periodontal Res. 2003, 38 (3): 255-261.
Cho EA, Dressier GR: TCF-4 binds beta-catenin and is expressed in distinct regions of the embryonic brain and limbs. Mech Dev. 1998, 77: 9-18. 10.1016/S0925-4773(98)00131-2.
Klein P, Schell H, Streitparth F, Heller M, Kassi JP, Kandziora F, Bragulla H, Haas NP, Duda GN: The initial phase of fracture healing is specifically sensitive to mechanical conditions. J Orthop Res. 2003, 21 (4): 662-669. 10.1016/S0736-0266(02)00259-0.
Lienau J, Schell H, Duda GN, Seebeck P, Muchow S, Bail HJ: Initial vascularization and tissue differentiation are influenced by fixation stability. J Orthop Res. 2005, 23 (3): 639-645. 10.1016/j.orthres.2004.09.006.
Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8 (3): 175-185.
Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.
Pruitt KD, Tatusova T, Maglott DR: NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2005, D501-D504. 33 Database
Ashburner M, Ball C, Blake J, Botstein D, Butler H, Cherry J, Davis A, Dolinski K, Dwight S, Eppig J, Harris M, Hill D, Issel-Tarver L, Kasarskis A, Lewis S, Matese J, Richardson J, Ringwald M, Rubin G, Sherlock G: Gene Ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature Genetics. 2000, 25: 25-29. 10.1038/75556.
Blake JA, Richardson JE, Bult CJ, Kadin JA, Eppig JT: MGD: the Mouse Genome Database. Nucleic Acids Research. 2003, 31: 193-195. 10.1093/nar/gkg047.
Robinson PN, Wollstein A, Böhme U, Beattie B: Ontologizing gene-expression microarray data: characterizing clusters with Gene Ontology. Bioinformatics. 2004, 20 (6): 979-981. 10.1093/bioinformatics/bth040.
Adams MD, Kerlavage AR, Fleischmann RD, Fuldner RA, Bult CJ, Lee NH, Kirkness EF, Weinstock KG, Gocayne JD, White O: Initial assessment of human gene diversity and expression patterns based upon 83 million nucleotides of cDNA sequence. Nature. 1995, 377 (6547 Suppl): 3-174.
Schmitt AO, Specht T, Beckmann G, Dahl E, Pilarsky CP, Hinzmann B, Rosenthal A: Exhaustive mining of EST libraries for genes differentially expressed in normal and tumour tissues. Nucleic Acids Res. 1999, 27 (21): 4251-4260. 10.1093/nar/27.21.4251.
Benjamini Y, Hochberg Y: Controlling the false discovery rate – A practical and powerful approach to multiple testing. J Roy Stat Soc B. 1995, 57: 289-300.
Wang Y, Middleton F, Horton JA, Reichel L, Farnum CE, Damron TA: Microarray analysis of proliferative and hypertrophic growth plate zones identifies differentiation markers and signal pathways. Bone. 2004, 35 (6): 1273-1293. 10.1016/j.bone.2004.09.009.
Paraoanu LE, Weiss B, Robitzki AA, Layer PG: Cytochrome-c oxidaseis one of several genes elevated in marginal retina of the chick embryo. Neuroscience. 2005, 132 (3): 665-672. 10.1016/j.neuroscience.2004.11.055.
Bland YS, Critchlow MA, Ashhurst DE: The expression of the fibrillar collagen genes during fracture healing: heterogeneity of the matrices and differentiation of the osteoprogenitor cells. Histochem J. 1999, 31 (12): 797-809. 10.1023/A:1003954104290.
Sato M, Yasui N, Nakase T, Kawahata H, Sugimoto M, Hirota S, Kitamura Y, Nomura S, Ochi T: Expression of bone matrix proteins mRNA during distraction osteogenesis. J Bone Miner Res. 1998, 13 (8): 1221-1231. 10.1359/jbmr.1918.104.22.1681.
Ou JH, Yen TS, Wang YF, Kam WK, Rutter WJ: Cloning and characterization of a human ribosomal protein gene with enhanced expression in fetal and neoplastic cells. Nucleic Acids Res. 1987, 15 (21): 8919-8934.
Liu X, Wu H, Byrne M, Krane S, Jaenisch R: Type III collagen is crucial for collagen I fibrillogenesis and for normal cardiovascular development. Proc Natl Acad Sci USA. 1997, 94 (5): 1852-1856. 10.1073/pnas.94.5.1852.
Brown TL, Fischer WC, Collins MD, De BK, Scott WJ: Identification of a 100-kDa phosphoprotein in developing murine embryos as elongation factor 2. Arch Biochem Biophys. 1994, 309: 105-110. 10.1006/abbi.1994.1091.
Gersdorff N, Müller M, Otto S, Poschadel R, Hübner S, Miosge N: Basement membrane composition in the early mouse embryo day 7. Dev Dyn. 2005, 233 (3): 1140-1148. 10.1002/dvdy.20425.
Armes N, Fried M: The genomic organization of the region containing the Drosophila melanogaster rpL7a (Surf-3) gene differs from those of the mammalian and avian Surfeit loci. Mol Cell Biol. 1995, 15 (5): 2367-2373.
Portman MA, Chen SH, Xiao Y, Ning XH: Maturational changes in gene expression for adenine nucleotide translocator isoforms and betaFl-ATPase in rabbit heart. Mol Genet Metab. 1999, 66 (2): 75-79. 10.1006/mgme.1998.2793.
Hwang KC, Cui XS, Park SP, Shin MR, Park SY, Kim EY, Kim NH: Identification of differentially regulated genes in bovine blastocysts using an annealing control primer system. Mol Reprod Dev. 2004, 69: 43-51. 10.1002/mrd.20156.
Hayman AR, Jones SJ, Boyde A, Foster D, Colledge WH, Carlton MB, Evans MJ, Cox TM: Mice lacking tartrate-resistant acid phosphatase (Acp 5) have disrupted endochondral ossification and mild osteopetrosis. Development. 1996, 122 (10): 3151-3162.
Rudolph U, Finegold MJ, Rich SS, Harriman GR, Srinivasan Y, Brabet P, Boulay G, Bradley A, Birnbaumer L: Ulcerative colitis and adenocarcinoma of the colon in G alpha i2-deficient mice. Nat Genet. 1995, 10 (2): 143-150.
Apte SS, Fukai N, Beier DR, Olsen BR: The matrix metalloproteinase-14 (MMP-14) gene is structurally distinct from other MMP genes and is co-expressed with the TIMP-2 gene during mouse embryogenesis. J Biol Chem. 1997, 272 (41): 25511-25517. 10.1074/jbc.272.41.25511.
Gabison EE, Hoang-Xuan T, Mauviel A, Menashi S: EMMPRIN/CD147, an MMP modulator in cancer, development and tissue repair. Biochimie. 2005, 87 (3–4): 361-368. 10.1016/j.biochi.2004.09.023.
Ishii M, Koike C, Igarashi A, Yamanaka K, Pan H, Higashi Y, Kawaguchi H, Sugiyama M, Kamata N, Iwata T, Matsubara T, Nakamura K, Kurihara H, Tsuji K, Kato Y: Molecular markers distinguish bone marrow mesenchymal stem cells from flbroblasts. Biochem Biophys Res Commun. 2005, 332: 297-303. 10.1016/j.bbrc.2005.04.118.
Zhao Y, Guan H, Liu SF, Wu RC, Wang Z: Overexpression of QM induces cell differentiation and mineralization in MC3T3-E1. BiolPharm Bull. 2005, 28 (8): 1371-1376. 10.1248/bpb.28.1371.
Lu ZH, Books JT, Ley TJ: YB-1 is important for late-stage embryonic development, optimal cellular stress responses, and the prevention of premature senescence. Mol Cell Biol. 2005, 25 (11): 4625-4637. 10.1128/MCB.25.11.4625-4637.2005.
Globus RK, Amblard D, Nishimura Y, Iwaniec UT, Kim JB, Almeida EAC, Damsky CD, Wronski TJ, van der Meulen MCH: Skeletal phenotype of growing transgenic mice that express a function-perturbing form of betal integrin in osteoblasts. Calcif Tissue Int. 2005, 76: 39-49. 10.1007/s00223-004-0309-4.
Keene DR, Sakai LY, Burgeson RE: Human bone contains type III collagen, type VI collagen, and fibrillin: type III collagen is present on specific fibers that may mediate attachment of tendons, ligaments, and periosteum to calcified bone cortex. J Histochem Cytochem. 1991, 39: 59-69.
Shi GP, Sukhova GK, Kuzuya M, Ye Q, Du J, Zhang Y, Pan JH, Lu ML, Cheng XW, Iguchi A, Perrey S, Lee AME, Chapman HA, Libby P: Deficiency of the cysteine protease cathepsin S impairs microvessel growth. Circ Res. 2003, 92 (5): 493-500. 10.1161/01.RES.0000060485.20318.96.
Angelastro JM, Töröcsik B, Greene LA: Nerve growth factor selectively regulates expression of transcripts encoding ribosomal proteins. BMC Neurosci. 2002, 3: 3-10.1186/1471-2202-3-3.
Menaa C, Devlin RD, Reddy SV, Gazitt Y, Choi SJ, Roodman GD: Annexin II increases osteoclast formation by stimulating the proliferation of osteoclast precursors in human marrow cultures. J Clin Invest. 1999, 103 (11): 1605-1613.
Wagner EF: Functions of AP1 (Fos/Jun) in bone development. Ann Rheum Dis. 2002, 61 (Suppl 2): ii40-ii42.
Raouf A, Ganss B, McMahon C, Vary C, Roughley PJ, Seth A: Lumican is a major proteoglycan component of the bone matrix. Matrix Biol. 2002, 21 (4): 361-367. 10.1016/S0945-053X(02)00027-6.
Hatori K, Sasano Y, Takahashi I, Kamakura S, Kagayama M, Sasaki K: Osteoblasts and osteocytes express MMP2 and -8 and TIMP1, -2, and -3 along with extracellular matrix molecules during appositional bone formation. Anat Rec A Discov Mol Cell Evol Biol. 2004, 277 (2): 262-271. 10.1002/ar.a.20007.
Gustafsson E, Aszodi A, Ortega N, Hunziker EB, Denker HW, Werb Z, Fassler R: Role of collagen type II and perlecan in skeletal development. Ann N Y Acad Sci. 2003, 995: 140-150.
Inada M, Wang Y, Byrne MH, Rahman MU, Miyaura C, López-Otín C, Krane SM: Critical roles for collagenase-3 (Mmpl3) in development of growth plate cartilage and in endochondral ossification. Proc Natl Acad Sci USA. 2004, 101 (49): 17192-17197. 10.1073/pnas.0407788101.
Weston AD, Hoffman LM, Underhill TM: Revisiting the role of retinoid signaling in skeletal development. Birth Defects Res C Embryo Today. 2003, 69 (2): 156-173. 10.1002/bdrc.10010.
Bennett JH, Carter DH, Alavi AL, Beresford JN, Walsh S: Patterns of integrin expression in a human mandibular explant model of osteoblast differentiation. Arch Oral Biol. 2001, 46 (3): 229-238. 10.1016/S0003-9969(00)00114-X.
Takahara K, Lyons GE, Greenspan DS: Bone morphogenetic protein-1 and a mammalian tolloid homologue (mTld) are encoded by alternatively spliced transcripts which are differentially expressed in some tissues. J Biol Chem. 1994, 269 (51): 32572-32578.
Duarte WR, Shibata T, Takenaga K, Takahashi E, Kubota K, Ohya K, Ishikawa I, Yamauchi M, Kasugai S: S100A4: a novel negative regulator of mineralization and osteoblast differentiation. J Bone Miner Res. 2003, 18 (3): 493-501. 10.1359/jbmr.2003.18.3.493.
Maeda Y, Tsuji K, Nifuji A, Noda M: Inhibitory helix-loop-helix transcription factors Idl/Id3 promote bone formation in vivo. J Cell Biochem. 2004, 93 (2): 337-344. 10.1002/jcb.20154.
Helseth DL, Veis A: Cathepsin D-mediated processing of procollagen: lysosomal enzyme involvement in secretory processing of procollagen. Proc Natl Acad Sci USA. 1984, 81 (11): 3302-3306. 10.1073/pnas.81.11.3302.
The authors gratefully acknowledge support by the Deutsche Forschungsgemeinschaft (KF102/2). RR is supported by the BMBF (NGFN-2, grant 01GR0414). The authors would like to thank Birol Köysüren and Janina Thiel for expert technical assistance.
JH performed and contributed to the design of the biological experiments with the help of ACS and VS. HK and AJP performed the EST sequencing under the guidance of RR. SAH analyzed and clustered the ESTs using GeneNest. SB and PNR were responsible for the bioinformatic analysis. JL and HS performed experiments on the sheep model under the guidance of GND. SM conceived and designed the study. PNR wrote the manuscript with the help of SB and JH. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: A gzip-compressed FASTA file containing sequences of the 19087 contigs described in this study. (GZ 4 MB)
Additional file 2: A complete list of sheep genes upregulated in early fracture healing as determined by EST analysis with links to human homologs in UniGene and to the GeneNest database. Additionally, RT-PCR primer sequences and detailed results are shown. (HTML 74 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Hecht, J., Kuhl, H., Haas, S.A. et al. Gene identification and analysis of transcripts differentially regulated in fracture healing by EST sequencing in the domestic sheep. BMC Genomics 7, 172 (2006). https://doi.org/10.1186/1471-2164-7-172
- Fracture Healing
- Bone Development
- Intramembranous Ossification
- Enchondral Ossification
- Putative Identity