Genome-wide association meta-analyses combining multiple risk phenotypes provide insights into the genetic architecture of cutaneous melanoma susceptibility

Most genetic susceptibility to cutaneous melanoma remains to be discovered. Meta-analysis genome-wide association study (GWAS) of 36,760 cases of melanoma (67% newly genotyped) and 375,188 controls identified 54 significant (P < 5 × 10-8) loci with 68 independent single nucleotide polymorphisms. Analysis of risk estimates across geographical regions and host factors suggests the acral melanoma subtype is uniquely unrelated to pigmentation. Combining this meta-analysis with GWAS of nevus count and hair color, and transcriptome association approaches, uncovered 31 potential secondary loci for a total of 85 cutaneous melanoma susceptibility loci. These findings provide insights into cutaneous melanoma genetic architecture, reinforcing the importance of nevogenesis, pigmentation and telomere maintenance, together with identifying potential new pathways for cutaneous melanoma pathogenesis.

Genetic susceptibility markers for a breast-colorectal cancer phenotype: Exploratory results from genome-wide association studies

BACKGROUND

Clustering of breast and colorectal cancer has been observed within some families and cannot be explained by chance or known high-risk mutations in major susceptibility genes. Potential shared genetic susceptibility between breast and colorectal cancer, not explained by high-penetrance genes, has been postulated. We hypothesized that yet undiscovered genetic variants predispose to a breast-colorectal cancer phenotype.

METHODS

To identify variants associated with a breast-colorectal cancer phenotype, we analyzed genome-wide association study (GWAS) data from cases and controls that met the following criteria: cases (n = 985) were women with breast cancer who had one or more first- or second-degree relatives with colorectal cancer, men/women with colorectal cancer who had one or more first- or second-degree relatives with breast cancer, and women diagnosed with both breast and colorectal cancer. Controls (n = 1769), were unrelated, breast and colorectal cancer-free, and age- and sex- frequency-matched to cases. After imputation, 6,220,060 variants were analyzed using the discovery set and variants associated with the breast-colorectal cancer phenotype at P<5.0E-04 (n = 549, at 60 loci) were analyzed for replication (n = 293 cases and 2,103 controls).

RESULTS

Multiple correlated SNPs in intron 1 of the ROBO1 gene were suggestively associated with the breast-colorectal cancer phenotype in the discovery and replication data (most significant; rs7430339, Pdiscovery = 1.2E-04; rs7429100, Preplication = 2.8E-03). In meta-analysis of the discovery and replication data, the most significant association remained at rs7429100 (P = 1.84E-06).

CONCLUSION

The results of this exploratory analysis did not find clear evidence for a susceptibility locus with a pleiotropic effect on hereditary breast and colorectal cancer risk, although the suggestive association of genetic variation in the region of ROBO1, a potential tumor suppressor gene, merits further investigation.

Gene-level association analysis of systemic sclerosis: A comparison of African-Americans and White populations

Gene-level analysis of ImmunoChip or genome-wide association studies (GWAS) data has not been previously reported for systemic sclerosis (SSc, scleroderma). The objective of this study was to analyze genetic susceptibility loci in SSc at the gene level and to determine if the detected associations were shared in African-American and White populations, using data from ImmunoChip and GWAS genotyping studies. The White sample included 1833 cases and 3466 controls (956 cases and 2741 controls from the US and 877 cases and 725 controls from Spain) and the African American sample, 291 cases and 260 controls. In both Whites and African Americans, we performed a gene-level analysis that integrates association statistics in a gene possibly harboring multiple SNPs with weak effect on disease risk, using Versatile Gene-based Association Study (VEGAS) software. The SNP-level analysis was performed using PLINK v.1.07. We identified 4 novel candidate genes (STAT1, FCGR2C, NIPSNAP3B, and SCT) significantly associated and 4 genes (SERBP1, PINX1, TMEM175 and EXOC2) suggestively associated with SSc in the gene level analysis in White patients. As an exploratory analysis we compared the results on Whites with those from African Americans. Of previously established susceptibility genes identified in Whites, only TNFAIP3 was significant at the nominal level (p = 6.13x10-3) in African Americans in the gene-level analysis of the ImmunoChip data. Among the top suggestive novel genes identified in Whites based on the ImmunoChip data, FCGR2C and PINX1 were only nominally significant in African Americans (p = 0.016 and p = 0.028, respectively), while among the top novel genes identified in the gene-level analysis in African Americans, UNC5C (p = 5.57x10-4) and CLEC16A (p = 0.0463) were also nominally significant in Whites. We also present the gene-level analysis of SSc clinical and autoantibody phenotypes among Whites. Our findings need to be validated by independent studies, particularly due to the limited sample size of African Americans.

Association Between Telomere Length and Risk of Cancer and Non-Neoplastic Diseases: A Mendelian Randomization Study

IMPORTANCE

The causal direction and magnitude of the association between telomere length and incidence of cancer and non-neoplastic diseases is uncertain owing to the susceptibility of observational studies to confounding and reverse causation.

OBJECTIVE

To conduct a Mendelian randomization study, using germline genetic variants as instrumental variables, to appraise the causal relevance of telomere length for risk of cancer and non-neoplastic diseases.

DATA SOURCES

Genomewide association studies (GWAS) published up to January 15, 2015.

STUDY SELECTION

GWAS of noncommunicable diseases that assayed germline genetic variation and did not select cohort or control participants on the basis of preexisting diseases. Of 163 GWAS of noncommunicable diseases identified, summary data from 103 were available.

DATA EXTRACTION AND SYNTHESIS

Summary association statistics for single nucleotide polymorphisms (SNPs) that are strongly associated with telomere length in the general population.

MAIN OUTCOMES AND MEASURES

Odds ratios (ORs) and 95% confidence intervals (CIs) for disease per standard deviation (SD) higher telomere length due to germline genetic variation.

RESULTS

Summary data were available for 35 cancers and 48 non-neoplastic diseases, corresponding to 420 081 cases (median cases, 2526 per disease) and 1 093 105 controls (median, 6789 per disease). Increased telomere length due to germline genetic variation was generally associated with increased risk for site-specific cancers. The strongest associations (ORs [95% CIs] per 1-SD change in genetically increased telomere length) were observed for glioma, 5.27 (3.15-8.81); serous low-malignant-potential ovarian cancer, 4.35 (2.39-7.94); lung adenocarcinoma, 3.19 (2.40-4.22); neuroblastoma, 2.98 (1.92-4.62); bladder cancer, 2.19 (1.32-3.66); melanoma, 1.87 (1.55-2.26); testicular cancer, 1.76 (1.02-3.04); kidney cancer, 1.55 (1.08-2.23); and endometrial cancer, 1.31 (1.07-1.61). Associations were stronger for rarer cancers and at tissue sites with lower rates of stem cell division. There was generally little evidence of association between genetically increased telomere length and risk of psychiatric, autoimmune, inflammatory, diabetic, and other non-neoplastic diseases, except for coronary heart disease (OR, 0.78 [95% CI, 0.67-0.90]), abdominal aortic aneurysm (OR, 0.63 [95% CI, 0.49-0.81]), celiac disease (OR, 0.42 [95% CI, 0.28-0.61]) and interstitial lung disease (OR, 0.09 [95% CI, 0.05-0.15]).

CONCLUSIONS AND RELEVANCE

It is likely that longer telomeres increase risk for several cancers but reduce risk for some non-neoplastic diseases, including cardiovascular diseases.

A comprehensive genome-wide analysis of melanoma Breslow thickness identifies interaction between CDC42 and SCIN genetic variants

Breslow thickness (BT) is a major prognostic factor of cutaneous melanoma (CM), the most fatal skin cancer. The genetic component of BT has only been explored by candidate gene studies with inconsistent results. Our objective was to uncover the genetic factors underlying BT using an hypothesis-free genome-wide approach. Our analysis strategy integrated a genome-wide association study (GWAS) of single nucleotide polymorphisms (SNPs) for BT followed by pathway analysis of GWAS outcomes using the gene-set enrichment analysis (GSEA) method and epistasis analysis within BT-associated pathways. This strategy was applied to two large CM datasets with Hapmap3-imputed SNP data: the French MELARISK study for discovery (966 cases) and the MD Anderson Cancer Center study (1,546 cases) for replication. While no marginal effect of individual SNPs was revealed through GWAS, three pathways, defined by gene ontology (GO) categories were significantly enriched in genes associated with BT (false discovery rate ≤5% in both studies): hormone activity, cytokine activity and myeloid cell differentiation. Epistasis analysis, within each significant GO, identified a statistically significant interaction between CDC42 and SCIN SNPs (pmeta-int =2.2 × 10(-6) , which met the overall multiple-testing corrected threshold of 2.5 × 10(-6) ). These two SNPs (and proxies) are strongly associated with CDC42 and SCIN gene expression levels and map to regulatory elements in skin cells. This interaction has important biological relevance since CDC42 and SCIN proteins have opposite effects in actin cytoskeleton organization and dynamics, a key mechanism underlying melanoma cell migration and invasion.

Mutations of HNRNPA0 and WIF1 predispose members of a large family to multiple cancers

We studied a large family that presented a strong familial susceptibility to multiple early onset cancers including prostate, breast, colon, and several other uncommon cancers. Through targeted gene, linkage, and whole genome sequencing analyses, we show that the presence of a variant in the regulatory region of HNRNPA0 associated with elevated cancer incidence in this family (Hazard ratio = 7.20, p = 0.0004). Whole genome sequencing identified a second rare protein changing mutation of WIF1 that interacted with the HNRNPA0 variant resulting in extremely high risk for cancer in carriers of mutations in both genes (p = 1.98 × 10(-13)). Analysis of downstream targets of the mutations in these two genes showed that the HNRNPA0 mutation affected expression patterns in the PI3 kinase and ERK/MAPK signaling pathways, while the WIF1 variant influenced expression of genes that play a role in NAD biosynthesis. This is a first report of variation in HNRNPA0 influencing common cancers or of a striking interaction between rare variants coexisting in an extended pedigree and jointly affecting cancer risk.

Genetic susceptibility loci of idiopathic interstitial pneumonia do not represent risk for systemic sclerosis: a case control study in Caucasian patients

BACKGROUND

Systemic sclerosis (SSc)-related interstitial lung disease (ILD) has phenotypic similarities to lung involvement in idiopathic interstitial pneumonia (IIP). We aimed to assess whether genetic susceptibility loci recently identified in the large IIP genome-wide association studies (GWASs) were also risk loci for SSc overall or severity of ILD in SSc.

METHODS

A total of 2571 SSc patients and 4500 healthy controls were investigated from the US discovery GWAS and additional US replication cohorts. Thirteen IIP-related selected single nucleotide polymorphisms (SNPs) were genotyped and analyzed for their association with SSc.

RESULTS

We found an association of SSc with the SNP rs6793295 in the LRRC34 gene (OR = 1.14, CI 95 % 1.03 to 1.25, p value = 0.009) and rs11191865 in the OBFC1 gene (OR = 1.09, CI 95 % 1.00 to 1.19, p value = 0.043) in the discovery cohort. Additionally, rs7934606 in MUC2 (OR = 1.24, CI 95 % 1.01 to 1.52, p value = 0.037) was associated with SSc-ILD defined by imaging. However, these associations failed to replicate in the validation cohort. Furthermore, SNPs rs2076295 in DSP (β = -2.29, CI 95 % -3.85 to -0.74, p value = 0.004) rs17690703 in SPPL2C (β = 2.04, CI 95 % 0.21 to 3.88, p value = 0.029) and rs1981997 in MAPT (β = 2.26, CI 95 % 0.35 to 4.17, p value = 0.02) were associated with percent predicted forced vital capacity (FVC%) even after adjusting for the anti-topoisomerase (ATA)-positive subset. However, these associations also did not replicate in the validation cohort.

CONCLUSIONS

Our results add new evidence that SSc and SSc-related ILD are genetically distinct from IIP, although they share phenotypic similarities.

Genome-wide meta-analysis identifies five new susceptibility loci for cutaneous malignant melanoma

Thirteen common susceptibility loci have been reproducibly associated with cutaneous malignant melanoma (CMM). We report the results of an international 2-stage meta-analysis of CMM genome-wide association studies (GWAS). This meta-analysis combines 11 GWAS (5 previously unpublished) and a further three stage 2 data sets, totaling 15,990 CMM cases and 26,409 controls. Five loci not previously associated with CMM risk reached genome-wide significance (P < 5 × 10(-8)), as did 2 previously reported but unreplicated loci and all 13 established loci. Newly associated SNPs fall within putative melanocyte regulatory elements, and bioinformatic and expression quantitative trait locus (eQTL) data highlight candidate genes in the associated regions, including one involved in telomere biology.

Integrated pathway and epistasis analysis reveals interactive effect of genetic variants at TERF1 and AFAP1L2 loci on melanoma risk

Genome-wide association studies (GWASs) have characterized 13 loci associated with melanoma, which only account for a small part of melanoma risk. To identify new genes with too small an effect to be detected individually but which collectively influence melanoma risk and/or show interactive effects, we used a two-step analysis strategy including pathway analysis of genome-wide SNP data, in a first step, and epistasis analysis within significant pathways, in a second step. Pathway analysis, using the gene-set enrichment analysis (GSEA) approach and the gene ontology (GO) database, was applied to the outcomes of MELARISK (3,976 subjects) and MDACC (2,827 subjects) GWASs. Cross-gene SNP-SNP interaction analysis within melanoma-associated GOs was performed using the INTERSNP software. Five GO categories were significantly enriched in genes associated with melanoma (false discovery rate ≤ 5% in both studies): response to light stimulus, regulation of mitotic cell cycle, induction of programmed cell death, cytokine activity and oxidative phosphorylation. Epistasis analysis, within each of the five significant GOs, showed significant evidence for interaction for one SNP pair at TERF1 and AFAP1L2 loci (pmeta-int  = 2.0 × 10(-7) , which met both the pathway and overall multiple-testing corrected thresholds that are equal to 9.8 × 10(-7) and 2.0 × 10(-7) , respectively) and suggestive evidence for another pair involving correlated SNPs at the same loci (pmeta-int  = 3.6 × 10(-6) ). This interaction has important biological relevance given the key role of TERF1 in telomere biology and the reported physical interaction between TERF1 and AFAP1L2 proteins. This finding brings a novel piece of evidence for the emerging role of telomere dysfunction into melanoma development.

Rare variants of large effect in BRCA2 and CHEK2 affect risk of lung cancer

We conducted imputation to the 1000 Genomes Project of four genome-wide association studies of lung cancer in populations of European ancestry (11,348 cases and 15,861 controls) and genotyped an additional 10,246 cases and 38,295 controls for follow-up. We identified large-effect genome-wide associations for squamous lung cancer with the rare variants BRCA2 p.Lys3326X (rs11571833, odds ratio (OR) = 2.47, P = 4.74 × 10(-20)) and CHEK2 p.Ile157Thr (rs17879961, OR = 0.38, P = 1.27 × 10(-13)). We also showed an association between common variation at 3q28 (TP63, rs13314271, OR = 1.13, P = 7.22 × 10(-10)) and lung adenocarcinoma that had been previously reported only in Asians. These findings provide further evidence for inherited genetic susceptibility to lung cancer and its biological basis. Additionally, our analysis demonstrates that imputation can identify rare disease-causing variants with substantive effects on cancer risk from preexisting genome-wide association study data.

Association between putative functional variants in the PSMB9 gene and risk of melanoma--re-analysis of published melanoma genome-wide association studies

To mine possibly hidden causal single-nucleotide polymorphisms (SNPs) of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome-wide association study (GWAS) data set with 1804 melanoma cases and 1026 cancer-free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in two other GWAS data sets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL data set (P = 0.013 and 0.004, respectively), but failed in validation using the Australian data set. Genotype-phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression level of PSMB9. Further experiments revealed that SNP rs2071480, which is in high LD with rs1351383 and rs2127675, may have a weak effect on the promoter activity of PSMB9. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility.

Association between functional polymorphisms in genes involved in the MAPK signaling pathways and cutaneous melanoma risk

Genome-wide association studies (GWASs) have mainly focused on top significant single nucleotide polymorphisms (SNPs), most of which did not have clear biological functions but were just surrogates for unknown causal variants. Studying SNPs with modest association and putative functions in biologically plausible pathways has become one complementary approach to GWASs. To unravel the key roles of mitogen-activated protein kinase (MAPK) pathways in cutaneous melanoma (CM) risk, we re-evaluated the associations between 47 818 SNPs in 280 MAPK genes and CM risk using our published GWAS dataset with 1804 CM cases and 1026 controls. We initially found 105 SNPs with P ≤ 0.001, more than expected by chance, 26 of which were predicted to be putatively functional SNPs. The risk associations with 16 SNPs around DUSP14 (rs1051849) and a previous reported melanoma locus MAFF/PLA2G6 (proxy SNP rs4608623) were replicated in the GenoMEL dataset (P < 0.01) but failed in the Australian dataset. Meta-analysis showed that rs1051849 in the 3' untranslated regions of DUSP14 was associated with a reduced risk of melanoma (odds ratio = 0.89, 95% confidence interval: 0.82-0.96, P = 0.003, false discovery rate = 0.056). Further genotype-phenotype correlation analysis using the 90 HapMap lymphoblastoid cell lines from Caucasians showed significant correlations between two SNPs (rs1051849 and rs4608623) and messenger RNA expression levels of DUSP14 and MAFF (P = 0.025 and P = 0.010, respectively). Gene-based tests also revealed significant SNPs were over-represented in MAFF, PLA2G6, DUSP14 and other 16 genes. Our results suggest that functional SNPs in MAPK pathways may contribute to CM risk. Further studies are warranted to validate our findings.

Whole-genome detection of disease-associated deletions or excess homozygosity in a case-control study of rheumatoid arthritis

Unlike genome-wide association studies, few comprehensive studies of copy number variation's contribution to complex human disease susceptibility have been performed. Copy number variations are abundant in humans and represent one of the least well-studied classes of genetic variants; in addition, known rheumatoid arthritis susceptibility loci explain only a portion of familial clustering. Therefore, we performed a genome-wide study of association between deletion or excess homozygosity and rheumatoid arthritis using high-density 550 K SNP genotype data from a genome-wide association study. We used a genome-wide statistical method that we recently developed to test each contiguous SNP locus between 868 cases and 1194 controls to detect excess homozygosity or deletion variants that influence susceptibility. Our method is designed to detect statistically significant evidence of deletions or homozygosity at individual SNPs for SNP-by-SNP analyses and to combine the information among neighboring SNPs for cluster analyses. In addition to successfully detecting the known deletion variants on major histocompatibility complex, we identified 4.3 and 28 kb clusters on chromosomes 10p and 13q, respectively, which were significant at a Bonferroni-type-corrected 0.05 nominal significant level. Independently, we performed analyses using PennCNV, an algorithm for identifying and cataloging copy numbers for individuals based on a hidden Markov model, and identified cases and controls that had chromosomal segments with copy number <2. Using Fisher's exact test for comparing the numbers of cases and controls with copy number <2 per SNP, we identified 26 significant SNPs (protective; more controls than cases) aggregating on chromosome 14 with P-values <10(-8).

Cell cycle-related genes as modifiers of age of onset of colorectal cancer in Lynch syndrome: a large-scale study in non-Hispanic white patients

Heterogeneity in age of onset of colorectal cancer in individuals with mutations in DNA mismatch repair genes (Lynch syndrome) suggests the influence of other lifestyle and genetic modifiers. We hypothesized that genes regulating the cell cycle influence the observed heterogeneity as cell cycle-related genes respond to DNA damage by arresting the cell cycle to provide time for repair and induce transcription of genes that facilitate repair. We examined the association of 1456 single nucleotide polymorphisms (SNPs) in 128 cell cycle-related genes and 31 DNA repair-related genes in 485 non-Hispanic white participants with Lynch syndrome to determine whether there are SNPs associated with age of onset of colorectal cancer. Genotyping was performed on an Illumina GoldenGate platform, and data were analyzed using Kaplan-Meier survival analysis, Cox regression analysis and classification and regression tree (CART) methods. Ten SNPs were independently significant in a multivariable Cox proportional hazards regression model after correcting for multiple comparisons (P < 5 × 10(-4)). Furthermore, risk modeling using CART analysis defined combinations of genotypes for these SNPs with which subjects could be classified into low-risk, moderate-risk and high-risk groups that had median ages of colorectal cancer onset of 63, 50 and 42 years, respectively. The age-associated risk of colorectal cancer in the high-risk group was more than four times the risk in the low-risk group (hazard ratio = 4.67, 95% CI = 3.16-6.92). The additional genetic markers identified may help in refining risk groups for more tailored screening and follow-up of non-Hispanic white patients with Lynch syndrome.

Model-free tests for genetic linkage

This unit covers statistical methods of linkage analysis that do not require the assumption of a detailed genetic model, as is required for standard lod score analysis. The unit has been updated to include the latest methods in sib-pair analysis, including updates to using the software program SIBPAL as well as the relative-pair analysis software applications GENEHUNTER, GENEHUNTER PLUS, and Merlin.

Abstract LB-335: Associations of functional polymorphisms in genes involved in the MAPK signaling pathways and cutaneous melanoma risk

Cutaneous melanoma (CM) is the most aggressive form of skin cancers, and its incidence is increasing annually in Caucasian populations. Previous linkage studies have identified two high-penetrance genes that influence CM risk, including cyclin-dependent kinase inhibitor 2A (CDKN2A) and cyclin-dependent kinase 4 (CDK4). Recent genome-wide association studies (GWASs) on CM have successfully identified several low-penetrance loci on CM susceptibility. However, most GWASs have mainly focused on top significant single nucleotide polymorphisms (SNPs) that could reach the required genome-wide significance level, most of which did not have clear biological functions but were just surrogates for causal variants. To overcome these limitations in GWASs, several complementary approaches have been proposed recently, such as pathway-based analysis and integration analysis of association results with gene expression. These applications have successfully revealed several new cancer susceptibility genes and pathways. Thus, studying SNPs with moderate significance and putative functions in a specific pathway might be helpful to identify SNPs with a relative small effect size and provide additional insight into genetic mechanisms of cancer. Because MAPK pathways play critical roles in development and progression of CM, in this study, we reevaluated the associations between SNPs in 280 MAPK genes and melanoma risk in the discovery GWAS dataset for 1804 CM cases and 1026 cancer-free controls. We found that 34 SNPs with a significance level of 0.001. We further performed in silico prediction using SNPinfo and identified 11 SNPs with putative functions or having a high LD with functional SNPs, five of which were selected to be validated in two other GWAS datasets. A meta-analysis showed that two SNPs (rs1051849 in DUSP14 and rs4608623 near MAFF/PLA2G6) were associated with reduced CM risk (P = 0.0002 for rs1051849 and P = 0.00007 for rs4608623). Further correlation analysis of gene expression levels and SNP genotypes from 90 HapMap CEU cell lines revealed that these two SNPs had significant effects on mRNA expression of DUSP14 and MAFF, respectively (P = 0.025 for rs1051849 and P = 0.010 for rs4608623). We also performed gene-based tests using the hypergeometric method and the Versatile Gene-Based Test for Genome-wide Association (VEGAS) and found that significant SNPs were over-represented in 9 MAPK genes including DUSP14 and MAFF. Our results indicated that genetic variants in DUSP14 and MAFF might be associated with melanoma risk by influencing mRNA expression levels. Further studies are warranted to validate our findings.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr LB-335. doi:1538-7445.AM2012-LB-335

Genome-wide association study identifies novel loci predisposing to cutaneous melanoma

We performed a multistage genome-wide association study of melanoma. In a discovery cohort of 1804 melanoma cases and 1026 controls, we identified loci at chromosomes 15q13.1 (HERC2/OCA2 region) and 16q24.3 (MC1R) regions that reached genome-wide significance within this study and also found strong evidence for genetic effects on susceptibility to melanoma from markers on chromosome 9p21.3 in the p16/ARF region and on chromosome 1q21.3 (ARNT/LASS2/ANXA9 region). The most significant single-nucleotide polymorphisms (SNPs) in the 15q13.1 locus (rs1129038 and rs12913832) lie within a genomic region that has profound effects on eye and skin color; notably, 50% of variability in eye color is associated with variation in the SNP rs12913832. Because eye and skin colors vary across European populations, we further evaluated the associations of the significant SNPs after carefully adjusting for European substructure. We also evaluated the top 10 most significant SNPs by using data from three other genome-wide scans. Additional in silico data provided replication of the findings from the most significant region on chromosome 1q21.3 rs7412746 (P = 6 × 10(-10)). Together, these data identified several candidate genes for additional studies to identify causal variants predisposing to increased risk for developing melanoma.

Novel genetic variants in the chromosome 5p15.33 region associate with lung cancer risk

Chromosome 5p15.33 has been identified by genome-wide association studies as one of the regions that associate with lung cancer risk. A few single-nucleotide polymorphisms (SNPs) in the telomerase reverse transcriptase (TERT) and cleft lip and palate transmembrane 1-like (CLPTM1L) genes located in this region have shown consistent associations. We performed dense genotyping of SNPs in this region to refine the previously reported association signals for lung cancer risk. Two hundred and fifteen SNPs were genotyped on an Illumina iSelect panel, in a hospital-based case-control study of 1681 lung cancer cases and 1235 unaffected controls. Association was tested using unconditional logistic regression, while adjusting for age, sex and pack-years smoked. Furthermore, since many of the SNPs were in linkage disequilibrium (LD), haplotype blocks were constructed, from which tagging SNPs at an r(2) threshold of ≥0.95 were included in a stepwise forward selection logistic regression model. Of the 215 SNPs, 69 were significant at P < 0.05 in univariate analysis; of these, 35 SNPs meeting the r(2) threshold were included in the multiple logistic regression model. Two SNPs, rs370348 (odds ratio = 0.76, P = 1.6 × 10(-6)) and rs4975538 (odds ratio = 1.18, P = 0.005), significantly associated with risk in the overall sample. Among ever smokers, rs4975615 (odds ratio = 0.75, P = 1.2 × 10(-4)) and rs4975538 (odds ratio = 1.26, P = 0.002) were significant, whereas among never-smokers, rs451360 (odds ratio = 0.62, P = 7.6 × 10(-5)) was significant. We refined the consistent association signal in this region, allowing for the considerable LD between SNPs and identified four novel SNPs that were independently and significantly associated with lung cancer risk. Results of these analyses strongly suggest effects on risk from several loci in the TERT/CLPTM1L region.

An analysis of single nucleotide polymorphisms of 125 DNA repair genes in the Texas genome-wide association study of lung cancer with a replication for the XRCC4 SNPs

DNA repair genes are important for maintaining genomic stability and limiting carcinogenesis. We analyzed all single nucleotide polymorphisms (SNPs) of 125 DNA repair genes covered by the Illumina HumanHap300 (v1.1) BeadChips in a previously conducted genome-wide association study (GWAS) of 1154 lung cancer cases and 1137 controls and replicated the top-hits of XRCC4 SNPs in an independent set of 597 cases and 611 controls in Texas populations. We found that six of 20 XRCC4 SNPs were associated with a decreased risk of lung cancer with a P-value of 0.01 or lower in the discovery dataset, of which the most significant SNP was rs10040363 (P for allelic test=4.89 x 10⁻⁴). Moreover, the data in this region allowed us to impute a potentially functional SNP rs2075685 (imputed P for allelic test=1.3 x 10⁻³). A luciferase reporter assay demonstrated that the rs2075685G>T change in the XRCC4 promoter increased expression of the gene. In the replication study of rs10040363, rs1478486, rs9293329, and rs2075685, however, only rs10040363 achieved a borderline association with a decreased risk of lung cancer in a dominant model (adjusted OR=0.80, 95% CI=0.62-1.03 and P=0.079). In the final combined analysis of both the Texas GWAS discovery and replication datasets, the strength of the association was increased for rs10040363 (adjusted OR=0.77, 95% CI=0.66-0.89, P(dominant)=5 x 10⁻⁴ and P for trend=5 x 10⁻⁴) and rs1478486 (adjusted OR=0.82, 95% CI=0.71-0.94, P(dominant)=6 x 10⁻³ and P for trend=3.5 x 10⁻³). Finally, we conducted a meta-analysis of these XRCC4 SNPs with available data from published GWA studies of lung cancer with a total of 12,312 cases and 47,921 controls, in which none of these XRCC4 SNPs was associated with lung cancer risk. It appeared that rs2075685, although associated with increased expression of a reporter gene and lung cancer risk in the Texas populations, did not have an effect on lung cancer risk in other populations. This study underscores the importance of replication using published data in larger populations.

Association of a novel functional promoter variant (rs2075533 C>T) in the apoptosis gene TNFSF8 with risk of lung cancer--a finding from Texas lung cancer genome-wide association study

Published genome-wide association studies (GWASs) have identified few variants in the known biological pathways involved in lung cancer etiology. To mine the possibly hidden causal single nucleotide polymorphisms (SNPs), we explored all SNPs in the extrinsic apoptosis pathway from our published GWAS dataset for 1154 lung cancer cases and 1137 cancer-free controls. In an initial association analysis of 611 tagSNPs in 41 apoptosis-related genes, we identified only 10 tagSNPs associated with lung cancer risk with a P value<10(-2), including four tagSNPs in DAPK1 and three tagSNPs in TNFSF8. Unlike DAPK1 SNPs, TNFSF8 rs2181033 tagged other four predicted functional but untyped SNPs (rs776576, rs776577, rs31813148 and rs2075533) in the promoter region. Therefore, we further tested binding affinity of these four SNPs by performing the electrophoretic mobility shift assay. We found that only rs2075533T allele modified levels of nuclear proteins bound to DNA, leading to significantly decreased expression of luciferase reporter constructs by 5- to -10-fold in H1299, HeLa and HCT116 cell lines compared with the C allele. We also performed a replication study of the untyped rs2075533 in an independent Texas population but did not confirm the protective effect. We further performed a mini meta-analysis for SNPs of TNFSF8 obtained from other four published lung cancer GWASs with 12  214 cases and 47  721 controls, and we found that only rs3181366 (r2=0.69 with the untyped rs2075533) was associated to lung cancer risk (P=0.008). Our findings suggest a possible role of novel TNFSF8 variants in susceptibility to lung cancer.

Genome-wide association study in alopecia areata implicates both innate and adaptive immunity

Alopecia areata (AA) is among the most highly prevalent human autoimmune diseases, leading to disfiguring hair loss due to the collapse of immune privilege of the hair follicle and subsequent autoimmune attack. The genetic basis of AA is largely unknown. We undertook a genome-wide association study (GWAS) in a sample of 1,054 cases and 3,278 controls and identified 139 single nucleotide polymorphisms that are significantly associated with AA (P <or= 5 x 10(-7)). Here we show an association with genomic regions containing several genes controlling the activation and proliferation of regulatory T cells (T(reg) cells), cytotoxic T lymphocyte-associated antigen 4 (CTLA4), interleukin (IL)-2/IL-21, IL-2 receptor A (IL-2RA; CD25) and Eos (also known as Ikaros family zinc finger 4; IKZF4), as well as the human leukocyte antigen (HLA) region. We also find association evidence for regions containing genes expressed in the hair follicle itself (PRDX5 and STX17). A region of strong association resides within the ULBP (cytomegalovirus UL16-binding protein) gene cluster on chromosome 6q25.1, encoding activating ligands of the natural killer cell receptor NKG2D that have not previously been implicated in an autoimmune disease. By probing the role of ULBP3 in disease pathogenesis, we also show that its expression in lesional scalp from patients with AA is markedly upregulated in the hair follicle dermal sheath during active disease. This study provides evidence for the involvement of both innate and acquired immunity in the pathogenesis of AA. We have defined the genetic underpinnings of AA, placing it within the context of shared pathways among autoimmune diseases, and implicating a novel disease mechanism, the upregulation of ULBP ligands, in triggering autoimmunity.

Detection of disease-associated deletions in case-control studies using SNP genotypes with application to rheumatoid arthritis

Genomic deletions have long been known to play a causative role in microdeletion syndromes. Recent whole-genome genetic studies have shown that deletions can increase the risk for several psychiatric disorders, suggesting that genomic deletions play an important role in the genetic basis of complex traits. However, the association between genomic deletions and common, complex diseases has not yet been systematically investigated in gene mapping studies. Likelihood-based statistical methods for identifying disease-associated deletions have recently been developed for familial studies of parent-offspring trios. The purpose of this study is to develop statistical approaches for detecting genomic deletions associated with complex disease in case-control studies. Our methods are designed to be used with dense single nucleotide polymorphism (SNP) genotypes to detect deletions in large-scale or whole-genome genetic studies. As more and more SNP genotype data for genome-wide association studies become available, development of sophisticated statistical approaches will be needed that use these data. Our proposed statistical methods are designed to be used in SNP-by-SNP analyses and in cluster analyses based on combined evidence from multiple SNPs. We found that these methods are useful for detecting disease-associated deletions and are robust in the presence of linkage disequilibrium using simulated SNP data sets. Furthermore, we applied the proposed statistical methods to SNP genotype data of chromosome 6p for 868 rheumatoid arthritis patients and 1,197 controls from the North American Rheumatoid Arthritis Consortium. We detected disease-associated deletions within the region of human leukocyte antigen in which genomic deletions were previously discovered in rheumatoid arthritis patients.

Comparison of genome-wide single-nucleotide polymorphism linkage analyses in Caucasian and Hispanic NARAC families

We performed linkage analysis on families with rheumatoid arthritis, stratifying by ethnic origin. We compared results using either Kong and Cox nonparametric LOD scores or MOD score analysis using the software GeneHunter MODSCORE. We first applied SNPLINK to remove markers showing excess linkage disequilibrium from the SNPs in the Illumina IV SNP Linkage panel. In this analysis there were 659 self-reported Caucasian families and 29 self-reported Hispanic families in the NARAC collection. Chromosome 19 yielded MOD scores > 3.00 in the Hispanic group, while chromosomes 2, 6, 7, 11, and XY had MOD scores > 3.00 in the Caucasian group. We performed simulation studies to evaluate the empirical distribution of the MOD score for autosomal loci separately in Hispanics and Caucasians. Results showed genome-wide significant evidence for linkage in Caucasians for chromosomes 2q and 6p, but no significant evidence for any linkages in the Hispanics, including little evidence for linkage to chromosome 6p in this group. An examination of the difference of phenotypes in two ethnic groups suggested significantly earlier mean age of onset, higher percentage of anti-cyclic citrullinated peptide positive people, and lower percentage of affected people carrying shared epitopes in Hispanics than those in Caucasians. A larger sample size of the Hispanic group is needed to identify linkage regions.

Linkage analyses of rheumatoid arthritis and related quantitative phenotypes: the GAW15 experience

The group that formed on the theme of linkage analyses of rheumatoid arthritis RA and related phenotypes (Group 10) in the Genetic Analysis Workshop 15 comprised 18 sets of investigators. Two data sets were available: one was a real set provided by the North American Rheumatoid Arthritis Consortium and collaborators in Canada, France (European Consortium Of Rheumatoid Arthritis Families) and the UK; the other was a simulated data set modelled after the real data set. Whereas a majority of the investigators analyzed the RA affection status as a binary phenotype, a few contributions considered data on correlated quantitative traits such as anti-cyclic citrullinated peptide and rheumatoid factor-immunoglobulin M. The different investigators applied a wide spectrum of linkage methods. As expected, most methods could identify the human leukocyfeantigen region on chromosome 6 as a major genetic factor for RA. In addition, some novel chromosomal regions provided significant evidence of linkage in multiple contributions in the group. In this report, we discuss the different strategies explored by the different investigators with the common goal of improving the power to detect linkage.

TRAF1-C5 as a risk locus for rheumatoid arthritis--a genomewide study

BACKGROUND

Rheumatoid arthritis has a complex mode of inheritance. Although HLA-DRB1 and PTPN22 are well-established susceptibility loci, other genes that confer a modest level of risk have been identified recently. We carried out a genomewide association analysis to identify additional genetic loci associated with an increased risk of rheumatoid arthritis.

METHODS

We genotyped 317,503 single-nucleotide polymorphisms (SNPs) in a combined case-control study of 1522 case subjects with rheumatoid arthritis and 1850 matched control subjects. The patients were seropositive for autoantibodies against cyclic citrullinated peptide (CCP). We obtained samples from two data sets, the North American Rheumatoid Arthritis Consortium (NARAC) and the Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA). Results from NARAC and EIRA for 297,086 SNPs that passed quality-control filters were combined with the use of Cochran-Mantel-Haenszel stratified analysis. SNPs showing a significant association with disease (P<1x10(-8)) were genotyped in an independent set of case subjects with anti-CCP-positive rheumatoid arthritis (485 from NARAC and 512 from EIRA) and in control subjects (1282 from NARAC and 495 from EIRA).

RESULTS

We observed associations between disease and variants in the major-histocompatibility-complex locus, in PTPN22, and in a SNP (rs3761847) on chromosome 9 for all samples tested, the latter with an odds ratio of 1.32 (95% confidence interval, 1.23 to 1.42; P=4x10(-14)). The SNP is in linkage disequilibrium with two genes relevant to chronic inflammation: TRAF1 (encoding tumor necrosis factor receptor-associated factor 1) and C5 (encoding complement component 5).

CONCLUSIONS

A common genetic variant at the TRAF1-C5 locus on chromosome 9 is associated with an increased risk of anti-CCP-positive rheumatoid arthritis.

STAT4 and the risk of rheumatoid arthritis and systemic lupus erythematosus

BACKGROUND

Rheumatoid arthritis is a chronic inflammatory disease with a substantial genetic component. Susceptibility to disease has been linked with a region on chromosome 2q.

METHODS

We tested single-nucleotide polymorphisms (SNPs) in and around 13 candidate genes within the previously linked chromosome 2q region for association with rheumatoid arthritis. We then performed fine mapping of the STAT1-STAT4 region in a total of 1620 case patients with established rheumatoid arthritis and 2635 controls, all from North America. Implicated SNPs were further tested in an independent case-control series of 1529 patients with early rheumatoid arthritis and 881 controls, all from Sweden, and in a total of 1039 case patients and 1248 controls from three series of patients with systemic lupus erythematosus.

RESULTS

A SNP haplotype in the third intron of STAT4 was associated with susceptibility to both rheumatoid arthritis and systemic lupus erythematosus. The minor alleles of the haplotype-defining SNPs were present in 27% of chromosomes of patients with established rheumatoid arthritis, as compared with 22% of those of controls (for the SNP rs7574865, P=2.81x10(-7); odds ratio for having the risk allele in chromosomes of patients vs. those of controls, 1.32). The association was replicated in Swedish patients with recent-onset rheumatoid arthritis (P=0.02) and matched controls. The haplotype marked by rs7574865 was strongly associated with lupus, being present on 31% of chromosomes of case patients and 22% of those of controls (P=1.87x10(-9); odds ratio for having the risk allele in chromosomes of patients vs. those of controls, 1.55). Homozygosity of the risk allele, as compared with absence of the allele, was associated with a more than doubled risk for lupus and a 60% increased risk for rheumatoid arthritis.

CONCLUSIONS

A haplotype of STAT4 is associated with increased risk for both rheumatoid arthritis and systemic lupus erythematosus, suggesting a shared pathway for these illnesses.

Dissecting the heterogeneity of rheumatoid arthritis through linkage analysis of quantitative traits

OBJECTIVE

To dissect the heterogeneity of rheumatoid arthritis (RA) through linkage analysis of quantitative traits, specifically, IgM rheumatoid factor (IgM-RF) and anti-cyclic citrullinated peptide (anti-CCP) autoantibody titers.

METHODS

Subjects, 1,002 RA patients from 491 multiplex families recruited by the North American RA Consortium, were typed for 379 microsatellite markers. Anti-CCP titers were determined based on a second-generation enzyme-linked immunosorbent assay, and IgM-RF levels were quantified by immunonephelometry. We used the Merlin statistical package to perform nonparametric quantitative trait linkage analysis.

RESULTS

For each of the quantitative traits, evidence of linkage, with logarithm of odds (LOD) scores of >1.0, was found in 9 regions. For both traits, the strongest evidence of linkage was for marker D6S1629 on chromosome 6p (LOD 14.02 for anti-CCP and LOD 12.09 for RF). Six other regions with LOD scores of >1.0 overlapped between the 2 traits, on chromosomes 1p21.1, 5q15, 8p23.1, 16p12.1, 16q23.1, and 18q21.31. Evidence of linkage to anti-CCP titer but not to RF titer was found in 2 regions (chromosomes 9p21.3 and 10q21.1), and evidence of linkage to RF titer but not to anti-CCP titer was found in 2 regions (chromosomes 5p15.2 and 1q42.3). Several covariates were significantly associated with 1 or both traits, and linkage analysis exploring the covariate effects revealed striking effects of sex in modulating linkage signals for several chromosomal regions. For example, sex had a striking impact on the linkage results for both quantitative traits on chromosome 6p (P = 0.0007 for anti-CCP titer and P = 0.0012 for RF titer), suggesting a sex-HLA region interaction.

CONCLUSION

Analysis of quantitative components of RA is a promising approach for dissecting the genetic heterogeneity of this complex disorder. These results highlight the potential importance of sex or other covariates that may modulate some of the genetic effects that influence the risk of specific disease manifestations.

High-density SNP analysis of 642 Caucasian families with rheumatoid arthritis identifies two new linkage regions on 11p12 and 2q33

We have completed a genome wide linkage scan using >5700 informative single-nucleotide polymorphism (SNP) markers (Illumina IV SNP linkage panel) in 642 Caucasian families containing affected sibling pairs with rheumatoid arthritis (RA), ascertained by the North American Rheumatoid Arthritis Consortium. The results show striking new evidence of linkage at chromosomes 2q33 and 11p12 with logarithm of odds (LOD) scores of 3.52 and 3.09, respectively. In addition to a strong and broad linkage interval surrounding the major histocompatibility complex (LOD>16), regions with LOD>2.5 were observed on chromosomes 5 and 10. Additional linkage evidence (LOD scores between 1.46 and 2.35) was also observed on chromosomes 4, 7, 12, 16 and 18. This new evidence for multiple regions of genetic linkage is partly explained by the significantly increased information content of the Illumina IV SNP linkage panel (75.6%) compared with a standard microsatellite linkage panel utilized previously (mean 52.6%). Stratified analyses according to whether or not the sibling pair members showed elevated anticyclic citrullinated peptide titers indicates significant variation in evidence for linkage among strata on chromosomes 4, 5, 6 and 7. Overall, these new linkage data should reinvigorate efforts to utilize positional information to identify susceptibility genes for RA.

Genome-wide linkage scan for spontaneous DZ twinning

In humans, spontaneous DZ twinning is known to have a genetic basis. A prior investigation in the Flemish and Dutch population showed that the phenotype of 'having DZ twins' was consistent with an autosomal monogenic dominant model, with a gene frequency of 3.5% and a female-specific lifetime penetrance of 10%. Recessive, X-linked, polygenic and sporadic models were rejected. This study reports on a genome-wide scan of 14 Flemish families containing 57 mothers of spontaneous DZ twins. Two-point linkage analysis using the autosomal dominant model showed nine chromosomal regions with a LOD score around 1. After multipoint linkage analysis, including heterogeneity, three chromosomes continued to give high LOD scores. These regions were further haplotyped with additional markers at 1 cM distance. The multipoint analysis was not in favour of linkage of the DZ twinning trait in most candidate genes and other regions (LOD score < -2) under the genetic model of autosomal dominance. To further evaluate the evidence for linkage given some uncertainty about the correct mode of inheritance of twinning susceptibility other models of inheritance were tested. Results of this analysis showed all models gave highest LOD scores under dominant models. If heterogeneity among the families is taken into account, the peaks that were observed on chromosome 2, 7 and 18 could well contain a potential gene contributing to DZ twinning. These results give suggestive evidence that the mode of inheritance of DZ twinning is probably more complex than was originally expected.

Genome-wide meta-analysis for rheumatoid arthritis

Meta-analysis is being increasingly used as a tool for integrating data from different studies of complex phenotypes, because the power of any one study to identify causal loci is limited. We applied a novel meta-analytical approach (Loesgen et al. in Genet Epidemiol 21(Suppl 1):S142-S147, 2001) in compiling results from four studies of rheumatoid arthritis in Caucasians including two studies from NARAC (Jawaheer et al. in Am J Hum Genet 68:927-936, 2001; Jawaheer et al. in Arthritis Rheum 48:906-916, 2003), one study from the UK (MacKay et al. in Arthritis Rheum 46:632-639, 2001) and one from France (Cornelis et al. in Proc Natl Acad Sci USA 95:10746-10750, 1998). For each study, we obtained NPL scores by performing interval mapping (2 cM intervals) using GeneHunter2 (Kruglyak et al. in Am J Hum Genet 58:1347-1363, 1996; Markianos et al. in Am J Hum Genet 68:963-977, 2001). The marker maps differed among the three consortium groups, therefore, the marker maps were aligned after the interval mapping was completed and the NPL scores that were within 1 cM of each other were combined using the method of Loesgen et al. (Genet Epidemiol 21(Suppl 1):S142-S147, 2001) by calculating the weighted average of the NPL score. This approach avoids some problems in analysis encountered by using GeneHunter2 when some markers in the sample are not genotyped. This procedure provided marginal evidence (P<0.05) of linkage on chromosome 1, 2, 5 and 18, strong evidence (P<0.01) on chromosomes 8 and 16, and overwhelming evidence in the HLA region of chromosome 6.

Screening the genome for rheumatoid arthritis susceptibility genes: a replication study and combined analysis of 512 multicase families

OBJECTIVE

A number of non-HLA loci that have shown evidence (P < 0.05) for linkage with rheumatoid arthritis (RA) have been previously identified. The present study attempts to confirm these findings.

METHODS

We performed a second genome-wide screen of 256 new multicase RA families recruited from across the United States by the North American Rheumatoid Arthritis Consortium. Affected sibling pair analysis on the new data set was performed using SIBPAL. We subsequently combined our first and second data sets in an attempt to enhance the evidence for linkages in a larger sample size. We also evaluated the impact of covariates on the support for linkage, using LODPAL.

RESULTS

Evidence of linkage at 1p13 (D1S1631), 6p21.3 (the HLA complex), and 18q21 (D18S858) (P < 0.05) was replicated in this independent data set. In addition, there was new evidence for linkage at 9p22 (D9S1121 [P = 0.001]) and 10q21 (D10S1221 [P = 0.0002] and D10S1225 [P = 0.0038]) in the current data set. The combined analysis of both data sets (512 families) showed evidence for linkage at the level of P < 0.005 at 1p13 (D1S1631), 1q43 (D1S235), 6q21 (D6S2410), 10q21 (D10S1221), 12q12 (D12S398), 17p13 (D17S1298), and 18q21 (D18S858). Linkage at HLA was also confirmed (P < 5 x 10(-12)). Inclusion of DRB1*04 as a covariate significantly increased the probability of linkage on chromosome 6. In addition, some linkages on chromosome 1 showed improved significance when modeling DRB1*04 or rheumatoid factor positivity as covariates.

CONCLUSION

These results provide a rational basis for pursuing high-density linkage and association studies of RA in several regions outside of the HLA region, particularly on chromosomes 1p, 1q, and 18q.

The mouse beta B1-crystallin promoter: strict regulation of lens fiber cell specificity

Previous studies have shown that the chicken beta B1-crystallin promoter (-434/+30) contains all of the signals necessary to specifically direct high level expression of heterologous genes to the lens fiber cells of mice. In the present study, the mouse beta B1-crystallin gene was cloned, and its regulation was investigated to further elucidate the mechanisms controlling lens fiber cell-specific gene expression. Phylogenetic footprinting analysis of the 5' flanking sequence from the mouse, rat, human and chicken beta B1-crystallin genes identified several known and putative functional cis elements including the PL2 element which is required for lens-specific expression of the chicken beta B1 promoter. Surprisingly, however, all six mouse beta B1-crystallin/CAT constructs tested (-1493/+44, -1493/+30, -870/+30, -250/+30, -135/+30 and -98/+30) were inactive in three different mammalian lens-derived cell lines while only the -870/+30 and -98/+30 constructs were active in chicken primary patched lens epithelial cells. In contrast, the chicken beta B1-crystallin promoter (-434/+30) was transcriptionally active in all lens-derived cells tested. Transgenic mice harboring a mouse beta B1-crystallin -1493/+44 CAT construct did express the transgene specifically in lens fiber cells, however, at lower levels than that previously reported for a chicken -434/+30 CAT construct. These data suggest that, as in other crystallin genes, the regulatory signals controlling lens fiber cell-specific expression are conserved between chicken and mouse. However, the inability of the mouse beta B1-crystallin promoter to function in mammalian lens-derived cultured cells implies that this gene has acquired additional cis-regulatory elements to ensure lens fiber cell specificity.

A genomewide screen in multiplex rheumatoid arthritis families suggests genetic overlap with other autoimmune diseases

Rheumatoid arthritis (RA) is an autoimmune/inflammatory disorder with a complex genetic component. We report the first major genomewide screen of multiplex families with RA gathered in the United States. The North American Rheumatoid Arthritis Consortium, using well-defined clinical criteria, has collected 257 families containing 301 affected sibling pairs with RA. A genome screen for allele sharing was performed, using 379 microsatellite markers. A nonparametric analysis using SIBPAL confirmed linkage of the HLA locus to RA (P < .00005), with lambdaHLA = 1.79. However, the analysis also revealed a number of non-HLA loci on chromosomes 1 (D1S235), 4 (D4S1647), 12 (D12S373), 16 (D16S403), and 17 (D17S1301), with evidence for linkage at a significance level of P<.005. Analysis of X-linked markers using the MLOD method from ASPEX also suggests linkage to the telomeric marker DXS6807. Stratifying the families into white or seropositive subgroups revealed some additional markers that showed improvement in significance over the full data set. Several of the regions that showed evidence for nominal significance (P < .05) in our data set had previously been implicated in RA (D16S516 and D17S1301) or in other diseases of an autoimmune nature, including systemic lupus erythematosus (D1S235), inflammatory bowel disease (D4S1647, D5S1462, and D16S516), multiple sclerosis (D12S1052), and ankylosing spondylitis (D16S516). Therefore, genes in the HLA complex play a major role in RA susceptibility, but several other regions also contribute significantly to overall genetic risk.