Genome-wide Association and Population Genetic Analysis of C-Reactive Protein in African American and Hispanic American Women

C-reactive protein (CRP) is a systemic inflammation marker that predicts future cardiovascular risk. CRP levels are higher in African Americans and Hispanic Americans than in European Americans, but the genetic determinants of CRP in these admixed United States minority populations are largely unknown. We performed genome-wide association studies (GWASs) of 8,280 African American (AA) and 3,548 Hispanic American (HA) postmenopausal women from the Women''s Health Initiative SNP Health Association Resource. We discovered and validated a CRP-associated variant of triggering receptors expressed by myeloid cells 2 (TREM2) in chromosomal region 6p21 (p = 10⁻¹⁰). The TREM2 variant associated with higher CRP is common in Africa but rare in other ancestral populations. In AA women, the CRP region in 1q23 contained a strong admixture association signal (p = 10⁻¹⁷), which appears to be related to several independent CRP-associated alleles; the strongest of these is present only in African ancestral populations and is associated with higher CRP. Of the other genomic loci previously associated with CRP through GWASs of European populations, most loci (LEPR, IL1RN, IL6R, GCKR, NLRP3, HNF1A, HNF4A, and APOC1) showed consistent patterns of association with CRP in AA and HA women. In summary, we have identified a common TREM2 variant associated with CRP in United States minority populations. The genetic architecture underlying the CRP phenotype in AA women is complex and involves genetic variants shared across populations, as well as variants specific to populations of African descent.     

Genetic determinants of lipid traits in diverse populations from the population architecture using genomics and epidemiology (PAGE) study

For the past five years, genome-wide association studies (GWAS) have identified hundreds of common variants associated with human diseases and traits, including high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglyceride (TG) levels. Approximately 95 loci associated with lipid levels have been identified primarily among populations of European ancestry. The Population Architecture using Genomics and Epidemiology (PAGE) study was established in 2008 to characterize GWAS-identified variants in diverse population-based studies. We genotyped 49 GWAS-identified SNPs associated with one or more lipid traits in at least two PAGE studies and across six racial/ethnic groups. We performed a meta-analysis testing for SNP associations with fasting HDL-C, LDL-C, and ln(TG) levels in self-identified European American (~20,000), African American (~9,000), American Indian (~6,000), Mexican American/Hispanic (~2,500), Japanese/East Asian (~690), and Pacific Islander/Native Hawaiian (~175) adults, regardless of lipid-lowering medication use. We replicated 55 of 60 (92%) SNP associations tested in European Americans at p<0.05. Despite sufficient power, we were unable to replicate ABCA1 rs4149268 and rs1883025, CETP rs1864163, and TTC39B rs471364 previously associated with HDL-C and MAFB rs6102059 previously associated with LDL-C. Based on significance (p<0.05) and consistent direction of effect, a majority of replicated genotype-phentoype associations for HDL-C, LDL-C, and ln(TG) in European Americans generalized to African Americans (48%, 61%, and 57%), American Indians (45%, 64%, and 77%), and Mexican Americans/Hispanics (57%, 56%, and 86%). Overall, 16 associations generalized across all three populations. For the associations that did not generalize, differences in effect sizes, allele frequencies, and linkage disequilibrium offer clues to the next generation of association studies for these traits.     

Leveraging genetic ancestry to study health disparities

Research to understand human genomic variation and its implications in health has great potential to contribute in the reduction of health disparities. Biological anthropology can play important roles in genomics and health disparities research using a biocultural approach. This paper argues that racial/ethnic categories should not be used as a surrogate for sociocultural factors or global genomic clusters in biomedical research or clinical settings, because of the high genetic heterogeneity that exists within traditional racial/ethnic groups. Genetic ancestry is used to show variation in ancestral genomic contributions to recently admixed populations in the United States, such as African Americans and Hispanic/Latino Americans. Genetic ancestry estimates are also used to examine the relationship between ancestry-related biological and sociocultural factors affecting health disparities. To localize areas of genomes that contribute to health disparities, admixture mapping and genome-wide association studies (GWAS) are often used. Recent GWAS have identified many genetic variants that are highly differentiated among human populations that are associated with disease risk. Some of these are population-specific variants. Many of these variants may impact disease risk and help explain a portion of the difference in disease burden among racial/ethnic groups. Genetic ancestry is also of particular interest in precision medicine and disparities in drug efficacy and outcomes. By using genetic ancestry, we can learn about potential biological differences that may contribute to the heterogeneity observed across self-reported racial groups.     

Genetic background of patients from a university medical center in Manhattan: implications for personalized medicine

Background

The rapid progress currently being made in genomic science has created interest in potential clinical applications; however, formal translational research has been limited thus far. Studies of population genetics have demonstrated substantial variation in allele frequencies and haplotype structure at loci of medical relevance and the genetic background of patient cohorts may often be complex.

Methods and Findings

To describe the heterogeneity in an unselected clinical sample we used the Affymetrix 6.0 gene array chip to genotype self-identified European Americans (N = 326), African Americans (N = 324) and Hispanics (N = 327) from the medical practice of Mount Sinai Medical Center in Manhattan, NY. Additional data from US minority groups and Brazil were used for external comparison. Substantial variation in ancestral origin was observed for both African Americans and Hispanics; data from the latter group overlapped with both Mexican Americans and Brazilians in the external data sets. A pooled analysis of the African Americans and Hispanics from NY demonstrated a broad continuum of ancestral origin making classification by race/ethnicity uninformative. Selected loci harboring variants associated with medical traits and drug response confirmed substantial within- and between-group heterogeneity.

Conclusion

As a consequence of these complementary levels of heterogeneity group labels offered no guidance at the individual level. These findings demonstrate the complexity involved in clinical translation of the results from genome-wide association studies and suggest that in the genomic era conventional racial/ethnic labels are of little value.     

Extent of linkage disequilibrium between the androgen receptor gene CAG and GGC repeats in human populations: implications for prostate cancer risk

While studies have implicated alleles at the CAG and GGC trinucleotide repeats of the androgen receptor gene with high-grade, aggressive prostate cancer disease, little is known about the normal range of variation for these two loci, which are separated by about 1.1 kb. More importantly, few data exist on the extent of linkage disequilibrium (LD) between the two loci in different human populations. Here we present data on CAG and GGC allelic variation and LD in six diverse populations. Alleles at the CAG and GGC repeat loci of the androgen receptor were typed in over 1000 chromosomes from Africa, Asia, and North America. Levels of linkage disequilibrium between the two loci were compared between populations. Haplotype variation and diversity were estimated for each population. Our results reveal that populations of African descent possess significantly shorter alleles for the two loci than non-African populations (P<0.0001). Allelic diversity for both markers was higher among African Americans than any other population, including indigenous Africans from Sierra Leone and Nigeria. Analysis of molecular variance revealed that approx. 20% of CAG and GGC repeat variance could be attributed to differences between the populations. All non-African populations possessed the same common haplotype while the three populations of African descent possessed three divergent common haplotypes. Significant LD was observed in our sample of healthy African Americans. The LD observed in the African American population may be due to several reasons; recent migration of African Americans from diverse rural communities following urbanization, recurrent gene flow from diverse West African populations, and admixture with European Americans. This study represents the largest genotyping effort to be performed on the two androgen receptor trinucleotide repeat loci in diverse human populations.     

Genetic meta-analysis of 15,901 African Americans identifies variation in EXOC3L1 is associated with HDL concentration

Meta-analyses of European populations has successfully identified genetic variants in over 150 loci associated with lipid levels, but results from additional ethnicities remain limited. Previously, we reported two novel lipid loci identified in a sample of 7,657 African Americans using a gene-centric array including 50,000 SNPs in 2,100 candidate genes. Initial discovery and follow-up of signals with P < 10⁻⁵ in additional African American samples confirmed CD36 and ICAM1. Using an additional 8,244 African American female samples from the Women’s Health Initiative SNP Health Association Resource genome-wide association study dataset, we further examined the previous meta-analyses results by attempting to replicate 20 additional putative lipid signals with P < 10⁻⁴. Replication confirmed rs868213, located in a splice donor region of exocyst complex component 3-like 1 (EXOC3L1) as a novel signal for HDL (additive allelic effect β = 0.02; P = 1.4 × 10⁻⁸; meta-analyses of discovery and replication). EXOC3L1 is strongly expressed in vascular endothelium and forms part of the exocyst complex, a key facilitator of the trafficking of lipid receptors. Increasing sample sizes for genetic studies in nonEuropean populations will continue to improve our understanding of lipid metabolism.     

β2-Adrenergic receptor polymorphisms: pharmacogenetic response to bronchodilator among African American asthmatics

Beta2-adrenergic receptor (beta2AR) gene polymorphisms have been reported to be associated with various asthma-related traits in different racial/ethnic populations. However, it is unknown whether beta2AR genetic variants are associated with asthma in African Americans. In this study, we have examined whether there is association between beta2AR genetic variants and asthma in African Americans. We have recruited 264 African American asthmatic subjects and 176 matched healthy controls participating in the Study of African Americans, Asthma, Genes and Environments (SAGE). We genotyped seven known and recently identified beta2AR SNP variants, then tested genotype and haplotype association of asthma-related traits with the beta2AR SNPs in our African American cohort with adjustment of confounding effect due to admixture background and environmental risk factors. We found a significant association of the SNP -47 (Arg-19Cys) polymorphism with DeltaFEF(25-75), a measure of bronchodilator drug responsiveness, in African American asthmatics after correction for multiple testing (P = 0.001). We did not observe association of the SNP +46 (Arg16Gly) variant with asthma disease diagnosis and asthma-related phenotypes. In contrast to previous results between the Arg16Gly variant and traits related to bronchodilator responsiveness, our results indicate that the Arg-19Cys polymorphism in beta upstream peptide may play an important role in bronchodilator drug responsiveness in African American subjects. Our findings highlight the importance of investigating genetic risk factors for asthma in different populations.     

Generalization and Dilution of Association Results from European GWAS in Populations of Non-European Ancestry: The PAGE Study

The vast majority of genome-wide association study (GWAS) findings reported to date are from populations with European Ancestry (EA), and it is not yet clear how broadly the genetic associations described will generalize to populations of diverse ancestry. The Population Architecture Using Genomics and Epidemiology (PAGE) study is a consortium of multi-ancestry, population-based studies formed with the objective of refining our understanding of the genetic architecture of common traits emerging from GWAS. In the present analysis of five common diseases and traits, including body mass index, type 2 diabetes, and lipid levels, we compare direction and magnitude of effects for GWAS-identified variants in multiple non-EA populations against EA findings. We demonstrate that, in all populations analyzed, a significant majority of GWAS-identified variants have allelic associations in the same direction as in EA, with none showing a statistically significant effect in the opposite direction, after adjustment for multiple testing. However, 25% of tagSNPs identified in EA GWAS have significantly different effect sizes in at least one non-EA population, and these differential effects were most frequent in African Americans where all differential effects were diluted toward the null. We demonstrate that differential LD between tagSNPs and functional variants within populations contributes significantly to dilute effect sizes in this population. Although most variants identified from GWAS in EA populations generalize to all non-EA populations assessed, genetic models derived from GWAS findings in EA may generate spurious results in non-EA populations due to differential effect sizes. Regardless of the origin of the differential effects, caution should be exercised in applying any genetic risk prediction model based on tagSNPs outside of the ancestry group in which it was derived. Models based directly on functional variation may generalize more robustly, but the identification of functional variants remains challenging.     

Genome-wide association study of coronary heart disease and its risk factors in 8,090 African Americans: the NHLBI CARe Project

Coronary heart disease (CHD) is the leading cause of mortality in African Americans. To identify common genetic polymorphisms associated with CHD and its risk factors (LDL- and HDL-cholesterol (LDL-C and HDL-C), hypertension, smoking, and type-2 diabetes) in individuals of African ancestry, we performed a genome-wide association study (GWAS) in 8,090 African Americans from five population-based cohorts. We replicated 17 loci previously associated with CHD or its risk factors in Caucasians. For five of these regions (CHD: CDKN2A/CDKN2B; HDL-C: FADS1-3, PLTP, LPL, and ABCA1), we could leverage the distinct linkage disequilibrium (LD) patterns in African Americans to identify DNA polymorphisms more strongly associated with the phenotypes than the previously reported index SNPs found in Caucasian populations. We also developed a new approach for association testing in admixed populations that uses allelic and local ancestry variation. Using this method, we discovered several loci that would have been missed using the basic allelic and global ancestry information only. Our conclusions suggest that no major loci uniquely explain the high prevalence of CHD in African Americans. Our project has developed resources and methods that address both admixture- and SNP-association to maximize power for genetic discovery in even larger African-American consortia.     

Genome scans of facial features in East Africans and cross-population comparisons reveal novel associations

Facial morphology is highly variable, both within and among human populations, and a sizable portion of this variation is attributable to genetics. Previous genome scans have revealed more than 100 genetic loci associated with different aspects of normal-range facial variation. Most of these loci have been detected in Europeans, with few studies focusing on other ancestral groups. Consequently, the degree to which facial traits share a common genetic basis across diverse sets of humans remains largely unknown. We therefore investigated the genetic basis of facial morphology in an East African cohort. We applied an open-ended data-driven phenotyping approach to a sample of 2,595 3D facial images collected on Tanzanian children. This approach segments the face into hierarchically arranged, multivariate features that capture the shape variation after adjusting for age, sex, height, weight, facial size and population stratification. Genome scans of these multivariate shape phenotypes revealed significant (p < 2.5 × 10⁻⁸) signals at 20 loci, which were enriched for active chromatin elements in human cranial neural crest cells and embryonic craniofacial tissue, consistent with an early developmental origin of the facial variation. Two of these associations were in highly conserved regions showing craniofacial-specific enhancer activity during embryological development (5q31.1 and 12q21.31). Six of the 20 loci surpassed a stricter threshold accounting for multiple phenotypes with study-wide significance (p < 6.25 × 10⁻¹⁰). Cross-population comparisons indicated 10 association signals were shared with Europeans (seven sharing the same associated SNP), and facilitated fine-mapping of causal variants at previously reported loci. Taken together, these results may point to both shared and population-specific components to the genetic architecture of facial variation.     

Genetic risk factors for BMI and obesity in an ethnically diverse population: Results from the population architecture using genomics and epidemiology (PAGE) study

Objective:

Several genome–wide association studies (GWAS) have demonstrated that common genetic variants contribute to obesity. However, studies of this complex trait have focused on ancestrally European populations, despite the high prevalence of obesity in some minority groups.

Design and Methods:

As part of the “Population Architecture using Genomics and Epidemiology (PAGE)” Consortium, we investigated the association between 13 GWAS‐identified single‐nucleotide polymorphisms (SNPs) and BMI and obesity in 69,775 subjects, including 6,149 American Indians, 15,415 African‐Americans, 2,438 East Asians, 7,346 Hispanics, 604 Pacific Islanders, and 37,823 European Americans. For the BMI‐increasing allele of each SNP, we calculated β coefficients using linear regression (for BMI) and risk estimates using logistic regression (for obesity defined as BMI ≥ 30) followed by fixed‐effects meta‐analysis to combine results across PAGE sites. Analyses stratified by racial/ethnic group assumed an additive genetic model and were adjusted for age, sex, and current smoking. We defined “replicating SNPs” (in European Americans) and “generalizing SNPs” (in other racial/ethnic groups) as those associated with an allele frequency‐specific increase in BMI.

Results:

By this definition, we replicated 9/13 SNP associations (5 out of 8 loci) in European Americans. We also generalized 8/13 SNP associations (5/8 loci) in East Asians, 7/13 (5/8 loci) in African Americans, 6/13 (4/8 loci) in Hispanics, 5/8 in Pacific Islanders (5/8 loci), and 5/9 (4/8 loci) in American Indians.

Conclusion:

Linkage disequilibrium patterns suggest that tagSNPs selected for European Americans may not adequately tag causal variants in other ancestry groups. Accordingly, fine‐mapping in large samples is needed to comprehensively explore these loci in diverse populations.     

Genetic architecture of lipid traits changes over time and differs by race: Princeton Lipid Follow-up Study

Dyslipidemia is a major risk factor for CVD. Previous studies on lipid heritability have largely focused on white populations assessed after the obesity epidemic. Given secular trends and racial differences in lipid levels, this study explored whether lipid heritability is consistent across time and between races. African American and white nuclear families had fasting lipids measured in the 1970s and 22–30 years later. Heritability was estimated, and bivariate analyses between visits were conducted by race using variance components analysis. A total of 1,454 individuals (age 14.1/40.6 for offspring/parents at baseline; 39.6/66.5 at follow-up) in 373 families (286 white, 87 African American) were included. Lipid trait heritabilities were typically stronger during the 1970s than the 2000s. At baseline, additive genetic variation for LDL was significantly lower in African Americans than whites (P = 0.015). Shared genetic contribution to lipid variability over time was significant in both whites (all P < 0.0001) and African Americans (P ≤ 0.05 for total, LDL, and HDL cholesterol). African American families demonstrated shared environmental contributions to lipid variation over time (all P ≤ 0.05). Lower heritability, lower LDL genetic variance, and durable environmental effects across the obesity epidemic in African American families suggest race-specific approaches are needed to clarify the genetic etiology of lipids.     

Genetic variation at twentythree microsatellite loci in sixteen human populations

We have analysed genetic variation at 23 microsatellite loci in a global sample of 16 ethnically and geographically diverse human populations. On the basis of their ancestral heritage and geographic locations, the studied populations can be divided into five major groups, viz. African, Caucasian, Asian Mongoloid, American Indian and Pacific Islander. With respect to the distribution of alleles at the 23 loci, large variability exists among the examined populations. However, with the exception of the American Indians and the Pacific Islanders, populations within a continental group show a greater degree of similarity. Phylogenetic analyses based on allele frequencies at the examined loci show that the first split of the present-day human populations had occurred between the Africans and all of the non-African populations, lending support to an African origin of modern human populations. Gene diversity analyses show that the coefficient of gene diversity estimated from the 23 loci is, in general, larger for populations that have remained isolated and probably of smaller effective sizes, such as the American Indians and the Pacific Islanders. These analyses also demonstrate that the component of total gene diversity, which is attributed to variation between groups of populations, is significantly larger than that among populations within each group. The empirical data presented in this work and their analyses reaffirm that evolutionary histories and the extent of genetic variation among human populations can be studied using microsatellite loci.     

Self-reported ethnicity, genetic structure and the impact of population stratification in a multiethnic study

It is well-known that population substructure may lead to confounding in case–control association studies. Here, we examined genetic structure in a large racially and ethnically diverse sample consisting of five ethnic groups of the Multiethnic Cohort study (African Americans, Japanese Americans, Latinos, European Americans and Native Hawaiians) using 2,509 SNPs distributed across the genome. Principal component analysis on 6,213 study participants, 18 Native Americans and 11 HapMap III populations revealed four important principal components (PCs): the first two separated Asians, Europeans and Africans, and the third and fourth corresponded to Native American and Native Hawaiian (Polynesian) ancestry, respectively. Individual ethnic composition derived from self-reported parental information matched well to genetic ancestry for Japanese and European Americans. STRUCTURE-estimated individual ancestral proportions for African Americans and Latinos are consistent with previous reports. We quantified the East Asian (mean 27%), European (mean 27%) and Polynesian (mean 46%) ancestral proportions for the first time, to our knowledge, for Native Hawaiians. Simulations based on realistic settings of case–control studies nested in the Multiethnic Cohort found that the effect of population stratification was modest and readily corrected by adjusting for race/ethnicity or by adjusting for top PCs derived from all SNPs or from ancestry informative markers; the power of these approaches was similar when averaged across causal variants simulated based on allele frequencies of the 2,509 genotyped markers. The bias may be large in case-only analysis of gene by gene interactions but it can be corrected by top PCs derived from all SNPs.     

Racial admixture and its impact on BMI and blood pressure in African and Mexican Americans

Admixed populations such as African Americans and Hispanic Americans present both challenges and opportunities in genetic epidemiologic research. Because of variation in admixture levels among individuals, case-control association studies may be subject to stratification bias. On the other hand, admixed populations also present special opportunities both for examining the role of genetic and environmental factors for observed racial/ethnic differences, and for possibly mapping alleles that contribute to such differences. Here we examined the distribution and relationship of individual admixture (IA) estimates with BMI and three measures of blood pressure in two admixed populations in the NHLBI Family Blood Pressure Program (FBPP): African Americans and Mexican Americans. For the African Americans, we observed modest but significant differences in average African IA among four recruitment sites. We observed a slight excess of African IA among hypertensives compared to normotensives, and a positive (non-significant) regression of African IA on blood pressure in untreated participants. Within Mexican Americans, we found no difference in average IA between hypertensives and normotensives, but a positive (marginally significant) regression of African IA on diastolic blood pressure. We also observed a significant positive regression of Caucasian IA (and negative regression of Native American IA) on BMI. Our results are suggestive of genetic differences between Africans and non-Africans that influence blood pressure, but such effects are likely to be modest compared to environmental ones. Excess obesity among Native Americans compared to whites is not consistent with a simple genetic explanation.     

Analysis of FTO gene variants with measures of obesity and glucose homeostasis in the IRAS Family Study

Multiple studies have identified FTO gene variants associated with measures of adiposity in European-derived populations. The objective of the study was to determine whether FTO variants were associated with adiposity, including visceral and subcutaneous adipose tissue (VAT, SAT), and glucose homeostasis measures in the Insulin Resistance Atherosclerosis Family Study (IRASFS). A total of 27 SNPs in FTO intron 1, including SNPs prominent in the literature (rs9939609, rs8050136, rs1121980, rs17817449, rs1421085, and rs3751812), were genotyped in 1,424 Hispanic Americans and 604 African Americans. Multiple SNPs were associated with BMI and SAT (P values ranging from 0.001 to 0.033), and trending or associated with waist circumference (P values ranging from 0.008 to 0.099) in the Hispanic Americans. No association was observed with VAT, illustrating that FTO variants are associated with overall fat mass instead of specific fat depots. For the glucose homeostasis measures, variants were associated with fasting insulin but, consistent with other studies, after BMI adjustment, no evidence of association remained. The lack of association of FTO SNPs with insulin sensitivity is consistent with the lack of association with VAT, since these traits are strongly correlated. In the African Americans, only rs8050136 and rs9939609 were associated with BMI and WAIST (P values of 0.011 and 0.034), and associated or trending towards association with SAT (P values of 0.038 and 0.058). These results confirm that FTO variants are associated with adiposity measures, predisposing individuals to obesity by increasing overall fat mass in Hispanic Americans and to a lesser degree in African Americans.     

New approaches to disease mapping in admixed populations

Admixed populations such as African Americans and Hispanic Americans are often medically underserved and bear a disproportionately high burden of disease. Owing to the diversity of their genomes, these populations have both advantages and disadvantages for genetic studies of complex phenotypes. Advances in statistical methodologies that can infer genetic contributions from ancestral populations may yield new insights into the aetiology of disease and may contribute to the applicability of genomic medicine to these admixed population groups.     

Leveraging Multi-ethnic Evidence for Mapping Complex Traits in Minority Populations: An Empirical Bayes Approach

Elucidating the genetic basis of complex traits and diseases in non-European populations is particularly challenging because US minority populations have been under-represented in genetic association studies. We developed an empirical Bayes approach named XPEB (cross-population empirical Bayes), designed to improve the power for mapping complex-trait-associated loci in a minority population by exploiting information from genome-wide association studies (GWASs) from another ethnic population. Taking as input summary statistics from two GWASs—a target GWAS from an ethnic minority population of primary interest and an auxiliary base GWAS (such as a larger GWAS in Europeans)—our XPEB approach reprioritizes SNPs in the target population to compute local false-discovery rates. We demonstrated, through simulations, that whenever the base GWAS harbors relevant information, XPEB gains efficiency. Moreover, XPEB has the ability to discard irrelevant auxiliary information, providing a safeguard against inflated false-discovery rates due to genetic heterogeneity between populations. Applied to a blood-lipids study in African Americans, XPEB more than quadrupled the discoveries from the conventional approach, which used a target GWAS alone, bringing the number of significant loci from 14 to 65. Thus, XPEB offers a flexible framework for mapping complex traits in minority populations.     

Testing assumptions of the Gilbert and Gill method for assessing ancestry using the femur subtrochanteric shape

In 1990, Gilbert and Gill proposed a simple metric technique using femoral subtrochanteric anteroposterior and mediolateral diaphyseal diameters for discriminating between Native American and American Black and White femora in medicolegal and bioarchaeological contexts. However, there are several inherent assumptions in the method that may affect its validity. The assumptions include minimal sexual dimorphism, temporal and geographical homogeneity within populations, and that differences between populations in femoral subtrochanteric size and shape are primarily due to genetic variation. In this study, these assumptions are tested using femora from seven populations (African, American Black, American White, Australian, Native American, Hispanic, and Polynesian). The results indicate that sexual dimorphism and geographical and temporal heterogeneity in proximal femur diaphyseal shape within Native Americans are not great enough to significantly affect the validity of the Gilbert and Gill method (GGM). Variation between populations is most likely due to combined genetic and environmental factors, but differences in proximal femur shape between Native Americans and American Blacks/Whites are sufficient to allow accurate discrimination between these groups. Caution, however, must be taken during investigations where populations other than Native Americans or American Blacks/Whites are present, and therefore, the GGM may have limited forensic anthropological application in many parts of the world.