* Institute of Anthropology and Human Genetics, Division of Molecular Genetics, University of Tübingen, Tübingen, Germany; Institut für Ur- und Frühgeschichte, Abteilung Ältere Urgeschichte und Quartärökologie, University of Tübingen, Tübingen, Germany;
Institute of Organic Chemistry, University of Tübingen, Tübingen, Germany;
Institute of Pathology, Division of Palaeopathology, Academic Teaching Hospital München-Bogenhausen, München, Germany; || Institute of Human Genetics, University of Münster, Münster, Germany; and ¶ Department of Zoology, Natural History Museums and Botanical Garden, University of Oslo, Oslo, Norway
Correspondence: E-mail: bachmann{at}nhm.uio.no.
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Key Words: ancient DNA extracts ancient nuclear DNA diagenesis extract-induced mutations manganese PCR errors
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Achondroplasia (ACH; MIM 100800) is the most common form of short-limb skeletal disorder and occurs with a frequency of 1:15,000. The disease has an autosomal dominant Mendelian inheritance with complete penetrance. The phenotype results from the disruption of the continuous process of enchondral ossification (Stanescu, Stanescu, and Maroteaux 1990). The critical chromosome region for ACH was localized to 4p16.3 by linkage analyses (Le Merrer et al. 1994; Velinov et al. 1994). In this chromosomal region, the gene for the fibroblast growth factor receptor 3 (FGFR3) described by Keegan et al. (1991) is associated with ACH, hypochondroplasia (HCH), and thanatophoric dysplasia (TD). It is noteworthy that FGFR3 was originally considered a candidate gene for the Huntington disease (Thompson et al. 1991). The vast majority of ACH cases are caused by a GA transition (97%) at cDNA position 1138 causing a G380R amino acid substitution; the remaining 3% are caused by a G
C transversion at the same position resulting in the same amino acid substitution of G380R (Shiang et al. 1994; Rousseau et al. 1994; Wilkin et al. 1998).
Here we report on the identification and evaluation of diagnostic mutations for achondroplasia in a mummified specimen from Semerchet, Egypt, as well as in a Merovingian skeleton from Kirchheim, Germany, both of which had short stature. A series of mock controls (i.e., modern human blood and prehistoric bone samples from human individuals not affected by ACH) and negative controls were included to monitor sequence heterogeneity in PCR products and to guarantee purity of the assays, respectively. Theoretically, the dominant mode of inheritance of ACH accompanied by its complete penetrance will allow an unequivocal identification, because any contamination of DNA extracts with DNA bearing the diagnostic ACH point mutation must stem from short-statured individuals with achondroplasia, which is very unlikely.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Amino Acid Racemization
The amino acid content of the samples was determined quantitatively by enantiomer labeling (Frank, Nicholson, and Bayer 1978). The degree of racemization of alanine, aspartic acid, leucine, phenylalanine, and serine was determined according to Gerhardt and Nicholson (1994). For this purpose, 1.0 mg of dry pulverized compacta was hydrolyzed in 6 N DCl in D2O. Ethyl ester/TFA derivatives of the released amino acids were subsequently analyzed by GC/MS (Agilent 6890/5973) with SIM-detection on a 20 m x 0.25 mm fused silica capillary (30% Lipodex E / 70% PS255, film thickness of 0.13 µ). The degree of racemization of serine and phenylalanine indicated that bone glue treatment can be excluded as a source for exogeneous contaminating DNA (Nicholson et al. 2002).
DNA Extraction
Extraction of DNA from bone samples was done as described in Pusch and Scholz (1997). Modern human control DNA was isolated from white blood cells (Miller, Dykes, and Polesky 1988).
PCR Cloning and Sequencing
A 164 bp stretch of the FGFR3 gene (position 11211284 of GenBank entry M58051) was amplified using the primers ACHF and ACHR of Shiang et al. (1994). Alternatively, the primer pair ACHvF 5'-GTGTATGCAGGCATCCTCAG-3' (position 11531172 in M58051) and ACHR was used to amplify a shorter 132 bp stretch of the FGFR3 gene. The high fidelity Pfu polymerase (Stratagene, La Jolla, Calif.) with proofreading property was used for amplifications of ancient DNA templates and the spiking experiments. Taq polymerase (Roche, Basil, Switzerland) was used for control amplifications of modern template DNA and for the screening of buffers for contamination.
Spiking experiments were done by adding 5 µl DNA extract from an Ursus spelaeus (Pleistocene) sample to 50 ng modern human template DNA.
Cloning of PCR products into plasmid vectors followed the temperature cycle-ligation protocol of Pusch, Schmitt, and Blin (1997). DNA sequencing was performed on an automated DNA sequencer using BigDye chemistry (Applied Biosystems, Foster City, Calif.). All sequences are provided as Supplementary Material online.
Uracil N-glycosylase (UNG) Treatment
The DNA extractions from the phenotypically healthy Neresheim 1 individual were treated with UNG from Escherichia coli as recently described (Hofreiter et al. 2001). To test for a UNG-mediated reduction of deaminated sites at cDNA position 1138 of the FGFR3 gene the obtained clones of Neresheim 1 were screened with SfcI. SfcI allows the differentiation between wild type alleles and mutated sequences that are altered by the predominant GA transition.
Software
The Lasergene/Seqscape packages were employed for sequence analyses and ClustalX (Thompson et al. 1997) was used for multiple sequence alignments.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Amplification of the FGFR3 Sequence
A 164 bp segment of the FGFR3 gene (11211284 of GenBank entry M58051) was amplified by using the primer pair ACHFACHR. Amplicons of the expected size were obtained for six mock controls and the Semerchet sample but the other five samples failed even in repeated rounds of PCR (table 3).
|
A reduction of the size of the target sequence to 132 bp (11531284 of GenBank entry M58051) by using the primer pair ACHvFACHR yielded amplicons also in the case of the Kirchheim sample, which was refractory to PCR amplification in the previous experiments. The diagnostic transition was detected in the Semerchet and Kirchheim samples, but again false-positives were found in the four phenotypically healthy mock controls Neresheim 1 and 2, Warburg 2, and al Buhais A (see Supplementary Material online).
When looking at the sequences one sees that replacing primer ACHF with ACHvF to amplify the 132 bp fragment increased the mutation rate for some template DNAs, e.g., for the Neresheim 1 sample from 1.0 x 102 to 1.6 x 102 per position (ratio 0.62), for the Neresheim 2 sample from 4.5 x 103 to 2.3 x 102 per position (ratio 0.19), and for the al Buhais A sample from 6.7 x 103 to 1.8 x 102 per position (ratio 0.37). However, the specific mutation rate at cDNA position 1138 is almost the same, i.e., 57% versus 60% of clones for Neresheim 1, 42.8% versus 37.5% of clones for Neresheim 2, and 33.3% versus 50% for al Buhais A. Therefore, we conclude that the length of the amplified fragment and possible primer effects on the mutation rate per position are not relevant for the frequency of false-positive GA transitions at cDNA position 1138.
Quantitation of Ancient DNA Molecules
Sequence alterations that occur in only a few sequences of a cloned PCR product may result from misincorporation of nucleotides at various stages of a PCR. Degradation and damage of ancient DNA such as lesions (Hansen et al. 2001) or miscoding bases (Höss et al. 1996) are likely to increase the frequency of such misincorporations. It is also known that PCR amplifications that start from less than 1,000 template molecules tend to yield inconsistent results (Handt et al. 1996).
As an example, we determined the number of initial template molecules in the Semerchet sample. The 164 bp segment of the FGFR3 gene (primer pair ACHFACHR) was amplified quantitatively from the Semerchet extract in the presence of defined numbers of competitor molecules that had 55 bp deleted (fig. 1). It turned out that 5 µl Semerchet extract amplified almost as good as 10 initial competitor molecules. Thus, the number of template molecules is 12 per µl Semerchet extract. This means that the PCR experiments on the Semerchet extracts certainly started from less than 20 template molecules.
|
Spiking Experiments
As suggested by Pusch and Bachmann (2004), the observed GA transition at FGFR3 cDNA position 1138 in the false-positive samples as well as some other sequence alterations may be introduced during PCR as a consequence of the presence of ancient DNA extracts in the reaction mixture. We tested this hypothesis by adding DNA extracts from Pleistocene Ursus spelaeus bones to 50 ng modern human DNA template from a phenotypically healthy individual and amplified the 132 bp FGFR3 fragment. The diagnostic G
A transition at FGFR3 cDNA position 1138 was detected in five out of 10 clones (50%) as well as other sequence alterations. In a replication experiment 22 clones were screened by SfcI cleavage for the diagnostic G
A transition and five clones (22.7%) were positive as determined by the correct length of the resulting DNA fragments after electrophoresis (table 3). The unspiked modern human template DNA yielded the wild type sequence without alterations. Furthermore, the control amplification using the Ursus spelaeus DNA extract yielded no 132 bp FGFR3 amplicon, i.e., contamination with human DNA was not detectable.
Amplification of the 132 bp FGFR3 Fragment Under the Presence of MnCl2
Pusch and Bachmann (2004) suggested that agents coextracted with DNA such as multivalent metal ions that accumulate postmortem during diagenesis may induce mutations under PCR. In particular, manganese in combination with magnesium is known to induce mutations (e.g., Svetlov and Cooper 1998; Kunichika, Hashimoto, and Imoto 2002). We amplified 50 ng of the same modern DNA template that was used for the spiking experiment under the presence of 0.25 mM MnCl2 and sequenced 10 clones. We obtained basically the same pattern of sequence alteration as in the spiking experiments. The diagnostic GA transition at FGFR3 cDNA position 1138 was detected in four out of 10 clones (40%) as well as other sequence alterations. In two clones position 1138 was altered by a G
T transversion (see Supplementary Material online).
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
There might be various reasons to explain the occurrence of the false-positive GA transition at cDNA position 1138 in the FGFR3 gene. Here, three lines of arguments will be discussed.
Post Mortem Damage of Template DNA
After the death of an organism, natural degeneration of bio-/macromolecules begins and the speed of the decay depends on the particular conditions of the environment. Moisture, temperature, and soil chemistry are important parameters in this context. Under certain conditions of preservation with low to moderate levels of hydrolysis and oxidation, various sized fragments of the DNA (in particular those with a low molecular weight) may survive in a bone. High salt concentrations at neutral pH and/or fast dehydration/mummification are expected to favor DNA survival (Pääbo and Wilson 1991; Lindahl 1993; Burger et al. 1999; Ovchinnikov et al. 2000). In addition to extensive enzymatic and physical DNA fragmentation, the most important route of decay for hydrated DNA is depurination (Lindahl 2000), leading, in part, to nicked double strands (i.e., partially single-stranded DNA; Pusch, Giddings, and Scholz 1998; Di Bernardo et al. 2002). It was suggested that sequence alterations may originate from such gaps or lesions in the DNA template (Hansen et al. 2001).
In our study, false-positive diagnostic FGFR3 1138 GA transitions occurred frequently. The application of UNG was recently proposed to reduce or exclude artificial type 2 (G
A/C
T) transitions (Hofreiter et al. 2001). UNG is expected to cleave specifically uracil bases in the ancient DNA to generate apyrimidinic sites (Pu and Struhl 1992). Such affected/altered DNA strands are no longer suitable template molecules for subsequent PCR experiments because they will fall apart during the first denaturation step. A treatment of the Neresheim 1 extract did not eliminate the occurrence of the false-positive diagnostic FGFR3 1138 G
A alterations. It should be noted that the first denaturation step of a PCR will also inactivate the UNG. UNG treatment can, therefore, not rule out the occurrence of sequence alterations via the deoxyuracil pathway that are introduced after the UNG treatment, thus during PCR.
Low Number of Template Molecules, PCR Jumping, and Contamination of Samples
Quantitation experiments gave an estimate of 12 template molecules per µl for the Semerchet sample. In contrast, 35 different sequences (out of 50 clones analyzed) were generated from the Semerchet extract (132 bp PCR product). This indicates that miscoding lesions in the ancient template DNA cannot account for the observed sequence heterogeneity. Most of the heterogeneity of PCR products must, therefore, be artificially induced de novo. One may assume that this effect is due to low template numbers, a general problem when subjecting ancient DNA extracts to PCR. However, PCR experiments on diluted modern DNA templates always yielded the expected authentic sequences. We therefore conclude that the observed sequence heterogeneity is independent of the number of template molecules. PCR jumping (Pääbo 1989) as a possible source for the observed sequence variation can be ruled out through the spiking experiments (i.e., cave bear extracts added to contemporary human DNA template). The obtained sequences are without any doubt human and there is no indication that damaged cave bear DNA was involved in the amplification process. The Ursus spelaeus sample itself did not yield any amplicon and was, therefore, not contaminated with human DNA. Further support comes from the amplification of the 132 bp stretch of the FGFR3 gene from the same contemporary human template DNA that was used for the spiking experiment under the presence of MnCl2 (see below). Contamination through previously cloned fragments may also be a possible explanation for the occurrence of the false-positive GA transition. However, this could be ruled out because all amplicons from modern blood controls of healthy individuals (no diagnostic G
A transition in the FGFR3 gene at cDNA position 1138 observed) and some of the mock controls never yielded the diagnostic ACH mutation. Moreover, the water and blank extraction controls never yielded PCR products.
PCR-Induced Sequence Alterations
Pusch and Bachmann (2004) showed that ancient DNA extracts can induce mutations in a nonrandom fashion. They have amplified a 148 bp stretch of the mitochondrial HVRI from contemporary human template DNA in PCR reactions spiked with ancient DNA extracts. In total, 34 positions of a 103 bp alignment were affected and most mutations occurred repeatedly in independent PCR amplifications. The spiking experiments presented here (i.e., cave bear extracts added to contemporary human DNA template) support the hypothesis that the phenomenon of extract-induced mutations is responsible for the unreliable typing of achondroplasia. The obtained sequencesalthough derived from a healthy individualdisplayed frequently (31.3%) the diagnostic GA transition at cDNA position 1138 as well as some other mutations. This indicates that the observed sequence alterations were introduced during PCR amplification. In other words, it was possible to generate human ACH sequences from modern human nucleic acids due to the presence of ancient cave bear extracts in the PCR. We therefore have to conclude that ancient DNA extracts contribute currently unidentified components to the PCR mixture which either alter template molecules directly or reduce the fidelity of the Pfu polymerase. Pusch and Bachmann (2004) have suggested that multivalent metal ions such as manganese that accumulate postmortem during diagenesis may induce mutations under PCR. PCR amplifications of contemporary human template DNA under the presence of MnCl2 indicate that cDNA position 1138 of the FGFR3 gene is indeed particularly sensitive to mutagenesis. In 40% of the obtained sequences the diagnostic G
A transition was observed at this particular position (in addition, 20% of the obtained sequences had a G
T transversion). It is still not understood why particular nucleotides are more prone to erroneous misincorporation of nucleotides than others. The frequent occurrence of the diagnostic ACH mutation in false-positive ancient samples may be the result of processes affecting CG dinucleotides. For contemporary nucleic acids it is known that such CG dinucleotides are preferential targets for spontaneous point mutations. They account for one-third of single-site mutations observed in inherited human diseases and are among the major types of miscoding lesions in the genome of living human cells (Cooper and Youssoufian 1988; Rideout et al. 1990).
To summarize, the experiments presented here show that the GA transitions at position 1138 of the FGFR3 gene, which is diagnostic for achondroplasia, occurred artificially in amplicons from ancient DNA extracts of phenotypically healthy individuals and render a reliable molecular typing of ancient samples for the disease impossible. We conclude that this is due to extract-induced mutations under PCR as described by Pusch and Bachmann (2004). They present the analyses of 547 sequences from cloned amplicons of a 148 bp stretch of the mitochondrial HVRI from contemporary human template DNA generated in spiked PCR reactions. The authors observed that extract-induced mutations occurred in a nonrandom fashion in independent PCR amplifications and in total 34 positions of a 103 bp alignment were affected. It is noteworthy that 15.7% of all sequences analyzed by Pusch and Bachmann (2004) differed from the closest human match in GenBank by nine substitutions and three gaps but shared seven out of 11 mutations that are in combination characteristic for the Neandertal sequence AF011222 (Krings et al. 1997). Their data might challenge the authenticity of several published sequences. However, showing that a combination of nucleotide substitutions can be generated artificially does not necessarily imply that similar published sequences are not authentic. In the example of Neandertal sequences, it might be difficult to tell apart authentic nucleotide substitutions from extract-induced mutations, because no contemporary Neandertal control DNA exists. Theoretically, ACH offers a more straightforward test system. ACH is a skeletal disorder with an autosomal dominant Mendelian inheritance and complete penetrance and the resulting phenotype cannot be overlooked. This means that the diagnostic G
A transition at position 1138 of the FGFR3 gene (or a G
C transversion at the same position) is expected in short-statured individuals but must not occur in phenotypically healthy individuals. The observed false-positives are, therefore, unambiguously artifacts. However, since this particular mutation can occur as a false-positive there is no reason to conclude that the occurrence of the diagnostic G
A transitions at position 1138 of FGFR3 in PCR products of DNA extracts from short-statured individuals is authentic. If so, we were biased to believe authenticity, because the sequence data meet the phenotype. Thus, the unreliable molecular typing of ancient bones for ACH provides an example that even in a test system with phenotypic control proof of authenticity might be impossible.
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Burger, J., S. Hummel, B. Herrmann, and W. Henke 1999. DNA preservation: A microsatellite-DNA study on ancient skeletal remains. Electrophoresis 20:17221728.[CrossRef][ISI][Medline]
Cooper, D. N., and H. Youssoufian. 1988. The CpG dinucleotide and human genetic disease. Hum. Genet. 78:151155.[ISI][Medline]
Di Bernardo, G., S. Del Gaudio, M. Cammarota, U. Galderisi, A. Cascino, and M. Cipollaro. 2002. Enzymatic repair of selected cross-linked homoduplex molecules enhances nuclear gene rescue from Pompeii and Herculaneum remains. Nucleic Acids Res. 30:e16.
Frank, H., G. J. Nicholson, and E. Bayer. 1978. Enantiomer labelling, a method for the quantitative analysis of amino acids. J. Chromatogr. 167:187196.[CrossRef][Medline]
Gerhardt, J., and G. J. Nicholson. 1994. Unambiguous determination of the optical purity of peptides via GC-MS. Pp. 241243 in R. S. Hodges and J. A. Smith, eds. Proceedings of the thirteenth American peptide symposium. Peptides chemistry, structure and biology. Escom, Leiden, The Netherlands.
Handt, O., M. Krings, R. H. Ward, and S. Pääbo. 1996. The retrieval of ancient human DNA sequences. Am. J. Hum. Genet. 59:368376.[ISI][Medline]
Hansen, A. J., E. Willerslev, C. Wiuf, T. Mourier, and P. Arctander. 2001. Statistical evidence dor miscoding lesions in ancient DNA templates. Mol. Biol. Evol. 18:262265.
Hofreiter, M., V. Jaenicke, D. Serre, A. van Haeseler, and S. Pääbo. 2001. DNA sequences from multiple amplifications reveal artifacts induced by cytosine deamination in ancient DNA. Nucleic Acids Res. 29:47934799.
Höss, M., P. Jaruga, T. H. Zastawny, M. Dizdaroglu, and S. Pääbo. 1996. DNA damage and DNA sequence retrieval from ancient tissues. Nucleic Acids Res. 7:13041307.[CrossRef]
Keegan, K., D. E. Johnson, L. T. Williams, and M. J. Hayman. 1991. Isolation of an additional member of the fibroblast growth factor receptor family, FGFR-3. Proc. Natl. Acad. Sci. USA 88:10951099.[Abstract]
Krings, M., A. Stone, R. W. Schmitz, H. Krainitzki, M. Stoneking, and S. Pääbo. 1997. Neandertal DNA sequences and the origin of modern humans. Cell 90:1930.[ISI][Medline]
Kunichika, K., H. Hashimoto, and T. Imoto. 2002. Robustness of hen lysozyme monitored by random mutations. Protein Eng. 15:805809.[CrossRef][ISI][Medline]
Le Merrer, M., F. Rousseau, L. Legeai-Mallet et al. (11 co-authors). 1994. A gene for achondroplasia-hypochondroplasia maps to chromosome 4p. Nat. Genet. 6:318321.[ISI][Medline]
Lindahl, T. 1993. Instability and decay of the primary structure of DNA. Nature 362:709715.[CrossRef][ISI][Medline]
. 2000. Fossil DNA. Curr. Biol. 10:R616.[CrossRef][ISI][Medline]
Miller, S. A., D. D. Dykes, and H. F. Polesky. 1988. A simple salting out procedure for extracting DNA from nucleated cells. Nucleic Acids Res. 16:101102.
Nerlich, A. G., C. J. Haas, A. Zink, U. Szeimies, and H. G. Hagedorn. 1997. Molecular evidence for tuberculosis in an ancient Egyptian mummy. Lancet 350:1404.[Medline]
Nicholson, G. J., J. Tomiuk, A. Czarnetzki, L. Bachmann, and C. M. Pusch. 2002. Detection of bone glue treatment as a major source of contamination in ancient DNA analyses. Am. J. Physiol. Anthropol. 118:117120.[CrossRef][ISI][Medline]
Ovchinnikov, I. V., A. Gotherstrom, G. P. Romanova, V. M. Kharitonov, K. Liden, and W. Goodwin. 2000. Molecular analysis of Neanderthal DNA from the northern Caucasus. Nature 404:453454.[CrossRef][ISI][Medline]
Pääbo, S. 1989. Ancient DNA: extraction, characterization, molecular cloning, and enzymatic amplification. Proc. Natl. Acad. Sci. USA 86:19391943.[Abstract]
Pääbo, S., and A. C. Wilson. 1991. Miocene DNA sequencesa dream come true? Curr. Biol. 1:4546.[CrossRef]
Poinar, H. N., M. Höss, J. L. Bada, and S. Pääbo. 1996. Amino acid racemization and the preservation of ancient DNA. Science 272:864866.[Abstract]
Pu, W. T., and K. Struhl. 1992. Uracil interference, a rapid and general method for defining protein-DNA interactions involving the 5-methyl group of thymines: the GCN4-DNA complex. Nucleic Acids Res. 20:771775.[Abstract]
Pusch, C., and L. Bachmann. 2004. Spiking of contemporary human template DNA with ancient DNA extracts induces mutations under PCR and generates non-authentic sequences. Mol. Biol. Evol. 21:957964.
Pusch, C., H. Schmitt, and N. Blin. 1997. Increased cloning efficiency by cycle restriction-ligation CRL. Trends Genet 13:164.[CrossRef]
Pusch, C., and M. Scholz. 1997. DNA extraction from ancient human bones via enzymatic treatment. Trends Genet. 13:417.[CrossRef][ISI][Medline]
Pusch, C. M., I. Giddings, and M. Scholz, 1998. Repair of degraded duplex DNA from prehistoric samples using Escherichia coli DNA polymerase I and T4 DNA ligase. Nucleic Acids Res. 26:857859.
Rideout, W. M. 3rd, G. A. Coetzee, A. F. Olumi, and P. A. Jones. 1990. 5-Methylcytosine as an endogenous mutagen in the human LDL receptor and p53 genes. Science 249:12881290.[ISI][Medline]
Rousseau, F., J. Bonaventure, L. Legeai-Mallet, A. Pelet, J.-M. Rozet, P. Maroteaux, M. Le Merrer, and A. Munnich. 1994. Mutations in the gene encoding fibroblast growth factor receptor-3 in achondroplasia. Nature 371:252254.[CrossRef][ISI][Medline]
Salo, W. L., A. C. Aufderheide, J. Buikstra, and T. A. Holcomb. 1994. Identification of Mycobacterium tuberculosis DNA in a pre-Columbian Peruvian mummy. Proc. Natl. Acad. Sci. USA 91:20912094.[Abstract]
Shiang, R., L. M. Thompson, Y. Z. Zhu, D. M. Church, T. J. Fielder, M. Bocian, S. T. Winokur, and J. J. Wasmuth. 1994. Mutations in the transmembrane domain of FGFR3 cause the most common genetic form of dwarfism, achondroplasia. Cell 78:335342.[ISI][Medline]
Stanescu, R., V. Stanescu, and P. Maroteaux. 1990. Homozygous achondroplasia: morphologic and biochemical study of cartilage. Am. J. Med. Genet. 37:412421.[ISI][Medline]
Svetlov, V., and T. G. Cooper. 1998. Efficient PCR-based random mutagenesis of sub-genic (100bp) DNA fragments. Yeast 14:8991.[CrossRef][ISI][Medline]
Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 1997. The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25:48764882.
Thompson, L. M., S. Plummer, M. Schalling, M. R. Altherr, J. F. Gusella, D. E. Housman, and J. J. Wasmuth. 1991. A gene encoding a fibroblast growth factor receptor isolated from the Huntington disease gene region of human chromosome 4. Genomics 11:11331142.[ISI][Medline]
Velinov, M., S. A. Slaugenhaupt, I. Stoilov, C. I. Scott, Jr., J. F. Gusella, and P. Tsipouras. 1994. The gene for achondroplasia maps to the telomeric region of chromosome 4p. Nat. Genet. 6:314317.[ISI][Medline]
Wilkin, D. J., J. K. Szabo, R. Cameron et al. (12 co-authors) 1998. Mutations in fibroblast growth-factor receptor 3 in sporadic cases of achondroplasia occur exclusively on the paternally derived chromosome. Am. J. Hum. Genet. 63:711716.[CrossRef][ISI][Medline]