Ancient Rapanui genomes reveal resilience and pre-European contact with the Americas.

Rapa Nui (also known as Easter Island) is one of the most isolated inhabited places in the world. It has captured the imagination of many owing to its archaeological record, which includes iconic megalithic statues called moai1. Two prominent contentions have arisen from the extensive study of Rapa Nui. First, the history of the Rapanui has been presented as a warning tale of resource overexploitation that would have culminated in a major population collapse-the 'ecocide' theory2-4. Second, the possibility of trans-Pacific voyages to the Americas pre-dating European contact is still debated5-7. Here, to address these questions, we reconstructed the genomic history of the Rapanui on the basis of 15 ancient Rapanui individuals that we radiocarbon dated (1670-1950 CE) and whole-genome sequenced (0.4-25.6×). We find that these individuals are Polynesian in origin and most closely related to present-day Rapanui, a finding that will contribute to repatriation efforts. Through effective population size reconstructions and extensive population genetics simulations, we reject a scenario involving a severe population bottleneck during the 1600s, as proposed by the ecocide theory. Furthermore, the ancient and present-day Rapanui carry similar proportions of Native American admixture (about 10%). Using a Bayesian approach integrating genetic and radiocarbon dates, we estimate that this admixture event occurred about 1250-1430 CE.

Identification of Sociodemographic, Clinical, and Genetic Factors to Aid Alaska Native and American Indian People to Successfully Quit Smoking.

INTRODUCTION: Alaska Native and American Indian (ANAI) people have a smoking prevalence of 23%. Nicotine metabolite ratio (NMR) and genetic testing may enable tailored selection of tobacco cessation medication. AIMS AND METHODS: The purpose of this study was to evaluate the relative contributions of NMR, cessation medication, demographics, and tobacco use history to cessation. Participants were recruited into an observational cohort study consisting of a baseline visit prior to their quit date and 6-week follow-up. Demographic and tobacco use surveys and blood, urine, and breath samples were collected at each visit. Electronic health records were queried for cessation medications. NMR was categorized into slow or normal nicotine metabolism phenotypes (<0.31 and ≥ 0.31, respectively). The main outcome was cessation at 6 weeks. Analyses consisted of descriptive statistics, medication and phenotype concordance, and estimates of relative risk (RR) of quitting. RESULTS: We enrolled 151 ANAI adults who smoked cigarettes daily. Two-thirds had normal nicotine metabolism phenotype. Retrospective medication and phenotype concordance was 39%. The overall quit rate was 25%. No demographic factors or tobacco use history were associated with quit success. Varenicline and bupropion increased the likelihood of quitting (RR = 2.93 [1.42, 6.03] and RR = 2.52 [1.12, 5.64], respectively) compared to nicotine replacement therapy. Non-optimal medication and phenotype concordance decreased likelihood of quit success (RR = 0.44 [0.22, 0.91]) compared to optimal concordance. CONCLUSIONS: This exploratory study found associations between quit success and tobacco cessation medication as well as medication and phenotype concordance. Additional research is needed to assess use of NMR for treatment selection among ANAI people. IMPLICATIONS: These results broadly support additional community-engaged research to improve medication and phenotype concordance in tribal health settings. Such future research on implementing meditcation and phenotype concordance holds promise to improve expectations, quit success, and health outcomes amongst individuals attempting to quit smoking.

Native American ancestry and breast cancer risk in Colombian and Mexican women: ruling out potential confounding through ancestry-informative markers.

BACKGROUND: Latin American and Hispanic women are less likely to develop breast cancer (BC) than women of European descent. Observational studies have found an inverse relationship between the individual proportion of Native American ancestry and BC risk. Here, we use ancestry-informative markers to rule out potential confounding of this relationship, estimating the confounder-free effect of Native American ancestry on BC risk. METHODS AND STUDY POPULATION: We used the informativeness for assignment measure to select robust instrumental variables for the individual proportion of Native American ancestry. We then conducted separate Mendelian randomization (MR) analyses based on 1401 Colombian women, most of them from the central Andean regions of Cundinamarca and Huila, and 1366 Mexican women from Mexico City, Monterrey and Veracruz, supplemented by sensitivity and stratified analyses. RESULTS: The proportion of Colombian Native American ancestry showed a putatively causal protective effect on BC risk (inverse variance-weighted odds ratio [OR] = 0.974 per 1% increase in ancestry proportion, 95% confidence interval [CI] 0.970-0.978, p = 3.1 × 10-40). The corresponding OR for Mexican Native American ancestry was 0.988 (95% CI 0.987-0.990, p = 1.4 × 10-44). Stratified analyses revealed a stronger association between Native American ancestry and familial BC (Colombian women: OR = 0.958, 95% CI 0.952-0.964; Mexican women: OR = 0.973, 95% CI 0.969-0.978), and stronger protective effects on oestrogen receptor (ER)-positive BC than on ER-negative and triple-negative BC. CONCLUSIONS: The present results point to an unconfounded protective effect of Native American ancestry on BC risk in both Colombian and Mexican women which appears to be stronger for familial and ER-positive BC. These findings provide a rationale for personalised prevention programmes that take genetic ancestry into account, as well as for future admixture mapping studies.

Genetic Architecture of Dilated Cardiomyopathy in Individuals of African and European Ancestry.

Importance: Black patients with dilated cardiomyopathy (DCM) have increased familial risk and worse outcomes than White patients, but most DCM genetic data are from White patients. Objective: To compare the rare variant genetic architecture of DCM by genomic ancestry within a diverse population of patients with DCM. Design: Cross-sectional study enrolling patients with DCM who self-identified as non-Hispanic Black, Hispanic, or non-Hispanic White from June 7, 2016, to March 15, 2020, at 25 US advanced heart failure programs. Variants in 36 DCM genes were adjudicated as pathogenic, likely pathogenic, or of uncertain significance. Exposure: Presence of DCM. Main Outcomes and Measures: Variants in DCM genes classified as pathogenic/likely pathogenic/uncertain significance and clinically actionable (pathogenic/likely pathogenic). Results: A total of 505, 667, and 26 patients with DCM of predominantly African, European, or Native American genomic ancestry, respectively, were included. Compared with patients of European ancestry, a lower percentage of patients of African ancestry had clinically actionable variants (8.2% [95% CI, 5.2%-11.1%] vs 25.5% [95% CI, 21.3%-29.6%]), reflecting the lower odds of a clinically actionable variant for those with any pathogenic variant/likely pathogenic variant/variant of uncertain significance (odds ratio, 0.25 [95% CI, 0.17-0.37]). On average, patients of African ancestry had fewer clinically actionable variants in TTN (difference, -0.09 [95% CI, -0.14 to -0.05]) and other genes with predicted loss of function as a disease-causing mechanism (difference, -0.06 [95% CI, -0.11 to -0.02]). However, the number of pathogenic variants/likely pathogenic variants/variants of uncertain significance was more comparable between ancestry groups (difference, -0.07 [95% CI, -0.22 to 0.09]) due to a larger number of non-TTN non-predicted loss of function variants of uncertain significance, mostly missense, in patients of African ancestry (difference, 0.15 [95% CI, 0.00-0.30]). Published clinical case-based evidence supporting pathogenicity was less available for variants found only in patients of African ancestry (P < .001). Conclusion and Relevance: Patients of African ancestry with DCM were less likely to have clinically actionable variants in DCM genes than those of European ancestry due to differences in genetic architecture and a lack of representation of African ancestry in clinical data sets.

Ancestry-driven metabolite variation provides insights into disease states in admixed populations.

BACKGROUND: Metabolic pathways are related to physiological functions and disease states and are influenced by genetic variation and environmental factors. Hispanics/Latino individuals have ancestry-derived genomic regions (local ancestry) from their recent admixture that have been less characterized for associations with metabolite abundance and disease risk. METHODS: We performed admixture mapping of 640 circulating metabolites in 3887 Hispanic/Latino individuals from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). Metabolites were quantified in fasting serum through non-targeted mass spectrometry (MS) analysis using ultra-performance liquid chromatography-MS/MS. Replication was performed in 1856 nonoverlapping HCHS/SOL participants with metabolomic data. RESULTS: By leveraging local ancestry, this study identified significant ancestry-enriched associations for 78 circulating metabolites at 484 independent regions, including 116 novel metabolite-genomic region associations that replicated in an independent sample. Among the main findings, we identified Native American enriched genomic regions at chromosomes 11 and 15, mapping to FADS1/FADS2 and LIPC, respectively, associated with reduced long-chain polyunsaturated fatty acid metabolites implicated in metabolic and inflammatory pathways. An African-derived genomic region at chromosome 2 was associated with N-acetylated amino acid metabolites. This region, mapped to ALMS1, is associated with chronic kidney disease, a disease that disproportionately burdens individuals of African descent. CONCLUSIONS: Our findings provide important insights into differences in metabolite quantities related to ancestry in admixed populations including metabolites related to regulation of lipid polyunsaturated fatty acids and N-acetylated amino acids, which may have implications for common diseases in populations.

Native American genetic ancestry and pigmentation allele contributions to skin color in a Caribbean population.

Our interest in the genetic basis of skin color variation between populations led us to seek a Native American population with genetically African admixture but low frequency of European light skin alleles. Analysis of 458 genomes from individuals residing in the Kalinago Territory of the Commonwealth of Dominica showed approximately 55% Native American, 32% African, and 12% European genetic ancestry, the highest Native American genetic ancestry among Caribbean populations to date. Skin pigmentation ranged from 20 to 80 melanin units, averaging 46. Three albino individuals were determined to be homozygous for a causative multi-nucleotide polymorphism OCA2NW273KV contained within a haplotype of African origin; its allele frequency was 0.03 and single allele effect size was -8 melanin units. Derived allele frequencies of SLC24A5A111T and SLC45A2L374F were 0.14 and 0.06, with single allele effect sizes of -6 and -4, respectively. Native American genetic ancestry by itself reduced pigmentation by more than 20 melanin units (range 24-29). The responsible hypopigmenting genetic variants remain to be identified, since none of the published polymorphisms predicted in prior literature to affect skin color in Native Americans caused detectable hypopigmentation in the Kalinago.

The variation in skin colour of modern humans is a product of thousands of years of natural selection. All human ancestry can be traced back to African populations, which were dark-skinned to protect them from the intense UV rays of the sun. Over time, humans spread to other parts of the world, and people in the northern latitudes with lower UV developed lighter skin through natural selection. This was likely driven by a need for vitamin D, which requires UV rays for production. Separate genetic mechanisms were involved in the evolution of lighter skin in each of the two main branches of human migration: the European branch (which includes peoples on the Indian subcontinent and Europe) and the East Asian branch (which includes East Asia and the Americas). A variant of the gene SLC24A5 is the primary contributor to lighter skin colour in the European branch, but a corresponding variant driving light skin colour evolution in the East Asian branch remains to be identified. One obstacle to finding such variants is the high prevalence of European ancestry in most people groups, which makes it difficult to separate the influence of European genes from those of other populations. To overcome this issue, Ang et al. studied a population that had a high proportion of Native American and African ancestors, but a relatively small proportion of European ancestors, the Kalinago people. The Kalinago live on the island of Dominica, one of the last Caribbean islands to be colonised by Europeans. Ang et al. were able to collect hundreds of skin pigmentation measurements and DNA samples of the Kalinago, to trace the effect of Native American ancestry on skin colour. Genetic analysis confirmed their oral history records of primarily Native American (55%) ­ one of the highest of any Caribbean population studied to date ­ compared with African (32%) and European (12%) ancestries. Native American ancestry had the highest effect on pigmentation and reduced it by more than 20 melanin units, while the European mutations in the genes SLC24A5 and SLC45A2 and an African gene variant for albinism only contributed 5, 4 and 8 melanin units, respectively. However, none of the so far published gene candidates responsible for skin lightening in Native Americans caused a detectable effect. Therefore, the gene responsible for lighter skin in Native Americans/East Asians has yet to be identified. The work of Ang et al. represents an important step in deciphering the genetic basis of lighter skin colour in Native Americans or East Asians. A better understanding of the genetics of skin pigmentation may help to identify why, for example, East Asians are less susceptible to melanoma than Europeans, despite both having a lighter skin colour. It may also further acceptance of how variations in human skin tones are the result of human migration, random genetic variation, and natural selection for pigmentation in different solar environments.

Bayesian inference of admixture graphs on Native American and Arctic populations.

Admixture graphs are mathematical structures that describe the ancestry of populations in terms of divergence and merging (admixing) of ancestral populations as a graph. An admixture graph consists of a graph topology, branch lengths, and admixture proportions. The branch lengths and admixture proportions can be estimated using numerous numerical optimization methods, but inferring the topology involves a combinatorial search for which no polynomial algorithm is known. In this paper, we present a reversible jump MCMC algorithm for sampling high-probability admixture graphs and show that this approach works well both as a heuristic search for a single best-fitting graph and for summarizing shared features extracted from posterior samples of graphs. We apply the method to 11 Native American and Siberian populations and exploit the shared structure of high-probability graphs to characterize the relationship between Saqqaq, Inuit, Koryaks, and Athabascans. Our analyses show that the Saqqaq is not a good proxy for the previously identified gene flow from Arctic people into the Na-Dene speaking Athabascans.

An ancient founder mutation located between ROBO1 and ROBO2 is responsible for increased microtia risk in Amerindigenous populations.

Microtia is a congenital malformation that encompasses mild hypoplasia to complete loss of the external ear, or pinna. Although the contribution of genetic variation and environmental factors to microtia remains elusive, Amerindigenous populations have the highest reported incidence. Here, using both transmission disequilibrium tests and association studies in microtia trios (parents and affected child) and microtia cohorts enrolled in Latin America, we map an ∼10-kb microtia locus (odds ratio = 4.7; P = 6.78e-18) to the intergenic region between Roundabout 1 (ROBO1) and Roundabout 2 (ROBO2) (chr3: 78546526 to 78555137). While alleles at the microtia locus significantly increase the risk of microtia, their penetrance is low (<1%). We demonstrate that the microtia locus contains a polymorphic complex repeat element that is expanded in affected individuals. The locus is located near a chromatin loop region that regulates ROBO1 and ROBO2 expression in induced pluripotent stem cell­derived neural crest cells. Furthermore, we use single nuclear RNA sequencing to demonstrate ROBO1 and ROBO2 expression in both fibroblasts and chondrocytes of the mature human pinna. Because the microtia allele is enriched in Amerindigenous populations and is shared by some East Asian subjects with craniofacial malformations, we propose that both populations share a mutation that arose in a common ancestor prior to the ancient migration of Eurasian populations into the Americas and that the high incidence of microtia among Amerindigenous populations reflects the population bottleneck that occurred during the migration out of Eurasia.

Genetic Distribution of Five Spinocerebellar Ataxia Microsatellite Loci in Mexican Native American Populations and Its Impact on Contemporary Mestizo Populations.

Spinocerebellar ataxias (SCAs) conform a heterogeneous group of neurodegenerative disorders with autosomal dominant inheritance. Five of the most frequent SCAs are caused by a CAG repeat expansion in the exons of specific genes. The SCAs incidence and the distribution of polymorphic CAG alleles vary among populations and ethnicities. Thus, characterization of the genetic architecture of ethnically diverse populations, which have undergone recent admixture and demographic events, could facilitate the identification of genetic risk factors. Owing to the great ethnic diversity of the Mexican population, this study aimed to analyze the allele frequencies of five SCA microsatellite loci (SCA1, SCA2, SCA3, SCA6, and SCA7) in eleven Mexican Native American (MNA) populations. Data from the literature were used to compare the allelic distribution of SCA loci with worldwide populations. The SCA loci allelic frequencies evidenced a certain genetic homogeneity in the MNA populations, except for Mayans, who exhibited distinctive genetic profiles. Neither pathological nor large normal alleles were found in MNA populations, except for the SCA2 pre-mutated allele in the Zapotec population. Collectively, our findings demonstrated the contribution of the MNA ancestry in shaping the genetic structure of contemporary Mexican Mestizo populations. Our results also suggest that Native American ancestry has no impact on the origin of SCAs in the Mexican population. Instead, the acquisition of pathological SCA alleles could be associated with European migration.

The Mitochondrial DNA Landscape of Modern Mexico.

Mexico is a rich source for anthropological and population genetic studies with high diversity in ethnic and linguistic groups. The country witnessed the rise and fall of major civilizations, including the Maya and Aztec, but resulting from European colonization, the population landscape has dramatically changed. Today, the majority of Mexicans do not identify themselves as Indigenous but as admixed, and appear to have very little in common with their pre-Columbian predecessors. However, when the maternally inherited mitochondrial (mt)DNA is investigated in the modern Mexican population, this is not the case. Control region sequences of 2021 samples deriving from all over the country revealed an overwhelming Indigenous American legacy, with almost 90% of mtDNAs belonging to the four major pan-American haplogroups A2, B2, C1, and D1. This finding supports a very low European contribution to the Mexican gene pool by female colonizers and confirms the effectiveness of employing uniparental markers as a tool to reconstruct a country's history. In addition, the distinct frequency and dispersal patterns of Indigenous American and West Eurasian clades highlight the benefit such large and country-wide databases provide for studying the impact of colonialism from a female perspective and population stratification. The importance of geographical database subsets not only for forensic application is clearly demonstrated.

Population diversity of three variants of the SLC47A2 gene (MATE2-K transporter) in Mexican Mestizos and Native Americans.

BACKGROUND: MATE2-K is an efflux transporter protein of organic cation expressed mainly in the kidney and encoded by the SLC47A2 gene. Different variants of this gene have shown an impact on the pharmacokinetics of various drugs, including metformin, which represents one of the most widely used drugs in treating type 2 diabetes. The SLC47A2 gene variants have been scarcely studied in Mexican populations, especially in Native American groups. For this reason, we analyzed the distribution of the variants rs12943590, rs35263947, and rs9900497 within the SLC47A2 gene in 173 Native Americans (Tarahumara, Huichol, Maya, Puerépecha) and 182 Mestizos (admixed) individuals from Mexico. METHODS AND RESULTS: Genotypes were determined through TaqMan probes (qPCR). The Hardy-Weinberg agreement was confirmed for all three SLC47A2 gene variants in all the Mexican populations analyzed. When worldwide populations were included for comparison purposes, for alleles and genotypes a relative interpopulation homogeneity was observed for rs35263947 (T allele; range 23.3-51.1%) and rs9900497 (T allele; range 18.6-40.9%). Conversely, heterogeneity was evident for rs12943590 (A allele, range 22.1-59.1%), where the most differentiated population was the Huichol, with high frequencies of the risk genotype associated with decreased response to metformin treatment (A/A = 40.9%). CONCLUSIONS: Although the SLC47A2 gene variants allow predicting favorable response to the metformin treatment in Mexican populations, the probable high frequency of ineffectiveness should be discarded in Huichols.

How natural selection shapes genetic differentiation in the MHC region: A case study with Native Americans.

The Human Leukocyte Antigen (HLA) loci are extremely well documented targets of balancing selection, yet few studies have explored how selection affects population differentiation at these loci. In the present study we investigate genetic differentiation at HLA genes by comparing differentiation at microsatellites distributed genomewide to those in the MHC region. Our study uses a sample of 494 individuals from 30 human populations, 28 of which are Native Americans, all of whom were typed for genomewide and MHC region microsatellites. We find greater differentiation in the MHC than in the remainder of the genome (FST-MHC = 0.130 and FST-Genomic = 0.087), and use a permutation approach to show that this difference is statistically significant, and not accounted for by confounding factors. This finding lies in the opposite direction to the expectation that balancing selection reduces population differentiation. We interpret our findings as evidence that selection favors different sets of alleles in distinct localities, leading to increased differentiation. Thus, balancing selection at HLA genes simultaneously increases intra-population polymorphism and inter-population differentiation in Native Americans.

Population structure and relatedness estimates in a Mexican sample.

Population stratification (PS) is a confounding factor in genome-wide association studies (GWASs) and also an interesting process itself. Latin American populations have mixed genetic ancestry, which may account for PS. We have analyzed the relatedness, by means of the identity-by-descent (IBD) estimations, in a sample of 1805 individuals and 1.006.703 autosomal mutations from a case-control study of colorectal cancer in Mexico. When using the recommended protocol for quality control assessment, 402 should have been removed due to relatedness. Our purpose was to analyze this value in the context of an admixed population. For that aim, we reanalyzed the sample using two software designed for admixed populations, obtaining estimates of 110 and 70 related individuals to remove. The results showed that the first estimation of relatedness was an effect of the higher Native American contribution in part of the data samples, being a confounding factor for IBD estimations. We conclude in the importance of considering PS and genetic ancestry in order to avoid spurious results, not only in GWAS but also in relatedness analysis.

Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

There has been limited study of Native American whole genome diversity to date, which impairs effective implementation of personalized medicine and a detailed description of its demographic history. Here we report high coverage whole genome sequencing of 76 unrelated individuals, from 27 indigenous groups across Mexico, with more than 97% average Native American ancestry. On average, each individual has 3.26 million Single Nucleotide Variants and short indels, that together comprise a catalog of 9,737,152 variants, 44,118 of which are novel. We report 497 common Single Nucleotide Variants (with allele frequency > 5%) mapped to drug responses and 316,577 in enhancer or promoter elements; interestingly we found some of these enhancer variants in PPARG, a nuclear receptor involved in highly prevalent health problems in Mexican population, such as obesity, diabetes, and insulin resistance. By detecting signals of positive selection we report 24 enriched key pathways under selection, most of them related to immune mechanisms. No missense variants in ACE2, the receptor responsible for the entry of the SARS CoV-2 virus, were found in any individual. Population genomics and phylogenetic analyses demonstrated stratification in a Northern-Central-Southern axis, with major substructure in the Central region. The Seri, a northern group with the most genetic divergence in our study, showed a distinctive genomic context with the most novel variants, and the most population specific genotypes. Genome-wide analysis showed that the average haplotype blocks are longer in Native Mexicans than in other world populations. With this dataset we describe previously undetected population level variation in Native Mexicans, helping to reduce the gap in genomic data representation of such groups.

Archaeogenomic distinctiveness of the Isthmo-Colombian area.

The recently enriched genomic history of Indigenous groups in the Americas is still meager concerning continental Central America. Here, we report ten pre-Hispanic (plus two early colonial) genomes and 84 genome-wide profiles from seven groups presently living in Panama. Our analyses reveal that pre-Hispanic demographic events contributed to the extensive genetic structure currently seen in the area, which is also characterized by a distinctive Isthmo-Colombian Indigenous component. This component drives these populations on a specific variability axis and derives from the local admixture of different ancestries of northern North American origin(s). Two of these ancestries were differentially associated to Pleistocene Indigenous groups that also moved into South America, leaving heterogenous genetic footprints. An additional Pleistocene ancestry was brought by a still unsampled population of the Isthmus (UPopI) that remained restricted to the Isthmian area, expanded locally during the early Holocene, and left genomic traces up to the present day.

Dissecting the role of Amerindian genetic ancestry and the ApoE ε4 allele on Alzheimer disease in an admixed Peruvian population.

Alzheimer disease (AD) is the leading cause of dementia in the elderly and occurs in all ethnic and racial groups. The apolipoprotein E (ApoE) ε4 is the most significant genetic risk factor for late-onset AD and shows the strongest effect among East Asian populations followed by non-Hispanic white populations and has a relatively lower effect in African descent populations. Admixture analysis in the African American and Puerto Rican populations showed that the variation in ε4 risk is correlated with the genetic ancestral background local to the ApoE gene. Native American populations are substantially underrepresented in AD genetic studies. The Peruvian population with up to ~80 of Amerindian (AI) ancestry provides a unique opportunity to assess the role of AI ancestry in AD. In this study, we assess the effect of the ApoE ε4 allele on AD in the Peruvian population. A total of 79 AD cases and 128 unrelated cognitive healthy controls from Peruvian population were included in the study. Genome-wide genotyping was performed using the Illumina Global screening array v2.0. Global ancestry and local ancestry analyses were assessed. The effect of the ApoE ε4 allele on AD was tested using a logistic regression model by adjusting for age, gender, and population substructure (first 3 principal components). Results showed that the genetic ancestry surrounding the ApoE gene is predominantly AI (60.6%) and the ε4 allele is significantly associated with increased risk of AD in the Peruvian population (odds ratio = 5.02, confidence interval: 2.3-12.5, p-value = 2e-4). Our results showed that the risk for AD from ApoE ε4 in Peruvians is higher than we have observed in non-Hispanic white populations. Given the high admixture of AI ancestry in the Peruvian population, it suggests that the AI genetic ancestry local to the ApoE gene is contributing to a strong risk for AD in ε4 carriers. Our data also support the findings of an interaction between the genetic risk allele ApoE ε4 and the ancestral backgrounds located around the genomic region of ApoE gene.

Novel DNA methylation signatures of tobacco smoking with trans-ethnic effects.

BACKGROUND: Smoking remains one of the leading preventable causes of death. Smoking leaves a strong signature on the blood methylome as shown in multiple studies using the Infinium HumanMethylation450 BeadChip. Here, we explore novel blood methylation smoking signals on the Illumina MethylationEPIC BeadChip (EPIC) array, which also targets novel CpG-sites in enhancers. METHOD: A smoking-methylation meta-analysis was carried out using EPIC DNA methylation profiles in 1407 blood samples from four UK population-based cohorts, including the MRC National Survey for Health and Development (NSHD) or 1946 British birth cohort, the National Child Development Study (NCDS) or 1958 birth cohort, the 1970 British Cohort Study (BCS70), and the TwinsUK cohort (TwinsUK). The overall discovery sample included 269 current, 497 former, and 643 never smokers. Replication was pursued in 3425 trans-ethnic samples, including 2325 American Indian individuals participating in the Strong Heart Study (SHS) in 1989-1991 and 1100 African-American participants in the Genetic Epidemiology Network of Arteriopathy Study (GENOA). RESULTS: Altogether 952 CpG-sites in 500 genes were differentially methylated between smokers and never smokers after Bonferroni correction. There were 526 novel smoking-associated CpG-sites only profiled by the EPIC array, of which 486 (92%) replicated in a meta-analysis of the American Indian and African-American samples. Novel CpG sites mapped both to genes containing previously identified smoking-methylation signals and to 80 novel genes not previously linked to smoking, with the strongest novel signal in SLAMF7. Comparison of former versus never smokers identified that 37 of these sites were persistently differentially methylated after cessation, where 16 represented novel signals only profiled by the EPIC array. We observed a depletion of smoking-associated signals in CpG islands and an enrichment in enhancer regions, consistent with previous results. CONCLUSION: This study identified novel smoking-associated signals as possible biomarkers of exposure to smoking and may help improve our understanding of smoking-related disease risk.

Mitochondrial DNA analysis of the putative skeletal remains of Sieur de Marle: Genetic support for anthropological assessment of biogeographic ancestry.

In 1932, seven burials were discovered on a Texas plantation that was originally the site of a 17th-century Caddo Indian village. Of the seven excavated graves, one set of remains (an adult male) was notably buried in a manner inconsistent with traditional Caddoan burial practices and has long been purported to be the remains of Sieur de Marle (a member of the French explorer La Salle's last expedition). Diary accounts of La Salle's expedition scribe report that Sieur de Marle died along a river near an Indian village during a trek to Canada to find help for colonists left behind at the ill-fated Fort St. Louis. Additionally, two lead projectiles recovered from the grave were ballistically analyzed and determined to be consistent with ammunition used in 17th-century weaponry. In the 1980s, anthropologists requested access to the remains for study, but the skull was missing. Cranial measurements recorded in 1940 and 1962 (by two independent anthropologists) were used to investigate the ancestry of this individual; and the Giles-Elliot (G-E) discriminant function was calculated to be 18.1, within the Anglo-European range. Dietary isotope testing on non-cranial skeletal elements determined that this unknown male's diet was rich in animal/marine protein sources, which differs appreciably from Caddo Indian populations of that time period. In order to genetically assess this individual's biogeographic ancestry and to provide further support that this individual is of European descent, mitochondrial DNA (mtDNA) sequencing was performed using the Applied Biosystems™ Precision ID mtDNA Whole Genome Panel. mtDNA sequencing of multiple sections from two different long bones yielded compiled results consistent with either Haplogroup H or R, both predominantly European mtDNA haplogroups. Further anthropological calculations were conducted using cranial measurements, FORDISC™ software, and discriminant function analysis. Two-way, four-way, and multigroup discriminant function analyses further classify this set of unidentified remains as being White (European) in origin, with posterior probabilities of 0.999, 0.881 and 0.986, respectively. Combined with historical records of Sieur de Marle's death, as well as overlays of historical and contemporary maps which demonstrate that the plantation site aligns with Joutel's diary accounts of de Marle's burial, these collective results support that these remains are of a European male and may possibly belong to this prominent member of La Salle's expedition team.

Maternal cardiometabolic factors and genetic ancestry influence epigenetic aging of the placenta.

Disruption of physiological aging of the placenta can lead to pregnancy complications and increased risk for cardiometabolic diseases during childhood and adulthood. Maternal metabolic and genetic factors need to operate in concert with placental development for optimal pregnancy outcome. However, it is unknown whether maternal cardiometabolic status and genetic ancestry contribute to differences in placental epigenetic age acceleration (PAA). We investigated whether maternal prepregnancy obesity, gestational weight gain (GWG), blood pressure, and genetic ancestry influence PAA. Among 301 pregnant women from 4 race/ethnic groups who provided placenta samples at delivery as part of the National Institute of Child Health and Human Development Fetal Growth Studies, placental DNA methylation age was estimated using 62 CpGs known to predict placental aging. PAA was defined to be the difference between placental DNA methylation age and gestational age at birth. Percentage of genetic ancestries was estimated using genotype data. We found that a 1 kg/week increase in GWG was associated with up to 1.71 (95% CI: -3.11, -0.32) week lower PAA. Offspring Native American ancestry and African ancestry were associated, respectively, with higher and lower PAA among Hispanics, and maternal East Asian ancestry was associated with lower PAA among Asians (p < 0.05). Among mothers with a male offspring, blood pressure was associated with lower PAA across all three trimesters (p < 0.05), prepregnancy obesity compared to normal weight was associated with 1.24 (95% CI: -2.24, -0.25) week lower PAA. In summary, we observed that maternal cardiometabolic factors and genetic ancestry influence placental epigenetic aging and some of these influences may be male offspring-specific.