RFMix: a discriminative modeling approach for rapid and robust local-ancestry inference

Meta-analysis of genome-wide association studies of stable warfarin dose in patients of African ancestry

ABSTRACT

Warfarin dose requirements are highly variable because of clinical and genetic factors. Although genetic variants influencing warfarin dose have been identified in European and East Asian populations, more work is needed to identify African-specific genetic variants to help optimize warfarin dosing. We performed genome-wide association studies (GWASs) in 4 African cohorts from Uganda, South Africa, and Zimbabwe, totaling 989 warfarin-treated participants who reached stable dose and had international normalized ratios within therapeutic ranges. We also included 2 African American cohorts recruited by the International Warfarin Pharmacogenetics Consortium (n = 316) and the University of Alabama at Birmingham (n = 199). After the GWAS, we performed standard error-weighted meta-analyses and then conducted stepwise conditional analyses to account for known loci in chromosomes 10 and 16. The genome-wide significance threshold was set at P < 5 × 10-8. The meta-analysis, comprising 1504 participants, identified 242 significant SNPs across 3 genomic loci, with 99.6% of these located within known loci on chromosomes 10 (top SNP: rs58800757, P = 4.27 × 10-13) and 16 (top SNP: rs9925964, P = 9.97 × 10-16). Adjustment for the VKORC1 SNP -1639G>A revealed an additional locus on chromosome 2 (top SNPs rs116057875/rs115254730/rs115240773, P = 3.64 × 10-8), implicating the MALL gene, that could indirectly influence warfarin response through interactions with caveolin-1. In conclusion, we reaffirmed the importance of CYP2C9 and VKORC1 in influencing warfarin dose requirements, and identified a new locus (MALL), that still requires direct evidence of biological plausibility.

Joint genotype and ancestry analysis identify novel loci associated with atopic dermatitis in African American population

Atopic dermatitis (AD) is a chronic itchy inflammatory disease of the skin. Genetic studies have identified multiple risk factors linked to the disease; however, most of the studies have been derived from European and East Asian populations. The admixed African American (AA) genome may provide an opportunity to discovery ancestry-specific loci involved in AD susceptibility. Herein, we present joint analysis of ancestry and genotype effects followed by validation using differential gene expression analysis on AD using 726 AD-affected individuals and 999 non-AD control individuals from the AA population, genotyped using Multi-Ethnic Global Array (MEGA) followed by imputation using the Consortium on Asthma among African Ancestry Populations in the Americas (CAAPA) reference panel. The joint analysis identified two novel AD-susceptibility loci, rs2195989 in gene ANGPT1 (8q23.1) and rs62538818 in the intergenic region between genes LURAP1L and MPDZ (9p23). Admixture mapping (AM) results showed potential genomic inflation, and we implemented genomic control and identified five ancestry-of-origin loci with European ancestry effects. The multi-omics functional prioritization of variants in AM signals prioritized the loci SLAIN2, RNF39, and FOXA2. Genome-wide association study (GWAS) identified variants significantly associated with AD in the AA population, including SGK1 (rs113357522, odds ratio [OR] = 2.81), EFR3A (rs16904552, OR = 1.725), and MMP14 (rs911912, OR = 1.791). GWAS variants were common in the AA but rare in the European population, which suggests an African-ancestry-specific risk of AD. Four genes (ANGPT1, LURAP1L, EFR3A, and SGK1) were further validated using qPCR from AD and healthy skin. This study highlighted the importance of genetic studies on admixed populations, as well as local ancestry and genotype-ancestry joint effects to identify risk loci for AD.

Understanding admixture fractions: theory and estimation of gene-flow

Estimation of admixture proportions has become one of the most commonly used computational tools in population genomics. However, there is remarkably little population genetic theory on statistical properties of these variables. We develop theoretical results that can accurately predict means and variances of admixture proportions within a population using models with recombination and genetic drift. Based on established theory on measures of multilocus disequilibrium, we show that there is a set of recurrence relations that can be used to derive expectations for higher moments of the admixture proportions distribution. We obtain closed form solutions for some special cases. Using these results, we develop a method for estimating admixture parameters from estimated admixture proportions obtained from programs such as Structure or Admixture. We apply this method to HapMap 3 data and find that the population history of African Americans, as expected, is not best explained by a single admixture event between people of European and African ancestry. The model of constant gene flow starting at 8 generations and ending at 2 generations before present gives the best fit.

Increased frequency of repeat expansion mutations across different populations

Repeat expansion disorders (REDs) are a devastating group of predominantly neurological diseases. Together they are common, affecting 1 in 3,000 people worldwide with population-specific differences. However, prevalence estimates of REDs are hampered by heterogeneous clinical presentation, variable geographic distributions and technological limitations leading to underascertainment. Here, leveraging whole-genome sequencing data from 82,176 individuals from different populations, we found an overall disease allele frequency of REDs of 1 in 283 individuals. Modeling disease prevalence using genetic data, age at onset and survival, we show that the expected number of people with REDs would be two to three times higher than currently reported figures, indicating underdiagnosis and/or incomplete penetrance. While some REDs are population specific, for example, Huntington disease-like 2 in Africans, most REDs are represented in all broad genetic ancestries (that is, Europeans, Africans, Americans, East Asians and South Asians), challenging the notion that some REDs are found only in specific populations. These results have worldwide implications for local and global health communities in the diagnosis and counseling of REDs.

Genetic disease risks of under-represented founder populations in New York City

The detection of founder pathogenic variants, those observed in high frequency only in a group of individuals with increased inter-relatedness, can help improve delivery of health care for that community. We identified 16 groups with shared ancestry, based on genomic segments that are shared through identity by descent (IBD), in New York City using the genomic data of 25,366 residents from the All Of Us Research Program and the Mount Sinai Bio Me biobank. From these groups we defined 8 as founder populations, mostly communities currently under-represented in medical genomics research, such as Puerto Rican, Garifuna and Filipino/Pacific Islanders. The enrichment analysis of ClinVar pathogenic or likely pathogenic (P/LP) variants in each group identified 202 of these damaging variants across the 8 founder populations. We confirmed disease-causing variants previously reported to occur at increased frequencies in Ashkenazi Jewish and Puerto Rican genetic ancestry groups, but most of the damaging variants identified have not been previously associated with any such founder populations, and most of these founder populations have not been described to have increased prevalence of the associated rare disease. Twenty-five of 51 variants meeting Tier 2 clinical screening criteria (1/100 carrier frequency within these founder groups) have never previously been reported. We show how population structure studies can provide insights into rare diseases disproportionately affecting under-represented founder populations, delivering a health care benefit but also a potential source of stigmatization of these communities, who should be part of the decision-making about implementation into health care delivery.

Author Summary

It is well recognized that genomic studies have been biased towards individuals of European ancestry, and that obtaining medical insights for populations under-represented in medical genomics is crucial to achieve health equity. Here, we use genomic information to identify networks of individuals in New York City who are distinctively related to each other, allowing us to define populations with common genetic ancestry based on genetic similarities rather than by self-reported race or ethnicity. In our study of >25,000 New Yorkers, we identified eight highly-interrelated founder populations, with 202 likely disease-causing variants with increased frequencies in specific founder populations. Many of these population-specific variants are new discoveries, despite their high frequency in founder populations. Studying recent genetic ancestry can help reveal population-specific disease insights that can help with early diagnosis, carrier screening, and opportunities for targeted therapies that all help to reduce health disparities in genomic medicine.

Fast and accurate local ancestry inference with Recomb-Mix

The availability of large genotyped cohorts brings new opportunities for revealing the high-resolution genetic structure of admixed populations via local ancestry inference (LAI), the process of identifying the ancestry of each segment of an individual haplotype. Though current methods achieve high accuracy in standard cases, LAI is still challenging when reference populations are more similar (e.g., intra-continental), when the number of reference populations is too numerous, or when the admixture events are deep in time, all of which are increasingly unavoidable in large biobanks. Here, we present a new LAI method, Recomb-Mix. Recomb-Mix integrates the elements of existing methods of the site-based Li and Stephens model and introduces a new graph collapsing trick to simplify counting paths with the same ancestry label readout. Through comprehensive benchmarking on various simulated datasets, we show that Recomb-Mix is more accurate than existing methods in diverse sets of scenarios while being competitive in terms of resource efficiency. We expect that Recomb-Mix will be a useful method for advancing genetics studies of admixed populations.

Joint models reveal genetic architecture of pubertal stage transitions and their association with BMI in admixed Chilean population

Early or late pubertal onset can lead to disease in adulthood, including cancer, obesity, type 2 diabetes, metabolic disorders, bone fractures, and psychopathologies. Thus, knowing the age at which puberty is attained is crucial as it can serve as a risk factor for future diseases. Pubertal development is divided into five stages of sexual maturation in boys and girls according to the standardized Tanner scale. We performed genome-wide association studies (GWAS) on the "Growth and Obesity Chilean Cohort Study" cohort composed of admixed children with mainly European and Native American ancestry. Using joint models that integrate time-to-event data with longitudinal trajectories of body mass index (BMI), we identified genetic variants associated with phenotypic transitions between pairs of Tanner stages. We identified $42$ novel significant associations, most of them in boys. The GWAS on Tanner $3\rightarrow 4$ transition in boys captured an association peak around the growth-related genes LARS2 and LIMD1 genes, the former of which causes ovarian dysfunction when mutated. The associated variants are expression and splicing Quantitative Trait Loci regulating gene expression and alternative splicing in multiple tissues. Further, higher individual Native American genetic ancestry proportions predicted a significantly earlier puberty onset in boys but not in girls. Finally, the joint models identified a longitudinal BMI parameter significantly associated with several Tanner stages' transitions, confirming the association of BMI with pubertal timing.

A blended genome and exome sequencing method captures genetic variation in an unbiased, high-quality, and cost-effective manner

We deployed the Blended Genome Exome (BGE), a DNA library blending approach that generates low pass whole genome (1-4× mean depth) and deep whole exome (30-40× mean depth) data in a single sequencing run. This technology is cost-effective, empowers most genomic discoveries possible with deep whole genome sequencing, and provides an unbiased method to capture the diversity of common SNP variation across the globe. To evaluate this new technology at scale, we applied BGE to sequence >53,000 samples from the Populations Underrepresented in Mental Illness Associations Studies (PUMAS) Project, which included participants across African, African American, and Latin American populations. We evaluated the accuracy of BGE imputed genotypes against raw genotype calls from the Illumina Global Screening Array. All PUMAS cohorts hadR 2 concordance ≥95% among SNPs with MAF≥1%, and never fell below ≥90%R 2 for SNPs with MAF<1%. Furthermore, concordance rates among local ancestries within two recently admixed cohorts were consistent among SNPs with MAF≥1%, with only minor deviations in SNPs with MAF<1%. We also benchmarked the discovery capacity of BGE to access protein-coding copy number variants (CNVs) against deep whole genome data, finding that deletions and duplications spanning at least 3 exons had a positive predicted value of ~90%. Our results demonstrate BGE scalability and efficacy in capturing SNPs, indels, and CNVs in the human genome at 28% of the cost of deep whole-genome sequencing. BGE is poised to enhance access to genomic testing and empower genomic discoveries, particularly in underrepresented populations.

Gene expression and splicing QTL analysis of blood cells in African American participants from the Jackson Heart Study

Most gene expression and alternative splicing quantitative trait loci (eQTL/sQTL) studies have been biased toward European ancestry individuals. Here, we performed eQTL and sQTL analyses using TOPMed whole-genome sequencing-derived genotype data and RNA-sequencing data from stored peripheral blood mononuclear cells in 1,012 African American participants from the Jackson Heart Study (JHS). At a false discovery rate of 5%, we identified 17,630 unique eQTL credible sets covering 16,538 unique genes; and 24,525 unique sQTL credible sets covering 9,605 unique genes, with lead QTL at P < 5e-8. About 24% of independent eQTLs and independent sQTLs with a minor allele frequency > 1% in JHS were rare (minor allele frequency < 0.1%), and therefore unlikely to be detected, in European ancestry individuals. Finally, we created an open database, which is freely available online, allowing fast query and bulk download of our QTL results.

Genomic Exploration of Essential Hypertension in African-Brazilian Quilombo Populations: A Comprehensive Approach with Pedigree Analysis and Family-Based Association Studies

Background

Essential Hypertension (EH) is a global health issue, responsible for approximately 9.4 million deaths annually. Its prevalence varies by region, with genetic factors contributing 30-60% to blood pressure variation. Despite extensive research, the genetic complexity of EH remains largely unexplained. This study aimed to investigate the genetic basis of EH in African-derived individuals from partially isolated quilombo remnant populations in Vale do Ribeira (SP-Brazil).

Methods

Samples from 431 individuals (167 affected, 261 unaffected, 3 with unknown phenotype) were genotyped using a 650k SNP array. Global ancestry proportions were estimated at 47% African, 36% European, and 16% Native American. Additional data from 673 individuals were used to construct six pedigrees. Pedigrees were pruned, and three non-overlapping marker subpanels were created. We phased haplotypes and performed local ancestry analysis to account for admixture. We then conducted genome-wide linkage analysis (GWLA) and performed fine-mapping through family-based association studies (FBAS) on imputed data and through EH-related genes investigation.

Results

Linkage analysis identified 22 ROIs with LOD scores ranging from 1.45 to 3.03, encompassing 2363 genes. Fine-mapping identified 60 EH-related candidate genes and 118 suggestive or significant variants (FBAS). Among these, 14 genes, including PHGDH, S100A10, MFN2, and RYR2, were strongly associated with hypertension and harbors 29 SNPs.

Conclusions

Through a complementary approach - combining admixture-adjusted genome-wide linkage analysis based on Markov chain Monte Carlo (MCMC) methods, association studies on imputed data, and in silico investigations - genetic regions, variants, and candidate genes were identified, offering insights into the genetic etiology of EH in quilombo remnant populations.

Genome-wide association analysis and admixture mapping in a Puerto Rican cohort supports an Alzheimer disease risk locus on chromosome 12

Introduction

Hispanic/Latino populations are underrepresented in Alzheimer Disease (AD) genetic studies. Puerto Ricans (PR), a three-way admixed (European, African, and Amerindian) population is the second-largest Hispanic group in the continental US. We aimed to conduct a genome-wide association study (GWAS) and comprehensive analyses to identify novel AD susceptibility loci and characterize known AD genetic risk loci in the PR population.

Materials and methods

Our study included Whole Genome Sequencing (WGS) and phenotype data from 648 PR individuals (345 AD, 303 cognitively unimpaired). We used a generalized linear-mixed model adjusting for sex, age, population substructure, and genetic relationship matrix. To infer local ancestry, we merged the dataset with the HGDP/1000G reference panel. Subsequently, we conducted univariate admixture mapping (AM) analysis.

Results

We identified suggestive signals within the SLC38A1 and SCN8A genes on chromosome 12q13. This region overlaps with an area of linkage of AD in previous studies (12q13) in independent data sets further supporting. Univariate African AM analysis identified one suggestive ancestral block (p = 7.2×10-6) located in the same region. The ancestry-aware approach showed that this region has both European and African ancestral backgrounds and both contributing to the risk in this region. We also replicated 11 different known AD loci -including APOE- identified in mostly European studies, which is likely due to the high European background of the PR population.

Conclusion

PR GWAS and AM analysis identified a suggestive AD risk locus on chromosome 12, which includes the SLC38A1 and SCN8A genes. Our findings demonstrate the importance of designing GWAS and ancestry-aware approaches and including underrepresented populations in genetic studies of AD.

The Genomic and Cultural Diversity of the Inka Qhapaq Hucha Ceremony in Chile and Argentina

The South American archaeological record has ample evidence of the socio-cultural dynamism of human populations in the past. This has also been supported through the analysis of ancient genomes, by showing evidence of gene flow across the region. While the extent of these signals is yet to be tested, the growing number of ancient genomes allows for more fine-scaled hypotheses to be evaluated. In this study, we assessed the genetic diversity of individuals associated with the Inka ritual, Qhapaq hucha. As part of this ceremony, one or more individuals were buried with Inka and local-style offerings on mountain summits along the Andes, leaving a very distinctive record. Using paleogenomic tools, we analyzed three individuals: two newly generated genomes from El Plomo Mountain (Chile) and El Toro Mountain (Argentina), and a previously published genome from Argentina (Aconcagua Mountain). Our results reveal a complex demographic scenario with each of the individuals showing different genetic affinities. Furthermore, while two individuals showed genetic similarities with present-day and ancient populations from the southern region of the Inka empire, the third individual may have undertaken long-distance movement. The genetic diversity we observed between individuals from similar cultural contexts supports the highly diverse strategies Inka implemented while incorporating new territories. More broadly, this research contributes to our growing understanding of the population dynamics in the Andes by discussing the implications and temporality of population movements in the region.

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 a specific APOE transcript and functional elements associated with Alzheimer's disease

BACKGROUND

The APOE gene is the strongest genetic risk factor for late-onset Alzheimer's Disease (LOAD). However, the gene regulatory mechanisms at this locus remain incompletely characterized.

METHODS

To identify novel AD-linked functional elements within the APOE locus, we integrated SNP variants with multi-omics data from human postmortem brains including 2,179 RNA-seq samples from 3 brain regions and two ancestries (European and African), 667 DNA methylation samples, and ChIP-seq samples. Additionally, we plotted the expression trajectory of APOE transcripts in human brains during development.

RESULTS

We identified an AD-linked APOE transcript (jxn1.2.2) particularly observed in the dorsolateral prefrontal cortex (DLPFC). The APOE jxn1.2.2 transcript is associated with brain neuropathological features, cognitive impairment, and the presence of the APOE4 allele in DLPFC. We prioritized two independent functional SNPs (rs157580 and rs439401) significantly associated with jxn1.2.2 transcript abundance and DNA methylation levels. These SNPs are located within active chromatin regions and affect brain-related transcription factor-binding affinities. The two SNPs shared effects on the jxn1.2.2 transcript between European and African ethnic groups.

CONCLUSION

The novel APOE functional elements provide potential therapeutic targets with mechanistic insight into the disease etiology.

Biological Insights from Schizophrenia-associated Loci in Ancestral Populations

Large-scale genome-wide association studies of schizophrenia have uncovered hundreds of associated loci but with extremely limited representation of African diaspora populations. We surveyed electronic health records of 200,000 individuals of African ancestry in the Million Veteran and All of Us Research Programs, and, coupled with genotype-level data from four case-control studies, realized a combined sample size of 13,012 affected and 54,266 unaffected persons. Three genome-wide significant signals - near PLXNA4, PMAIP1, and TRPA1 - are the first to be independently identified in populations of predominantly African ancestry. Joint analyses of African, European, and East Asian ancestries across 86,981 cases and 303,771 controls, yielded 376 distinct autosomal loci, which were refined to 708 putatively causal variants via multi-ancestry fine-mapping. Utilizing single-cell functional genomic data from human brain tissue and two complementary approaches, transcriptome-wide association studies and enhancer-promoter contact mapping, we identified a consensus set of 94 genes across ancestries and pinpointed the specific cell types in which they act. We identified reproducible associations of schizophrenia polygenic risk scores with schizophrenia diagnoses and a range of other mental and physical health problems. Our study addresses a longstanding gap in the generalizability of research findings for schizophrenia across ancestral populations, underlining shared biological underpinnings of schizophrenia across global populations in the presence of broadly divergent risk allele frequencies.

Characterizing features affecting local ancestry inference performance in admixed populations

In recent years, significant efforts have been made to improve methods for genomic studies of admixed populations using Local Ancestry Inference (LAI). Accurate LAI is crucial to ensure downstream analyses reflect the genetic ancestry of research participants accurately. Here, we test analytic strategies for LAI to provide guidelines for optimal accuracy, focusing on admixed populations reflective of Latin America's primary continental ancestries - African (AFR), Amerindigenous (AMR), and European (EUR). Simulating LD-informed admixed haplotypes under a variety of 2 and 3-way admixture models, we implemented a standard LAI pipeline, testing three reference panel compositions to quantify their overall and ancestry-specific accuracy. We examined LAI miscall frequencies and true positive rates (TPR) across simulation models and continental ancestries. AMR tracts have notably reduced LAI accuracy as compared to EUR and AFR tracts in all comparisons, with TPR means for AMR ranging from 88-94%, EUR from 96-99% and AFR 98-99%. When LAI miscalls occurred, they most frequently erroneously called European ancestry in true Amerindigenous sites. Using a reference panel well-matched to the target population, even with a lower sample size, LAI produced true-positive estimates that were not statistically different from a high sample size but mismatched reference, while being more computationally efficient. While directly responsive to admixed Latin American cohort compositions, these trends are broadly useful for informing best practices for LAI across other admixed populations. Our findings reinforce the need for inclusion of more underrepresented populations in sequencing efforts to improve reference panels.

Ancestral Diversity in Pharmacogenomics Affects Treatment for Hispanic/Latine Populations With Inflammatory Bowel Disease

BACKGROUND AND AIMS

The prevalence of inflammatory bowel disease among Hispanic/Latine communities is increasing. Pharmacogenomic studies reveal genetic markers that influence treatment decisions. The aim of our study was to examine the frequency and impact of genetic polymorphisms on thiopurine-associated leukopenia (NUDT15, TPMT) and anti-tumor necrosis factor (TNF) immunogenicity (HLA-DQA1∗05) in a cohort of Hispanic patients of diverse ancestral backgrounds.

METHODS

We performed a multicenter, retrospective cohort study comprising 2225 Hispanic participants. We measured the frequency of variation affecting drug response in NUDT15, TPMT, and HLA-DQA1∗05; their ancestral origin (European, African, or Amerindian); and the rate of development of myelosuppression and immunogenicity to thiopurines and anti-TNFs, in exposed patients.

RESULTS

NUDT15 and TPMT variants were rare, except for rs116855232 in NUDT15, which was common only in alleles of Amerindian origin. All NUDT15 variant alleles were inherited on an Amerindian haplotype, and among the Amerindian allele subset, the variant frequency of NUDT15∗4 (rs147390019) was a remarkable 23% in patients with leukopenia but only 3% in patients without leukopenia. HLA-DQA1∗05 and its European tagging variant rs2097432 were common in alleles from all ancestral origins and demonstrated association with immunogenicity to anti-TNFs. However, rs2097432 was only correlated with HLA-DQA1∗05 in the European allele subset.

CONCLUSIONS

These findings indicate that NUDT15 testing should become standard clinical practice before prescribing thiopurines in individuals with Amerindian/Alaska Native ancestry, including Hispanic individuals. Additionally, rs2097432 should not be used as a surrogate for HLA-DQA1∗05 testing for diverse populations. Ultimately, incorporating ancestry in personalized therapeutic approaches is important for treatment of Hispanic patients with inflammatory bowel disease.

Local Ancestry Inference Based on Population-Specific Single-Nucleotide Polymorphisms-A Study of Admixed Populations in the 1000 Genomes Project

Human populations have interacted throughout history, and a considerable portion of modern human populations show evidence of admixture. Local ancestry inference (LAI) is focused on detecting the genetic ancestry of chromosomal segments in admixed individuals and has wide applications. In this work, we proposed a new LAI method based on population-specific single-nucleotide polymorphisms (SNPs) and applied it in the analysis of admixed populations in the 1000 Genomes Project (1KGP). Based on population-specific SNPs in a sliding window, we computed local ancestry information vectors, which are moment estimators of local ancestral proportions, for two haplotypes of an admixed individual and inferred the local ancestral origins. Then we used African (AFR), East Asian (EAS), European (EUR) and South Asian (SAS) populations from the 1KGP and indigenous American (AMR) populations from the Human Genome Diversity Project (HGDP) as reference populations and conducted the proposed LAI analysis on African American populations and American populations in the 1KGP. The results were compared with those obtained by RFMix, G-Nomix and FLARE. We demonstrated that the existence of alleles in a chromosomal region that are specific to a particular reference population and the absence of alleles specific to the other reference populations provide reasonable evidence for determining the ancestral origin of the region. Contemporary AFR, AMR and EUR populations approximate ancestral populations of the admixed populations well, and the results from RFMix, G-Nomix and FLARE largely agree with those from the Ancestral Spectrum Analyzer (ASA), in which the proposed method was implemented. When admixtures are ancient and contemporary reference populations do not satisfactorily approximate ancestral populations, the performances of RFMix, G-Nomix and FLARE deteriorate with increased error rates and fragmented chromosomal segments. In contrast, our method provides fair results.

Urofacial (Ochoa) syndrome with a founder pathogenic variant in the HPSE2 gene: a case report and mutation origin

Urofacial syndrome or Ochoa syndrome (UFS or UFOS) is a rare disease characterized by inverted facial expression and bladder dysfunction that was described for the first time in Colombia. It is an autosomal recessive pathology with mutations in the HPSE2 and LRIG2 genes. However, 16% of patients do not have any mutations associated with the syndrome. Despite the importance of neurobiology in its pathophysiology, there are no neurological, neuropsychological, or psychological studies in these patients. A 30-year-old male from Medellín, Colombia, with a significant perinatal history, was diagnosed with grade 4 hydronephrosis on his first ultrasound test. At 4 months of age, symptoms such as hypomimia, lagophthalmos, and recurrent urinary tract infections started to manifest. Imaging studies revealed urinary tract dilatation, vesicoureteral reflux, and a double collector system on his left side, which led to the diagnosis of UFS. Multiple procedures, including vesicostomy, ureterostomy, and enterocystoplasty, were performed. At 20 years of age, he achieved urinary sphincter control. Genetic analysis revealed a founder pathogenic variant, c.1516C > T (p.Arg506Ter), in the HPSE2 gene, which produces a truncated protein that lacks 86 amino acids. This variant is classified as pathogenic according to the ClinVar database for UFS. The mutation age is approximately 260-360 years, and the two alleles share a 7.2-7.4 Mb IBD segment. Moreover, we detected European local ancestry in the IBD segment, which is consistent with a Spanish introduction. Neurological examination, neuropsychological assessment, and psychological testing revealed no abnormalities, except for high stress levels. Clinical analysis of this patient revealed distorted facial expression and detrusor-sphincter dyssynergia, which are typical of patients with UFS. Genetic analysis revealed a pathogenic variant in the HPSE2 gene of European origin and a mutation age of 260-360 years. From a neurological, neuropsychological, and psychological (emotional and personality) perspective, the patient showed no signs or symptoms of clinical interest.

Association of genetic ancestry with molecular tumor profiles in colorectal cancer

BACKGROUND

There are known disparities in incidence and outcomes of colorectal cancer (CRC) by race and ethnicity. Some of these disparities may be mediated by molecular changes in tumors that occur at different rates across populations. Genetic ancestry is a measure complementary to race and ethnicity that can overcome missing data issues and better capture genetic similarity in admixed populations. We aimed to identify somatic mutations and tumor gene expression differences associated with both genetic ancestry and imputed race and ethnicity.

METHODS

Sequencing was performed with the Tempus xT NGS 648-gene panel and whole exome capture RNA-Seq for 8454 primarily late-stage CRC patients. Genetic ancestry proportions for five continental groups-Africa (AFR), American indigenous (AMR), East Asia (EAS), Europe (EUR), and South Asia (SAS)-were estimated using ancestry informative markers. To address data gaps, race and ethnicity categories were imputed, resulting in assignments for 952 Hispanic/Latino, 420 non-Hispanic (NH) Asian, 1061 NH Black, and 5763 NH White individuals. We assessed association of genetic ancestry proportions and imputed race and ethnicity categories with somatic mutations in relevant CRC genes and in 2608 expression profiles, as well as 1957 consensus molecular subtypes (CMS).

RESULTS

Increased AFR ancestry was associated with higher odds of somatic mutations in APC, KRAS, and PIK3CA and lower odds of BRAF mutations. Additionally, increased EAS ancestry was associated with lower odds of mutations in KRAS, EUR with higher odds in BRAF, and the Hispanic/Latino category with lower odds in BRAF. Greater AFR ancestry and the NH Black category were associated with higher rates of CMS3, while a higher proportion of Hispanic/Latino patients exhibited indeterminate CMS classifications.

CONCLUSIONS

Molecular differences in CRC tumor mutation frequencies and gene expression that may underlie observed differences by race and ethnicity were identified. The association of AFR ancestry with increased KRAS mutations aligns with higher CMS3 subtype rates in NH Black patients. The increase of indeterminate CMS in Hispanic/Latino patients suggests that subtype classification methods could benefit from enhanced patient diversity.

Genomic data reveal a north-south split and introgression history of blood fluke (Schistosoma haematobium) populations from across Africa

The human parasitic fluke, Schistosoma haematobium hybridizes with the livestock parasite S. bovis in the laboratory, but the extent of hybridization in nature is unclear. We analyzed 34.6 million single nucleotide variants in 162 samples from 18 African countries, revealing a sharp genetic discontinuity between northern and southern S. haematobium. We found no evidence for recent hybridization. Instead the data reveal admixture events that occurred 257-879 generations ago in northern S. haematobium populations. Fifteen introgressed S. bovis genes are approaching fixation in northern S. haematobium with four genes potentially driving adaptation. We identified 19 regions that were resistant to introgression; these were enriched on the sex chromosomes. These results (i) demonstrate strong barriers to gene flow between these species, (ii) indicate that hybridization may be less common than currently envisaged, but (iii) reveal profound genomic consequences of interspecific hybridization between schistosomes of medical and veterinary importance.

Genetic determinants and phenotypic consequences of blood T-cell proportions in 207,000 diverse individuals

T-cells play a critical role in multiple aspects of human health and disease. However, to date the genetic determinants of human T-cell abundance have not been studied at scale because assays quantifying T-cell abundance are not widely used in clinical or research settings. The complete blood count clinical assay quantifies lymphocyte abundance which includes T-cells, B-cells, and NK-cells. To address this gap, we directly estimate T-cell fractions from whole genome sequencing data in over 200,000 individuals from the multi-ethnic TOPMed and All of Us studies. We identified 27 loci associated with T-cell fraction. Interrogating electronic health records identified clinical phenotypes associated with T-cell fraction, including notable changes in T-cell proportions that were highly dynamic over the course of pregnancy. In summary, by estimating T-cell fraction, we obtained new insights into the genetic regulation of T-cells and identified disease consequences of T-cell fractions across the human phenome.

Assessing the limits of local ancestry inference from small reference panels

Admixture is a common biological phenomenon among populations of the same or different species. Identifying admixed tracts within individual genomes can provide valuable information to date admixture events, reconstruct ancestry-specific demographic histories, or detect adaptive introgression, genetic incompatibilities, as well as regions of the genomes affected by (associative-) overdominance. Although many local ancestry inference (LAI) methods have been developed in the last decade, their performance was accessed using large reference panels, which are rarely available for non-model organisms or ancient samples. Moreover, the demographic conditions for which LAI becomes unreliable have not been explicitly outlined. Here, we identify the demographic conditions for which local ancestries can be best estimated using very small reference panels. Furthermore, we compare the performance of two LAI methods (RFMix and MOSAIC) with the performance of a newly developed approach (simpLAI) that can be used even when reference populations consist of single individuals. Based on simulations of various demographic models, we also determine the limits of these LAI tools and propose post-painting filtering steps to reduce false-positive rates and improve the precision and accuracy of the inferred admixed tracts. Besides providing a guide for using LAI, our work shows that reasonable inferences can be obtained from a single diploid genome per reference under demographic conditions that are not uncommon among past human groups and non-model organisms.

Extended genome-wide association study employing the African genome resources panel identifies novel susceptibility loci for Alzheimer's disease in individuals of African ancestry

INTRODUCTION

Despite a two-fold risk, individuals of African ancestry have been underrepresented in Alzheimer's disease (AD) genomics efforts.

METHODS

Genome-wide association studies (GWAS) of 2,903 AD cases and 6,265 controls of African ancestry. Within-dataset results were meta-analyzed, followed by functional genomics analyses.

RESULTS

A novel AD-risk locus was identified in MPDZ on chromosome (chr) 9p23 (rs141610415, MAF = 0.002, p = 3.68×10-9). Two additional novel common and nine rare loci were identified with suggestive associations (P < 9×10-7). Comparison of association and linkage disequilibrium (LD) patterns between datasets with higher and lower degrees of African ancestry showed differential association patterns at chr12q23.2 (ASCL1), suggesting that this association is modulated by regional origin of local African ancestry.

DISCUSSION

These analyses identified novel AD-associated loci in individuals of African ancestry and suggest that degree of African ancestry modulates some associations. Increased sample sets covering as much African genetic diversity as possible will be critical to identify additional loci and deconvolute local genetic ancestry effects.

HIGHLIGHTS

Genetic ancestry significantly impacts risk of Alzheimer's Disease (AD). Although individuals of African ancestry are twice as likely to develop AD, they are vastly underrepresented in AD genomics studies. The Alzheimer's Disease Genetics Consortium has previously identified 16 common and rare genetic loci associated with AD in African American individuals. The current analyses significantly expand this effort by increasing the sample size and extending ancestral diversity by including populations from continental Africa. Single variant meta-analysis identified a novel genome-wide significant AD-risk locus in individuals of African ancestry at the MPDZ gene, and 11 additional novel loci with suggestive genome-wide significance at p < 9×10-7. Comparison of African American datasets with samples of higher degree of African ancestry demonstrated differing patterns of association and linkage disequilibrium at one of these loci, suggesting that degree and/or geographic origin of African ancestry modulates the effect at this locus. These findings illustrate the importance of increasing number and ancestral diversity of African ancestry samples in AD genomics studies to fully disentangle the genetic architecture underlying AD, and yield more effective ancestry-informed genetic screening tools and therapeutic interventions.

Exploring the relationship between admixture and genetic susceptibility to attention deficit hyperactivity disorder in two Latin American cohorts

Contemporary research on the genomics of Attention Deficit Hyperactivity Disorder (ADHD) often underrepresents admixed populations of diverse genomic ancestries, such as Latin Americans. This study explores the relationship between admixture and genetic associations for ADHD in Colombian and Mexican cohorts. Some 546 participants in two groups, ADHD and Control, were genotyped with Infinium PsychArray®. Global ancestry levels were estimated using overall admixture proportions and principal component analysis, while local ancestry was determined using a method to estimate ancestral components along the genome. Genome-wide association analysis (GWAS) was conducted to identify significant associations. Differences between Colombia and Mexico were evaluated using appropriate statistical tests. 354 Single-nucleotide polymorphisms (SNPs) and Single-nucleotide variants (SNVs) related to some genes and intergenic regions exhibited suggestive significance (p-value < 5*10e-5) in the GWAS. None of the variants revealed genome-wide significance (p-value < 5*10e-8). The study identified a significant relationship between risk SNPs and the European component of admixture, notably observed in the LOC105379109 gene. Despite differences in risk association loci, such as FOXP2, our findings suggest a possible homogeneity in genetic variation's impact on ADHD between Colombian and Mexican populations. Current reference datasets for ADHD predominantly consist of samples with high European ancestry, underscoring the need for further research to enhance the representation of reference populations and improve the identification of ADHD risk traits in Latin Americans.

RFMix-reader: Accelerated reading and processing for local ancestry studies

Motivation

Local ancestry inference is a powerful technique in genetics, revealing population history and the genetic basis of diseases. It is particularly valuable for improving eQTL discovery and fine-mapping in admixed populations. Despite the widespread use of the RFMix software for local ancestry inference, large-scale genomic studies face challenges of high memory consumption and processing times when handling RFMix output files.

Results

Here, I present RFMix-reader , a new Python-based parsing software, designed to streamline the analysis of large-scale local ancestry datasets. This software prioritizes computational eiciency and memory optimization, leveraging GPUs when available for additional speed boosts. By overcoming these data processing hurdles, RFMix-reader empowers researchers to unlock the full potential of local ancestry data for understanding human health and health disparities.

Availability

RFMix-reader is freely available on PyPI at https://pypi.org/project/rfmix-reader/ , implemented in Python 3, and supported on Linux, Windows, and Mac OS.

Contact

KynonJade.Benjamin@libd.org.

Supplementary information

Supplementary data are available at https://rfmix-reader.readthedocs.io/en/latest/ .

Genetic Ancestry and Self-Reported "Skin Color/Race" in the Urban Admixed Population of São Paulo City, Brazil

Epidemiological studies frequently classify groups based on phenotypes like self-reported skin color/race, which inaccurately represent genetic ancestry and may lead to misclassification, particularly among individuals of multiracial backgrounds. This study aimed to characterize both global and local genome-wide genetic ancestries and to assess their relationship with self-reported skin color/race in an admixed population of Sao Paulo city. We analyzed 226,346 single-nucleotide polymorphisms from 841 individuals participating in the population-based ISA-Nutrition study. Our findings confirmed the admixed nature of the population, demonstrating substantial European, significant Sub-Saharan African, and minor Native American ancestries, irrespective of skin color. A correlation was observed between global genetic ancestry and self-reported color-race, which was more evident in the extreme proportions of African and European ancestries. Individuals with higher African ancestry tended to identify as Black, those with higher European ancestry tended to identify as White, and individuals with higher Native American ancestry were more likely to self-identify as Mixed, a group with diverse ancestral compositions. However, at the individual level, this correlation was notably weak, and no deviations were observed for specific regions throughout the individual's genome. Our findings emphasize the significance of accurately defining and thoroughly analyzing race and ancestry, especially within admixed populations.

Molecular and clinical characterization of a founder mutation causing G6PC3 deficiency

G6PC3 deficiency is a monogenic immunometabolic disorder that causes syndromic congenital neutropenia. Patients display heterogeneous extra-hematological manifestations, contributing to delayed diagnosis. Here, we investigated the origin and functional consequence of the G6PC3 c.210delC variant found in patients of Mexican origin. Based on the shared haplotypes amongst carriers of the c.210delC mutation, we estimated that this variant originated from a founder effect in a common ancestor. Furthermore, by ancestry analysis, we concluded that it originated in the indigenous Mexican population. At the protein level, we showed that this frameshift mutation leads to an aberrant protein expression in overexpression and patient-derived cells. G6PC3 pathology is driven by the intracellular accumulation of the metabolite 1,5-anhydroglucitol-6-phosphate (1,5-AG6P) that inhibits glycolysis. We characterized how the variant c.210delC impacts glycolysis by performing extracellular flux assays on patient-derived cells. When treated with 1,5-anhydroglucitol (1,5-AG), the precursor to 1,5-AG6P, patient-derived cells exhibited markedly reduced engagement of glycolysis. Finally, we compared the clinical presentation of patients with the mutation c.210delC and all other G6PC3 deficient patients reported in the literature to date, and we found that c.210delC carriers display all prominent clinical features observed in prior G6PC3 deficient patients. In conclusion, G6PC3 c.210delC is a loss-of-function mutation that arose from a founder effect in the indigenous Mexican population. These findings may facilitate the diagnosis of additional patients in this geographical area. Moreover, the in vitro 1,5-AG-dependent functional assay used in our study could be employed to assess the pathogenicity of additional G6PC3 variants.

Increased frequency of repeat expansion mutations across different populations

Repeat expansion disorders (REDs) are a devastating group of predominantly neurological diseases. Together they are common, affecting 1 in 3,000 people worldwide with population-specific differences. However, prevalence estimates of REDs are hampered by heterogeneous clinical presentation, variable geographic distributions, and technological limitations leading to under-ascertainment. Here, leveraging whole genome sequencing data from 82,176 individuals from different populations, we found an overall disease allele frequency of REDs of 1 in 283 individuals. Modelling disease prevalence using genetic data, age at onset and survival, we show that the expected number of people with REDs would be two to three times higher than currently reported figures, indicating under-diagnosis and/or incomplete penetrance. While some REDs are population-specific, e.g. Huntington disease-like 2 in Africans, most REDs are represented in all broad genetic ancestries (i.e. Europeans, Africans, Americans, East Asians, and South Asians), challenging the notion that some REDs are found only in specific populations. These results have worldwide implications for local and global health communities in the diagnosis and counselling of REDs.

A genome-wide association study of alloimmunization in the TOPMed OMG-SCD cohort identifies a locus on chromosome 12

BACKGROUND

Red cell alloimmunization after exposure to donor red cells is a very common complication of transfusion for patients with sickle cell disease (SCD), resulting frequently in accelerated donor red blood cell destruction. Patients show substantial differences in their predisposition to alloimmunization, and genetic variability is one proposed component. Although several genetic association studies have been conducted for alloimmunization, the results have been inconsistent, and the genetic determinants of alloimmunization remain largely unknown.

STUDY DESIGN AND METHODS

We performed a genome-wide association study (GWAS) in 236 African American (AA) SCD patients from the Outcome Modifying Genes in Sickle Cell Disease (OMG-SCD) cohort, which is part of Trans-Omics for Precision Medicine (TOPMed), with whole-genome sequencing data available. We also performed sensitivity analyses adjusting for different sets of covariates and applied different sample grouping strategies based on the number of alloantibodies patients developed.

RESULTS

We identified one genome-wide significant locus on chr12 (p = 3.1e-9) with no evidence of genomic inflation (lambda = 1.003). Further leveraging QTL evidence from GTEx whole blood and/or Jackson Heart Study PBMC RNA-Seq data, we identified a number of potential genes, such as ARHGAP9, STAT6, and ATP23, that may be driving the association signal. We also discovered some suggestive loci using different analysis strategies.

DISCUSSION

We call for the community to collect additional alloantibody information within SCD cohorts to further the understanding of the genetic basis of alloimmunization in order to improve transfusion outcomes.

Ancestral Origins and Admixture History of Kazakhs

Kazakh people, like many other populations that settled in Central Asia, demonstrate an array of mixed anthropological features of East Eurasian (EEA) and West Eurasian (WEA) populations, indicating a possible scenario of biological admixture between already differentiated EEA and WEA populations. However, their complex biological origin, genomic makeup, and genetic interaction with surrounding populations are not well understood. To decipher their genetic structure and population history, we conducted, to our knowledge, the first whole-genome sequencing study of Kazakhs residing in Xinjiang (KZK). We demonstrated that KZK derived their ancestries from 4 ancestral source populations: East Asian (∼39.7%), West Asian (∼28.6%), Siberian (∼23.6%), and South Asian (∼8.1%). The recognizable interactions of EEA and WEA ancestries in Kazakhs were dated back to the 15th century BCE. Kazakhs were genetically distinctive from the Uyghurs in terms of their overall genomic makeup, although the 2 populations were closely related in genetics, and both showed a substantial admixture of western and eastern peoples. Notably, we identified a considerable sex-biased admixture, with an excess of western males and eastern females contributing to the KZK gene pool. We further identified a set of genes that showed remarkable differentiation in KZK from the surrounding populations, including those associated with skin color (SLC24A5, OCA2), essential hypertension (HLA-DQB1), hypertension (MTHFR, SLC35F3), and neuron development (CNTNAP2). These results advance our understanding of the complex history of contacts between Western and Eastern Eurasians, especially those living or along the old Silk Road.

Global and local ancestry estimation in a captive baboon colony

The last couple of decades have highlighted the importance of studying hybridization, particularly among primate species, as it allows us to better understand our own evolutionary trajectory. Here, we report on genetic ancestry estimates using dense, full genome data from 881 olive (Papio anubus), yellow (Papio cynocephalus), or olive-yellow crossed captive baboons from the Southwest National Primate Research Center. We calculated global and local ancestry information, imputed low coverage genomes (n = 830) to improve marker quality, and updated the genetic resources of baboons available to assist future studies. We found evidence of historical admixture in some putatively purebred animals and identified errors within the Southwest National Primate Research Center pedigree. We also compared the outputs between two different phasing and imputation pipelines along with two different global ancestry estimation software. There was good agreement between the global ancestry estimation software, with R2 > 0.88, while evidence of phase switch errors increased depending on what phasing and imputation pipeline was used. We also generated updated genetic maps and created a concise set of ancestry informative markers (n = 1,747) to accurately obtain global ancestry estimates.

Inference of Locus-Specific Population Mixtures from Linked Genome-Wide Allele Frequencies

Admixture between populations and species is common in nature. Since the influx of new genetic material might be either facilitated or hindered by selection, variation in mixture proportions along the genome is expected in organisms undergoing recombination. Various graph-based models have been developed to better understand these evolutionary dynamics of population splits and mixtures. However, current models assume a single mixture rate for the entire genome and do not explicitly account for linkage. Here, we introduce TreeSwirl, a novel method for inferring branch lengths and locus-specific mixture proportions by using genome-wide allele frequency data, assuming that the admixture graph is known or has been inferred. TreeSwirl builds upon TreeMix that uses Gaussian processes to estimate the presence of gene flow between diverged populations. However, in contrast to TreeMix, our model infers locus-specific mixture proportions employing a hidden Markov model that accounts for linkage. Through simulated data, we demonstrate that TreeSwirl can accurately estimate locus-specific mixture proportions and handle complex demographic scenarios. It also outperforms related D- and f-statistics in terms of accuracy and sensitivity to detect introgressed loci.

Admixture mapping of cognitive function in diverse Hispanic and Latino adults: Results from the Hispanic Community Health Study/Study of Latinos

INTRODUCTION

We conducted admixture mapping and fine-mapping analyses to identify ancestry-of-origin loci influencing cognitive abilities.

METHODS

We estimated the association of local ancestry intervals across the genome with five neurocognitive measures in 7140 diverse Hispanic and Latino adults (mean age 55 years). We prioritized genetic variants in associated loci and tested them for replication in four independent cohorts.

RESULTS

We identified nine local ancestry-associated regions for the five neurocognitive measures. There was strong biological support for the observed associations to cognitive function at all loci and there was statistical evidence of independent replication at 4q12, 9p22.1, and 13q12.13.

DISCUSSION

Our study identified multiple novel loci harboring genes implicated in cognitive functioning and dementia, and uncovered ancestry-relevant genetic variants. It adds to our understanding of the genetic architecture of cognitive function in Hispanic and Latino adults and demonstrates the power of admixture mapping to discover unique haplotypes influencing cognitive function, complementing genome-wide association studies.

HIGHLIGHTS

We identified nine ancestry-of-origin chromosomal regions associated with five neurocognitive traits. In each associated region, we identified single nucleotide polymorphisms (SNPs) that explained, at least in part, the admixture signal and were tested for replication in independent samples of Black, non-Hispanic White, and Hispanic/Latino adults with the same or similar neurocognitive tests. Statistical evidence of independent replication of the prioritized SNPs was observed for three of the nine associations, at chr4q12, chr9p22.1, and chr13q12.13. At all loci, there was strong biological support for the observed associations to cognitive function and dementia, prioritizing genes such as KIT, implicated in autophagic clearance of neurotoxic proteins and on mast cell and microglial-mediated inflammation; SLC24A2, implicated in synaptic plasticity associated with learning and memory; and MTMR6, implicated in phosphoinositide lipids metabolism.

A harmonized public resource of deeply sequenced diverse human genomes

Underrepresented populations are often excluded from genomic studies owing in part to a lack of resources supporting their analyses. The 1000 Genomes Project (1kGP) and Human Genome Diversity Project (HGDP), which have recently been sequenced to high coverage, are valuable genomic resources because of the global diversity they capture and their open data sharing policies. Here, we harmonized a high-quality set of 4094 whole genomes from 80 populations in the HGDP and 1kGP with data from the Genome Aggregation Database (gnomAD) and identified over 153 million high-quality SNVs, indels, and SVs. We performed a detailed ancestry analysis of this cohort, characterizing population structure and patterns of admixture across populations, analyzing site frequency spectra, and measuring variant counts at global and subcontinental levels. We also show substantial added value from this data set compared with the prior versions of the component resources, typically combined via liftOver and variant intersection; for example, we catalog millions of new genetic variants, mostly rare, compared with previous releases. In addition to unrestricted individual-level public release, we provide detailed tutorials for conducting many of the most common quality-control steps and analyses with these data in a scalable cloud-computing environment and publicly release this new phased joint callset for use as a haplotype resource in phasing and imputation pipelines. This jointly called reference panel will serve as a key resource to support research of diverse ancestry populations.

Germline variation contributes to false negatives in CRISPR-based experiments with varying burden across ancestries

Reducing disparities is vital for equitable access to precision treatments in cancer. Socioenvironmental factors are a major driver of disparities, but differences in genetic variation likely also contribute. The impact of genetic ancestry on prioritization of cancer targets in drug discovery pipelines has not been systematically explored due to the absence of pre-clinical data at the appropriate scale. Here, we analyze data from 611 genome-scale CRISPR/Cas9 viability experiments in human cell line models to identify ancestry-associated genetic dependencies essential for cell survival. Surprisingly, we find that most putative associations between ancestry and dependency arise from artifacts related to germline variants. Our analysis suggests that for 1.2-2.5% of guides, germline variants in sgRNA targeting sequences reduce cutting by the CRISPR/Cas9 nuclease, disproportionately affecting cell models derived from individuals of recent African descent. We propose three approaches to mitigate this experimental bias, enabling the scientific community to address these disparities.

Latin American Trans-ancestry INitiative for OCD genomics (LATINO): Study protocol

Obsessive-compulsive disorder (OCD) is a debilitating psychiatric disorder. Worldwide, its prevalence is ~2% and its etiology is mostly unknown. Identifying biological factors contributing to OCD will elucidate underlying mechanisms and might contribute to improved treatment outcomes. Genomic studies of OCD are beginning to reveal long-sought risk loci, but >95% of the cases currently in analysis are of homogenous European ancestry. If not addressed, this Eurocentric bias will result in OCD genomic findings being more accurate for individuals of European ancestry than other ancestries, thereby contributing to health disparities in potential future applications of genomics. In this study protocol paper, we describe the Latin American Trans-ancestry INitiative for OCD genomics (LATINO, https://www.latinostudy.org). LATINO is a new network of investigators from across Latin America, the United States, and Canada who have begun to collect DNA and clinical data from 5000 richly phenotyped OCD cases of Latin American ancestry in a culturally sensitive and ethical manner. In this project, we will utilize trans-ancestry genomic analyses to accelerate the identification of OCD risk loci, fine-map putative causal variants, and improve the performance of polygenic risk scores in diverse populations. We will also capitalize on rich clinical data to examine the genetics of treatment response, biologically plausible OCD subtypes, and symptom dimensions. Additionally, LATINO will help elucidate the diversity of the clinical presentations of OCD across cultures through various trainings developed and offered in collaboration with Latin American investigators. We believe this study will advance the important goal of global mental health discovery and equity.

Analysis of gene expression in the postmortem brain of neurotypical Black Americans reveals contributions of genetic ancestry

Ancestral differences in genomic variation affect the regulation of gene expression; however, most gene expression studies have been limited to European ancestry samples or adjusted to identify ancestry-independent associations. Here, we instead examined the impact of genetic ancestry on gene expression and DNA methylation in the postmortem brain tissue of admixed Black American neurotypical individuals to identify ancestry-dependent and ancestry-independent contributions. Ancestry-associated differentially expressed genes (DEGs), transcripts and gene networks, while notably not implicating neurons, are enriched for genes related to the immune response and vascular tissue and explain up to 26% of heritability for ischemic stroke, 27% of heritability for Parkinson disease and 30% of heritability for Alzheimer's disease. Ancestry-associated DEGs also show general enrichment for the heritability of diverse immune-related traits but depletion for psychiatric-related traits. We also compared Black and non-Hispanic white Americans, confirming most ancestry-associated DEGs. Our results delineate the extent to which genetic ancestry affects differences in gene expression in the human brain and the implications for brain illness risk.

Data Harmonization Guidelines to Combine Multi-platform Genomic Data from Admixed Populations and Boost Power in Genome-Wide Association Studies

Data harmonization involves combining data from multiple independent sources and processing the data to produce one uniform dataset. Merging separate genotypes or whole-genome sequencing datasets has been proposed as a strategy to increase the statistical power of association tests by increasing the effective sample size. However, data harmonization is not a widely adopted strategy due to the difficulties with merging data (including confounding produced by batch effects and population stratification). Detailed data harmonization protocols are scarce and are often conflicting. Moreover, data harmonization protocols that accommodate samples of admixed ancestry are practically non-existent. Existing data harmonization procedures must be modified to ensure the heterogeneous ancestry of admixed individuals is incorporated into additional downstream analyses without confounding results. Here, we propose a set of guidelines for merging multi-platform genetic data from admixed samples that can be adopted by any investigator with elementary bioinformatics experience. We have applied these guidelines to aggregate 1544 tuberculosis (TB) case-control samples from six separate in-house datasets and conducted a genome-wide association study (GWAS) of TB susceptibility. The GWAS performed on the merged dataset had improved power over analyzing the datasets individually and produced summary statistics free from bias introduced by batch effects and population stratification. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Processing separate datasets comprising array genotype data Alternate Protocol 1: Processing separate datasets comprising array genotype and whole-genome sequencing data Alternate Protocol 2: Performing imputation using a local reference panel Basic Protocol 2: Merging separate datasets Basic Protocol 3: Ancestry inference using ADMIXTURE and RFMix Basic Protocol 4: Batch effect correction using pseudo-case-control comparisons.

Recent selection and introgression facilitated high-altitude adaptation in cattle

During the past 3000 years, cattle on the Qinghai-Xizang Plateau have developed adaptive phenotypes under the selective pressure of hypoxia, ultraviolet (UV) radiation, and extreme cold. The genetic mechanism underlying this rapid adaptation is not yet well understood. Here, we present whole-genome resequencing data for 258 cattle from 32 cattle breeds/populations, including 89 Tibetan cattle representing eight populations distributed at altitudes ranging from 3400 m to 4300 m. Our genomic analysis revealed that Tibetan cattle exhibited a continuous phylogeographic cline from the East Asian taurine to the South Asian indicine ancestries. We found that recently selected genes in Tibetan cattle were related to body size (HMGA2 and NCAPG) and energy expenditure (DUOXA2). We identified signals of sympatric introgression from yak into Tibetan cattle at different altitudes, covering 0.64%-3.26% of their genomes, which included introgressed genes responsible for hypoxia response (EGLN1), cold adaptation (LRP11), DNA damage repair (LATS1), and UV radiation resistance (GNPAT). We observed that introgressed yak alleles were associated with noncoding variants, including those in present EGLN1. In Tibetan cattle, three yak introgressed SNPs in the EGLN1 promoter region reduced the expression of EGLN1, suggesting that these genomic variants enhance hypoxia tolerance. Taken together, our results indicated complex adaptation processes in Tibetan cattle, where recently selected genes and introgressed yak alleles jointly facilitated rapid adaptation to high-altitude environments.

Genetic determinants of serum bilirubin using inferred native American gene variants in Chilean adolescents

Gene variants in the UGT1A1 gene are strongly associated with circulating bilirubin levels in several populations, as well as other variants of modest effect across the genome. However, the effects of such variants are unknown regarding the Native American ancestry of the admixed Latino population. Our objective was to assess the Native American genetic determinants of serum bilirubin in Chilean admixed adolescents using the local ancestry deconvolution approach. We measured total serum bilirubin levels in 707 adolescents of the Chilean Growth and Obesity Cohort Study (GOCS) and performed high-density genotyping using the Illumina-MEGA array (>1.7 million genotypes). We constructed a local ancestry reference panel with participants from the 1000 Genomes Project, the Human Genome Diversity Project, and our GOCS cohort. Then, we inferred and isolated haplotype tracts of Native American, European, or African origin to perform genome-wide association studies. In the whole cohort, the rs887829 variant and others near UGT1A1 were the unique signals achieving genome-wide statistical significance (b = 0.30; p = 3.34 × 10-57). After applying deconvolution methods, we found that significance is also maintained in Native American (b = 0.35; p = 3.29 × 10-17) and European (b = 0.28; p = 1.14 × 10-23) ancestry components. The rs887829 variant explained a higher percentage of the variance of bilirubin in the Native American (37.6%) compared to European ancestry (28.4%). In Native American ancestry, carriers of the TT genotype of this variant averaged 4-fold higher bilirubinemia compared to the CC genotype (p = 2.82 × 10-12). We showed for the first time that UGT1A1 variants are the primary determinant of bilirubin levels in Native American ancestry, confirming its pan-ethnic relevance. Our study illustrates the general value of the local ancestry deconvolution approach to assessing isolated ancestry effects in admixed populations.

scAI-SNP: a method for inferring ancestry from single-cell data

Collaborative efforts, such as the Human Cell Atlas, are rapidly accumulating large amounts of single-cell data. To ensure that single-cell atlases are representative of human genetic diversity, we need to determine the ancestry of the donors from whom single-cell data are generated. Self-reporting of race and ethnicity, although important, can be biased and is not always available for the datasets already collected. Here, we introduce scAI-SNP, a tool to infer ancestry directly from single-cell genomics data. To train scAI-SNP, we identified 4.5 million ancestry-informative single-nucleotide polymorphisms (SNPs) in the 1000 Genomes Project dataset across 3201 individuals from 26 population groups. For a query single-cell data set, scAI-SNP uses these ancestry-informative SNPs to compute the contribution of each of the 26 population groups to the ancestry of the donor from whom the cells were obtained. Using diverse single-cell data sets with matched whole-genome sequencing data, we show that scAI-SNP is robust to the sparsity of single-cell data, can accurately and consistently infer ancestry from samples derived from diverse types of tissues and cancer cells, and can be applied to different modalities of single-cell profiling assays, such as single-cell RNA-seq and single-cell ATAC-seq. Finally, we argue that ensuring that single-cell atlases represent diverse ancestry, ideally alongside race and ethnicity, is ultimately important for improved and equitable health outcomes by accounting for human diversity.

Molecular and clinical characterization of a founder mutation causing G6PC3 deficiency

Background

G6PC3 deficiency is a rare genetic disorder that causes syndromic congenital neutropenia. It is driven by the intracellular accumulation of a metabolite named 1,5-anhydroglucitol-6-phosphate (1,5-AG6P) that inhibits glycolysis. Patients display heterogeneous extra-hematological manifestations, contributing to delayed diagnosis.

Objective

The G6PC3 c.210delC variant has been identified in patients of Mexican origin. We set out to study the origin and functional consequence of this mutation. Furthermore, we sought to characterize the clinical phenotypes caused by it.

Methods

Using whole-genome sequencing data, we conducted haplotype analysis to estimate the age of this allele and traced its ancestral origin. We examined how this mutation affected G6PC3 protein expression and performed extracellular flux assays on patient-derived cells to characterize how this mutation impacts glycolysis. Finally, we compared the clinical presentations of patients with the c.210delC mutation relative to other G6PC3 deficient patients published to date.

Results

Based on the length of haplotypes shared amongst ten carriers of the G6PC3 c.210delC mutation, we estimated that this variant originated in a common ancestor of indigenous American origin. The mutation causes a frameshift that introduces a premature stop codon, leading to a complete loss of G6PC3 protein expression. When treated with 1,5-anhydroglucitol (1,5-AG), the precursor to 1,5-AG6P, patient-derived cells exhibited markedly reduced engagement of glycolysis. Clinically, c.210delC carriers display all the clinical features of syndromic severe congenital neutropenia type 4 observed in prior reports of G6PC3 deficiency.

Conclusion

The G6PC3 c.210delC is a loss-of-function mutation that arose from a founder effect in the indigenous Mexican population. These findings may facilitate the diagnosis of additional patients in this geographical area. Moreover, the in vitro 1,5-AG-dependent functional assay used in our study could be employed to assess the pathogenicity of additional G6PC3 variants.

Characterizing substructure via mixture modeling in large-scale genetic summary statistics

Genetic summary data are broadly accessible and highly useful including for risk prediction, causal inference, fine mapping, and incorporation of external controls. However, collapsing individual-level data into groups masks intra- and inter-sample heterogeneity, leading to confounding, reduced power, and bias. Ultimately, unaccounted substructure limits summary data usability, especially for understudied or admixed populations. Here, we present Summix2, a comprehensive set of methods and software based on a computationally efficient mixture model to estimate and adjust for substructure in genetic summary data. In extensive simulations and application to public data, Summix2 characterizes finer-scale population structure, identifies ascertainment bias, and identifies potential regions of selection due to local substructure deviation. Summix2 increases the robust use of diverse publicly available summary data resulting in improved and more equitable research.

Phasing millions of samples achieves near perfect accuracy, enabling parent-of-origin classification of variants

Haplotype phasing, the process of determining which genetic variants are physically located on the same chromosome, is crucial for various genetic analyses. In this study, we first benchmark SHAPEIT and Beagle, two state-of-the-art phasing methods, on two large datasets: > 8 million diverse, research-consented 23andMe, Inc. customers and the UK Biobank (UKB). We find that both perform exceptionally well. Beagle's median switch error rate (SER) (after excluding single SNP switches) in white British trios from UKB is 0.026% compared to 0.00% for European ancestry 23andMe research participants; 55.6% of European ancestry 23andMe research participants have zero non-single SNP switches, compared to 42.4% of white British trios. South Asian ancestry 23andMe research participants have the highest median SER amongst the 23andMe populations, but it is still remarkably low at 0.46%. We also investigate the relationship between identity-by-descent (IBD) and SER, finding that switch errors tend to occur in regions of little or no IBD segment coverage. SHAPEIT and Beagle excel at 'intra-chromosomal' phasing, but lack the ability to phase across chromosomes, motivating us to develop an inter-chromosomal phasing method, called HAPTIC ( HAP lotype TI ling and C lustering), that assigns paternal and maternal variants discretely genome-wide. Our approach uses identity-by-descent (IBD) segments to phase blocks of variants on different chromosomes. HAPTIC represents the segments a focal individual shares with their relatives as nodes in a signed graph and performs bipartite clustering on the signed graph using spectral clustering. We test HAPTIC on 1022 UKB trios, yielding a median phase error of 0.08% in regions covered by IBD segments (33.5% of sites). We also ran HAPTIC in the 23andMe database and found a median phase error rate (the rate of mismatching alleles between the inferred and true phase) of 0.92% in Europeans (93.8% of sites) and 0.09% in admixed Africans (92.7% of sites). HAPTIC's precision depends heavily on data from relatives, so will increase as datasets grow larger and more diverse. HAPTIC enables analyses that require the parent-of-origin of variants, such as association studies and ancestry inference of untyped parents.

Global and Local Ancestry and its Importance: A Review

The fastest way to significantly change the composition of a population is through admixture, an evolutionary mechanism. In animal breeding history, genetic admixture has provided both short-term and long-term advantages by utilizing the phenomenon of complementarity and heterosis in several traits and genetic diversity, respectively. The traditional method of admixture analysis by pedigree records has now been replaced greatly by genome-wide marker data that enables more precise estimations. Among these markers, SNPs have been the popular choice since they are cost-effective, not so laborious, and automation of genotyping is easy. Certain markers can suggest the possibility of a population's origin from a sample of DNA where the source individual is unknown or unwilling to disclose their lineage, which are called Ancestry-Informative Markers (AIMs). Revealing admixture level at the locus-specific level is termed as local ancestry and can be exploited to identify signs of recent selective response and can account for genetic drift. Considering the importance of genetic admixture and local ancestry, in this mini-review, both concepts are illustrated, encompassing basics, their estimation/identification methods, tools/software used and their applications.

MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Understanding the genomic heterogeneity of North African Imazighen: from broad to microgeographical perspectives

The strategic location of North Africa has led to cultural and demographic shifts, shaping its genetic structure. Historical migrations brought different genetic components that are evident in present-day North African genomes, along with autochthonous components. The Imazighen (plural of Amazigh) are believed to be the descendants of autochthonous North Africans and speak various Amazigh languages, which belong to the Afro-Asiatic language family. However, the arrival of different human groups, especially during the Arab conquest, caused cultural and linguistic changes in local populations, increasing their heterogeneity. We aim to characterize the genetic structure of the region, using the largest Amazigh dataset to date and other reference samples. Our findings indicate microgeographical genetic heterogeneity among Amazigh populations, modeled by various admixture waves and different effective population sizes. A first admixture wave is detected group-wide around the twelfth century, whereas a second wave appears in some Amazigh groups around the nineteenth century. These events involved populations with higher genetic ancestry from south of the Sahara compared to the current North Africans. A plausible explanation would be the historical trans-Saharan slave trade, which lasted from the Roman times to the nineteenth century. Furthermore, our investigation shows that assortative mating in North Africa has been rare.

Single-nucleotide variant calling in single-cell sequencing data with Monopogen

Single-cell omics technologies enable molecular characterization of diverse cell types and states, but how the resulting transcriptional and epigenetic profiles depend on the cell's genetic background remains understudied. We describe Monopogen, a computational tool to detect single-nucleotide variants (SNVs) from single-cell sequencing data. Monopogen leverages linkage disequilibrium from external reference panels to identify germline SNVs and detects putative somatic SNVs using allele cosegregating patterns at the cell population level. It can identify 100 K to 3 M germline SNVs achieving a genotyping accuracy of 95%, together with hundreds of putative somatic SNVs. Monopogen-derived genotypes enable global and local ancestry inference and identification of admixed samples. It identifies variants associated with cardiomyocyte metabolic levels and epigenomic programs. It also improves putative somatic SNV detection that enables clonal lineage tracing in primary human clonal hematopoiesis. Monopogen brings together population genetics, cell lineage tracing and single-cell omics to uncover genetic determinants of cellular processes.

Genomic signatures of strawberry domestication and diversification

Cultivated strawberry (Fragaria × ananassa) has a brief history of less than 300 yr, beginning with the hybridization of octoploids Fragaria chiloensis and Fragaria virginiana. Here we explored the genomic signatures of early domestication and subsequent diversification for different climates using whole-genome sequences of 289 wild, heirloom, and modern varieties from two major breeding programs in the United States. Four nonadmixed wild octoploid populations were identified, with recurrent introgression among the sympatric populations. The proportion of F. virginiana ancestry increased by 20% in modern varieties over initial hybrids, and the proportion of F. chiloensis subsp. pacifica rose from 0% to 3.4%. Effective population size rapidly declined during early breeding. Meanwhile, divergent selection for distinct environments reshaped wild allelic origins in 21 out of 28 chromosomes. Overlapping divergent selective sweeps in natural and domesticated populations revealed 16 convergent genomic signatures that may be important for climatic adaptation. Despite 20 breeding cycles since initial hybridization, more than half of loci underlying yield and fruit size are still not under artificial selection. These insights add clarity to the domestication and breeding history of what is now the most widely cultivated fruit in the world.