Impact Factor 3.789

Frontiers reaches 6.4 on Journal Impact Factors

Original Research ARTICLE

Front. Genet., 01 June 2016 | https://doi.org/10.3389/fgene.2016.00089

The Genetic Contribution of West-African Ancestry to Protection against Central Obesity in African-American Men but Not Women: Results from the ARIC and MESA Studies

Yann C. Klimentidis1*, Amit Arora1, Jin Zhou1, Rick Kittles2 and David B. Allison3,4
  • 1Department of Epidemiology and Biostatistics, Mel and Enid Zuckerman College of Public Health, University of Arizona, Tucson, AZ, USA
  • 2Department of Surgery, College of Medicine, University of Arizona, Tucson, AZ, USA
  • 3Nutrition and Obesity Research Center, University of Alabama at Birmingham, Birmingham, AL, USA
  • 4Office of Energetics, University of Alabama at Birmingham, Birmingham, AL, USA

Over 80% of African-American (AA) women are overweight or obese. A large racial disparity between AA and European-Americans (EA) in obesity rates exists among women, but curiously not among men. Although socio-economic and/or cultural factors may partly account for this race-by-sex interaction, the potential involvement of genetic factors has not yet been investigated. Among 2814 self-identified AA in the Atherosclerosis Risk in Communities study, we estimated each individual's degree of West-African genetic ancestry using 3437 ancestry informative markers. We then tested whether sex modifies the association between West-African genetic ancestry and body mass index (BMI), waist-circumference (WC), and waist-to-hip ratio (WHR), adjusting for income and education levels, and examined associations of ancestry with the phenotypes separately in males and females. We replicated our findings in the Multi-Ethnic Study of Atherosclerosis (n = 1611 AA). In both studies, we find that West-African ancestry is negatively associated with obesity, especially central obesity, among AA men, but not among AA women (pinteraction = 4.14 × 10−5 in pooled analysis of WHR). In conclusion, our results suggest that the combination of male gender and West-African genetic ancestry is associated with protection against central adiposity, and suggest that the large racial disparity that exists among women, but not men, may be at least partly attributed to genetic factors.

Introduction

In the United States, African Americans (AA) exhibit higher rates of overweight and obesity. However, this trend appears to be mainly driven by very high rates among AA women. According to data from the Centers for Disease Control and the National Center for Health Statistics, the proportions of overweight or obese AA and EA women are 82.2 and 60.9%, respectively, whereas for men, the proportions are 70.7 and 73.2%, respectively (CDC/NCHS, 2013). A similar trend exists for rates of type-2 diabetes, whereby the race disparity is larger among women than among men (Brancati et al., 2000; Signorello et al., 2007). This puzzling pattern of a greater racial disparity only among women for obesity has not yet been explored from a genetic perspective.

Previous studies have shown that West-African ancestry (WAA) is associated with a greater body mass index (BMI) (Fernandez et al., 2003; Tang et al., 2006; Halder et al., 2012; Nassir et al., 2012), and greater whole-body and trunk fat mass (Ochs-Balcom et al., 2013). On the other hand, several studies have shown that WAA and AA race are associated with a lower body fat percentage, including visceral fat (Yanovski et al., 1996; Gower et al., 1999; Beasley et al., 2009; Cardel et al., 2011). However, some of these studies examined only women, or only BMI, or did not examine sex differences in the association between WAA and obesity phenotypes. Therefore, little is known as to whether this pattern differs among men and women, potentially resulting in the striking racial disparity observed among women, but not among men. A recent study by Goonesekera et al. (2015) reported a stronger (positive) association between WAA and BMI, percent body fat, and waist circumference (WC) among women compared to men in a cohort study of Boston residents, and a stronger positive association between WAA and BMI among less educated women. These findings are consistent with the epidemiological pattern described above, and deserving of further examination.

Based on these previous studies, it is possible that WAA, representing a collection of many individual genetic variants, predisposes to less central/visceral adiposity, and that sex modifies this relationship through biological and/or environmental factors. To help further understand if genetic factors may play a role in the large racial disparity among women, and the relatively small racial disparity among men, we examined the interaction of WAA with sex on BMI, WC, and waist-to-hip ratio (WHR) in two cohorts of AA. We also examined whether the association of WAA with these phenotypes differs according to income, education, or menopausal status.

Methods

Studies

We used data from the Atherosclerosis Risk in Communities (ARIC) study (The Aric Investigators, 1989) and from the Multi-Ethnic Study of Atherosclerosis (MESA) study (Bild et al., 2002) (n = 2814 and 1611 AA, respectively). Briefly, ARIC is a prospective multi-center cohort study of EA and AA, aged 45–64, to examine risk factors for atherosclerosis. The MESA study is also a prospective multi-center cohort study of men and women, aged 45–85, of diverse ethnic and racial backgrounds, to investigate the risk factors for cardiovascular disease. We only included self-identified AA in this study. Datasets were obtained from the database of Genotypes and Phenotypes (dbGaP). We obtained IRB approval from the University of Arizona.

Genotypes and Ancestry Estimation

Both ARIC and MESA participants were genotyped with the Affymetrix Genome-Wide Human SNP Array 6.0 (Affymetrix, Santa Clara, CA, USA). To estimate WAA, we used a panel of 3437 SNPs developed by Tandon et al. which are especially informative for West African vs. European ancestry among AA, are spaced throughout the genome, are not in linkage disequilibrium with each other in the ancestral populations, and are available on the Affymetrix 6.0 array (Tandon et al., 2011). With the genotypes of these SNPs, we performed a principal components analysis (PCA) on the scaled and centered genotypes, using the prcomp function in R Statistical Software (R Development Core Team, 2011). To confirm that PC1 of this analysis corresponds to WAA, we used a supervised Bayesian approach with STRUCTURE v2.3 software (Pritchard et al., 2000) with genotypes at the same set of SNPs for HapMap reference individuals (Yoruba in Ibadan, Nigeria: n = 176, and Utah residents with ancestry from northern and western Europe: n = 174) (The International HapMap Consortium, 2003) as putative W. African and European parental populations. We used K = 2, representing two founding populations, and a burn-in length of 20,000 for 80,000 repetitions. There was a very strong relationship between PC1 and WAA as determined by the supervised ancestry analysis (r2≥ 0.96).

Phenotypes, Income, Education, and Menopausal Status

BMI was calculated as weight (kg) divided by squared height (m2). Waist measurements were taken at the umbilicus to the nearest 0.1 cm in MESA, and nearest 1 cm in ARIC. Hip circumference was measured from the largest diameter of the hip, to the nearest 0.1 cm in MESA, and nearest 1 cm in ARIC. WHR was calculated by dividing waist circumference by hip circumference. Income was recoded to correspond to yearly household/family income in US dollars. As yearly income intervals were presented in the questionnaires, we chose midpoints of these intervals. In both datasets, education was recoded to correspond to years of education (using interval midpoints, if needed). Menopausal status was assessed as a binary variable in both cohorts using the response to a question asking if women had gone through menopause or not.

Statistical Analyses

Differences in WAA, phenotypes, and other characteristics were compared across cohorts and in men vs. women using t-tests. The difference in percentage of women having undergone menopause between the two cohorts was tested using a chi-squared test. We examined the distributions of these three phenotypes as well as the distributions of WAA, income, and education to assess normality. We then used linear regression to test the association of WAA with each phenotype among all, male, and female AA, in each study and in the combined dataset. We log transformed BMI and WC to conform to the assumptions of residual normality. We included age, sex, income, education, and WAA as covariates. In the combined analyses, we additionally included a two-level categorical “study” variable (ARIC/MESA) in order to account for variation across studies. The interaction of WAA with sex was modeled as the product of these two variables and included as a predictor in the linear regression model. To test whether the association of WAA with the obesity-related phenotypes differs according to income, education or menopausal status, we tested the interaction of WAA with each phenotype separately in men and/or women, including in the model all other covariates. Statistical analyses were carried out with R Statistical software, version 2.15.2 (R Development Core Team, 2011).

Results

Association of WAA with Obesity- and Metabolic-Related Phenotypes

Characteristics of AA participants in MESA and ARIC are shown in Table 1. In addition to mean age being higher in MESA, levels of income, education, BMI, and WC are also higher in MESA. The distribution of WAA in both ARIC and MESA are shown in Figure 1. As shown in Table 1, WAA is significantly higher among ARIC participants (88%) than among MESA participants (80%). However, we do not find any difference in WAA between AA men and women in either MESA or ARIC (p = 0.81 and 0.93, respectively, see Figure 1). The main effect models for ARIC and MESA and for each of the three phenotypes are shown in Table 2. Notably, among AA in ARIC, WAA is positively and significantly associated with BMI in ARIC AA (p = 0.004), but not MESA AA (see Table 2). In general, we observe a trend of a positive association between WAA and BMI, no consistent association for WC, and a negative association of WAA with WHR. We also observe that both income and education are negatively associated with all three phenotypes, and these associations are much stronger in women, and nearly absent in men (data not shown).

TABLE 1
www.frontiersin.org

Table 1. Characteristics of male and female African-American participants in ARIC and MESA.

FIGURE 1
www.frontiersin.org

Figure 1. Histograms comparing the distribution of WAA among men and women in ARIC and MESA.

TABLE 2
www.frontiersin.org

Table 2. Main-effect models of association with anthropometric phenotypes.

Figure 2 shows the sex-stratified analysis of WAA with BMI, WC, and WHR, as well as the respective interaction p-values for the interaction of sex with WAA. The interaction of sex with WAA is most pronounced for WHR (pinteraction = 4.14 × 10−5), whereby we observe a negative association of WAA with WHR among men. Among women, there is essentially no association between WAA and WHR. The interaction trend is similar but slightly weaker for WC (pinteraction = 1.44 × 10−4) and BMI (pinteraction = 5.35 × 10−4). Among women, there is a trend of a positive association of WAA with BMI and WC. Finally, we observed a trend in which the overall association of WAA with BMI and WC appears to be more positive in ARIC than in MESA, which is corroborated by results of the main-effects model shown in Table 2.

FIGURE 2
www.frontiersin.org

Figure 2. Sex-stratified analysis of association between WAA and BMI, WC, and WHR in ARIC, MESA and in combined dataset. Error bars represent 95% confidence intervals around the beta coefficients.

In the pooled dataset, we do not find any evidence for statistically significant interactions of WAA with either income or education among men or women, or of WAA with menopausal status among women (see Table 3).

TABLE 3
www.frontiersin.org

Table 3. Interaction of WAA with income, education, and menopausal status on three obesity-related phenotypes in men and/or women.

Discussion

We used genetic admixture/ancestry as a tool to investigate whether there may be a genetic basis to the large racial disparity in obesity that exists among women, but not among men. Our findings suggest that WAA is associated with protection against central obesity among AA men, but not among AA women. The pattern may explain to a great extent the much greater racial disparity among women (or alternatively, the much greater sex disparity among AA, as compared to EA). One interpretation is that AAs are exposed to environmental and/or cultural factors that predispose them to greater obesity than EAs. Possibly, some of the genes that are inherited as part of their West-African ancestry are protective against obesity, thereby “canceling out” the obesifying effects of environment/culture, but only in men. Another interpretation is that genetic protection is afforded to all individuals of African descent, but this protection is overwhelmed by cultural and/or other factors in women. For example, previous studies have found a tendency for a larger ideal body size among AA women as compared to White women (e.g., Rucker and Cash, 1992; Becker et al., 1999). In addition, evidence for greater obesity risk among women than men from socio-economic factors has been noted in a cross-country study (Wells et al., 2012), and in several studies examining gender and race gaps in obesity in the United States (Zhang and Wang, 2004; Robinson et al., 2009; Seamans et al., 2015).

Our findings are quite consistent across the two cohorts examined. However, we did observe a trend of a stronger positive association of WAA with BMI and WC in ARIC than in MESA in both genders. This difference could be attributed to differences in age, WAA, or other characteristics that differ between these two cohorts, as shown in Table 1. The pattern of stronger positive association among women however, holds up in both cohorts for all three phenotypes, and is quite consistent across both cohorts for WHR.

Previous studies have shown that AA have lower levels of central adiposity than EAs and other ethnic/racial groups (Yanovski et al., 1996; Beasley et al., 2009), and that WAA is associated with lower central adiposity (Cardel et al., 2011). However, some previous studies were limited to female participants, or did not examine sex differences in the association of WAA with obesity phenotypes. Here, we confirm the Goonesekera et al. finding that the association of WAA with BMI and WC is stronger among women. Our results suggest that this trend may be driven by the negative association of WAA with WC and WHR among men.

This study is strengthened by the use of two cohorts of AA with relatively large sample sizes, the use of multiple obesity-related measurements, and the availability of thousands of SNPs to more accurately estimate WAA. Limitations include the lack of more refined and metabolically relevant measurements such as amount of visceral adiposity, and the absence of other potentially important socio-economic confounders such as discrimination, and access to health-care. Regarding the latter limitation, it is difficult to tease genetic factors apart from environmental factors. Thus, it could be that variation in WAA is also capturing variation in environmental, social, cultural factors beyond those that we have included in our models as covariates.

This study suggests that there are specific genetic variants and physiological mechanism(s) that differ among West African and European populations, which when in the context of male sex, result in lower central adiposity. Further research, likely requiring larger sample sizes, is needed to identify these putative specific genetic factors and their corresponding physiological mechanisms. The identification and understanding of these mechanisms could lead to improved prevention and treatment approaches. It will also be important to also consider and examine non-genetic factors that may underlie the greater racial disparity among women.

Author Contributions

YK and DA contributed to the conception and design of this work. YK and AA analyzed the data. YK, JZ, RK, AA and DA contributed to the interpretation of the data, and to the writing and revising of the manuscript.

Funding

YK and AA are supported by a National Institutes of Health (NIH) grant: K01DK095032. DA is supported by NIH grant P30DK056336. The opinions expressed are those of the authors and do not necessarily represent those of the NIH or any other organization.

Conflict of Interest Statement

Dr. DA has received, anticipates, or has had financial interests with the Frontiers Foundation; University of Wisconsin; University of Arizona; Paul Weiss, Wharton and Garrison LLP; and Sage Publications.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank the participants and organizers of the ARIC and MESA studies. Data from these studies was obtained from dbGaP through accession numbers: phs000280.v2.p1 and phs000209.v10.p2. Atherosclerosis Risk in Communities: The Atherosclerosis Risk in Communities Study is carried out as a collaborative study supported by National Heart, Lung, and Blood Institute contracts (HHSN268201100005C, HHSN26 8201100006C, HHSN268201100007C, HHSN268201100008C, HHSN26820110 0009C, HHSN268201100010C, HHSN268201100011C, and HHSN26820110001 2C). Funding for GENEVA was provided by National Human Genome Research Institute grant U01HG004402 (E. Boerwinkle). The authors thank the staff and participants of the ARIC study for their important contributions. Multi-Ethnic Study of Atherosclerosis: MESA and the MESA SHARe project are conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with MESA investigators. Support for MESA is provided by contracts N01-HC-95159, N01-HC-95160, N01-HC-95161, N01-HC-95162, N01-HC-95163, N01-HC-95164, N01-HC-95165, N01-HC-95166, N01-HC-95167, N01-HC-95168, N01-HC-95169 and CTSA UL1-RR-024156. MESA Family is conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with MESA investigators. Support is provided by grants and contracts R01HL071051, R01HL071205, R01HL071250, R01HL071251, R01HL071258, R01HL071259, UL1-RR-025005, by the National Center for Research Resources, Grant UL1RR033176, and the National Center for Advancing Translational Sciences, Grant UL1TR000124. This manuscript was not prepared in collaboration with MESA investigators and does not necessarily reflect the opinions or views of MESA, or the NHLBI. Funding for SHARe genotyping was provided by NHLBI Contract N02-HL-64278. Genotyping was performed at Affymetrix (Santa Clara, California, USA) and the Broad Institute of Harvard and MIT (Boston, Massachusetts, USA) using the Affymetrix Genome-Wide Human SNP Array 6.0.

References

Beasley, L. E., Koster, A., Newman, A. B., Javaid, M. K., Ferrucci, L., Kritchevsky, S. B., et al. (2009). Inflammation and race and gender differences in computerized tomography-measured adipose depots. Obesity (Silver Spring) 17, 1062–1069. doi: 10.1038/oby.2008.627

PubMed Abstract | CrossRef Full Text | Google Scholar

Becker, D. M., Yanek, L. R., Koffman, D. M., and Bronner, Y. C. (1999). Body image preferences among urban African Americans and whites from low income communities. Ethn. Dis. 9, 377–386.

PubMed Abstract | Google Scholar

Bild, D. E., Bluemke, D. A., Burke, G. L., Detrano, R., Diez Roux, A. V., Folsom, A. R., et al. (2002). Multi-ethnic study of atherosclerosis: objectives and design. Am. J. Epidemiol. 156, 871–881. doi: 10.1093/aje/kwf113

PubMed Abstract | CrossRef Full Text | Google Scholar

Brancati, F. L., Kao, W. H., Folsom, A. R., Watson, R. L., and Szklo, M. (2000). Incident type 2 diabetes mellitus in african american and white adults: the atherosclerosis risk in communities study. JAMA 283, 2253–2259. doi: 10.1001/jama.283.17.2253

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardel, M., Higgins, P. B., Willig, A. L., Keita, A. D., Casazza, K., Gower, B. A., et al. (2011). African genetic admixture is associated with body composition and fat distribution in a cross-sectional study of children. Int. J. Obes. 35, 60–65. doi: 10.1038/ijo.2010.203

PubMed Abstract | CrossRef Full Text | Google Scholar

CDC/NCHS, I. (2013). National Health and Nutrition Examination Survey [Internet]. Avaliable online at: http://www.cdc.gov/nchs/data/hus/2013/069.pdf (Accessed May 10, 2016).

Fernandez, J. R., Shriver, M. D., Beasley, T. M., Rafla-Demetrious, N., Parra, E., Albu, J., et al. (2003). Association of African genetic admixture with resting metabolic rate and obesity among women. Obes. Res. 11, 904–911. doi: 10.1038/oby.2003.124

PubMed Abstract | CrossRef Full Text | Google Scholar

Goonesekera, S. D., Fang, S. C., Piccolo, R. S., Florez, J. C., and McKinlay, J. B. (2015). Biogeographic ancestry is associated with higher total body adiposity among African-American Females: the boston area community health survey. PLoS ONE 10:e0122808. doi: 10.1371/journal.pone.0122808

PubMed Abstract | CrossRef Full Text | Google Scholar

Gower, B. A., Nagy, T. R., and Goran, M. I. (1999). Visceral fat, insulin sensitivity, and lipids in prepubertal children. Diabetes 48, 1515–1521.

PubMed Abstract | Google Scholar

Halder, I., Kip, K. E., Mulukutla, S. R., Aiyer, A. N., Marroquin, O. C., Huggins, G. S., et al. (2012). Biogeographic ancestry, self-identified race, and admixture-phenotype associations in the heart SCORE Study. Am. J. Epidemiol. 176, 146–155. doi: 10.1093/aje/kwr518

PubMed Abstract | CrossRef Full Text | Google Scholar

Nassir, R., Qi, L., Kosoy, R., Garcia, L., Allison, M., Ochs-Balcom, H. M., et al. (2012). Relationship between adiposity and admixture in African-American and Hispanic-American women. Int. J. Obes. 36, 304–313. doi: 10.1038/ijo.2011.84

PubMed Abstract | CrossRef Full Text | Google Scholar

Ochs-Balcom, H. M., Preus, L., Wactawski-Wende, J., Nie, J., Johnson, N. A., Zakharia, F., et al. (2013). Association of DXA-derived bone mineral density and fat mass with african ancestry. J. Clin. Endocrinol. Metab. 98, E713–E717. doi: 10.1210/jc.2012-3921

PubMed Abstract | CrossRef Full Text | Google Scholar

Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics. 155, 945–959. Available online at: http://www.genetics.org/content/155/2/945

PubMed Abstract

R Development Core Team (2011). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Available online at: http://www.rproject.org

Robinson, W. R., Gordon-Larsen, P., Kaufman, J. S., Suchindran, C. M., and Stevens, J. (2009). The female-male disparity in obesity prevalence among black American young adults: contributions of sociodemographic characteristics of the childhood family. Am. J. Clin. Nutr. 89, 1204–1212. doi: 10.3945/ajcn.2007.25751

PubMed Abstract | CrossRef Full Text | Google Scholar

Rucker, C., and Cash, T. (1992). Body images, body-size perceptions, and eating behaviors among African-American and White College women. Int. J. Eat. Disord. 12, 291–299.

Google Scholar

Seamans, M. J., Robinson, W. R., Thorpe, R. J. Jr., Cole, S. R., and LaVeist, T. A. (2015). Exploring racial differences in the obesity gender gap. Ann. Epidemiol. 25, 420–425. doi: 10.1016/j.annepidem.2015.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Signorello, L. B., Schlundt, D. G., Cohen, S. S., Steinwandel, M. D., Buchowski, M. S., McLaughlin, J. K., et al. (2007). Comparing diabetes prevalence between African Americans and whites of similar socioeconomic status. Am. J. Public Health 97, 2260–2267. doi: 10.2105/AJPH.2006.094482

PubMed Abstract | CrossRef Full Text | Google Scholar

Tandon, A., Patterson, N., and Reich, D. (2011). Ancestry informative marker panels for African Americans based on subsets of commercially available SNP arrays. Genet. Epidemiol. 35, 80–83. doi: 10.1002/gepi.20550

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, H., Jorgenson, E., Gadde, M., Kardia, S. L., Rao, D. C., Zhu, X., et al. (2006). Racial admixture and its impact on BMI and blood pressure in African and Mexican Americans. Hum. Genet. 119, 624–633. doi: 10.1007/s00439-006-0175-4

PubMed Abstract | CrossRef Full Text | Google Scholar

The Aric Investigators (1989). The Atherosclerosis Risk in communities (ARIC) study: design and objectives. the ARIC investigators. Am. J. Epidemiol. 129, 687–702.

PubMed Abstract

The International HapMap Consortium (2003). The International HapMap Project. Nature 426, 789–796. doi: 10.1038/nature02168

PubMed Abstract | CrossRef Full Text

Wells, J. C. K., Marphatia, A. A., Cole, T. J., and McCoy, D. (2012). Associations of economic and gender inequality with global obesity prevalence: understanding the female excess. Soc. Sci. Med. 75, 482–490. doi: 10.1016/j.socscimed.2012.03.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Yanovski, J. A., Yanovski, S. Z., Filmer, K. M., Hubbard, V. S., Avila, N., Lewis, B., et al. (1996). Differences in body composition of black and white girls. Am J. Clin. Nutr. 64, 833–839.

PubMed Abstract | Google Scholar

Zhang, Q., and Wang, Y. (2004). Socioeconomic inequality of obesity in the United States: do gender, age, and ethnicity matter? Soc. Sci. Med. 58, 1171–1180. doi: 10.1016/S0277-9536(03)00288-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: African-American, genetic admixture, genetic ancestry, obesity, central obesity, waist-to-hip ratio, health disparities, sex differences

Citation: Klimentidis YC, Arora A, Zhou J, Kittles R and Allison DB (2016) The Genetic Contribution of West-African Ancestry to Protection against Central Obesity in African-American Men but Not Women: Results from the ARIC and MESA Studies. Front. Genet. 7:89. doi: 10.3389/fgene.2016.00089

Received: 07 February 2016; Accepted: 01 May 2016;
Published: 01 June 2016.

Edited by:

Geoff Hayes, Northwestern University, USA

Reviewed by:

Tesfaye B. Mersha, Cincinnati Children's Hospital Medical Center, USA
Amy Rebecca Bentley, National Human Genome Research Institute, USA
Anne E. Sumner, National Institutes of Health, USA

Copyright © 2016 Klimentidis, Arora, Zhou, Kittles and Allison. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yann C. Klimentidis, yann@email.arizona.edu