Mapping a gene for rheumatoid arthritis on chromosome 18q21

William Tapper; Andrew Collins; Newton E Morton

doi:10.1186/1753-6561-1-S1-S18

Volume 1 Supplement 1

Genetic Analysis Workshop 15: Gene Expression Analysis and Approaches to Detecting Multiple Functional Loci

Proceedings
Open Access

Mapping a gene for rheumatoid arthritis on chromosome 18q21

William Tapper¹Email author,
Andrew Collins¹ and
Newton E Morton¹

BMC Proceedings20071(Suppl 1):S18

https://doi.org/10.1186/1753-6561-1-S1-S18

Published: 18 December 2007

Abstract

Although single chi-square analysis of the North American Rheumatoid Arthritis Consortium (NARAC) data identifies many single-nucleotide polymorphisms (SNPs) with p-values less than 0.05, none remain significant after Bonferroni correction. In contrast, CHROMSCAN evades heavy Bonferroni correction and auto-correlation between SNPs by using composite likelihood to model association across all markers in a region and permutation to assess significance. Analysis by CHROMSCAN identifies a 36-kb interval that includes the most significant SNP (msSNP) observed in a 10-Mb target suggested by linkage. Unexpectedly, stratification by gender and age of onset shows that association evidence comes almost entirely from females with age of onset less than 40. Combining evidence from a meta-analysis of linkage studies and three subsets of the NARAC data provides significant evidence for a determinant of rheumatoid arthritis in a 36-kb interval and illustrates the principle that estimates of location and its information are more powerful than estimates of p-values alone.

Keywords

Association Mapping
Causal Locus
Exponential Decline
Complex Inheritance
North American Rheumatoid Arthritis Consortium

Background

Initially, linkage mapping dealt with rare and highly penetrant genes. Without cytogenetic assignment, the preferred strategy was segregation analysis to determine all relevant parameters except recombination, followed by linkage analysis to determine recombination frequency [1]. Complex inheritance with uncertain segregation parameters proved much more difficult, giving rise to many unconfirmed claims based on microsatellites and leading to meta-analysis without point locations [2]. The HapMap project provides dense SNPs that can be used to localize causal loci with or without pedigrees. This procedure, called association mapping, revolutionized identification of disease genes. Recent developments of linkage disequilibrium units (LDU), composite likelihood, control of auto-correlation, and meta-analysis are incorporated into the CHROMSCAN program [3, 4] to increase its precision for association mapping. Here we use these methods to establish the location and weight of evidence for a gene predisposing to rheumatoid arthritis.

Methods

Data preparation

The data, provided by NARAC (North American Rheumatoid Arthritis Consortium) consist of 2300 single-nucleotide polymorphisms (SNPs) in a 10-Mb region of 18q21 with linkage evidence in U.S. and French scans [5]. Illumina genotyped these markers in 460 cases and 460 controls, matched for age and gender, from New York. The genotypic data for controls were screened and 7 SNPs with $χ_{1}^{2}$ > 10 for the Hardy-Weinberg test [6] were removed, leaving 2293 to be analyzed. CHROMSCAN requires SNPs to be located on both physical and LDU scales. Physical locations were taken from build 35 of the human genome sequence. Unlike physical maps, study-specific and various LDU maps are available, corresponding to the four HapMap samples separately and combined (CEU, CHB, JPT, YRI, and cosmopolitan). The LDU map with the highest SNP density and population attributes closest to the experimental data should be optimal. We therefore used LDU locations relative to the CEU HapMap data with a density of 1 SNP per 863 bp compared to 1 SNP per 4139 bp in the NARAC data. We also used the kilobase map to determine the robustness and power of LDU maps compared with physical maps.

LDU map construction

The theory for constructing LDU maps has been described [7]. Briefly, the LDU distance for the i^th SNP interval is given by ε_id_i, where ε_idescribes the exponential decline of association with physical distance d_iin kb. Values of ε_iare estimated by composite likelihood that fits the Malecot model [8] to multiple pairwise diplotype data. The Malecot equation, given by $ρ = (1 - L) M e^{- \sum ε_{i} d_{i}} + L$ , uses additional parameters to describes association at the last major bottleneck (M), and residual association at large distance (L) to predict rho (ρ), the probability of association.

Association mapping

The CHROMSCAN program [3] uses a model similar to LDU maps except the exponential term is replaced by εΔ(S_i- S) to estimate the location (S) of a disease gene, where S_iis the location of the i^thmarker in kilobases or LDU. The Kronecker Δ is used for map direction and assures a correct sign, with Δ = 1 if S_i≥ S or -1 if S_i<S. To calculate the expected association with distance, z_i, the model becomes $z_{i} = (1 - L) M e^{- ε Δ (S_{i} - S)} + L$ , where M is diminished by complex inheritance and L is the association at large distance. The observed association ${\overset{⌢}{z}}_{i}$ is determined by a 2 × 2 table between affection status and the two alleles of each SNP to give $| \begin{matrix} a & b \\ c & d \end{matrix} |$ and ${\overset{⌢}{z}}_{i} = \frac{(a d - b c)}{(a + b) (b + d)}$ , where ad - bc ≥ 0 and b ≤ c is ensured by rearrangement of columns and rows [9]. Given the observed associations ${\hat{z}}_{i}$ , the Malecot parameters are estimated iteratively using composite likelihood, which evades a heavy Bonferroni correction by combining information over all loci within a region as $Λ = {\sum K_{i} ({\hat{z}}_{i} - z_{i})}^{2}$ , where ${\hat{z}}_{i}$ and z_iare the observed and expected association values, respectively, at the i^th SNP. Their squared difference is weighted by information (K_i) which is estimated as: $K_{i} = χ_{1}^{2} / {\hat{z}}_{i}^{2}$ , where $χ_{1}^{2}$ is the Pearson $χ_{1}^{2}$ from the 2 × 2 table.

Sub-hypotheses of the Malecot model are used to test for a causal polymorphism. Model A, which estimates none of the parameters and uses M = 0 with predicted L [10], is taken as the null hypothesis H₀ in which there is no association between affection status and SNPs. Model D estimates M, S, and L. Therefore the Λ_A - Λ_D comparison tests for a disease determinant at location S. For both models, ε is fixed to 1 for the LDU map and to a value of ε determined from pairwise marker-by-marker association data for the kilobase map. In order to account for autocorrelation between SNPs as a result of LD, the significance of evidence is determined by a rank-based permutation test [3].

Three separate analyses of the data were performed by CHROMSCAN. The first is a preliminary screen of the entire 10-Mb bin, which is divided into 18 nonoverlapping regions, each with at least 30 SNPs and covering at least 10 LDUs. To determine accurate levels of significance, the number of permutation replicates must approach the actual level of significance so that interpolation of the variance under H₁ is reliable. To minimize computation time, the initial analysis was restricted to 100 replicates. Significant regions identified by the initial screen were re-analyzed separately using 1000 and 5000 replicates in order to verify convergence. To demonstrate the power of LDU maps, this analysis was repeated using the kilobase map and two estimates of the exponential decline ε derived from the significant region and the 10-Mb region [11]. The risk for rheumatoid arthritis is elevated in females, especially with late onset (≥35–≤60) [12]. Our third analysis therefore stratified cases into three groups corresponding to males, females with onset ≤39, and females with onset ≥40. The partition of females around an onset age of 40 was chosen to give approximately equal numbers of 'early' and 'late' onset cases. Unaffected controls for the three groups were all males (with similar age and total number of individuals as affected males), and females divided by current age to give similar total numbers of individuals as cases, respectively. This analysis was restricted to significant regions from the initial screen and used 5000 replicates.

Results

Association mapping

Single chi-square analyses of the 10-Mb region identifies 125 SNPs with p < 0.05, none of which reach significance after Bonferroni correction (0.05/2293). The initial screen by CHROMSCAN divides the 18q21 bin into 18 nonoverlapping regions. Although the most significant SNP (msSNP, rs3745064) occurs in region 6, the next msSNP in region 11 is deceptively close in terms of significance, and several other regions contain suggestive SNPs (Table 1). In contrast, the composite likelihood approach, which models association across all markers in a region, identifies region 6 as the only significant region (p = 0.01259). The intensive screen of region 6 identified a large increase in significance between 100 and 1000 replicates, which is attributed to the relationship between number of replicates and significance, while the small decrease in significance between 1000 and 5000 replicates suggests that convergence has been achieved (Table 2). These analyses estimate a causal locus (S) at 53308 kb.

Table 1

Regions screened with 100 replicates

Region	No. SNPs	First kb	Last kb	$χ_{1}^{2}$ for msSNP	Composite likelihood
					$χ_{1}^{2}$	LOD₁	p-Values
1	256	48896	49928	6.58	0.05948	0.01291	0.80701
2	238	49935	51009	11.18	0.00691	0.00150	0.93358
3	197	51010	51919	10.83	2.22538	0.48323	0.13584
4	132	51926	52398	4.42	0.05000	0.01086	0.82274
5	241	52398	53262	5.90	0.00011	0.00002	0.99158
6	113	53262	53596	12.25	6.22830	1.35246	0.01259
7	74	53599	53919	6.84	2.95054	0.64070	0.08592
8	38	53920	54006	1.75	0.78229	0.16987	0.37637
9	103	54008	54448	6.64	1.71241	0.37185	0.19074
10	49	54496	54974	5.23	0.50780	0.11027	0.47592
11	70	54977	55285	12.06	0.34910	0.07581	0.55438
12	60	55286	55381	10.89	0.00011	0.00002	0.99147
13	89	55382	55634	4.24	0.13215	0.02870	0.71588
14	29	55641	55838	2.42	0.05010	0.01088	0.82258
15	134	55839	56399	5.66	0.97584	0.21190	0.32321
16	209	56408	57248	3.87	0.00166	0.00036	0.96744
17	55	57253	57545	9.51	1.51557	0.32910	0.21834
18	200	57548	58415	3.86	1.58390	0.34394	0.20826

Table 2

Intensive screening of region 6

No. replicates	Map^a	ε	S		K	SE	$χ_{1}^{2}$	LOD₁	p
			kb	LDU
100	LDU	1.000	53308.370	1195.126	41.233	0.1557	6.2283	1.3525	0.0126
1000	LDU	1.000	53308.370	1195.126	44.329	0.1502	9.6876	2.1036	0.0019
5000	LDU	1.000	53308.370	1195.126	36.611	0.1589	8.3245	1.8076	0.0039
1000	kb	0.021	53306.868	1195.126	0.0195	7.1529	7.2178	1.5673	0.0072
1000	kb	0.031	53307.748	1195.126	0.0367	5.2220	7.2910	1.5832	0.0069

^aValues of information (K) and standard error (SE) are relative to S_LDU and S_kb, indicated by the map column, when using the LDU and kb map, respectively.

The CHROMSCAN analysis of region 6 was repeated using the kilobase map so that its performance can be compared with the LDU map. Using a kilobase map requires specification of the exponential decline ε [11]. Two values of ε, corresponding to the 10 Mb interval (0.021) or region 6 alone (0.031), were investigated. Despite the large difference between ε values for the kilobase map, the significance level and location were almost identical. However, the ratios of $χ_{1}^{2}$ indicate that the kilobase maps have a relative efficiency of 75% compared with an LDU map at 1000 replicates (Table 2).

Because King et al. [12] demonstrated that the risk for rheumatoid arthritis is elevated in females, especially with late onset, we stratified cases into three groups according to sex and age of onset. The effect of this stratification is highly suggestive despite its crudeness (Table 3) and small sample sizes. Females with onset ≤39 account for most of the association. The other two classes give such small chi-square values that they would undoubtedly be assigned to other regions if the partition test had not been restricted to region 6 on the pooled evidence. However, when considering region 6 alone, there is remarkable agreement between point estimates for 'early' and 'late' onset females and those from males. At this time it is impossible to say whether this consistency is caused by imperfectly divided onset groups or a small effect at late age.

Table 3

Stratification by gender and age of onset (5000 replicates)

Sex	Age of onset	Case/Control		S		K	SE	$χ_{1}^{2}$	LOD₁	p
				kb	LDU
M	-	95	95	53308.370	1195.126	3.13	0.57	0.0299	0.007	0.8624
F	≤39	193	192	53308.370	1195.126	47.3	0.15	11.844	2.575	0.0006
F	≥40	170	173	53262.482	1194.474	0.04	5.18	0.0183	0.004	0.8589

Linkage

Choi et al. [13] reported a meta-analysis of four linkage studies with microsatellites in a 10-Mb bin of chromosome 18. The results from this study were reported as p-values without estimates of location or standard errors. Without this information, the power for meta-analysis is reduced because the sum of two $χ_{1}^{2}$ values must be converted back to $χ_{1}^{2}$ and LOD₁ instead of weighting estimates of location by their information. Perhaps because of this inefficiency, the combined LOD₁ from this meta-analysis is 1.542, well below the conventional value of 3 for asserting significance. The corresponding p-value in large-sample theory is 0.007714, providing strong but inconclusive evidence for localization in the 18q21 region. Despite its limitations, linkage contributes evidence that should not be ignored.

Joint significance of linkage and association

The simplest meta-analysis is based on n independent samples, the i^th of which contributes a P_ivalue that on the null hypothesis is uniformly distributed. Then -2 ln P_iwould be distributed as

χ_{2}^{2}

, with

χ_{2 n}^{2} = - 2 \sum \ln P_{i}

. This is the only test applicable to data that do not provide an estimate of location S_iand information K_i, but has three disadvantages; first, equal weight is given to samples with different standard errors; second, there is no test of homogeneity; and third, there is no point estimate to become more precise as n increases. As a consequence, much information is lost. Accepting these limitations and assuming accuracy of the P estimates, Table 4 shows that combining pooled association with linkage provides suggestive evidence to assign a gene for rheumatoid arthritis to the 18q21.31 interval. The LOD₁ with no Bonferroni correction is 2.676 for linkage and pooled association. When location and information weight are available, the evidence for association is combined by determination of the difference between

Σ χ_{1}^{2}

with n degrees of freedom and

Σ K_{i} {({\hat{S}}_{i} - \bar{S})}^{2}

, which tests for heterogeneity with n - 1 degrees freedom where

\bar{S} = Σ K_{i} S_{i} / Σ K_{i}

. When the stratified association samples are combined in this manner, the heterogeneity test is negligible. As expected, power is increased when pooled with linkage (LOD₁ = 3.401, p = 0.000076). Even with conservative adjustment of the p-value to account for the 18 regions tested by association (18*0.000569), and despite strong although not formally significant, evidence from linkage for at least one causal gene in the 18 regions, the meta-analysis is supportive (LOD₁ = 2.327, p = 0.001062). We conclude that evidence for region 6 is probative, with linkage and association both providing critical evidence despite lack of a point estimate and information weight for linkage.

Table 4

Meta-analysis of association (5000 replicates) and linkage

Study	Source	$χ_{1}^{2}$	LOD₁	p
1	Choi et al. [13]	7.101	1.542	0.007714
2	Association (region 6)	8.3245	1.808	0.003916
Pooled 1 and 2		12.326	2.676	0.000447
Association (stratified) using P		7.0896	1.539	0.007762
Association (stratified) using S, K		11.875	2.579	0.000569
Pooled with 1		15.662	3.401	0.000076
Association (stratified) using S, K and adjusted		6.5944	1.432	0.010242
Pooled with 1		10.718	2.327	0.001062

Discussion

This application demonstrates that CHROMSCAN is a powerful approach for gene mapping in complex inheritance, which is applicable to meta-analysis. Obvious extensions include identification of a causal locus and more precise definition of the phenotype associated with it. The 95% confidence interval, given by S ± 1.96 (SE), covers 36 kb between 53296 and 53332 kb and includes the msSNP rs3745064. Although no described genes are within this region, it does include four human mRNAs from GenBank: CR590917, AK021217, AK124558, and BC01314, all to the left of point estimate (S). Of these, CR590917 appears to be the most interesting because it is expressed within T cells and could therefore conceivably affect risk for rheumatoid arthritis. Finally, geneid [14] and Genscan [15] predict a similar gene, which is the closest annotated sequence to the point estimate (S). However, nothing is known about the function of this gene and its reliability is questionable. The fascinating directions revealed by these findings have yet to be explored. Ultimately, interaction with other contributing loci and environmental factors will be recognized and, more importantly, locus-specific treatment will be found.

Recent papers testify to growing interest in meta-analysis, looking backward to linkage rather than forward to association mapping. Rank permutation provides a valid significance test, but the genome search meta-analysis (GSMA) that uses regional assignment with arbitrary weights cannot give a reliable estimate of effect and therefore has low power for estimating point location and detecting heterogeneity [16, 17]. Most of the few papers on association mapping assume family data rarely feasible for diseases of late onset and are restricted to single markers without composite likelihood to estimate both location S and its information K. One manuscript presented in GAW15 that used meta-analysis without those estimates failed to detect the strong signal on chromosome 18q demonstrated by composite likelihood [18].

Declarations

Acknowledgements

This article has been published as part of BMC Proceedings Volume 1 Supplement 1, 2007: Genetic Analysis Workshop 15: Gene Expression Analysis and Approaches to Detecting Multiple Functional Loci. The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/1?issue=S1.

Authors’ Affiliations

(1)

Human Genetics Division, University of Southampton, Southampton General Hospital, Southampton, UK

References

Shuman S: Structure, mechanism, and evolution of the mRNA capping apparatus. Prog Nucleic Acid Res Mol Biol. 2000, 66: 1-40.View ArticleGoogle Scholar
Morton NE: Sequential tests for the detection of linkage. Am J Hum Genet. 1955, 7: 277-318.PubMed CentralPubMedGoogle Scholar
Levinson DF, Levinson MD, Segurado R, Lewis CM: Genome scan meta-analysis of schizophrenia and bipolar disorder. Part I: Methods and power analysis. Am J Hum Genet. 2003, 73: 17-33. 10.1086/376548.View ArticlePubMed CentralPubMedGoogle Scholar
Morton NE, Maniatis N, Zhang W, Ennis S, Collins A: Genome scanning by composite likelihood. Am J Hum Genet. 2007, 80: 19-28. 10.1086/510401.View ArticlePubMed CentralPubMedGoogle Scholar
CHROMSCAN. [http://www.som.soton.ac.uk/research/geneticsdiv/epidemiology/chromscan/]
Amos CI, Chen WV, Lee A, Li W, Kern M, Lundsten R, Batliwalla F, Wener M, Remmers E, Kastner DA, Chrisiwell LA, Seldin MF, Gregersen PK: High-density SNP analysis of 642 Caucasian families with rheumatoid arthritis identifies two new linkage regions in 11p12 and 2q33. Genes Immun. 2006, 7: 277-286. 10.1038/sj.gene.6364295.View ArticlePubMedGoogle Scholar
Gomes I, Collins A, Lonjou C, Thomas NS, Wilkinson J, Watson M, Morton N: Hardy-Weinberg quality control. Ann Hum Genet. 1999, 63: 535-538. 10.1046/j.1469-1809.1999.6360535.x.View ArticlePubMedGoogle Scholar
Maniatis N, Collins A, Xu CF, McCarthy LC, Hewett DR, Tapper W, Ennis S, Ke X, Morton NE: The first linkage disequilibrium (LD) maps: delineation of hot and cold blocks by diplotype analysis. Proc Natl Acad Sci USA. 2002, 99: 2228-2233. 10.1073/pnas.042680999.View ArticlePubMed CentralPubMedGoogle Scholar
Collins A, Morton NE: Mapping a disease locus by allelic association. Proc Natl Acad Sci USA. 1998, 95: 1741-1745. 10.1073/pnas.95.4.1741.View ArticlePubMed CentralPubMedGoogle Scholar
Maniatis N, Morton NE, Gibson J, Xu CF, Hosking LK, Collins A: The optimal measure of linkage disequilibrium reduces error in association mapping of affection status. Hum Mol Genet. 2005, 14: 145-153. 10.1093/hmg/ddi019.View ArticlePubMedGoogle Scholar
Morton NE, Zhang W, Taillon-Miller P, Ennis S, Kwok PY, Collins A: The optimal measure of allelic association. Proc Natl Acad Sci USA. 2001, 98: 5217-5221. 10.1073/pnas.091062198.View ArticlePubMed CentralPubMedGoogle Scholar
Lau W, Kuo TY, Tapper W, Cox S, Collins A: Exploiting large scale computing to construct high resolution linkage disequilibrium maps of the human genome. Bioinformatics. 2007, 23: 517-519. 10.1093/bioinformatics/btl615.View ArticlePubMedGoogle Scholar
King RA, Rotter JI, Motulsky AG: The Genetic Basis of Common Disease. 1992, New York: Oxford University Press, 598-599.Google Scholar
Choi SJ, Rho YH, Ji JD, Song GG, Lie YH: Genome scan meta-analysis of rheumatoid arthritis. Rheumatology. 2006, 45: 166-170. 10.1093/rheumatology/kei128.View ArticlePubMedGoogle Scholar
Blanco E, Parra G, Guigó R: Using geneid to identify genes. Current Protocols in Bioinformatics. Edited by: Baxevanis AD, Davison DB. 2002, New York: John Wiley & Sons Inc, 1-26.Google Scholar
Burge C, Karlin S: Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997, 268: 78-94. 10.1006/jmbi.1997.0951.View ArticlePubMedGoogle Scholar
Zintzaras E, Kitsios G: Identification of chromosomal regions linked to premature myocardial infarction: a meta-analysis of whole-genome searches. J Hum Genet. 2006, 51: 1015-1021. 10.1007/s10038-006-0053-x.View ArticlePubMedGoogle Scholar
Lewis CM, Levinson DF: Testing for genetic heterogeneity in the genome search meta-analysis method. Genet Epidemiol. 2006, 30: 348-355. 10.1002/gepi.20149.View ArticlePubMedGoogle Scholar
Segurado R, Hamshere ML, Glaser B, Nikolov I, Moskvina V, Holmans P: Combining linkage datasets for meta-analysis and mega-analysis: the GAW15 rheumatoid arthritis data set. BMC Genet. 2007, xx: xxx-Google Scholar

Copyright

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Genetic Analysis Workshop 15: Gene Expression Analysis and Approaches to Detecting Multiple Functional Loci

Mapping a gene for rheumatoid arthritis on chromosome 18q21

Abstract

Keywords

Background

Methods

Data preparation

LDU map construction

Association mapping

Results

Association mapping

Linkage

Joint significance of linkage and association

Discussion

Declarations

Acknowledgements

Authors’ Affiliations

References

Copyright

BMC Proceedings

Contact us