A type I interferon transcriptional signature precedes autoimmunity in children genetically at risk for type 1 diabetes.

Diagnosis of the autoimmune disease type 1 diabetes (T1D) is preceded by the appearance of circulating autoantibodies to pancreatic islets. However, almost nothing is known about events leading to this islet autoimmunity. Previous epidemiological and genetic data have associated viral infections and antiviral type I interferon (IFN) immune response genes with T1D. Here, we first used DNA microarray analysis to identify IFN-ß-inducible genes in vitro and then used this set of genes to define an IFN-inducible transcriptional signature in peripheral blood mononuclear cells from a group of active systemic lupus erythematosus patients (n = 25). Using this predefined set of 225 IFN signature genes, we investigated the expression of the signature in cohorts of healthy controls (n = 87), patients with T1D (n = 64), and a large longitudinal birth cohort of children genetically predisposed to T1D (n = 109; 454 microarrayed samples). Expression of the IFN signature was increased in genetically predisposed children before the development of autoantibodies (P = 0.0012) but not in patients with established T1D. Upregulation of IFN-inducible genes was transient, temporally associated with a recent history of upper respiratory tract infections (P = 0.0064), and marked by increased expression of SIGLEC-1 (CD169), a lectin-like receptor expressed on CD14(+) monocytes. DNA variation in IFN-inducible genes altered T1D risk (P = 0.007), as exemplified by IFIH1, one of the genes in our IFN signature for which increased expression is a known risk factor for disease. These findings identify transient increased expression of type I IFN genes in preclinical diabetes as a risk factor for autoimmunity in children with a genetic predisposition to T1D.

Seven new loci associated with age-related macular degeneration.

Age-related macular degeneration (AMD) is a common cause of blindness in older individuals. To accelerate the understanding of AMD biology and help design new therapies, we executed a collaborative genome-wide association study, including >17,100 advanced AMD cases and >60,000 controls of European and Asian ancestry. We identified 19 loci associated at P < 5 × 10(-8). These loci show enrichment for genes involved in the regulation of complement activity, lipid metabolism, extracellular matrix remodeling and angiogenesis. Our results include seven loci with associations reaching P < 5 × 10(-8) for the first time, near the genes COL8A1-FILIP1L, IER3-DDR1, SLC16A8, TGFBR1, RAD51B, ADAMTS9 and B3GALTL. A genetic risk score combining SNP genotypes from all loci showed similar ability to distinguish cases and controls in all samples examined. Our findings provide new directions for biological, genetic and therapeutic studies of AMD.

Genetically distinct subsets within ANCA-associated vasculitis.

BACKGROUND: Antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis is a severe condition encompassing two major syndromes: granulomatosis with polyangiitis (formerly known as Wegener's granulomatosis) and microscopic polyangiitis. Its cause is unknown, and there is debate about whether it is a single disease entity and what role ANCA plays in its pathogenesis. We investigated its genetic basis. METHODS: A genomewide association study was performed in a discovery cohort of 1233 U.K. patients with ANCA-associated vasculitis and 5884 controls and was replicated in 1454 Northern European case patients and 1666 controls. Quality control, population stratification, and statistical analyses were performed according to standard criteria. RESULTS: We found both major-histocompatibility-complex (MHC) and non-MHC associations with ANCA-associated vasculitis and also that granulomatosis with polyangiitis and microscopic polyangiitis were genetically distinct. The strongest genetic associations were with the antigenic specificity of ANCA, not with the clinical syndrome. Anti-proteinase 3 ANCA was associated with HLA-DP and the genes encoding α(1)-antitrypsin (SERPINA1) and proteinase 3 (PRTN3) (P=6.2×10(-89), P=5.6×10(-12,) and P=2.6×10(-7), respectively). Anti-myeloperoxidase ANCA was associated with HLA-DQ (P=2.1×10(-8)). CONCLUSIONS: This study confirms that the pathogenesis of ANCA-associated vasculitis has a genetic component, shows genetic distinctions between granulomatosis with polyangiitis and microscopic polyangiitis that are associated with ANCA specificity, and suggests that the response against the autoantigen proteinase 3 is a central pathogenic feature of proteinase 3 ANCA-associated vasculitis. These data provide preliminary support for the concept that proteinase 3 ANCA-associated vasculitis and myeloperoxidase ANCA-associated vasculitis are distinct autoimmune syndromes. (Funded by the British Heart Foundation and others.).

Genome-wide association study of age-related macular degeneration identifies associated variants in the TNXB-FKBPL-NOTCH4 region of chromosome 6p21.3.

Age-related macular degeneration (AMD) is a leading cause of visual loss in Western populations. Susceptibility is influenced by age, environmental and genetic factors. Known genetic risk loci do not account for all the heritability. We therefore carried out a genome-wide association study of AMD in the UK population with 893 cases of advanced AMD and 2199 controls. This showed an association with the well-established AMD risk loci ARMS2 (age-related maculopathy susceptibility 2)-HTRA1 (HtrA serine peptidase 1) (P =2.7 × 10(-72)), CFH (complement factor H) (P =2.3 × 10(-47)), C2 (complement component 2)-CFB (complement factor B) (P =5.2 × 10(-9)), C3 (complement component 3) (P =2.2 × 10(-3)) and CFI (P =3.6 × 10(-3)) and with more recently reported risk loci at VEGFA (P =1.2 × 10(-3)) and LIPC (hepatic lipase) (P =0.04). Using a replication sample of 1411 advanced AMD cases and 1431 examined controls, we confirmed a novel association between AMD and single-nucleotide polymorphisms on chromosome 6p21.3 at TNXB (tenascin XB)-FKBPL (FK506 binding protein like) [rs12153855/rs9391734; discovery P =4.3 × 10(-7), replication P =3.0 × 10(-4), combined P =1.3 × 10(-9), odds ratio (OR) = 1.4, 95% confidence interval (CI) = 1.3-1.6] and the neighbouring gene NOTCH4 (Notch 4) (rs2071277; discovery P =3.2 × 10(-8), replication P =3.8 × 10(-5), combined P =2.0 × 10(-11), OR = 1.3, 95% CI = 1.2-1.4). These associations remained significant in conditional analyses which included the adjacent C2-CFB locus. TNXB, FKBPL and NOTCH4 are all plausible AMD susceptibility genes, but further research will be needed to identify the causal variants and determine whether any of these genes are involved in the pathogenesis of AMD.

Statistical colocalization of monocyte gene expression and genetic risk variants for type 1 diabetes.

One mechanism by which disease-associated DNA variation can alter disease risk is altering gene expression. However, linkage disequilibrium (LD) between variants, mostly single-nucleotide polymorphisms (SNPs), means it is not sufficient to show that a particular variant associates with both disease and expression, as there could be two distinct causal variants in LD. Here, we describe a formal statistical test of colocalization and apply it to type 1 diabetes (T1D)-associated regions identified mostly through genome-wide association studies and expression quantitative trait loci (eQTLs) discovered in a recently determined large monocyte expression data set from the Gutenberg Health Study (1370 individuals), with confirmation sought in an additional data set from the Cardiogenics Transcriptome Study (558 individuals). We excluded 39 out of 60 overlapping eQTLs in 49 T1D regions from possible colocalization and identified 21 coincident eQTLs, representing 21 genes in 14 distinct T1D regions. Our results reflect the importance of monocyte (and their derivatives, macrophage and dendritic cell) gene expression in human T1D and support the candidacy of several genes as causal factors in autoimmune pancreatic beta-cell destruction, including AFF3, CD226, CLECL1, DEXI, FKRP, PRKD2, RNLS, SMARCE1 and SUOX, in addition to the recently described GPR183 (EBI2) gene.

Age-related macular degeneration: the importance of family history as a risk factor.

BACKGROUND: Family history is considered a risk factor for age-related macular degeneration (AMD). With the advent of effective therapy for the disease, the importance of family history merits further investigation. This study quantifies the risk associated with family history, first, by a case-control study of reported family history and, second, by examining the siblings of AMD cases. METHODS: The authors recruited cases with advanced AMD, spouses and siblings. All subjects were carefully phenotyped. Clinical findings in the siblings were compared with spouses. Information about family history was collected. The ORs for reported family history of AMD were calculated. Analyses were adjusted for age, smoking and genotype. RESULTS: 495 AMD cases, 259 spouses and 171 siblings were recruited. The OR for AMD was 27.8 (CI 3.8 to 203.0; p=0.001) with a reported family history of an affected parent and 12.0 (CI 3.7 to 38.6; p<0.0001) with a history of an affected sibling. ORs adjusted for age and smoking were higher. Examination of siblings confirmed their increased risk with 23% affected by AMD and an OR of 10.8 (4.5 to 25.8; p<0.0001). Adjusting for age increased the OR to 16.1 (6.2 to 41.8). CONCLUSION: The risk of AMD is greatly increased by having an affected first-degree relative. Those at risk need to be made aware of this and AMD patients should advise siblings and children to seek prompt ophthalmological advice if they develop visual symptoms of distortion or reduced vision.

Genetic variation in complement regulators and susceptibility to age-related macular degeneration.

OBJECTIVES: Age-related macular degeneration (AMD) is the commonest cause of blindness in Western populations. Risk is influenced by age, genetic and environmental factors. Complement activation appears to be important in the pathogenesis and associations have been found between AMD and genetic variations in complement regulators such as complement factor H. We therefore investigated other complement regulators for association with AMD. METHODS: We carried out a case-control study to test for association between AMD and single nucleotide polymorphisms (SNPs) spanning the genes encoding complement factor P (CFP, properdin), CD46 (membrane cofactor protein, MCP), CD55 (decay accelerating factor, DAF) and CD59 (protectin). All cases and controls were examined by an ophthalmologist and had independent grading of fundus photographs to confirm their disease status. RESULTS: 20 SNPs were genotyped in 446 cases and 262 controls. For two SNPs with p-values approaching significance additional subjects were genotyped to increase the numbers to 622 cases and 359 controls. There was no evidence of association between AMD and any of the SNPs typed in CFP, CD46, CD55 or CD59. CONCLUSIONS: In a case-control sample that has shown the well established associations between AMD and variants in CFH, CFB and C3 there was absence of association with SNPs in CFP, CD46, CD55 and CD59. This suggests that these are not important susceptibility genes for AMD.

Copy number, linkage disequilibrium and disease association in the FCGR locus.

The response of a leukocyte to immune complexes (ICs) is modulated by receptors for the Fc region of IgG (FcgammaRs), and alterations in their affinity or function have been associated with risk of autoimmune diseases, including systemic lupus erythematosus (SLE). The low-affinity FcgammaR genomic locus is complex, containing regions of copy number variation (CNV) which can alter receptor expression and leukocyte responses to IgG. Combined paralogue ratio tests (PRTs) were used to distinguish three intervals within the FCGR locus which undergo CNV, and to determine FCGR gene copy number (CN). There were significant differences in FCGR3B and FCGR3A CNV profiles between Caucasian, East Asian and Kenyan populations. A previously noted association of low FCGR3B CN with SLE in Caucasians was supported [OR = 1.57 (1.08-2.27), P = 0.018], and replicated in Chinese [OR = 1.65 (1.25-2.18), P = 4 x 10(-4)]. There was no association of FCGR3B CNV with vasculitis, nor with malarial or bacterial infection. Linkage disequilibrium (LD) between multi-allelic FCGR3B CNV and SLE-associated SNPs in the FCGR locus was defined for the first time. Despite LD between FCGR3B CNV and a variant in FcgammaRIIB (I232T) which abolishes inhibitory function, both reduced CN of FCGR3B and homozygosity of the FcgammaRIIB-232T allele were individually strongly associated with SLE risk. Thus CN of FCGR3B, which controls IC responses and uptake by neutrophils, and variations in FCGR2B, which controls factors such as antibody production and macrophage activation, are important in SLE pathogenesis. Further interpretations of contributions to pathogenesis by FcgammaRs must be made in the context of LD involving CNV regions.

Genome-wide analysis of allelic expression imbalance in human primary cells by high-throughput transcriptome resequencing.

Many disease-associated variants identified by genome-wide association (GWA) studies are expected to regulate gene expression. Allele-specific expression (ASE) quantifies transcription from both haplotypes using individuals heterozygous at tested SNPs. We performed deep human transcriptome-wide resequencing (RNA-seq) for ASE analysis and expression quantitative trait locus discovery. We resequenced double poly(A)-selected RNA from primary CD4(+) T cells (n = 4 individuals, both activated and untreated conditions) and developed tools for paired-end RNA-seq alignment and ASE analysis. We generated an average of 20 million uniquely mapping 45 base reads per sample. We obtained sufficient read depth to test 1371 unique transcripts for ASE. Multiple biases inflate the false discovery rate which we estimate to be approximately 50% for random SNPs. However, after controlling for these biases and considering the subset of SNPs that pass HapMap QC, 4.6% of heterozygous SNP-sample pairs show evidence of imbalance (P < 0.001). We validated four findings by both bacterial cloning and Sanger sequencing assays. We also found convincing evidence for allelic imbalance at multiple reporter exonic SNPs in CD6 for two samples heterozygous at the multiple sclerosis-associated variant rs17824933, linking GWA findings with variation in gene expression. Finally, we show in CD4(+) T cells from a further individual that high-throughput sequencing of genomic DNA and RNA-seq following enrichment for targeted gene sequences by sequence capture methods offers an unbiased means to increase the read depth for transcripts of interest, and therefore a method to investigate the regulatory role of many disease-associated genetic variants.

The imprinted DLK1-MEG3 gene region on chromosome 14q32.2 alters susceptibility to type 1 diabetes.

Genome-wide association (GWA) studies to map common disease susceptibility loci have been hugely successful, with over 300 reproducibly associated loci reported to date. However, these studies have not yet provided convincing evidence for any susceptibility locus subject to parent-of-origin effects. Using imputation to extend existing GWA datasets, we have obtained robust evidence at rs941576 for paternally inherited risk of type 1 diabetes (T1D; ratio of allelic effects for paternal versus maternal transmissions = 0.75; 95% confidence interval (CI) = 0.71-0.79). This marker is in the imprinted region of chromosome 14q32.2, which contains the functional candidate gene DLK1. Our meta-analysis also provided support at genome-wide significance for a T1D locus at chromosome 19p13.2. The highest association was at marker rs2304256 (odds ratio (OR) = 0.86; 95%CI = 0.82-0.90) in the TYK2 gene, which has previously been associated with systemic lupus erythematosus and multiple sclerosis.

Sex chromosomes and genetic association studies.

Although the literature concerning statistical testing for genotype-phenotype association in family-based and population-based studies is very extensive, until recently the sex chromosomes have received little attention. Here it is shown that the X chromosome in particular presents special problems with respect to efficient analysis of mixed-sex population studies, and as a result of X inactivation. This paper reviews recent developments in approaching these problems.

Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes.

Type 1 diabetes (T1D) is a common autoimmune disorder that arises from the action of multiple genetic and environmental risk factors. We report the findings of a genome-wide association study of T1D, combined in a meta-analysis with two previously published studies. The total sample set included 7,514 cases and 9,045 reference samples. Forty-one distinct genomic locations provided evidence for association with T1D in the meta-analysis (P < 10(-6)). After excluding previously reported associations, we further tested 27 regions in an independent set of 4,267 cases, 4,463 controls and 2,319 affected sib-pair (ASP) families. Of these, 18 regions were replicated (P < 0.01; overall P < 5 × 10(-8)) and 4 additional regions provided nominal evidence of replication (P < 0.05). The many new candidate genes suggested by these results include IL10, IL19, IL20, GLIS3, CD69 and IL27.

Statistical independence of the colocalized association signals for type 1 diabetes and RPS26 gene expression on chromosome 12q13.

Following the recent success of genome-wide association studies in uncovering disease-associated genetic variants, the next challenge is to understand how these variants affect downstream pathways. The most proximal trait to a disease-associated variant, most commonly a single nucleotide polymorphism (SNP), is differential gene expression due to the cis effect of SNP alleles on transcription, translation, and/or splicing gene expression quantitative trait loci (eQTL). Several genome-wide SNP-gene expression association studies have already provided convincing evidence of widespread association of eQTLs. As a consequence, some eQTL associations are found in the same genomic region as a disease variant, either as a coincidence or a causal relationship. Cis-regulation of RPS26 gene expression and a type 1 diabetes (T1D) susceptibility locus have been colocalized to the 12q13 genomic region. A recent study has also suggested RPS26 as the most likely susceptibility gene for T1D in this genomic region. However, it is still not clear whether this colocalization is the result of chance alone or if RPS26 expression is directly correlated with T1D susceptibility, and therefore, potentially causal. Here, we derive and apply a statistical test of this hypothesis. We conclude that RPS26 expression is unlikely to be the molecular trait responsible for T1D susceptibility at this locus, at least not in a direct, linear connection.

Shared and distinct genetic variants in type 1 diabetes and celiac disease.

BACKGROUND: Two inflammatory disorders, type 1 diabetes and celiac disease, cosegregate in populations, suggesting a common genetic origin. Since both diseases are associated with the HLA class II genes on chromosome 6p21, we tested whether non-HLA loci are shared. METHODS: We evaluated the association between type 1 diabetes and eight loci related to the risk of celiac disease by genotyping and statistical analyses of DNA samples from 8064 patients with type 1 diabetes, 9339 control subjects, and 2828 families providing 3064 parent-child trios (consisting of an affected child and both biologic parents). We also investigated 18 loci associated with type 1 diabetes in 2560 patients with celiac disease and 9339 control subjects. RESULTS: Three celiac disease loci--RGS1 on chromosome 1q31, IL18RAP on chromosome 2q12, and TAGAP on chromosome 6q25--were associated with type 1 diabetes (P<1.00x10(-4)). The 32-bp insertion-deletion variant on chromosome 3p21 was newly identified as a type 1 diabetes locus (P=1.81x10(-8)) and was also associated with celiac disease, along with PTPN2 on chromosome 18p11 and CTLA4 on chromosome 2q33, bringing the total number of loci with evidence of a shared association to seven, including SH2B3 on chromosome 12q24. The effects of the IL18RAP and TAGAP alleles confer protection in type 1 diabetes and susceptibility in celiac disease. Loci with distinct effects in the two diseases included INS on chromosome 11p15, IL2RA on chromosome 10p15, and PTPN22 on chromosome 1p13 in type 1 diabetes and IL12A on 3q25 and LPP on 3q28 in celiac disease. CONCLUSIONS: A genetic susceptibility to both type 1 diabetes and celiac disease shares common alleles. These data suggest that common biologic mechanisms, such as autoimmunity-related tissue damage and intolerance to dietary antigens, may be etiologic features of both diseases.

PTPN22 Trp620 explains the association of chromosome 1p13 with type 1 diabetes and shows a statistical interaction with HLA class II genotypes.

OBJECTIVE: The disease association of the common 1858C>T Arg620Trp (rs2476601) nonsynonymous single nucleotide polymorphism (SNP) of protein tyrosine phosphatase; nonreceptor type 22 (PTPN22) on chromosome 1p13 has been confirmed in type 1 diabetes and also in other autoimmune diseases, including rheumatoid arthritis and Graves' disease. Some studies have reported additional associated SNPs independent of rs2476601/Trp(620), suggesting that it may not be the sole causal variant in the region and that the relative risk of rs2476601/Trp(620) is greater in lower risk by HLA class II genotypes than in the highest risk class II risk category. RESEARCH DESIGN AND METHODS: We resequenced PTPN22 and used these and other data to provide >150 SNPs to evaluate the association of the PTPN22 gene and its flanking chromosome region with type 1 diabetes in a minimum of 2,000 case subjects and 2,400 control subjects. RESULTS: Due to linkage disequilibrium, we were unable to distinguish between rs2476601/Trp(620) (P = 2.11 x10(-87)) and rs6679677 (P = 3.21 x10(-87)), an intergenic SNP between the genes putative homeodomain transcription factor 1 and round spermatid basic protein 1. None of the previously reported disease-associated SNPs proved to be independent of rs2476601/Trp(620). We did not detect any interaction with age at diagnosis or sex. However, we found that rs2476601/Trp(620) has a higher relative risk in type 1 diabetic case subjects carrying lower risk HLA class II genotypes than in those carrying higher risk ones (P = 1.36 x 10(-4) in a test of interaction). CONCLUSIONS: In our datasets, there was no evidence for allelic heterogeneity at the PTPN22 locus in type 1 diabetes, indicating that the SNP rs2476601/Trp(620) remains the best candidate in this chromosome region in European populations. The heterogeneity of rs2476601/Trp(620) disease risk by HLA class II genotype is consistent with previous studies, and the joint effect of the two loci is still greater in the high-risk group.

Localization of type 1 diabetes susceptibility to the MHC class I genes HLA-B and HLA-A.

The major histocompatibility complex (MHC) on chromosome 6 is associated with susceptibility to more common diseases than any other region of the human genome, including almost all disorders classified as autoimmune. In type 1 diabetes the major genetic susceptibility determinants have been mapped to the MHC class II genes HLA-DQB1 and HLA-DRB1 (refs 1-3), but these genes cannot completely explain the association between type 1 diabetes and the MHC region. Owing to the region's extreme gene density, the multiplicity of disease-associated alleles, strong associations between alleles, limited genotyping capability, and inadequate statistical approaches and sample sizes, which, and how many, loci within the MHC determine susceptibility remains unclear. Here, in several large type 1 diabetes data sets, we analyse a combined total of 1,729 polymorphisms, and apply statistical methods-recursive partitioning and regression-to pinpoint disease susceptibility to the MHC class I genes HLA-B and HLA-A (risk ratios >1.5; P(combined) = 2.01 x 10(-19) and 2.35 x 10(-13), respectively) in addition to the established associations of the MHC class II genes. Other loci with smaller and/or rarer effects might also be involved, but to find these, future searches must take into account both the HLA class II and class I genes and use even larger samples. Taken together with previous studies, we conclude that MHC-class-I-mediated events, principally involving HLA-B*39, contribute to the aetiology of type 1 diabetes.