Genetic Susceptibility Loci in Rheumatoid Arthritis Establish Transcriptional Regulatory Networks with Other Genes

Linkage studies have identified the human leukocyte antigen (HLA)-DRB1 as a putative rheumatoid arthritis (RA) susceptibility locus (SL). Nevertheless, it was estimated that its contribution was partial, suggesting that other non-HLA genes may play a role in RA susceptibility. To test this hypothesis, we conducted microarray transcription profiling of peripheral blood mononuclear cells in 15 RA patients and analyzed the data, using bioinformatics programs (significance analysis of microarrays method and GeneNetwork), which allowed us to determine the differentially expressed genes and to reconstruct transcriptional networks. The patients were grouped according to disease features or treatment with tumor necrosis factor blocker. Transcriptional networks that were reconstructed allowed us to identify the interactions occurring between RA SL and other genes, for example, HLA-DRB1 interacting with FNDC3A (fibronectin type III domain containing 3A). Given that fibronectin fragments can stimulate mediators of matrix and cartilage destruction in RA, this interaction is of special interest and may contribute to a clearer understanding of the functional role of HLA-DRB1 in RA pathogenesis.

Six novel susceptibility Loci for early-onset androgenetic alopecia and their unexpected association with common diseases.

Androgenetic alopecia (AGA) is a highly heritable condition and the most common form of hair loss in humans. Susceptibility loci have been described on the X chromosome and chromosome 20, but these loci explain a minority of its heritable variance. We conducted a large-scale meta-analysis of seven genome-wide association studies for early-onset AGA in 12,806 individuals of European ancestry. While replicating the two AGA loci on the X chromosome and chromosome 20, six novel susceptibility loci reached genome-wide significance (p = 2.62×10(-9)-1.01×10(-12)). Unexpectedly, we identified a risk allele at 17q21.31 that was recently associated with Parkinson's disease (PD) at a genome-wide significant level. We then tested the association between early-onset AGA and the risk of PD in a cross-sectional analysis of 568 PD cases and 7,664 controls. Early-onset AGA cases had significantly increased odds of subsequent PD (OR = 1.28, 95% confidence interval: 1.06-1.55, p = 8.9×10(-3)). Further, the AGA susceptibility alleles at the 17q21.31 locus are on the H1 haplotype, which is under negative selection in Europeans and has been linked to decreased fertility. Combining the risk alleles of six novel and two established susceptibility loci, we created a genotype risk score and tested its association with AGA in an additional sample. Individuals in the highest risk quartile of a genotype score had an approximately six-fold increased risk of early-onset AGA [odds ratio (OR) = 5.78, p = 1.4×10(-88)]. Our results highlight unexpected associations between early-onset AGA, Parkinson's disease, and decreased fertility, providing important insights into the pathophysiology of these conditions.

Genome-wide enrichment analysis between endometriosis and obesity-related traits reveals novel susceptibility loci.

Endometriosis is a chronic inflammatory condition in women that results in pelvic pain and subfertility, and has been associated with decreased body mass index (BMI). Genetic variants contributing to the heritable component have started to emerge from genome-wide association studies (GWAS), although the majority remain unknown. Unexpectedly, we observed an intergenic locus on 7p15.2 that was genome-wide significantly associated with both endometriosis and fat distribution (waist-to-hip ratio adjusted for BMI; WHRadjBMI) in an independent meta-GWAS of European ancestry individuals. This led us to investigate the potential overlap in genetic variants underlying the aetiology of endometriosis, WHRadjBMI and BMI using GWAS data. Our analyses demonstrated significant enrichment of common variants between fat distribution and endometriosis (P = 3.7 × 10(-3)), which was stronger when we restricted the investigation to more severe (Stage B) cases (P = 4.5 × 10(-4)). However, no genetic enrichment was observed between endometriosis and BMI (P = 0.79). In addition to 7p15.2, we identify four more variants with statistically significant evidence of involvement in both endometriosis and WHRadjBMI (in/near KIFAP3, CAB39L, WNT4, GRB14); two of these, KIFAP3 and CAB39L, are novel associations for both traits. KIFAP3, WNT4 and 7p15.2 are associated with the WNT signalling pathway; formal pathway analysis confirmed a statistically significant (P = 6.41 × 10(-4)) overrepresentation of shared associations in developmental processes/WNT signalling between the two traits. Our results demonstrate an example of potential biological pleiotropy that was hitherto unknown, and represent an opportunity for functional follow-up of loci and further cross-phenotype comparisons to assess how fat distribution and endometriosis pathogenesis research fields can inform each other.

Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci

We undertook a meta-analysis of six Crohn's disease genome-wide association studies (GWAS) comprising 6,333 affected individuals (cases) and 15,056 controls and followed up the top association signals in 15,694 cases, 14,026 controls and 414 parent-offspring trios. We identified 30 new susceptibility loci meeting genome-wide significance (P < 5 × 10 ? ? ). A series of in silico analyses highlighted particular genes within these loci and, together with manual curation, implicated functionally interesting candidate genes including SMAD3, ERAP2, IL10, IL2RA, TYK2, FUT2, DNMT3A, DENND1B, BACH2 and TAGAP. Combined with previously confirmed loci, these results identify 71 distinct loci with genome-wide significant evidence for association with Crohn's disease.

Meta-analysis identifies seven susceptibility loci involved in the atopic march.

Eczema often precedes the development of asthma in a disease course called the 'atopic march'. To unravel the genes underlying this characteristic pattern of allergic disease, we conduct a multi-stage genome-wide association study on infantile eczema followed by childhood asthma in 12 populations including 2,428 cases and 17,034 controls. Here we report two novel loci specific for the combined eczema plus asthma phenotype, which are associated with allergic disease for the first time; rs9357733 located in EFHC1 on chromosome 6p12.3 (OR 1.27; P=2.1 × 10(-8)) and rs993226 between TMTC2 and SLC6A15 on chromosome 12q21.3 (OR 1.58; P=5.3 × 10(-9)). Additional susceptibility loci identified at genome-wide significance are FLG (1q21.3), IL4/KIF3A (5q31.1), AP5B1/OVOL1 (11q13.1), C11orf30/LRRC32 (11q13.5) and IKZF3 (17q21). We show that predominantly eczema loci increase the risk for the atopic march. Our findings suggest that eczema may play an important role in the development of asthma after eczema.

Identification of novel susceptibility loci for inflammatory bowel disease on chromosomes 1p, 3q, and 4q: Evidence for epistasis between 1p and IBD1

The idiopathic inflammatory bowel diseases, Crohn’s disease (CD) and ulcerative colitis (UC), are chronic, frequently disabling diseases of the intestines. Segregation analyses, twin concordance, and ethnic differences in familial risks have established that CD and UC are complex, non-Mendelian, related genetic disorders. We performed a genome-wide screen using 377 autosomal markers, on 297 CD, UC, or mixed relative pairs from 174 families, 37% Ashkenazim. We observed evidence for linkage at 3q for all families (multipoint logarithm of the odds score (MLod) = 2.29, P = 5.7 × 10−4), with greatest significance for non-Ashkenazim Caucasians (MLod = 3.39, P = 3.92 × 10−5), and at chromosome 1p (MLod = 2.65, P = 2.4 × 10−4) for all families. In a limited subset of mixed families (containing one member with CD and another with UC), evidence for linkage was observed on chromosome 4q (MLod = 2.76, P = 1.9 × 10−4), especially among Ashkenazim. There was confirmatory evidence for a CD locus, overlapping IBD1, in the pericentromeric region of chromosome 16 (MLod = 1.69, P = 2.6 × 10−3), particularly among Ashkenazim (MLod = 1.51, P = 7.8 × 10−3); however, positive MLod scores were observed over a very broad region of chromosome 16. Furthermore, evidence for epistasis between IBD1 and chromosome 1p was observed. Thirteen additional loci demonstrated nominal (MLod > 1.0, P < 0.016) evidence for linkage. This screen provides strong evidence that there are several major susceptibility loci contributing to the genetic risk for CD and UC.

Genome scan for teratogen-induced clefting susceptibility loci in the mouse: Evidence of both allelic and locus heterogeneity distinguishing cleft lip and cleft palate

Nonsyndromic clefting of the lip and palate in humans has a highly complex etiology, with both multiple genetic loci and exposure to teratogens influencing susceptibility. Previous studies using mouse models have examined only very small portions of the genome. Here we report the findings of a genome-wide search for susceptibility genes for teratogen-induced clefting in the AXB and BXA set of recombinant inbred mouse strains. We compare results obtained using phenytoin (which induces cleft lip) and 6-aminonicotinamide (which induces cleft palate). We use a new statistical approach based on logistic regression suitable for these categorical data to identify several chromosomal regions as possible locations of clefting susceptibility loci, and we review candidate genes located within each region. Because cleft lip and cleft palate do not frequently co-aggregate in human families and because these structures arise semi-independently during development, these disorders are usually considered to be distinct in etiology. Our data, however, implicate several of the same chromosomal regions for both forms of clefting when teratogen-induced. Furthermore, different parental strain alleles are usually associated with clefting of the lip versus that of the palate (i.e., allelic heterogeneity). Because several other chromosomal regions are associated with only one form of clefting, locus heterogeneity also appears to be involved. Our findings in this mouse model suggest several priority areas for evaluation in human epidemiological studies.

Localization of non-Mhc collagen-induced arthritis susceptibility loci in DBA/1j mice

One approach to understanding common human diseases is to determine the genetic defects responsible for similar diseases in animal models and place those defective genes in their corresponding biochemical pathways. Our laboratory is working with an animal model for human rheumatoid arthritis called collagen-induced arthritis (CIA). We are particularly interested in determining the location of disease-predisposing loci. To that end, we performed experiments to localize susceptibility loci for CIA in an F2 cross between the highly susceptible mouse strain DBA/1j and the highly resistant mouse strain SWR/j. Specifically, a quantitative trait locus analysis was performed to localize regions of the mouse genome responsible for susceptibility/severity to CIA. One susceptibility locus, Cia1 in the major histocompatibility locus, had been identified previously. Two additional loci were detected in our analysis that contribute to CIA severity (Cia2, Cia3) on chromosomes 2 and 6. A third locus was detected that contributes to the age of onset of the disease. This locus (Cia4) was located on chromosome 2 and was linked to the same region as Cia2. Determining the identity of these loci may provide insights into the etiology of human rheumatoid arthritis.

Identification of rat susceptibility loci for adjuvant-oil-induced arthritis

One intradermal injection of incomplete Freund’s adjuvant-oil induces a T cell-mediated inflammatory joint disease in DA rats. Susceptibility genes for oil-induced arthritis (OIA) are located both within and outside the major histocompatibility complex (MHC, Oia1). We have searched for disease-linked non-MHC loci in an F2 intercross between DA rats and MHC-identical but arthritis-resistant LEW.1AV1 rats. A genome-wide scan with microsatellite markers revealed two major chromosome regions that control disease incidence and severity: Oia2 on chromosome 4 (P = 4 × 10−13) and Oia3 on chromosome 10 (P = 1 × 10−6). All animals homozygous for DA alleles at both loci developed severe arthritis, whereas all those homozygous for LEW.1AV1 alleles were resistant. These results have general implications for situations where nonspecific activation of the immune system (e.g., incomplete Freund’s adjuvant-oil) causes inflammation and disease, either alone or in conjunction with specific antigens. They may also provide clues to the etiology of inflammatory diseases in humans, including rheumatoid arthritis.

Localization of atherosclerosis susceptibility loci to chromosomes 4 and 6 using the Ldlr knockout mouse model

Atherosclerosis is a complex disease resulting from the interaction of multiple genes. We have used the Ldlr knockout mouse model in an interspecific genetic cross to map atherosclerosis susceptibility loci. A total of 174 (MOLF/Ei × B6.129S7-Ldlrtm1Her) × C57BL/6J-Ldlrtm1Her backcross mice, homozygous for the Ldlr null allele, were fed a Western-type diet for 3 months and then killed for quantification of aortic lesions. A genome scan was carried out by using DNA pools and microsatellite markers spaced at ≈18-centimorgan intervals. Quantitative trait locus analysis of individual backcross mice confirmed linkages to chromosomes 4 (Athsq1, logarithm of odds = 6.2) and 6 (Athsq2, logarithm of odds = 6.7). Athsq1 affected lesions in females only whereas Athsq2 affected both sexes. Among females, the loci accounted for ≈50% of the total variance of lesion area. The susceptible allele at Athsq1 was derived from the MOLF/Ei genome whereas the susceptible allele at Athsq2 was derived from C57BL/6J. Inheritance of susceptible alleles at both loci conferred a 2-fold difference in lesion area, suggesting an additive effect of Athsq1 and Athsq2. No associations were observed between the quantitative trait loci and levels of plasma total cholesterol, high density lipoprotein cholesterol, non-high density lipoprotein cholesterol, insulin, or body weight. We provide strong evidence for complex inheritance of atherosclerosis in mice with elevated plasma low density lipoprotein cholesterol and show a major influence of nonlipoprotein-related factors on disease susceptibility. Athsq1 and Athsq2 represent candidate susceptibility loci for human atherosclerosis, most likely residing on chromosomes 1p36–32 and 12p13–12, respectively.

Novel multiple sclerosis susceptibility loci implicated in epigenetic regulation.

We conducted a genome-wide association study (GWAS) on multiple sclerosis (MS) susceptibility in German cohorts with 4888 cases and 10,395 controls. In addition to associations within the major histocompatibility complex (MHC) region, 15 non-MHC loci reached genome-wide significance. Four of these loci are novel MS susceptibility loci. They map to the genes L3MBTL3, MAZ, ERG, and SHMT1. The lead variant at SHMT1 was replicated in an independent Sardinian cohort. Products of the genes L3MBTL3, MAZ, and ERG play important roles in immune cell regulation. SHMT1 encodes a serine hydroxymethyltransferase catalyzing the transfer of a carbon unit to the folate cycle. This reaction is required for regulation of methylation homeostasis, which is important for establishment and maintenance of epigenetic signatures. Our GWAS approach in a defined population with limited genetic substructure detected associations not found in larger, more heterogeneous cohorts, thus providing new clues regarding MS pathogenesis.

A simple method to localise pleiotropic susceptibility loci using univariate linkage analyses of correlated traits

Univariate linkage analysis is used routinely to localise genes for human complex traits. Often, many traits are analysed but the significance of linkage for each trait is not corrected for multiple trait testing, which increases the experiment-wise type-I error rate. In addition, univariate analyses do not realise the full power provided by multivariate data sets. Multivariate linkage is the ideal solution but it is computationally intensive, so genome-wide analysis and evaluation of empirical significance are often prohibitive. We describe two simple methods that efficiently alleviate these caveats by combining P-values from multiple univariate linkage analyses. The first method estimates empirical pointwise and genome-wide significance between one trait and one marker when multiple traits have been tested. It is as robust as an appropriate Bonferroni adjustment, with the advantage that no assumptions are required about the number of independent tests performed. The second method estimates the significance of linkage between multiple traits and one marker and, therefore, it can be used to localise regions that harbour pleiotropic quantitative trait loci (QTL). We show that this method has greater power than individual univariate analyses to detect a pleiotropic QTL across different situations. In addition, when traits are moderately correlated and the QTL influences all traits, it can outperform formal multivariate VC analysis. This approach is computationally feasible for any number of traits and was not affected by the residual correlation between traits. We illustrate the utility of our approach with a genome scan of three asthma traits measured in families with a twin proband.

Familial Paget's disease of bone: Nonlinkage to the PDB1 and PDB2 loci on chromosomes 6p and 18q in a large pedigree

Paget's disease of bone is a common condition characterized by bone pain, deformity, pathological fracture, and an increased incidence of osteosarcoma. Genetic factors play a role in the pathogenesis of Paget's disease but the molecular basis remains largely unknown. Susceptibility loci for Paget's disease of bone have been mapped to chromosome 6p21.3 (PDB1) and 18q121.1-q22 (PDB2) in different pedigrees, We have identified a large pedigree of over 250 individuals with 49 informative individuals affected with Paget's disease of bone; 31 of whom are available for genotypic analysis. The disease is inherited as an autosomal dominant trait in the pedigree with high penetrance by the sixth decade. Linkage analysis has been performed with markers at PDB1; these data show significant exclusion of linkage with log,, of the odds ratio (LOD) scores < -2 in this region. Linkage analysis of microsatellite markers from the PDB2 region has excluded linkage with this region, with a 30 cM exclusion region (LOD score < -2.0) centered on D18S42, These data confirm the genetic heterogeneity of Paget's disease of bone. Our hypothesis is that a novel susceptibility gene relevant to the pathogenesis of Paget's disease of bone lies elsewhere in the genome in the affected members of this pedigree and will be identified using a microsatellite genomewide scan followed by positional cloning.

Common variants in the ATP2B1 gene are associated with susceptibility to hypertension: the Japanese Millennium Genome Project.

Hypertension is one of the most common complex genetic disorders. We have described previously 38 single nucleotide polymorphisms (SNPs) with suggestive association with hypertension in Japanese individuals. In this study we extend our previous findings by analyzing a large sample of Japanese individuals (n=14 105) for the most associated SNPs. We also conducted replication analyses in Japanese of susceptibility loci for hypertension identified recently from genome-wide association studies of European ancestries. Association analysis revealed significant association of the ATP2B1 rs2070759 polymorphism with hypertension (P=5.3×10(-5); allelic odds ratio: 1.17 [95% CI: 1.09 to 1.26]). Additional SNPs in ATP2B1 were subsequently genotyped, and the most significant association was with rs11105378 (odds ratio: 1.31 [95% CI: 1.21 to 1.42]; P=4.1×10(-11)). Association of rs11105378 with hypertension was cross-validated by replication analysis with the Global Blood Pressure Genetics consortium data set (odds ratio: 1.13 [95% CI: 1.05 to 1.21]; P=5.9×10(-4)). Mean adjusted systolic blood pressure was highly significantly associated with the same SNP in a meta-analysis with individuals of European descent (P=1.4×10(-18)). ATP2B1 mRNA expression levels in umbilical artery smooth muscle cells were found to be significantly different among rs11105378 genotypes. Seven SNPs discovered in published genome-wide association studies were also genotyped in the Japanese population. In the combined analysis with replicated 3 genes, FGF5 rs1458038, CYP17A1, rs1004467, and CSK rs1378942, odds ratio of the highest risk group was 2.27 (95% CI: 1.65 to 3.12; P=4.6×10(-7)) compared with the lower risk group. In summary, this study confirmed common genetic variation in ATP2B1, as well as FGF5, CYP17A1, and CSK, to be associated with blood pressure levels and risk of hypertension.

In vitro whole-genome analysis identifies a susceptibility locus for HIV-1.

Advances in large-scale analysis of human genomic variability provide unprecedented opportunities to study the genetic basis of susceptibility to infectious agents. We report here the use of an in vitro system for the identification of a locus on HSA8q24.3 associated with cellular susceptibility to HIV-1. This locus was mapped through quantitative linkage analysis using cell lines from multigeneration families, validated in vitro, and followed up by two independent association studies in HIV-positive individuals. Single nucleotide polymorphism rs2572886, which is associated with cellular susceptibility to HIV-1 in lymphoblastoid B cells and in primary T cells, was also associated with accelerated disease progression in one of two cohorts of HIV-1-infected patients. Biological analysis suggests a role of the rs2572886 region in the regulation of the LY6 family of glycosyl-phosphatidyl-inositol (GPI)-anchored proteins. Genetic analysis of in vitro cellular phenotypes provides an attractive approach for the discovery of susceptibility loci to infectious agents.