Mark Seielstad

Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci

To identify new genetic risk factors for rheumatoid arthritis, we conducted a genome-wide association study meta-analysis of 5,539 autoantibody-positive individuals with rheumatoid arthritis (cases) and 20,169 controls of European descent, followed by replication in an independent set of 6,768 rheumatoid arthritis cases and 8,806 controls. Of 34 SNPs selected for replication, 7 new rheumatoid arthritis risk alleles were identified at genome-wide significance (P < 5 × 10−8) in an analysis of all 41,282 samples. The associated SNPs are near genes of known immune function, including IL6ST, SPRED2, RBPJ, CCR6, IRF5 and PXK. We also refined associations at two established rheumatoid arthritis risk loci (IL2RA and CCL21) and confirmed the association at AFF3. These new associations bring the total number of confirmed rheumatoid arthritis risk loci to 31 among individuals of European ancestry. An additional 11 SNPs replicated at P < 0.05, many of which are validated autoimmune risk alleles, suggesting that most represent genuine rheumatoid arthritis risk alleles.

Patterns of linkage disequilibrium in the human genome

Particular alleles at neighbouring loci tend to be co-inherited. For tightly linked loci, this might lead to associations between alleles in the population — a property known as linkage disequilibrium (LD). LD has recently become the focus of intense study in the hope that it might facilitate the mapping of complex disease loci through whole-genome association studies. This approach depends crucially on the patterns of LD in the human genome. In this review, we draw on empirical studies in humans and Drosophila, as well as simulation studies, to assess the current state of knowledge about patterns of LD, and consider the implications for the use of LD as a mapping tool.

Y chromosome sequence variation and the history of human populations

Binary polymorphisms associated with the non-recombining region of the human Y chromosome (NRY) preserve the paternal genetic legacy of our species that has persisted to the present, permitting inference of human evolution, population affinity and demographic history. We used denaturing high-performance liquid chromatography (DHPLC; ref. 2) to identify 160 of the 166 bi-allelic and 1 tri-allelic site that formed a parsimonious genealogy of 116 haplotypes, several of which display distinct population affinities based on the analysis of 1062 globally representative individuals. A minority of contemporary East Africans and Khoisan represent the descendants of the most ancestral patrilineages of anatomically modern humans that left Africa between 35,000 and 89,000 years ago.

Genetic evidence for a higher female migration rate in humans

Mitochondrial DNA and the Y chromosome have been used extensively in the study of modern human origins and other phylogenetic questions, but not in the context of their sex-specific modes of transmission. mtDNA is transmitted exclusively by females, whereas the Y chromosome is passed only among males. As a result, differences in the reproductive output or migration rate of males and females will influence the geographic patterns and relative level of genetic diversity on the Y chromosome, autosomes and mtDNA (ref. 1). We have found that Y chromosome variants tend to be more localized geographically than those of mtDNA and the autosomes. The fraction of variation within human populations for Y chromosome single nucleotide polymorphisms (SNPs) is 35.5%, versus 80–85% for the autosomes and mtDNA (refs 6, 7, 8). A higher female than male migration rate (via patrilocality, the tendency for a wife to move into her husbands natal household) explains most of this discrepancy, because diverse Y chromosomes would enter a population at a lower rate than mtDNA or the autosomes. Polygyny may also contribute, but the reduction of variation within populations that we measure for the Y chromosome, relative to the autosomes and mitochondrial DNA, is of such magnitude that differences in the effective population sizes of the sexes alone are insufficient to produce the observation.

Common variants at CD40 and other loci confer risk of rheumatoid arthritis

To identify rheumatoid arthritis risk loci in European populations, we conducted a meta-analysis of two published genome-wide association (GWA) studies totaling 3,393 cases and 12,462 controls. We genotyped 31 top-ranked SNPs not previously associated with rheumatoid arthritis in an independent replication of 3,929 autoantibody-positive rheumatoid arthritis cases and 5,807 matched controls from eight separate collections. We identified a common variant at the CD40 gene locus (rs4810485, P = 0.0032 replication, P = 8.2 × 10−9 overall, OR = 0.87). Along with other associations near TRAF1 (refs. 2,3) and TNFAIP3 (refs. 4,5), this implies a central role for the CD40 signaling pathway in rheumatoid arthritis pathogenesis. We also identified association at the CCL21 gene locus (rs2812378, P = 0.00097 replication, P = 2.8 × 10−7 overall), a gene involved in lymphocyte trafficking. Finally, we identified evidence of association at four additional gene loci: MMEL1-TNFRSF14 (rs3890745, P = 0.0035 replication, P = 1.1 × 10−7 overall), CDK6 (rs42041, P = 0.010 replication, P = 4.0 × 10−6 overall), PRKCQ (rs4750316, P = 0.0078 replication, P = 4.4 × 10−6 overall), and KIF5A-PIP4K2C (rs1678542, P = 0.0026 replication, P = 8.8 × 10−8 overall).

Genome-wide association identifies multiple ulcerative colitis susceptibility loci

Ulcerative colitis is a chronic, relapsing inflammatory condition of the gastrointestinal tract with a complex genetic and environmental etiology. In an effort to identify genetic variation underlying ulcerative colitis risk, we present two distinct genome-wide association studies of ulcerative colitis and their joint analysis with a previously published scan, comprising, in aggregate, 2,693 individuals with ulcerative colitis and 6,791 control subjects. Fifty-nine SNPs from 14 independent loci attained an association significance of P < 10−5. Seven of these loci exceeded genome-wide significance (P < 5 × 10−8). After testing an independent cohort of 2,009 cases of ulcerative colitis and 1,580 controls, we identified 13 loci that were significantly associated with ulcerative colitis (P < 5 × 10−8), including the immunoglobulin receptor gene FCGR2A, 5p15, 2p16 and ORMDL3 (orosomucoid1-like 3). We confirmed association with 14 previously identified ulcerative colitis susceptibility loci, and an analysis of acknowledged Crohns disease loci showed that roughly half of the known Crohns disease associations are shared with ulcerative colitis. These data implicate approximately 30 loci in ulcerative colitis, thereby providing insight into disease pathogenesis.

The Distribution of Human Genetic Diversity: A Comparison of Mitochondrial, Autosomal, and Y-Chromosome Data

We report a comparison of worldwide genetic variation among 255 individuals by using autosomal, mitochondrial, and Y-chromosome polymorphisms. Variation is assessed by use of 30 autosomal restriction-site polymorphisms (RSPs), 60 autosomal short-tandem-repeat polymorphisms (STRPs), 13 Alu-insertion polymorphisms and one LINE-1 element, 611 bp of mitochondrial control-region sequence, and 10 Y-chromosome polymorphisms. Analysis of these data reveals substantial congruity among this diverse array of genetic systems. With the exception of the autosomal RSPs, in which an ascertainment bias exists, all systems show greater gene diversity in Africans than in either Europeans or Asians. Africans also have the largest total number of alleles, as well as the largest number of unique alleles, for most systems. GST values are 11%-18% for the autosomal systems and are two to three times higher for the mtDNA sequence and Y-chromosome RSPs. This difference is expected because of the lower effective population size of mtDNA and Y chromosomes. A lower value is seen for Y-chromosome STRs, reflecting a relative lack of continental population structure, as a result of rapid mutation and genetic drift. Africa has higher GST values than does either Europe or Asia for all systems except the Y-chromosome STRs and Alus. All systems except the Y-chromosome STRs show less variation between populations within continents than between continents. These results are reassuring in their consistency and offer broad support for an African origin of modern human populations.

Meta-analysis of genome-wide association studies identifies eight new loci for type 2 diabetes in east Asians.

We conducted a three-stage genetic study to identify susceptibility loci for type 2 diabetes (T2D) in east Asian populations. We followed our stage 1 meta-analysis of eight T2D genome-wide association studies (6,952 cases with T2D and 11,865 controls) with a stage 2 in silico replication analysis (5,843 cases and 4,574 controls) and a stage 3 de novo replication analysis (12,284 cases and 13,172 controls). The combined analysis identified eight new T2D loci reaching genome-wide significance, which mapped in or near GLIS3, PEPD, FITM2-R3HDML-HNF4A, KCNK16, MAEA, GCC1-PAX4, PSMD6 and ZFAND3. GLIS3, which is involved in pancreatic beta cell development and insulin gene expression, is known for its association with fasting glucose levels. The evidence of an association with T2D for PEPD and HNF4A has been shown in previous studies. KCNK16 may regulate glucose-dependent insulin secretion in the pancreas. These findings, derived from an east Asian population, provide new perspectives on the etiology of T2D.

Mapping Human Genetic Diversity in Asia

Patterns of Early Migration In order to gain insight into various migrations that must have happened during movement of early humans into Asia and the subsequent populating of the largest continent on Earth, the HUGO Pan-Asian SNP Consortium (p. 1541) analyzed genetic variation in almost 2000 individuals representing 73 Asian and two non-Asian populations. The results suggest that there may have been a single major migration of people into Asia and a subsequent south-to-north migration across the continent. While most populations from the same linguistic group tend to cluster together in terms of relatedness, several do not, clustering instead with their geographic neighbors, suggesting either substantial recent mixing among the populations or language replacement. Furthermore, data from indigenous Taiwanese populations appear to be inconsistent with the idea of a Taiwan homeland for Austronesian populations. Genetic analyses of Asian peoples suggest that the continent was populated through a single migration event. Asia harbors substantial cultural and linguistic diversity, but the geographic structure of genetic variation across the continent remains enigmatic. Here we report a large-scale survey of autosomal variation from a broad geographic sample of Asian human populations. Our results show that genetic ancestry is strongly correlated with linguistic affiliations as well as geography. Most populations show relatedness within ethnic/linguistic groups, despite prevalent gene flow among populations. More than 90% of East Asian (EA) haplotypes could be found in either Southeast Asian (SEA) or Central-South Asian (CSA) populations and show clinal structure with haplotype diversity decreasing from south to north. Furthermore, 50% of EA haplotypes were found in SEA only and 5% were found in CSA only, indicating that SEA was a major geographic source of EA populations.

Meta-analysis of genome-wide association studies identifies common variants associated with blood pressure variation in East Asians

We conducted a meta-analysis of genome-wide association studies of systolic (SBP) and diastolic (DBP) blood pressure in 19,608 subjects of east Asian ancestry from the AGEN-BP consortium followed up with de novo genotyping (n = 10,518) and further replication (n = 20,247) in east Asian samples. We identified genome-wide significant (P < 5 × 10−8) associations with SBP or DBP, which included variants at four new loci (ST7L-CAPZA1, FIGN-GRB14, ENPEP and NPR3) and a newly discovered variant near TBX3. Among the five newly discovered variants, we obtained significant replication in the independent samples for all of these loci except NPR3. We also confirmed seven loci previously identified in populations of European descent. Moreover, at 12q24.13 near ALDH2, we observed strong association signals (P = 7.9 × 10−31 and P = 1.3 × 10−35 for SBP and DBP, respectively) with ethnic specificity. These findings provide new insights into blood pressure regulation and potential targets for intervention.