Faculty Bibliography

For admixture mapping studies in Mexican Americans (MAM), we define a genomewide single-nucleotide-polymorphism (SNP) panel that can distinguish between chromosomal segments of Amerindian (AMI) or European (EUR) ancestry. These studies used genotypes for >400,000 SNPs, defined in EUR and both Pima and Mayan AMI, to define a set of ancestry-informative markers (AIMs). The use of two AMI populations was necessary to remove a subset of SNPs that distinguished genotypes of only one AMI subgroup from EUR genotypes. The AIMs set contained 8,144 SNPs separated by a minimum of 50 kb with only three intermarker intervals >1 Mb and had EUR/AMI FST values >0.30 (mean FST = 0.48) and Mayan/Pima FST values <0.05 (mean FST < 0.01). Analysis of a subset of these SNP AIMs suggested that this panel may also distinguish ancestry between EUR and other disparate AMI groups, including Quechuan from South America. We show, using realistic simulation parameters that are based on our analyses of MAM genotyping results, that this panel of SNP AIMs provides good power for detecting disease-associated chromosomal segments for genes with modest ethnicity risk ratios. A reduced set of 5,287 SNP AIMs captured almost the same admixture mapping information, but smaller SNP sets showed substantial drop-off in admixture mapping information and power. The results will enable studies of type 2 diabetes, rheumatoid arthritis, and other diseases among which epidemiological studies suggest differences in the distribution of ancestry-associated susceptibility

Gene-gene and gene-environment interactions are key features in the development of rheumatoid arthritis (RA) and other complex diseases. The aim of this study was to use and compare three different definitions of interaction between the two major genetic risk factors of RA--the HLA-DRB1 shared epitope (SE) alleles and the PTPN22 R620W allele--in three large case-control studies: the Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA) study, the North American RA Consortium (NARAC) study, and the Dutch Leiden Early Arthritis Clinic study (in total, 1,977 cases and 2,405 controls). The EIRA study was also used to analyze interactions between smoking and the two genes. 'Interaction' was defined either as a departure from additivity, as interaction in a multiplicative model, or in terms of linkage disequilibrium--for example, deviation from independence of penetrance of two unlinked loci. Consistent interaction, defined as departure from additivity, between HLA-DRB1 SE alleles and the A allele of PTPN22 R620W was seen in all three studies regarding anti-CCP-positive RA. Testing for multiplicative interactions demonstrated an interaction between the two genes only when the three studies were pooled. The linkage disequilibrium approach indicated a gene-gene interaction in EIRA and NARAC, as well as in the pooled analysis. No interaction was seen between smoking and PTPN22 R620W. A new pattern of interactions is described between the two major known genetic risk factors and the major environmental risk factor concerning the risk of developing anti-CCP-positive RA. The data extend the basis for a pathogenetic hypothesis for RA involving genetic and environmental factors. The study also raises and illustrates principal questions concerning ways to define interactions in complex diseases

Ig class switch recombination (CSR) and somatic hypermutation serve to diversify antibody responses and are orchestrated by the activity of activation-induced cytidine deaminase and many proteins involved in DNA repair and genome surveillance. Msh5, a gene encoded in the central MHC class III region, and its obligate heterodimerization partner Msh4 have a critical role in regulating meiotic homologous recombination and have not been implicated in CSR. Here, we show that MRL/lpr mice carrying a congenic H-2(b/b) MHC interval exhibit several abnormalities regarding CSR, including a profound deficiency of IgG3 in most mice and long microhomologies at Ig switch (S) joints. We found that Msh5 is expressed at low levels on the H-2(b) haplotype and, importantly, a similar long S joint microhomology phenotype was observed in both Msh5 and Msh4-null mice. We also present evidence that genetic variation in MSH5 is associated with IgA deficiency and common variable immune deficiency (CVID) in humans. One of the human MSH5 alleles identified contains two nonsynonymous polymorphisms, and the variant protein encoded by this allele shows impaired binding to MSH4. Similar to the mice, Ig S joints from CVID and IgA deficiency patients carrying disease-associated MSH5 alleles show increased donor/acceptor microhomology, involving pentameric DNA repeat sequences and lower mutation rates than controls. Our findings suggest that Msh4/5 heterodimers contribute to CSR and support a model whereby Msh4/5 promotes the resolution of DNA breaks with low or no terminal microhomology by a classical nonhomologous end-joining mechanism while possibly suppressing an alternative microhomology-mediated pathway

We tested the hypothesis that de novo copy number variation (CNV) is associated with autism spectrum disorders (ASDs). We performed comparative genomic hybridization (CGH) on the genomic DNA of patients and unaffected subjects to detect copy number variants not present in their respective parents. Candidate genomic regions were validated by higher-resolution CGH, paternity testing, cytogenetics, fluorescence in situ hybridization, and microsatellite genotyping. Confirmed de novo CNVs were significantly associated with autism (P = 0.0005). Such CNVs were identified in 12 out of 118 (10%) of patients with sporadic autism, in 2 out of 77 (3%) of patients with an affected first-degree relative, and in 2 out of 196 (1%) of controls. Most de novo CNVs were smaller than microscopic resolution. Affected genomic regions were highly heterogeneous and included mutations of single genes. These findings establish de novo germline mutation as a more significant risk factor for ASD than previously recognized

Systematic genome-wide studies to map genomic regions associated with human diseases are becoming more practical. Increasingly, efforts will be focused on the identification of the specific functional variants responsible for the disease. The challenges of identifying causal variants include the need for complete ascertainment of genetic variants and the need to consider the possibility of multiple causal alleles. We recently reported that risk of systemic lupus erythematosus (SLE) is strongly associated with a common SNP in IFN regulatory factor 5 (IRF5), and that this variant altered spicing in a way that might provide a functional explanation for the reproducible association to SLE risk. Here, by resequencing and genotyping in patients with SLE, we find evidence for three functional alleles of IRF5: the previously described exon 1B splice site variant, a 30-bp in-frame insertion/deletion variant of exon 6 that alters a proline-, glutamic acid-, serine- and threonine-rich domain region, and a variant in a conserved polyA+ signal sequence that alters the length of the 3' UTR and stability of IRF5 mRNAs. Haplotypes of these three variants define at least three distinct levels of risk to SLE. Understanding how combinations of variants influence IRF5 function may offer etiological and therapeutic insights in SLE; more generally, IRF5 and SLE illustrates how multiple common variants of the same gene can together influence risk of common disease

High Mobility Group Box-1 (HMGB1) is a cytokine implicated in the pathogenesis of rheumatoid arthritis (RA) and other inflammatory diseases. The cholinergic anti-inflammatory pathway, a vagus nerve-dependent mechanism, inhibits HMGB1 release in experimental disease models. Here, we examine the relationship between vagus nerve activity and HMGB1 in patients with RA. We compared RR interval variability, an index of cardiac vagal modulation, HMGB1 and hsCRP serum levels, and disease activity scores in thirteen RA patients and eleven age- and sex-matched controls. In RA patients, serum levels of HMGB1 and hsCRP were elevated as compared with controls (HMGB1=71 ng/mL [45-99] vs. 18 ng/mL [0-40], P<0.0001; hsCRP=14.5 mg/L [0.7-59] vs. 1 mg/L [0.4-2.9], P<0.001). RR interval variability in RA patients was significantly decreased as compared with controls (HF=38 msec2 [14-80] vs. 288 msec2 [38-364], P<0.0001; rMSSD=20.9+/-9.79 msec, 52.6+/-35.3 msec, P<0.01). HMGB1 levels and RR interval variability were significantly related (rho=-0.49, P<0.01). HMGB1 serum levels significantly correlated with disease activity scores (DAS-28) in patients with RA (P=0.004). The study design does not enable a determination of causality, but the results are consistent with the hypothesis that decreased cholinergic anti-inflammatory pathway activity is associated with increased HMGB1 levels in patients with RA

Recent studies have shown increased expression of interferon (IFN)-regulated genes in the peripheral blood cells of patients with systemic lupus erythematosus. A similar interferon signature has been observed in affected muscle tissue from patients with dermatomyositis (DM), but it has not yet been determined if this signature extends to the peripheral blood in DM. We performed global gene expression profiling of peripheral blood cells from adult and juvenile DM patients and healthy controls. Several interesting groups of genes were differentially expressed in DM, including genes with immune function, and others that function in muscle or are involved in mitochondrial/oxidative phosphorylation. Investigation of type I IFN-regulated transcripts revealed a striking interferon signature present in most DM patients studied. Levels of type I IFN-regulated proteins were also elevated in DM serum samples. Furthermore, both the transcript and serum protein IFN signatures were associated with disease activity. These data suggest that the IFN signature may be a useful marker for DM disease activity, and that sampling peripheral blood may be a more practical alternative to muscle biopsy for measuring this signature

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

ABSTRACT : 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