Enhanced statistical tests for GWAS in admixed populations: assessment using African Americans from CARe and a Breast Cancer Consortium

Improved multi-ancestry fine-mapping identifies cis-regulatory variants underlying molecular traits and disease risk

Multi-ancestry statistical fine-mapping of cis-molecular quantitative trait loci (cis-molQTL) aims to improve the precision of distinguishing causal cis-molQTLs from tagging variants. However, existing approaches fail to reflect shared genetic architectures. To solve this limitation, we present the Sum of Shared Single Effects (SuShiE) model, which leverages LD heterogeneity to improve fine-mapping precision, infer cross-ancestry effect size correlations, and estimate ancestry-specific expression prediction weights. We apply SuShiE to mRNA expression measured in PBMCs (n=956) and LCLs (n=814) together with plasma protein levels (n=854) from individuals of diverse ancestries in the TOPMed MESA and GENOA studies. We find SuShiE fine-maps cis-molQTLs for 16% more genes compared with baselines while prioritizing fewer variants with greater functional enrichment. SuShiE infers highly consistent cis-molQTL architectures across ancestries on average; however, we also find evidence of heterogeneity at genes with predicted loss-of-function intolerance, suggesting that environmental interactions may partially explain differences in cis-molQTL effect sizes across ancestries. Lastly, we leverage estimated cis-molQTL effect-sizes to perform individual-level TWAS and PWAS on six white blood cell-related traits in AOU Biobank individuals (n=86k), and identify 44 more genes compared with baselines, further highlighting its benefits in identifying genes relevant for complex disease risk. Overall, SuShiE provides new insights into the cis-genetic architecture of molecular traits.

Admix-kit: an integrated toolkit and pipeline for genetic analyses of admixed populations

SUMMARY

Admixed populations, with their unique and diverse genetic backgrounds, are often underrepresented in genetic studies. This oversight not only limits our understanding but also exacerbates existing health disparities. One major barrier has been the lack of efficient tools tailored for the special challenges of genetic studies of admixed populations. Here, we present admix-kit, an integrated toolkit and pipeline for genetic analyses of admixed populations. Admix-kit implements a suite of methods to facilitate genotype and phenotype simulation, association testing, genetic architecture inference, and polygenic scoring in admixed populations.

AVAILABILITY AND IMPLEMENTATION

Admix-kit package is open-source and available at https://github.com/KangchengHou/admix-kit. Additionally, users can use the pipeline designed for admixed genotype simulation available at https://github.com/UW-GAC/admix-kit_workflow.

Phase-free local ancestry inference mitigates the impact of switch errors on phase-based methods

Local ancestry inference (LAI) is an indispensable component of a variety of analyses in medical and population genetics, from admixture mapping to characterizing demographic history. However, the accuracy of LAI depends on a number of factors such as phase quality (for phase-based LAI methods), time since admixture of the population under study, and other factors. Here we present an empirical analysis of four LAI methods using simulated individuals of mixed African and European ancestry, examining the impact of variable phase quality and a range of demographic scenarios. We found that regardless of phasing options, calls from LAI methods that operate on unphased genotypes (phase-free LAI) have 2.6-4.6% higher Pearson correlation with the ground truth than methods that operate on phased genotypes (phase-based LAI). Applying the TRACTOR phase-correction algorithm led to modest improvements in phase-based LAI, but despite this, the Pearson correlation of phase-free LAI remained 2.4-3.8% higher than phase-corrected phase-based approaches (considering the best performing methods in each category). Phase-free and phase-based LAI accuracy differences can dramatically impact downstream analyses: estimates of the time since admixture using phase-based LAI tracts are upwardly biased by ≈10 generations using our highest quality phased data but have virtually no bias using phase-free LAI calls. Our study underscores the strong dependence of phase-based LAI accuracy on phase quality and highlights the merits of LAI approaches that analyze unphased genetic data.

Joint multi-ancestry and admixed GWAS reveals the complex genetics behind human cranial vault shape

The cranial vault in humans is highly variable, clinically relevant, and heritable, yet its genetic architecture remains poorly understood. Here, we conduct a joint multi-ancestry and admixed multivariate genome-wide association study on 3D cranial vault shape extracted from magnetic resonance images of 6772 children from the ABCD study cohort yielding 30 genome-wide significant loci. Follow-up analyses indicate that these loci overlap with genomic risk loci for sagittal craniosynostosis, show elevated activity cranial neural crest cells, are enriched for processes related to skeletal development, and are shared with the face and brain. We present supporting evidence of regional localization for several of the identified genes based on expression patterns in the cranial vault bones of E15.5 mice. Overall, our study provides a comprehensive overview of the genetics underlying normal-range cranial vault shape and its relevance for understanding modern human craniofacial diversity and the etiology of congenital malformations.

Universal genome-wide association studies: Powerful joint ancestry and association testing

The vast majority of human populations and individuals have mixed ancestry. Consequently, adjustment for locus-specific ancestry is essential for genetic association studies. To empower association studies for all populations, it is necessary to integrate effects of locus-specific ancestry and genotype. We developed a joint test of ancestry and association that can be performed with summary statistics, is independent of study design, can take advantage of locus-specific ancestry effects to boost power in association testing, and can utilize association effects to fine map admixture peaks. We illustrate the test using the association between serum triglycerides and LPL. By combining data from African Americans, European Americans, and West Africans, we identify three conditionally independent variants with varying amounts of ancestrally differentiated allele frequencies. Using out-of-sample data, we demonstrate improved prediction achievable by accounting for multiple causal variants and locus-specific ancestry effects at a single locus.

Polygenic prediction across populations is influenced by ancestry, genetic architecture, and methodology

Polygenic risk scores (PRSs) developed from multi-ancestry genome-wide association studies (GWASs), PRSmulti, hold promise for improving PRS accuracy and generalizability across populations. To establish best practices for leveraging the increasing diversity of genomic studies, we investigated how various factors affect the performance of PRSmulti compared with PRSs constructed from single-ancestry GWASs (PRSsingle). Through extensive simulations and empirical analyses, we showed that PRSmulti overall outperformed PRSsingle in understudied populations, except when the understudied population represented a small proportion of the multi-ancestry GWAS. Furthermore, integrating PRSs based on local ancestry-informed GWASs and large-scale, European-based PRSs improved predictive performance in understudied African populations, especially for less polygenic traits with large-effect ancestry-enriched variants. Our work highlights the importance of diversifying genomic studies to achieve equitable PRS performance across ancestral populations and provides guidance for developing PRSs from multiple studies.

Admix-kit: An Integrated Toolkit and Pipeline for Genetic Analyses of Admixed Populations

Admixed populations, with their unique and diverse genetic backgrounds, are often underrepresented in genetic studies. This oversight not only limits our understanding but also exacerbates existing health disparities. One major barrier has been the lack of efficient tools tailored for the special challenges of genetic study of admixed populations. Here, we present admix-kit, an integrated toolkit and pipeline for genetic analyses of admixed populations. Admix-kit implements a suite of methods to facilitate genotype and phenotype simulation, association testing, genetic architecture inference, and polygenic scoring in admixed populations.

Strategies for the Genomic Analysis of Admixed Populations

Admixed populations constitute a large portion of global human genetic diversity, yet they are often left out of genomics analyses. This exclusion is problematic, as it leads to disparities in the understanding of the genetic structure and history of diverse cohorts and the performance of genomic medicine across populations. Admixed populations have particular statistical challenges, as they inherit genomic segments from multiple source populations-the primary reason they have historically been excluded from genetic studies. In recent years, however, an increasing number of statistical methods and software tools have been developed to account for and leverage admixture in the context of genomics analyses. Here, we provide a survey of such computational strategies for the informed consideration of admixture to allow for the well-calibrated inclusion of mixed ancestry populations in large-scale genomics studies, and we detail persisting gaps in existing tools.

Impact of cross-ancestry genetic architecture on GWASs in admixed populations

Genome-wide association studies (GWASs) have identified thousands of variants for disease risk. These studies have predominantly been conducted in individuals of European ancestries, which raises questions about their transferability to individuals of other ancestries. Of particular interest are admixed populations, usually defined as populations with recent ancestry from two or more continental sources. Admixed genomes contain segments of distinct ancestries that vary in composition across individuals in the population, allowing for the same allele to induce risk for disease on different ancestral backgrounds. This mosaicism raises unique challenges for GWASs in admixed populations, such as the need to correctly adjust for population stratification. In this work we quantify the impact of differences in estimated allelic effect sizes for risk variants between ancestry backgrounds on association statistics. Specifically, while the possibility of estimated allelic effect-size heterogeneity by ancestry (HetLanc) can be modeled when performing a GWAS in admixed populations, the extent of HetLanc needed to overcome the penalty from an additional degree of freedom in the association statistic has not been thoroughly quantified. Using extensive simulations of admixed genotypes and phenotypes, we find that controlling for and conditioning effect sizes on local ancestry can reduce statistical power by up to 72%. This finding is especially pronounced in the presence of allele frequency differentiation. We replicate simulation results using 4,327 African-European admixed genomes from the UK Biobank for 12 traits to find that for most significant SNPs, HetLanc is not large enough for GWASs to benefit from modeling heterogeneity in this way.

Causal effects on complex traits are similar for common variants across segments of different continental ancestries within admixed individuals

Individuals of admixed ancestries (for example, African Americans) inherit a mosaic of ancestry segments (local ancestry) originating from multiple continental ancestral populations. This offers the unique opportunity of investigating the similarity of genetic effects on traits across ancestries within the same population. Here we introduce an approach to estimate correlation of causal genetic effects (radmix) across local ancestries and analyze 38 complex traits in African-European admixed individuals (N = 53,001) to observe very high correlations (meta-analysis radmix = 0.95, 95% credible interval 0.93-0.97), much higher than correlation of causal effects across continental ancestries. We replicate our results using regression-based methods from marginal genome-wide association study summary statistics. We also report realistic scenarios where regression-based methods yield inflated heterogeneity-by-ancestry due to ancestry-specific tagging of causal effects, and/or polygenicity. Our results motivate genetic analyses that assume minimal heterogeneity in causal effects by ancestry, with implications for the inclusion of ancestry-diverse individuals in studies.

Local Ancestry-Informed Candidate Pathway Analysis of Warfarin Stable Dose in Latino Populations

Accuracy of warfarin dose prediction algorithms may be improved by including data from diverse populations in genetic studies of dose variability. Here, we surveyed single nucleotide polymorphisms in vitamin K-related genetic pathways for association with warfarin dose requirements in two admixed Latino populations in standard-principal component adjusted and contemporary-local ancestry adjusted regression models. A total of five variants from vitamin K-related genes/pathways were associated with warfarin dose in both cohorts (P < 0.0125) in standard models. Local ancestry-adjusted analysis unveiled 35 associated variants with absolute effects ranging from β = 9.04 ( ±2.23) to 39.18 ( ±10.89) per ancestral allele in the discovery cohort and β = 6.47 (± 2.02) to 17.82 (± 6.83) in the replication cohort. Importantly, we demonstrate the technical validity of the Tractor model in cohorts with admixed ancestry from three founder populations and bring attention to the technical hurdles obstructing the inclusion of diverse, especially admixed, populations in pharmacogenomic research.

Impact of cross-ancestry genetic architecture on GWAS in admixed populations

Genome-wide association studies (GWAS) have identified thousands of variants for disease risk. These studies have predominantly been conducted in individuals of European ancestries, which raises questions about their transferability to individuals of other ancestries. Of particular interest are admixed populations, usually defined as populations with recent ancestry from two or more continental sources. Admixed genomes contain segments of distinct ancestries that vary in composition across individuals in the population, allowing for the same allele to induce risk for disease on different ancestral backgrounds. This mosaicism raises unique challenges for GWAS in admixed populations, such as the need to correctly adjust for population stratification to balance type I error with statistical power. In this work we quantify the impact of differences in estimated allelic effect sizes for risk variants between ancestry backgrounds on association statistics. Specifically, while the possibility of estimated allelic effect-size heterogeneity by ancestry (HetLanc) can be modeled when performing GWAS in admixed populations, the extent of HetLanc needed to overcome the penalty from an additional degree of freedom in the association statistic has not been thoroughly quantified. Using extensive simulations of admixed genotypes and phenotypes we find that modeling HetLanc in its absence reduces statistical power by up to 72%. This finding is especially pronounced in the presence of allele frequency differentiation. We replicate simulation results using 4,327 African-European admixed genomes from the UK Biobank for 12 traits to find that for most significant SNPs HetLanc is not large enough for GWAS to benefit from modeling heterogeneity.

Including diverse and admixed populations in genetic epidemiology research

The inclusion of ancestrally diverse participants in genetic studies can lead to new discoveries and is important to ensure equitable health care benefit from research advances. Here, members of the Ethical, Legal, Social, Implications (ELSI) committee of the International Genetic Epidemiology Society (IGES) offer perspectives on methods and analysis tools for the conduct of inclusive genetic epidemiology research, with a focus on admixed and ancestrally diverse populations in support of reproducible research practices. We emphasize the importance of distinguishing socially defined population categorizations from genetic ancestry in the design, analysis, reporting, and interpretation of genetic epidemiology research findings. Finally, we discuss the current state of genomic resources used in genetic association studies, functional interpretation, and clinical and public health translation of genomic findings with respect to diverse populations.

Allelic imbalance of chromatin accessibility in cancer identifies candidate causal risk variants and their mechanisms

While many germline cancer risk variants have been identified through genome-wide association studies (GWAS), the mechanisms by which these variants operate remain largely unknown. Here we used 406 cancer ATAC-Seq samples across 23 cancer types to identify 7,262 germline allele-specific accessibility QTLs (as-aQTLs). Cancer as-aQTLs had stronger enrichment for cancer risk heritability (up to 145 fold) than any other functional annotation across seven cancer GWAS. Most cancer as-aQTLs directly altered transcription factor (TF) motifs and exhibited differential TF binding and gene expression in functional screens. To connect as-aQTLs to putative risk mechanisms, we introduced the regulome-wide associations study (RWAS). RWAS identified genetically associated accessible peaks at >70% of known breast and prostate loci and discovered new risk loci in all examined cancer types. Integrating as-aQTL discovery, motif analysis and RWAS identified candidate causal regulatory elements and their probable upstream regulators. Our work establishes cancer as-aQTLs and RWAS analysis as powerful tools to study the genetic architecture of cancer risk.

Ethnic Diversity and Warfarin Pharmacogenomics

Warfarin has remained the most commonly prescribed vitamin K oral anticoagulant worldwide since its approval in 1954. Dosing challenges including having a narrow therapeutic window and a wide interpatient variability in dosing requirements have contributed to making it the most studied drug in terms of genotype-phenotype relationships. However, most of these studies have been conducted in Whites or Asians which means the current pharmacogenomics evidence-base does not reflect ethnic diversity. Due to differences in minor allele frequencies of key genetic variants, studies conducted in Whites/Asians may not be applicable to underrepresented populations such as Blacks, Hispanics/Latinos, American Indians/Alaska Natives and Native Hawaiians/other Pacific Islanders. This may exacerbate health inequalities when Whites/Asians have better anticoagulation profiles due to the existence of validated pharmacogenomic dosing algorithms which fail to perform similarly in the underrepresented populations. To examine the extent to which individual races/ethnicities are represented in the existing body of pharmacogenomic evidence, we review evidence pertaining to published pharmacogenomic dosing algorithms, including clinical utility studies, cost-effectiveness studies and clinical implementation guidelines that have been published in the warfarin field.

Human genetic admixture through the lens of population genomics

Over the past 50 years, geneticists have made great strides in understanding how our species' evolutionary history gave rise to current patterns of human genetic diversity classically summarized by Lewontin in his 1972 paper, ‘The Apportionment of Human Diversity’. One evolutionary process that requires special attention in both population genetics and statistical genetics is admixture: gene flow between two or more previously separated source populations to form a new admixed population. The admixture process introduces ancestry-based structure into patterns of genetic variation within and between populations, which in turn influences the inference of demographic histories, identification of genetic targets of selection and prediction of complex traits. In this review, we outline some challenges for admixture population genetics, including limitations of applying methods designed for populations without recent admixture to the study of admixed populations. We highlight recent studies and methodological advances that aim to overcome such challenges, leveraging genomic signatures of admixture that occurred in the past tens of generations to gain insights into human history, natural selection and complex trait architecture.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

Haplotypes of single cancer driver genes and their local ancestry in a highly admixed long-lived population of Northeast Brazil

Admixed populations have not been examined in detail in cancer genetic studies. Here, we inferred the local ancestry of cancer-associated single nucleotide polymorphisms (SNPs) and haplotypes of a highly admixed Brazilian population. SNP array was used to genotype 73 unrelated individuals aged 80-102 years. Local ancestry inference was performed by merging genotyped regions with phase three data from the 1000 Genomes Project Consortium using RFmix. The average ancestry tract length was 9.12-81.71 megabases. Strong linkage disequilibrium was detected in 48 haplotypes containing 35 SNPs in 10 cancer driver genes. All together, 19 risk and eight protective alleles were identified in 23 out of 48 haplotypes. Homozygous individuals were mainly of European ancestry, whereas heterozygotes had at least one Native American and one African ancestry tract. Native-American ancestry for homozygous individuals with risk alleles for HNF1B, CDH1, and BRCA1 was inferred for the first time. Results indicated that analysis of SNP polymorphism in the present admixed population has a high potential to identify new ancestry-associated alleles and haplotypes that modify cancer susceptibility differentially in distinct human populations. Future case-control studies with populations with a complex history of admixture could help elucidate ancestry-associated biological differences in cancer incidence and therapeutic outcomes.

Admixed Populations Improve Power for Variant Discovery and Portability in Genome-Wide Association Studies

Genome-wide association studies (GWAS) are primarily conducted in single-ancestry settings. The low transferability of results has limited our understanding of human genetic architecture across a range of complex traits. In contrast to homogeneous populations, admixed populations provide an opportunity to capture genetic architecture contributed from multiple source populations and thus improve statistical power. Here, we provide a mechanistic simulation framework to investigate the statistical power and transferability of GWAS under directional polygenic selection or varying divergence. We focus on a two-way admixed population and show that GWAS in admixed populations can be enriched for power in discovery by up to 2-fold compared to the ancestral populations under similar sample size. Moreover, higher accuracy of cross-population polygenic score estimates is also observed if variants and weights are trained in the admixed group rather than in the ancestral groups. Common variant associations are also more likely to replicate if first discovered in the admixed group and then transferred to an ancestral population, than the other way around (across 50 iterations with 1,000 causal SNPs, training on 10,000 individuals, testing on 1,000 in each population, p = 3.78e-6, 6.19e-101, ∼0 for F_ST = 0.2, 0.5, 0.8, respectively). While some of these F_ST values may appear extreme, we demonstrate that they are found across the entire phenome in the GWAS catalog. This framework demonstrates that investigation of admixed populations harbors significant advantages over GWAS in single-ancestry cohorts for uncovering the genetic architecture of traits and will improve downstream applications such as personalized medicine across diverse populations.

Tractor uses local ancestry to enable the inclusion of admixed individuals in GWAS and to boost power

Admixed populations are routinely excluded from genomic studies due to concerns over population structure. Here, we present a statistical framework and software package, Tractor, to facilitate the inclusion of admixed individuals in association studies by leveraging local ancestry. We test Tractor with simulated and empirical two-way admixed African-European cohorts. Tractor generates accurate ancestry-specific effect-size estimates and P values, can boost genome-wide association study (GWAS) power and improves the resolution of association signals. Using a local ancestry-aware regression model, we replicate known hits for blood lipids, discover novel hits missed by standard GWAS and localize signals closer to putative causal variants.

Cohort Profile: ZOE 2.0-A Community-Based Genetic Epidemiologic Study of Early Childhood Oral Health

Early childhood caries (ECC) is an aggressive form of dental caries occurring in the first five years of life. Despite its prevalence and consequences, little progress has been made in its prevention and even less is known about individuals' susceptibility or genomic risk factors. The genome-wide association study (GWAS) of ECC ("ZOE 2.0") is a community-based, multi-ethnic, cross-sectional, genetic epidemiologic study seeking to address this knowledge gap. This paper describes the study's design, the cohort's demographic profile, data domains, and key oral health outcomes. Between 2016 and 2019, the study enrolled 8059 3-5-year-old children attending public preschools in North Carolina, United States. Participants resided in 86 of the state's 100 counties and racial/ethnic minorities predominated-for example, 48% (n = 3872) were African American, 22% white, and 20% (n = 1611) were Hispanic/Latino. Seventy-nine percent (n = 6404) of participants underwent clinical dental examinations yielding ECC outcome measures-ECC (defined at the established caries lesion threshold) prevalence was 54% and the mean number of decayed, missing, filled surfaces due to caries was eight. Nearly all (98%) examined children provided sufficient DNA from saliva for genotyping. The cohort's community-based nature and rich data offer excellent opportunities for addressing important clinical, epidemiologic, and biological questions in early childhood.

Impact of admixture and ancestry on eQTL analysis and GWAS colocalization in GTEx

BACKGROUND

Population structure among study subjects may confound genetic association studies, and lack of proper correction can lead to spurious findings. The Genotype-Tissue Expression (GTEx) project largely contains individuals of European ancestry, but the v8 release also includes up to 15% of individuals of non-European ancestry. Assessing ancestry-based adjustments in GTEx improves portability of this research across populations and further characterizes the impact of population structure on GWAS colocalization.

RESULTS

Here, we identify a subset of 117 individuals in GTEx (v8) with a high degree of population admixture and estimate genome-wide local ancestry. We perform genome-wide cis-eQTL mapping using admixed samples in seven tissues, adjusted by either global or local ancestry. Consistent with previous work, we observe improved power with local ancestry adjustment. At loci where the two adjustments produce different lead variants, we observe 31 loci (0.02%) where a significant colocalization is called only with one eQTL ancestry adjustment method. Notably, both adjustments produce similar numbers of significant colocalizations within each of two different colocalization methods, COLOC and FINEMAP. Finally, we identify a small subset of eQTL-associated variants highly correlated with local ancestry, providing a resource to enhance functional follow-up.

CONCLUSIONS

We provide a local ancestry map for admixed individuals in the GTEx v8 release and describe the impact of ancestry and admixture on gene expression, eQTLs, and GWAS colocalization. While the majority of the results are concordant between local and global ancestry-based adjustments, we identify distinct advantages and disadvantages to each approach.

Mixed-model admixture mapping identifies smoking-dependent loci of lung function in African Americans

Admixture mapping has led to the discovery of many genes associated with differential disease risk by ancestry, highlighting the importance of ancestry-based approaches to association studies. However, the potential of admixture mapping in deciphering the interplay between genes and environment exposures has been seldom explored. Here we performed a genome-wide screening of local ancestry-smoking interactions for five spirometric lung function phenotypes in 3300 African Americans from the COPDGene study. To account for population structure and outcome heterogeneity across exposure groups, we developed a multi-component linear mixed model for mapping gene-environment interactions and empirically showed its robustness and increased power. When applied to the COPDGene study, our approach identified two 11p15.2-3 and 2q37 loci, exhibiting local ancestry-smoking interactions at genome-wide significant level, which would have been missed by standard single-nucleotide polymorphism analyses. These two loci harbor the PARVA and RAB17 genes previously recognized to be involved in smoking behavior. Overall, our study provides the first evidence for potential synergistic effects between African ancestry and smoking on pulmonary function, and underlines the importance of ethnic diversity in genetic studies.

Identifying Rare Variant Associations in Admixed Populations

An admixed population and its ancestral populations bear different burdens of a complex disease. The ancestral populations may have different haplotypes of deleterious alleles and thus ancestry-gene interaction can influence disease risk in the admixed population. Among admixed individuals, deleterious haplotypes and their ancestries are dependent and can provide non-redundant association information. Herein we propose a local ancestry boosted sum test (LABST) for identifying chromosomal blocks that harbor rare variants but have no ancestry switches. For such a stable ancestral block, our LABST exploits ancestry-gene interaction and the number of rare alleles therein. Under the null of no genetic association, the test statistic asymptotically follows a chi-square distribution with one degree of freedom (1-df). Our LABST properly controlled type I error rates under extensive simulations, suggesting that the asymptotic approximation was accurate for the null distribution of the test statistic. In terms of power for identifying rare variant associations, our LABST uniformly outperformed several famed methods under four important modes of disease genetics over a large range of relative risks. In conclusion, exploiting ancestry-gene interaction can boost statistical power for rare variant association mapping in admixed populations.

A robust and powerful two-step testing procedure for local ancestry adjusted allelic association analysis in admixed populations

Genetic association studies in admixed populations allow us to gain deeper understanding of the genetic architecture of human diseases and traits. However, population stratification, complicated linkage disequilibrium (LD) patterns, and the complex interplay of allelic and ancestry effects on phenotypic traits pose challenges in such analyses. These issues may lead to detecting spurious associations and/or result in reduced statistical power. Fortunately, if handled appropriately, these same challenges provide unique opportunities for gene mapping. To address these challenges and to take these opportunities, we propose a robust and powerful two-step testing procedure Local Ancestry Adjusted Allelic (LAAA) association. In the first step, LAAA robustly captures associations due to allelic effect, ancestry effect, and interaction effect, allowing detection of effect heterogeneity across ancestral populations. In the second step, LAAA identifies the source of association, namely allelic, ancestry, or the combination. By jointly modeling allele, local ancestry, and ancestry-specific allelic effects, LAAA is highly powerful in capturing the presence of interaction between ancestry and allele effect. We evaluated the validity and statistical power of LAAA through simulations over a broad spectrum of scenarios. We further illustrated its usefulness by application to the Candidate Gene Association Resource (CARe) African American participants for association with hemoglobin levels. We were able to replicate independent groups' previously identified loci that would have been missed in CARe without joint testing. Moreover, the loci, for which LAAA detected potential effect heterogeneity, were replicated among African Americans from the Women's Health Initiative study. LAAA is freely available at https://yunliweb.its.unc.edu/LAAA.

Genetic Ancestry for Sleep Research: Leveraging Health Inequalities to Identify Causal Genetic Variants

Recent evidence has highlighted the health inequalities in sleep behaviors and sleep disorders that adversely affect outcomes in select populations, including African-American and Hispanic-American subjects. Race-related sleep health inequalities are ascribed to differences in multilevel and interlinked health determinants, such as sociodemographic factors, health behaviors, and biology. African-American and Hispanic-American subjects are admixed populations whose genetic inheritance combines two or more ancestral populations originating from different continents. Racial inequalities in admixed populations can be parsed into relevant groups of mediating factors (environmental vs genetic) with the use of measures of genetic ancestry, including the proportion of an individual's genetic makeup that comes from each of the major ancestral continental populations. This review describes sleep health inequalities in African-American and Hispanic-American subjects and considers the potential utility of ancestry studies to exploit these differences to gain insight into the genetic underpinnings of these phenotypes. The inclusion of genetic approaches in future studies of admixed populations will allow greater understanding of the potential biological basis of race-related sleep health inequalities.

UBD modifies APOL1-induced kidney disease risk

People of recent African ancestry develop kidney disease at much higher rates than most other groups. Two specific coding variants in the Apolipoprotein-L1 gene APOL1 termed G1 and G2 are the causal drivers of much of this difference in risk, following a recessive pattern of inheritance. However, most individuals with a high-risk APOL1 genotype do not develop overt kidney disease, prompting interest in identifying those factors that interact with APOL1 We performed an admixture mapping study to identify genetic modifiers of APOL1-associated kidney disease. Individuals with two APOL1 risk alleles and focal segmental glomerulosclerosis (FSGS) have significantly increased African ancestry at the UBD (also known as FAT10) locus. UBD is a ubiquitin-like protein modifier that targets proteins for proteasomal degradation. African ancestry at the UBD locus correlates with lower levels of UBD expression. In cell-based experiments, the disease-associated APOL1 alleles (known as G1 and G2) lead to increased abundance of UBD mRNA but to decreased levels of UBD protein. UBD gene expression inversely correlates with G1 and G2 APOL1-mediated cell toxicity, as well as with levels of G1 and G2 APOL1 protein in cells. These studies support a model whereby inflammatory stimuli up-regulate both UBD and APOL1, which interact in a functionally important manner. UBD appears to mitigate APOL1-mediated toxicity by targeting it for destruction. Thus, genetically encoded differences in UBD and UBD expression appear to modify the APOL1-associated kidney phenotype.

Critical Issues in the Inclusion of Genetic and Epigenetic Information in Prevention and Intervention Trials

Human genetic research in the past decade has generated a wealth of data from the genome-wide association scan era, much of which is catalogued and freely available. These data will typically test the relationship between a single nucleotide variant or polymorphism (SNP) and some outcome, disease, or trait. Ongoing investigations will yield a similar wealth of data regarding epigenetic phenomena. These data will typically test the relationship between DNA methylation at a single genomic location/region and some outcome. Most of these findings will be the result of cross-sectional investigations typically using ascertained cases and controls. Consequently, most methodological consideration focuses on methods appropriate for simple case-control comparisons. It is expected that a growing number of investigators with longitudinal experimental prevention or intervention cohorts will also measure genetic and epigenetic indicators as part of their investigations, harvesting the wealth of information generated by the genome-wide association study (GWAS) era to allow for targeted hypothesis testing in the next generation of prevention and intervention trials. Herein, we discuss appropriate quality control and statistical modelling of genetic, polygenic, and epigenetic measures in longitudinal models. We specifically discuss quality control, population stratification, genotype imputation, pathway approaches, and proper modelling of an interaction between a specific genetic variant and an environment variable (GxE interaction).

Properties of global- and local-ancestry adjustments in genetic association tests in admixed populations

Population substructure can lead to confounding in tests for genetic association, and failure to adjust properly can result in spurious findings. Here we address this issue of confounding by considering the impact of global ancestry (average ancestry across the genome) and local ancestry (ancestry at a specific chromosomal location) on regression parameters and relative power in ancestry-adjusted and -unadjusted models. We examine theoretical expectations under different scenarios for population substructure; applying different regression models, verifying and generalizing using simulations, and exploring the findings in real-world admixed populations. We show that admixture does not lead to confounding when the trait locus is tested directly in a single admixed population. However, if there is more complex population structure or a marker locus in linkage disequilibrium (LD) with the trait locus is tested, both global and local ancestry can be confounders. Additionally, we show the genotype parameters of adjusted and unadjusted models all provide tests for LD between the marker and trait locus, but in different contexts. The local ancestry adjusted model tests for LD in the ancestral populations, while tests using the unadjusted and the global ancestry adjusted models depend on LD in the admixed population(s), which may be enriched due to different ancestral allele frequencies. Practically, this implies that global-ancestry adjustment should be used for screening, but local-ancestry adjustment may better inform fine mapping and provide better effect estimates at trait loci.

Genome-wide association study of coronary artery calcified atherosclerotic plaque in African Americans with type 2 diabetes

Background

Coronary artery calcified atherosclerotic plaque (CAC) predicts cardiovascular disease (CVD). Despite exposure to more severe conventional CVD risk factors, African Americans (AAs) are less likely to develop CAC, and when they do, have markedly lower levels than European Americans. Genetic factors likely contribute to the observed ethnic differences. To identify genes associated with CAC in AAs with type 2 diabetes (T2D), a genome-wide association study (GWAS) was performed using the Illumina 5 M chip in 691 African American-Diabetes Heart Study participants (AA-DHS), with replication in 205 Jackson Heart Study (JHS) participants with T2D. Genetic association tests were performed on the genotyped and 1000 Genomes-imputed markers separately for each study, and combined in a meta-analysis.

Results

Single nucleotide polymorphisms (SNPs), rs11353135 (2q22.1), rs16879003 (6p22.3), rs5014012, rs58071836 and rs10244825 (all on chromosome 7), rs10918777 (9q31.2), rs13331874 (16p13.3) and rs4459623 (18q12.1) were associated with presence and/or quantity of CAC in the AA-DHS and JHS, with meta-analysis p-values ≤8.0 × 10⁻⁷. The strongest result in AA-DHS alone was rs6491315 in the 13q32.1 region (parameter estimate (SE) = −1.14 (0.20); p-value = 9.1 × 10⁻⁹). This GWAS peak replicated a previously reported AA-DHS CAC admixture signal (rs7492028, LOD score 2.8).

Conclusions

Genetic association between SNPs on chromosomes 2, 6, 7, 9, 16 and 18 and CAC were detected in AAs with T2D from AA-DHS and replicated in the JHS. These data support a role for genetic variation on these chromosomes as contributors to CAC in AAs with T2D, as well as to variation in CAC between populations of African and European ancestry.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0572-9) contains supplementary material, which is available to authorized users.

Genome-wide admixture and association study of subclinical atherosclerosis in the Women's Interagency HIV Study (WIHS)

Cardiovascular disease (CVD) is a major comorbidity among HIV-infected individuals. Common carotid artery intima-media thickness (cCIMT) is a valid and reliable subclinical measure of atherosclerosis and is known to predict CVD. We performed genome-wide association (GWA) and admixture analysis among 682 HIV-positive and 288 HIV-negative Black, non-Hispanic women from the Women’s Interagency HIV study (WIHS) cohort using a combined and stratified analysis approach. We found some suggestive associations but none of the SNPs reached genome-wide statistical significance in our GWAS analysis. The top GWAS SNPs were rs2280828 in the region intergenic to mediator complex subunit 30 and exostosin glycosyltransferase 1 (MED30 | EXT1) among all women, rs2907092 in the catenin delta 2 (CTNND2) gene among HIV-positive women, and rs7529733 in the region intergenic to family with sequence similarity 5, member C and regulator of G-protein signaling 18 (FAM5C | RGS18) genes among HIV-negative women. The most significant local European ancestry associations were in the region intergenic to the zinc finger and SCAN domain containing 5D gene and NADH: ubiquinone oxidoreductase complex assembly factor 1 (ZSCAN5D | NDUF1) pseudogene on chromosome 19 among all women, in the region intergenic to vomeronasal 1 receptor 6 pseudogene and zinc finger protein 845 (VN1R6P | ZNF845) gene on chromosome 19 among HIV-positive women, and in the region intergenic to the SEC23-interacting protein and phosphatidic acid phosphatase type 2 domain containing 1A (SEC23IP | PPAPDC1A) genes located on chromosome 10 among HIV-negative women. A number of previously identified SNP associations with cCIMT were also observed and included rs2572204 in the ryanodine receptor 3 (RYR3) and an admixture region in the secretion-regulating guanine nucleotide exchange factor (SERGEF) gene. We report several SNPs and gene regions in the GWAS and admixture analysis, some of which are common across HIV-positive and HIV-negative women as demonstrated using meta-analysis, and also across the two analytic approaches (i.e., GWA and admixture). These findings suggest that local European ancestry plays an important role in genetic associations of cCIMT among black women from WIHS along with other environmental factors that are related to CVD and may also be triggered by HIV. These findings warrant confirmation in independent samples.

Population Stratification in Genetic Association Studies

Population stratification (PS) is a primary consideration in studies of genetic determinants of human traits. Failure to control for PS may lead to confounding, causing a study to fail for lack of significant results, or resources to be wasted following false-positive signals. Here, historical and current approaches for addressing PS when performing genetic association studies in human populations are reviewed. Methods for detecting the presence of PS, including global and local ancestry methods, are described. Also described are approaches for accounting for PS when calculating association statistics, such that measures of association are not confounded. Many traits are being examined for the first time in minority populations, which may inherently feature PS. © 2017 by John Wiley & Sons, Inc.

Overview of Admixture Mapping

Admixture mapping is a powerful method of gene mapping for diseases or traits that show differential risk by ancestry. Admixture mapping has been applied most often to Americans who trace ancestry to various combinations of Native Americans, Europeans, and West Africans. Recent developments in admixture mapping include improvements in methods and the reference data needed to make inferences about ancestry, as well as extensions of the mapping approach in the framework of linear mixed models. In this unit, the key concepts of admixture mapping are outlined. Several approaches for inferring local ancestry are described, and strategies for performing admixture mapping depending on the study design are provided. Finally, comparisons and contrasts between linkage analysis, association analysis, and admixture mapping are provided, with an emphasis on integrating admixture mapping and association testing. © 2017 by John Wiley & Sons, Inc.

Joint genotype- and ancestry-based genome-wide association studies in admixed populations

In genome-wide association studies (GWAS) genetic loci that influence complex traits are localized by inspecting associations between genotypes of genetic markers and the values of the trait of interest. On the other hand, admixture mapping, which is performed in case of populations consisting of a recent mix of two ancestral groups, relies on the ancestry information at each locus (locus-specific ancestry). Recently it has been proposed to jointly model genotype and locus-specific ancestry within the framework of single marker tests. Here, we extend this approach for population-based GWAS in the direction of multimarker models. A modified version of the Bayesian information criterion is developed for building a multilocus model that accounts for the differential correlation structure due to linkage disequilibrium (LD) and admixture LD. Simulation studies and a real data example illustrate the advantages of this new approach compared to single-marker analysis or modern model selection strategies based on separately analyzing genotype and ancestry data, as well as to single-marker analysis combining genotypic and ancestry information. Depending on the signal strength, our procedure automatically chooses whether genotypic or locus-specific ancestry markers are added to the model. This results in a good compromise between the power to detect causal mutations and the precision of their localization. The proposed method has been implemented in R and is available at http://www.math.uni.wroc.pl/~mbogdan/admixtures/.

Human Demographic History Impacts Genetic Risk Prediction across Diverse Populations

The vast majority of genome-wide association studies (GWASs) are performed in Europeans, and their transferability to other populations is dependent on many factors (e.g., linkage disequilibrium, allele frequencies, genetic architecture). As medical genomics studies become increasingly large and diverse, gaining insights into population history and consequently the transferability of disease risk measurement is critical. Here, we disentangle recent population history in the widely used 1000 Genomes Project reference panel, with an emphasis on populations underrepresented in medical studies. To examine the transferability of single-ancestry GWASs, we used published summary statistics to calculate polygenic risk scores for eight well-studied phenotypes. We identify directional inconsistencies in all scores; for example, height is predicted to decrease with genetic distance from Europeans, despite robust anthropological evidence that West Africans are as tall as Europeans on average. To gain deeper quantitative insights into GWAS transferability, we developed a complex trait coalescent-based simulation framework considering effects of polygenicity, causal allele frequency divergence, and heritability. As expected, correlations between true and inferred risk are typically highest in the population from which summary statistics were derived. We demonstrate that scores inferred from European GWASs are biased by genetic drift in other populations even when choosing the same causal variants and that biases in any direction are possible and unpredictable. This work cautions that summarizing findings from large-scale GWASs may have limited portability to other populations using standard approaches and highlights the need for generalized risk prediction methods and the inclusion of more diverse individuals in medical genomics.

Human Demographic History Impacts Genetic Risk Prediction across Diverse Populations

The vast majority of genome-wide association studies (GWASs) are performed in Europeans, and their transferability to other populations is dependent on many factors (e.g., linkage disequilibrium, allele frequencies, genetic architecture). As medical genomics studies become increasingly large and diverse, gaining insights into population history and consequently the transferability of disease risk measurement is critical. Here, we disentangle recent population history in the widely used 1000 Genomes Project reference panel, with an emphasis on populations underrepresented in medical studies. To examine the transferability of single-ancestry GWASs, we used published summary statistics to calculate polygenic risk scores for eight well-studied phenotypes. We identify directional inconsistencies in all scores; for example, height is predicted to decrease with genetic distance from Europeans, despite robust anthropological evidence that West Africans are as tall as Europeans on average. To gain deeper quantitative insights into GWAS transferability, we developed a complex trait coalescent-based simulation framework considering effects of polygenicity, causal allele frequency divergence, and heritability. As expected, correlations between true and inferred risk are typically highest in the population from which summary statistics were derived. We demonstrate that scores inferred from European GWASs are biased by genetic drift in other populations even when choosing the same causal variants and that biases in any direction are possible and unpredictable. This work cautions that summarizing findings from large-scale GWASs may have limited portability to other populations using standard approaches and highlights the need for generalized risk prediction methods and the inclusion of more diverse individuals in medical genomics.

Evaluating genetic susceptibility to Staphylococcus aureus bacteremia in African Americans using admixture mapping

The incidence of Staphylococcus aureus bacteremia (SAB) is significantly higher in African American (AA) than in European-descended populations. We used admixture mapping (AM) to test the hypothesis that genomic variations with different frequencies in European and African ancestral genomes influence susceptibility to SAB in AAs. A total of 565 adult AAs (390 cases with SAB; 175 age-matched controls) were genotyped for AM analysis. A case-only admixture score and a mixed χ²(1df) score (MIX) to jointly evaluate both single-nucleotide polymorphism (SNP) and admixture association (P<5.00e-08) were computed using MIXSCORE. In addition, a permutation scheme was implemented to derive multiplicity adjusted P-values (genome-wide 0.05 significance threshold: P<9.46e-05). After empirical multiplicity adjustment, one region on chromosome 6 (52 SNPs, P=4.56e-05) in the HLA class II region was found to exhibit a genome-wide statistically significant increase in European ancestry. This region encodes genes involved in HLA-mediated immune response and these results provide additional evidence for genetic variation influencing HLA-mediated immunity, modulating susceptibility to SAB.

Mixed Model Association with Family-Biased Case-Control Ascertainment

Mixed models have become the tool of choice for genetic association studies; however, standard mixed model methods may be poorly calibrated or underpowered under family sampling bias and/or case-control ascertainment. Previously, we introduced a liability threshold-based mixed model association statistic (LTMLM) to address case-control ascertainment in unrelated samples. Here, we consider family-biased case-control ascertainment, where case and control subjects are ascertained non-randomly with respect to family relatedness. Previous work has shown that this type of ascertainment can severely bias heritability estimates; we show here that it also impacts mixed model association statistics. We introduce a family-based association statistic (LT-Fam) that is robust to this problem. Similar to LTMLM, LT-Fam is computed from posterior mean liabilities (PML) under a liability threshold model; however, LT-Fam uses published narrow-sense heritability estimates to avoid the problem of biased heritability estimation, enabling correct calibration. In simulations with family-biased case-control ascertainment, LT-Fam was correctly calibrated (average χ² = 1.00-1.02 for null SNPs), whereas the Armitage trend test (ATT), standard mixed model association (MLM), and case-control retrospective association test (CARAT) were mis-calibrated (e.g., average χ² = 0.50-1.22 for MLM, 0.89-2.65 for CARAT). LT-Fam also attained higher power than other methods in some settings. In 1,259 type 2 diabetes-affected case subjects and 5,765 control subjects from the CARe cohort, downsampled to induce family-biased ascertainment, LT-Fam was correctly calibrated whereas ATT, MLM, and CARAT were again mis-calibrated. Our results highlight the importance of modeling family sampling bias in case-control datasets with related samples.

Genome-wide methylation data mirror ancestry information

Background

Genetic data are known to harbor information about human demographics, and genotyping data are commonly used for capturing ancestry information by leveraging genome-wide differences between populations. In contrast, it is not clear to what extent population structure is captured by whole-genome DNA methylation data.

Results

We demonstrate, using three large-cohort 450K methylation array data sets, that ancestry information signal is mirrored in genome-wide DNA methylation data and that it can be further isolated more effectively by leveraging the correlation structure of CpGs with cis-located SNPs. Based on these insights, we propose a method, EPISTRUCTURE, for the inference of ancestry from methylation data, without the need for genotype data.

Conclusions

EPISTRUCTURE can be used to infer ancestry information of individuals based on their methylation data in the absence of corresponding genetic data. Although genetic data are often collected in epigenetic studies of large cohorts, these are typically not made publicly available, making the application of EPISTRUCTURE especially useful for anyone working on public data. Implementation of EPISTRUCTURE is available in GLINT, our recently released toolset for DNA methylation analysis at: http://glint-epigenetics.readthedocs.io.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0108-y) contains supplementary material, which is available to authorized users.

Genetic architecture of age-related cognitive decline in African Americans

OBJECTIVE

To identify genetic risk factors associated with susceptibility to age-related cognitive decline in African Americans (AAs).

METHODS

We performed a genome-wide association study (GWAS) and an admixture-mapping scan in 3,964 older AAs from 5 longitudinal cohorts; for each participant, we calculated a slope of an individual's global cognitive change from neuropsychological evaluations. We also performed a pathway-based analysis of the age-related cognitive decline GWAS.

RESULTS

We found no evidence to support the existence of a genomic region which has a strongly different contribution to age-related cognitive decline in African and European genomes. Known Alzheimer disease (AD) susceptibility variants in the ABCA7 and MS4A loci do influence this trait in AAs. Of interest, our pathway-based analyses returned statistically significant results highlighting a shared risk from lipid/metabolism and protein tyrosine signaling pathways between cognitive decline and AD, but the role of inflammatory pathways is polarized, being limited to AD susceptibility.

CONCLUSIONS

The genetic architecture of aging-related cognitive in AA individuals is largely similar to that of individuals of European descent. In both populations, we note a surprising lack of enrichment for immune pathways in the genetic risk for cognitive decline, despite strong enrichment of these pathways among genetic risk factors for AD.

Genetic Loci and Novel Discrimination Measures Associated with Blood Pressure Variation in African Americans Living in Tallahassee

Sequencing of the human genome and decades of genetic association and linkage studies have dramatically improved our understanding of the etiology of many diseases. However, the multiple causes of complex diseases are still not well understood, in part because genetic and sociocultural risk factors are not typically investigated concurrently. Hypertension is a leading risk factor for cardiovascular disease and afflicts more African Americans than any other racially defined group in the US. Few genetic loci for hypertension have been replicated across populations, which may reflect population-specific differences in genetic variants and/or inattention to relevant sociocultural factors. Discrimination is a salient sociocultural risk factor for poor health and has been associated with hypertension. Here we use a biocultural approach to study blood pressure (BP) variation in African Americans living in Tallahassee, Florida by genotyping over 30,000 single nucleotide polymorphisms (SNPs) and capturing experiences of discrimination using novel measures of unfair treatment of self and others (n = 157). We perform a joint admixture and genetic association analysis for BP that prioritizes regions of the genome with African ancestry. We only report significant SNPs that were confirmed through our simulation analyses, which were performed to determine the false positive rate. We identify eight significant SNPs in five genes that were previously associated with cardiovascular diseases. When we include measures of unfair treatment and test for interactions between SNPs and unfair treatment, we identify a new class of genes involved in multiple phenotypes including psychosocial distress and mood disorders. Our results suggest that inclusion of culturally relevant stress measures, like unfair treatment in African Americans, may reveal new genes and biological pathways relevant to the etiology of hypertension, and may also improve our understanding of the complexity of gene-environment interactions that underlie complex diseases.

Admixture in Latin America

Latin Americans arguably represent the largest recently admixed populations in the world. This reflects a history of massive settlement by immigrants (mostly Europeans and Africans) and their variable admixture with Natives, starting in 1492. This process resulted in the population of Latin America showing an extensive genetic and phenotypic diversity. Here we review how genetic analyses are being applied to examine the demographic history of this population, including patterns of mating, population structure and ancestry. The admixture history of Latin America, and the resulting extensive diversity of the region, represents a natural experiment offering an advantageous setting for genetic association studies. We review how recent analyses in Latin Americans are contributing to elucidating the genetic architecture of human complex traits.

Assigning breed origin to alleles in crossbred animals

BACKGROUND

For some species, animal production systems are based on the use of crossbreeding to take advantage of the increased performance of crossbred compared to purebred animals. Effects of single nucleotide polymorphisms (SNPs) may differ between purebred and crossbred animals for several reasons: (1) differences in linkage disequilibrium between SNP alleles and a quantitative trait locus; (2) differences in genetic backgrounds (e.g., dominance and epistatic interactions); and (3) differences in environmental conditions, which result in genotype-by-environment interactions. Thus, SNP effects may be breed-specific, which has led to the development of genomic evaluations for crossbred performance that take such effects into account. However, to estimate breed-specific effects, it is necessary to know breed origin of alleles in crossbred animals. Therefore, our aim was to develop an approach for assigning breed origin to alleles of crossbred animals (termed BOA) without information on pedigree and to study its accuracy by considering various factors, including distance between breeds.

RESULTS

The BOA approach consists of: (1) phasing genotypes of purebred and crossbred animals; (2) assigning breed origin to phased haplotypes; and (3) assigning breed origin to alleles of crossbred animals based on a library of assigned haplotypes, the breed composition of crossbred animals, and their SNP genotypes. The accuracy of allele assignments was determined for simulated datasets that include crosses between closely-related, distantly-related and unrelated breeds. Across these scenarios, the percentage of alleles of a crossbred animal that were correctly assigned to their breed origin was greater than 90 %, and increased with increasing distance between breeds, while the percentage of incorrectly assigned alleles was always less than 2 %. For the remaining alleles, i.e. 0 to 10 % of all alleles of a crossbred animal, breed origin could not be assigned.

CONCLUSIONS

The BOA approach accurately assigns breed origin to alleles of crossbred animals, even if their pedigree is not recorded.

Genome-wide scans reveal variants at EDAR predominantly affecting hair straightness in Han Chinese and Uyghur populations

Hair straightness/curliness is one of the most conspicuous features of human variation and is particularly diverse among populations. A recent genome-wide scan found common variants in the Trichohyalin (TCHH) gene that are associated with hair straightness in Europeans, but different genes might affect this phenotype in other populations. By sampling 2899 Han Chinese, we performed the first genome-wide scan of hair straightness in East Asians, and found EDAR (rs3827760) as the predominant gene (P = 4.67 × 10^-16), accounting for 3.66 % of the total variance. The candidate gene approach did not find further significant associations, suggesting that hair straightness may be affected by a large number of genes with subtle effects. Notably, genetic variants associated with hair straightness in Europeans are generally low in frequency in Han Chinese, and vice versa. To evaluate the relative contribution of these variants, we performed a second genome-wide scan in 709 samples from the Uyghur, an admixed population with both eastern and western Eurasian ancestries. In Uyghurs, both EDAR (rs3827760: P = 1.92 × 10^-12) and TCHH (rs11803731: P = 1.46 × 10^-3) are associated with hair straightness, but EDAR (OR 0.415) has a greater effect than TCHH (OR 0.575). We found no significant interaction between EDAR and TCHH (P = 0.645), suggesting that these two genes affect hair straightness through different mechanisms. Furthermore, haplotype analysis indicates that TCHH is not subject to selection. While EDAR is under strong selection in East Asia, it does not appear to be subject to selection after the admixture in Uyghurs. These suggest that hair straightness is unlikely a trait under selection.

Mitochondrial reactive oxygen species and inflammation: Molecular mechanisms, diseases and promising therapies

Over the last few decades, many different groups have been engaged in studies of new roles for mitochondria, particularly the coupling of alterations in the redox pathway with the inflammatory responses involved in different diseases, including Alzheimer's disease, Parkinson's disease, atherosclerosis, cerebral cavernous malformations, cystic fibrosis and cancer. Mitochondrial dysfunction is important in these pathological conditions, suggesting a pivotal role for mitochondria in the coordination of pro-inflammatory signaling from the cytosol and signaling from other subcellular organelles. In this regard, mitochondrial reactive oxygen species are emerging as perfect liaisons that can trigger the assembly and successive activation of large caspase-1- activating complexes known as inflammasomes. This review offers a glimpse into the mechanisms by which inflammasomes are activated by mitochondrial mechanisms, including reactive oxygen species production and mitochondrial Ca²⁺ uptake, and the roles they can play in several inflammatory pathologies.

Admixture mapping of serum vitamin D and parathyroid hormone concentrations in the African American-Diabetes Heart Study

Vitamin D and intact parathyroid hormone (iPTH) concentrations differ between individuals of African and European descent and may play a role in observed racial differences in bone mineral density (BMD). These findings suggest that mapping by admixture linkage disequilibrium (MALD) may be informative for identifying genetic variants contributing to these ethnic disparities. Admixture mapping was performed for serum 25-hydroxyvitamin D, 1,25-dihydroxyvitamin D, vitamin D-binding protein (VDBP), bioavailable vitamin D, and iPTH concentrations and computed tomography measured thoracic and lumbar vertebral volumetric BMD in 552 unrelated African Americans with type 2 diabetes from the African American-Diabetes Heart Study. Genotyping was performed using a custom Illumina ancestry informative marker (AIM) panel. For each AIM, the probability of inheriting 0, 1, or 2 copies of a European-derived allele was determined. Non-parametric linkage analysis was performed by testing for association between each AIM using these probabilities among phenotypes, accounting for global ancestry, age, and gender. Fine-mapping of MALD peaks was facilitated by genome-wide association study (GWAS) data. VDBP levels were significantly linked in proximity to the protein coding locus (rs7689609, LOD=11.05). Two loci exhibited significant linkage signals for 1,25-dihydroxyvitamin D on 13q21.2 (rs1622710, LOD=3.20) and 12q13.2 (rs11171526, LOD=3.10). iPTH was significantly linked on 9q31.3 (rs7854368, LOD=3.14). Fine-mapping with GWAS data revealed significant known (rs7041 with VDBP, P=1.38×10(-82)) and novel (rs12741813 and rs10863774 with VDBP, P<6.43×10(-5)) loci with plausible biological roles. Admixture mapping in combination with fine-mapping has focused efforts to identify loci contributing to ethnic differences in vitamin D-related traits.

Mapping the genomic architecture of adaptive traits with interspecific introgressive origin: a coalescent-based approach

Recent studies of eukaryotes including human and Neandertal, mice, and butterflies have highlighted the major role that interspecific introgression has played in adaptive trait evolution. A common question arises in each case: what is the genomic architecture of the introgressed traits? One common approach that can be used to address this question is association mapping, which looks for genotypic markers that have significant statistical association with a trait. It is well understood that sample relatedness can be a confounding factor in association mapping studies if not properly accounted for. Introgression and other evolutionary processes (e.g., incomplete lineage sorting) typically introduce variation among local genealogies, which can also differ from global sample structure measured across all genomic loci. In contrast, state-of-the-art association mapping methods assume fixed sample relatedness across the genome, which can lead to spurious inference. We therefore propose a new association mapping method called Coal-Map, which uses coalescent-based models to capture local genealogical variation alongside global sample structure. Using simulated and empirical data reflecting a range of evolutionary scenarios, we compare the performance of Coal-Map against EIGENSTRAT, a leading association mapping method in terms of its popularity, power, and type I error control. Our empirical data makes use of hundreds of mouse genomes for which adaptive interspecific introgression has recently been described. We found that Coal-Map's performance is comparable or better than EIGENSTRAT in terms of statistical power and false positive rate. Coal-Map's performance advantage was greatest on model conditions that most closely resembled empirically observed scenarios of adaptive introgression. These conditions had: (1) causal SNPs contained in one or a few introgressed genomic loci and (2) varying rates of gene flow - from high rates to very low rates where incomplete lineage sorting dominated as a primary cause of local genealogical variation.