Major Ancestry Research Articles

A Catalogue of Structural Variation across Ancestrally Diverse Asian Genomes

Structural variants (SVs) are significant contributors to inter-individual genetic variation associated with traits and diseases. Current SV studies using whole-genome sequencing (WGS) have a largely Eurocentric composition, with little known about SV diversity in other ancestries, particularly from Asia. Here, we present a WGS catalogue of 73,035 SVs from 8392 Singaporeans of East Asian, Southeast Asian and South Asian ancestries, of which ~65% (47,770 SVs) are novel. We show that Asian populations can be stratified by their global SV patterns and identified 42,239 novel SVs that are specific to Asian populations. 52% of these novel SVs are restricted to one of the three major ancestry groups studied (Indian, Chinese or Malay). We uncovered SVs affecting major clinically actionable loci. Lastly, by identifying SVs in linkage disequilibrium with single-nucleotide variants, we demonstrate the utility of our SV catalogue in the fine-mapping of Asian GWAS variants and identification of potential causative variants. These results augment our knowledge of structural variation across human populations, thereby reducing current ancestry biases in global references of genetic variation afflicting equity, diversity and inclusion in genetic research.

Genetic ancestry proportion influences risk of adverse events from tuberculosis treatment in Brazil.

Tuberculosis (TB) treatment is highly effective, but response to therapy can vary by geography, race, and ethnicity. We assessed for differences in TB treatment response in a representative and heterogeneous Brazilian population. We estimated genetic ancestry proportion according to major ancestry groups (African, European, and Amerindian ancestry) in the Regional Prospective Observational Research in Tuberculosis (RePORT)-Brazil cohort. RePORT-Brazil is an observational prospective cohort study of individuals with newly-diagnosed, culture-confirmed, pulmonary TB. TB treatment outcomes that were attributed to TB treatment included Grade 2 or higher adverse drug reaction (ADR), Grade 3 or higher ADR, hepatic ADR, and failure/recurrence. Ancestry proportion was the main predictor in logistic regression for each outcome, with adjustments for candidate confounders. There were 941 pulmonary TB patients included in this study. We observed a decreased risk of Grade 2+ ADR when African ancestry proportion increased by 10% (Odds Ratio [OR] 0.41, 95% Confidence Interval [CI] 0.20 -0.85) and an increased risk for Grade 2+ ADR with increasing European ancestry (OR 2.84, 95% CI 1.47 - 5.48). We then performed the same analysis adding HIV status as an interaction term. The decreased risk for Grade 2+ ADR seen for African ancestry proportion did not hold for persons living with HIV; we observed increased risk for Grade 2+ ADR with increasing African ancestry proportion. There were no associations with Amerindian ancestry or for other treatment outcomes. In this Brazilian TB cohort, toxicity risk was associated with African and European ancestry, divergent, and affected by HIV. #RePORT-Brazil Consortium members include: Aline Benjamin and Flavia M. Sant'Anna Instituto Nacional de Infectologia Evandro Chagas, Fiocruz, Rio de Janeiro, Brazil Jamile Garcia de Oliveira and João Marine Clínica de Saúde Rinaldo Delmare, Rio de Janeiro, Brazil Adriana Rezende and Anna Cristina Carvalho Secretaria de Saúde de Duque de Caxias, Rio de Janeiro, Brazil Michael Rocha and Betânia Nogueira Instituto Brasileiro para Investigação da Tuberculose, Fundação José Silveira, Salvador, Brazil Alexandra Brito and Renata Spener Fundação Medicina Tropical Dr. Heitor Vieira Dourado, Manaus, Brazil Megan Turner Vanderbilt University Medical Center, Nashville, USA.

IFNL4 genotype and other personal characteristics to predict response to 8-week sofosbuvir-based treatment for chronic hepatitis C

Shorter duration therapy for hepatitis C virus (HCV) infection might reduce treatment costs and increase the number of patients treated and cured. We determined factors associated with treatment response after an 8-week sofosbuvir-based therapy and developed a simple model to predict an individual's likelihood of treatment success. Among 2907 patients who received ledipasvir/sofosbuvir for 8 weeks, we determined failure rates by demographic and clinical characteristics, and IFNL4-∆G/TT genotype. We estimated the average IFNL4 genotype-related treatment failure rate in major ancestry groups by applying our IFNL4 genotype results to genotype distributions from reference populations. We created a treatment response model based on three personal characteristics. Overall, 4.4% of the patients failed treatment. We observed significantly lower failure rates for persons <50 years (1.6%), females (2.6%), those carrying the IFNL4-TT/TT genotype (1.8%), those with HCV RNA <5.8log10 copies/mL (2.0%) or HCV genotype-1B infection (2.6%). In a model based on ancestry, age and sex, the predicted probability of treatment failure ranged from 0.5% among females of East Asian ancestry <50 years of age to 8.5% among males of African ancestry age ≥65 years. Applying this algorithm at the point-of-care might facilitate HCV elimination, especially in low- and middle-income countries.

Folate Metabolism and Risk of Childhood Acute Lymphoblastic Leukemia: A Genetic Pathway Analysis from the Childhood Cancer and Leukemia International Consortium.

Prenatal folate supplementation has been consistently associated with a reduced risk of childhood acute lymphoblastic leukemia (ALL). Previous germline genetic studies examining the one carbon (folate) metabolism pathway were limited in sample size, scope, and population diversity and led to inconclusive results. We evaluated whether ∼2,900 single-nucleotide polymorphisms (SNP) within 46 candidate genes involved in the folate metabolism pathway influence the risk of childhood ALL, using genome-wide data from nine case-control studies in the Childhood Cancer and Leukemia International Consortium (n = 9,058 cases including 4,510 children of European ancestry, 3,018 Latinx, and 1,406 Asians, and 92,364 controls). Each study followed a standardized protocol for quality control and imputation of genome-wide data and summary statistics were meta-analyzed for all children combined and by major ancestry group using METAL software. None of the selected SNPs reached statistical significance, overall and for major ancestry groups (using adjusted Bonferroni P-value of 5 × 10-6 and less-stringent P-value of 3.5 × 10-5 accounting for the number of "independent" SNPs). None of the 10 top (nonsignificant) SNPs and corresponding genes overlapped across ancestry groups. This large meta-analysis of original data does not reveal associations between many common genetic variants in the folate metabolism pathway and childhood ALL in various ancestry groups. Genetic variants in the folate pathway alone do not appear to substantially influence childhood acute lymphoblastic leukemia risk. Other mechanisms such as gene-folate interaction, DNA methylation, or maternal genetic effects may explain the observed associations with self-reported prenatal folate intake.

Network of large pedigrees reveals social practices of Avar communities.

From AD 567-568, at the onset of the Avar period, populations from the Eurasian Steppe settled in the Carpathian Basin for approximately 250 years1. Extensive sampling for archaeogenomics (424 individuals) and isotopes, combined with archaeological, anthropological and historical contextualization of four Avar-period cemeteries, allowed for a detailed description of the genomic structure of these communities and their kinship and social practices. We present a set of large pedigrees, reconstructed using ancient DNA, spanning nine generations and comprising around 300 individuals. We uncover a strict patrilineal kinship system, in which patrilocality and female exogamy were the norm and multiple reproductive partnering and levirate unions were common. The absence of consanguinity indicates that this society maintained a detailed memory of ancestry over generations. These kinship practices correspond with previous evidence from historical sources and anthropological research on Eurasian Steppe societies2. Network analyses of identity-by-descent DNA connections suggest that social cohesion between communities was maintained via female exogamy. Finally, despite the absence of major ancestry shifts, the level of resolution of our analyses allowed us to detect genetic discontinuity caused by the replacement of a community at one of the sites. This was paralleled with changes in the archaeological record and was probably a result of local political realignment.

Multi-Ancestry Polygenic Risk Score for Coronary Heart Disease Based on an Ancestrally Diverse Genome-Wide Association Study and Population-Specific Optimization.

Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (CHD; PRSCHD) for 5 genetic ancestry groups. We derived ancestry-specific and multi-ancestry PRSCHD based on pruning and thresholding (PRSPT) and ancestry-based continuous shrinkage priors (PRSCSx) applied to summary statistics from the largest multi-ancestry genome-wide association study meta-analysis for CHD to date, including 1.1 million participants from 5 major genetic ancestry groups. Following training and optimization in the Million Veteran Program, we evaluated the best-performing PRSCHD in 176,988 individuals across 9 diverse cohorts. Multi-ancestry PRSPT and PRSCSx outperformed ancestry-specific PRSPT and PRSCSx across a range of tuning values. Two best-performing multi-ancestry PRSCHD (ie, PRSPTmult and PRSCSxmult) and 1 ancestry-specific (PRSCSxEUR) were taken forward for validation. PRSPTmult demonstrated the strongest association with CHD in individuals of South Asian ancestry and European ancestry (odds ratio per 1 SD [95% CI, 2.75 [2.41-3.14], 1.65 [1.59-1.72]), followed by East Asian ancestry (1.56 [1.50-1.61]), Hispanic/Latino ancestry (1.38 [1.24-1.54]), and African ancestry (1.16 [1.11-1.21]). PRSCSxmult showed the strongest associations in South Asian ancestry (2.67 [2.38-3.00]) and European ancestry (1.65 [1.59-1.71]), lower in East Asian ancestry (1.59 [1.54-1.64]), Hispanic/Latino ancestry (1.51 [1.35-1.69]), and the lowest in African ancestry (1.20 [1.15-1.26]). The use of summary statistics from a large multi-ancestry genome-wide meta-analysis improved the performance of PRSCHD in most ancestry groups compared with single-ancestry methods. Despite the use of one of the largest and most diverse sets of training and validation cohorts to date, improvement of predictive performance was limited in African ancestry. This highlights the need for larger genome-wide association study datasets of underrepresented populations to enhance the performance of PRSCHD.

Prevalence of FH-Causing Variants and Impact on LDL-C Concentration in European, South Asian, and African Ancestry Groups of the UK Biobank-Brief Report.

Familial hypercholesterolemia (FH) is a monogenic disease that causes high low-density lipoprotein cholesterol (LDL-C) and higher risk of premature coronary heart disease. The prevalence of FH-causing variants and their association with LDL-C in non-European populations remains largely unknown. Using DNA diagnosis in a population-based cohort, we aimed to estimate the prevalence of FH across 3 major ancestry groups in the United Kingdom. Principal component analysis was used to distinguish genetic ancestry in UK Biobank participants. Whole exome sequencing data were analyzed to provide a genetic diagnosis of FH. LDL-C concentrations were adjusted for statin use. Principal component analysis distinguished 140 439 European, 4067 South Asian, and 3906 African participants with lipid and whole exome sequencing data. There were significant differences between the 3 groups, including total and LDL-C concentrations, and prevalence and incidence of coronary heart disease. We identified 488, 18, and 15 participants of European, South Asian, and African ancestry carrying a likely pathogenic or pathogenic FH-variant. No statistical difference in the prevalence of an FH-causing variant was observed: 1 out of 288 (95% CI, 1/316-1/264) in European, 1 out of 260 (95% CI, 1/526-1/173) in African, and 1 out of 226 (95% CI, 1/419-1/155) in South Asian. Carriers of an FH-causing variant had significantly higher LDL-C concentration than noncarriers in every ancestry group. There was no difference in median (statin-use adjusted) LDL-C concentration in FH-variant carriers depending on their ancestry background. Self-reported statin use was nonsignificantly highest in FH-variant carriers of South Asian ancestry (55.6%), followed by African (40.0%) and European (33.8%; P=0.15). The prevalence of FH-causing variants in the UK Biobank is similar across the ancestry groups analyzed. Despite overall differences in lipid concentrations, FH-variant carriers across the 3 ancestry groups had similar LDL-C levels. In all ancestry groups, the proportion of FH-variant carriers treated with lipid-lowering therapy should be improved to reduce future risk of premature coronary heart disease.

Abstract 4260: Understanding genomic and social determinants of cancer immunotherapy outcome across ancestry

Abstract Compared with previous standards of care, the use of immune checkpoint inhibitors (ICI) has brought significant improvements in survival and quality of life for lung cancer patients. However, only a small proportion of these patients respond durably. People with different ancestries differ probabilistically in genetic factors, environmental exposures, and socio-economic conditions. Whether patients of different ancestry benefit equally from ICIs remains unclear. We studied the impact of genomic ancestry, tumor genomics, and social determinants of health (SDH) factors and factors that are impacted from SDH including recorded race/ethnicity, inferred low-income status from patient zip codes, exposure to smoking, and BMI on ICI response, defined by cancer progression-free survival (PFS, minimum 6 months FU), for non-small cell lung cancer (NSCLC) patients with MSK-IMPACT targeted panel sequencing. This FDA approved assay includes matched tumor-white blood cell sequencing to distinguish germline from somatic variants and has been applied to 1,802 NSCLC patients who received ICI treatment, including 81 and 117 patients with at least 80% of African (AFR) and East Asian (EAS) ancestry, respectively. Moreover, 173 samples were derived from admixed patients with more than one major ancestry. We first used a natural language processing (NLP) model to obtain PFS from free-text clinical notes. A multivariable cox proportional hazard model was then used to associate PFS with ancestry, race, smoking status, ICI drug regimen, PD-L1 status, disease stage, tumor mutational burden (TMB), inferred income, and BMI. Neither genetic ancestry nor self-reported race/ethnicity was associated with the PFS. Moreover, ICI drug regimen types, low-income status, and BMI were not associated with PFS in our cohort. TMB-high was associated with longer PFS across all ancestries, although TMB was lower in patients with EAS ancestry (Median 7.9 vs. 5.3 mut/Mb, p&amp;lt;0.001). These results suggest that the benefits of ICI extend across ancestry, race, and income lines in a single institution, arguing for more equitable patient access to these medications. We also show that TMB is a generalizable biomarker for ICI outcome across ancestries. However, more diverse patient populations are needed to understand whether there is ancestry-specificity in other ICI outcome biomarkers. Citation Format: Christopher J. Fong, Michele Waters, Karl Pichotta, Justin Jee, Devika R. Jutagir, David Ma, Tomin Perea-Chamblee, Susie Kim, Kanika Arora, Brooke Mastrogiacomo, Thinh Tran, Steven Maron, Mirella Altoe, Anisha Luthra, Joseph Kholodenko, Arfath Patha, Doori Rose, Michael F. Berger, Gregory J. Riely, Nikolaus Schultz, Sanna Goyert, Adam Schoenfeld, Francesca Gany, Jian Carrot-Zhang. Understanding genomic and social determinants of cancer immunotherapy outcome across ancestry. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4260.

Genome-wide allele and haplotype-sharing patterns suggested one unique Hmong–Mein-related lineage and biological adaptation history in Southwest China

BackgroundFine-scale genetic structure of ethnolinguistically diverse Chinese populations can fill the gap in the missing diversity and evolutionary landscape of East Asians, particularly for anthropologically informed Chinese minorities. Hmong–Mien (HM) people were one of the most significant indigenous populations in South China and Southeast Asia, which were suggested to be the descendants of the ancient Yangtze rice farmers based on linguistic and archeological evidence. However, their deep population history and biological adaptative features remained to be fully characterized.ObjectivesTo explore the evolutionary and adaptive characteristics of the Miao people, we genotyped genome-wide SNP data in Guizhou HM-speaking populations and merged it with modern and ancient reference populations via a comprehensive population genetic analysis and evolutionary admixture modeling.ResultsThe overall genetic admixture landscape of Guizhou Miao showed genetic differentiation between them and other linguistically diverse Guizhou populations. Admixture models further confirmed that Miao people derived their primary ancestry from geographically close Guangxi Gaohuahua people. The estimated identity by descent and effective population size confirmed a plausible population bottleneck, contributing to their unique genetic diversity and population structure patterns. We finally identified several natural selection candidate genes associated with several biological pathways.ConclusionsGuizhou Miao possessed a specific genetic structure and harbored a close genetic relationship with geographically close southern Chinese indigenous populations and Guangxi historical people. Miao people derived their major ancestry from geographically close Guangxi Gaohuahua people and experienced a plausible population bottleneck which contributed to the unique pattern of their genetic diversity and structure. Future ancient DNA from Shijiahe and Qujialing will provide new insights into the origin of the Miao people.

Development and Validation of a Breast Cancer Polygenic Risk Score on the Basis of Genetic Ancestry Composition.

Polygenic risk scores (PRSs) for breast cancer (BC) risk stratification have been developed primarily in women of European ancestry. Their application to women of non-European ancestry has lagged because of the lack of a formal approach to incorporate genetic ancestry and ancestry-dependent variant frequencies and effect sizes. Here, we propose a multiple-ancestry PRS (MA-PRS) that addresses these issues and may be useful in the development of equitable PRSs across other cancers and common diseases. Women referred for hereditary cancer testing were divided into consecutive cohorts for development (n = 189,230) and for independent validation (n = 89,126). Individual genetic composition as fractions of three reference ancestries (African, East Asian, and European) was determined from ancestry-informative single-nucleotide polymorphisms. The MA-PRS is a combination of three ancestry-specific PRSs on the basis of genetic ancestral composition. Stratification of risk was evaluated by multivariable logistic regression models controlling for family cancer history. Goodness-of-fit analysis compared expected with observed relative risks by quantiles of the MA-PRS distribution. In independent validation, the MA-PRS was significantly associated with BC risk in the full cohort (odds ratio, 1.43; 95% CI, 1.40 to 1.46; P = 8.6 × 10-308) and within each major ancestry. The top decile of the MA-PRS consistently identified patients with two-fold increased risk of developing BC. Goodness-of-fit tests showed that the MA-PRS was well calibrated and predicted BC risk accurately in the tails of the distribution for both European and non-European women. The MA-PRS uses genetic ancestral composition to expand the utility of polygenic risk prediction to non-European women. Inclusion of genetic ancestry in polygenic risk prediction presents an opportunity for more personalized treatment decisions for women of varying and mixed ancestries.

The Genetic Echo of the Tarim Mummies in Modern Central Asians.

The diversity of Central Asians has been shaped by multiple migrations and cultural diffusion. Although ancient DNA studies have revealed the demographic changes of the Central Asian since the Bronze Age, the contribution of the ancient populations to the modern Central Asian remains opaque. Herein, we performed high-coverage sequencing of 131 whole genomes of Indo-European-speaking Tajik and Turkic-speaking Kyrgyz populations to explore their genomic diversity and admixture history. By integrating the ancient DNA data, we revealed more details of the origins and admixture history of Central Asians. We found that the major ancestry of present-day Tajik populations can be traced back to the admixture of the Bronze Age Bactria–Margiana Archaeological Complex and Andronovo-related populations. Highland Tajik populations further received additional gene flow from the Tarim mummies, an isolated ancient North Eurasian–related population. The West Eurasian ancestry of Kyrgyz is mainly derived from Historical Era populations in Xinjiang of China. Furthermore, the recent admixture signals detected in both Tajik and Kyrgyz are ascribed to the expansions of Eastern Steppe nomadic pastoralists during the Historical Era.

The genomic history of southwestern Chinese populations demonstrated massive population migration and admixture among proto-Hmong-Mien speakers and incoming migrants.

Southwest China was the crossroad for the initial settler people of East Asia, which shows the highest diversity in languages and genetics. This region played a significant role in the formation of the genetic makeup of the proto-Hmong-Mien-speaking people and in the north-to-south human expansion during the Neolithic-to-historic transformation. Their genetic history covering migration events and the admixture processes still needs to be further explored. Therefore, in the current study, we have generated genome-wide data from three genomic aspects covering autosomal, mitochondrial and Y-chromosomal regions in 260 Hmong-Mien, Tibeto-Burman, and Sinitic people from 29 different southwestern Chinese groups, and further analyzed them with 2676 published modern and ancient Eurasian genomes. Here, we have noticed a new southwestern East Asian genetic cline composed of the Hmong-Mien-specific ancestry enrichedin modern Hmong and Pathen. This newly identified southern inland East Asian lineage contributed to a great extent of the gene pool in the modern southern East Asians. We also have observed genetic substructure among Hmong-Mien-speaking populations. The southern Hmong-Mien-speaking people showed more genetic affinity with modern Tai-Kadai/Austroasiatic people, while the northern Hmong-Mien speakers expressed a closer genetic connection with the Neolithic-to-modern northern East Asians. Moreover, southwestern Sinitic populations had a strong genomic affinity with the adjacent Hmong-Mien-speaking populations and the lowlander Tibeto-Burman-speaking populations, which suggested the large-scale genetic admixture occurred between them. Allele-sharing-based qpAdm/qpGraph results further confirmed that all included southwestern Chinese populations could be modeled as a mixed result of the major ancestry component from the northern millet farmers in the Yellow River basin and the minor ancestry component from the southern rice farmers in the Yangtze River basin. Usually, this newly identified Hmong-Mien-associated southern East Asian ancestry could improve our understanding of the full-scale genetic landscape of the evolutionary and admixture history of southwestern East Asians. Further ancient genomic studies from southeastern China are required to shed deeper light on our established phylogeny context.

Kouprey (Bos sauveli) genomes unveil polytomic origin of wild Asian Bos

SummaryThe evolution of the genera Bos and Bison, and the nature of gene flow between wild and domestic species, is poorly understood, with genomic data of wild species being limited. We generated two genomes from the likely extinct kouprey (Bos sauveli) and analyzed them alongside other Bos and Bison genomes. We found that B. sauveli possessed genomic signatures characteristic of an independent species closely related to Bos javanicus and Bos gaurus. We found evidence for extensive incomplete lineage sorting across the three species, consistent with a polytomic diversification of the major ancestry in the group, potentially followed by secondary gene flow. Finally, we detected significant gene flow from an unsampled Asian Bos-like source into East Asian zebu cattle, demonstrating both that the full genomic diversity and evolutionary history of the Bos complex has yet to be elucidated and that museum specimens and ancient DNA are valuable resources to do so.

Genomic Insight Into the Population Admixture History of Tungusic-Speaking Manchu People in Northeast China.

Manchu is the third-largest ethnic minority in China and has the largest population size among the Tungusic-speaking groups. However, the genetic origin and admixture history of the Manchu people are far from clear due to the sparse sampling and a limited number of markers genotyped. Here, we provided the first batch of genome-wide data of genotyping approximate 700,000 single-nucleotide polymorphisms (SNPs) in 93 Manchu individuals collected from northeast China. We merged the newly generated data with data of publicly available modern and ancient East Asians to comprehensively characterize the genetic diversity and fine-scale population structure, as well as explore the genetic origin and admixture history of northern Chinese Manchus. We applied both descriptive methods of ADMIXTURE, fineSTRUCTURE, FST, TreeMix, identity by decedent (IBD), principal component analysis (PCA), and qualitative f-statistics (f3, f4, qpAdm, and qpWave). We found that Liaoning Manchus have a close genetic relationship and significant admixture signal with northern Han Chinese, which is in line with the cluster patterns in the haplotype-based results. Additionally, the qpAdm-based admixture models showed that modern Manchu people were formed as major ancestry related to Yellow River farmers and minor ancestry linked to ancient populations from Amur River Bain, or others. In summary, the northeastern Chinese Manchu people in Liaoning were an exception to the coherent genetic structure of Tungusic-speaking populations, probably due to the large-scale population migrations and genetic admixtures in the past few hundred years.

Patterns of genetic connectedness between modern and medieval Estonian genomes reveal the origins of a major ancestry component of the Finnish population

The Finnish population is a unique example of a genetic isolate affected by a recent founder event. Previous studies have suggested that the ancestors of Finnic-speaking Finns and Estonians reached the circum-Baltic region by the 1st millennium BC. However, high linguistic similarity points to a more recent split of their languages. To study genetic connectedness between Finns and Estonians directly, we first assessed the efficacy of imputation of low-coverage ancient genomes by sequencing a medieval Estonian genome to high depth (23×) and evaluated the performance of its down-sampled replicas. We find that ancient genomes imputed from >0.1× coverage can be reliably used in principal-component analyses without projection. By searching for long shared allele intervals (LSAIs; similar to identity-by-descent segments) in unphased data for >143,000 present-day Estonians, 99 Finns, and 14 imputed ancient genomes from Estonia, we find unexpectedly high levels of individual connectedness between Estonians and Finns for the last eight centuries in contrast to their clear differentiation by allele frequencies. High levels of sharing of these segments between Estonians and Finns predate the demographic expansion and late settlement process of Finland. One plausible source of this extensive sharing is the 8th-10th centuries AD migration event from North Estonia to Finland that has been proposed to explain uniquely shared linguistic features between the Finnish language and the northern dialect of Estonian and shared Christianity-related loanwords from Slavic. These results suggest that LSAI detection provides a computationally tractable way to detect fine-scale structure in large cohorts.

Prediction of eye and hair pigmentation phenotypes using the HIrisPlex system in a Brazilian admixed population sample.

Human pigmentation is a complex trait, probably involving more than 100 genes. Predicting phenotypes using SNPs present in those genes is important for forensic purpose. For this, the HIrisPlex tool was developed for eye and hair color prediction, with both models achieving high accuracy among Europeans. Its evaluation in admixed populations is important, since they present a higher frequency of intermediate phenotypes, and HIrisPlex has demonstrated limitations in such predictions; therefore, the performance of this tool may be impaired in such populations. Here, we evaluate the set of 24 markers from the HIrisPlex system in 328 individuals from Ribeirão Preto (SP) region, predicting eye and hair color and comparing the predictions with their real phenotypes. We used the HaloPlex Target Enrichment System and MiSeq Personal Sequencer platform for massively parallel sequencing. The prediction of eye and hair color was accomplished by the HIrisPlex online tool, using the default prediction settings. Ancestry was estimated using the SNPforID 34-plex to observe if and how an individual's ancestry background would affect predictions in this admixed sample. Our sample presented major European ancestry (70.5%), followed by African (21.1%) and Native American/East Asian (8.4%). HIrisPlex presented an overall sensitivity of 0.691 for hair color prediction, with sensitivities ranging from 0.547 to 0.782. The lowest sensitivity was observed for individuals with black hair, who present a reduced European contribution (48.4%). For eye color prediction, the overall sensitivity was 0.741, with sensitivities higher than 0.85 for blue and brown eyes, although it failed in predicting intermediate eye color. Such struggle in predicting this phenotype category is in accordance with what has been seen in previous studies involving HIrisPlex. Individuals with brown eye color are more admixed, with European ancestry decreasing to 62.6%; notwithstanding that, sensitivity for brown eyes was almost 100%. Overall sensitivity increases to 0.791 when a 0.7 threshold is set, though 12.5% of the individuals become undefined. When combining eye and hair prediction, hit rates between 51.3 and 68.9% were achieved. Despite the difficulties with intermediate phenotypes, we have shown that HIrisPlex results can be very helpful when interpreted with caution.

Origins and genetic legacy of prehistoric dogs.

Dogs were the first domestic animal, but little is known about their population history and to what extent it was linked to humans. We sequenced 27 ancient dog genomes and found that all dogs share a common ancestry distinct from present-day wolves, with limited gene flow from wolves since domestication but substantial dog-to-wolf gene flow. By 11,000 years ago, at least five major ancestry lineages had diversified, demonstrating a deep genetic history of dogs during the Paleolithic. Coanalysis with human genomes reveals aspects of dog population history that mirror humans, including Levant-related ancestry in Africa and early agricultural Europe. Other aspects differ, including the impacts of steppe pastoralist expansions in West and East Eurasia and a near-complete turnover of Neolithic European dog ancestry.

Assessment of the potential role of natural selection in type 2 diabetes and related traits across human continental ancestry groups: comparison of phenotypic with genotypic divergence.

Prevalence of type 2 diabetes differs among human ancestry groups, and many hypotheses invoke differential natural selection to account for these differences. We sought to assess the potential role of differential natural selection across major continental ancestry groups for diabetes and related traits, by comparison of genetic and phenotypic differences. This was a cross-sectional comparison among 734 individuals from an urban sample (none of whom was more closely related to another than third-degree relatives), including 83 African Americans, 523 American Indians and 128 European Americans. Participants were not recruited based on diabetes status or other traits. BMI was calculated, and diabetes was diagnosed by a 75g oral glucose tolerance test. In those with normal glucose tolerance (n = 434), fasting insulin and 30min post-load insulin, adjusted for 30min glucose, were taken as measures of insulin resistance and secretion, respectively. Whole exome sequencing was performed, resulting in 97,388 common (minor allele frequency ≥ 5%) variants; the coancestry coefficient (FST) was calculated across all markers as a measure of genetic divergence among ancestry groups. The phenotypic divergence index (PST) was also calculated from the phenotypic differences and heritability (which was estimated from genetic relatedness calculated empirically across all markers in 761 American Indian participants prior to the exclusion of close relatives). Under evolutionary neutrality, the expectation is PST = FST, while for traits under differential selection PST is expected to be significantly greater than FST. A bootstrap procedure was used to test the hypothesis PST = FST. RESULTS: With adjustment for age and sex, prevalence of type 2 diabetes was 34.0% in American Indians, 12.4% in African Americans and 10.4% in European Americans (p = 2.9 × 10-10 for difference among groups). Mean BMI was 36.3, 33.4 and 33.0kg/m2, respectively (p =1.9 × 10-7). Mean fasting insulin was 63.8, 48.4 and 45.2pmol/l (p = 9.2 × 10-5), while mean 30min insulin was 559.8, 553.5 and 358.8pmol/l, respectively (p =5.7 × 10-8). FST across all markers was 0.130, while PST for liability to diabetes, adjusted for age and sex, was 0.149 (p = 0.35 for difference with FST). PST was 0.094 for BMI (p = 0.54), 0.095 for fasting insulin (p = 0.54) and 0.216 (p = 0.18) for 30min insulin. For type 2 diabetes and BMI, the maximum divergence between populations was observed between American Indians and European Americans (PST-MAX = 0.22, p = 0.37, and PST-MAX = 0.14, p = 0.61), which suggests that a relatively modest 22% or 14% of the genetic variance, respectively, can potentially be explained by differential selection (assuming the absence of neutral drift). These analyses suggest that while type 2 diabetes and related traits differ significantly among continental ancestry groups, the differences are consistent with neutral expectations based on heritability and genetic distances. While these analyses do not exclude a modest role for natural selection, they do not support the hypothesis that differential natural selection is necessary to explain the phenotypic differences among these ancestry groups. Graphical abstract.

Common homozygosity for predicted loss-of-function variants reveals both redundant and advantageous effects of dispensable human genes

Humans homozygous or hemizygous for variants predicted to cause a loss of function (LoF) of the corresponding protein do not necessarily present with overt clinical phenotypes. We report here 190 autosomal genes with 207 predicted LoF variants, for which the frequency of homozygous individuals exceeds 1% in at least one human population from five major ancestry groups. No such genes were identified on the X and Y chromosomes. Manual curation revealed that 28 variants (15%) had been misannotated as LoF. Of the 179 remaining variants in 166 genes, only 11 alleles in 11 genes had previously been confirmed experimentally to be LoF. The set of 166 dispensable genes was enriched in olfactory receptor genes (41 genes). The 41 dispensable olfactory receptor genes displayed a relaxation of selective constraints similar to that observed for other olfactory receptor genes. The 125 dispensable nonolfactory receptor genes also displayed a relaxation of selective constraints consistent with greater redundancy. Sixty-two of these 125 genes were found to be dispensable in at least three human populations, suggesting possible evolution toward pseudogenes. Of the 179 LoF variants, 68 could be tested for two neutrality statistics, and 8 displayed robust signals of positive selection. These latter variants included a known FUT2 variant that confers resistance to intestinal viruses, and an APOL3 variant involved in resistance to parasitic infections. Overall, the identification of 166 genes for which a sizeable proportion of humans are homozygous for predicted LoF alleles reveals both redundancies and advantages of such deficiencies for human survival.

An ancestry informative marker panel design for individual ancestry estimation of Hispanic population using whole exome sequencing data

BackgroundEuropeans and American Indians were major genetic ancestry of Hispanics in the U.S. These ancestral groups have markedly different incidence rates and outcomes in many types of cancers. Therefore, the genetic admixture may cause biased genetic association study with cancer susceptibility variants specifically in Hispanics. For example, the incidence rate of liver cancer has been shown with substantial disparity between Hispanic, Asian and non-Hispanic white populations. Currently, ancestry informative marker (AIM) panels have been widely utilized with up to a few hundred ancestry-informative single nucleotide polymorphisms (SNPs) to infer ancestry admixture. Notably, current available AIMs are predominantly located in intron and intergenic regions, while the whole exome sequencing (WES) protocols commonly used in translational research and clinical practice do not cover these markers. Thus, it remains challenging to accurately determine a patient’s admixture proportion without additional DNA testing.ResultsIn this study we designed an unique AIM panel that infers 3-way genetic admixture from three distinct and selective continental populations (African (AFR), European (EUR), and East Asian (EAS)) within evolutionarily conserved exonic regions. Initially, about 1 million exonic SNPs from selective three populations in the 1000 Genomes Project were trimmed by their linkage disequilibrium (LD), restricted to biallelic variants, and finally we optimized to an AIM panel with 250 SNP markers, or the UT-AIM250 panel, using their ancestral informativeness statistics. Comparing to published AIM panels, UT-AIM250 performed better accuracy when we tested with three ancestral populations (accuracy: 0.995 ± 0.012 for AFR, 0.997 ± 0.007 for EUR, and 0.994 ± 0.012 for EAS). We further demonstrated the performance of the UT-AIM250 panel to admixed American (AMR) samples of the 1000 Genomes Project and obtained similar results (AFR, 0.085 ± 0.098; EUR, 0.665 ± 0.182; and EAS, 0.250 ± 0.205) to previously published AIM panels (Phillips-AIM34: AFR, 0.096 ± 0.127, EUR, 0.575 ± 0.290, and EAS, 0.330 ± 0.315; Wei-AIM278: AFR, 0.070 ± 0.096, EUR, 0.537 ± 0.267, and EAS, 0.393 ± 0.300). Subsequently, we applied the UT-AIM250 panel to a clinical dataset of 26 self-reported Hispanic patients in South Texas with hepatocellular carcinoma (HCC). We estimated the admixture proportions using WES data of adjacent non-cancer liver tissues (AFR, 0.065 ± 0.043; EUR, 0.594 ± 0.150; and EAS, 0.341 ± 0.160). Similar admixture proportions were identified from corresponding tumor tissues. In addition, we estimated admixture proportions of The Cancer Genome Atlas (TCGA) collection of hepatocellular carcinoma (TCGA-LIHC) samples (376 patients) using the UT-AIM250 panel. The panel obtained consistent admixture proportions from tumor and matched normal tissues, identified 3 possible incorrectly reported race/ethnicity, and/or provided race/ethnicity determination if necessary.ConclusionsHere we demonstrated the feasibility of using evolutionarily conserved exonic regions to infer admixture proportions and provided a robust and reliable control for sample collection or patient stratification for genetic analysis. R implementation of UT-AIM250 is available at https://github.com/chenlabgccri/UT-AIM250.