Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428.

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Joint multi-ancestry and admixed GWAS reveals the complex genetics behind human cranial vault shape.

The cranial vault in humans is highly variable, clinically relevant, and heritable, yet its genetic architecture remains poorly understood. Here, we conduct a joint multi-ancestry and admixed multivariate genome-wide association study on 3D cranial vault shape extracted from magnetic resonance images of 6772 children from the ABCD study cohort yielding 30 genome-wide significant loci. Follow-up analyses indicate that these loci overlap with genomic risk loci for sagittal craniosynostosis, show elevated activity cranial neural crest cells, are enriched for processes related to skeletal development, and are shared with the face and brain. We present supporting evidence of regional localization for several of the identified genes based on expression patterns in the cranial vault bones of E15.5 mice. Overall, our study provides a comprehensive overview of the genetics underlying normal-range cranial vault shape and its relevance for understanding modern human craniofacial diversity and the etiology of congenital malformations.

ADHD and its neurocognitive substrates: A two sample Mendelian randomization study.

Attention-deficit/hyperactivity disorder (ADHD) is associated with a wide array of neural and cognitive features, and other psychiatric disorders, identified mainly through cross-sectional associations studies. However, it is unclear if the disorder is causally associated with these neurocognitive features. Here, we applied a two-sample bidirectional Mendelian randomization (MR) study to summary GWAS data to explore the presence and direction of a causal effect between ADHD and a range of neurocognitive features and other psychiatric disorders. The inverse variance weighted method was used in the main analysis, and two MR methods (MR-Egger, weighted median) were used for robustness checks. We found that genetic risk for ADHD was causally associated with a decreased area of lateral orbitofrontal cortex. Conversely, we found that brain volume and some features of intrinsic functional connectivity had causal effects on ADHD risk. Bidirectional causal links were found between ADHD and adult general intelligence, as well as depression and autistic spectrum disorders. Such work highlights the important ties between ADHD and general cognitive ability, and suggest some neural features, previously merely associated with the disorder, may play a causal role in its pathogenesis.

TBX3 and EFNA4 Variant in a Family with Ulnar-Mammary Syndrome and Sagittal Craniosynostosis.

Ulnar-mammary syndrome (UMS) is a rare, autosomal dominant disorder characterized by anomalies affecting the limbs, apocrine glands, dentition, and genital development. This syndrome is caused by haploinsufficiency in the T-Box3 gene (TBX3), with considerable variability in the clinical phenotype being observed even within families. We describe a one-year-old female with unilateral, postaxial polydactyly, and bilateral fifth fingernail duplication. Next-generation sequencing revealed a novel, likely pathogenic, variant predicted to affect the canonical splice site in intron 3 of the TBX3 gene (c.804 + 1G > A, IVS3 + 1G > A). This variant was inherited from the proband's father who was also diagnosed with UMS with the additional clinical finding of congenital, sagittal craniosynostosis. Subsequent whole genome analysis in the proband's father detected a variant in the EFNA4 gene (c.178C > T, p.His60Tyr), which has only been reported to be associated with sagittal craniosynostosis in one patient prior to this report but reported in other cranial suture synostosis. The findings in this family extend the genotypic spectrum of UMS, as well as the phenotypic spectrum of EFNA4-related craniosynostosis.

Targeted Sequencing of Candidate Regions Associated with Sagittal and Metopic Nonsyndromic Craniosynostosis.

Craniosynostosis (CS) is a major birth defect in which one or more skull sutures fuse prematurely. We previously performed a genome-wide association study (GWAS) for sagittal non-syndromic CS (sNCS), identifying associations downstream from BMP2 on 20p12.3 and intronic to BBS9 on 7p14.3; analyses of imputed variants in DLG1 on 3q29 were also genome-wide significant. We followed this work with a GWAS for metopic non-syndromic NCS (mNCS), discovering a significant association intronic to BMP7 on 20q13.31. In the current study, we sequenced the associated regions on 3q29, 7p14.3, and 20p12.3, including two candidate genes (BMP2 and BMPER) near some of these regions in 83 sNCS child-parent trios, and sequenced regions on 7p14.3 and 20q13.2-q13.32 in 80 mNCS child-parent trios. These child-parent trios were selected from the original GWAS cohorts if the probands carried at least one copy of the top associated GWAS variant (rs1884302 C allele for sNCS; rs6127972 T allele for mNCS). Many of the variants sequenced in these targeted regions are strongly predicted to be within binding sites for transcription factors involved in craniofacial development or bone morphogenesis. Variants enriched in more than one trio and predicted to be damaging to gene function are prioritized for functional studies.

Gain-of-function variants and overexpression of RUNX2 in patients with nonsyndromic midline craniosynostosis.

Craniosynostosis (CS), the premature fusion of one or more cranial sutures, is a relatively common congenital anomaly, occurring in 3-5 per 10,000 live births. Nonsyndromic CS (NCS) accounts for up to 80% of all CS cases, yet the genetic factors contributing to the disorder remain largely unknown. The RUNX2 gene, encoding a transcription factor critical for bone and skull development, is a well known CS candidate gene, as copy number variations of this gene locus have been found in patients with syndromic craniosynostosis. In the present study, we aimed to characterize RUNX2 to better understand its role in the genetic etiology and in the molecular mechanisms underlying midline suture ossification in NCS. We report four nonsynonymous variants, one intronic variant and one 18 bp in-frame deletion in RUNX2 not found in our study control population. Significant difference in allele frequency (AF) for the deletion variant RUNX2 p.Ala84-Ala89del (ClinVar 257,095; dbSNP rs11498192) was observed in our sagittal NCS cohort when compared to the general population (P = 1.28 × 10-6), suggesting a possible role in the etiology of NCS. Dual-luciferase assays showed that three of four tested RUNX2 variants conferred a gain-of-function effect on RUNX2, further suggesting their putative pathogenicity in the tested NCS cases. Downregulation of RUNX2 expression was observed in prematurely ossified midline sutures. Metopic sites showed significant downregulation of promoter 1-specific isoforms compared to sagittal sites. Suture-derived mesenchymal stromal cells showed an increased expression of RUNX2 over matched unfused suture derived cells. This demonstrates that RUNX2, and particularly the distal promoter 1-isoform group, are overexpressed in the osteogenic precursors within the pathological suture sites.

A genome-wide association study implicates the BMP7 locus as a risk factor for nonsyndromic metopic craniosynostosis.

Our previous genome-wide association study (GWAS) for sagittal nonsyndromic craniosynostosis (sNCS) provided important insights into the genetics of midline CS. In this study, we performed a GWAS for a second midline NCS, metopic NCS (mNCS), using 215 non-Hispanic white case-parent triads. We identified six variants with genome-wide significance (P ≤ 5 × 10-8): rs781716 (P = 4.71 × 10-9; odds ratio [OR] = 2.44) intronic to SPRY3; rs6127972 (P = 4.41 × 10-8; OR = 2.17) intronic to BMP7; rs62590971 (P = 6.22 × 10-9; OR = 0.34), located ~ 155 kb upstream from TGIF2LX; and rs2522623, rs2573826, and rs2754857, all intronic to PCDH11X (P = 1.76 × 10-8, OR = 0.45; P = 3.31 × 10-8, OR = 0.45; P = 1.09 × 10-8, OR = 0.44, respectively). We performed a replication study of these variants using an independent non-Hispanic white sample of 194 unrelated mNCS cases and 333 unaffected controls; only the association for rs6127972 (P = 0.004, OR = 1.45; meta-analysis P = 1.27 × 10-8, OR = 1.74) was replicated. Our meta-analysis examining single nucleotide polymorphisms common to both our mNCS and sNCS studies showed the strongest association for rs6127972 (P = 1.16 × 10-6). Our imputation analysis identified a linkage disequilibrium block encompassing rs6127972, which contained an enhancer overlapping a CTCF transcription factor binding site (chr20:55,798,821-55,798,917) that was significantly hypomethylated in mesenchymal stem cells derived from fused metopic compared to open sutures from the same probands. This study provides additional insights into genetic factors in midline CS.

Genetic associations with childhood brain growth, defined in two longitudinal cohorts.

Genome-wide association studies (GWASs) are unraveling the genetics of adult brain neuroanatomy as measured by cross-sectional anatomic magnetic resonance imaging (aMRI). However, the genetic mechanisms that shape childhood brain development are, as yet, largely unexplored. In this study we identify common genetic variants associated with childhood brain development as defined by longitudinal aMRI. Genome-wide single nucleotide polymorphism (SNP) data were determined in two cohorts: one enriched for attention-deficit/hyperactivity disorder (ADHD) (LONG cohort: 458 participants; 119 with ADHD) and the other from a population-based cohort (Generation R: 257 participants). The growth of the brain's major regions (cerebral cortex, white matter, basal ganglia, and cerebellum) and one region of interest (the right lateral prefrontal cortex) were defined on all individuals from two aMRIs, and a GWAS and a pathway analysis were performed. In addition, association between polygenic risk for ADHD and brain growth was determined for the LONG cohort. For white matter growth, GWAS meta-analysis identified a genome-wide significant intergenic SNP (rs12386571, P = 9.09 × 10-9 ), near AKR1B10. This gene is part of the aldo-keto reductase superfamily and shows neural expression. No enrichment of neural pathways was detected and polygenic risk for ADHD was not associated with the brain growth phenotypes in the LONG cohort that was enriched for the diagnosis of ADHD. The study illustrates the use of a novel brain growth phenotype defined in vivo for further study.

A variant associated with sagittal nonsyndromic craniosynostosis alters the regulatory function of a non-coding element.

Craniosynostosis presents either as a nonsyndromic congenital anomaly or as a finding in nearly 200 genetic syndromes. Our previous genome-wide association study of sagittal nonsyndromic craniosynostosis identified associations with variants downstream from BMP2 and intronic in BBS9. Because no coding variants in BMP2 were identified, we hypothesized that conserved non-coding regulatory elements may alter BMP2 expression. In order to identify and characterize noncoding regulatory elements near BMP2, two conserved noncoding regions near the associated region on chromosome 20 were tested for regulatory activity with a Renilla luciferase assay. For a 711 base pair noncoding fragment encompassing the most strongly associated variant, rs1884302, the luciferase assay showed that the risk allele (C) of rs1884302 drives higher expression of the reporter than the common allele (T). When this same DNA fragment was tested in zebrafish transgenesis studies, a strikingly different expression pattern of the green fluorescent reporter was observed depending on whether the transgenic fish had the risk (C) or the common (T) allele at rs1884302. The in vitro results suggest that altered BMP2 regulatory function at rs1884302 may contribute to the etiology of sagittal nonsyndromic craniosynostosis. The in vivo results indicate that differences in regulatory activity depend on the presence of a C or T allele at rs1884302.

Type I error rates of rare single nucleotide variants are inflated in tests of association with non-normally distributed traits using simple linear regression methods.

In this study, the effects of (a) the minor allele frequency of the single nucleotide variant (SNV), (b) the degree of departure from normality of the trait, and (c) the position of the SNVs on type I error rates were investigated in the Genetic Analysis Workshop (GAW) 19 whole exome sequence data. To test the distribution of the type I error rate, 5 simulated traits were considered: standard normal and gamma distributed traits; 2 transformed versions of the gamma trait (log10 and rank-based inverse normal transformations); and trait Q1 provided by GAW 19. Each trait was tested with 313,340 SNVs. Tests of association were performed with simple linear regression and average type I error rates were determined for minor allele frequency classes. Rare SNVs (minor allele frequency < 0.05) showed inflated type I error rates for non-normally distributed traits that increased as the minor allele frequency decreased. The inflation of average type I error rates increased as the significance threshold decreased. Normally distributed traits did not show inflated type I error rates with respect to the minor allele frequency for rare SNVs. There was no consistent effect of transformation on the uniformity of the distribution of the location of SNVs with a type I error.

Evaluation of IRX Genes and Conserved Noncoding Elements in a Region on 5p13.3 Linked to Families with Familial Idiopathic Scoliosis and Kyphosis.

Because of genetic heterogeneity present in idiopathic scoliosis, we previously defined clinical subsets (a priori) from a sample of families with idiopathic scoliosis to find genes involved with spinal curvature. Previous genome-wide linkage analysis of seven families with at least two individuals with kyphoscoliosis found linkage (P-value = 0.002) in a 3.5-Mb region on 5p13.3 containing only three known genes, IRX1, IRX2, and IRX4 In this study, the exons of IRX1, IRX2, and IRX4, the conserved noncoding elements in the region, and the exons of a nonprotein coding RNA, LOC285577, were sequenced. No functional sequence variants were identified. An intrafamilial test of association found several associated noncoding single nucleotide variants. The strongest association was with rs12517904 (P = 0.00004), located 6.5 kb downstream from IRX1 In one family, the genotypes of nine variants differed from the reference allele in all individuals with kyphoscoliosis, and two of three individuals with scoliosis, but did not differ from the reference allele in all other genotyped individuals. One of these variants, rs117273909, was located in a conserved noncoding region that functions as an enhancer in mice. To test whether the variant allele at rs117273909 had an effect on enhancer activity, zebrafish transgenesis was performed with overlapping fragments of 198 and 687 bp containing either the wild type or the variant allele. Our data suggests that this region acts as a regulatory element; however, its size and target gene(s) need to be identified to determine its role in idiopathic scoliosis.

Sequencing of the TBX6 Gene in Families with Familial Idiopathic Scoliosis.

STUDY DESIGN: A hypothesis-driven study was conducted in a familial cohort to determine the potential association between variants within the TBX6 gene and Familial Idiopathic Scoliosis (FIS). OBJECTIVE: To determine if variants within exons of the TBX6 gene segregate with the FIS phenotype within a sample of families with FIS. SUMMARY OF BACKGROUND DATA: Idiopathic Scoliosis (IS) is a structural curvature of the spine whose underlying genetic etiology has not been established. IS has been reported to occur at a higher rate than expected in family members of individuals with congenital scoliosis (CS), suggesting that the two diseases might have a shared etiology. The TBX6 gene on chromosome 16p, essential to somite development, has been associated with CS in a Chinese population. Previous studies have identified linkage to this locus in families with FIS, and specifically with rs8060511, located in an intron of the TBX6 gene. METHODS: Parent-offspring trios from 11 families (13 trios, 42 individuals) with FIS were selected for Sanger sequencing of the TBX6 gene. Trios were selected from a large population of families with FIS in which a genome-wide scan had resulted in linkage to 16p. RESULTS: Sequencing analyses of the subset of families resulted in the identification of five coding variants. Three of the five variants were novel; the remaining two variants were previously characterized and account for 90% of the observed variants in these trios. In all cases, there was no correlation between transmission of the TBX6 variant allele and FIS phenotype. However, an analysis of regulatory markers in osteoblasts showed that rs8060511 is in a putative enhancer element. CONCLUSIONS: Although this study did not identify any TBX6 coding variants that segregate with FIS, we identified a variant that is located in a potential TBX6 enhancer element. Therefore, further investigation of the region is needed.

CHD7 gene polymorphisms and familial idiopathic scoliosis.

STUDY DESIGN: Model-independent linkage analysis and tests of association were performed for 22 single nucleotide polymorphisms in the CHD7 gene in 244 families of European descent with familial idiopathic scoliosis (FIS). OBJECTIVE: To replicate an association between FIS and the CHD7 gene on 8q12.2 in an independent sample of families of European descent. SUMMARY OF BACKGROUND DATA: The CHD7 gene on chromosome 8, responsible for the CHARGE syndrome, was previously associated with FIS in an independent study that included 52 families of European descent. METHODS: Model-independent linkage analysis and intrafamilial tests of association were performed on the degree of lateral curvature considered as a qualitative trait (with thresholds of ≥10°, ≥15°, ≥20°, and ≥30°) and as a quantitative trait (degree of lateral curvature). Results from the tests of associations from this study and the previous study were combined in a weighted meta-analysis. RESULTS: No significant results (P < 0.01) were found for linkage analysis or tests of association between genetic variants of the CHD7 and FIS in this study, failing to replicate the findings from the previous study. Furthermore, no significant results (P < 0.01) were found from meta-analysis of the results from the tests of association from this sample and from the previous sample. CONCLUSION: No association between the 22 genotyped single nucleotide polymorphisms in the CHD7 gene and FIS within this study sample was found, failing to replicate the earlier findings. Further investigation of the CHD7 gene and its potential association to FIS may be required. LEVEL OF EVIDENCE: N/A.

Intra-familial tests of association between familial idiopathic scoliosis and linked regions on 9q31.3-q34.3 and 16p12.3-q22.2.

OBJECTIVE: Custom genotyping of markers in families with familial idiopathic scoliosis were used to fine-map candidate regions on chromosomes 9 and 16 in order to identify candidate genes that contribute to this disorder and prioritize them for next-generation sequence analysis. METHODS: Candidate regions on 9q and 16p-16q, previously identified as linked to familial idiopathic scoliosis in a study of 202 families, were genotyped with a high-density map of single nucleotide polymorphisms. Tests of linkage for fine-mapping and intra-familial tests of association, including tiled regression, were performed on scoliosis as both a qualitative and quantitative trait. RESULTS AND CONCLUSIONS: Nominally significant linkage results were found for markers in both candidate regions. Results from intra-familial tests of association and tiled regression corroborated the linkage findings and identified possible candidate genes suitable for follow-up with next-generation sequencing in these same families. Candidate genes that met our prioritization criteria included FAM129B and CERCAM on chromosome 9 and SYT1, GNAO1, and CDH3 on chromosome 16.

A genome-wide association study identifies susceptibility loci for nonsyndromic sagittal craniosynostosis near BMP2 and within BBS9.

Sagittal craniosynostosis is the most common form of craniosynostosis, affecting approximately one in 5,000 newborns. We conducted, to our knowledge, the first genome-wide association study for nonsyndromic sagittal craniosynostosis (sNSC) using 130 non-Hispanic case-parent trios of European ancestry (NHW). We found robust associations in a 120-kb region downstream of BMP2 flanked by rs1884302 (P = 1.13 × 10(-14), odds ratio (OR) = 4.58) and rs6140226 (P = 3.40 × 10(-11), OR = 0.24) and within a 167-kb region of BBS9 between rs10262453 (P = 1.61 × 10(-10), OR = 0.19) and rs17724206 (P = 1.50 × 10(-8), OR = 0.22). We replicated the associations to both loci (rs1884302, P = 4.39 × 10(-31) and rs10262453, P = 3.50 × 10(-14)) in an independent NHW population of 172 unrelated probands with sNSC and 548 controls. Both BMP2 and BBS9 are genes with roles in skeletal development that warrant functional studies to further understand the etiology of sNSC.

Old lessons learned anew: family-based methods for detecting genes responsible for quantitative and qualitative traits in the Genetic Analysis Workshop 17 mini-exome sequence data.

Family-based study designs are again becoming popular as new next-generation sequencing technologies make whole-exome and whole-genome sequencing projects economically and temporally feasible. Here we evaluate the statistical properties of linkage analyses and family-based tests of association for the Genetic Analysis Workshop 17 mini-exome sequence data. Based on our results, the linkage methods using relative pairs or nuclear families had low power, with the best results coming from variance components linkage analysis in nuclear families and Elston-Stewart model-based linkage analysis in extended pedigrees. For family-based tests of association, both ASSOC and ROMP performed well for genes with large effects, but ROMP had the advantage of not requiring parental genotypes in the analysis. For the linkage analyses we conclude that genome-wide significance levels appear to control type I error well but that "suggestive" significance levels do not. Methods that make use of the extended pedigrees are well powered to detect major loci segregating in the families even when there is substantial genetic heterogeneity and the trait is mainly polygenic. However, large numbers of such pedigrees will be necessary to detect all major loci. The family-based tests of association found the same major loci as the linkage analyses and detected low-frequency loci with moderate effect sizes, but control of type I error was not as stringent.

Identification of susceptibility loci for scoliosis in FIS families with triple curves.

The triple curve pattern (three lateral curvatures of equal severity) has been recognized as a distinct and unique clinical subtype of scoliosis. As part of a large study of familial idiopathic scoliosis (FIS), a subset of five families with a triple curve pattern (at least one member of each family having a triple curve) was evaluated to determine if this curve pattern was linked to any of the markers previously genotyped as part of the STRP-based previous linkage screen. Model independent linkage analysis (SIBPAL, v4.5) of the initial genomic screen identified candidate regions on chromosomes 6 and 10 when FIS was analyzed both as qualitative and quantitative traits in single- and multipoint linkage analyses. Additional fine mapping analyses of this subgroup with SNPs corroborated the findings in these regions (P < 0.001). These regions have been previously linked to FIS, however, this is the first time these regions have been implicated in a clinically well-defined subgroup and may suggest a unique genetic etiology for the formation of a triple curve.

Males with familial idiopathic scoliosis: a distinct phenotypic subgroup.

STUDY DESIGN: Statistical analysis of genomic screening and fine mapping data. OBJECTIVE: The goals of this study were to analyze a region on chromosome 17 and to identify specific genetic determinants within this region linked to familial idiopathic scoliosis (FIS) in a subgroup of families in which affected males have undergone surgery. SUMMARY OF BACKGROUND DATA: The high prevalence and variability of FIS is indicative of genetic heterogeneity. To localize genes related to scoliosis, identification of groups of families with common clinical characteristics is a strategy that reduces genetic heterogeneity. Two independent studies have implicated a region on chromosome 17 as related to FIS. METHODS: With approval of the Institutional Review Board, the initial study population consisted of 202 families (1198 individuals), each of which had 2 or more affected individuals; 17 of those families had an affected male who had surgery. Individuals underwent genomic screening and subsequent fine mapping. Results were obtained using model-independent linkage analysis, with scoliosis set as a qualitative and as a quantitative trait, as implemented in SIBPAL (S.A.G.E., v4.5). The level of significance was set at P < or = 0.05. RESULTS: The initial study population had significant results at markers d17s975 and d17s2196. Analyses of a subgroup of families with males having undergone surgery using a customized single nucleotide polymorphism panel resulted in increased significance of this region. CONCLUSION: The data confirm a previously reported genetic locus on chromosome 17 as statistically significant in the etiology of FIS within a subgroup of families in which an affected male had spinal surgery.

Idiopathic scoliosis: identification of candidate regions on chromosome 19p13.

STUDY DESIGN: We performed genomic screening, statistical linkage analysis, and fine mapping of 202 families with at least 2 individuals with idiopathic scoliosis. OBJECTIVE: To identify regions on chromosome 19p13 statistically linked to the phenotypic expression of idiopathic scoliosis. SUMMARY OF BACKGROUND DATA: Idiopathic scoliosis is a common structural curvature of the spine affecting otherwise healthy children. Presently, no clear consensus exists regarding the underlying abnormality or genetic determinants of this disease. METHODS: Model-independent linkage analysis of qualitative and quantitative traits related to scoliosis was used to screen genotyping data from 391 markers in 202 families (1198 individuals). Subsets of families with probands having a curve > or = 30 degrees were dichotomized based on the most likely mode of inheritance for each family (autosomal dominant or X-linked dominant). Fine mapping was performed to show linkage to candidate regions on chromosome 19. RESULTS: When the threshold of disease was set at a curvature of > or = 30 degrees, qualitative linkage analysis revealed significant results at 2 successive markers on chromosome 19. CONCLUSION: The data confirm a previously reported genetic locus on chromosome 19 as potentially significant in the etiology of idiopathic scoliosis.

Lack of association between the aggrecan gene and familial idiopathic scoliosis.

STUDY DESIGN: A study was conducted to determine the potential association between a specific aggrecan gene polymorphism and familial idiopathic scoliosis (FIS). OBJECTIVES: To determine the allelic distribution of the exon 12 polymorphism within a sample of families with FIS. SUMMARY OF BACKGROUND DATA: FIS is a structural curvature of the spine where the underlying genetic etiology has not been established. The aggrecan locus has been linked to multiple skeletal disorders. A polymorphism, within the aggrecan gene, was previously reported to be associated with curve severity in individuals with scoliosis. METHODS: Fifty-eight families with FIS were genotyped for the aggrecan exon 12 polymorphism using a polymerase chain reaction method. Model-independent sib-pair linkage analyses and tests of association were performed to analyze the genetic effects of the exon 12 polymorphism. RESULTS: Linkage analyses of a genomic screen performed on a subgroup of 48 families with a most likely to be X-linked dominant mode of inheritance of FIS showed marginally significant results on chromosome 15q25-26 (P < 0.05). The overall distribution of the alleles was consistent with previously reported literature; no evidence of association and marginal significance of linkage was found between the polymorphism and FIS or the degree of lateral curvature. CONCLUSIONS: Despite the negative association reported here, further investigation of the gene and its potential association to FIS is required.