The genetic basis of spondyloarthritis: SPARTAN/IGAS 2009

A joint meeting was held in July 2009 in Houston, Texas, of members of the Spondyloarthritis Research and Therapy Network (SPARTAN), founded in 2003 to promote research, education, and treatment of ankylosing spondylitis (AS) and related forms of spondyloarthritis (SpA), and members of International Genetics of AS (IGAS), founded in 2003 to encourage and coordinate studies internationally in the genetics of AS. The general topic was the genetic basis of SpA, with presentations on the future of human genetic studies; microbes, SpA, and innate immunity; susceptibility of AS to the major histocompatibility complex (MHC) and non-MHC; and individual discussions of the genetics of psoriasis and psoriatic arthritis, uveitis, inflammatory bowel disease, and enteropathic arthritis. Summaries of those discussions are presented.

Genome-wide association studies and musculoskeletal diseases

Bone and joint diseases are major causes of morbidity and mortality worldwide, and their prevalence is increasing as the average population age increases. Most common musculoskeletal diseases show significant heritability, and few have treatments that prevent disease or can induce true treatment-free, disease-free remission. Furthermore, despite valiant efforts of hypothesis-driven research, our understanding of the etiopathogenesis of these conditions is, with few exceptions, at best moderate. Therefore, there has been a long-standing interest in genetics research in musculoskeletal disease as a hypothesis-free method for investigating disease etiopathogenesis. Important contributions have been made through the identification of monogenic causes of disease, but the holy grail of human genetics research has been the identification of the genes responsible for common diseases. The development of genome-wide association (GWA) studies has revolutionized this field, and led to an explosion in the number of genes identified that are definitely involved in musculoskeletal disease pathogenesis. However, this approach will not identify all common disease genes, and although the current progress is exciting and proves the potential of this research discipline, other approaches will be required to identify many of the types of genetic variation likely to be involved.

Mutations in STIL, Encoding a Pericentriolar and Centrosomal Protein, Cause Primary Microcephaly

Primary microcephaly (MCPH) is an autosomal-recessive congenital disorder characterized by smaller-than-normal brain size and mental retardation. MCPH is genetically heterogeneous with six known loci: MCPH1-MCPH6. We report mapping of a novel locus, MCPH7, to chromosome 1p32.3-p33 between markers D1S2797 and D1S417, corresponding to a physical distance of 8.39 Mb. Heterogeneity analysis of 24 families previously excluded from linkage to the six known MCPH loci suggested linkage of five families (20.83%) to the MCPH7 locus. In addition, four families were excluded from linkage to the MCPH7 locus as well as all of the six previously known loci, whereas the remaining 15 families could not be conclusively excluded or included. The combined maximum two-point LOD score for the linked families was 5.96 at marker D1S386 at theta = 0.0. The combined multipoint LOD score was 6.97 between markers D1S2797 and D1S417. Previously, mutations in four genes, MCPH1, CDK5RAP2, ASPM, and CENPJ, that code for centrosomal proteins have been shown to cause this disorder. Three different homozygous mutations in STIL, which codes for a pericentriolar and centrosomal protein, were identified in patients from three of the five families linked to the MCPH7 locus; all are predicted to truncate the STIL protein. Further, another recently ascertained family was homozygous for the same mutation as one of the original families. There was no evidence for a common haplotype. These results suggest that the centrosome and its associated structures are important in the control of neurogenesis in the developing human brain.

Low birth weight in MZ twins discordant for birth weight is associated with shorter telomere length and lower IQ, but not anxiety/depression in later life

Shorter telomere length (TL) has found to be associated with lower birth weight and with lower cognitive ability and psychiatric disorders. However, the direction of causation of these associations and the extent to which they are genetically or environmentally mediated are unclear. Within-pair comparisons of monozygotic (MZ) and dizygotic (DZ) twins can throw light on these questions. We investigated correlations of within pair differences in telomere length, IQ, and anxiety/depression in an initial sample from Brisbane (242 MZ pairs, 245 DZ same sex (DZSS) pairs) and in replication samples from Amsterdam (514 MZ pairs, 233 DZSS pairs) and Melbourne (19 pairs selected for extreme high or low birth weight difference). Intra-pair differences of birth weight and telomere length were significantly correlated in MZ twins, but not in DZSS twins. Greater intra-pair differences of telomere length were observed in the 10% of MZ twins with the greatest difference in birth weight compared to the bottom 90% in both samples and also in the Melbourne sample. Intra-pair differences of telomere length and IQ, but not of TL and anxiety/depression, were correlated in MZ twins, and to a smaller extent in DZSS twins. Our findings suggest that the same prenatal effects that reduce birth weight also influence telomere length in MZ twins. The association between telomere length and IQ is partly driven by the same prenatal effects that decrease birth weight.

Independent replication and meta-analysis for endometriosis risk loci

Endometriosis is a complex disease that affects 6-10% of women in their reproductive years and 20-50% of women with infertility. Genome-wide and candidate-gene association studies for endometriosis have identified 10 independent risk loci, and of these, nine (rs7521902, rs13394619, rs4141819, rs6542095, rs1519761, rs7739264, rs12700667, rs1537377, and rs10859871) are polymorphic in European populations. Here we investigate the replication of nine SNP loci in 998 laparoscopically and histologically confirmed endometriosis cases and 783 disease-free controls from Belgium. SNPs rs7521902, rs13394619, and rs6542095 show nominally significant (p <.05) associations with endometriosis, while the directions of effect for seven SNPs are consistent with the original reports. Association of rs6542095 at the IL1A locus with 'All' (p =.066) and 'Grade-B' (p =.01) endometriosis is noteworthy because this is the first successful replication in an independent population. Meta-analysis with the published results yields genome-wide significant evidence for rs7521902, rs13394619, rs6542095, rs12700667, rs7739264, and rs1537377. Notably, three coding variants in GREB1 (near rs13394619) and CDKN2B-AS1 (near rs1537377) also showed nominally significant associations with endometriosis. Overall, this study provides important replication in a uniquely characterized independent population, and indicates that the majority of the original genome-wide association findings are not due to chance alone. © The Author(s) 2015.

Inference of the genetic architecture underlying BMI and height with the use of 20,240 sibling pairs

Evidence that complex traits are highly polygenic has been presented by population-based genome-wide association studies (GWASs) through the identification of many significant variants, as well as by family-based de novo sequencing studies indicating that several traits have a large mutational target size. Here, using a third study design, we show results consistent with extreme polygenicity for body mass index (BMI) and height. On a sample of 20,240 siblings (from 9,570 nuclear families), we used a within-family method to obtain narrow-sense heritability estimates of 0.42 (SE = 0.17, p = 0.01) and 0.69 (SE = 0.14, p = 6 x 10(-)(7)) for BMI and height, respectively, after adjusting for covariates. The genomic inflation factors from locus-specific linkage analysis were 1.69 (SE = 0.21, p = 0.04) for BMI and 2.18 (SE = 0.21, p = 2 x 10(-10)) for height. This inflation is free of confounding and congruent with polygenicity, consistent with observations of ever-increasing genomic-inflation factors from GWASs with large sample sizes, implying that those signals are due to true genetic signals across the genome rather than population stratification. We also demonstrate that the distribution of the observed test statistics is consistent with both rare and common variants underlying a polygenic architecture and that previous reports of linkage signals in complex traits are probably a consequence of polygenic architecture rather than the segregation of variants with large effects. The convergent empirical evidence from GWASs, de novo studies, and within-family segregation implies that family-based sequencing studies for complex traits require very large sample sizes because the effects of causal variants are small on average.

Genome-wide linkage and association analyses implicate FASN in predisposition to Uterine Leiomyomata

Uterine leiomyomata (UL), the most prevalent pelvic tumors in women of reproductive age, pose a major public health problem given their high frequency, associated morbidities, and most common indication for hysterectomies. A genetic component to UL predisposition is supported by analyses of ethnic predisposition, twin studies, and familial aggregation. A genome-wide SNP linkage panel was genotyped and analyzed in 261 white UL-affected sister-pair families from the Finding Genes for Fibroids study. Two significant linkage regions were detected in 10p11 (LOD = 4.15) and 3p21 (LOD = 3.73), and five additional linkage regions were identified with LOD scores > 2.00 in 2q37, 5p13, 11p15, 12q14, and 17q25. Genome-wide association studies were performed in two independent cohorts of white women, and a meta-analysis was conducted. One SNP (rs4247357) was identified with a p value (p = 3.05 x 10(-8)) that reached genome-wide significance (odds ratio = 1.299). The candidate SNP is under a linkage peak and in a block of linkage disequilibrium in 17q25.3, which spans fatty acid synthase (FASN), coiled-coil-domain-containing 57 (CCDC57), and solute-carrier family 16, member 3 (SLC16A3). By tissue microarray immunohistochemistry, we found elevated (3-fold) FAS levels in UL-affected tissue compared to matched myometrial tissue. FAS transcripts and/or protein levels are upregulated in various neoplasms and implicated in tumor cell survival. FASN represents the initial UL risk allele identified in white women by a genome-wide, unbiased approach and opens a path to management and potential therapeutic intervention.

Meta-analysis of telomere length in 19 713 subjects reveals high heritability, stronger maternal inheritance and a paternal age effect

Telomere length (TL) has been associated with aging and mortality, but individual differences are also influenced by genetic factors, with previous studies reporting heritability estimates ranging from 34 to 82%. Here we investigate the heritability, mode of inheritance and the influence of parental age at birth on TL in six large, independent cohort studies with a total of 19 713 participants. The meta-analysis estimate of TL heritability was 0.70 (95% CI 0.64–0.76) and is based on a pattern of results that is highly similar for twins and other family members. We observed a stronger mother–offspring (r=0.42; P-value=3.60 × 10−61) than father–offspring correlation (r=0.33; P-value=7.01 × 10−5), and a significant positive association with paternal age at offspring birth (β=0.005; P-value=7.01 × 10−5). Interestingly, a significant and quite substantial correlation in TL between spouses (r=0.25; P-value=2.82 × 10−30) was seen, which appeared stronger in older spouse pairs (mean age ≥55 years; r=0.31; P-value=4.27 × 10−23) than in younger pairs (mean age<55 years; r=0.20; P-value=3.24 × 10−10). In summary, we find a high and very consistent heritability estimate for TL, evidence for a maternal inheritance component and a positive association with paternal age.

Identification of the Meckel syndrome gene (MKS1) exposes a novel ciliopathy

Meckel syndrome (MKS, MIM 249000) is an autosomal recessive developmental disorder causing death in utero or shortly after birth. The hallmarks of the disease are cystic kidney dysplasia and fibrotic changes of the liver, occipital encephalocele with or without hydrocephalus and polydactyly. Other anomalies frequently seen in the patients are incomplete development of the male genitalia, club feet and cleft lip or palate. The clinical picture has been well characterized in the literature while the molecular pathology underlying the disease has remained unclear until now. In this study we identified the first MKS gene by utilizing the disease haplotypes in Finnish MKS families linked to the MKS1 locus on chromosome 17q23 (MKS1) locus. Subsequently, the genetic heterogeneity of MKS was established in the Finnish families. Mutations in at least four different genes can cause MKS. These genes have been mapped to the chromosomes 17q23 (MKS1), 11q13 (MKS2), 8q22 (MKS3) and 9q33 (MKS4). Two of these genes have been identified so far: The MKS1 gene (this work) and the MKS3 gene. The identified MKS1 gene was initially a novel human gene which is conserved among species. We found three different MKS mutations, one of them being the Finnish founder mutation. The information available from MKS1 orthologs in other species convinced us that the MKS1 gene is required for normal ciliogenesis. Defects of the cilial system in other human diseases and model organisms actually cause phenotypic features similar to those seen in MKS patients. The MKS3 (TMEM67) gene encodes a transmembrane protein and the gene maps to the syntenic Wpk locus in the rat, which is a model with polycystic kidney disease, agenesis of the corpus callosum and hydrocephalus. The available information from these two genes suggest that MKS1 would encode a structural component of the centriole required for normal ciliary functions, and MKS3 would be a transmembrane component most likely required for normal ciliary sensory signaling. The MKS4 locus was localized to chromosme 9q32-33 in this study by using an inbred Finnish family with two affected and two healthy children. This fourth locus contains TRIM32 gene, which is associated to another well characterized human ciliopathy, Bardet Biedl syndrome (BBS). Future studies should identify the MKS4 gene on chromosome 9q and confirm if there are more than two genes causing MKS Finnish families. The research on critical signaling pathways in organogenesis have shown that both Wnt and Hedgehog pathways are dependent on functional cilia. The MKS gene products will serve as excellent model molecules for more detailed studies of the functional role of cilia in organogenesis in more detail.

Genomic inflation factors under polygenic inheritance

Population structure, including population stratification and cryptic relatedness, can cause spurious associations in genome-wide association studies (GWAS). Usually, the scaled median or mean test statistic for association calculated from multiple single-nucleotide-polymorphisms across the genome is used to assess such effects, and 'genomic control' can be applied subsequently to adjust test statistics at individual loci by a genomic inflation factor. Published GWAS have clearly shown that there are many loci underlying genetic variation for a wide range of complex diseases and traits, implying that a substantial proportion of the genome should show inflation of the test statistic. Here, we show by theory, simulation and analysis of data that in the absence of population structure and other technical artefacts, but in the presence of polygenic inheritance, substantial genomic inflation is expected. Its magnitude depends on sample size, heritability, linkage disequilibrium structure and the number of causal variants. Our predictions are consistent with empirical observations on height in independent samples of ~4000 and ~133,000 individuals.

Genome-wide association study identifies a locus at 7p15.2 associated with endometriosis

Endometriosis is a common gynecological disease associated with pelvic pain and subfertility. We conducted a genome-wide association study (GWAS) in 3,194 individuals with surgically confirmed endometriosis (cases) and 7,060 controls from Australia and the UK. Polygenic predictive modeling showed significantly increased genetic loading among 1,364 cases with moderate to severe endometriosis. The strongest association signal was on 7p15.2 (rs12700667) for 'all' endometriosis (P = 2.6 x 10(-)(7), odds ratio (OR) = 1.22, 95% CI 1.13-1.32) and for moderate to severe disease (P = 1.5 x 10(-)(9), OR = 1.38, 95% CI 1.24-1.53). We replicated rs12700667 in an independent cohort from the United States of 2,392 self-reported, surgically confirmed endometriosis cases and 2,271 controls (P = 1.2 x 10(-)(3), OR = 1.17, 95% CI 1.06-1.28), resulting in a genome-wide significant P value of 1.4 x 10(-)(9) (OR = 1.20, 95% CI 1.13-1.27) for 'all' endometriosis in our combined datasets of 5,586 cases and 9,331 controls. rs12700667 is located in an intergenic region upstream of the plausible candidate genes NFE2L3 and HOXA10.

Variation in BMPR1B, TGFRB1 and BMPR2 and control of dizygotic twinning

Genes in the TGF9 signaling pathway play important roles in the regulation of ovarian follicle growth and ovulation rate. Mutations in three genes in this pathway, growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and the bone morphogenetic protein receptor B 1 (BMPRB1), influence dizygotic (DZ) twinning rates in sheep. To date, only variants in GDF9 and BMP15, but not their receptors transforming growth factor ss receptor 1 (TGFBR1), bone morphogenetic protein receptor 2 (BMPR2) and BMPR1B, have been investigated with respect to their roles in human DZ twinning. We screened for rare and novel variants in TGFBR1, BMPR2 and BMPR1B in mothers of dizygotic twins (MODZT) from twin-dense families, and assessed association between genotyped and imputed variants and DZ twinning in another large sample of MODZT. Three novel variants were found: a deep intronic variant in BMPR2, and one intronic and one non-synonymous exonic variant in BMPRB1 which would result in the replacement of glutamine by glutamic acid at amino acid position 294 (p.Gln294Glu). None of these variants were predicted to have major impacts on gene function. However, the p.Gln294Glu variant changes the same amino acid as a sheep BMPR1B functional variant and may have functional consequences. Six BMPR1B variants were marginally associated with DZ twinning in the larger case-control sample, but these were no longer significant once multiple testing was taken into account. Our results suggest that variation in the TGF9 signaling pathway type II receptors has limited effects on DZ twinning rates in humans.

Meta-analysis of genome-wide association for migraine in six population-based European cohorts

Migraine is a common neurological disorder with a genetically complex background. This paper describes a meta-analysis of genome-wide association (GWA) studies on migraine, performed by the Dutch-Icelandic migraine genetics (DICE) consortium, which brings together six population-based European migraine cohorts with a total sample size of 10,980 individuals (2446 cases and 8534 controls). A total of 32 SNPs showed marginal evidence for association at a P-value<10(-5). The best result was obtained for SNP rs9908234, which had a P-value of 8.00 x 10(-8). This top SNP is located in the nerve growth factor receptor (NGFR) gene. However, this SNP did not replicate in three cohorts from the Netherlands and Australia. Of the other 31 SNPs, 18 SNPs were tested in two replication cohorts, but none replicated. In addition, we explored previously identified candidate genes in the meta-analysis data set. This revealed a modest gene-based significant association between migraine and the metadherin (MTDH) gene, previously identified in the first clinic-based GWA study (GWAS) for migraine (Bonferroni-corrected gene-based P-value=0.026). This finding is consistent with the involvement of the glutamate pathway in migraine. Additional research is necessary to further confirm the involvement of glutamate.

Association between ORMDL3, IL1RL1 and a deletion on chromosome 17q21 with asthma risk in Australia

Genome-wide association studies followed by replication provide a powerful approach to map genetic risk factors for asthma. We sought to search for new variants associated with asthma and attempt to replicate the association with four loci reported previously (ORMDL3, PDE4D, DENND1B and IL1RL1). Genome-wide association analyses of individual single nucleotide polymorphisms (SNPs), rare copy number variants (CNVs) and overall CNV burden were carried out in 986 asthma cases and 1846 asthma-free controls from Australia. The most-associated locus in the SNP analysis was ORMDL3 (rs6503525, P = 4.8 x 10(-)(7)). Five other loci were associated with P < 10(-)(5), most notably the chemokine CXC motif ligand 14 (CXCL14) gene (rs31263, P = 7.8 x 10(-)(6)). We found no evidence for association with the specific risk variants reported recently for PDE4D, DENND1B and ILR1L1. However, a variant in IL1RL1 that is in low linkage disequilibrium with that reported previously was associated with asthma risk after accounting for all variants tested (rs10197862, gene wide P = 0.01). This association replicated convincingly in an independent cohort (P = 2.4 x 10(-)(4)). A 300-kb deletion on chromosome 17q21 was associated with asthma risk, but this did not reach experiment-wide significance. Asthma cases and controls had comparable CNV rates, length and number of genes affected by deletions or duplications. In conclusion, we confirm the association between asthma risk and variants in ORMDL3 and identify a novel risk variant in IL1RL1. Follow-up of the 17q21 deletion in larger cohorts is warranted.

A variant in LIN28B is associated with 2D:4D finger-length ratio, a putative retrospective biomarker of prenatal testosterone exposure

The ratio of the lengths of an individual's second to fourth digit (2D:4D) is commonly used as a noninvasive retrospective biomarker for prenatal androgen exposure. In order to identify the genetic determinants of 2D:4D, we applied a genome-wide association approach to 1507 11-year-old children from the Avon Longitudinal Study of Parents and Children (ALSPAC) in whom 2D:4D ratio had been measured, as well as a sample of 1382 12- to 16-year-olds from the Brisbane Adolescent Twin Study. A meta-analysis of the two scans identified a single variant in the LIN28B gene that was strongly associated with 2D:4D (rs314277: p = 4.1 x 10(-8)) and was subsequently independently replicated in an additional 3659 children from the ALSPAC cohort (p = 1.53 x 10(-6)). The minor allele of the rs314277 variant has previously been linked to increased height and delayed age at menarche, but in our study it was associated with increased 2D:4D in the direction opposite to that of previous reports on the correlation between 2D:4D and age at menarche. Our findings call into question the validity of 2D:4D as a simplistic retrospective biomarker for prenatal testosterone exposure.