Finnish HLA studies confirm the increased risk conferred by HLA-B27 homozygosity in ankylosing spondylitis

Objective: To determine the influence of HLA-B27 homozygosity and HLA-DRB1 alleles in the susceptibility to, and severity of, ankylosing spondylitis in a Finnish population. Methods: 673 individuals from 261 families with ankylosing spondylitis were genotyped for HLA-DRB1 alleles and HLA-B27 heterozygosity/ homozygosity. The frequencies of HLA-B27 homozygotes in probands from these families were compared with the expected number of HLA-B27 homozygotes in controls under Hardy-Weinberg equilibrium (HWE). The effect of HLA-DRB1 alleles was assessed using a logistic regression procedure conditioned on HLA-B27 and case-control analysis. Results: HLA-B27 was detected in 93% of cases of ankylosing spondylitis. An overrepresentation of HLA-B27 homozygotes was noted in ankylosing spondylitis (11%) compared with the expected number of HLA-B27 homozygotes under HWE (4%) (odds ratio (OR) = 3.3 (95% confidence interval, 1.6 to 6.8), p = 0.002). HLA-B27 homozygosity was marginally associated with reduced BASDAI (HLA-B27 homozygotes, 4.5 (1.6); HLA-B27 heterozygotes, 5.4 (1.8) (mean (SD)), p = 0.05). Acute anterior uveitis (AAU) was present in significantly more HLA-B27 positive cases (50%) than HLA-B27 negative cases (16%) (OR = 5.4 (1.7 to 17), p<0.004). HLA-B27 positive cases had a lower average age of symptom onset (26.7 (8.0) years) compared with HLA-B27 negative cases (35.7 (11.2) years) (p<0.0001). Conclusions: HLA-627 homozygosity is associated with a moderately increased risk of ankylosing spondylitis compared with HLA-β27 heterozygosity. HLA-B27 positive cases had an earlier age of onset of ankylosing spondylitis than HLA-B27 negative cases and were more likely to develop AAU. HLA-DRB1 alleles may influence the age of symptom onset of ankylosing spondylitis.

A genome-wide screen for susceptibility loci in ankylosing spondylitis

Objective. To localize the regions containing genes that determine susceptibility to ankylosing spondylitis (AS). Methods. One hundred five white British families with 121 affected sibling pairs with AS were recruited, largely from the Royal National Hospital for Rheumatic Diseases AS database. A genome-wide linkage screen was undertaken using 254 highly polymorphic microsatellite markers from the Medical Research Council (UK) (MRC) set. The major histocompatibility complex (MHC) region was studied more intensively using 5 microsatellites lying within the HLA class III region and HLA-DRB1 typing. The Analyze package was used for 2-point analysis, and GeneHunter for multipoint analysis. Results. When only the MRC set was considered, 11 markers in 7 regions achieved a P value of ≤0.01. The maximum logarithm of odds score obtained was 3.8 (P = 1.4 x 10-5) using marker D6S273, which lies in the HLA class III region. A further marker used in mapping of the MHC class III region achieved a LOD score of 8.1 (P = 1 x 10-9). Nine of 118 affected sibling pairs (7.6%) did not share parental haplotypes identical by descent across the MHC, suggesting that only 31% of the susceptibility to AS is coded by genes linked to the MHC. The maximum non-MHC LOD score obtained was 2.6 (P = 0.0003) for marker D16S422. Conclusion. The results of this study confirm the strong linkage of the MHC with AS, and provide suggestive evidence regarding the presence and location of non-MHC genes influencing susceptibility to the disease.

The effect of HLA-DR on susceptibility to rheumatoid arthritis is influenced by the associated lymphotoxin α-tumor necrosis factor haplotype

Objective. HLA-DRB1, a major genetic determinant of susceptibility to rheumatoid arthritis (RA), is located within 1,000 kb of the gene encoding tumor necrosis factor (TNF). Because certain HLA-DRB1*04 subtypes increase susceptibility to RA, investigation of the role of the TNF gene is complicated by linkage disequilibrium (LD) between TNF and DRB1 alleles. By adequately controlling for this LD, we aimed to investigate the presence of additional major histocompatibility complex (MHC) susceptibility genes. Methods. We identified 274 HLA-DRB1*04-positive cases of RA and 271 HLA-DRB1*04-positive population controls. Each subject was typed for 6 single-nucleotide polymorphisms within a 4.5-kb region encompassing TNF and lymphotoxin a (LTA). LTA-TNF haplotypes in these unrelated individuals were determined using a combination of family data and the PHASE software program. Results. Significant differences in LTA-TNF haplotype frequencies were observed between different subtypes of HLA-DRB1*04. The LTA-TNF haplotypes observed were very restricted, with only 4 haplotypes constituting 81% of all haplotypes present. Among individuals carrying DRB1*0401, the LTA-TNF 2 haplotype was significantly underrepresented in cases compared with controls (odds ratio 0.5 [95% confidence interval 0.3-0.8], P = 0.007), while in those with DRB1*0404, the opposite effect was observed (P = 0.007). Conclusion. These findings suggest that the MHC contains genetic elements outside the LTA-TNF region that modify the effect of HLA-DRB1 on susceptibility to RA.

The effect of LRP5 polymorphisms on bone mineral density is apparent in childhood

Bone mass acquired during childhood is the primary determinant of adult bone mineral density (BMD) and osteoporosis risk. Bone accrual is subject to genetic influences. Activating and inactivating LRP5 gene mutations elicit extreme bone phenotypes, while more common LRP5 polymorphisms are associated with normal variation of BMD. Our aim was to test the hypothesis that LRP5 gene polymorphisms influence bone mass acquisition during childhood. The association between LRP5 gene polymorphisms and bone size and mineralization was examined in 819 unrelated British Caucasian children (n = 429 boys) aged 9 years. Height, weight, pubertal status (where available), total-body and spinal bone area, bone mineral content (BMC), BMD, and area-adjusted BMC (aBMC) were assessed. Dual-energy X-ray absorptiometry (DXA)-gene associations were assessed by linear regression, with adjustment for age, gender, pubertal status, and body size parameters. There were 140, 79, 12, and 2 girls who achieved Tanner stages I-IV, respectively, and 179 and 32 boys who achieved Tanner stages I and II, respectively. The rs2306862 (N740N) coding polymorphism in exon 10 of the LRP5 gene was associated with spinal BMD and aBMC (each P = 0.01) and total-body BMD and aBMC (P = 0.04 and 0.03, respectively). Adjusting for pubertal stage strengthened associations between this polymorphism and spinal BMD and aBMC (P = 0.01 and 0.002, respectively). Individuals homozygous for the T allele had greater spinal BMD and aBMC scores than those homozygous for the C allele. A dose effect was apparent as the mean spinal BMD and aBMC of heterozygous TC individuals were intermediate between those of their TT and CC counterparts. The N740N polymorphism in exon 10 of LRP5 was associated with spinal BMD and aBMC in pre- and early pubertal children. These results indicate that LRP5 influences volumetric bone density in childhood, possibly through effects on trabecular bone formation.

PTHR1 polymorphisms influence BMD variation through effects on the growing skeleton

We investigated whether polymorphisms in PTHR1 are associated with bone mineral density (BMD), to determine whether the association of this gene with BMD was due to effects on attainment of peak bone mass or effects on subsequent bone loss. The PTHR1 gene, including its 14 exons, their exon-intron boundaries, and 1,500 bp of its promoter region, was screened for polymorphisms by denaturing high-performance liquid chromatography (dHPLC) and sequencing in 36 osteoporotic cases. Eleven single-nucleotide polymorphisms (SNPs), one tetranucleotide repeat, and one tetranucleotide deletion were identified. A cohort of 634 families, including 1,236 men (39%) and 1,926 women (61%) ascertained with probands with low BMD (Z< -2.0) and the Children in Focus subset of the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort (785 unrelated individuals, mean age 118 months), were genotyped for the five most informative SNPs (minor allele frequency >5%) and the tetranucleotide repeat. In our osteoporosis families, association was noted between lumbar spine BMD and alleles of a known functional tetranucleotide repeat (U4) in the PTHR1 promoter region (P = 0.042) and between two and three marker haplotypes of PTHR1 polymorphisms with lumbar spine, femoral neck, and total hip BMD (P = 0.021-0.047). This association was restricted to the youngest tertile of the population (age 16-39 years, P = 0.013-0.048). A similar association was found for the ALSPAC cohort: two marker haplotypes of SNPs A48609T and C52813T were associated with height (P = 0.006) and total body less head BMD (P = 0.02), corrected for age and gender, confirming the family findings. These findings suggest a role for PTHR1 variation in determining peak BMD.

Major histocompatibility complex associations of ankylosing spondylitis are complex and involve further epistasis with ERAP1

Ankylosing spondylitis (AS) is a common, highly heritable, inflammatory arthritis for which HLA-B*27 is the major genetic risk factor, although its role in the aetiology of AS remains elusive. To better understand the genetic basis of the MHC susceptibility loci, we genotyped 7,264 MHC SNPs in 22,647 AS cases and controls of European descent. We impute SNPs, classical HLA alleles and amino-acid residues within HLA proteins, and tested these for association to AS status. Here we show that in addition to effects due to HLA-B*27 alleles, several other HLA-B alleles also affect susceptibility. After controlling for the associated haplotypes in HLA-B, we observe independent associations with variants in the HLA-A, HLA-DPB1 and HLA-DRB1 loci. We also demonstrate that the ERAP1 SNP rs30187 association is not restricted only to carriers of HLA-B*27 but also found in HLA-B*40:01 carriers independently of HLA-B*27 genotype.

Genetic insights into common pathways and complex relationships among iImmune-mediated diseases

Shared aetiopathogenic factors among immune-mediated diseases have long been suggested by their co-familiality and co-occurrence, and molecular support has been provided by analysis of human leukocyte antigen (HLA) haplotypes and genome-wide association studies. The interrelationships can now be better appreciated following the genotyping of large immune disease sample sets on a shared SNP array: the 'Immunochip'. Here, we systematically analyse loci shared among major immune-mediated diseases. This reveals that several diseases share multiple susceptibility loci, but there are many nuances. The most associated variant at a given locus frequently differs and, even when shared, the same allele often has opposite associations. Interestingly, risk alleles conferring the largest effect sizes are usually disease-specific. These factors help to explain why early evidence of extensive 'sharing' is not always reflected in epidemiological overlap. © 2013 Macmillan Publishers Limited. All rights reserved.

Common variants in the HLA-DRB1-HLA-DQA1 HLA class II region are associated with susceptibility to visceral leishmaniasis

To identify susceptibility loci for visceral leishmaniasis, we undertook genome-wide association studies in two populations: 989 cases and 1,089 controls from India and 357 cases in 308 Brazilian families (1,970 individuals). The HLA-DRB1-HLA-DQA1 locus was the only region to show strong evidence of association in both populations. Replication at this region was undertaken in a second Indian population comprising 941 cases and 990 controls, and combined analysis across the three cohorts for rs9271858 at this locus showed P combined = 2.76 × 10 -17 and odds ratio (OR) = 1.41, 95% confidence interval (CI) = 1.30-1.52. A conditional analysis provided evidence for multiple associations within the HLA-DRB1-HLA-DQA1 region, and a model in which risk differed between three groups of haplotypes better explained the signal and was significant in the Indian discovery and replication cohorts. In conclusion, the HLA-DRB1-HLA-DQA1 HLA class II region contributes to visceral leishmaniasis susceptibility in India and Brazil, suggesting shared genetic risk factors for visceral leishmaniasis that cross the epidemiological divides of geography and parasite species. © 2013 Nature America, Inc. All rights reserved.

Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci

Ankylosing spondylitis is a common, highly heritable inflammatory arthritis affecting primarily the spine and pelvis. In addition to HLA-B*27 alleles, 12 loci have previously been identified that are associated with ankylosing spondylitis in populations of European ancestry, and 2 associated loci have been identified in Asians. In this study, we used the Illumina Immunochip microarray to perform a case-control association study involving 10,619 individuals with ankylosing spondylitis (cases) and 15,145 controls. We identified 13 new risk loci and 12 additional ankylosing spondylitis-associated haplotypes at 11 loci. Two ankylosing spondylitis-associated regions have now been identified encoding four aminopeptidases that are involved in peptide processing before major histocompatibility complex (MHC) class I presentation. Protective variants at two of these loci are associated both with reduced aminopeptidase function and with MHC class I cell surface expression.

Genome-wide association study identifies five new schizophrenia loci

We examined the role of common genetic variation in schizophrenia in a genome-wide association study of substantial size: a stage 1 discovery sample of 21,856 individuals of European ancestry and a stage 2 replication sample of 29,839 independent subjects. The combined stage 1 and 2 analysis yielded genome-wide significant associations with schizophrenia for seven loci, five of which are new (1p21.3, 2q32.3, 8p23.2, 8q21.3 and 10q24.32-q24.33) and two of which have been previously implicated (6p21.32-p22.1 and 18q21.2). The strongest new finding (P = 1.6 × 10 -11) was with rs1625579 within an intron of a putative primary transcript for MIR137 (microRNA 137), a known regulator of neuronal development. Four other schizophrenia loci achieving genome-wide significance contain predicted targets of MIR137, suggesting MIR137-mediated dysregulation as a previously unknown etiologic mechanism in schizophrenia. In a joint analysis with a bipolar disorder sample (16,374 affected individuals and 14,044 controls), three loci reached genome-wide significance: CACNA1C (rs4765905, P = 7.0 × 10 -9), ANK3 (rs10994359, P = 2.5 × 10 -8) and the ITIH3-ITIH4 region (rs2239547, P = 7.8 × 10 -9).

The multiple sclerosis whole blood mRNA transcriptome and genetic associations indicate dysregulation of specific T cell pathways in pathogenesis

Multiple sclerosis (MS) is an autoimmune disease with a genetic component, caused at least in part by aberrant lymphocyte activity. The whole blood mRNA transcriptome was measured for 99 untreated MS patients: 43 primary progressive MS, 20 secondary progressive MS, 36 relapsing remitting MS and 45 age-matched healthy controls. The ANZgene Multiple Sclerosis Genetics Consortium genotyped more than 300 000 SNPs for 115 of these samples. Transcription from genes on translational regulation, oxidative phosphorylation, immune synapse and antigen presentation pathways was markedly increased in all forms of MS. Expression of genes tagging T cells was also upregulated (P < 10-12) in MS. A T cell gene signature predicts disease state with a concordance index of 0.79 with age and gender as co-variables, but the signature is not associated with clinical course or disability. The ANZgene genome wide association screen identified two novel regions with genome wide significance: one encoding the T cell co-stimulatory molecule, CD40; the other a region on chromosome 12q13-14. The CD40 haplotype associated with increased MS susceptibility has decreased gene expression in MS (P < 0.0007). The second MS susceptibility region includes 17 genes on 12q13-14 in tight linkage disequilibrium. Of these, only 13 are expressed in leukocytes, and of these the expression of one, FAM119B, is much lower in the susceptibility haplotype (P tdthomlt; 10-14). Overall, these data indicate dysregulation of T cells can be detected in the whole blood of untreated MS patients, and supports targeting of activated T cells in therapy for all forms of MS.

Partial epilepsy syndrome in a Gypsy family linked to 5q31.3-q32

PURPOSE The restricted genetic diversity and homogeneous molecular basis of Mendelian disorders in isolated founder populations have rarely been explored in epilepsy research. Our long-term goal is to explore the genetic basis of epilepsies in one such population, the Gypsies. The aim of this report is the clinical and genetic characterization of a Gypsy family with a partial epilepsy syndrome. METHODS Clinical information was collected using semistructured interviews with affected subjects and informants. At least one interictal electroencephalography (EEG) recording was performed for each patient and previous data obtained from records. Neuroimaging included structural magnetic resonance imaging (MRI). Linkage and haplotype analysis was performed using the Illumina IVb Linkage Panel, supplemented with highly informative microsatellites in linked regions and Affymetrix SNP 5.0 array data. RESULTS We observed an early-onset partial epilepsy syndrome with seizure semiology strongly suggestive of temporal lobe epilepsy (TLE), with mild intellectual deficit co-occurring in a large proportion of the patients. Psychiatric morbidity was common in the extended pedigree but did not cosegregate with epilepsy. Linkage analysis definitively excluded previously reported loci, and identified a novel locus on 5q31.3-q32 with an logarithm of the odds (LOD) score of 3 corresponding to the expected maximum in this family. DISCUSSION The syndrome can be classified as familial temporal lobe epilepsy (FTLE) or possibly a new syndrome with mild intellectual deficit. The linked 5q region does not contain any ion channel-encoding genes and is thus likely to contribute new knowledge about epilepsy pathogenesis. Identification of the mutation in this family and in additional patients will define the full phenotypic spectrum.

Fine mapping of the MHC Class III region demonstrates association of AIF1 and rheumatoid arthritis

Objective. The heritability of RA has been estimated to be ∼55%, of which the MHC contributes about one-third. HLA-DRB1 alleles are strongly associated with RA, but it is likely that significant non-DRB1 MHC genetic susceptibility factors are involved. Previously, we identified two three-marker haplotypes in a 106-kb region in the MHC class III region immediately centromeric to TNF, which are strongly associated with RA on HLA-DRB1*0404 haplotypes. In the present study, we aimed to refine these associations further using a combination of genotyping and gene expression studies. Methods. Thirty-nine nucleotide polymorphisms (SNPs) were genotyped in 95 DRB1*0404 carrying unrelated RA cases, 125 DRB1*0404 - carrying healthy controls and 87 parent-case trio RA families in which the affected child carried HLA-DRB1*04. Quantitative RT-PCR was used to assess the expression of the positional candidate MHC class III genes APOM, BAT2, BAT3, BAT4, BAT5, AIF1, C6orf47, CSNK2β and LY6G5C, and the housekeeper genes, hypoxanthine-guanine phosphoribosyltransferase (HPRT) and β2-microglobulin (B2M) in 31 RA cases and 21 ethnically, age- and sex-matched healthy controls. Synovial membrane specimens from RA, PsA and OA cases were stained by an indirect immunoperoxidase technique using a mouse-anti-human AIF1 monoclonal antibody. Results. Association was observed between RA and single markers or two marker haplotypes involving AIF1, BAT3 and CSNK. AIF1 was also significantly overexpressed in RA mononuclear cells (1.5- to 1.9-fold difference, P = 0.02 vs HPRT, P = 0.002 vs B2M). AIF1 protein was clearly expressed by synovial macrophages in all the inflammatory synovial samples in contrast to the non-inflammatory OA samples. Conclusions. The results of the genotyping and expression studies presented here suggest a role for AIF1 in both the aetiology and pathogenesis of RA.

Estrogen receptor α regulates area-adjusted bone mineral content in late pubertal girls

Context: Whether the action of estrogen in skeletal development depends on estrogen receptor α as encoded by the ESR1 gene is unknown. Objectives: The aim of this study was to establish whether the gain in area-adjusted bone mineral content (ABMC) in girls occurs in late puberty and to examine whether the magnitude of this gain is related to ESR1 polymorphisms. Design: We conducted a cross-sectional analysis. Setting: The study involved the Avon Longitudinal Study of Parents and Children (ALSPAC), a population-based prospective study. Participants: Participants included 3097 11-yr-olds with DNA samples, dual x-ray absorptiometry measurements, and pubertal stage information. Outcomes: Outcome measures included separate prespecified analyses in boys and girls of the relationship between ABMC derived from total body dual x-ray absorptiometry scans and Tanner stage and of the interaction between ABMC, Tanner stage, and ESR1 polymorphisms. Results: Total body less head and spinal ABMC were higher in girls in Tanner stages 4 and 5, compared with those in Tanner stages 1, 2, and 3. In contrast, height increased throughout puberty. No differences were observed in ABMC according to Tanner stage in boys. For rs2234693 (PvuII) and rs9340799 (XbaI) polymorphisms, differences in spinal ABMC in late puberty were 2-fold greater in girls who were homozygous for the C and G alleles, respectively (P = 0.001). For rs7757956, the difference in total body less head ABMC in late puberty was 50% less in individuals homozygous or heterozygous for the A allele (P = 0.006). Conclusions: Gains in ABMC in late pubertal girls are strongly associated with ESR1 polymorphisms, suggesting that estrogen contributes to this process via an estrogen receptor α-dependent pathway.

Non-B27 MHC associations of ankylosing spondylitis

Ankylosing spondylitis (AS) has been associated with human leukocyte antigen (HLA)-B27 for over 30 years; however, the mechanism of action has remained elusive. Although many studies have reported associations between AS and other genes in the major histocompatibility complex (MHC) in AS, no conclusive results have emerged. To investigate the contribution of non-B27 MHC genes to AS, a large cohort of AS families and controls were B27 typed and genotyped across the region. Interrogation of the data identified a region of 270kb, lying from 31952649 to 32221738 base pairs from the p-telomere of chromosome 6 and containing 23 genes, which is likely to include genes involved with susceptibility to AS.