- Research article
- Open Access
A major QTL controls susceptibility to spinal curvature in the curveback guppy
BMC Geneticsvolume 12, Article number: 16 (2011)
Understanding the genetic basis of heritable spinal curvature would benefit medicine and aquaculture. Heritable spinal curvature among otherwise healthy children (i.e. Idiopathic Scoliosis and Scheuermann kyphosis) accounts for more than 80% of all spinal curvatures and imposes a substantial healthcare cost through bracing, hospitalizations, surgery, and chronic back pain. In aquaculture, the prevalence of heritable spinal curvature can reach as high as 80% of a stock, and thus imposes a substantial cost through production losses. The genetic basis of heritable spinal curvature is unknown and so the objective of this work is to identify quantitative trait loci (QTL) affecting heritable spinal curvature in the curveback guppy. Prior work with curveback has demonstrated phenotypic parallels to human idiopathic-type scoliosis, suggesting shared biological pathways for the deformity.
A major effect QTL that acts in a recessive manner and accounts for curve susceptibility was detected in an initial mapping cross on LG 14. In a second cross, we confirmed this susceptibility locus and fine mapped it to a 5 cM region that explains 82.6% of the total phenotypic variance.
We identify a major QTL that controls susceptibility to curvature. This locus contains over 100 genes, including MTNR1B, a candidate gene for human idiopathic scoliosis. The identification of genes associated with heritable spinal curvature in the curveback guppy has the potential to elucidate the biological basis of spinal curvature among humans and economically important teleosts.
Idiopathic-type spinal curvature is a heritable condition that occurs during development without other structural malformations that would be indicative of a congenital defect. The deformity is observed among humans and teleosts, but not among quadrupedal animals . The phenotype is a consequence of genetic, biomechanical, and environmental factors that affect the spine during development.
Heritable idiopathic-type spinal curvature among humans
Idiopathic-type spinal curvature is the primary defect in the Idiopathic Scoliosis syndrome (IS) and Scheuermann kyphosis, and is also associated with other heritable disorders such as Prader-Willi syndrome and Turner syndrome. The extent to which etiological factors for curvature are shared between these developmental syndromes is unknown. Scheuermann kyphosis and IS may share genetic factors, based on familial clustering of the two pathologies , and/or common pathological processes [3, 4]. However, the syndromes differ morphologically; IS occurs in all three planes of the body and Scheuermann kyphosis is primarily a sagittal defect.
The global prevalence of IS among children is 0.5-10% (reviewed in ), and the prevalence for Scheuermann kyphosis is 4-10% . Neither the genetic basis nor the biological processes involved in the aetiology of either of these idiopathic-type spinal curvatures are known. The poor understanding of causative factors is due to trait complexity, including a high degree of phenotypic variability coinciding with growth, and a historic lack of a genetic and developmental animal model.
Recent work with the guppy (Poecilia reticulata) curveback lineage has demonstrated morphological and developmental similarities to the human IS syndrome and Scheuermann kyphosis [7–10]. Considering that teleosts share most developmental pathways, physiological mechanisms and organ systems with humans, it is likely that there are shared biological pathways for spinal curvature.
Heritable spinal curvature among teleosts
Spinal curvature is the most common morphological deformity among teleosts and can be caused by a variety of influences including nutritional, environmental, and genetic factors (reviewed in ). Heritable spinal curvature has been noted among laboratory (e.g. swordtail, guppy, medaka) and aquaculture (e.g. seabream, salmon, trout, sea bass) stocks. Among aquaculture stocks, spinal curvature can detract from the value of fish, and has prevented aquaculture from achieving optimal production [11–16]. The incidence of spinal curvature among aquaculture stocks ranges from 1-81%, depending on the species [16, 17]. Although environmental factors have been associated with spinal curvature, several lines of evidence indicate that a genetic component contributes to a significant portion of these deformities in teleosts [13, 16, 18–20]. However, no genes have as yet been identified that contribute to this deformity, even though such identification would enable aquaculture populations to be screened for alleles causing deformity, which could then be eliminated via marker assisted selection.
The curveback guppy lineage has a phenotype that has been extensively characterized, consisting of a heritable, idiopathic-type curve that has a propensity to increase in magnitude with growth . Guppy offspring are born with a fully ossified skeleton after ~3 weeks of gestation (guppies are live-bearers), and in the curveback lineage, curvature begins after birth and is generally stable by sexual maturity (approximately a month after birth). The curveback phenotype occurs primarily on the sagittal plane of the fish as an anterior lordosis and a posterior kyphosis, both of variable magnitude. In addition, some individuals exhibit coronal deviation of the spine (Figure1).
As described for human idiopathic-type scoliosis syndromes, there is extensive variability in the curveback population for curve magnitude, as well as for the propensity for a curve to progress, and for the rate of progression. Moreover, our initial survey of inheritance for curvature suggests that it is a complex trait under the control of multiple genes . Therefore, we used quantitative trait locus (QTL) mapping to identify genomic regions that are associated with curvature. Genes that are associated with idiopathic-type scoliosis in curveback will suggest candidate genes and biological pathways that may be involved in the human and other teleost curve phenotypes. Here we describe the mapping of chromosomal regions associated with spinal curvature in the curveback guppy.
Detection of QTL for idiopathic-type spinal curvature in curveback
QTL were detected using mapping progeny generated from Cross 1. From this cross 30 F1 offspring were produced, none of which were curved, indicating that curvature is a recessive trait. In the 129 backcross (BC) progeny, 48% had curves of variable magnitude whereas in the 286 F2 progeny, 22% showed curves of varying magnitude. The percent curved fish among BC and F2 progeny suggests that curvature is caused by a major recessive locus that is inherited in a Mendelian recessive manner. However, variation for curve magnitude among curved progeny suggests the presence of modifier alleles. The distribution of curve magnitude among BC and F2 offspring is given in Table1.
Of the 376 SNPs tested on BC offspring of Cross 1, 168 markers (44.6%) were polymorphic, 151 of which were assigned to 23 guppy linkage groups by Tripathi, et al., 2009  (Additional file 1). With the exception of LG 3, which has one marker, the number of markers per LG ranges from 3-14 and the size of intervals between markers ranges from 0.1-19.57 cM (Table2). Interval mapping identified one major effect QTL on LG14 that has a significant association with curvature (LOD 9.4; F = 57.66, chromosome-wide 1% significance threshold F = 9.03). The estimated location of this QTL is near Marker 0294 on LG14 (15.5 cM) . The allele substitution effect estimated to be 1.39 (0.17 SE). At Marker 0294 the alleles are fixed in the parents and 94% of the curved individuals and 25% of non-curved individuals were homozygous for the allele of the curved parent, suggesting that the QTL contains a gene or genes for curve predisposition. Linkage group 3 was omitted from the interval analysis because it only contained one marker, which is insufficient for interval mapping. One-way ANOVAs showed no association between the marker on LG3 and curve magnitude [F(1,82) = 1.31, p = 0.26], or between the markers that could not be assigned to a linkage group and curve magnitude (results not shown). We did not detect any other QTL with the single or two QTL model.
Confirmation and fine-mapping of QTL on LG 14
From Cross 2, 27 F1 were produced, none of which were curved. From the F1, four intercrosses were made. To localize the major QTL on LG14, we genotyped 175 F2 from Cross 2 with 12 guppy map markers between 0 and 24 cM , including 7 that were not informative in Cross 1 (Additional file 1). The parents of Cross 2 were homozygous for all markers. Interval analysis using genotype data from both crosses and the additional markers estimated the QTL location to be at 13 cM on LG14 (LOD 19.53), with a 95% confidence interval estimated by bootstrapping to be 9 cM  (Figure2). Calculation of the 95% confidence interval based on the two-LOD support interval yielded a slightly larger confidence interval, suggesting the region from 11 to 23 cM. Using the genotypes for Marker 0289 at 13 cM, this QTL explains approximately 82.6% of the phenotypic variation in Cross 2, and 54% of the difference between the parents.
No curved individuals were observed in the F1 populations of Cross 1 and Cross 2, suggesting that curve predisposition is recessive. Recessive inheritance was also confirmed by estimates of the additive and dominance effects from Cross 2, which were equal (1.08 in both cases). Therefore, curved individuals that were not homozygous for the allele from the curved parent at a given marker excluded that marker from the QTL region. In this way, we were able to narrow the QTL region further than by interval mapping (which does not assume that curved individuals are homozygous at the QTL). In Cross 2, Marker 0289 (at 13 cM) is homozygous for the curved allele in all curved individuals (and all non-curved individuals are either homozygous for the wild-type allele or are heterozygous), whereas flanking Markers 0635 and 0381(at 9.35 and 14.5 cM) are heterozygous in one or more curved individuals and so delimit the QTL region (Figure3). Thus, the QTL has been narrowed to a 5 cM region.
To identify genomic regions that are associated with idiopathic-type spinal curvature in the curveback guppy, we screened 84 backcross progeny with 367 guppy-specific markers to identify a major locus that accounts for curve susceptibility. In a second mapping cross between curveback and a more outbred lineage, the susceptibility locus was confirmed. The 57 non-curved F1 from mapping Crosses 1 and 2, estimates of allele effects from our mapping crosses, as well as evidence from the curveback pedigree (Gorman, et al., 2007: F1 of the original pedigree cross are non-curved), confirm that the inheritance of curve susceptibility in curveback is recessive.
Localization of the susceptibility QTL on LG14 using additional F2 mapping cross progeny from the second cross identified a 5 cM interval containing a marker whose genotypes explain 100% of the susceptibility for curvature (i.e., comparing curved vs. non-curved individuals), 82.6% of the phenotypic variation in curve severity (i.e, qualitative score 0-4), and 54% of the difference between the parents. The estimated difference in curvature between curved homozygotes and heterozygotes at the LG14 QTL is slightly higher in Cross 2 (2.16) than in Cross 1 (1.39), but this difference could be due to different genetic backgrounds in the two crosses.
Based on synteny between the guppy and medaka genomes, we estimate that there are 100-150 genes in the QTL interval. Synteny between the guppy and medaka genomes has been established by Tripathi et al., (2009), using the same guppy markers as this QTL study (Additional file 1). Moreover, to confirm synteny for LG 14, we blasted all guppy LG14 EST and BAC end clone sequences (available on: http://www.ncbi.nlm.nih.gov/pubmed/) to the medaka genome using the BLAST algorithm (Ensembl release 57: http://www.ensembl.org). We found that the gene order for 17 out of 20 markers is conserved between guppy and medaka on LG 14, demonstrating chromosome-wide synteny between these species. Therefore, we searched the medaka region corresponding to the 5 cM Guppy QTL for candidate genes named in human studies of IS. We found that the melatonin 1B receptor gene (MTNR1B), which has been implicated as a candidate for curve predisposition in human IS, is contained in the QTL .
The importance in identifying genes involved in spinal curvature
Spinal curvature is a prevalent and costly deformity among humans and teleosts. The estimated annual cost of treating children hospitalized with idiopathic scoliosis (IS) in the United States alone is over $3 billion. This cost estimate does not consider Scheuermann kyphosis, or adults with idiopathic-type spinal curvature who suffer from chronic back pain, contributing to the estimated $849 billion cost of treatment and lost wages associated with musculoskeletal disease . Among teleosts, spinal column deformities reduce total production in the aquaculture industry substantially . In contrast to humans, teleost curve phenotypes are less well characterized; although heritable curves are acknowledged to account for many cases among cultured stocks, whether these cases are from congenital defects or are idiopathic-type is often not known.
Despite the prevalence and impact of this type of deformity, the genetic architecture and specific genes involved are unknown. The current view is that human idiopathic-type scoliosis is a complex genetic disorder with multiple genes segregating in the population exhibiting complex genotype by environment interactions [25–30]. In aquaculture stocks, inheritance for spinal curvature has been described as Mendelian recessive or dominant, as well as polygenic, depending on how well the phenotype is characterized and what stock is considered (reviewed in , ). The guppy curveback phenotype has been extensively characterized so that the lineage can be applied as a model for understanding the biological context of heritable spinal curvature [1, 7–9]. Future studies can use approaches to map QTL affecting shape based on digital photos [31, 32], rather than the qualitative scale used in the present study. The identification of QTL in this study is a first step in understanding the genetics of this type of deformity and will lead to the identification of biological pathways associated with spinal integrity.
We detected a major QTL associated with idiopathic-type spinal curvature in the curveback guppy, and identify this as a susceptibility locus. The locus acts in a recessive manner and accounts for a large portion of phenotypic variance. Fine mapping of the susceptibility locus has identified a 5 cM region that contains over 100 genes among which is the human candidate for IS, MTNR1B. The identification of genes for spinal curvature in curveback will determine the molecular pathway(s) leading to spinal curvature in the guppy, and these findings will provide insights into the etiopathogenesis of spinal curvature in humans and other teleost species.
Study population and phenotypic evaluation
The guppy, P. reticulata, is a live-bearing teleost native to streams in northeast South America that has been used as a genetic model since the 1920's. The curveback lineage was established from laboratory guppies that are derived from a wild population caught in Central Cumaná, Venezuela and raised under standardized conditions since 2000. The curveback lineage originated from a curved male crossed to a non-curved, unrelated female in 2003, followed by full-sib mating and backcrossing. Breeding pairs are maintained in 4L plastic aquaria, and offspring are separated into individual 600 ml plastic containers after birth. Laboratory fish are kept under standardized conditions as described in Gorman et al., (2007), in compliance with protocols approved by the Simon Fraser University Animal Care Facility and the Canadian Council on Animal Care (#763B-05).
The guppy spine is visible without magnification. Fish were scored from the side (for sagittal curvature) and above (for coronal curvature) while in a plastic view tank 4"long × 2"wide × 3"high. Since all affected individuals exhibit sagittal curvature, the degree of lordosis is used as a standard for comparison of curve magnitude. Evaluation of curve magnitude is based on scores defined by Gorman et al., (2007) in which 0 denotes no curvature and scores of 1 though 4 reflect increasing curve magnitudes. Fish were scored for curve magnitude after birth and then once a month until 3 months of age (sexual maturity is at approximately one month). To maintain consistency and precision in scores, the same individual evaluated all fish.
At 3 months past birth fish were culled and photographed on a light table with a digital camera (Toshiba PDR-3310, NYC, USA) under 3X magnification. As described in Gorman and Breden (2010), all fish were positioned on their side and photographs were taken from above so that the camera looked down on the sagittal profile of the fish. The qualitative adult score 0-4 was used for analysis . Adult fish were frozen in buffered (EDTA) ethanol for genotyping. DNA was extracted from tissue using the DNeasy Blood and Tissue kit and protocol (Qiagen), and diluted as per specifications from Sequenom Inc. (see below).
Crosses used for mapping QTL
To detect QTL associated with curvature, mapping Cross 1 was generated from an extremely curved male (Pmale1) from the curveback lineage mated with a non-curved and unrelated female from a lineage having no prior incidence of spinal curvature (Pfemale1). The maternal lineage is a mixture of laboratory fish derived from Central Cumaná, Venezuela and a pet store fancy guppy stock, and has been inbred (via sib mating) for 7 generations as part of another experiment. According to the qualitative scale for curve magnitude (), Pmale1 had a magnitude of 4 (the highest curve magnitude), and Pfemale1 had a score of 0 (non-curved). From Cross 1 30 F1 offspring were produced, none of which were curved. Two non-curved F1 daughters were subsequently backcrossed to their curved father to obtain a total of 129 backcross (BC) mapping progeny. In addition, two non-curved F1 progeny were crossed to create 286 F2s that were used to estimate the distribution of curve magnitude, but were not used for genome-wide mapping.
For confirmation of the QTL detected in Cross1 and for fine mapping, additional offspring were generated from a second mapping cross. Cross 2 is an extremely curved curveback female (Pfemale2) crossed to a non-curved male from the QUL lineage (Pmale2). The QUL lineage is derived from the lower Quare River in Trinidad and was chosen as the outcross lineage because it exhibits a high rate of polymorphism in molecular markers compared to the Central Cumaná strain  from which the curveback lineage was derived. The Pfemale2 is derived from the same curveback lineage as Pmale1, but is the result of further generations of inbreeding. From Cross 2, 27 non-curved F1offspring were produced, from which we made four crosses to generate a total of 259 F2 mapping progeny.
Marker selection and QTL mapping
The guppy genome consists of 23 haploid chromosomes , with a genetic map distance estimated to be 899 cM, with no linkage group (LG) longer than 58 cM . The 2 parentals, 9 F1 males and females (including those used in production of BC and F2 offspring), and 53 non-curved and 31 curved BC offspring, were genotyped by Sequenom, Inc., using 189 EST-based (expressed sequence tag) and 178 genomic guppy-specific SNP (single nucleotide polymorphism) markers . These SNP markers are optimized for high-throughput genotyping by Sequenom Inc. using mass spectrometry. Markers were chosen based on their map positions published in Tripathi et al., (2009), with the goal of creating a map with 10 cM resolution.
To detect QTL that are associated with curvature in Cross 1, and to confirm the location of the major QTL in Cross 2, we performed interval mapping using GridQTL at http://www.gridqtl.org.uk/[35, 36]. Grid QTL evaluates the statistical support for the presence of a QTL at positions throughout the genome using least squares regression, which has been shown to be robust for use with non-normal and discrete characters, including binary traits [37, 38]. This program is able to analyze outbred populations; although the parental lines used in this study have been inbred for several generations, they are not homozygous at all marker loci. Because the Cross 1 parents were heterozygous at numerous markers, and because of the smaller sample sizes for the mapping crosses, we used genetic distances from Tripathi et al. (2009).
We first conducted analyses under the single locus model, and then under the two QTL model. For putative QTL, the substitution effect (difference between curved parent allele homozygotes and heterozygotes) was estimated in Cross 1 and additive and dominance effects were estimated in Cross 2. Sex was included as a fixed effect, and we also tested for an interaction between putative QTL and sex. To determine significance thresholds, one thousand permutations per chromosome (linkage group) were carried out for each trait to determine the distribution of the F-statistic under the null hypothesis that no QTL was segregating on that chromosome . For markers that could not be assigned to a linkage group for interval analysis, association with curvature was tested using individual marker genotypes and a one-way ANOVA. We used a significance threshold of 0.05 for all ANOVAs, which were conducted using JMP statistical software for MacOSX, Version 7.0, SAS Institute, INC., Cary, NC, USA.
Fine-mapping of a major effect QTL
The initial genome screen of Cross 1 detected a major effect QTL on LG14 (see Results) and so we genotyped 9 additional SNP markers in this region  by either direct sequencing of PCR products or PCR based RFLP  (Additional file 1 for marker numbers). However, further localization of the QTL was limited due to a lack of markers polymorphic between the parental strains in Cross 1. Using the F2s from Cross 2, we genotyped 175 individuals at markers on LG 14, including 7 markers that were not polymorphic in Cross 1  (see supplemental data), by either direct sequencing of PCR products or PCR- based RFLP .
Gorman KF, Breden F: Idiopathic-type scoliosis is not exclusive to bipedalism. Med. Hyp. 2009, 72 (3): 348-52. 10.1016/j.mehy.2008.09.052.
Axenovich TI, Zaidman AM, Zorkoltseva IV, Kalashnikova EV, Borodin PM: Segregation analysis of Scheuermann disease in ninety families from Siberia. Am J Med Genet. 2001, 100 (4): 275-9. 10.1002/ajmg.1290.
Roth M: Idiopathic scoliosis and Scheuermann's disease: essentially identical manifestations of neuro-vertebral growth disproportion. Radiol Diagn (Berl). 1981, 22 (3): 380-91.
Deacon P, Berkin CR, Dickson RA: Combined idiopathic kyphosis and scoliosis. An analysis of the lateral spinal curvatures associated with Scheuermann's disease. J Bone Joint Surg Br. 1985, 67 (2): 189-92.
Asher MA, Burton DC: Adolescent idiopathic scoliosis: natural history and long term treatment effects. Scoliosis. 2006, 1: 1-10. 10.1186/1748-7161-1-2.
Papagelopoulos PJ, Klassen RA, Peterson HA, Dekiituski MB: Surgical treatment of Scheuermann's disease with segtnental compression instrumentation. Clin Orthop Relat Res. 2001, 386: 139-149. 10.1097/00003086-200105000-00018.
Gorman KF, Tredwell SJ, Breden F: The mutant guppy syndrome curveback as a model for human heritable spinal curvature. Spine. 2007, 32 (7): 735-741. 10.1097/01.brs.0000259081.40354.e2.
Gorman KF, Breden F: Disproportionate body lengths correlate with idiopathic-type spinal curvature in the curveback guppy. Spine. 2010, 35 (5): 511-516. 10.1097/BRS.0b013e3181b52baf.
Gorman KF, Handrigan GR, Jin G, Wallis R, Breden F: Structural and micro-anatomical changes in vertebrae associated with idiopathic-type spinal curvature in the curveback guppy model. Scoliosis. 2010, 5: 10-10.1186/1748-7161-5-10.
Gorman KF, Breden F: Teleosts as models for human vertebral stability and deformity. Comparative Biochemistry and Physiology, Part C. 2007, 145: 28-38.
Fraser MR, Anderson TA, de Nys R: Ontogenic development of the spine and spinal deformities in larval barramundi (Lates calcarifer) culture. Aquaculture. 2004, 242: 697-711. 10.1016/j.aquaculture.2004.09.018.
Andrades JA, Becerra J, Fernándes-Llebrez P: Skeletal deformities in larval, juvenile and adult stages of cultured gilthead sea bream (Sparus aurata L.). Aquaculture. 1996, 141: 1-11. 10.1016/0044-8486(95)01226-5.
Afonso JM, Montero D, Robaina L, Astorga N, Izquierdo MS, Gines R: Association of a lordosis-scoliosis-kypohsis deformity in gilthead seabream (Sparus aurata) with family structure. Fish Physiology and Biochemistry. 2000, 22: 159-163. 10.1023/A:1007811702624.
Koumoundouros G, Maingot E, Divanach P, Kentouri M: Kyphosis in reared sea bass (Dicentrarchus labrax L.): ontogeny and effects on mortality. Aquaculture. 2002, 209 (1): 49-58. 10.1016/S0044-8486(01)00821-3.
Silverstone AM, Hammell L: Spinal deformities in farmed Atlantic salmon. Can Vet J. 2002, 43 (10): 782-784.
Bardon A, Vandeputte M, Dupont-Nivet M, Chavanne H, Haffray P, Vergnet A, Chatain B: What is the heritable component of spinal deformities in the European sea bass (Dicentrarchus labrax)?. Aquaculture. 2009, 294 (3-4): 194-201. 10.1016/j.aquaculture.2009.06.018.
Eissa AE, Moustafa M, El-Husseiny IN, S. Saeid S, Saleh O, Borhan T: Identification of some skeletal deformities in freshwater teleosts raised in Egyptian aquaculture. Chemoshpere. 2009, 77 (3): 419-425. 10.1016/j.chemosphere.2009.06.050.
Purdom CE: Genetics and Fish Breeding. 1993, Chapman & Hall, London
Divanach P, Boglione C, Menu B, Koumoundouros G, Kentouri M, Cataudella S: Abnormalities in finfish mariculture: an overview of the problem, causes and solutions. Seabass and Seabream Culture: Problems and Prospects. Edited by: Chatain, B., Saroglia, M., Ž . Sweetman, J., Lavens, P. 1996, European Aquaculture Society, Oostende, Belgium, 45-66.
Evans ML, Neff BD: Non-additive genetic effects contribute to larval spinal deformity in two populations of Chinook salmon (Oncorhynchus tshawytscha). Aquaculture. 2009, 296: 169-173. 10.1016/j.aquaculture.2009.08.018.
Tripathi N, Hoffmann M, Willing EM, Lanz C, Weigel D, Dreyer C: Genetic linkage map of the guppy, Poecilia reticulata, and quantitative trait loci analysis of male colour variation. Proc. R. Soc. B. 2009, 276: 2195-2208. 10.1098/rspb.2008.1930.
Visscher PM, Thompson R, Haley CS: Confidence intervals in QTL mapping by bootstrapping. Genetics. 1996, 143: 1013-1020.
Qiu XS, Tang NL, Yeung HY, Lee KM, Hung VW, Ng BK, Ma SL, Kwok RH, Qin L, Qiu Y, Cheng JC: Melatonin receptor 1B (MTNR1B) gene polymorphism is associated with the occurrence of adolescent idiopathic scoliosis. Spine. 2007, 32 (16): 1748-53. 10.1097/BRS.0b013e3180b9f0ff.
American Academy of Orthopedic Surgeons: United States Bone and Joint Decade: The Burden of Musculoskeletal Diseases in the United States. 2008, Chap. 3 pages 63-64 and Tables 3.3 and 3.4 Rosemont, IL
Miller NH, Marosy B, Justice CM, Novak SM, Tang EY, Boyce P, Pettengil J, Doheny KF, Pugh EW, Wilson AF: Linkage analysis of genetic loci for kyphoscoliosis on chromosomes 5p13, 13q13.3, and 13q32. American Journal of Medical Genetics Part A. 2006, 140A: 1059-1068. 10.1002/ajmg.a.31211.
Ogilvie JW, Braun J, Argyle V, Nelson L, Meade M, Ward K: The search for idiopathic scoliosis genes. Spine. 2006, 31 (6): 679-681. 10.1097/01.brs.0000202527.25356.90.
Cheung KMC, Wang T, Qui GX, Luk KDK: Recent advances in the aetiology of adolescent idiopathic scoliosis. International Orthopaedics (SICOT). 2008, 32 (6): 729-34. 10.1007/s00264-007-0393-y.
Cheng JCY, Tang NLS, Yeung HY, Miller N: Genetic association of complex traits using idiopathic scoliosis as an example. Clinical Orthopaedics and Related Research. 2007, 462: 38-44. 10.1097/BLO.0b013e3180d09dcc.
Xia CW, Qiu Y, Sun X, Qiu XS, Wang SF, Zhu ZZ, Zhu F: Vitamin D receptor gene polymorphisms in female adolescent idiopathic scoliosis patients. Zhonghua Yi Xue Za Zhi. 2007, 21 (87): 1465-9.
Wise CA, Gao X, Shoemaker S, Gordon D, Herring JA: Understanding genetic factors in idiopathic scoliosis, a complex disease of childhood. Curr Genomics. 2008, 9: 51-9. 10.2174/138920208783884874.
Klingenberg CP, Leamy LJ: Quantitative genetics of geometric shape in the mouse mandible. Evolution. 2001, 55: 2342-2352.
Fu GF, Berg A, Das K, Li JH, Li R, Wu RL: A statistical model for mapping morphological shape. Theoretical Biology and Medical Modeling. 2010, 7: 28-10.1186/1742-4682-7-28.
Alexander H, Breden F: Sexual isolation and extreme morphological divergence in the Cumaná guppy: a possible case of incipient speciation. Journal of Evolutionary Biology. 2004, 17: 1238-1254. 10.1111/j.1420-9101.2004.00788.x.
Winge O: The location of eighteen genes in Lebistes reticulates. J. Genetics. 1922, 12: 137-144. 10.1007/BF02983077.
Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM: QTL Express: mapping quantitative trait loci in simple and complex pedigrees. Bioinformatics. 2002, 18: 339-340. 10.1093/bioinformatics/18.2.339.
Seaton G, Hernandez J, Grunchec JA, White I, Allen J, De Koning DJ, Wei W, Berry D, Haley C, Knott S: GridQTL: A Grid Portal for QTL Mapping of Compute Intensive Datasets. Proceedings of the 8th World Congress on Genetics. Applied to Livestock Production, August 13-18. 2006, . Belo Horizonte, Brazil
Rebai A: Comparison of methods for regression interval mapping in QTL analysis with non-normal traits. Genet. Res. 1997, 69: 69-74. 10.1017/S0016672396002558. 1997
Visscher PM, Haley CS, Knott SA: Mapping QTLs for binary traits in backcross and F2 populations. Genet. Res. 1996, 68: 55-63. 10.1017/S0016672300033887.
Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138 (3): 963-971.
Ota M, Fukushima H, Kulski JK, Inoko H: Single nucleotide polymorphism detection by polymerase chain reaction-restriction fragment length polymorphism. Nature Protocols. 2007, 2: 2857-2864. 10.1038/nprot.2007.407.
We acknowledge the reviewers of this paper, and Corey Watson at Simon Fraser University for their insights and suggestions. This study received support from the Natural Sciences and Engineering Research Council of Canada, and the US National Institute of Health (#1R21AR053730-01A2). The SFU University Animal Care Committee approved this research (project #763B).
KFG and FB participated in the project conception and oversight. KFG was also in charge of mapping crosses and phenotyping. JP and RA were in charge of DNA extraction, and JP performed all PCR and RFLP procedures and sequence analysis for Cross 2. CD and DW provided the guppy specific markers. KFG, FB, and JKC performed statistical analyses. KFG and JKC drafted the manuscript and FB contributed to revisions for the final draft. All authors have read and approved the final manuscript.