Faculty Bibliography

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

MOTIVATION: Genetic association analysis is based on statistical correlations which do not assign any cause-to-effect arrows between the two correlated variables. Normally, such assignment of cause and effect label is not necessary in genetic analysis since genes are always the cause and phenotypes are always the effect. However, among intermediate phenotypes and biomarkers, assigning cause and effect becomes meaningful, and causal inference can be useful. RESULTS: We show that causal inference is possible by an example in a study of rheumatoid arthritis. With the help of genotypic information, the shared epitope, the causal relationship between two biomarkers related to the disease, anti-cyclic citrullinated peptide (anti-CCP) and rheumatoid factor (RF) has been established. We emphasize the fact that third variable must be a genotype to be able to resolve potential ambiguities in causal inference. Two non-trivial conclusions have been reached by the causal inference: (1) anti-CCP is a cause of RF and (2) it is unlikely that a third confounding factor contributes to both anti-CCP and RF

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

Systemic lupus erythematosus (SLE) is a complex autoimmune disease characterized by activation of the type I interferon (IFN) pathway. Here we convincingly replicate association of the IFN regulatory factor 5 (IRF5) rs2004640 T allele with SLE in four independent case-control cohorts (P = 4.4 x 10(-16)) and by family-based transmission disequilibrium test analysis (P = 0.0006). The rs2004640 T allele creates a 5' donor splice site in an alternate exon 1 of IRF5, allowing expression of several unique IRF5 isoforms. We also identify an independent cis-acting variant associated with elevated expression of IRF5 and linked to the exon 1B splice site. Haplotypes carrying the variant associated with elevated expression and lacking the exon 1B donor site do not confer risk of SLE. Thus, a common IRF5 haplotype driving elevated expression of multiple unique isoforms of IRF5 is an important genetic risk factor for SLE, establishing a causal role for type I IFN pathway genes in human autoimmunity

OBJECTIVE: HLA is the most strongly associated locus in rheumatoid arthritis (RA), accounting for up to one-third of the genetic contribution. Conditioning on the effect of true disease loci such as HLA can lead to increased power to detect effects at other loci and, in addition, allows investigation of the underlying disease models, including interactions. The aim of this study was to detect susceptibility loci for RA by conditioning on HLA in a large sample of affected sibling pairs (ASPs) and to test for evidence of interaction between novel loci and HLA. METHODS: Genotype data from 3 whole-genome linkage scans for RA in a US population and a UK population were pooled, resulting in a combined data set of 886 ASPs. This pooling of data increased the power to detect loci showing low levels of heterogeneity. Nonparametric linkage analysis was performed to identify regions of interest. Joint 2-locus analysis was then performed for HLA and each of the loci that demonstrated evidence of linkage in the 886 ASPs. RESULTS: Evidence for linkage was most significant at HLA (P = 4 x 10(-16)), with 7 non-HLA loci showing some evidence for linkage (P = 0.05-0.003). Joint modeling of these loci with HLA provided evidence for linkage at a genome-wide significance level for loci on 6q (P = 2.7 x 10(-6)) and 16p (P = 2 x 10(-4)). CONCLUSION: These data provide the most convincing evidence to date that 6q and 16p harbor susceptibility genes. In addition, these loci may interact with HLA, facilitating the search for candidate genes within this region.

Human autoimmune diseases are well suited for the application of gene expression profiling. Sampling of blood cells and target tissues has already revealed many important pathways contributing to this spectrum of disorders, and many commonalities are emerging. For instance, clinically distinct diseases such as systemic lupus erythematosus, Sjogren's syndrome, dermatomyositis, and psoriasis all show evidence for dysregulation of the type I interferon pathway. These data suggest that autoimmune diseases will eventually be categorized at the level of gene expression. This work has led to advances in our understanding of disease pathogenesis and in the future promises to facilitate assessments of disease activity and improve targeting of therapies. Here, we review the literature on gene profiling in human autoimmune diseases and provide perspective on the current state of the art

OBJECTIVE: To test for association of the ectonucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1) K121Q polymorphism with body mass index (BMI) and diabetes in a large sample of Caucasians and African-Americans by selectively genotyping individuals at the extremes of the phenotypic distribution. SUBJECTS: Subsets comprising the extremes of the BMI distribution (10th-20th and above the 90th BMI percentile for Caucasians and between the 10th-30th and above the 80th percentile for African-Americans) from a group of 10,260 Caucasian and 2268 African-American adults participating in New York Cancer Project were studied. METHODS: Subjects were genotyped for the ENPP1 K121Q polymorphism by pyrosequencing and tested for association with BMI and diabetes by regression analysis. RESULTS: Regression analysis with BMI as the dependent variable demonstrated a significant association (P = 0.02) of genotype at K121Q with BMI, with no significant race-by-genotype interaction (P = 0.30). Compared with Q/Q or Q/K individuals, the K/K individuals had a BMI approximately 1.3 kg/m2 higher, without effects of age, gender or race. By logistic regression analysis, the K121Q alleles had no significant effect on diabetes status (P = 0.37) in obese subjects. CONCLUSION: In both Caucasians and African-Americans, the K121 polymorphism in ENPP1 was associated with increased BMI, but not with diabetes

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