Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Osteopontin and Systemic Lupus Erythematosus Association: A Probable Gene-Gender Interaction

  • Shizhong Han,

    Affiliation Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America

  • Joel M. Guthridge,

    Affiliation Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America

  • Isaac T. W. Harley,

    Affiliation Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America

  • Andrea L. Sestak,

    Affiliation Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America

  • Xana Kim-Howard,

    Affiliation Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America

  • Kenneth M. Kaufman,

    Affiliations Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America, United States Department of Veterans Affairs Medical Center, Oklahoma City, Oklahoma, United States of America, Department of Medicine, University of Oklahoma Health Sciences Center, Oklahoma City, Oklahoma, United States of America

  • Bahram Namjou,

    Affiliation Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America

  • Harshal Deshmukh,

    Affiliation Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America

  • Gail Bruner,

    Affiliation Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America

  • Luis R. Espinoza,

    Affiliation Louisiana State University Health Science Center, New Orleans, Louisiana, United States of America

  • Gary S. Gilkeson,

    Affiliation Medical University of South Carolina, Charleston, South Carolina, United States of America

  • John B. Harley,

    Affiliations Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America, United States Department of Veterans Affairs Medical Center, Oklahoma City, Oklahoma, United States of America, Department of Medicine, University of Oklahoma Health Sciences Center, Oklahoma City, Oklahoma, United States of America

  • Judith A. James,

    Affiliations Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America, Department of Medicine, University of Oklahoma Health Sciences Center, Oklahoma City, Oklahoma, United States of America

  • Swapan K. Nath

    swapan-nath@omrf.ouhsc.edu

    Affiliations Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma, United States of America, Department of Medicine, University of Oklahoma Health Sciences Center, Oklahoma City, Oklahoma, United States of America

Abstract

Osteopontin (SPP1) is an important bone matrix mediator found to have key roles in inflammation and immunity. SPP1 genetic polymorphisms and increased osteopontin protein levels have been reported to be associated with SLE in small patient collections. The present study evaluates association between SPP1 polymorphisms and SLE in a large cohort of 1141 unrelated SLE patients [707 European-American (EA) and 434 African-American (AA)], and 2009 unrelated controls (1309 EA and 700 AA). Population-based case-control association analyses were performed. To control for potential population stratification, admixture adjusted logistic regression, genomic control (GC), structured association (STRAT), and principal components analysis (PCA) were applied. Combined analysis of 2 ethnic groups, showed the minor allele of 2 SNPs (rs1126616T and rs9138C) significantly associated with higher risk of SLE in males (P = 0.0005, OR = 1.73, 95% CI = 1.28–2.33), but not in females. Indeed, significant gene-gender interactions in the 2 SNPs, rs1126772 and rs9138, were detected (P = 0.001 and P = 0.0006, respectively). Further, haplotype analysis identified rs1126616T-rs1126772A-rs9138C which demonstrated significant association with SLE in general (P = 0.02, OR = 1.30, 95%CI 1.08–1.57), especially in males (P = 0.0003, OR = 2.42, 95%CI 1.51–3.89). Subgroup analysis with single SNPs and haplotypes also identified a similar pattern of gender-specific association in AA and EA. GC, STRAT, and PCA results within each group showed consistent associations. Our data suggest SPP1 is associated with SLE, and this association is especially stronger in males. To our knowledge, this report serves as the first association of a specific autosomal gene with human male lupus.

Introduction

Systemic Lupus Erythematosus (SLE) is a prototypic human autoimmune disease characterized by impaired T cell responses, dysregulated B cell activation, hyperactive B cells and autoantibody production leading to inflammation and potential end-organ damage. While the etiology of SLE remains complex, genetic factors are known to be important in the pathogenesis of SLE [1], [2]. The current collection of genetic information suggests that SLE susceptibility arises from specific combinations of multiple gene-gene and gene-environment interactions. Among the genetic factors believed to influence SLE susceptibility, the major histocompatibility complex (MHC) alleles show the most significant association, but these do not explain the total genetic background of the disease. Importantly, several recent studies show that non-HLA genes play a role in SLE development [3][7]. Recently, several lines of evidence suggest that secreted phosphoprotein 1 (SPP1) located at 4q22, also called osteopontin and early T-lymphocyte activation 1, may have a role in the pathogenesis of SLE as well as other autoimmune disorders.

SPP1 plays a key role in bone biology and has recently found to also be important in regulating inflammation and immunity. The immunologic functions of SPP1 include enhancing the proinflammatory Th1 cell response and inhibiting the Th2 responses [8][9]. In addition, some studies have suggested that SPP1 plays a role in the survival of activated T cells by inducing apoptosis, while others have demonstrated the essential role of an intracellular form of SPP1 in the production of interferon-alpha by plasmacytoid dendritic cells [10], [11]. Humans with SLE and autoimmune prone mice over express osteopontin suggesting that abnormal expression of this protein may participate in SLE disease pathogenesis [12], [13]. Further, polymorphic osteopontin alleles have been implicated in the development of a mouse model of lupus [14].

SNPs in the SPP1 gene have also been reported to be associated with human SLE, adding further support to the role of this gene in SLE pathogenesis [15]. A significant association between rs11226616 and SLE was first demonstrated in a small North American Caucasian cohort study [15]. Two SNPs (rs1126772 and rs9138) in the 3′ UTR in the SPP1 gene were associated with high levels of SPP1 and elevated risk of developing autoimmune/lymphoproliferative syndrome (ALPS), a disorder which leads to an autoimmune pattern similar to lupus prone strains of mice [16]. The same group later showed significant associations between SLE and 2 SPP1 SNPs (rs7687316 and rs9138) in an Italian population [17]. This information prompted us to test association between SPP1 polymorphisms and SLE in a large, multi-ethnic collection.

Results

Marker information, minor allele frequency and the statistical significance for allele distributions between cases and controls are presented in Table 1. The only significant difference in allele distribution was observed in the combined male-female group for rs6840362, which showed a significant difference allele distribution in EA (P = 0.015). However, significant differences were evident in the male subgroup. This finding was especially strong in EA males, where 3 SNPs (rs1126616, rs1126772 and rs9138) showed significant differences in allele distributions. Similarly, in AA males, 2 SNPs (rs1126616 and rs9138) demonstrated significant differences in allele distribution.

thumbnail
Table 1. Marker information and minor allele frequency in African-American and European-American sample.

https://doi.org/10.1371/journal.pone.0001757.t001

We evaluated the association of each polymorphism with SLE, adjusting for the admixture proportion utilizing logistic regression under the multiplicative genetic model for minor alleles by combined analysis. Considering the gender effect in allele distributions and possible race specific effect in disease susceptibility, subgroup analysis stratified by gender and race were also performed. Table 2 shows the association results in detail. Briefly, 2 SNPs' minor alleles (rs1126616T and rs9138C) showed significant associations with SLE in AA and EA combined males, but not in females, both of which conferred a high risk of SLE (P = 0.0005, OR = 1.73, 95%CI 1.28–2.33). Indeed, significant gene-gender interactions in the 2 SNPs, rs1126772 and rs9138, were detected (P = 0.001, P = 0.0006, respectively). Subgroup analysis by race revealed the same trend in AA and EA. For AA, 3 SNPs (rs11728697, rs1126616, and rs9138) showed significant association in males only (P = 0.02, P = 0.027, P = 0.027, respectively). In EA males, 3 SNPs (rs1126616, rs1126772, and rs9138) also showed significant association (P = 0.003, P = 0.028, P = 0.003). Furthermore, to exclude the false positive association which can arise from hidden population substructure, we utilized GC, STRAT, and PCA to verify the association result in each population. All associations remained consistent with the admixture adjusted logistic regression analysis (Table 3).

thumbnail
Table 2. Associations of SNPs with SLE by logistic regression under multiplicative genetic model

https://doi.org/10.1371/journal.pone.0001757.t002

thumbnail
Table 3. Associations of SNPs with SLE by genomic control (GC), structured association test (STRAT), and principle components analysis (PCA).

https://doi.org/10.1371/journal.pone.0001757.t003

To assess the effect of any particular clinical feature on the genetic association, we have performed a subgroup analysis classified by eight clinical characteristics (Cutaneous manifestations, arthritis, serositis, renal involvement, neuro-psychiatric manifestations, hematological features, anti-dsDNA, antinuclear antibody) available for males (91 male SLE) and females (754 female SLE) separately. Although we did not find any significant association in females, we did observed evidence of association in males for some of the clinical features; however, these associations are not stronger than the overall male-specific analysis. Therefore, there seems to be no evidence that the overall effect is dominated by a particular subset.

Haplotype analysis was performed to further evaluate the role of SPP1 in SLE susceptibility. We conducted haplotype reconstruction and linkage disequilibrium analysis incorporating all 11 SNPs in the SPP1 gene. The LD map and each SNP associated P value are depicted in Fig. 1. Notably, rs1126616 and rs9138 are almost in complete LD in both AA (Fig. 1A) and EA (Fig. 1B), which explained their similar behavior in the association with SLE. Initially, we included all selected 11 SNPs for haplotype analysis in combined data and each race specific population adjusted by admixture proportion. In accordance with the results of single SNP analysis, no significant global association was detected either in the combined analysis or in each specific population. The same trend remained in the female subgroup. However, significant haplotype association was observed in the male subgroup for the combined, AA and EA (P = 0.001, P = 0.003, P = 0.028, respectively).

thumbnail
Figure 1. Linkage disequilibrium (LD) structure and association P value in AA.

LD is calculated from the genotype data in AA healthy control population. Fig. 1B: LD structure and association P value in EA. LD is calculated from the genotype data in EA healthy control population.

https://doi.org/10.1371/journal.pone.0001757.g001

To further explore the haplotype effect in male subgroup, we performed conditional analysis to detect if there was a specific SNP or a subset of SNPs which can explain the global haplotype association. Interestingly, the last 3 SNPs (rs1126616, rs1126772, and rs9138) explained the whole association for each significant group. When we performed haplotype association analysis conditional on these 3 SNPs, all significant associations disappear in the combined AA and EA subgroup. We next focused on these 3 SNPs for haplotype analysis. Detailed results are presented in Table 4. Although single SNP analysis revealed only marginal significant association, it is notable that in AA males the haplotype analysis demonstrated a much stronger association (P = 0.002). Specifically, the haplotype TAC confers a high risk of SLE in AA males (P = 0.001, OR = 3.37, 95%CI 1.65–6.92). In EA males, the TAC haplotype is not significant, but there is a trend in the same direction towards association with SLE (P = 0.079, OR = 1.84, 95%CI 0.96–3.53). Furthermore, in the combined male subgroup, the haplotype TAC significantly increases the risk of SLE (P = 0.0004, OR = 2.42, 95%CI 1.51–3.89). All the association P-values were also verified by 10,000 permutations. Accordingly, interaction analysis between haplotype and gender also showed significant interactions in AA, EA, and the combined analysis (P = 0.018, P = 0.02, P = 0.0018, respectively).

thumbnail
Table 4. Haplotypic analysis of SPP1 polymorphisms (rs1126616, rs1126772, and rs9138) and SLE risk.

https://doi.org/10.1371/journal.pone.0001757.t004

Discussion

Our study confirms the previously reported genetic association with SLE and presents additional support in a large multiethnic cohort. Previous studies have suggested that increased SPP1 plasma concentration, as a result of increased gene/protein expression and local production, was associated with SLE [13]. Therefore, SPP1 is a reasonable candidate gene for SLE susceptibility. In the first report, Forton et al found that a silent polymorphism (rs1126616) in exon 7 was significantly associated with SLE [15]. Subsequently, 2 SNPs in the 5′ (rs7687316) and 3′ (rs9138) ends of the SPP1 gene were reported to contribute to SLE susceptibility [16]. Although our combined male-female results were marginally significant, associations with SLE were found for haplotype analysis of 3 SNPs (rs1126616, rs1126772, and rs9138). Since significant interactions were detected between gender and the SPP1 SNPs, haplotypes were analyzed separately for males and females. Our unique findings focus on stronger associations found in male SLE patients from the combined analysis of samples from 2 different ethnic populations. Haplotype analyses revealed that the last 3 SNPs in 3′UTR explain the global association in males and supports the hypothesis that the causal variant of SPP1 might be near 3′UTR, which could affect the expression level of SPP1.

SLE is at least nine times more prevalent in female compared to male subjects; however, the underlying cause of the gender effect has not been clearly established. There is evidence that gender-specific genetic effects exist, both from the many differences in animal models of lupus and from other previous work in humans [18][21]. In a recent study by our group [22], a SLE linkage at 13q32 was identified and replicated by restricting the samples to a relatively homogenous population of African-American families containing at least one male affected. Therefore, taken all the evidences together, the evolving data suggest unique clinical courses and genetic predispositions in male lupus. Further study is warranted to see if differences in osteopontin regulation are involved in these differences.

Our haplotype analyses findings could also be significant due to possible functionality of the three 3′UTR SNPs. Recent years have seen an increased appreciation for the importance of post-transcriptional regulation in eukaryotic organisms [23]. The untranslated region at the 3′ end of a gene (3′UTR) is involved in regulating gene expression at both the pre-mRNA level and the mature mRNA level. In the former, the 3′UTR plays a central role in mRNA 3′ processing and polyadenylation, whereas in the latter cis-elements in the 3′UTR are bound by trans-acting factors which modulate mRNA stability, nuclear export, subcellular localization, and translation efficiency [24][27]. Micro-RNA target sites in 3′UTRs have been shown to be highly conserved [28]. Polymorphisms in microRNA target sites within the 3′UTR may influence gene expression in complex phenotypes, such as lupus [29]. 3′UTRs are also preferred sites of cis-encoded natural antisense transcripts [30]. A search of the PolyA_DB2 database shows that there is supporting cDNA/EST evidence for 2 transcripts with different polyadenylation sites in the SPP1 gene [31]. In addition, the PolymiRTS (polymorphic microRNA Target Site) database revealed evidence that suggests that SNP rs1126772, which shows association with lupus in this report, is predicted to disrupt a non-conserved microRNA target site [32]. Two 3′UTR SNPs, rs1126893, and rs2853754, which are not tested in the present association study, are predicted to disrupt two evolutionarily conserved microRNA target sites. Interestingly, the rs1126893 is in complete LD and within the same haplotype block with our associated SNPs (rs1126772 and rs9138 at the 3′UTR) based on Hapmap data of CEPH families. In fact, the distance between rs1126893 and rs2853754 is only 30 bp, the alternative allele for SNP rs2853754 is predicted to create microRNA target sites for 2 microRNAs. A more extensive analysis of these potential mechanisms for the 3′UTR control of SPP1 gene expression will be required to determine if these predicted sites may be responsible for this genetic association.

However, we acknowledge that the current study suffers from some shortcomings. First, the samples size for the male SLE patients is small (121 males); the association results in male subgroup analysis might arise by chance. However, the possibility of false positives arising from the multiple testing problems should be small since the reported P value is significant enough to survive even the most conservative Bonferroni correction. Second, none of the previously published functional studies were performed to evaluate the role of SPP1 in males and females separately. While this experiment is beyond the scope of the present study, we hope to assess SPP1 gene and protein expression in male and female lupus patients in future. It is worth noting that sex-specific differences in SPP1 gene expression have been observed in rats [33], [34]. Third, an independent replication study should be performed to further verify the male association. Like many other association studies with complex phenotype, no matter how large, validation in a second cohort is needed.

In summary, our data suggest the SPP1 gene might be associated with the development of SLE in general and especially in males. To our knowledge, this report serves as the first description of a gender-specific human lupus genetic association. Further in-depth molecular, genetic, and functional studies should improve our understanding of the disease, and hopefully will provide a more accurate diagnostic algorithm and improved genetic counseling and management strategies.

Materials and Methods

Patients and controls

Genomic DNA samples from SLE patients and control subjects were collected after obtaining written, informed consent from 1141 unrelated SLE patients [707 European-American (EA), 434 African-American (AA)], including 121 males (31 AA male patients and 90 EA male patients), and 2009 unrelated controls (1309 EA, 700 AA), including 598 males (225 AA male controls and 373 EA male controls). Coded DNA samples were obtained from the Lupus Family Registry and Repository (LFRR: http://lupus.omrf.org/). The study was approved by the Institutional Review Boards (IRB) at the Oklahoma Medical Research Foundation and University of Oklahoma Health Sciences Center. All patients utilized in the study met SLE classification based upon the revised criteria of the American College of Rheumatology (ACR) [35].

Genotyping

Initially, we evaluated 32 SNPs in the SPP1 gene on the Illumina SNP platform (San Diego, CA) at the University of Texas Southwestern Microarray Core Facility (Dallas, TX) using standard methods described in more detail (www.illumina.com). A summary of the genotyping information for all the markers is shown in supplementary Table S1. We included 11 SNPs in the final analysis in which the controls in Hardy-Weinberg equilibrium, and a minor allele frequency of >1% in both the populations. Marker information, genotype, and allele distributions of the 11 SNPs in cases and controls are summarized in supplementary Table S2. Of the SNPs that failed quality control, 13 SNPs were monomorphic in 2 studied populations, 6 SNPs had minor allele frequencies which were below 1% at least in one population, 2 SNPs failed genotyping, and 1 SNP deviated significantly from HWE in the black control group.

To know how much of the genetic variation of the SPP1 gene that is captured by the 11 successfully genotyped SNPs, we downloaded the tagged SNPs of SPP1 from the CEU population from Hapmap database (http://www.hapmap.org) using Tagger Pairwise method, in which MAF and LD (r-squared) cutoffs were set to 0.05 and 0.8, respectively. We found 4 tagged SNPs, rs6840362, rs9138, rs2853749, and rs11728697, which can capture all other 10 genotyped SNPs by Hapmap project across the gene. Actually, all these 4 tag SNPs were included in our analyzed marker list. Therefore, all the known genetic variation of the SPP1 gene should be captured by the 11 SNPs we analyzed, especially for the European population.

Statistical analysis

Departures from Hardy–Weinberg equilibrium in cases and controls were tested for SNPs with a Pearson Chi-square test. Allele frequencies in SLE cases were compared to those in control subjects using the Chi-square test or Fisher's exact probability test, where appropriate. Statistical evaluations for testing genetic effects were performed using multivariate logistic regression analysis with adjustments for gender and admixture proportions under multiplicative genetic model. Interaction analysis was performed by introducing the interaction term into the logistic regression model including admixture proportion as a covariate. Statistical significance was obtained by the likelihood ratio test comparing the models with and without the interaction term. Haplotype-based association analysis and conditional analysis were performed by the WHAP program [36].

To control for possible confounding due to population stratification, a panel of 221 ancestry informative markers (AIMs) was genotyped using the same Illumina SNP platform for an ongoing association study project in these LFRR samples. Frequencies of the AIMS in the ancestry population European and African are shown in supplementary Table S3. These AIMs were selected based on the criteria of large allele frequency differences (20% or greater allele frequency difference) between 2 ancestral populations, HWE in ancestral populations and all 221 AIMs were separated by al least 1 cM to minimize the possibility of strong LD between AIMs. We used the STRUCTURE program to estimate the admixture proportion for each individual in AA and EA separately [37]. The log likelihood of each analysis at varying number of population groups (k) was estimated from the average of 3 independent runs (20,000 burn in and 30,000 iterations). As expected, the results favored a two-ancestry population model in both AA and EA. The average proportion of European ancestry was 0.17 in AA samples and 0.99 in EA samples. We included the European ancestry proportion in each individual as a covariate in the logistic regression model to control the population stratification in subgroup and combined analysis. Additionally, genomic control (GC) [38], structured association (STRAT) [39], and principal components analysis (PCA) [40] were applied to control for population stratification in race specific analysis.

Acknowledgments

These authors are grateful to the patients and their families for their cooperation and blood samples.

Author Contributions

Conceived and designed the experiments: SN JG IH AS JJ. Performed the experiments: KK. Analyzed the data: SN SH XK BN HD. Contributed reagents/materials/analysis tools: GG SN GB LE JH. Wrote the paper: SN SH JG JJ.

References

  1. 1. Nath SK, Kilpatrick J, Harley JB (2004) Genetics of human systemic lupus erythematosus: the emerging picture. Curr Opin Immunol 16: 794–800.
  2. 2. Sestak AL, Nath SK, Sawalha AH, Harley JB (2007) Current status of lupus genetics. Arthritis Res Ther 9: 210.
  3. 3. Prokunina L, Castillejo-Lopez C, Oberg F, Gunnarsson I, Berg L, et al. (2002) A regulatory polymorphism in PDCD1 is associated with susceptibility to systemic lupus erythematosus in humans. Nat Genet 32: 666–669.
  4. 4. Kyogoku C, Langefeld CD, Ortmann WA, Lee A, Selby S, et al. (2004) Genetic association of the R620W polymorphism of protein tyrosine phosphatase PTPN22 with human SLE. Am J Hum Genet 75: 504–507.
  5. 5. Kaufman KM, Kelly JA, Herring BJ, Adler AJ, Glenn SB, et al. (2006) Evaluation of the genetic association of the PTPN22 R620W polymorphism in familial and sporadic systemic lupus erythematosus. Arthritis Rheum 54: 2533–2540.
  6. 6. Graham DS, Wong AK, McHugh NJ, Whittaker JC, Vyse TJ (2006) Evidence for unique association signals in SLE at the CD28-CTLA4-ICOS locus in a family-based study. Hum Mol Genet 15: 3195–3205.
  7. 7. Graham RR, Kyogoku C, Sigurdsson S, Vlasova IA, Davies LR, et al. (2007) Three functional variants of IFN regulatory factor 5 IRF5; define risk and protective haplotypes for human lupus. Proc Natl Acad Sci 104: 6758–6763.
  8. 8. O'Regan AW, Nau GJ, Chupp GL, Berman JS (2000) Osteopontin (Eta-1) in cell-mediated immunity: teaching an old dog new tricks. Immunol Today 21: 475–478.
  9. 9. Ashkar S, Weber GF, Panoutsakopoulou V, Sanchirico ME, Jansson M, et al. (2000) Eta-1(osteopontin): an early component of type-1(cell-mediated) immunity. Science 287: 860–864.
  10. 10. Hur EM, Youssef S, Haws ME, Zhang SY, Sobel RA, et al. (2007) Osteopontin-induced relapse and progression of autoimmune brain disease through enhanced survival of activated T cells. Nat Immunol 8: 74–83.
  11. 11. Shinohara ML, Lu L, Bu J, Werneck MB, Kobayashi KS, et al. (2006) Osteopontin expression is essential for interferon-alpha production by plasmacytoid dendritic cells. Nat Immunol 75: 498–506.
  12. 12. Katagiri Y, Mori K, Hara T, Tanaka K, Murakami M, et al. (1995) Functional analysis of the osteopontin molecule. Ann N Y Acad Sci 760: 371–374.
  13. 13. Wong CK, Lit LC, Tam LS, Li EK, Lam CW (2005) Elevation of plasma osteopontin concentration is correlated with disease activity in patients with systemic lupus erythematosus. Rheumatology (Oxford) 44: 602–606.
  14. 14. Miyazaki T, Ono M, Qu WM, Zhang MC, Mori S, et al. (2005) Implication of allelic polymorphism of osteopontin in the development of lupus nephritis in MRL/lpr mice. Eur J Immunol 355: 1510–1520.
  15. 15. Forton AC, Petri MA, Goldman D, Sullivan KE (2002) An osteopontin (SPP1) polymorphism is associated with systemic lupus erythematosus. Hum Mutat 19: 459.
  16. 16. Chiocchetti A, Indelicato M, Bensi T, Mesturini R, Giordano M, et al. (2004) High levels of osteopontin associated with polymorphisms in its gene are a risk factor for development of autoimmunity/lymphoproliferation. Blood 103: 1376–1382.
  17. 17. D'Alfonso S, Barizzone N, Giordano M, Chiocchetti A, Magnani C, et al. (2005) Two single-nucleotide polymorphisms in the 5′ and 3′ ends of the osteopontin gene contribute to susceptibility to systemic lupus erythematosus. Arthritis Rheum 52: 539–547.
  18. 18. Mok CC, Lau CS, Chan TM, Wong RWS (1999) Clinical characteristics and outcome of southern Chinese males with systemic lupus erythematosus. Lupus 8: 188–196.
  19. 19. Molina JF, Drenkard C, Molina J, Cardiel MH, Uribe O, et al. (1996) Systemic lupus erytematosus in males: a study of 107 Latin American patients. Medicine (Baltimore) 75: 124–130.
  20. 20. Stein CM, Olson JM, Gray-McGuire C, Bruner GR, Harley JB, et al. (2002) Increased prevalence of renal disease in systemic lupus erythematosus families with affected male relatives. Arthritis Rheum 46: 428–435.
  21. 21. Arbuckle MR, James JA, Dennis GJ, Rubertone MV, McClain MT, et al. (2003) Rapid clinical progression to diagnosis among African-American men with systemic lupus erythematosus. Lupus 12: 99–106.
  22. 22. Xing C, Gray-McGuire C, Kelly JA, Garriott P, Bukulmez H, et al. (2005) Genetic linkage of systemic lupus erythematosus to 13q32 in African American families with affected male members. Hum Genet 118: 309–321.
  23. 23. Hieronymus H, Silver PA (2004) A systems view of mRNP biology. Genes Dev 18: 2845–2860.
  24. 24. Zhao J, Hyman L, Moore C (1999) Formation of mRNA 3′ ends in eukaryotes: mechanism, regulation, and interrelationships with other steps in mRNA synthesis. Microbiol Mol Biol Rev 63: 405–45.
  25. 25. Conne B, Stutz A, Vassalli JD (2000) The 3′ untranslated region of messenger RNA: a molecular “hotspot' for pathology? Nat Med 6: 637–641.
  26. 26. Mignone F, Gissi C, Liuni S, Pesole G (2002) Untranslated regions of mRNAs. Genome Biol 3: REVIEWS004.
  27. 27. Chabanon H, Mickelburgh I, Hesketh J (2004) Zipcodes and postage stamps: mRNA localization signals and their trans-acting binding proteins. Brief Funct Genomic Proteomic 3: 240–256.
  28. 28. Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, et al. (2005) Systematic discovery of regulatory motifs in human promoters and 3′ UTRs by comparison of several mammals. Nature 434: 338–345.
  29. 29. Martin MM, Buckenberger JA, Jiang J, Malana GE, Nuovo GJ, et al. (2007) The human angiotensin II type 1 receptor +1166 A/C polymorphism attenuates microRNA-155 binding. J Biol Chem 282: 24262–24269.
  30. 30. Sun M, Hurst LD, Carmichael GG, Chen J (2005) Evidence for a preferential targeting of 3′-UTRs by cis-encoded natural antisense transcripts. Nucleci Acids Res 33: 5533–5543.
  31. 31. Zhang H, Hu J, Recce M, Tian B (2005) PolyA_DB: a database for mammalian mRNA polyadenylation. Nucleic Acids Res 33: 201–212.
  32. 32. Bao L, Zhou M, Wu L, Lu L, Godowitz D, et al. (2007) PolymiRTS Database: linking polymorphisms in microRNA target sites with complex traits. Nucleic Acids Res 35: Database Issue D51–D54.
  33. 33. Cowley BD, Ricardo SD, Nagao S, Diamond JR (2001) Increased renal expression of monocyte chemoattractant protein-1 and osteopontin in ADPKD in rats. Kidney International 60: 2087–2096.
  34. 34. Matkovits T, Christakos S (1995) Variable in Vivo Regulation of Rat Vitamin D-Dependent Genes Osteopontin, Ca,Mg-Adenosine Triphosphatase, and 25-Hydroxyvitamin D3 24-Hydroxylase: Implications for Differing Mechanisms of Regulation and Involvement of Multiple Factors. Endocrinology 136: 3971–3982.
  35. 35. Hochberg MC (1997) Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum 40: 1725.
  36. 36. Purcell S, Daly MJ, Sham PC (2007) WHAP: haplotype-based association analysis. Bioinformatics 23: 255–256.
  37. 37. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959.
  38. 38. Devlin B, Roeder K (1999) Genomic control for association studies. Biometrics 55: 997–1004.
  39. 39. Pritchard JK, Stephens M, Rosenberg NA, Donnelly P (2000) Association mapping in structured populations. Am J Hum Genet 67: 170–181.
  40. 40. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, et al. (2006) Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet 38: 904–909.