Comparing the effect of imputation reference panel composition in four distinct Latin American cohorts.

Genome-wide association studies have been useful in identifying genetic risk factors for various phenotypes. These studies rely on imputation and many existing panels are largely composed of individuals of European ancestry, resulting in lower levels of imputation quality in underrepresented populations. We aim to analyze how the composition of imputation reference panels affects imputation quality in four target Latin American cohorts. We compared imputation quality for chromosomes 7 and X when altering the imputation reference panel by: 1) increasing the number of Latin American individuals; 2) excluding either Latin American, African, or European individuals, or 3) increasing the Indigenous American (IA) admixture proportions of included Latin Americans. We found that increasing the number of Latin Americans in the reference panel improved imputation quality in the four populations; however, there were differences between chromosomes 7 and X in some cohorts. Excluding Latin Americans from analysis resulted in worse imputation quality in every cohort, while differential effects were seen when excluding Europeans and Africans between and within cohorts and between chromosomes 7 and X. Finally, increasing IA-like admixture proportions in the reference panel increased imputation quality at different levels in different populations. The difference in results between populations and chromosomes suggests that existing and future reference panels containing Latin American individuals are likely to perform differently in different Latin American populations.

Pharmacogenetic variability of tuberculosis biomarkers in native and mestizo Peruvian populations.

In Peru, 29 292 people were diagnosed with tuberculosis in 2022. Although tuberculosis treatments are effective, 3.4%-13% are associated with significant adverse drug reactions, with drug-induced liver injury (DILI) considered the most predominant. Among the first-line antituberculosis drugs, isoniazid is the main drug responsible for the appearance of DILI. In liver, isoniazid (INH) is metabolized by N-acetyltransferase-2 (NAT2) and cytochrome P450 2E1 (CYP2E1). Limited information exists on genetic risk factors associated with the presence of DILI to antituberculosis drugs in Latin America, and even less is known about these factors in the native and mestizo Peruvian population. The aim of this study was to determine the prevalence of NAT2 and CYP2E1 genotypes in native and mestizo population. An analytical cross-sectional analysis was performed using genetic data from mestizo population in Lima and native participants from south of Peru. NAT2 metabolizer was determined as fast, intermediate and slow, and CYP2E1 genotypes were classified as c1/c1, c1/c2 and c2/c2, from molecular tests and bioinformatic analyses. Of the 472 participants, 36 and 6 NAT2 haplotypes were identified in the mestizo and native population, respectively. In mestizo population, the most frequent NAT2*5B and NAT2*7B haplotypes were associated with DILI risk; while in natives, NAT2*5G and NAT2*13A haplotypes were associated with decreased risk of DILI. For CYP2E1, c1/c1 and c1/c2 genotypes are the most frequent in natives and mestizos, respectively. The linkage disequilibrium of NAT2 single nucleotide polymorphisms (SNPs) was estimated, detecting a block between all SNPs natives. In addition, a block between rs1801280 and rs1799929 for NAT2 was detected in mestizos. Despite the limitations of a secondary study, it was possible to report associations between NAT2 and CYP2E alleles with Peruvian native and mestizo by prevalence ratios. The results of this study will help the development of new therapeutic strategies for a Tuberculosis efficient control between populations.

Genome-wide study investigating effector genes and polygenic prediction for kidney function in persons with ancestry from Africa and the Americas.

Chronic kidney disease is a leading cause of death and disability globally and impacts individuals of African ancestry (AFR) or with ancestry in the Americas (AMS) who are under-represented in genome-wide association studies (GWASs) of kidney function. To address this bias, we conducted a large meta-analysis of GWASs of estimated glomerular filtration rate (eGFR) in 145,732 AFR and AMS individuals. We identified 41 loci at genome-wide significance (p < 5 × 10-8), of which two have not been previously reported in any ancestry group. We integrated fine-mapped loci with epigenomic and transcriptomic resources to highlight potential effector genes relevant to kidney physiology and disease, and reveal key regulatory elements and pathways involved in renal function and development. We demonstrate the varying but increased predictive power offered by a multi-ancestry polygenic score for eGFR and highlight the importance of population diversity in GWASs and multi-omics resources to enhance opportunities for clinical translation for all.

Unappreciated subcontinental admixture in Europeans and European Americans and implications for genetic epidemiology studies.

European-ancestry populations are recognized as stratified but not as admixed, implying that residual confounding by locus-specific ancestry can affect studies of association, polygenic adaptation, and polygenic risk scores. We integrate individual-level genome-wide data from ~19,000 European-ancestry individuals across 79 European populations and five European American cohorts. We generate a new reference panel that captures ancestral diversity missed by both the 1000 Genomes and Human Genome Diversity Projects. Both Europeans and European Americans are admixed at the subcontinental level, with admixture dates differing among subgroups of European Americans. After adjustment for both genome-wide and locus-specific ancestry, associations between a highly differentiated variant in LCT (rs4988235) and height or LDL-cholesterol were confirmed to be false positives whereas the association between LCT and body mass index was genuine. We provide formal evidence of subcontinental admixture in individuals with European ancestry, which, if not properly accounted for, can produce spurious results in genetic epidemiology studies.

Pharmacogenetics: ethnicity, treatment and health in Latin American populations.

The Ibero-American Network of Pharmacogenetics and Pharmacogenomics (RIBEF) studies Latin American populations to benefit from the implementation of personalized medicine. Since 2006, it has studied ethnicity to apply pharmacogenetics knowledge in autochthonous populations of Latin America, considering ancestral medicine. The meeting 'Pharmacogenetics: ethnicity, Treatment and Health in Latin American Populations' was held in Mexico City, Mexico, and presented the relevance of RIBEF collaboration with Latin American researchers and the governments of Mexico, Spain and the Autonomous Community of Extremadura. The results of 17 years of uninterrupted work by RIBEF, the Declaration of Mérida/T'Hó and the call for the Dr José María Cantú Award for studies focused on the pharmacogenetics of native populations in Latin America were presented.

X-Chromosome Association Study in Latin American Cohorts Identifies New Loci in Parkinson's Disease.

BACKGROUND: Sex differences in Parkinson's disease (PD) risk are well-known. However, the role of sex chromosomes in the development and progression of PD is still unclear. OBJECTIVE: The objective of this study was to perform the first X-chromosome-wide association study for PD risk in a Latin American cohort. METHODS: We used data from three admixed cohorts: (1) Latin American Research consortium on the Genetics of Parkinson's Disease (n = 1504) as discover cohort, and (2) Latino cohort from International Parkinson Disease Genomics Consortium (n = 155) and (3) Bambui Aging cohort (n = 1442) as replication cohorts. We also developed an X-chromosome framework specifically designed for admixed populations. RESULTS: We identified eight linkage disequilibrium regions associated with PD. We replicated one of these regions (top variant rs525496; discovery odds ratio [95% confidence interval]: 0.60 [0.478-0.77], P = 3.13 × 10-5 replication odds ratio: 0.60 [0.37-0.98], P = 0.04). rs5525496 is associated with multiple expression quantitative trait loci in brain and non-brain tissues, including RAB9B, H2BFM, TSMB15B, and GLRA4, but colocalization analysis suggests that rs5525496 may not mediate risk by expression of these genes. We also replicated a previous X-chromosome-wide association study finding (rs28602900), showing that this variant is associated with PD in non-European populations. CONCLUSIONS: Our results reinforce the importance of including X-chromosome and diverse populations in genetic studies. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Clinical genomics and precision medicine.

Precision Medicine emerges from the genomic paradigm of health and disease. For precise molecular diagnoses of genetic diseases, we must analyze the Whole Exome (WES) or the Whole Genome (WGS). By not needing exon capture, WGS is more powerful to detect single nucleotide variants and copy number variants. In healthy individuals, we can observe monogenic highly penetrant variants, which may be causally responsible for diseases, and also susceptibility variants, associated with common polygenic diseases. But there is the major problem of penetrance. Thus, there is the question of whether it is worthwhile to perform WGS in all healthy individuals as a step towards Precision Medicine. The genetic architecture of disease is consistent with the fact that they are all polygenic. Moreover, ancestry adds another layer of complexity. We are now capable of obtaining Polygenic Risk Scores for all complex diseases using only data from new generation sequencing. Yet, review of available evidence does not at present favor the idea that WGS analyses are sufficiently developed to allow reliable predictions of the risk components for monogenic and polygenic hereditary diseases in healthy individuals. Probably, it is still better for WGS to remain reserved for the diagnosis of pathogenic variants of Mendelian diseases.

GSTT1/GSTM1 Genotype and Anti-Tuberculosis Drug-Induced Hepatotoxicity in Peruvian Patients.

In Peru, 24,581 people were diagnosed with tuberculosis (TB) in 2020. Although TB treatments are effective, 3.4-13% are associated with significant adverse drug reactions (ADRs), with drug-induced liver injury (DILI) considered the most predominant. Among the first-line antituberculosis drugs, isoniazid (INH) is the main drug responsible for the appearance of DILI. In the liver, INH is metabolized by the enzymes N-acetyltransferase-2 (NAT2), cytochrome P450 2E1 (CYP2E1), and glutathione S-transferase (GST) with two isoforms, GSTT1 and GSTM1. Based on previous studies, we hypothesized that interactions between the GSTT1 and GSTM1 null genotypes induce DILI in TB patients. In this cross-sectional study of 377 participants who completed their anti-TB treatment, we genotyped by revealing the presence or absence of 215- and 480-bp bands of GSTM1 and GSTT1, respectively. We found that the prevalence of the GSTM1 genotype was 52.79% and 47.21% for presence and null, respectively, and for GSTT1 it was 69.76% and 30.24% for presence and null, respectively. Neither genotype was prevalent in the patients who developed DILI (n = 16). We did not confirm our hypothesis; however, we found that the combination of GSTM1 present genotype, GSTT1 null genotype, fast NAT2 acetylators, and CYP2E1 c1/c1 genotype had a significant risk for the development of ADR (OR 11; p = 0.017; 95% CI: (0.54-186.35)). We propose that the presence of the GSTM1 present genotype, GSTT1 null genotype, fast NAT2 acetylators, and CYP2E1 c1/c1 genotype in the Peruvian population could be considered a risk factor for the development of ADR due to therapeutic drug intake.

Genetic analyses of the electrocardiographic QT interval and its components identify additional loci and pathways.

The QT interval is an electrocardiographic measure representing the sum of ventricular depolarization and repolarization, estimated by QRS duration and JT interval, respectively. QT interval abnormalities are associated with potentially fatal ventricular arrhythmia. Using genome-wide multi-ancestry analyses (>250,000 individuals) we identify 177, 156 and 121 independent loci for QT, JT and QRS, respectively, including a male-specific X-chromosome locus. Using gene-based rare-variant methods, we identify associations with Mendelian disease genes. Enrichments are observed in established pathways for QT and JT, and previously unreported genes indicated in insulin-receptor signalling and cardiac energy metabolism. In contrast for QRS, connective tissue components and processes for cell growth and extracellular matrix interactions are significantly enriched. We demonstrate polygenic risk score associations with atrial fibrillation, conduction disease and sudden cardiac death. Prioritization of druggable genes highlight potential therapeutic targets for arrhythmia. Together, these results substantially advance our understanding of the genetic architecture of ventricular depolarization and repolarization.

Global and local ancestry modulate APOE association with Alzheimer's neuropathology and cognitive outcomes in an admixed sample.

Dementia is more prevalent in Blacks than in Whites, likely due to a combination of environmental and biological factors. Paradoxically, clinical studies suggest an attenuation of APOE ε4 risk of dementia in African ancestry (AFR), but a dearth of neuropathological data preclude the interpretation of the biological factors underlying these findings, including the association between APOE ε4 risk and Alzheimer's disease (AD) pathology, the most frequent cause of dementia. We investigated the interaction between African ancestry, AD-related neuropathology, APOE genotype, and functional cognition in a postmortem sample of 400 individuals with a range of AD pathology severity and lack of comorbid neuropathology from a cohort of community-dwelling, admixed Brazilians. Increasing proportions of African ancestry (AFR) correlated with a lower burden of neuritic plaques (NP). However, for individuals with a severe burden of NP and neurofibrillary tangles (NFT), AFR proportion was associated with worse Clinical Dementia Rating sum of boxes (CDR-SOB). Among APOE ε4 carriers, the association between AFR proportion and CDR-SOB disappeared. APOE local ancestry inference of a subset of 309 individuals revealed that, in APOE ε4 noncarriers, non-European APOE background correlated with lower NP burden and, also, worse cognitive outcomes than European APOE when adjusting by NP burden. Finally, APOE ε4 was associated with worse AD neuropathological burden only in a European APOE background. APOE genotype and its association with AD neuropathology and clinical pattern are highly influenced by ancestry, with AFR associated with lower NP burden and attenuated APOE ε4 risk compared to European ancestry.

Identifying signatures of natural selection in Indian populations.

In this study, we present the results of a genome-wide scan for signatures of positive selection using data from four tribal groups (Kokana, Warli, Bhil, and Pawara) and two caste groups (Deshastha Brahmin and Kunbi Maratha) from West of the Maharashtra State In India, as well as two samples of South Asian ancestry from the 1KG project (Gujarati Indian from Houston, Texas and Indian Telugu from UK). We used an outlier approach based on different statistics, including PBS, xpEHH, iHS, CLR, Tajima's D, as well as two recently developed methods: Graph-aware Retrieval of Selective Sweeps (GRoSS) and Ascertained Sequentially Markovian Coalescent (ASMC). In order to minimize the risk of false positives, we selected regions that are outliers in all the samples included in the study using more than one method. We identified putative selection signals in 107 regions encompassing 434 genes. Many of the regions overlap with only one gene. The signals observed using microarray-based data are very consistent with our analyses using high-coverage sequencing data, as well as those identified with a novel coalescence-based method (ASMC). Importantly, at least 24 of these genomic regions have been identified in previous selection scans in South Asian populations or in other population groups. Our study highlights genomic regions that may have played a role in the adaptation of anatomically modern humans to novel environmental conditions after the out of Africa migration.

NAT2 and CYP2E1 polymorphisms and antituberculosis drug-induced hepatotoxicity in Peruvian patients.

BACKGROUND: In Peru, 32,970 people were diagnosed with tuberculosis (TB) in 2019. Although TB treatment is effective, 3.4%-13% is associated with significant adverse drug reactions (ADR), considering drug-induced liver injury (DILI) as the most prevalent. Among the first-line anti-TB drugs, isoniazid (INH) is primarily responsible for the occurrence of DILI. INH is metabolized in the liver by the enzymes N-acetyltransferase-2 (NAT2) and Cytochrome P450 2E1 (CYP2E1). Based on the previous studies, we hypothesized that the interactions between slow CYP2E1 genotype and NAT2 slow acetylators will induce DILI in TB patients. METHODS: In this cross-sectional study, all 377 participants completed their anti-TB treatment, and we genotyped SNPs: rs1041983, rs1801280, rs1799929, rs1799930, rs1208, and rs1799931 for NAT2 and rs3813867 and rs2031920 for CYP2E1. RESULTS: We found that rapid, intermediate, and slow NAT2 acetylator were 15%, 38%, and 47%, respectively, in the general population. Intermediate NAT2 acetylator is the least prevalent among patients with adverse reactions (p = 0.024). We did not confirm our hypothesis, however, we found that the combination of intermediate NAT2 acetylators and CYP2E1 c1/c1 genotype significantly protected (OR = 0.16; p = 0.049) against the development of DILI in our population. CONCLUSION: We propose that the presence of NAT2 intermediate and CYP2E1 c1/c1 genotype could help in therapeutic drug monitoring, and optimize its therapeutic benefits while minimizing its risk for side effects or toxicity.

NAToRA, a relatedness-pruning method to minimize the loss of dataset size in genetic and omics analyses.

Genetic and omics analyses frequently require independent observations, which is not guaranteed in real datasets. When relatedness cannot be accounted for, solutions involve removing related individuals (or observations) and, consequently, a reduction of available data. We developed a network-based relatedness-pruning method that minimizes dataset reduction while removing unwanted relationships in a dataset. It uses node degree centrality metric to identify highly connected nodes (or individuals) and implements heuristics that approximate the minimal reduction of a dataset to allow its application to complex datasets. When compared with two other popular population genetics methodologies (PLINK and KING), NAToRA shows the best combination of removing all relatives while keeping the largest possible number of individuals in all datasets tested and also, with similar effects on the allele frequency spectrum and Principal Component Analysis than PLINK and KING. NAToRA is freely available, both as a standalone tool that can be easily incorporated as part of a pipeline, and as a graphical web tool that allows visualization of the relatedness networks. NAToRA also accepts a variety of relationship metrics as input, which facilitates its use. We also release a genealogies simulator software used for different tests performed in this study.

Population genetics of PDE4B (phosphodiesterase-4B) in neglected Native Americans: Implications for cancer pharmacogenetics.

PDE4B (phosphodiesterase-4B) has an important role in cancer and in pharmacology of some disorders, such as inflammatory diseases. Remarkably in Native Americans, PDE4B variants are associated with acute lymphoblastic leukemia (ALL) relapse, as this gene modulates sensitivity of glucocorticoids used in ALL chemotherapy. PDE4B allele rs6683977.G, associated with genomic regions of Native American origin in US-Hispanics (admixed among Native Americans, Europeans, and Africans), increases ALL relapse risk, contributing to an association between Native American ancestry and ALL relapse that disappeared with an extra-phase of chemotherapy. This result insinuates that indigenous populations along the Americas may have high frequencies of rs6683977.G, but this has never been corroborated. We studied ancestry and PDE4B diversity in 951 healthy individuals from nine Latin American populations. In non-admixed Native American populations rs6683977.G has frequencies greater than 90%, is in linkage disequilibrium with other ALL relapse associated and regulatory variants in PDE4B-intron-7, conforming haplotypes showing their highest worldwide frequencies in Native Americans (>0.82). Our findings inform the discussion on the pertinence of an extra-phase of chemotherapy in Native American populations, and exemplifies how knowledge generated in US-Hispanics is relevant for their even more neglected and vulnerable Native American ancestors along the American continent.

Whole-genome sequencing of 1,171 elderly admixed individuals from São Paulo, Brazil.

As whole-genome sequencing (WGS) becomes the gold standard tool for studying population genomics and medical applications, data on diverse non-European and admixed individuals are still scarce. Here, we present a high-coverage WGS dataset of 1,171 highly admixed elderly Brazilians from a census-based cohort, providing over 76 million variants, of which ~2 million are absent from large public databases. WGS enables identification of ~2,000 previously undescribed mobile element insertions without previous description, nearly 5 Mb of genomic segments absent from the human genome reference, and over 140 alleles from HLA genes absent from public resources. We reclassify and curate pathogenicity assertions for nearly four hundred variants in genes associated with dominantly-inherited Mendelian disorders and calculate the incidence for selected recessive disorders, demonstrating the clinical usefulness of the present study. Finally, we observe that whole-genome and HLA imputation could be significantly improved compared to available datasets since rare variation represents the largest proportion of input from WGS. These results demonstrate that even smaller sample sizes of underrepresented populations bring relevant data for genomic studies, especially when exploring analyses allowed only by WGS.

Pharmacogenetics research in Brazil: a systematic review.

Aim: Pharmacogenomics (PGx) is a rising scientific area in many countries, such as Brazil. Objectives: To identify biomarkers, therapeutic areas, probe drugs and regions/ethnicities most studied in the country in order to guide future studies. Materials & methods: Systematic review of 1060 studies (from 1968 to 2020) comprising 80 genes, six probe drugs and 3,819,233 individuals. Results:MTHFR and HLA-A/B were the most studied genes and metoprolol and dextromethorphan the most studied probe drugs. Oncology was the most studied therapeutic area considering PGx biomarkers. The country's regions and ethnic groups were studied unevenly, with south/southeast and White people over-represented in respect to their demographic relevance, in detriment of the center-west/northeast/north and Black/mixed individuals. Conclusion: Many of the gaps and possible paths to be covered to reach even PGx data are pointed out by this review.

HLA-G 3'UTR haplotype frequencies in highland and lowland South Native American populations.

Human Leukocyte Antigen (HLA)-G participates in several biological processes, including reproduction, vascular remodeling, immune tolerance, and hypoxia response. HLA-G is a potential candidate gene for high altitude adaptation since its expression is modulated in both micro and macro environment under hypoxia and constant cellular stress. Besides the promoter region, the HLA-G 3'untranslated region (UTR) influences HLA-G expression patterns through several post-transcriptional mechanisms. Currently, the 3'UTR genetic diversity in terms of altitude adaptation of Native American populations is still unexplored, particularly at high altitude ecoregions. Here, we evaluated 288 Native Americans from 9 communities located in the Andes [highland (HL); ≥2,500 m (range = 2,838-4,433 m)] and 8 populations located in lowland (LL) regions [<2,500 m (range = 80-431 m); Amazonian tropical forest, Brazilian central plateau, and Chaco] of South America. In total, nine polymorphic sites and ten haplotypes were observed. The most frequent haplotypes (UTR-1, UTR-2, and UTR-3) accounted for âˆ¼ 77% of haplotypes found in LL, while in the HL, the same haplotypes reach âˆ¼ 93%. Also, a remarkable high frequency of putative ancestral UTR-5 haplotype was observed in LL (21.5%), while in HL UTR-2 reaches up to 47%. Further, UTR-2 frequency positively correlates with altitude-related variables, while a negative correlation for UTR-5 was observed. From an evolutionary perspective, we observed a tendency towards balancing selection in HL and LL populations thus suggesting that haplotypes of ancient and more derived alleles may have been co-opted for relatively recent adaptations such as those experienced by modern humans in the highland and lowland of South America. We also discuss how long-term balancing selection can be a reservoir of genetic variants that can be positively selected. Finally, our study provides some pieces of evidence that HLA-G 3'UTR haplotypes may have contributed to high altitude adaptation in the Andes.

Tracing the Distribution of European Lactase Persistence Genotypes Along the Americas.

In adulthood, the ability to digest lactose, the main sugar present in milk of mammals, is a phenotype (lactase persistence) observed in historically herder populations, mainly Northern Europeans, Eastern Africans, and Middle Eastern nomads. As the -13910∗T allele in the MCM6 gene is the most well-characterized allele responsible for the lactase persistence phenotype, the -13910C > T (rs4988235) polymorphism is commonly evaluated in lactase persistence studies. Lactase non-persistent adults may develop symptoms of lactose intolerance when consuming dairy products. In the Americas, there is no evidence of the consumption of these products until the arrival of Europeans. However, several American countries' dietary guidelines recommend consuming dairy for adequate human nutrition and health promotion. Considering the extensive use of dairy and the complex ancestry of Pan-American admixed populations, we studied the distribution of -13910C > T lactase persistence genotypes and its flanking haplotypes of European origin in 7,428 individuals from several Pan-American admixed populations. We found that the -13910∗T allele frequency in Pan-American admixed populations is directly correlated with allele frequency of the European sources. Moreover, we did not observe any overrepresentation of European haplotypes in the -13910C > T flanking region, suggesting no selective pressure after admixture in the Americas. Finally, considering the dominant effect of the -13910∗T allele, our results indicate that Pan-American admixed populations are likely to have higher frequency of lactose intolerance, suggesting that general dietary guidelines deserve further evaluation across the continent.

Human-SARS-CoV-2 interactome and human genetic diversity: TMPRSS2-rs2070788, associated with severe influenza, and its population genetics caveats in Native Americans.

For human/SARS-CoV-2 interactome genes ACE2, TMPRSS2 and BSG, there is a convincing evidence of association in Asians with influenza-induced SARS for TMPRSS2-rs2070788, tag-SNP of the eQTL rs383510. This case illustrates the importance of population genetics and of sequencing data in the design of genetic association studies in different human populations: the high linkage disequilibrium (LD) between rs2070788 and rs383510 is Asian-specific. Leveraging on a combination of genotyping and sequencing data for Native Americans (neglected in genetic studies), we show that while their frequencies of the Asian tag-SNP rs2070788 is, surprisingly, the highest worldwide, it is not in LD with the eQTL rs383510, that therefore, should be directly genotyped in genetic association studies of SARS in populations with Native American ancestry.