Ancestry informative marker sets for determining continental origin and admixture proportions in common populations in America

Replication of genetic loci for ages at menarche and menopause in the multi-ethnic Population Architecture using Genomics and Epidemiology (PAGE) study

STUDY QUESTION

Do genetic associations identified in genome-wide association studies (GWAS) of age at menarche (AM) and age at natural menopause (ANM) replicate in women of diverse race/ancestry from the Population Architecture using Genomics and Epidemiology (PAGE) Study?

SUMMARY ANSWER

We replicated GWAS reproductive trait single nucleotide polymorphisms (SNPs) in our European descent population and found that many SNPs were also associated with AM and ANM in populations of diverse ancestry.

WHAT IS KNOWN ALREADY

Menarche and menopause mark the reproductive lifespan in women and are important risk factors for chronic diseases including obesity, cardiovascular disease and cancer. Both events are believed to be influenced by environmental and genetic factors, and vary in populations differing by genetic ancestry and geography. Most genetic variants associated with these traits have been identified in GWAS of European-descent populations.

STUDY DESIGN, SIZE, DURATION

A total of 42 251 women of diverse ancestry from PAGE were included in cross-sectional analyses of AM and ANM.

MATERIALS, SETTING, METHODS

SNPs previously associated with ANM (n = 5 SNPs) and AM (n = 3 SNPs) in GWAS were genotyped in American Indians, African Americans, Asians, European Americans, Hispanics and Native Hawaiians. To test SNP associations with ANM or AM, we used linear regression models stratified by race/ethnicity and PAGE sub-study. Results were then combined in race-specific fixed effect meta-analyses for each outcome. For replication and generalization analyses, significance was defined at P < 0.01 for ANM analyses and P < 0.017 for AM analyses.

MAIN RESULTS AND THE ROLE OF CHANCE

We replicated findings for AM SNPs in the LIN28B locus and an intergenic region on 9q31 in European Americans. The LIN28B SNPs (rs314277 and rs314280) were also significantly associated with AM in Asians, but not in other race/ethnicity groups. Linkage disequilibrium (LD) patterns at this locus varied widely among the ancestral groups. With the exception of an intergenic SNP at 13q34, all ANM SNPs replicated in European Americans. Three were significantly associated with ANM in other race/ethnicity populations: rs2153157 (6p24.2/SYCP2L), rs365132 (5q35/UIMC1) and rs16991615 (20p12.3/MCM8). While rs1172822 (19q13/BRSK1) was not significant in the populations of non-European descent, effect sizes showed similar trends.

LIMITATIONS, REASONS FOR CAUTION

Lack of association for the GWAS SNPs in the non-European American groups may be due to differences in locus LD patterns between these groups and the European-descent populations included in the GWAS discovery studies; and in some cases, lower power may also contribute to non-significant findings.

WIDER IMPLICATIONS OF THE FINDINGS

The discovery of genetic variants associated with the reproductive traits provides an important opportunity to elucidate the biological mechanisms involved with normal variation and disorders of menarche and menopause. In this study we replicated most, but not all reported SNPs in European descent populations and examined the epidemiologic architecture of these early reported variants, describing their generalizability and effect size across differing ancestral populations. Such data will be increasingly important for prioritizing GWAS SNPs for follow-up in fine-mapping and resequencing studies, as well as in translational research.

Replication of genome wide association studies of alcohol dependence: support for association with variation in ADH1C

Genome-wide association studies (GWAS) have revealed many single nucleotide polymorphisms (SNPs) associated with complex traits. Although these studies frequently fail to identify statistically significant associations, the top association signals from GWAS may be enriched for true associations. We therefore investigated the association of alcohol dependence with 43 SNPs selected from association signals in the first two published GWAS of alcoholism. Our analysis of 808 alcohol-dependent cases and 1,248 controls provided evidence of association of alcohol dependence with SNP rs1614972 in the ADH1C gene (unadjusted p = 0.0017). Because the GWAS study that originally reported association of alcohol dependence with this SNP [1] included only men, we also performed analyses in sex-specific strata. The results suggest that this SNP has a similar effect in both sexes (men: OR (95%CI) = 0.80 (0.66, 0.95); women: OR (95%CI) = 0.83 (0.66, 1.03)). We also observed marginal evidence of association of the rs1614972 minor allele with lower alcohol consumption in the non-alcoholic controls (p = 0.081), and independently in the alcohol-dependent cases (p = 0.046). Despite a number of potential differences between the samples investigated by the prior GWAS and the current study, data presented here provide additional support for the association of SNP rs1614972 in ADH1C with alcohol dependence and extend this finding by demonstrating association with consumption levels in both non-alcoholic and alcohol-dependent populations. Further studies should investigate the association of other polymorphisms in this gene with alcohol dependence and related alcohol-use phenotypes.

From single-SNP to wide-locus: genome-wide association studies identifying functionally related genes and intragenic regions in small sample studies

BACKGROUND

Genome-wide association studies (GWAS) have had limited success when applied to complex diseases. Analyzing SNPs individually requires several large studies to integrate the often divergent results. In the presence of epistasis, multivariate approaches based on the linear model (including stepwise logistic regression) often have low sensitivity and generate an abundance of artifacts.

METHODS

Recent advances in distributed and parallel processing spurred methodological advances in nonparametric statistics. U-statistics for structured multivariate data (µStat) are not confounded by unrealistic assumptions (e.g., linearity, independence).

RESULTS

By incorporating knowledge about relationships between SNPs, µGWAS (GWAS based on µStat) can identify clusters of genes around biologically relevant pathways and pinpoint functionally relevant regions within these genes.

CONCLUSION

With this computational biostatistics approach increasing power and guarding against artifacts, personalized medicine and comparative effectiveness will advance while subgroup analyses of Phase III trials can now suggest risk factors for adverse events and novel directions for drug development.

ETHNOPRED: a novel machine learning method for accurate continental and sub-continental ancestry identification and population stratification correction

BACKGROUND

Population stratification is a systematic difference in allele frequencies between subpopulations. This can lead to spurious association findings in the case-control genome wide association studies (GWASs) used to identify single nucleotide polymorphisms (SNPs) associated with disease-linked phenotypes. Methods such as self-declared ancestry, ancestry informative markers, genomic control, structured association, and principal component analysis are used to assess and correct population stratification but each has limitations. We provide an alternative technique to address population stratification.

RESULTS

We propose a novel machine learning method, ETHNOPRED, which uses the genotype and ethnicity data from the HapMap project to learn ensembles of disjoint decision trees, capable of accurately predicting an individual's continental and sub-continental ancestry. To predict an individual's continental ancestry, ETHNOPRED produced an ensemble of 3 decision trees involving a total of 10 SNPs, with 10-fold cross validation accuracy of 100% using HapMap II dataset. We extended this model to involve 29 disjoint decision trees over 149 SNPs, and showed that this ensemble has an accuracy of ≥ 99.9%, even if some of those 149 SNP values were missing. On an independent dataset, predominantly of Caucasian origin, our continental classifier showed 96.8% accuracy and improved genomic control's λ from 1.22 to 1.11. We next used the HapMap III dataset to learn classifiers to distinguish European subpopulations (North-Western vs. Southern), East Asian subpopulations (Chinese vs. Japanese), African subpopulations (Eastern vs. Western), North American subpopulations (European vs. Chinese vs. African vs. Mexican vs. Indian), and Kenyan subpopulations (Luhya vs. Maasai). In these cases, ETHNOPRED produced ensembles of 3, 39, 21, 11, and 25 disjoint decision trees, respectively involving 31, 502, 526, 242 and 271 SNPs, with 10-fold cross validation accuracy of 86.5% ± 2.4%, 95.6% ± 3.9%, 95.6% ± 2.1%, 98.3% ± 2.0%, and 95.9% ± 1.5%. However, ETHNOPRED was unable to produce a classifier that can accurately distinguish Chinese in Beijing vs. Chinese in Denver.

CONCLUSIONS

ETHNOPRED is a novel technique for producing classifiers that can identify an individual's continental and sub-continental heritage, based on a small number of SNPs. We show that its learned classifiers are simple, cost-efficient, accurate, transparent, flexible, fast, applicable to large scale GWASs, and robust to missing values.

Meta-analysis of genetic polymorphisms in granulomatosis with polyangiitis (Wegener's) reveals shared susceptibility loci with rheumatoid arthritis

OBJECTIVE

To examine the association of previously identified autoimmune disease susceptibility loci with granulomatosis with polyangiitis (Wegener's) (GPA), and to determine whether the genetic susceptibility profiles of other autoimmune diseases are associated with those of GPA.

METHODS

Genetic data from 2 cohorts were meta-analyzed. Genotypes for 168 previously identified single-nucleotide polymorphisms (SNPs) associated with susceptibility to different autoimmune diseases were ascertained in a total of 880 patients with GPA and 1,969 control subjects of European descent. Single-marker associations were identified using additive logistic regression models. Associations of multiple SNPs with GPA were assessed using genetic risk scores based on susceptibility loci for Crohn's disease, type 1 diabetes, systemic lupus erythematosus, rheumatoid arthritis (RA), celiac disease, and ulcerative colitis. Adjustment for population substructure was performed in all analyses, using ancestry-informative markers and principal components analysis.

RESULTS

Genetic polymorphisms in CTLA4 were significantly associated with GPA in the single-marker meta-analysis (odds ratio [OR] 0.79, 95% confidence interval [95% CI] 0.70-0.89, P = 9.8 × 10(-5) ). The genetic risk score for RA susceptibility markers was significantly associated with GPA (OR 1.05 per 1-unit increase in genetic risk score, 95% CI 1.02-1.08, P = 5.1 × 10(-5) ).

CONCLUSION

RA and GPA may arise from a similar genetic predisposition. Aside from CTLA4, other loci previously found to be associated with common autoimmune diseases were not statistically significantly associated with GPA in this study.

Impact of genetic ancestry and sociodemographic status on the clinical expression of systemic lupus erythematosus in American Indian-European populations

OBJECTIVE

American Indian-Europeans, Asians, and African Americans have an excess morbidity from systemic lupus erythematosus (SLE) and a higher prevalence of lupus nephritis than do Caucasians. The aim of this study was to analyze the relationship between genetic ancestry and sociodemographic characteristics and clinical features in a large cohort of American Indian-European SLE patients.

METHODS

A total of 2,116 SLE patients of American Indian-European origin and 4,001 SLE patients of European descent for whom we had clinical data were included in the study. Genotyping of 253 continental ancestry-informative markers was performed on the Illumina platform. Structure and Admixture software were used to determine genetic ancestry proportions of each individual. Logistic regression was used to test the association between genetic ancestry and sociodemographic and clinical characteristics. Odds ratios (ORs) were calculated with 95% confidence intervals (95% CIs).

RESULTS

The average American Indian genetic ancestry of 2,116 SLE patients was 40.7%. American Indian genetic ancestry conferred increased risks of renal involvement (P < 0.0001, OR 3.50 [95% CI 2.63- 4.63]) and early age at onset (P < 0.0001). American Indian ancestry protected against photosensitivity (P < 0.0001, OR 0.58 [95% CI 0.44-0.76]), oral ulcers (P < 0.0001, OR 0.55 [95% CI 0.42-0.72]), and serositis (P < 0.0001, OR 0.56 [95% CI 0.41-0.75]) after adjustment for age, sex, and age at onset. However, age and sex had stronger effects than genetic ancestry on malar rash, discoid rash, arthritis, and neurologic involvement.

CONCLUSION

In general, American Indian genetic ancestry correlates with lower sociodemographic status and increases the risk of developing renal involvement and SLE at an earlier age.

Amerindian genetic ancestry protects against Alzheimer's disease

BACKGROUND

Alzheimer's disease (AD) is the most common form of dementia worldwide, and bears remarkable evidence for a differential prevalence among continental populations. In this scenario, estimating ancestry proportions in recently admixed populations is a strategy that can help increasing knowledge about the genetic structure of this complex trait.

AIM/METHODS

Our purpose was to assess mean ancestry estimates for the three main parental contributors to the Brazilian contingent (European, African and Amerindian) using a panel of 12 ancestry informative markers. Outpatients with the late-onset form of AD (n = 120) were compared for ancestry levels with non-cognitively impaired subjects (n = 412) in the Midwest Brazil, controlling for classic clinical, social and anthropometric risk factors.

RESULTS

Our findings show a 3-fold greater genetic Amerindian content among control subjects compared to AD patients (p < 0.001).

CONCLUSION

Our results suggest that the allelic architecture of Native Americans can confer protection against the onset of the disease.

Ancestral effect on HOMA-IR levels quantitated in an American population of Mexican origin

OBJECTIVE

An elevated insulin resistance index (homeostasis model assessment of insulin resistance [HOMA-IR]) is more commonly seen in the Mexican American population than in European populations. We report quantitative ancestral effects within a Mexican American population, and we correlate ancestral components with HOMA-IR.

RESEARCH DESIGN AND METHODS

We performed ancestral analysis in 1,551 participants of the Cameron County Hispanic Cohort by genotyping 103 ancestry-informative markers (AIMs). These AIMs allow determination of the percentage (0-100%) ancestry from three major continental populations, i.e., European, African, and Amerindian.

RESULTS

We observed that predominantly Amerindian ancestral components were associated with increased HOMA-IR (β = 0.124, P = 1.64 × 10(-7)). The correlation was more significant in males (Amerindian β = 0.165, P = 5.08 × 10(-7)) than in females (Amerindian β = 0.079, P = 0.019).

CONCLUSIONS

This unique study design demonstrates how genomic markers for quantitative ancestral information can be used in admixed populations to predict phenotypic traits such as insulin resistance.

Human leukocyte antigen profiles of latin american populations: differential admixture and its potential impact on hematopoietic stem cell transplantation

The outcome of hematopoietic stem cell transplantation (HSCT) is shaped by both clinical and genetic factors that determine its success. Genetic factors including human leukocyte antigen (HLA) and non-HLA genetic variants are believed to influence the risk of potentially fatal complications after the transplant. Moreover, ethnicity has been proposed as a factor modifying the risk of graft-versus-host disease. The populations of Latin America are a complex array of different admixture processes with varying degrees of ancestral population proportions that came in different migration waves. This complexity makes the study of genetic risks in this region complicated unless the extent of this variation is thoroughly characterized. In this study we compared the HLA-A and HLA-B allele group profiles for 31 Latin American populations and 61 ancestral populations from Iberia, Italy, Sub-Saharan Africa, and America. Results from population genetics comparisons show a wide variation in the HLA profiles from the Latin American populations that correlate with different admixture proportions. Populations in Latin America seem to be organized in at least three groups with (1) strong Amerindian admixture, (2) strong Caucasian component, and (3) a Caucasian-African gradient. These results imply that genetic risk assessment for HSCT in Latin America has to be adapted for different population subgroups rather than as a pan-Hispanic/Latino analysis.

First all-in-one diagnostic tool for DNA intelligence: genome-wide inference of biogeographic ancestry, appearance, relatedness, and sex with the Identitas v1 Forensic Chip

When a forensic DNA sample cannot be associated directly with a previously genotyped reference sample by standard short tandem repeat profiling, the investigation required for identifying perpetrators, victims, or missing persons can be both costly and time consuming. Here, we describe the outcome of a collaborative study using the Identitas Version 1 (v1) Forensic Chip, the first commercially available all-in-one tool dedicated to the concept of developing intelligence leads based on DNA. The chip allows parallel interrogation of 201,173 genome-wide autosomal, X-chromosomal, Y-chromosomal, and mitochondrial single nucleotide polymorphisms for inference of biogeographic ancestry, appearance, relatedness, and sex. The first assessment of the chip’s performance was carried out on 3,196 blinded DNA samples of varying quantities and qualities, covering a wide range of biogeographic origin and eye/hair coloration as well as variation in relatedness and sex. Overall, 95 % of the samples (N = 3,034) passed quality checks with an overall genotype call rate >90 % on variable numbers of available recorded trait information. Predictions of sex, direct match, and first to third degree relatedness were highly accurate. Chip-based predictions of biparental continental ancestry were on average ~94 % correct (further support provided by separately inferred patrilineal and matrilineal ancestry). Predictions of eye color were 85 % correct for brown and 70 % correct for blue eyes, and predictions of hair color were 72 % for brown, 63 % for blond, 58 % for black, and 48 % for red hair. From the 5 % of samples (N = 162) with <90 % call rate, 56 % yielded correct continental ancestry predictions while 7 % yielded sufficient genotypes to allow hair and eye color prediction. Our results demonstrate that the Identitas v1 Forensic Chip holds great promise for a wide range of applications including criminal investigations, missing person investigations, and for national security purposes.

Underlying genetic structure impacts the association between CYP2B6 polymorphisms and response to efavirenz and nevirapine

OBJECTIVE

CYP2B6 variation predicts pharmacokinetic characteristics of its substrates. Consideration for underlying genetic structure is critical to protect against spurious associations with the highly polymorphic CYP2B6 gene.

DESIGN

The effect of CYP2B6 variation on response to its substrates, nonnucleoside reverse transcriptase inhibitors (NNRTIs), was explored in the Women's Interagency HIV Study.

METHODS

Five putative functional polymorphisms were tested for associations with virologic suppression within 1 year after NNRTI initiation in women naive to antiretroviral agents (n = 91). Principal components were generated to control for population substructure. Logistic regression was used to test the joint effect of rs3745274 and rs28399499, which together indicate slow, intermediate, and extensive metabolizers.

RESULTS

Rs3745274 was significantly associated with virologic suppression [odds ratio = 3.61, 95% confidence interval (CI) 1.16-11.22, P trend = 0.03]; the remaining polymorphisms tested were not significantly associated with response. Women classified as intermediate and slow metabolizers were 2.90 (95% CI 0.79-12.28) and 13.44 (95% CI 1.66 to infinity) times as likely to achieve virologic suppression compared to extensive metabolizers after adjustment for principal components (P trend = 0.005). Failure to control for genetic ancestry resulted in substantial confounding of the relationship between the metabolizer phenotype and treatment response.

CONCLUSION

The CYP2B6 metabolizer phenotype was significantly associated with virologic response to NNRTIs; this relationship would have been masked by simple adjustment for self-reported ethnicity. Given the appreciable genetic heterogeneity that exists within self-reported ethnicity, these results exemplify the importance of characterizing underlying genetic structure in pharmacogenetic studies. Further follow-up of the CYP2B6 metabolizer phenotype is warranted, given the potential clinical importance of this finding.

Risk for myasthenia gravis maps to a (151) Pro→Ala change in TNIP1 and to human leukocyte antigen-B*08

OBJECTIVE

The objective of this study is to comprehensively define the genetic basis of early onset myasthenia gravis (EOMG).

METHODS

We have carried out a 2-stage genome-wide association study on a total of 649 North European EOMG patients. Cases were matched 1:4 with controls of European ancestry. We performed imputation and conditional analyses across the major histocompatibility complex, as well as in the top regions of association outside the human leukocyte antigen (HLA) region.

RESULTS

We observed the strongest association in the HLA class I region at rs7750641 (p = 1.2 × 10(-92) ; odds ratio [OR], 6.25). By imputation and conditional analyses, HLA-B*08 proves to be the major associated allele (p = 2.87 × 10(-113) ; OR, 6.41). In addition to the expected association with PTPN22 (rs2476601; OR, 1.71; p = 8.2 × 10(-10) ), an imputed coding variant (rs2233290) at position 151 (Pro→Ala) in the TNFAIP3-interacting protein 1, TNIP1, confers even stronger risk than PTPN22 (OR, 1.91; p = 3.2 × 10(-10) ).

INTERPRETATION

The association at TNIP1 in EOMG implies disease mechanisms involving ubiquitin-dependent dysregulation of NF-κB signaling. The localization of the major HLA signal to the HLA-B*08 allele suggests that CD8(+) T cells may play a key role in disease initiation or pathogenesis.

European genetic ancestry is associated with a decreased risk of lupus nephritis

OBJECTIVE

African Americans, East Asians, and Hispanics with systemic lupus erythematosus (SLE) are more likely to develop renal disease than are SLE patients of European descent. This study was undertaken to investigate whether European genetic ancestry protects against the development of lupus nephritis, with the aim of exploring the genetic and socioeconomic factors that might explain this effect.

METHODS

This was a cross-sectional study of SLE patients from a multiethnic case collection. Participants were genotyped for 126 single-nucleotide polymorphisms (SNPs) informative for ancestry. A subset of participants was also genotyped for 80 SNPs in 14 candidate genes for renal disease in SLE. Logistic regression was used to test the association between European ancestry and renal disease. Analyses were adjusted for continental ancestries, socioeconomic status (SES), and candidate genes.

RESULTS

Participants (n = 1,906) had, on average, 62.4% European, 15.8% African, 11.5% East Asian, 6.5% Amerindian, and 3.8% South Asian ancestry. Among the participants, 656 (34%) had renal disease. A 10% increase in the proportion of European ancestry estimated in each participant was associated with a 15% reduction in the odds of having renal disease, after adjustment for disease duration and sex (odds ratio 0.85, 95% confidence interval 0.82-0.87; P = 1.9 × 10(-30) ). Adjustment for other genetic ancestries, measures of SES, or SNPs in the genes most associated with renal disease (IRF5 [rs4728142], BLK [rs2736340], STAT4 [rs3024912], and HLA-DRB1*0301 and DRB1*1501) did not substantively alter this relationship.

CONCLUSION

European ancestry is protective against the development of renal disease in SLE, an effect that is independent of other genetic ancestries, candidate risk alleles, and socioeconomic factors.

Introducing the Forensic Research/Reference on Genetics knowledge base, FROG-kb

BACKGROUND

Online tools and databases based on multi-allelic short tandem repeat polymorphisms (STRPs) are actively used in forensic teaching, research, and investigations. The Fst value of each CODIS marker tends to be low across the populations of the world and most populations typically have all the common STRP alleles present diminishing the ability of these systems to discriminate ethnicity. Recently, considerable research is being conducted on single nucleotide polymorphisms (SNPs) to be considered for human identification and description. However, online tools and databases that can be used for forensic research and investigation are limited.

METHODS

The back end DBMS (Database Management System) for FROG-kb is Oracle version 10. The front end is implemented with specific code using technologies such as Java, Java Servlet, JSP, JQuery, and GoogleCharts.

RESULTS

We present an open access web application, FROG-kb (Forensic Research/Reference on Genetics-knowledge base, http://frog.med.yale.edu), that is useful for teaching and research relevant to forensics and can serve as a tool facilitating forensic practice. The underlying data for FROG-kb are provided by the already extensively used and referenced ALlele FREquency Database, ALFRED (http://alfred.med.yale.edu). In addition to displaying data in an organized manner, computational tools that use the underlying allele frequencies with user-provided data are implemented in FROG-kb. These tools are organized by the different published SNP/marker panels available. This web tool currently has implemented general functions possible for two types of SNP panels, individual identification and ancestry inference, and a prediction function specific to a phenotype informative panel for eye color.

CONCLUSION

The current online version of FROG-kb already provides new and useful functionality. We expect FROG-kb to grow and expand in capabilities and welcome input from the forensic community in identifying datasets and functionalities that will be most helpful and useful. Thus, the structure and functionality of FROG-kb will be revised in an ongoing process of improvement. This paper describes the state as of early June 2012.

PNPLA3 polymorphisms and liver aminotransferase levels in a Mexican American population

PURPOSE

This study examined genetic associations of patatin-like phospholipase domain containing 3 gene (PNPLA3) polymorphisms and liver aminotransferases in an extensively documented, randomly recruited Mexican American population at high risk of liver disease.

METHODS

Two single nucleotide polymorphisms (SNP) in the PNPLA3 gene (i.e., rs738409 and rs2281135) were genotyped in 1532 individuals. Population stratification was corrected by the genotyping of 103 ancestry informative markers (AIMs) for Mexican Americans.

RESULTS

Both PNPLA3 SNPs showed highly significant association with alanine aminotransferase (ALT) levels, but was also, in males, associated with aspartate aminotransferase (AST) levels. Haplotypic association test of the two SNPs suggested stronger genetic association with rs738409 than rs2281135. Obvious sex effects were observed: rs738409-sex interaction in ALT levels P = 8.37 x 10(-4); rs738409-sex interaction in AST levels P = 5.03 x 10(-3).

CONCLUSIONS

This population study highlights a sex-specific association of PNPLA3 polymorphisms and elevated liver enzymes in a population-based study, independent of common pathological factors of the metabolic syndrome. The strong genetic association found in women ≤ 50 years old, but not in women > 50 years old, suggests that sex hormones may mediate the sex effect.

Integrative analysis of single nucleotide polymorphisms and gene expression efficiently distinguishes samples from closely related ethnic populations

Background

Ancestry informative markers (AIMs) are a type of genetic marker that is informative for tracing the ancestral ethnicity of individuals. Application of AIMs has gained substantial attention in population genetics, forensic sciences, and medical genetics. Single nucleotide polymorphisms (SNPs), the materials of AIMs, are useful for classifying individuals from distinct continental origins but cannot discriminate individuals with subtle genetic differences from closely related ancestral lineages. Proof-of-principle studies have shown that gene expression (GE) also is a heritable human variation that exhibits differential intensity distributions among ethnic groups. GE supplies ethnic information supplemental to SNPs; this motivated us to integrate SNP and GE markers to construct AIM panels with a reduced number of required markers and provide high accuracy in ancestry inference. Few studies in the literature have considered GE in this aspect, and none have integrated SNP and GE markers to aid classification of samples from closely related ethnic populations.

Results

We integrated a forward variable selection procedure into flexible discriminant analysis to identify key SNP and/or GE markers with the highest cross-validation prediction accuracy. By analyzing genome-wide SNP and/or GE markers in 210 independent samples from four ethnic groups in the HapMap II Project, we found that average testing accuracies for a majority of classification analyses were quite high, except for SNP-only analyses that were performed to discern study samples containing individuals from two close Asian populations. The average testing accuracies ranged from 0.53 to 0.79 for SNP-only analyses and increased to around 0.90 when GE markers were integrated together with SNP markers for the classification of samples from closely related Asian populations. Compared to GE-only analyses, integrative analyses of SNP and GE markers showed comparable testing accuracies and a reduced number of selected markers in AIM panels.

Conclusions

Integrative analysis of SNP and GE markers provides high-accuracy and/or cost-effective classification results for assigning samples from closely related or distantly related ancestral lineages to their original ancestral populations. User-friendly BIASLESS (Biomarkers Identification and Samples Subdivision) software was developed as an efficient tool for selecting key SNP and/or GE markers and then building models for sample subdivision. BIASLESS was programmed in R and R-GUI and is available online at http://www.stat.sinica.edu.tw/hsinchou/genetics/prediction/BIASLESS.htm.

Risk variants for psoriasis vulgaris in a large case-control collection and association with clinical subphenotypes

Recent genome-wide association studies (GWASs) have identified >20 new loci associated with the susceptibility to psoriasis vulgaris (PsV) risk. We investigated the association of PsV and its main clinical subphenotypes with 32 loci having previous genome-wide evidence of association with PsV (P < 5e-8) or strong GWAS evidence (P < 5e-5 in discovery and P < 0.05 in replication sample) in a large cohort of PsV patients (n = 2005) and controls (n = 1497). We provide the first independent replication for COG6 (P = 0.00079) and SERPINB8 (P = 0.048) loci with PsV. In those patients having developed psoriatic arthritis (n = 955), we found, for the first time, a strong association with IFIH1 (P = 0.013). Analyses of clinically relevant PsV subtypes yielded a significant association of severity of cutaneous disease with variation at LCE3D locus (P = 0.0005) in PsV and nail involvement with IL1RN in purely cutaneous psoriasis (PsC, P = 0.007). In an exploratory analysis of epistasis, we replicated the previously described HLA-C-ERAP1 interaction with PsC. Our findings show that common genetic variants associated with a complex phenotype like PsV influence different subphenotypes of high clinical relevance.

Genetic ancestry-smoking interactions and lung function in African Americans: a cohort study

BACKGROUND

Smoking tobacco reduces lung function. African Americans have both lower lung function and decreased metabolism of tobacco smoke compared to European Americans. African ancestry is also associated with lower pulmonary function in African Americans. We aimed to determine whether African ancestry modifies the association between smoking and lung function and its rate of decline in African Americans.

METHODOLOGY/PRINCIPAL FINDINGS

We evaluated a prospective ongoing cohort of 1,281 African Americans participating in the Health, Aging, and Body Composition (Health ABC) Study initiated in 1997. We also examined an ongoing prospective cohort initiated in 1985 of 1,223 African Americans in the Coronary Artery Disease in Young Adults (CARDIA) Study. Pulmonary function and tobacco smoking exposure were measured at baseline and repeatedly over the follow-up period. Individual genetic ancestry proportions were estimated using ancestry informative markers selected to distinguish European and West African ancestry. African Americans with a high proportion of African ancestry had lower baseline forced expiratory volume in one second (FEV₁) per pack-year of smoking (-5.7 ml FEV₁/ smoking pack-year) compared with smokers with lower African ancestry (-4.6 ml in FEV₁/ smoking pack-year) (interaction P value  = 0.17). Longitudinal analyses revealed a suggestive interaction between smoking, and African ancestry on the rate of FEV(1) decline in Health ABC and independently replicated in CARDIA.

CONCLUSIONS/SIGNIFICANCE

African American individuals with a high proportion of African ancestry are at greater risk for losing lung function while smoking.

Genetic risk factors for thrombosis in systemic lupus erythematosus

OBJECTIVE

Thrombosis is a serious complication of systemic lupus erythematosus (SLE). We investigated whether genetic variants implicated in thrombosis pathways are associated with thrombosis among 2 ethnically diverse SLE cohorts.

METHODS

Our discovery cohort consisted of 1698 patients with SLE enrolled in the University of California, San Francisco, Lupus Genetics Project and our replication cohort included 1361 patients with SLE enrolled in the PROFILE cohort. Patients fulfilled American College of Rheumatology SLE criteria, and data relevant to thrombosis were available. Thirty-three single nucleotide polymorphisms (SNP) previously shown to be associated with risk of deep venous thrombosis in the general population or implicated in thrombosis pathways were genotyped and tested for association with thrombosis in bivariate allelic analyses. SNP with p < 0.1 in the bivariate analyses were further tested in multivariable logistic regression models adjusted for age, sex, disease duration, antiphospholipid antibody status, smoking, nephritis, and medications.

RESULTS

In the discovery cohort, 23% of patients with SLE experienced a thrombotic event. SNP in the following genes demonstrated association with thrombosis risk overall in the discovery or replication cohorts and were assessed using metaanalytic methods: factor V Leiden (FVL) rs6025 (OR 1.85, p = 0.02) and methylenetetrahydrofolate reductase (MTHFR) rs1801133 (OR 0.75, p = 0.04) in whites, and fibrinogen gamma (FGG) rs2066865 (OR 1.91, p = 0.01) in Hispanic Americans. SNP in these genes showed association with venous thrombosis risk in whites: MTHFR rs1801131 (OR 1.51, p = 0.01), MTHFR rs1801133 (OR 0.70, p = 0.04), FVL rs6025 (OR 2.69, p = 0.002), and FGG rs2066865 (OR 1.49, p = 0.02) in whites. A SNP in FGG rs2066865 (OR 2.19, p = 0.003) demonstrated association with arterial thrombosis risk in Hispanics.

CONCLUSION

Our results implicate specific genetic risk factors for thrombosis in patients with SLE and suggest that genetic risk for thrombosis differs across ethnic groups.

Design, conduct, and analysis of a multicenter, pharmacogenomic, biomarker study in matched patients with severe sepsis treated with or without drotrecogin Alfa (activated)

BACKGROUND

A genomic biomarker identifying patients likely to benefit from drotrecogin alfa (activated) (DAA) may be clinically useful as a companion diagnostic. This trial was designed to validate biomarkers (improved response polymorphisms (IRPs)). Each IRP (A and B) contains two single nucleotide polymorphisms that were associated with a differential DAA treatment effect.

METHODS

DAA is typically given to younger patients with greater disease severity; therefore, a well-matched control group is critical to this multicenter, retrospective, controlled, outcome-blinded, genotype-blinded trial. Within each center, DAA-treated patients will be matched to controls treated within 24 months of each other taking into account age, APACHE II, cardiovascular, respiratory, renal, and hematologic dysfunction, mechanical ventilation status, medical/surgical status, and infection site. A propensity score will estimate the probability that a patient would have received DAA given their baseline characteristics. Two-phase data transfer will ensure unbiased selection of matched controls. The first transfer will be for eligibility and matching data and the second transfer for outcomes and genotypic data. The primary analysis will compare the effect of DAA in IRP + and IRP - groups on in-hospital mortality through day 28.

DISCUSSION

A design-based approach matching DAA-free to DAA-treated patients in a multicenter study of patients who have severe sepsis and high risk of death will directly compare control to DAA-treated groups for mortality by genotype. Results, which should be available in 2012, may help to identify the group of patients who would benefit from DAA and may provide a model for future investigation of sepsis therapies.

Relationship between gallbladder surgery and ethnic admixture in African American and Hispanic American women

OBJECTIVES

The objective of this study was to investigate whether differences in admixture in African American and Hispanic American adult women are associated with differences in gallbladder surgery.

METHODS

Gallbladder surgery history on entry to the Women's Health Initiative's (WHI) study was used as a dichotomous outcome measure for this study. The proportion of European, sub-Saharan African, and Amerindian (AMI) admixture was estimated for 10,841 African American and 4,620 Hispanic American women in WHI using 92 ancestry informative markers. Logistic regression analyses assessed the relationship between admixture and gallbladder surgery in WHI women (enrollment at ages >50, mean age 61 years) with or without adjusting for multiple covariates, including measures of adiposity, parity, alcohol use, and education.

RESULTS

There was a significant positive association between AMI admixture and the frequency of gallbladder surgery in Hispanic Americans. The odds ratio (OR) and 95% confidence intervals (CIs) for AMI admixture group was OR=2.97, CI=2.01-4.38, P<10(-4). Although there were strong positive associations with parity and adiposity, and negative associations with alcohol consumption and education, accounting for these covariates did not remove the admixture association (OR=2.46, CI=1.62-3.73). In contrast, the effect of African admixture was nearly indistinguishable from that of the European admixture, both of which were protective in the Hispanic American group, and African admixture had a marginal association with decreased gallbladder surgery in the African American group. Measures of adiposity were associated with increased risk for gallbladder surgery and remained significant after accounting for admixture and each of the other covariates. Education level and alcohol use were associated with decreased risk for gallbladder disease. However, after accounting for the other covariates these variably remained significant.

CONCLUSIONS

AMI admixture is strongly associated with gallbladder surgery in women, even after adjustment for selected risk factors for cholelithiasis. Additional studies to ascertain the specific genetic risk factors underlying these associations are warranted.

Consistent directions of effect for established type 2 diabetes risk variants across populations: the population architecture using Genomics and Epidemiology (PAGE) Consortium

Common genetic risk variants for type 2 diabetes (T2D) have primarily been identified in populations of European and Asian ancestry. We tested whether the direction of association with 20 T2D risk variants generalizes across six major racial/ethnic groups in the U.S. as part of the Population Architecture using Genomics and Epidemiology Consortium (16,235 diabetes case and 46,122 control subjects of European American, African American, Hispanic, East Asian, American Indian, and Native Hawaiian ancestry). The percentage of positive (odds ratio [OR] >1 for putative risk allele) associations ranged from 69% in American Indians to 100% in European Americans. Of the nine variants where we observed significant heterogeneity of effect by racial/ethnic group (P(heterogeneity) < 0.05), eight were positively associated with risk (OR >1) in at least five groups. The marked directional consistency of association observed for most genetic variants across populations implies a shared functional common variant in each region. Fine-mapping of all loci will be required to reveal markers of risk that are important within and across populations.

Examination of ancestral informative markers and self-reported race with tumor characteristics of breast cancer among Black and White women

African American (AA) women have a higher mortality from breast cancer (BC) compared to European American (EA) women. This may be due to the higher proportion of AA women with tumors that are diagnosed at more advanced stages and are characterized as being estrogen receptor negative (ER-)/progesterone receptor negative (PR-). Our study sought to determine whether self-reported race and percent African ancestry were associated with BC tumor characteristics. In a multi-center, population-based case-control study of BC, we determined percent African ancestry using ancestry informative markers (AIM) among women self-reporting race as AA or Black. BC tumor characteristics were associated with self-reported race (including a 30 % reduction in ER+/PR+ tumors [95 % confidence interval [CI]: 0.6-0.9] and a 1.5-fold increased risk of high grade [95 % CI: 1.2-1.9] for AA women compared to EA women). AIMs among AA women were not associated with BC tumor characteristics (AA women with ≥95 % versus <80 % African ancestry, odds ratio [OR] = 1.0 for ER+/PR+ [95 % CI: 0.6-1.8] and OR = 0.9 for high-grade tumors [95 % CI: 0.6-1.4]). Similar findings were observed for BC stage. While BC subtypes were associated with self-reported race, BC subtypes were not associated with percent African ancestry. These study results suggest that subtle differences in percent African ancestry are less important than the overall presence of African ancestry in relation to BC tumor characteristics.

Genetic analysis of ancestry, admixture and selection in Bolivian and Totonac populations of the New World

BACKGROUND

Populations of the Americas were founded by early migrants from Asia, and some have experienced recent genetic admixture. To better characterize the native and non-native ancestry components in populations from the Americas, we analyzed 815,377 autosomal SNPs, mitochondrial hypervariable segments I and II, and 36 Y-chromosome STRs from 24 Mesoamerican Totonacs and 23 South American Bolivians.

RESULTS AND CONCLUSIONS

We analyzed common genomic regions from native Bolivian and Totonac populations to identify 324 highly predictive Native American ancestry informative markers (AIMs). As few as 40-50 of these AIMs perform nearly as well as large panels of random genome-wide SNPs for predicting and estimating Native American ancestry and admixture levels. These AIMs have greater New World vs. Old World specificity than previous AIMs sets. We identify highly-divergent New World SNPs that coincide with high-frequency haplotypes found at similar frequencies in all populations examined, including the HGDP Pima, Maya, Colombian, Karitiana, and Surui American populations. Some of these regions are potential candidates for positive selection. European admixture in the Bolivian sample is approximately 12%, though individual estimates range from 0-48%. We estimate that the admixture occurred ~360-384 years ago. Little evidence of European or African admixture was found in Totonac individuals. Bolivians with pre-Columbian mtDNA and Y-chromosome haplogroups had 5-30% autosomal European ancestry, demonstrating the limitations of Y-chromosome and mtDNA haplogroups and the need for autosomal ancestry informative markers for assessing ancestry in admixed populations.

Classical HLA-DRB1 and DPB1 alleles account for HLA associations with primary biliary cirrhosis

Susceptibility to primary biliary cirrhosis (PBC) is strongly associated with HLA region polymorphisms. To determine if associations can be explained by classical HLA determinants we studied Italian 676 cases and 1440 controls with genotyped with dense single nucleotide polymorphisms (SNPs) for which classical HLA alleles and amino acids were imputed. Although previous genome-wide association studies and our results show stronger SNP associations near DQB1, we demonstrate that the HLA signals can be attributed to classical DRB1 and DPB1 genes. Strong support for the predominant role of DRB1 is provided by our conditional analyses. We also demonstrate an independent association of DPB1. Specific HLA-DRB1 genes (*08, *11 and *14) account for most of the DRB1 association signal. Consistent with previous studies, DRB1*08 (p = 1.59 × 10⁻¹¹) was the strongest predisposing allele where as DRB1*11 (p = 1.42 × 10⁻¹⁰) was protective. Additionally DRB1*14 and the DPB1 association (DPB1*03:01) (p = 9.18 × 10⁻⁷) were predisposing risk alleles. No signal was observed in the HLA class 1 or class 3 regions. These findings better define the association of PBC with HLA and specifically support the role of classical HLA-DRB1 and DPB1 genes and alleles in susceptibility to PBC.

Relationship between diabetes risk and admixture in postmenopausal African-American and Hispanic-American women

AIMS/HYPOTHESIS

Type 2 diabetes is more prevalent in African-Americans (AFAs) and Hispanic-Americans (HAs) than in European-Americans. We assessed whether continental admixture was correlated with diabetes risk in these high-risk groups.

METHODS

We estimated the proportion of sub-Saharan African (AFR), Amerindian (AMI) and European admixture using 92 ancestry-informative marker genotypes in 16,476 AFA and HA women from the Women's Health Initiative. Cox regression models were used to examine the association between admixture and diabetes risk, with and without accounting for socioeconomic status (SES) and adiposity measurements.

RESULTS

AFR admixture was significantly associated with diabetes risk in AFA women when adjusting for entry age, neighbourhood SES and BMI or waist/hip ratio (WHR) (all p < 0.0001). In HA women, AMI admixture had significant associations with diabetes risk that remained significant after adjustment for SES and BMI (all p < 0.0005). In both AFAs and HAs, SES showed significant negative associations while BMI or WHR had significant positive associations with diabetes risk, with and without adjustment for genetic admixture.

CONCLUSIONS/INTERPRETATION

In AFAs, admixture, SES and BMI/WHR each independently contribute to diabetes risk after accounting for each of the other factors; in HAs, admixture, SES and BMI each independently contribute to diabetes risk after accounting for each of the other factors, whereas admixture is not significantly associated with diabetes risk after accounting for SES and WHR. The findings emphasise the importance of considering both genetic and environmental causes in the aetiology of type 2 diabetes.

Sparse principal component analysis for identifying ancestry-informative markers in genome-wide association studies

Genome-wide association studies (GWAS) routinely apply principal component analysis (PCA) to infer population structure within a sample to correct for confounding due to ancestry. GWAS implementation of PCA uses tens of thousands of single-nucleotide polymorphisms (SNPs) to infer structure, despite the fact that only a small fraction of such SNPs provides useful information on ancestry. The identification of this reduced set of ancestry-informative markers (AIMs) from a GWAS has practical value; for example, researchers can genotype the AIM set to correct for potential confounding due to ancestry in follow-up studies that utilize custom SNP or sequencing technology. We propose a novel technique to identify AIMs from genome-wide SNP data using sparse PCA. The procedure uses penalized regression methods to identify those SNPs in a genome-wide panel that significantly contribute to the principal components while encouraging SNPs that provide negligible loadings to vanish from the analysis. We found that sparse PCA leads to negligible loss of ancestry information compared to traditional PCA analysis of genome-wide SNP data. We further demonstrate the value of sparse PCA for AIM selection using real data from the International HapMap Project and a genomewide study of inflammatory bowel disease. We have implemented our approach in open-source R software for public use.

Estrogen-related genes and their contribution to racial differences in breast cancer risk

Racial differences in breast cancer risk, including the risks of hormone receptor subtypes of breast cancer, have been previously reported. We evaluated whether variation in genes related to estrogen metabolism (COMT, CYP1A1, CYP1B1, CYP17A1, CYP19A1, ESR1, GSTM1, GSTP1, GSTT1, HSD17B1, SULT1A1, and UGT1A1) contributes to breast cancer risk and/or racial differences in risk within the CARE study, a multi-centered, population-based case-control study of breast cancer. Genetic variation was assessed as single nucleotide polymorphisms (SNPs), haplotypes, and SNP-hormone therapy (HT) interactions within a subset of 1,644 cases and 1,451 controls, including 949 Black women (493 cases and 456 controls), sampled from the CARE study population. No appreciable associations with breast cancer risk were detected for single SNPs or haplotypes in women overall. We detected SNP-HT interactions in women overall within CYP1B1 (rs1800440; p (het) = 0.003) and within CYP17A1 (rs743572; p (het) = 0.009) in which never users of HT were at a decreased risk of breast cancer, while ever users were at a non-significant increased risk. When investigated among racial groups, we detected evidence of an SNP-HT interaction with CYP1B1 in White women (p value = 0.02) and with CYP17A1 in Black women (p value = 0.04). This analysis suggests that HT use may modify the effect of variation in estrogen-related genes on breast cancer risk, which may affect Black and White women to a different extent.

Associations between genetic variants in vitamin D metabolism and asthma characteristics in young African Americans: a pilot study

INTRODUCTION

Low vitamin D levels have been associated with asthma severity in children. Young, urban African Americans (AAs) have high rates of hypovitaminosis D and asthma. Our objective was to determine associations between variants in vitamin D metabolism genes and asthma characteristics in a pilot study of young urban AAs.

MATERIALS AND METHODS

Two urban AA cohorts of subjects aged 6 to 20 years (139 subjects with asthma and 74 subjects without asthma) were genotyped for 12 single nucleotide polymorphisms (SNPs) in 3 vitamin D metabolism genes: VDR (vitamin D receptor), CYP24A1 (cytochrome P450 vitamin D 24-hydroxylase), and CYP2R1 (cytochrome P450 vitamin D 25-hydroxylase). In a case-control analysis, SNPs were studied for associations with an asthma diagnosis. Within the asthmatic cohort, SNPs were analyzed for associations with quantitative asthma characteristics. All analyses were adjusted for age, sex, and body mass index percentile.

RESULTS

Only the CYP2R1 SNP rs10766197 homozygous minor genotype was associated with asthma (P = 0.044). CYP24A1 SNP rs2248137 was associated with lower vitamin D levels (P = 0.006). Within the asthma cohort, multiple significant associations between SNPs and asthma characteristics were identified; VDR SNP rs2228570 was associated with the higher nighttime asthma morbidity scores (P = 0.04), lower baseline spirometric measures (P < 0.05), 1 or more positive aeroallergen skin test (P = 0.003), and increased immunoglobulin E levels (P < 0.001).

DISCUSSION

This pilot study demonstrates that variants in vitamin D metabolism genes are associated with quantitative asthma characteristics in young, urban AAs. The collection of these associations provides evidence for the need for a large population-based study of vitamin D-relevant SNPs in this cohort.

Non-hispanic whites have higher risk for pulmonary impairment from pulmonary tuberculosis

BACKGROUND

Disparities in outcomes associated with race and ethnicity are well documented for many diseases and patient populations. Tuberculosis (TB) disproportionately affects economically disadvantaged, racial and ethnic minority populations. Pulmonary impairment after tuberculosis (PIAT) contributes heavily to the societal burden of TB. Individual impacts associated with PIAT may vary by race/ethnicity or socioeconomic status.

METHODS

We analyzed the pulmonary function of 320 prospectively identified patients with pulmonary tuberculosis who had completed at least 20 weeks standard anti-TB regimes by directly observed therapy. We compared frequency and severity of spirometry-defined PIAT in groups stratified by demographics, pulmonary risk factors, and race/ethnicity, and examined clinical correlates to pulmonary function deficits.

RESULTS

Pulmonary impairment after tuberculosis was identified in 71% of non-Hispanic Whites, 58% of non-Hispanic Blacks, 49% of Asians and 32% of Hispanics (p < 0.001). Predictors for PIAT varied between race/ethnicity. PIAT was evenly distributed across all levels of socioeconomic status suggesting that PIAT and socioeconomic status are not related. PIAT and its severity were significantly associated with abnormal chest x-ray, p < 0.0001. There was no association between race/ethnicity and time to beginning TB treatment, p = 0.978.

CONCLUSIONS

Despite controlling for cigarette smoking, socioeconomic status and time to beginning TB treatment, non-Hispanic White race/ethnicity remained an independent predictor for disproportionately frequent and severe pulmonary impairment after tuberculosis relative to other race/ethnic groups. Since race/ethnicity was self reported and that race is not a biological construct: these findings must be interpreted with caution. However, because race/ethnicity is a proxy for several other unmeasured host, pathogen or environment factors that may contribute to disparate health outcomes, these results are meant to suggest hypotheses for further research.

Evidence for gene-gene epistatic interactions among susceptibility loci for systemic lupus erythematosus

OBJECTIVE

Several confirmed genetic susceptibility loci for lupus have been described. To date, no clear evidence for genetic epistasis in lupus has been established. The aim of this study was to test for gene-gene interactions in a number of known lupus susceptibility loci.

METHODS

Eighteen single-nucleotide polymorphisms tagging independent and confirmed lupus susceptibility loci were genotyped in a set of 4,248 patients with lupus and 3,818 normal healthy control subjects of European descent. Epistasis was tested by a 2-step approach using both parametric and nonparametric methods. The false discovery rate (FDR) method was used to correct for multiple testing.

RESULTS

We detected and confirmed gene-gene interactions between the HLA region and CTLA4, IRF5, and ITGAM and between PDCD1 and IL21 in patients with lupus. The most significant interaction detected by parametric analysis was between rs3131379 in the HLA region and rs231775 in CTLA4 (interaction odds ratio 1.19, Z = 3.95, P = 7.8 × 10(-5) [FDR ≤0.05], P for multifactor dimensionality reduction = 5.9 × 10(-45)). Importantly, our data suggest that in patients with lupus, the presence of the HLA lupus risk alleles in rs1270942 and rs3131379 increases the odds of also carrying the lupus risk allele in IRF5 (rs2070197) by 17% and 16%, respectively (P = 0.0028 and P = 0.0047, respectively).

CONCLUSION

We provide evidence for gene-gene epistasis in systemic lupus erythematosus. These findings support a role for genetic interaction contributing to the complexity of lupus heritability.

Straightforward inference of ancestry and admixture proportions through ancestry-informative insertion deletion multiplexing

Ancestry-informative markers (AIMs) show high allele frequency divergence between different ancestral or geographically distant populations. These genetic markers are especially useful in inferring the likely ancestral origin of an individual or estimating the apportionment of ancestry components in admixed individuals or populations. The study of AIMs is of great interest in clinical genetics research, particularly to detect and correct for population substructure effects in case-control association studies, but also in population and forensic genetics studies. This work presents a set of 46 ancestry-informative insertion deletion polymorphisms selected to efficiently measure population admixture proportions of four different origins (African, European, East Asian and Native American). All markers are analyzed in short fragments (under 230 basepairs) through a single PCR followed by capillary electrophoresis (CE) allowing a very simple one tube PCR-to-CE approach. HGDP-CEPH diversity panel samples from the four groups, together with Oceanians, were genotyped to evaluate the efficiency of the assay in clustering populations from different continental origins and to establish reference databases. In addition, other populations from diverse geographic origins were tested using the HGDP-CEPH samples as reference data. The results revealed that the AIM-INDEL set developed is highly efficient at inferring the ancestry of individuals and provides good estimates of ancestry proportions at the population level. In conclusion, we have optimized the multiplexed genotyping of 46 AIM-INDELs in a simple and informative assay, enabling a more straightforward alternative to the commonly available AIM-SNP typing methods dependent on complex, multi-step protocols or implementation of large-scale genotyping technologies.

African ancestry is associated with asthma risk in African Americans

Background

Asthma is a common complex condition with clear racial and ethnic differences in both prevalence and severity. Asthma consultation rates, mortality, and severe symptoms are greatly increased in African descent populations of developed countries. African ancestry has been associated with asthma, total serum IgE and lower pulmonary function in African-admixed populations. To replicate previous findings, here we aimed to examine whether African ancestry was associated with asthma susceptibility in African Americans. In addition, we examined for the first time whether African ancestry was associated with asthma exacerbations.

Methodology/Principal Findings

After filtering for self-reported ancestry and genotype data quality, samples from 1,117 self-reported African-American individuals from New York and Baltimore (394 cases, 481 controls), and Chicago (321 cases followed for asthma exacerbations) were analyzed. Genetic ancestry was estimated based on ancestry informative markers (AIMs) selected for being highly divergent among European and West African populations (95 AIMs for New York and Baltimore, and 66 independent AIMs for Chicago). Among case-control samples, the mean African ancestry was significantly higher in asthmatics than in non-asthmatics (82.0±14.0% vs. 77.8±18.1%, mean difference 4.2% [95% confidence interval (CI):2.0–6.4], p<0.0001). This association remained significant after adjusting for potential confounders (odds ratio: 4.55, 95% CI: 1.69–12.29, p = 0.003). African ancestry failed to show an association with asthma exacerbations (p = 0.965) using a model based on longitudinal data of the number of exacerbations followed over 1.5 years.

Conclusions/Significance

These data replicate previous findings indicating that African ancestry constitutes a risk factor for asthma and suggest that elevated asthma rates in African Americans can be partially attributed to African genetic ancestry.

Indigenous American ancestry is associated with arsenic methylation efficiency in an admixed population of northwest Mexico

Many studies provide evidence relating lower human arsenic (As) methylation efficiency, represented by high percent urinary monomethylarsonic acid (MMA(V)), with several As-induced diseases, possibly due to the fact that MMA(V) serves as a proxy for MMA(III), the most toxic As metabolite. Some epidemiological studies suggested that indigenous Americans (AME) methylate As more efficiently; however, data supporting this have been equivocal. The aim of this study was to characterize the association between AME ancestry and As methylation efficiency using a panel of ancestry informative genetic markers to determine individual ancestry proportions in an admixed population (composed of two or more isolated ancestral populations) of 746 individuals environmentally exposed to As in northwest Mexico. Total urinary As (TAs) mean and range were 170.4 and 2.3-1053.5 μg/L, while percent AME (%AME) mean and range were 72.4 and 23-100. Adjusted (gender, age, AS3MT 7388/M287T haplotypes, body mass index [BMI], and TAs) multiple regression model showed that higher AME ancestry is significantly associated with lower percentage of urinary As excreted as MMA(V) (%uMMA) in this population (p < .01). Data also demonstrated a significant interaction between BMI and gender, indicating negative association between BMI and %uMMA, stronger in women than men (p < .01). Moreover, age and the AS3MT variants 7388 (intronic) and M287T (nonsynonymous) were also significantly associated with As methylation efficiency (p < .01). This study highlights the importance of BMI and indigenous American ancestry in some of the observed variability in As methylation efficiency, underscoring the need to be considered in epidemiology studies, particularly those carried out in admixed populations.

Identifying genetic variation for alcohol dependence

Researchers are using various strategies to identify the genes that may be associated with alcoholism. The initial efforts primarily relied on candidate gene and linkage studies; more recently, however, modern advances in genotyping have resulted in widespread use of genome-wide association studies for alcohol dependence. The key findings of the earlier studies were that variations (i.e., polymorphisms) in the DNA sequences of the genes encoding alcohol dehydrogenase 1B (i.e., the ADH1B gene), aldehyde dehydrogenase 2 (i.e., the ALDH2 gene), and other alcohol-metabolizing enzymes mediate the risk for alcoholism; moreover, these polymorphisms also have an impact on the risk of alcohol-related cancers, such as esophageal cancer. In addition, a gene encoding one of the receptors for the neurotransmitter γ-aminobutyric acid (GABA) known as GABRA2 seems to have a role in the development of alcohol dependence. Genome-wide association studies now offer a host of emerging opportunities, as well as challenges, for discovering the genetic etiology of alcohol dependence and for unveiling new treatment strategies.

The rs4774 CIITA missense variant is associated with risk of systemic lupus erythematosus

The major histocompatibility complex (MHC) class II transactivator gene (CIITA) encodes an important transcription factor required for human leukocyte antigens (HLA) class II MHC-restricted antigen presentation. MHC genes, including the HLA class II DRB1*03:01 allele, are strongly associated with systemic lupus erythematosus (SLE). Recently the rs4774 CIITA missense variant (+1632G/C) was reported to be associated with susceptibility to multiple sclerosis. In the current study, we investigated CIITA, DRB1*03:01 and risk of SLE using a multi-stage analysis. In stage 1, 9 CIITA variants were tested in 658 cases and 1363 controls (N=2021). In stage 2, rs4774 was tested in 684 cases and 2938 controls (N=3622). We also performed a meta-analysis of the pooled 1342 cases and 4301 controls (N=5643). In stage 1, rs4774(*)C was associated with SLE (odds ratio (OR)=1.24, 95% confidence interval (95% CI)=1.07-1.44, P=4.2 × 10(-3)). Similar results were observed in stage 2 (OR=1.16, 95% CI=1.02-1.33, P=8.5 × 10(-3)) and the meta-analysis of the combined data set (OR=1.20, 95% CI=1.09-1.33, P(meta)=2.5 × 10(-4)). In all three analyses, the strongest evidence for association between rs4774(*)C and SLE was present in individuals who carried at least one copy of DRB1*03:01 (P(meta)=1.9 × 10(-3)). Results support a role for CIITA in SLE, which appears to be stronger in the presence of DRB1*03:01.

Genetic associations with metabolic syndrome and its quantitative traits by race/ethnicity in the United States

BACKGROUND

Elevated insulin resistance (IR), triglycerides (TG), body mass index (BMI), and waist circumference (WC) are features of the metabolic syndrome. Although several single-nucleotide polymorphisms (SNPs) associated with these traits have been reported, no study has reported their risk allele frequencies and effect sizes among the major U.S. race/ethnic groups in a nationally representative sample.

METHODS

We compared the risk allele frequencies of eight SNPs previously associated with IR, TG, BMI, or WC by race/ethnicity (non-Hispanic white, non-Hispanic black, Mexican American) in 3,030 participants of the National Health and Nutrition Examination Study III (NHANES III). In regression models predicting IR, TG, BMI, WC, and metabolic syndrome, we tested whether the SNP effect sizes on these traits varied by race/ethnicity.

RESULTS

Risk allele frequencies varied by race/ethnicity for all eight loci (P<0.0001). The directionality of effects of the variants on IR, TG, WC, and BMI was generally consistent with previous observations and did not differ by race/ethnicity (P>0.001), although our study had low power for this test. No SNP predicted metabolic syndrome in any of the three groups (P>0.05).

CONCLUSIONS

The significance of racial/ethnic differences in risk allele frequencies merits consideration if genetic discoveries are to have clinical and public health applicability.

Type 2 diabetes risk variants and colorectal cancer risk: the Multiethnic Cohort and PAGE studies

BACKGROUND

Diabetes has been positively associated with the risk of colorectal cancer. This study investigated whether recently established risk variants for diabetes also have effects on colorectal cancer.

METHODS

19 single nucleotide repeats (SNPs) associated with type 2 diabetes in genome-wide association studies were tested in a case-control study of 2011 colorectal cancer cases and 6049 controls nested in the Multiethnic Cohort study as part of the Population Architecture using Genomics and Epidemiology (PAGE) initiative. ORs and 95% CIs were estimated by unconditional logistic regression to evaluate the association between SNPs and colorectal cancer risk, adjusting for age, sex and race/ethnicity. Permutation testing was conducted to correct for multiple hypothesis testing.

RESULTS

Four type 2 diabetes SNPs were associated with colorectal cancer risk: rs7578597 (THADA), rs864745 (JAZF1), rs5219 (KCNJ11) and rs7961581 (TSPAN8, LGR5). The strongest association was for the rs7578597 (THADA) Thr1187Ala missense polymorphism (P(trend)=0.004 adjusted for multiple testing), with the high risk allele for colorectal cancer being the low risk allele for diabetes. Similar patterns of associations were seen with further adjustment for diabetes status and body mass index. The association of diabetes status with colorectal cancer risk was somewhat weakened after adjustment for these SNPs.

CONCLUSION

The findings suggest that diabetes risk variants also influence colorectal cancer susceptibility, possibly through mechanisms different from those for diabetes.

Nuclear receptor-mediated induction of CYP450 by antiretrovirals: functional consequences of NR1I2 (PXR) polymorphisms and differential prevalence in whites and sub-Saharan Africans

BACKGROUND

Antiretroviral therapy including HIV protease inhibitors and nonnucleoside reverse transcriptase inhibitors can both inhibit and induce expression of cytochrome P450s, potentially leading to drug interactions. However, information is lacking on the impact of genetic polymorphism on this interaction.

METHODS

This study examines the prevalence of 33 polymorphisms in NR1I2 (pregnane X receptor [PXR]), CYP3A4, and CYP2B6 in 1013 white and sub-Saharan African patients with HIV; explores the inductive ability of 16 antiretrovirals on CYP3A4 and CYP2B6 promoter activity through nuclear receptors PXR and constitutive androstane receptor (CAR); and evaluates the influence of naturally occurring PXR genetic variants on antiretroviral activation.

RESULTS

Seventeen polymorphisms were present at different frequencies between the two ethnicities. Darunavir, fosamprenavir, lopinavir, nelfinavir, tipranavir, efavirenz, and abacavir increased CYP3A4 and/or CYP2B6 promoter activity, some through constitutive androstane receptor but mainly through PXR. Addition of low-dose ritonavir enhanced levels of CYP promoter activity for several protease inhibitors. Some PXR variants displayed lower fosamprenavir- and lopinavir-induced CYP3A4 promoter activity than the PXR reference sequence, whereas efavirenz and nelfinavir induction was unchanged.

CONCLUSIONS

The presence of NR1I2 polymorphisms can alter the induction of CYP3A4 and CYP2B6 promoter activity, potentially adding to the unpredictable nature of antiretroviral drug interactions. These polymorphisms differ in prevalence between whites and sub-Saharan Africans.

Identification of novel genetic susceptibility loci in African American lupus patients in a candidate gene association study

OBJECTIVE

Candidate gene and genome-wide association studies have identified several disease susceptibility loci in lupus patients. These studies have largely been performed in lupus patients who are Asian or of European ancestry. This study was undertaken to examine whether some of these same susceptibility loci increase lupus risk in African American individuals.

METHODS

Single-nucleotide polymorphisms tagging 15 independent lupus susceptibility loci were genotyped in a set of 1,724 lupus patients and 2,024 healthy controls of African American descent. The loci examined included PTPN22, FCGR2A, TNFSF4, STAT4, CTLA4, PDCD1, PXK, BANK1, MSH5 (HLA region), CFB (HLA region), C8orf13-BLK region, MBL2, KIAA1542, ITGAM, and MECP2/IRAK1.

RESULTS

We found the first evidence of genetic association between lupus in African American patients and 5 susceptibility loci (C8orf13-BLK, BANK1, TNFSF4, KIAA1542, and CTLA4; P = 8.0 × 10⁻⁶, P = 1.9 × 10⁻⁵, P = 5.7 × 10⁻⁵, P = 0.00099, and P = 0.0045, respectively). Further, we confirmed the genetic association between lupus and 5 additional lupus susceptibility loci (ITGAM, MSH5, CFB, STAT4, and FCGR2A; P = 7.5 × 10⁻¹¹, P = 5.2 × 10⁻⁸, P = 8.7 × 10⁻⁷ , P = 0.0058, and P = 0.0070, respectively), and provided evidence, for the first time, of genome-wide significance for the association between lupus in African American patients and ITGAM and MSH5 (HLA region).

CONCLUSION

These findings provide evidence of novel genetic susceptibility loci for lupus in African Americans and demonstrate that the majority of lupus susceptibility loci examined confer lupus risk across multiple ethnicities.

BRCA1 and BRCA2 mutations among ovarian cancer patients from Colombia

OBJECTIVE

The contribution of BRCA1 and BRCA2 mutations to ovarian cancer in Colombia has not yet been explored. Five founder mutations have been identified in two previous studies of breast cancer patients in the Bogota region [1,2]. It is important that the frequency of mutations be established among unselected cases of ovarian cancer in order to estimate the genetic burden of this cancer in Colombia and to plan genetic and preventive services.

METHODS

We enrolled 100 unselected women with ovarian cancer from the Bogota region, and from northern and southern central regions of Colombia. A detailed family history was obtained from each patient and a blood sample was processed for DNA analysis. DNA quality was adequate for BRCA testing for 96 women. Mutations in BRCA1 and BRCA2 were sought using a Hispanic BRCA mutation testing panel. All mutations were confirmed by direct sequencing.

RESULTS

Fifteen mutations were identified (two in BRCA2 and thirteen in BRCA1) representing 15.6% of the total (95% CI: 7.8% to 21.3%). Among the 15 mutation-positive families there were nine breast-ovarian cancer families, one gastric cancer family, one prostate cancer family, three uterine cancer families, and one family with no history of cancer. A single founder mutation in BRCA1 (3450del4) was seen in 11 patients.

CONCLUSION

In summary, BRCA1 founder mutations are common in Colombian women with ovarian cancer. Approximately 11.5% of all ovarian cancer cases in the Bogota region are attributable to a single BRCA1 founder mutation.

ALFRED: an allele frequency resource for research and teaching

ALFRED (http://alfred.med.yale.edu) is a free, web accessible, curated compilation of allele frequency data on DNA sequence polymorphisms in anthropologically defined human populations. Currently, ALFRED has allele frequency tables on over 663 400 polymorphic sites; 170 of them have frequency tables for more than 100 different population samples. In ALFRED, a population may have multiple samples with each ‘sample’ consisting of many individuals on which an allele frequency is based. There are 3566 population samples from 710 different populations with allele frequency tables on at least one polymorphism. Fifty of those population samples have allele frequency data for over 650 000 polymorphisms. Records also have active links to relevant resources (dbSNP, PharmGKB, OMIM, Ethnologue, etc.). The flexible search options and data display and download capabilities available through the web interface allow easy access to the large quantity of high-quality data in ALFRED.

Phenotypic associations of genetic susceptibility loci in systemic lupus erythematosus

OBJECTIVE

Systemic lupus erythematosus is a clinically heterogeneous autoimmune disease. A number of genetic loci that increase lupus susceptibility have been established. This study examines if these genetic loci also contribute to the clinical heterogeneity in lupus.

MATERIALS AND METHODS

4001 European-derived, 1547 Hispanic, 1590 African-American and 1191 Asian lupus patients were genotyped for 16 confirmed lupus susceptibility loci. Ancestry informative markers were genotyped to calculate and adjust for admixture. The association between the risk allele in each locus was determined and compared in patients with and without the various clinical manifestations included in the ACR criteria.

RESULTS

Renal disorder was significantly correlated with the lupus risk allele in ITGAM (p=5.0 × 10(-6), OR 1.25, 95% CI 1.12 to 1.35) and in TNFSF4 (p=0.0013, OR 1.14, 95% CI 1.07 to 1.25). Other significant findings include the association between risk alleles in FCGR2A and malar rash (p=0.0031, OR 1.11, 95% CI 1.17 to 1.33), ITGAM and discoid rash (p=0.0020, OR 1.20, 95% CI 1.06 to 1.33), STAT4 and protection from oral ulcers (p=0.0027, OR 0.89, 95% CI 0.83 to 0.96) and IL21 and haematological disorder (p=0.0027, OR 1.13, 95% CI 1.04 to 1.22). All these associations are significant with a false discovery rate of <0.05 and pass the significance threshold using Bonferroni correction for multiple testing.

CONCLUSION

Signifi cant associations were found between clinical manifestations and the FCGR2A, ITGAM, STAT4, TNSF4 and IL21 genes. The findings suggest that genetic profiling might be a useful tool to predict disease manifestations in lupus patients in the future.

The Next PAGE in understanding complex traits: design for the analysis of Population Architecture Using Genetics and Epidemiology (PAGE) Study

Genetic studies have identified thousands of variants associated with complex traits. However, most association studies are limited to populations of European descent and a single phenotype. The Population Architecture using Genomics and Epidemiology (PAGE) Study was initiated in 2008 by the National Human Genome Research Institute to investigate the epidemiologic architecture of well-replicated genetic variants associated with complex diseases in several large, ethnically diverse population-based studies. Combining DNA samples and hundreds of phenotypes from multiple cohorts, PAGE is well-suited to address generalization of associations and variability of effects in diverse populations; identify genetic and environmental modifiers; evaluate disease subtypes, intermediate phenotypes, and biomarkers; and investigate associations with novel phenotypes. PAGE investigators harmonize phenotypes across studies where possible and perform coordinated cohort-specific analyses and meta-analyses. PAGE researchers are genotyping thousands of genetic variants in up to 121,000 DNA samples from African-American, white, Hispanic/Latino, Asian/Pacific Islander, and American Indian participants. Initial analyses will focus on single nucleotide polymorphisms (SNPs) associated with obesity, lipids, cardiovascular disease, type 2 diabetes, inflammation, various cancers, and related biomarkers. PAGE SNPs are also assessed for pleiotropy using the “phenome-wide association study” approach, testing each SNP for associations with hundreds of phenotypes. PAGE data will be deposited into the National Center for Biotechnology Information's Database of Genotypes and Phenotypes and made available via a custom browser.

Race, ancestry, and development of food-allergen sensitization in early childhood

OBJECTIVE

We examined whether the risk of food-allergen sensitization varied according to self-identified race or genetic ancestry.

METHODS

We studied 1104 children (mean age: 2.7 years) from an urban multiethnic birth cohort. Food sensitization was defined as specific immunoglobulin E (sIgE) levels of ≥ 0.35 kilo-units of allergen (kUA)/L for any of 8 common food allergens. Multivariate logistic regression analyses were used to evaluate the associations of self-identified race and genetic ancestry with food sensitization. Analyses also examined associations with numbers of food sensitizations (0, 1 or 2, and ≥ 3 foods) and with logarithmically transformed allergen sIgE levels.

RESULTS

In this predominantly minority cohort (60.9% black and 22.5% Hispanic), 35.5% of subjects exhibited food sensitizations. In multivariate models, both self-reported black race (odds ratio [OR]: 2.34 [95% confidence interval [CI]: 1.24-4.44]) and African ancestry (in 10% increments; OR: 1.07 [95% CI: 1.02-1.14]) were associated with food sensitization. Self-reported black race (OR: 3.76 [95% CI: 1.09-12.97]) and African ancestry (OR: 1.19 [95% CI: 1.07-1.32]) were associated with a high number (≥ 3) of food sensitizations. African ancestry was associated with increased odds of peanut sIgE levels of ≥ 5 kUA/L (OR: 1.25 [95% CI: 1.01-1.52]). Similar ancestry associations were seen for egg sIgE levels of ≥ 2 kUA/L (OR: 1.13 [95% CI: 1.01-1.27]) and milk sIgE levels of ≥ 5 kUA/L (OR: 1.24 [95% CI: 0.94-1.63]), although findings were not significant for milk.

CONCLUSIONS

Black children were more likely to be sensitized to food allergens and were sensitized to more foods. African ancestry was associated with peanut sensitization.

Ancestral heterogeneity in a biethnic stroke population

To test for and characterize heterogeneity in ancestral contributions to individuals among a population of Mexican American (MA) and non-Hispanic white (NHW) stroke/transient ischemic attack (TIA) cases, data from a community-based stroke surveillance study in south Texas were used. Strokes/TIA cases were identified (2004-2006) with a random sample asked to provide blood. Race-ethnicity was self-reported. Thirty-three ancestry informative markers were genotyped and individual genetic admixture estimated using maximum likelihood methods. Three hypotheses were tested for each MA using likelihood ratio tests: (1) H(0) : μi = 0 (100% Native American), (2) H(0) : μi = 1.00 (100% European), (3) H(0) : μi = 0.59 (average European). Among 154 self-identified MAs, estimated European ancestry varied from 0.26 to 0.98, with an average of 0.59 (SE = 0.014). We rejected hypothesis 1 for every MA and rejected hypothesis 2 for all but two MAs. We rejected hypothesis 3 for 40 MAs (20 < 59%, 20 > 59%). Among 84 self-identified NHWs, the estimated fraction of European ancestry ranged from 0.83 to 1.0, with an average of 0.97 (SE = 0.014). Self-identified MAs, and to a lesser extent NHWs, from an established bi-ethnic community were heterogeneous with respect to genetic admixture. Researchers should not use simple race-ethnic categories as proxies for homogeneous genetic populations when conducting gene mapping and disease association studies in multi-ethnic populations.

African ancestry, early life exposures, and respiratory morbidity in early childhood

BACKGROUND

Racial disparities persist in early childhood wheezing and cannot be completely explained by known risk factors.

OBJECTIVE

To evaluate the associations of genetic ancestry and self-identified race with early childhood recurrent wheezing, accounting for socio-economic status (SES) and early life exposures.

METHODS

We studied 1034 children in an urban, multi-racial, prospective birth cohort. Multivariate logistic regression was used to evaluate the association of genetic ancestry as opposed to self-identified race with recurrent wheezing (>3 episodes). Sequential models accounted for demographic, socio-economic factors and early life risk factors. Genetic ancestry, estimated using 150 ancestry informative markers, was expressed in deciles.

RESULTS

Approximately 6.1% of subjects (mean age 3.1 years) experienced recurrent wheezing. After accounting for SES and demographic factors, African ancestry (OR: 1.16, 95% CI: 1.02-1.31) was significantly associated with recurrent wheezing. By self-reported race, hispanic subjects had a borderline decrease in risk of wheeze compared with African Americans (OR: 0.44, 95% CI: 0.19-1.00), whereas white subjects (OR: 0.46, 95% CI: 0.14-1.57) did not have. After further adjustment for known confounders and early life exposures, both African (OR: 1.19, 95% CI: 1.05-1.34) and European ancestry (OR: 0.84, 95% CI: 0.74-0.94) retained a significant association with recurrent wheezing, as compared with self-identified race (OR(whites) : 0.31, 95% CI: 0.09-1.14; OR(hispanic) : 0.47, 95% CI: 0.20-1.08). There were no significant interactions between ancestry and early life factors on recurrent wheezing.

CONCLUSIONS AND CLINICAL RELEVANCE

In contrast to self-identified race, African ancestry remained a significant, independent predictor of early childhood wheezing after accounting for early life and other known risk factors associated with lung function changes and asthma. Genetic ancestry may be a powerful way to evaluate wheezing disparities and a proxy for differentially distributed genetic and early life risk factors associated with childhood recurrent wheezing.