- Research article
- Open Access
Canine fibroblast growth factor receptor 3 sequence is conserved across dogs of divergent skeletal size
BMC Genetics volume 9, Article number: 67 (2008)
Fibroblast growth factor receptor 3 (FGFR3) is expressed in the growth plate of endochondral bones and serves as a negative regulator of linear bone elongation. Activating mutations severely limit bone growth, resulting in dwarfism, while inactivating mutations significantly enhance bone elongation and overall skeletal size. Domesticated dogs exhibit the greatest skeletal size diversity of any species and, given the regulatory role of FGFR3 on growth plate proliferation, we asked whether sequence differences in FGFR3 could account for some of the size differences.
All exons, the promoter region, and 60 bp of the 3' flanking region of the canine FGFR3 gene were sequenced for nine different dog breeds representing a spectrum of skeletal size. The resultant sequences were compared to the reference Boxer genome sequence.
There was no variation in sequence for any FGFR3 exons, promoter region, or 3' flanking sequence across all breeds evaluated.
The results suggest that, regardless of domestication selection pressure to develop breeds having extreme differences in skeletal size, the FGFR3 gene is conserved. This implies a critical role for this gene in normal skeletal integrity and indicates that other genes account for size variability in dogs.
Fibroblast growth factor receptor 3 (FGFR3) is a membrane-bound tyrosine kinase receptor that regulates cellular proliferation within the growth plate of long bones. In bone elongation, FGFR3 serves to limit the proliferative activity of the chondrocytes while promoting differentiation and contributing to the mineralization at the chondro-osseus junction . Mice and humans with a gain of function mutation have impaired bone elongation, resulting in achondroplasia [2–4], whereas loss of function mutations result in skeletal overgrowth and severe appendicular abnormalities [5–7].
In mice and humans, the FGFR3 gene contains 19 exons  that are differentially spliced to form two distinct isoforms that differ in tissue expression, ligand binding affinity, and cellular response . The two FGFR3 isoforms, IIIb and IIIc, are generated by alternative RNA splicing of exon 7 to either exon 8 or 9. A third isoform, IIIa has been described for other members of the FGFR family but not for FGFR3 . The alternative splicing creates receptors with distinct extracellular binding domains with different ligand binding specificities and differential expression: FGFR3 IIIb is expressed in epithelial cells and FGFR3 IIIc is expressed in mesenchymal-derived cells [6, 11]. Mice lacking FGFR3 IIIc display significant skeletal overgrowth, exaggerated limb growth, distorted growth plates indicative of elevated proliferation, and diminished mineralization . Mice lacking the FGFR3 IIIb isoform do not exhibit those skeletal phenotypes  indicating the IIIc form as critical for normal skeletal development.
Mutations in the ligand binding domain, the transmembrane domain, or the tyrosine kinase domains have all been associated with constitutive activation and short stature . In contrast, genetically engineered mice that fail to express functional FGFR3 exhibit extreme skeletal overgrowth [12, 13]. In sheep, a naturally-occurring mutation in FGFR3 causes inactivation of a kinase domain and results in similar excessive growth . Furthermore, sheep heterozygous for the naturally-occurring loss of FGFR3 function mutation exhibit enhanced frame size without the detrimental skeletal effects associated with two mutant alleles . The reported FGFR3 mutations in humans, mice and sheep arise predominantly from point mutations in the coding regions [16, 17]. Depending on the location of a single nucleotide polymorphism (SNP) in FGFR3, "graded activation" of the gene occurs, generating a spectrum of skeletal size . For example, a SNP in the transmembrane domain of FGFR3 commonly results in achondrodysplasia, a specific skeletal disorder that results in short stature and disproportionately short limbs . However, point mutations in the second tyrosine kinase domain of FGFR3 can result in lethal thanatophoric dysplasia or severe achondroplasia with developmental delay and acanthosis nigricans [17, 18].
As a consequence of domestication and selective breeding, dogs exhibit the greatest range of skeletal size diversity in any single species. Given the pivotal role of FGFR3 in modulating skeletal size, we asked whether SNPs within the FGFR3 gene could account for height variations seen in the three different Poodle varieties specifically selected on wither height: Toy (less than or equal to 25.4 cm), Miniature (between 25.4 and 38.1 cm), and Standard (greater than 38.1 cm). In addition, the FGFR3 gene was sequenced in several dog breeds defined as chondrodysplastic by Stockard . These breeds exhibit dwarfism as a fixed trait with the entire population of a given breed sharing the same mutation and exhibiting the same altered limb morphology [20, 21]. Similarities between human and canine phenotypes suggest a potential role for FGFR3 in different dog breeds displaying dwarfism.
The Ensembl FGFR3 gene model was used to define putative exons for the canine FGFR3 gene  including the FGFR3 IIIb and FGFR3 IIIc isoforms. The primer pairs for the promoter, the 19 exons comprising FGFR3 IIIb and FGFR3 IIIc, and 60 bases at the 3' end of the FGFR3 gene (Table 1) were developed based on the homology of the canine Ensembl reference Boxer sequence and the published human genome sequence. In many cases the amplicon for a given exon included adjacent intronic sequences. A predicted exon defined in the Ensembl FGFR3 IIIc isoform model, located between exon 4 and 5, was excluded from sequencing because of complete lack of homology to other FGFR3 cDNA transcript sequences (human, mouse, and cow).
All amplicons were sequenced for three to eight dogs of each of the three Poodle varieties, and, representing chondrodysplastic breeds, a single Basset hound, Dachshund, Pekingese, Cardigan Welsh Corgi, and Pembroke Welsh Corgi were sequenced. A Nova Scotia Duck Tolling Retriever and Chesapeake Bay Retriever, both affected with short stature, curved limb skeletal dysplasia, were also sequenced as was a Dalmatian that exhibited no skeletal abnormalities. All dogs sequenced were registered purebreds with the American Kennel Club with the Poodles unrelated at the grandparent level. The Nova Scotia Duck Tolling Retriever and Chesapeake Bay Retriever were diagnosed as affected with chondrodysplasia by x-ray. Buccal swab-derived DNA was extracted for each Poodle as previously described  and then subjected to PCR amplification using primers specific for the regions of interest. The 20 μL PCR reactions contained 0.5 U Amplitaq Gold DNA Polymerase (Applied Biosystems, Foster City, CA, USA), 200 μm dNTPs (Applied Biosystems), 1.5 mm MgCl2 (Applied Biosystems), 1 × Geneamp PCR Gold buffer (Applied Biosystems); 0.2 μm of each forward and reverse oligonucleotide primer (Operon Biotechnologies, Huntsville, AL, USA), and 10–20 ng template genomic DNA. Reaction mixtures were subjected to a thermal cycling program of 94°C for 12 min, followed by 35 cycles of 94°C for 30 s, 45 s at the annealing temperature (see Table 1), and 72°C for 45 s and a final elongation stage of 72°C for 10 min were completed. Blood samples for the other breeds were collected from patients of the University of California, Davis, Veterinary Medical Teaching Hospital. DNA was extracted using the QiaAmp Blood Mini Kit (Qiagen) and PCR amplified as described. The amplified products for each region for each dog were gel purified with a gel-purification kit (Qiagen, Valencia, Calif.) and sequenced at an automated DNA sequencing facility (Davis Sequencing Inc, Davis, Calif.). Alternatively, some amplified products were purified with the QiaQuick PCR Purification Kit (Qiagen), sequencing reactions carried out using the Big Dye Terminator v3.1 cycle sequencing kit (Applied Biosystems), and products processed on an ABI 3100 Genetic Analyzer.
The derived sequences for each exon and the flanking promoter and 3'end were aligned with the reference Boxer sequence (Vector NTI Advance, version 10.3.0, Invitrogen Corp, Carlsbad Calif). Any potential SNPs in the sequences were confirmed or rejected by sequencing the complementary strand as well as sequencing additional dogs of the breed that exhibited the SNP. The 699 bp upstream of exon 1 was analyzed for the presence of promoter elements using the Genomatix  software tools Gene2Promoter and MatInspector .
Remarkably, Poodles of all three varieties that differed in skeletal size had fully conserved FGFR3 sequence across all 19 exons, the promoter region, and the downstream 3' flanking sequence. Looking across breeds that exhibit extreme skeletal morphologic differences also revealed 100% exon sequence conservation relative to the published reference Boxer sequence. As predicted by Ensembl, there are 17 introns for the canine FGFR3 gene. The primers used for sequencing the exons also included complete intron sequencing for ten of those introns and partial coverage for the remaining seven. Based upon the predicted sizes of the introns, partial sequence coverage ranged from a high of ~99% of the intronic region sequenced to a low of just over 15% for the largest introns. The only sequence variation observed relative to the reference Boxer sequence, across all dogs evaluated, was in intronic regions. Single SNPs were identified within introns 1, 8, 10, 12, and 16 for a total of 5 SNPs (Table 2). The identified intronic SNPs were sporadic and not present in all the dogs exhibiting a common phenotype.
Analysis of the 699 bp region upstream of exon 1 using the Gene2Promoter software indicated that no classical promoter elements were apparent although MatInspector identified several possible transcription factor elements within the promoter region sequenced. The most notable transcription factor elements included two CCAAT box promoter elements, a GA box, three cAMP response elements, three GC box Sp1-binding sites, two RXR heterodimer binding sites, and two estrogen response elements. Additionally, two Prox1 binding sites were identified in the intron between exons 1 and 2.
In addition to differential skeletal size observed in sheep carrying mutations in FGFR3, greater than 95% of human dwarfism achondrodysplasia cases are caused by a gain-of-function mutation in FGFR3 . Thus, the FGFR3 gene is the foremost candidate gene for regulating the diversity of skeletal size observed in dog breeds. The lack of sequence divergence in the FGFR3 gene or its regulatory regions across all the dogs and dog breeds assessed in this study suggests that selection pressure to maintain this precise sequence remains very strong. The conservation of this gene persists even in the face of human breeding schemes and breed development as seen by the absence of sequence mutations within the coding region of the gene across all the diverse dog breeds analyzed. Rare and only minor sequence divergence was observed for intronic regions. The intronic SNPs identified were deemed neutral as they were not associated with a common phenotype or height category. Due to the homozygosity present in breeds exhibiting fixed traits, it was unnecessary to sequence additional individuals from each breed.
Achondroplasia in humans is most frequently due to constitutive activation mutations within the transmembrane domain of FGFR3 . However, when Martinez et al.,  sequenced only the FGFR3 transmembrane domain, they found no difference in that sequence when comparing achondroplastic dogs to a non-achondroplastic dog, the German shepherd. Likewise, dwarfism in Dexter cattle is not due to mutations in the transmembrane domain of FGFR3 . These findings, coupled with the current results, imply the essentiality of the FGFR3 gene in normal skeletal development. Spontaneous inactivating or genetically engineered null mutations of the FGFR3 gene result in appendicular abnormalities, enhanced proliferation of growth plate chondrocytes, and reduced cortical bone thickness [6, 28]. Furthermore, lifespan is reduced in genetically FGFR3 null mice .
The FGFR3 protein and mRNA show a high degree of sequence homology across divergent mammalian species. For example, a pairwise alignment of the dog, human, and mouse FGFR3 protein sequences exhibits a 96% homology for consensus of amino acid pairs among all ungapped positions in the IIIb isoform. Similarly in the dog FGFR3 IIIc isoform, a pairwise alignment to the human and mouse protein sequence results in 95% and 94% homology, respectively. As expected, when the sequenced dog FGFR3 genomic exons were compared with published dog ESTs (mRNA) they were identical and when aligned with human and mouse ESTs exhibited homologies similar to human and mouse protein alignments as described above .
Expression of FGFR3 has been proposed to be the essential link to coordinate growth plate cell proliferation with the osteogenesis and mineralization necessary to structurally support the elongating skeleton . Disruption of this coordination could jeopardize the integrity of the skeleton in terrestrial species. The pathology in animals with gain of function mutations likewise exhibit dysregulation between cellular proliferation at the growth plate and ossification . Thus, the selection pressure to maintain a coordinated balance between bone elongation and ossification apparently outweighed any human selection pressure to develop dog breeds with distinct morphologies.
Recent reports have identified a QTL association between IGF-I and skeletal size in the dog , reiterating the association between IGF-I and growth in livestock [32–34]. The IGF-I gene has a more expansive set of target tissues when compared to the restrictive targets of FGFR3. The systemic effects of IGF-I, with its multitude of actions, would be subject to less stringent selection pressure than a gene that balances structural support with overall size. Therefore, genetic change in the IGF-I locus could be tolerated whereas sequence change in FGFR3 would not be permitted.
Although a classical promoter was not identified by the sequence analysis of the dog FGFR3 promoter region, many genes utilize diverse core promoters that vary from the standard TATA-dependent promoters (reviewed in ) including the murine FGFR3 . The in silico FGFR3 promoter analysis did underscore the similarities and differences between dogs, mice, and humans. A pairwise alignment of the dog FGFR3 promoter region with mouse, and human resulted in a poor identity (52% and 49%, respectively). Nevertheless, the key functional transcriptional enhancer elements characterized in the mouse and human FGFR3 gene, such as cAMP response elements, SP1, and Prox1 binding sites [36, 37], were identified in the canine FGFR3 promoter region. This suggests common regulation of the expression and function of the FGFR3 gene. Interestingly, potential estrogen binding domains were identified within the dog FGFR3 promoter region. Estrogen is known to play a major role in growth plate arrest and ossification . Although the mechanism used in growth plate closure is unknown, the presence of estrogen binding domains suggests a mechanism for estrogen regulation of chondrocyte activity through FGFR3 during skeletal maturation .
The present study revealed homogeneity across many different breeds of dogs, indicating the essentiality of the FGFR3 gene product in maintaining skeletal integrity. Deviations from the consensus sequence appear to be quickly subjected to purifying selection in purebred dog breeds such that they are not maintained.
L'Hôte CG, Knowles MA: Cell responses to FGFR3 signalling: growth, differentiation and apoptosis. Exp Cell Res. 2005, 304: 417-31. 10.1016/j.yexcr.2004.11.012.
Wang Y, Spatz MK, Kannan K, Hayk H, Avivi A, Gorivodsky M, Pines M, Yayon A, Lonai P, Givol D: A mouse model for achondroplasia produced by targeting fibroblast growth factor receptor 3. Proc Natl Acad Sci USA. 1999, 96: 4455-60. 10.1073/pnas.96.8.4455.
Li C, Chen L, Iwata T, Kitagawa M, Fu XY, Deng CX: A Lys644Glu substitution in fibroblast growth factor receptor 3 (FGFR3) causes dwarfism in mice by activation of STATs and ink4 cell cycle inhibitors. Hum Mol Genet. 1999, 8: 35-44. 10.1093/hmg/8.1.35.
Horton WA, Hall JG, Hecht JT: Achondroplasia. Lancet. 2007, 370: 162-72. 10.1016/S0140-6736(07)61090-3.
Brodie SG, Deng CX: Mouse models orthologous to FGFR3-related skeletal dysplasias. Pediatr Pathol Mol Med. 2003, 22: 87-103. 10.1080/15227950307697.
Eswarakumar VP, Schlessinger J: Skeletal overgrowth is mediated by deficiency in a specific isoform of fibroblast growth factor receptor 3. Proc Natl Acad Sci USA. 2007, 104: 3937-42. 10.1073/pnas.0700012104.
Toydemir RM, Brassington AE, Bayrak-Toydemir P, Krakowiak PA, Jorde LB, Whitby FG, Longo N, Viskochil DH, Carey JC, Bamshad MJ: A novel mutation in FGFR3 causes camptodactyly, tall stature, and hearing loss (CATSHL) syndrome. Am J Hum Genet. 2006, 79: 935-41. 10.1086/508433.
Wüchner C, Hilbert K, Zabel B, Winterpacht A: Human fibroblast growth factor receptor 3 gene (FGFR3): genomic sequence and primer set information for gene analysis. Hum Genet. 1997, 100: 215-9. 10.1007/s004390050493.
Ornitz DM, Xu J, Colvin JS, McEwen DG, MacArthur CA, Coulier F, Gao G, Goldfarb M: Receptor Specificity of the fibroblast growth factor receptor family. J Biol Chem. 1996, 271: 15292-15297. 10.1074/jbc.271.25.15292.
Perez-Castro AV, Wilson J, Altherr MR: Genomic organization of the mouse fibroblast growth factor receptor 3 (Fgfr3) gene. Genomics. 1995, 30 (2): 157-62. 10.1006/geno.1995.9890.
Zhang X, Ibrahimi OA, Olsen SK, Umemori H, Mohammadi M, Ornitz DM: Receptor specificity of the fibroblast growth factor family. The complete mammalian FGF family. J Biol Chem. 2006, 281: 15694-700. 10.1074/jbc.M601252200.
Deng CX, Wynshaw-Boris A, Zhou F, Kuo A, Leder P: Fibroblast growth factor receptor 3 is a negative regulator of bone growth. Cell. 1996, 84 (6): 911-21. 10.1016/S0092-8674(00)81069-7.
Colvin JS, Bohne BA, Harding GW, McEwen DG, Ornitz DM: Skeletal overgrowth and deafness in mice lacking fibroblast growth factor receptor 3. Nat Genet. 1996, 12 (4): 390-7. 10.1038/ng0496-390.
Beever JE, Smit MA, Meyers SN, Hadfield TS, Bottema C, Albretsen J, Cockett NE: A single-base change in the tyrosine kinase II domain of ovine FGFR3 causes hereditary chondrodysplasia in sheep. Anim Genet. 2006, 37: 66-71. 10.1111/j.1365-2052.2005.01398.x.
Smith LB, Dally MR, Sainz RD, Rodrigue KL, Oberbauer AM: Enhanced skeletal growth of sheep heterozygous for an inactivated fibroblast growth factor receptor 3. J Anim Sci. 2006, 84: 2942-9. 10.2527/jas.2006-255.
Kannan K, Givol D: FGF receptor mutations: dimerization syndromes, cell growth suppression, and animal models. IUBMB Life. 2000, 49: 197-205.
Naski MC, Wang Q, Xu J, Ornitz DM: Graded activation of fibroblast growth factor receptor 3 by mutations causing achondroplasia and thanatophoric dysplasia. Nat Genet. 1996, 13 (2): 233-7. 10.1038/ng0696-233.
Shiang R, Thompson LM, Zhu YZ, Church DM, Fielder TJ, Bocian M, Winokur ST, Wasmuth JJ: Mutations in the transmembrane domain of FGFR3 cause the most common genetic form of dwarfism, achondroplasia. Cell. 1994, 78 (2): 335-42. 10.1016/0092-8674(94)90302-6.
Stockard CR: In the genetic and endocrinic basis for differences in form and behavior. The Wistar Institute of Anatomy and Biology, Philadelphia. 1941, Section II: 45-140.
Young AE, Bannasch DL: SNPs in the promoter regions of the canine RMRP and SHOX genes are not associated with canine chondrodysplasia. Anim Biotech. 2008, 19: 1-5. 10.1080/10495390701638328.
Pollinger JP, Bustamante CD, Fledel-Alon A, Schmutz S, Gray MM, Wayne RK: Selective sweep mapping of genes with large phenotypic effects. Genome Res. 2005, 15 (12): 1809-19. 10.1101/gr.4374505.
Ensembl Canis Familiaris FGFR3 gene website. [http://www.ensembl.org/Canis_familiaris/geneview?gene=ENSCAFG00000014993]
Oberbauer AM, Grossman DI, Eggleston ML, Irion DN, Schaffer AL, Pedersen NC, Belanger JM: Alternatives to blood as a source of DNA for large-scale scanning studies of canine genome linkages. Vet Res Commun. 2003, 27: 27-38. 10.1023/A:1022006606796.
Genomatix website. [http://www.genomatix.de]
Cartharius K, Frech K, Grote K, Klocke B, Haltmeier M, Klingenhoff A, Frisch M, Bayerlein M, Werner T: MatInspector and beyond: promoter analysis based on transcription factor binding sites. Bioinformatics. 2005, 21: 2933-42. 10.1093/bioinformatics/bti473.
Martínez S, Valdés J, Alonso RA: Achondroplastic dog breeds have no mutations in the transmembrane domain of the FGFR-3 gene. Can J Vet Res. 2000, 64: 243-5.
Usha AP, Lester DH, Williams JL: Dwarfism in Dexter cattle is not caused by the mutations in FGFR3 responsible for achondroplasia in humans. Anim Genet. 1997, 28: 55-7. 10.1111/j.1365-2052.1997.00068.x.
Vanek J, Alstad A, Berg I, Misek A, Moore B, Limesand W: Spider syndrome in lambs: A clinical and postmortem analysis. Vet Med. 1986, 81: 663-668.
Valverde-Franco G, Liu H, Davidson D, Chai S, Valderrama-Carvajal H, Goltzman D, Ornitz DM, Henderson JE: Defective bone mineralization and osteopenia in young adult FGFR3 -/- mice. Hum Mol Genet. 2004, 13: 271-84. 10.1093/hmg/ddh034.
Basic Local Alignment Search Tool (BLAST) website. [http://blast.ncbi.nlm.nih.gov/Blast.cgi]
Sutter NB, Bustamante CD, Chase K, Gray MM, Zhao K, Zhu L, Padhukasahasram B, Karlins E, Davis S, Jones PG, Quignon P, Johnson GS, Parker HG, Fretwell N, Mosher DS, Lawler DF, Satyaraj E, Nordborg M, Lark KG, Wayne RK, Ostrander EA: A single IGF1 allele is a major determinant of small size in dogs. Science. 2007, 316: 112-115. 10.1126/science.1137045.
Curi R, de Oliveira H, Silveira A, Lopes C: Association between IGF-I, IGF-IR and GHRH gene polymorphisms and growth and carcass traits in beef cattle. Livestock Prod Sci. 2005, 94: 159-167. 10.1016/j.livprodsci.2004.10.009.
Casas-Carrillo E, Prill-Adams A, Price SG, Clutter AC, Kirkpatrick BW: Relationship of growth hormone and insulin-like growth factor-1 genotypes with growth and carcass traits in swine. Anim Genet. 1997, 28: 88-93. 10.1111/j.1365-2052.1997.00086.x.
Moody DE, Pomp D, Newman S, MacNeil MD: Characterization of DNA polymorphisms in three populations of Hereford cattle and their associations with growth and maternal EPD in line 1 herefords. J Anim Sci. 1996, 74: 1784-93.
Smale ST, Kadonaga JT: The RNA polymerase II core promoter. Annu Rev Biochem. 2003, 72: 449-79. 10.1146/annurev.biochem.72.121801.161520.
McEwen DG, Green RP, Naski MC, Towler DA, Ornitz DM: Fibroblast growth factor receptor 3 gene transcription is suppressed by cyclic adenosine 3',5'-monophosphate. Identification of a chondrocytic regulatory element. J Biol Chem. 1999, 274 (43): 30934-42. 10.1074/jbc.274.43.30934.
Shin JW, Min M, Larrieu-Lahargue F, Canron X, Kunstfeld R, Nguyen L, Henderson JE, Bikfalvi A, Detmar M, Hong YK: Prox1 promotes lineage-specific expression of fibroblast growth factor (FGF) receptor-3 in lymphatic endothelium: a role for FGF signaling in lymphangiogenesis. Mol Biol Cell. 2006, 17 (2): 576-84. 10.1091/mbc.E05-04-0368.
Sävendahl L: Hormonal regulation of growth plate cartilage. Horm Res. 2005, 64 (Suppl 2): 94-7. 10.1159/000087764.
Perry RJ, Farquharson C, Ahmed SF: The role of sex steroids in controlling pubertal growth. Clin Endocrinol (Oxf). 2008, 68 (1): 4-15.
This work was supported by California Agricultural Experiment Station Funds (CA-D*-ASC-5256-AH) with partial funding also provided by the Center for Companion Animal Health, University of California, Davis. LBS received financial assistance from Scholarships funded by the Ford Family Foundation and the endowment of G. Kirk Swingle.
LBS participated in carrying out the study, sequencing coordination, sequence alignment and in drafting the manuscript. DLB participated in conceiving the study and in the study design. AEY participated in carrying out the study, sequencing, and sequence alignment. DIG and JMB participated in sequencing and sequence alignment. AMO participated in conceiving the study, study design, coordination, and in drafting the manuscript.
About this article
Cite this article
Smith, L.B., Bannasch, D.L., Young, A.E. et al. Canine fibroblast growth factor receptor 3 sequence is conserved across dogs of divergent skeletal size. BMC Genet 9, 67 (2008). https://doi.org/10.1186/1471-2156-9-67
- Growth Plate
- Skeletal Size
- Thanatophoric Dysplasia
- Bone Elongation